Télécharger gdef2.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : gdef2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ********************************************************************
  5. * *
  6. * TEST GDEF2 *
  7. * *
  8. * CISAILLEMENT PUR EN GRANDES DEFORMATIONS PLASTIQUES *
  9. * ( MODELE DE PLASTICITE ISOTROPE ) *
  10. * On compare avec les valeurs obtenues a une solution *
  11. * analytique *
  12. * *
  13. ********************************************************************
  14. *
  15. OPTI ECHO 1;
  16. TITRE ' TEST GDEF2 ' ;
  17. GRAPH = 'N';
  18. OPTI DIME 2 ELEM QUA4 MODE PLAN CONT;
  19. *
  20. *==================================================================*
  21. * geometrie : maillage *
  22. *==================================================================*
  23. *
  24. P1 = 0 0 ;
  25. P2 = 1. 0. ;
  26. P3 = 1. 1. ;
  27. P4 = 0. 1. ;
  28. L1 = P1 D 1 P2 ;
  29. L2 = P2 D 1 P3;
  30. L3 = P3 D 1 P4 ;
  31. L4 = P4 D 1 P1 ;
  32. SU = DALLER L1 L2 L3 L4 ;
  33. SI (NEG GRAPH 'N');
  34. TRAC SU ;
  35. FINSI;
  36. *
  37. *==================================================================*
  38. * modele - materiau - caracteristiques en Pa *
  39. * conditions aux limites *
  40. *==================================================================*
  41. *
  42. MO = MODE SU MECANIQUE ELASTIQUE PLASTIQUE ISOTROPE
  43. epsilon jaumann ;
  44. YU = 1000. ; XNU = 0.3 ;
  45. LTRAC = EVOL MANU 'EPS' ( PROG 0. 1.E-3 10.001 )
  46. 'SIG' ( PROG 0. 1. 30.91026920);
  47. LEPS1 = LTRAC extr absc ;
  48. LSIG1 = LTRAC extr ordo ;
  49. lsig1 = lsig1 enle 1 ;
  50. leps1 = (leps1 enle 1) - (lsig1 / yu) ;
  51. LECRO = evol vert manu eps leps1 sig lsig1 ;
  52. si (neg graph 'N')
  53. dess (ltrac et lecro) titr ' Courbes de traction et d ecrouissage (vert)' ;
  54. fins ;
  55. MA = MATE MO YOUN YU NU XNU ECRO LECRO ;
  56. BL1 = BLOQ DEPL P1 ;
  57. BL2 = BLOQ DEPL P2 ;
  58. BL3 = BLOQ UX L3 ;
  59. BL4 = BLOQ UY L3 ;
  60. FF = DEPI BL3 1. ;
  61. *
  62. *==================================================================*
  63. * Initialisation de la table pour appel a PASAPAS *
  64. *==================================================================*
  65. *
  66. TAB1 = TABLE;
  67. tab1.grands_deplacements=vrai;
  68. TAB1 . BLOCAGES_MECANIQUES= BL1 ET BL2 ET BL3 ET BL4 ;
  69. TAB1 . MODELE = MO;
  70. TAB1 . CARACTERISTIQUES = MA ;
  71. EV1 = EVOL MANU X (PROG 0 5 10 15 20 )
  72. Y (PROG 0 5 10 5 0 );
  73. TAB1 . CHARGEMENT = CHAR MECA FF EV1 ;
  74. TAB1 . TEMPS_CALCULES = PROG 0.1 PAS 0.1 5.;
  75. PASAPAS TAB1 ;
  76. *
  77. *==================================================================*
  78. * affichage des resultats *
  79. *==================================================================*
  80. *
  81. MESS ' RESULTATS TEST GDEF2 ';
  82. MESS ' -------------------- ';
  83. MESS ' On cherche le cisaillement ';
  84. MESS ' Il est compare a une solution analytique ';
  85. *
  86. *----------- r{sultats --------------------------------------------------
  87. *
  88. AA = TAB1 . DEPLACEMENTS ;
  89. BB = TAB1 . CONTRAINTES ;
  90. N = DIME AA ;
  91. SP12= PROG ;
  92. SPTH12= PROG ;
  93. EP = PROG ;
  94. EP = TAB1 . TEMPS_CALCULES ;
  95. REPE BOU (N - 1 ) ;
  96. CSX = CHAN CHPO MO ( BB . &BOU ) ;
  97. DX = EXTR ( AA . &BOU ) UX P3 ;
  98. S12 = EXTR CSX SMXY P1 ;
  99. SP12 = SP12 ET ( PROG S12 ) ;
  100. SDX = 0.577350269 + DX ;
  101. SPTH12= SPTH12 ET ( PROG SDX ) ;
  102. FIN BOU ;
  103. EV12 = EVOL ROSE MANU 'DEPLA' EP 'SIG12' SP12 ;
  104. FVTH12 = EVOL TURQ MANU 'DEPLA' EP 'SIG12' SPTH12 ;
  105. SI (NEG GRAPH 'N');
  106. DESS ( FVTH12 ET EV12) ;
  107. FINSI;
  108. *
  109. * COMPARAISON SOL ANALYTIQUE
  110. * ecart
  111. *
  112. ER12 =(MAXI ABS ( SP12 - SPTH12 ))/(MAXI ABS SPTH12);
  113. *
  114. MESS ' Ecart relatif en contraintes sig12 :' ER12 ;
  115. *
  116. *
  117. *==================================================================*
  118. * code fonctionnement *
  119. *==================================================================*
  120. *
  121. * L'ecart maximum entre valeurs de reference et calculee doit etre
  122. * inferieur a 5.d-2
  123. *
  124. SAUT 2 LIGNE;
  125. SI (ER12 <EG 0.05 ) ;
  126. ERRE 0;
  127. SINON;
  128. ERRE 5;
  129. FINSI;
  130. *
  131. *
  132. FIN;
  133.  
  134.  
  135.  
  136.  
  137.  
  138.  
  139.  
  140.  
  141.  
  142.  
  143.  
  144.  
  145.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales