Télécharger dragon.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : dragon.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. *-----------------------------------------------
  6. *
  7. * Modèle de comportement : DRAGON
  8. *
  9. * 'MODE' maillage 'ELASTIQUE' 'ISOTROPE'
  10. * 'PLASTIQUE_ENDOM' 'DRAGON' ;
  11. *
  12. * Validation en traction et en compression
  13. * sur un élément
  14. *
  15. *-----------------------------------------------
  16. *
  17. opti echo 0 dime 3 elem cub8 ;
  18. *
  19. *-----------------------------------------------
  20. * Comportement en traction
  21. *----------------------------------------------
  22. *
  23. l1 = 1. ;
  24. p1 = 0. 0. 0. ;
  25. p2 = l1 0. 0. ;
  26. p3 = l1 l1 0. ;
  27. p4 = 0. l1 0. ;
  28. d1 = d p1 1 p2 ;
  29. d2 = d p2 1 p3 ;
  30. d3 = d p3 1 p4 ;
  31. d4 = d p4 1 p1 ;
  32. sur1 = 'SURF' ( d1 'ET' d2 'ET' d3 'ET' d4 ) 'PLAN' ;
  33. p5 = 0. 0. l1 ;
  34. p6 = l1 0. l1 ;
  35. p7 = l1 l1 l1 ;
  36. p8 = 0. l1 l1 ;
  37. d5 = d p5 1 p6 ;
  38. d6 = d p6 1 p7 ;
  39. d7 = d p7 1 p8 ;
  40. d8 = d p8 1 p5 ;
  41. sur2 = 'SURF' ( d5 'ET' d6 'ET' d7 'ET' d8 ) 'PLAN' ;
  42. mail0 = 'VOLU' sur1 1 sur2 ;
  43. 'TITR' ' Maillage ' ;
  44. *'TRAC' mail0 ;
  45. *
  46. * Modèle
  47. *
  48. youn0 = .45E5 ;
  49. nu0 = .32 ;
  50. g0 = -0.10800E+03 ;
  51. alfa0 = 0.20000E+04 ;
  52. beta0 = -0.16513E+05 ;
  53. c00 = 0.00009E+01 ;
  54. c10 = 0.54000E+00 ;
  55. b0 = 0.00000E+00 ;
  56. mod0 = 'MODE' mail0 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  57. 'PLASTIQUE_ENDOM' 'DRAGON' ;
  58. mat0 = 'MATE' mod0 'YOUN' youn0 'NU' nu0 'g' g0 'ALFA' alfa0
  59. 'BETA' beta0 'C0' c00 'C1' c10 'B' b0 ;
  60. *
  61. * Temps du calcul
  62. *
  63. tfin0 = 100. ;
  64. dt0 = 10. ;
  65. *
  66. * Condition aux limites
  67. *
  68. cl1 = 'BLOQ' sur1 'UZ' ;
  69. cl2 = 'BLOQ' sur2 'UZ' ;
  70. cl3 = 'BLOQ' ( p1 'ET' p4 'ET' p5 'ET' p8 ) 'UX' ;
  71. cl4 = 'BLOQ' ( p1 'ET' p2 'ET' p5 'ET' p6 ) 'UY' ;
  72. cl0 = cl1 'ET' cl2 'ET' cl3 'ET' cl4 ;
  73. *
  74. * Chargement ( déplacement imposé )
  75. *
  76. depmax0 = 1.E-3 ;
  77. dep0 = 'DEPI' cl2 depmax0 ;
  78. ev0 = 'EVOL' 'MANU' temps ( 'PROG' 0. tfin0 ) y ( 'PROG' 0. 1. ) ;
  79. cha0 = 'CHAR' 'DIMP' dep0 ev0 ;
  80. *
  81. * Appel à PASAPAS
  82. *
  83. ta1 = 'TABLE' ;
  84. ta1 .'MODELE' = mod0 ;
  85. ta1 .'CARACTERISTIQUES' = mat0 ;
  86. ta1 .'BLOCAGES_MECANIQUES' = cl0 ;
  87. ta1 .'CHARGEMENT' = cha0 ;
  88. ta1 .'TEMPS_CALCULES' = 'PROG' 0. 'PAS' dt0 tfin0 ;
  89. *
  90. TMASAU=table;
  91. ta1 . 'MES_SAUVEGARDES'=TMASAU;
  92. ta1 . 'PRECISION'=1e-12;
  93. ta1 . 'PRECISINTER'=1e-12;
  94. TMASAU .'DEFTO'=VRAI;
  95. TMASAU .'DEFIN'=VRAI;
  96. PASAPAS ta1 ;
  97. *
  98. * Post traitement
  99. *
  100. dim0 = ( 'DIME' ( ta1 .'TEMPS' ) ) - 1 ;
  101. sigm0 = ta1 .'CONTRAINTES'.dim0 ;
  102. var0 = ta1. 'VARIABLES_INTERNES'.dim0 ;
  103. eps0 = ta1. 'DEFORMATIONS_INELASTIQUES'.dim0 ;
  104. depl0 = ta1 .'DEPLACEMENTS'.dim0 ;
  105. eps1 = 'EPSI' depl0 mod0 ;
  106. *
  107. * Erreur
  108. *
  109. err0 = 1.E-3 ;
  110. smxx1 = 'MAXI' ( 'EXCO' sigm0 'SMXX' ) ;
  111. smyy1 = 'MAXI' ( 'EXCO' sigm0 'SMYY' ) ;
  112. smzz1 = 'MAXI' ( 'EXCO' sigm0 'SMZZ' ) ;
  113. smxy1 = 'MAXI' ( 'EXCO' sigm0 'SMXY' ) ;
  114. smxz1 = 'MAXI' ( 'EXCO' sigm0 'SMXZ' ) ;
  115. smyz1 = 'MAXI' ( 'EXCO' sigm0 'SMYZ' ) ;
  116. *
  117. end0 = 'MAXI' ( 'EXCO' var0 'EPSE' ) ;
  118. end1 = 'MAXI' ( 'EXCO' var0 'END1' ) ;
  119. end2 = 'MAXI' ( 'EXCO' var0 'END2' ) ;
  120. end3 = 'MAXI' ( 'EXCO' var0 'END3' ) ;
  121. end4 = 'MAXI' ( 'EXCO' var0 'END4' ) ;
  122. end5 = 'MAXI' ( 'EXCO' var0 'END5' ) ;
  123. end6 = 'MAXI' ( 'EXCO' var0 'END6' ) ;
  124. *
  125. epxx1 = 'MAXI' ( 'EXCO' eps0 'EIXX' ) ;
  126. epyy1 = 'MAXI' ( 'EXCO' eps0 'EIYY' ) ;
  127. epzz1 = 'MAXI' ( 'EXCO' eps0 'EIZZ' ) ;
  128. epxy1 = 'MAXI' ( 'EXCO' eps0 'GIXY' ) ;
  129. epxz1 = 'MAXI' ( 'EXCO' eps0 'GIXZ' ) ;
  130. epyz1 = 'MAXI' ( 'EXCO' eps0 'GIYZ' ) ;
  131. *
  132. err1 = 'ABS' ( smxx1 - 0.0 ) ;
  133. 'SI' ( err1 '>' err0 ) ;
  134. 'ERRE' 5 ;
  135. 'FINSI' ;
  136. err1 = 'ABS' ( smyy1 - 0.0 ) ;
  137. 'SI' ( err1 '>' err0 ) ;
  138. 'ERRE' 5 ;
  139. 'FINSI' ;
  140. err1 = 'ABS' ( smzz1 - 20.875 ) ;
  141. 'SI' ( err1 '>' ( err0 * 20.875 ) ) ;
  142. 'ERRE' 5 ;
  143. 'FINSI' ;
  144. err1 = 'ABS' ( smxy1 - 0.0 ) ;
  145. 'SI' ( err1 '>' err0 ) ;
  146. 'ERRE' 5 ;
  147. 'FINSI' ;
  148. err1 = 'ABS' ( smxz1 - 0.0 ) ;
  149. 'SI' ( err1 '>' err0 ) ;
  150. 'ERRE' 5 ;
  151. 'FINSI' ;
  152. err1 = 'ABS' ( smyz1 - 0.0 ) ;
  153. 'SI' ( err1 '>' err0 ) ;
  154. 'ERRE' 5 ;
  155. 'FINSI' ;
  156. err1 = 'ABS' ( end0 - 0.1398 ) ;
  157. 'SI' ( err1 '>' err0 ) ;
  158. 'ERRE' 5 ;
  159. 'FINSI' ;
  160. err1 = 'ABS' ( end1 - 0.0 ) ;
  161. 'SI' ( err1 '>' err0 ) ;
  162. 'ERRE' 5 ;
  163. 'FINSI' ;
  164. err1 = 'ABS' ( end2 - 0.0 ) ;
  165. 'SI' ( err1 '>' err0 ) ;
  166. 'ERRE' 5 ;
  167. 'FINSI' ;
  168. err1 = 'ABS' ( end3 - 0.1398 ) ;
  169. 'SI' ( err1 '>' err0 ) ;
  170. 'ERRE' 5 ;
  171. 'FINSI' ;
  172. err1 = 'ABS' ( end4 - 0.0 ) ;
  173. 'SI' ( err1 '>' err0 ) ;
  174. 'ERRE' 5 ;
  175. 'FINSI' ;
  176. err1 = 'ABS' ( end5 - 0.0 ) ;
  177. 'SI' ( err1 '>' err0 ) ;
  178. 'ERRE' 5 ;
  179. 'FINSI' ;
  180. err1 = 'ABS' ( end6 - 0.0 ) ;
  181. 'SI' ( err1 '>' err0 )4/spaj>;
  182. 'ERRE' 5 ;
  183. 'FINSI' ;
  184. err1 = 'ABS' ( epxx1 + 3.23E-4 ) ;
  185. 'SI' ( err1 '>' err0 ) ;
  186. 'ERRE' 5 ;
  187. 'FINSI' ;
  188. err1 = 'ABS' ( epyy1 + 3.23E-4 ) ;
  189. 'SI' ( err1 '>' err0 ) ;
  190. 'ERRE' 5 ;
  191. 'FINSI' ;
  192. err1 = 'ABS' ( epzz1 - 1.E-3 ) ;
  193. 'SI' ( err1 '>' err0 ) ;
  194. 'ERRE' 5 ;
  195. 'FINSI' ;
  196. err1 = 'ABS' ( epxy1 - 0.0 ) ;
  197. 'SI' ( err1 '>' err0 ) ;
  198. 'ERRE' 5 ;
  199. 'FINSI' ;
  200. err1 = 'ABS' ( epxz1 - 0.0 ) ;
  201. 'SI' ( err1 '>' err0 ) ;
  202. 'ERRE' 5 ;
  203. 'FINSI' ;
  204. err1 = 'ABS' ( epyz1 - 0.0 ) ;
  205. 'SI' ( err1 '>' err0 ) ;
  206. 'ERRE' 5 ;
  207. 'FINSI' ;
  208. *
  209. *
  210. *-----------------------------------------------
  211. * Comportement en compression
  212. *-----------------------------------------------
  213. *
  214. * Chargement ( déplacement imposé )
  215. *
  216. depmax0 = -1.E-3 ;
  217. dep0 = 'DEPI' cl2 depmax0 ;
  218. ev0 = 'EVOL' 'MANU' temps ( 'PROG' 0. tfin0 ) y ( 'PROG' 0. 1. ) ;
  219. cha0 = 'CHAR' 'DIMP' dep0 ev0 ;
  220. *
  221. * Appel à PASAPAS
  222. *
  223. ta1 = 'TABLE' ;
  224. ta1 .'MODELE' = mod0 ;
  225. ta1 .'CARACTERISTIQUES' = mat0 ;
  226. ta1 .'BLOCAGES_MECANIQUES' = cl0 ;
  227. ta1 .'CHARGEMENT' = cha0 ;
  228. ta1 .'TEMPS_CALCULES' = 'PROG' 0. 'PAS' dt0 tfin0 ;
  229. *
  230. TMASAU=table;
  231. ta1 . 'MES_SAUVEGARDES'=TMASAU;
  232. ta1 . 'PRECISION'=1e-12;
  233. ta1 . 'PRECISINTER'=1e-12;
  234. TMASAU .'DEFTO'=VRAI;
  235. TMASAU .'DEFIN'=VRAI;
  236. PASAPAS ta1 ;
  237. *
  238. * Post traitement
  239. *
  240. dim0 = ( 'DIME' ( ta1 .'TEMPS' ) ) - 1 ;
  241. sigm0 = ta1 .'CONTRAINTES'.dim0 ;
  242. var0 = ta1. 'VARIABLES_INTERNES'.dim0 ;
  243. eps0 = ta1. 'DEFORMATIONS_INELASTIQUES'.dim0 ;
  244. depl0 = ta1 .'DEPLACEMENTS'.dim0 ;
  245. eps1 = 'EPSI' depl0 mod0 ;
  246. *
  247. * Erreur
  248. *
  249. err0 = 1.E-3 ;
  250. smxx1 = 'MAXI' ( 'EXCO' sigm0 'SMXX' ) ;
  251. smyy1 = 'MAXI' ( 'EXCO' sigm0 'SMYY' ) ;
  252. smzz1 = 'MAXI' ( 'EXCO' sigm0 'SMZZ' ) ;
  253. smxy1 = 'MAXI' ( 'EXCO' sigm0 'SMXY' ) ;
  254. smxz1 = 'MAXI' ( 'EXCO' sigm0 'SMXZ' ) ;
  255. smyz1 = 'MAXI' ( 'EXCO' sigm0 'SMYZ' ) ;
  256. *
  257. end0 = 'MAXI' ( 'EXCO' var0 'EPSE' ) ;
  258. end1 = 'MAXI' ( 'EXCO' var0 'END1' ) ;
  259. end2 = 'MAXI' ( 'EXCO' var0 'END2' ) ;
  260. end3 = 'MAXI' ( 'EXCO' var0 'END3' ) ;
  261. end4 = 'MAXI' ( 'EXCO' var0 'END4' ) ;
  262. end5 = 'MAXI' ( 'EXCO' var0 'END5' ) ;
  263. end6 = 'MAXI' ( 'EXCO' var0 'END6' ) ;
  264. *
  265. epxx1 = 'MAXI' ( 'EXCO' eps0 'EIXX' ) ;
  266. epyy1 = 'MAXI' ( 'EXCO' eps0 'EIYY' ) ;
  267. epzz1 = 'MAXI' ( 'EXCO' eps0 'EIZZ' ) ;
  268. epxy1 = 'MAXI' ( 'EXCO' eps0 'GIXY' ) ;
  269. epxz1 = 'MAXI' ( 'EXCO' eps0 'GIXZ' ) ;
  270. epyz1 = 'MAXI' ( 'EXCO' eps0 'GIYZ' ) ;
  271. *
  272. err1 = 'ABS' ( smxx1 + 0.0 ) ;
  273. 'SI' ( err1 '>' err0 ) ;
  274. 'ERRE' 5 ;
  275. 'FINSI' ;
  276. err1 = 'ABS' ( smyy1 + 0.0 ) ;
  277. 'SI' ( err1 '>' err0 ) ;
  278. 'ERRE' 5 ;
  279. 'FINSI' ;
  280. err1 = 'ABS' ( smzz1 + 4.19E1 ) ;
  281. 'SI' ( err1 '>' ( err0 * 41.9 ) ) ;
  282. 'ERRE' 5 ;
  283. 'FINSI' ;
  284. err1 = 'ABS' ( smxy1 - 0.0 ) ;
  285. 'SI' ( err1 '>' err0 ) ;
  286. 'ERRE' 5 ;
  287. 'FINSI' ;
  288. err1 = 'ABS' ( smxz1 - 0.0 ) ;
  289. 'SI' ( err1 '>' err0 ) ;
  290. 'ERRE' 5 ;
  291. 'FINSI' ;
  292. err1 = 'ABS' ( smyz1 - 0.0 ) ;
  293. 'SI' ( err1 '>' err0 ) ;
  294. 'ERRE' 5 ;
  295. 'FINSI' ;
  296. err1 = 'ABS' ( end0 - 7.23E-2 ) ;
  297. 'SI' ( err1 '>' err0 ) ;
  298. 'ERRE' 5 ;
  299. 'FINSI' ;
  300. err1 = 'ABS' ( end1 - 3.619E-2 ) ;
  301. 'SI' ( err1 '>' err0 ) ;
  302. 'ERRE' 5 ;
  • 'FINSI' ;
  • err1 = 'ABS' ( end2 - 3.619E-2 ) ;
  • 'SI' ( err1 '>' err0 ) ;
  • 'ERRE' 5 ;
  • 'FINSI' ;
  • err1 = 'ABS' ( end3 - 0.0 ) ;
  • 'SI' ( err1 '>' err0 ) ;
  • 'ERRE' 5 ;
  • 'FINSI' ;
  • err1 = 'ABS' ( end4 - 0.0 ) ;
  • 'SI' ( err1 '>' err0 ) ;
  • 'ERRE' 5 ;
  • 'FINSI' ;
  • err1 = 'ABS' ( end5 - 0.0 ) ;
  • 'SI' ( err1 '>' err0 ) ;
  • 'ERRE' 5 ;
  • 'FINSI' ;
  • err1 = 'ABS' ( end6 - 0.0 ) ;
  • 'SI' ( err1 '>' err0 ) ;
  • 'ERRE' 5 ;
  • 'FINSI' ;
  • err1 = 'ABS' ( epxx1 - 3.702E-4 ) ;
  • 'SI' ( err1 '>' err0 ) ;
  • 'ERRE' 5 ;
  • 'FINSI' ;
  • err1 = 'ABS' ( epyy1 - 3.702E-4 ) ;
  • 'SI' ( err1 '>' err0 ) ;
  • 'ERRE' 5 ;
  • 'FINSI' ;
  • err1 = 'ABS' ( epzz1 + 1.E-3 ) ;
  • 'SI' ( err1 '>' err0 ) ;
  • 'ERRE' 5 ;
  • 'FINSI' ;
  • err1 = 'ABS' ( epxy1 - 0.0 ) ;
  • 'SI' ( err1 '>' err0 ) ;
  • 'ERRE' 5 ;
  • 'FINSI' ;
  • err1 = 'ABS' ( epxz1 - 0.0 ) ;
  • 'SI' ( err1 '>' err0 ) ;
  • 'ERRE' 5 ;
  • 'FINSI' ;
  • err1 = 'ABS' ( epyz1 - 0.0 ) ;
  • 'SI' ( err1 '>' err0 ) ;
  • 'ERRE' 5 ;
  • 'FINSI' ;
  • *
  • 'FIN' ;
  •  
  •  
  •  
  •  
  •  
  •  
  •  
  •  
  •  
  •  
  •  
  • © Cast3M 2003 - Tous droits réservés.
    Mentions légales