Télécharger compar_syco_plast.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : compar_syco_plast.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *******************************************************
  5. **-------------------------------------------------------------------
  6. MESS ' COMPARAISON ' ;
  7. MESS ' CALCUL DYNAMIQUE COMPORTEMENT VISCOPLASTIQUE SYCO1 ';
  8. MESS ' avec parametre p de la loi SYCO1 valant 1 ';
  9. MESS ' CALCUL DYNAMIQUE COMPORTEMENT ELASTOPLASTIQUE ';
  10. MESS ' en 2D defo plane ' ;
  11. MESS ' Unites: N Pa m kg s' ;
  12. MESS '' ;
  13. **--------------------------------------------------------------------
  14. *** Options de calcul
  15. opti dime 2 elem QUA4 mode plan defo;
  16.  
  17. *******************************************************
  18. *** maillage
  19. pa1 = 0. 0.;
  20. pa2 = 5. 0.;
  21. pa3 = 5. 3.;
  22. pa4 = 0. 3.;
  23. li12 = pa1 droi 1 pa2;
  24. li23 = pa2 droi 1 pa3;
  25. li34 = pa3 droi 1 pa4;
  26. li41 = pa4 droi 1 pa1;
  27. su1 = dall li12 li23 li34 li41;
  28. elim 0.001 su1;
  29. *
  30.  
  31. **********************************************************
  32. *** MODELE & MATERIAU ***
  33. **********************************************************
  34.  
  35. Epsprog0 = ('PROG'
  36. 0. 0.01 0.1 0.5 );
  37. Sigprog0 = ('PROG'
  38. 0. 500. 1000. 1200.) * 1.e6;
  39. lcxx = 'EVOL' 'MANU' 'Epsilon' Epsprog0 'Sigma' Sigprog0 ;
  40. *
  41. ab1 = extr lcxx 'ABSC';
  42. ord1 = extr lcxx 'ORDO';
  43. eey1 = (extr ord1 2)/(extr ab1 2);
  44. *
  45. nu0 = 0.3 ;
  46. rho0= 7800.;
  47. **********************************************************
  48. *** Modele & materiau ***
  49. *
  50. * elements standards
  51. mod1 = 'MODELISER' su1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  52. 'VISCOPLASTIQUE' 'SYCO1' 'QUA4';
  53. mat1 = 'MATERIAU' mod1 'YOUN' eey1 'NU' nu0 'RHO' rho0
  54. 'TRAC' lcxx 'PSYC' 1. 'DSYC' 1.d8 ;
  55. *
  56.  
  57. **********************************************************
  58. *** CL et DEPLACEMENTS IMPOSES ***
  59. **********************************************************
  60.  
  61. cl1 = bloq (pa1 et pa4) 'UX';
  62. cl2 = bloq (pa1 et pa2) 'UY';
  63. cl3 = bloq (pa2 et pa3) 'UX';
  64. cltot = (cl1 'ET' cl2 'ET' cl3 ) ;
  65. *
  66. *... chargement en deplacement impose
  67. umax = 1.;
  68. uimp = depi cl3 umax;
  69.  
  70. *
  71. xx1 = prog 0. 1.;
  72. yy1 = prog 0. 50.;
  73. evochar = evol manu 'temps' xx1 'charge' yy1;
  74. char1 = char meca uimp evochar;
  75. *
  76. *
  77. *=======================================================*
  78. * *
  79. * RESOLUTION *
  80. * *
  81. *=======================================================*
  82.  
  83. xx0 = prog 0. 1.15e-3 1.25e-2;
  84. *
  85. tab1 = table;
  86. tab1.dynamique = vrai;
  87. tab1.modele = mod1;
  88. tab1.caracteristiques = mat1;
  89. tab1.blocages_mecaniques = cltot;
  90. tab1.chargement = char1;
  91. tab1.temps_calcules = xx0;
  92. tab1.precision = 1.e-5;
  93. pasapas tab1;
  94. *
  95. *post traitement
  96. sig1 = tab1.contraintes. 2;
  97. ut1 = tab1.deplacements. 2;
  98. vt1 = tab1.vitesses. 2;
  99. tt1 = tab1.temps. 2;
  100. var1 = tab1.variables_internes. 2;
  101. *
  102. *
  103. val1 = extr var1 'EPSE' 1 1 3;
  104. val2 = extr var1 'VP' 1 1 3;
  105. val3 = extr ut1 'UY' pa3;
  106.  
  107. *******************************************************
  108. *
  109. *******************************************************
  110. **-------------------------------------------------------------------
  111. MESS '' ;
  112. MESS ' CALCUL DYNAMIQUE COMPORTEMENT ELASTOPLASTIQUE ';
  113. MESS ' en 2D defo plane ' ;
  114. MESS ' Unites: N Pa m kg s' ;
  115. MESS '' ;
  116. **--------------------------------------------------------------------
  117. *** Options de calcul
  118. opti dime 2 elem QUA4 mode plan defo;
  119.  
  120. *******************************************************
  121. *** maillage
  122. pa1 = 0. 0.;
  123. pa2 = 5. 0.;
  124. pa3 = 5. 3.;
  125. pa4 = 0. 3.;
  126. li12 = pa1 droi 1 pa2;
  127. li23 = pa2 droi 1 pa3;
  128. li34 = pa3 droi 1 pa4;
  129. li41 = pa4 droi 1 pa1;
  130. su1 = dall li12 li23 li34 li41;
  131. elim 0.001 su1;
  132. *
  133.  
  134. **********************************************************
  135. *** MODELE & MATERIAU ***
  136. **********************************************************
  137.  
  138. Epsprog0 = ('PROG'
  139. 0. 0.01 0.1 0.5 );
  140. Sigprog0 = ('PROG'
  141. 0. 500. 1000. 1200.) * 1.e6;
  142. lcxx = 'EVOL' 'MANU' 'Epsilon' Epsprog0 'Sigma' Sigprog0 ;
  143. *
  144. ab1 = extr lcxx 'ABSC';
  145. ord1 = extr lcxx 'ORDO';
  146. eey1 = (extr ord1 2)/(extr ab1 2);
  147. *
  148. nu0 = 0.3 ;
  149. rho0= 7800.;
  150. **********************************************************
  151. *** Modele & materiau ***
  152. *
  153. * elements standards
  154. *
  155. mod1 = 'MODELISER' su1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  156. 'PLASTIQUE' ;
  157. mat1 = 'MATERIAU' mod1 'YOUN' eey1 'NU' nu0 'RHO' rho0
  158. 'TRAC' lcxx ;
  159.  
  160. **********************************************************
  161. *** CL et DEPLACEMENTS IMPOSES ***
  162. **********************************************************
  163.  
  164. cl1 = bloq (pa1 et pa4) 'UX';
  165. cl2 = bloq (pa1 et pa2) 'UY';
  166. cl3 = bloq (pa2 et pa3) 'UX';
  167. cltot = (cl1 'ET' cl2 'ET' cl3 ) ;
  168. *
  169. *... chargement en deplacement impose
  170. umax = 1.;
  171. uimp = depi cl3 umax;
  172.  
  173. *
  174. xx1 = prog 0. 1.;
  175. yy1 = prog 0. 50.;
  176. evochar = evol manu 'temps' xx1 'charge' yy1;
  177. char1 = char meca uimp evochar;
  178. *
  179. *=======================================================*
  180. * *
  181. * RESOLUTION *
  182. * *
  183. *=======================================================*
  184.  
  185. xx0 = prog 0. 1.15e-3 1.25e-2;
  186. *
  187. tab2 = table;
  188. tab2.dynamique = vrai;
  189. tab2.modele = mod1;
  190. tab2.caracteristiques = mat1;
  191. tab2.blocages_mecaniques = cltot;
  192. tab2.chargement = char1;
  193. tab2.temps_calcules = xx0;
  194. tab2.precision = 1.e-5;
  195. pasapas tab2;
  196. *
  197. *post traitement
  198. sig1 = tab2.contraintes. 2;
  199. ut1 = tab2.deplacements. 2;
  200. vt1 = tab2.vitesses. 2;
  201. tt1 = tab2.temps. 2;
  202. var1 = tab2.variables_internes. 2;
  203. *
  204. *
  205. xal1 = extr var1 'EPSE' 1 1 3;
  206. xal3 = extr ut1 'UY' pa3;
  207.  
  208. aa1 = ABS(val1 - xal1);
  209. aa3 = ABS(val3 - xal3);
  210. *
  211. si (aa1 >EG 1.e-8);
  212. erre 5;
  213. finsi;
  214. si (aa3 >EG 1.e-8);
  215. erre 5;
  216. finsi;
  217. fin;
  218.  
  219.  
  220.  
  221.  
  222.  
  223.  
  224.  

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