Télécharger chaboche1.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : chaboche1.dgibi
  2. * section : Mecanique Viscoplastique
  3. *----------------------------------------------------------------------*
  4. * CHABOCHE1.DGIBI *
  5. *----------------------------------------------------------------------*
  6. *
  7. * Objet :
  8. * -------
  9. *
  10. * Test de validation d'une loi de comportement de materiau.
  11. * Loi de comportement elastoviscoplastique de Chaboche.
  12. *
  13. *
  14. * Description :
  15. * -------------
  16. * Essai de traction-compression.
  17. *
  18. * Type de calcul : Mecanique
  19. * Mode de calcul : 2D contraintes planes
  20. * Type d'element : QUA8 BBAR
  21. * Chargement : Deplacement impose
  22. * Reference : solution comparee a la loi ONERA
  23. * Objectif : Ecart relatif sur la contrainte SMZZ < 1.e-4
  24. *
  25. * Reference :
  26. * -----------
  27. * L&M : J. Lemaitre, J.-L. Chaboche, "Mecanique des materiaux solides",
  28. * Dunod, 2e edition, 1996.
  29. *
  30. *----------------------------------------------------------------------*
  31. *
  32. * Pour afficher les traces, mettre IGRAPH a vrai :
  33. IGRAPH = faux ;
  34. *
  35. *------------------------ Maillage (1 element) ------------------------*
  36. *
  37. opti dime 2 elem qua8 mode plan cont ;
  38.  
  39. ne1 = 1 ;
  40. l1 = (0 0) droi ne1 (1. 0.) ;
  41. S1 = l1 tran ne1 (0. 1.) ;
  42.  
  43. l2 = s1 cote 2 ;
  44. l3 = S1 cote 3 ;
  45. l4 = S1 cote 4 ;
  46.  
  47. si IGRAPH ;
  48. trac qual S1 titr ' Maillage ' ;
  49. fins ;
  50. *
  51. *---------------------- Modele & Carateristiques ----------------------*
  52. *
  53. * Par rapport aux parametres du modele presentes dans le L&M, p. 310 :
  54. * - CLM : coefficient C d'ecrouissage cinematique du L&M
  55. * - Gam : coefficient gamma d'ecrouissage cinematique du L&M
  56. * les coefficients d'ecrouissage du modele de Cast3M sont :
  57. * - C = Gam
  58. * - A = CLM / Gam
  59. * Valeurs des parametres du modele pour acier 316L a 20 degC (L&M, p. 319) :
  60. SIGY1 = 82.e6 ;
  61. N1 = 24. ;
  62. K1 = 151.e6 ;
  63. CLM = 162400.e6 ;
  64. Gam = 2800. ;
  65. A1 = CLM / Gam ;
  66. C1 = Gam ;
  67. Q1 = 142.e6 - SIGY1 ;
  68. B1 = 8. ;
  69.  
  70. * Module d'Young, coeff. de Poisson et epaisseur (DIM3) :
  71. YM1 = 200.e9 ;
  72. NU1 = 0.3 ;
  73. E1 = 1.e-3 ;
  74.  
  75. * Modeles & Caracteristiques :
  76. mo1 = mode s1 mecanique elastique isotrope viscoplastique chaboche bbar ;
  77. ma1 = mate mo1 'YOUN' YM1 'NU' NU1 'SIGY' SIGY1 'DIM3' E1
  78. 'K' K1 'N' N1 'A' A1 'C' C1 'B' B1 'Q' Q1 ;
  79.  
  80. mo2 = 'MODE' S1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 'VISCOPLASTIQUE' 'ONERA' 'BBAR' ;
  81. ma2 = 'MATE' mo2 'YOUN' YM1 'NU' NU1 'KK' SIGY1 'DIM3' E1
  82. 'K0' K1 'N' N1 'B' B1
  83. 'ALFK' 0. 'ALFR' 1. 'ALF' 0.
  84. 'A1' A1 'C1' C1 'BET1' 0. 'R1' 0.
  85. 'A2' 0. 'C2' 0. 'BET2' 0. 'R2' 0.
  86. 'PHI' 1. 'GAMA' 0. 'M' 0.
  87. 'QMAX' 1. 'QSTA' 0. 'MU' 0. 'ETA' 0. ;
  88.  
  89. *-------------------------- CL & Chargement ---------------------------*
  90. *
  91. cl1 = bloq uy l1 ;
  92. cl2 = bloq ux l4 ;
  93. cl3 = bloq uy l3 ;
  94.  
  95. u1 = 0.002 ;
  96. f1 = depi cl3 u1 ;
  97.  
  98. xt1 = 100. ;
  99. ev1 = evol manu temp (xt1 * (prog 0. 1. 2. 3. 4. 6.))
  100. (prog 0. 1. 0. -1. 0. 2.) ;
  101. cgu1 = char dimp f1 ev1 ;
  102.  
  103. si IGRAPH ;
  104. trac qual (S1 et (l1 coul roug) et (l3 coul vert) et (l4 coul bleu))
  105. titr ' CL : Uy=0 en L1 (rouge), Ux=0 en L4 (bleu) et Uy impose sur L3 (vert)' ;
  106. dess (u1*ev1) titre ' Evolution temporelle du deplacement impose ';
  107. fins ;
  108.  
  109. *------------------------- Resolution PASAPAS -------------------------*
  110. *
  111. * Modele de Chaboche :
  112. ta1 = tabl ;
  113. ta1.modele = mo1 ;
  114. ta1.caracteristiques = ma1 ;
  115. ta1.blocages_mecaniques = cl1 et cl2 et cl3 ;
  116. ta1.chargement = cgu1 ;
  117. ta1.temps_calcules = (prog 0. pas 0.1 6.) * xt1 ;
  118. *ta1.temps_calcules = prog 0. 1.0 ;
  119. ta1.hypothese_deformation = 'LINEAIRE' ;
  120.  
  121. ta1.mes_sauvegardes = tabl ;
  122. ta1.mes_sauvegardes.defin = vrai ;
  123.  
  124. pasapas ta1 ;
  125.  
  126. * Modele ONERA :
  127. ta2 = tabl ;
  128. ta2.modele = mo2 ;
  129. ta2.caracteristiques = ma2 ;
  130. ta2.blocages_mecaniques = cl1 et cl2 et cl3 ;
  131. ta2.chargement = cgu1 ;
  132. ta2.temps_calcules = ta1.temps_calcules ;
  133.  
  134. ta2.mes_sauvegardes = tabl ;
  135. ta2.mes_sauvegardes.defin = vrai ;
  136.  
  137. chq0 = manu chml mo2 epse 0. qqq Q1 stresses ;
  138. ta2.variables_internes = tabl ;
  139. ta2.variables_internes.0 = chq0 ;
  140.  
  141. pasapas ta2 ;
  142.  
  143. *-------------------- Postraitement des resultats ---------------------*
  144.  
  145. nb1 = dime ta1.temps ;
  146. ltps1 = prog ;
  147. leps1 = prog ;
  148. lsig1 = prog ;
  149. lsig2 = prog ;
  150. i1 = 0 ;
  151. repe bo1 nb1 ;
  152. tpsi1 = ta1.temps.i1 ;
  153. depi1 = ta1.deplacements.i1 ;
  154. epsi1 = (epsi mo1 ma1 depi1) exco epyy scal maxi ;
  155. sigi1 = ta1.contraintes.i1 exco smyy scal maxi ;
  156. sigi2 = ta2.contraintes.i1 exco smyy scal maxi ;
  157. ltps1 = ltps1 et (prog tpsi1) ;
  158. leps1 = leps1 et (prog epsi1) ;
  159. lsig1 = lsig1 et (prog sigi1) ;
  160. lsig2 = lsig2 et (prog sigi2) ;
  161. i1 = i1 + 1 ;
  162. fin bo1 ;
  163. lsig1 = 1.e-6 * lsig1 ;
  164. lsig2 = 1.e-6 * lsig2 ;
  165. evezz1 = evol manu temps ltps1 epzz leps1 ;
  166. evszz1 = evol roug manu temps ltps1 smzz lsig1 ;
  167. evszz2 = evol vert manu temps ltps1 smzz lsig2 ;
  168. evsve1 = evol roug manu epzz leps1 smzz lsig1 ;
  169. evsve2 = evol vert manu epzz leps1 smzz lsig2 ;
  170.  
  171. * Affichages :
  172. si IGRAPH ;
  173. tl1 = table ;
  174. tl1 . 1 = 'NOLI MARQ CROI' ;
  175. tl1 . titre = table ;
  176. tl1 . titre . 1 = 'ONERA' ;
  177. tl1 . titre . 2 = 'CHABOCHE' ;
  178. dess (evszz1 et evszz2) lege tl1
  179. titre ' Evolution temporelle de la contrainte de traction ' ;
  180. dess (evsve1 et evsve2) lege tl1 titre ' Courbe de traction ' ;
  181. fins ;
  182.  
  183. *------------------------- Test de validation -------------------------*
  184. *
  185. err0 = (maxi abs (lsig2 - lsig1)) / (maxi abs lsig2) ;
  186.  
  187. opti echo 0 ;
  188. mess ' Comparaison solution calculee / reference ' ;
  189. mess ' ----------------------------------------- ' ;
  190. mess ' Max.|SMZZ(cal.) - SMZZ(Ref.)| / Max.|SMZZ(cal.)| =' err0 ;
  191. opti echo 1 ;
  192.  
  193. si (err0 > 1.e-4) ;
  194. erre 5 ;
  195. fins ;
  196.  
  197. fin ;
  198.  
  199.  
  200.  
  201.  
  202.  
  203.  

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