Télécharger guionnet_tra.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : guionnet_tra.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Guionnet_tra.dgibi: Jeux de données *
  5. * ---------------------------------------- *
  6. * *
  7. **************************************************
  8. * *
  9. * TEST DE VALIDATION D'UNE LOI DE COMPORTEMENT *
  10. * -------------------------------------------- *
  11. * DE MATERIAU *
  12. * ----------- *
  13. * *
  14. * LOI DE COMPORTEMENT VISCOPLASTIQUE DE : *
  15. * GUIONNET *
  16. * *
  17. * COMPARAISON DE CALCULS SUR DES ELEMENTS: *
  18. * - COQUE EPAISSE ( MFR=5 ) *
  19. * - MASSIF 3D ( MFR=1 ) *
  20. * *
  21. * STRUCTURE : UNE PLAQUE CARREE *
  22. * - COTE : 1. M *
  23. * - EPAISSEUR : 0.01 M *
  24. * *
  25. * CHARGEMENT : DEPLACEMENTS IMPOSES *
  26. * ESSAI DE TRACTION *
  27. * *
  28. **************************************************
  29. *opti echo 0 dime 3 elem cub8 ;
  30. opti dime 3 elem cub8 ;
  31. *
  32. * Maillage
  33. *
  34. l1 = 1. ;
  35. l2 = 1. ;
  36. ep0 = .01 ;
  37. p1 = 0. 0. 0. ;
  38. p2 = ep0 0. 0. ;
  39. p3 = ep0 l2 0. ;
  40. p4 = 0. l2 0. ;
  41. p10 = 0. 0. l1 ;
  42. p20 = ep0 0. l1 ;
  43. p30 = ep0 l2 l1 ;
  44. p40 = 0. l2 l1 ;
  45. li1 = d p1 1 p2 ;
  46. li2 = d p2 1 p3 ;
  47. li3 = d p3 1 p4 ;
  48. li4 = d p4 1 p1 ;
  49. li10 = d p10 1 p20 ;
  50. li20 = d p20 1 p30 ;
  51. li30 = d p30 1 p40 ;
  52. li40 = d p40 1 p10 ;
  53. sur1 = 'DALL' li1 li2 li3 li4 ;
  54. sur10 = 'DALL' li10 li20 li30 li40 ;
  55. mail1 = sur1 'VOLU' 1 sur10 ;
  56. *
  57. * Modele
  58. *
  59. youn0 = 2.1E11 ;
  60. nu0 = .3 ;
  61. mod0 = 'MODE' mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 'VISCOPLASTIQUE' 'GUIONNET' ;
  62. mat0 = MATE mod0 'YOUN' 192500. 'NU ' 0.3 'N ' 24 'K ' 200. 'CK ' 0. 'R0 ' 25. 'CD ' 0. 'M ' 0.8 'A ' 24000. 'A1 ' 3.3E5 'C ' 180. 'C1 ' 40. 'G ' 0. 'R ' 0. 'NN ' 0.005 'C2 ' 40. 'Q ' 1.466 'G1 ' 0. 'R1 ' 0. 'BETA' 1.1 'C0 ' 180. 'P1M0' 0.00333 ;
  63. * * Conditions aux limites *
  64. cl1 = 'BLOQ' sur1 'UZ' ;
  65. cl2 = 'BLOQ' ( li1 'ET' li10 ) 'UY' ;
  66. cl3 = 'BLOQ' ( li2 'ET' li20 ) 'UX' ;
  67. cl4 = 'BLOQ' sur10 'UZ' ;
  68. cl0 = cl1 'ET' cl2 'ET' cl3 'ET' cl4 ;
  69. *
  70. * Chargement
  71. *
  72. depmax = .1 ;
  73. dep0 = 'DEPI' cl4 depmax ;
  74. ev0 = 'EVOL' 'MANU' temps ( 'PROG' 0. 1000. ) y ( 'PROG' 0. 1. ) ;
  75. cha0 = 'CHAR' 'MECA' dep0 ev0 ;
  76. *
  77. dt0 = 10. ;
  78. tfin0 = 100. ;
  79. pro0 = 'PROG' 0. 'PAS' dt0 tfin0 ;
  80. *
  81. * Appel à PASAPAS
  82. *
  83. ta1 = 'TABLE' ;
  84. ta1 .'MODELE' = mod0 ;
  85. ta1 .'CARACTERISTIQUES' = mat0 ;
  86. ta1 .'CHARGEMENT' = cha0 ;
  87. ta1 .'BLOCAGES_MECANIQUES' = cl0 ;
  88. ta1 .'TEMPS_CALCULES' = pro0 ;
  89. *
  90. PASAPAS ta1 ;
  91. *
  92. * Post traitement
  93. *
  94. dim0 = 'DIME' ta1 .'TEMPS' ;
  95. *
  96. * Résultats calculés
  97. *
  98. reac0 = ta1 .'REACTIONS'.( dim0 - 1 ) ;
  99. depl0 = ta1 .'DEPLACEMENTS'.( dim0 - 1 ) ;
  100. temp0 = ta1 .'TEMPS'.( dim0 - 1 ) ;
  101. sigm0 = ta1 .'CONTRAINTES'.(dim0 - 1 ) ;
  102. eps0 = 'EPSI' mod0 depl0 ;
  103. *
  104. defo0 = 'DEFO' mail1 depl0 0. ;
  105. defo1 = 'DEFO' mail1 depl0 10. 'ROUG' ;
  106. *titr ' deformee finale' ;
  107. *trace ( defo0 'ET' defo1 ) ;
  108. *
  109. sigm0 = 'CHAN' 'CHPO' sigm0 mod0 ;
  110. *titr ' Contraintes suivant X' ;
  111. *trace ( 'EXCO' sigm0 'SMXX' ) mail1 ;
  112. *titr ' Contraintes suivant Y' ;
  113. *trace ( 'EXCO' sigm0 'SMYY' ) mail1 ;
  114. *titr ' Contraintes suivant Z' ;
  115. *trace ( 'EXCO' sigm0 'SMZZ' ) mail1 ;
  116. *titr ' Contraintes d cisaillement suivant YZ' ;
  117. *trace ( 'EXCO' sigm0 'SMYZ' ) mail1 ;
  118. reac0 = 'REDU' reac0 sur10 ;
  119. reac0 = 'EXCO' reac0 'FZ' 'SCAL' ;
  120. chp1 = 'MANU' 'CHPO' sur10 1 'SCAL' 1. ;
  121. reac1 = 'XTY' reac0 chp1 ( 'MOTS' 'SCAL' ) ( 'MOTS' 'SCAL' ) ;
  122. *
  123. *opti echo 0 dime 3 elem qua8 ;
  124. opti dime 3 elem qua8 ;
  125. *
  126. * Maillage
  127. *
  128. p1 = 0. 0. 0. ;
  129. p4 = 0. l2 0. ;
  130. p10 = 0. 0. l1 ;
  131. p40 = 0. l2 l1 ;
  132. li1 = d p1 1 p4 ;
  133. li2 = d p4 1 p40 ;
  134. li3 = d p40 1 p10 ;
  135. li4 = d p10 1 p1 ;
  136. mail1 = 'DALL' li1 li2 li3 li4 ;
  137. *titr ' Maillage' ;
  138. *trace mail1 ;
  139. *
  140. * Modele
  141. *
  142. mod0 = 'MODE' mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 'VISCOPLASTIQUE' 'GUIONNET' 'COQ8' ;
  143. mat0 = MATE mod0 'YOUN' 192500. 'NU ' 0.3 'N ' 24 'K ' 200. 'CK ' 0. 'A ' 24000. 'R0 ' 25. 'CD ' 0. 'M ' 0.8 'A1 ' 3.3E5 'C ' 180. 'C1 ' 40. 'G ' 0. 'R ' 0. 'NN ' 0.005 'C2 ' 40. 'Q ' 1.466 'G1 ' 0. 'R1 ' 0. 'BETA' 1.1 'C0 ' 180. 'P1M0' 0.00333 ;
  144. car0 = 'CARA' mod0 'EPAI' ep0 ;
  145. *
  146. * Conditions aux limites
  147. *
  148. cl1 = 'BLOQ' li1 'UZ' ;
  149. cl2 = 'BLOQ' li4 'UY' ;
  150. cl3 = 'BLOQ' mail1 'UX' ;
  151. cl4 = 'BLOQ' li3 'UZ' ;
  152. cl0 = cl1 'ET' cl2 'ET' cl3 'ET' cl4 ;
  153. *
  154. * Chargement
  155. *
  156. dep0 = 'DEPI' cl4 depmax ;
  157. ev0 = 'EVOL' 'MANU' temps ( 'PROG' 0. 1000. ) y ( 'PROG' 0. 1. ) ;
  158. cha0 = 'CHAR' 'DIMP' dep0 ev0 ;
  159. *
  160. * Appel à PASAPAS
  161. *
  162. ta1 = 'TABLE' ;
  163. ta1 .'MODELE' = mod0 ;
  164. ta1 .'CARACTERISTIQUES' = mat0 'ET' car0 ;
  165. ta1 .'CHARGEMENT' = cha0 ;
  166. ta1 .'BLOCAGES_MECANIQUES' = cl0 ;
  167. ta1 .'TEMPS_CALCULES' = pro0 ;
  168. *
  169. PASAPAS ta1 ;
  170. *
  171. * Post traitement
  172. *
  173. dim0 = 'DIME' ta1 .'TEMPS' ;
  174. *
  175. * Résultats calculés
  176. *
  177. reac0 = ta1 .'REACTIONS'.( dim0 - 1 ) ;
  178. depl0 = ta1 .'DEPLACEMENTS'.( dim0 - 1 ) ;
  179. temp0 = ta1 .'TEMPS'.( dim0 - 1 ) ;
  180. sigm0 = ta1 .'CONTRAINTES'.( dim0 - 1 ) ;
  181. eps0 = 'EPSI' mod0 depl0 car0 ;
  182. *
  183. defo0 = 'DEFO' mail1 depl0 0. ;
  184. defo1 = 'DEFO' mail1 depl0 10. 'ROUG' ;
  185. *titr ' deformee finale' ;
  186. *trace ( defo0 'ET' defo1 ) ;
  187. *
  188. sigm0 = 'CHAN' 'CHPO' sigm0 mod0 ;
  189. *titr ' Contraintes suivant Y' ;
  190. *trace ( 'EXCO' sigm0 'SMSS' ) mail1 ;
  191. *titr ' Contraintes suivant Z' ;
  192. *trace ( 'EXCO' sigm0 'SMTT' ) mail1 ;
  193. *titr ' Contraintes d cisaillement suivant YZ' ;
  194. *trace ( 'EXCO' sigm0 'SMST' ) mail1 ;
  195. reac0 = 'REDU' reac0 li3 ;
  196. reac0 = 'EXCO' reac0 'FZ' 'SCAL' ;
  197. chp1 = 'MANU' 'CHPO' li3 1 'SCAL' 1. ;
  198. reac2 = 'XTY' reac0 chp1 ( 'MOTS' 'SCAL' ) ( 'MOTS' 'SCAL' ) ;
  199. *
  200. 'MESS' ' Réaction au sommet de la structure , calcul massif =' reac1 ;
  201. 'MESS' ' Réaction au sommet de la structure , calcul coque épaisse =' reac2 ;
  202. 'SI' ( reac1 '>' 0. ) ;
  203. err1 = ( 'ABS' ( reac1 -reac2 ) ) / reac1 ;
  204. err1 = err1 * 100. ;
  205. 'SINON' ;
  206. err1 = 0. ;
  207. 'FINSI' ;
  208. 'MESS' ' Erreur obtenue en % =' err1 ;
  209. *
  210. 'SI' ( err1 '>' 0.03 ) ;
  211. 'ERRE' 5 ;
  212. 'SINON' ;
  213. 'ERRE' 0 ;
  214. 'FINSI' ;
  215. *
  216. 'FIN' ;
  217.  
  218.  
  219.  
  220.  
  221.  
  222.  
  223.  
  224.  
  225.  
  226.  
  227.  

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