Télécharger la plaquette Cast3M


fr en

Télécharger norton_tra2.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  norton_tra2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. **************************************************
  5. *                                                *
  6. *  TEST DE VALIDATION D'UNE LOI DE COMPORTEMENT  *
  7. *  --------------------------------------------  *
  8. *              DE MATERIAU                       *
  9. *              -----------                       *
  10. *                                                *
  11. *  LOI DE COMPORTEMENT DE FLUAGE DE :            *
  12. *            NORTON                              *
  13. *                                                *
  14. *  COMPARAISON DE CALCULS SUR DES ELEMENTS:      *
  15. *     - COQUE EPAISSE ( MFR=5 )                  *
  16. *     - COQUE MINCE   ( MFR=3 )                  *
  17. *     - MASSIF 3D     ( MFR=1 )                  *
  18. *                                                *
  19. *  STRUCTURE : UNE PLAQUE CARREE                 *
  20. *     - COTE : 1. M                              *
  21. *     - EPAISSEUR : 0.01  M                      *
  22. *                                                * 
  23. *  CHARGEMENT : FORCES IMPOSEES                  *
  24. *               ESSAI DE TRACTION                *
  25. *                                                *
  26. ************************************************** 
  27. * 
  28. opti echo 0 dime 3 elem cub8 ;  
  29. *  
  30. *  Maillage  
  31. *  
  32. l1 = 1. ; 
  33. l2 = 1. ; 
  34. ep0 = .01 ; 
  35. p1 = 0. 0. 0. ; 
  36. p2 = ep0 0. 0. ; 
  37. p3 = ep0 l2 0. ; 
  38. p4 = 0. l2 0. ; 
  39. p10 = 0. 0. l1 ;
  40. p20 = ep0 0. l1 ;
  41. p30 = ep0 l2 l1 ;
  42. p40 = 0. l2 l1 ; 
  43. li1 = d p1 1 p2 ; 
  44. li2 = d p2 1 p3 ; 
  45. li3 = d p3 1 p4 ;  
  46. li4 = d p4 1 p1 ; 
  47. li10 = d p10 1 p20 ;
  48. li20 = d p20 1 p30 ;
  49. li30 = d p30 1 p40 ; 
  50. li40 = d p40 1 p10 ;
  51. sur1 = 'DALL' li1 li2 li3 li4 ; 
  52. sur10 = 'DALL' li10 li20 li30 li40 ; 
  53. mail1 = sur1 'VOLU' 1 sur10 ; 
  54. *  
  55. *  Modele 
  56. *  
  57. AF1 = 1.E-45 ; 
  58. AF2 = 5. ; 
  59. AF3 = 1.2 ; 
  60. SMAX = 2.1E8 ; 
  61. youn0 = 2.1E11 ; 
  62. nu0 = .3 ; 
  63. mod0 = 'MODE' mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 
  64.        'FLUAGE' 'NORTON' ; 
  65. mat0 = 'MATE' mod0 'YOUN' youn0 'NU' nu0 'AF1' AF1 
  66.        'AF2' AF2 'AF3' AF3 'SMAX' SMAX ; 
  67. *  
  68. *  Conditions aux limites 
  69. *  
  70. cl1 = 'BLOQ' sur1 'UZ' ; 
  71. cl2 = 'BLOQ' ( li1 'ET' li10 ) 'UY' ; 
  72. cl3 = 'BLOQ' ( li4 'ET' li40 ) 'UX' ; 
  73. cl0 = cl1 'ET' cl2 'ET' cl3 ; 
  74. *  
  75. *  Chargement 
  76. *  
  77. depmax = 5.E5 ;  
  78. dep0 = 'FORCE' sur10 'FX' 0. 'FY' 0. 'FZ' depmax ;  
  79. ev0 = 'EVOL' 'MANU' temps ( 'PROG' 0. 1.E-2 1000. )
  80. y ( 'PROG' 0. 1. 1. ) ; 
  81. cha0 = 'CHAR' 'MECA' dep0 ev0 ; 
  82. *    
  83. dt0 = .0500 ;
  84. dt1 = 2. * dt0 ;
  85. dt2 = 2. * dt1 ;
  86. dt3 = 5. * dt2 ;
  87. tfin0 = 40. ; 
  88. pro0 = ( 'PROG' 0. 'PAS' 1.E-2 dt0 ) 'ET'
  89. ( 'PROG' ( 2. * dt0 ) 'PAS' 5. 50. ) ;
  90. *  
  91. *  Appel à PASAPAS 
  92. *  
  93. ta1 = 'TABLE' ; 
  94. ta1 .'MODELE' = mod0 ; 
  95. ta1 .'CARACTERISTIQUES' = mat0 ; 
  96. ta1 .'CHARGEMENT' = cha0 ; 
  97. ta1 .'BLOCAGES_MECANIQUES' = cl0 ;
  98. ta1 .'TEMPS_CALCULES' = pro0 ;  
  99. *  
  100. TMASAU=table;
  101. ta1 . 'MES_SAUVEGARDES'=TMASAU;
  102. TMASAU .'DEFTO'=VRAI;
  103. TMASAU .'DEFIN'=VRAI;
  104. PASAPAS ta1 ;  
  105. * 
  106. *  Post traitement  
  107. *  
  108. dim0 = 'DIME' ta1 .'TEMPS' ; 
  109. * 
  110. *  Résultats calculés  
  111. * 
  112. def0 = ta1 .'DEFORMATIONS_INELASTIQUES'.( dim0 - 1 ) ;  
  113. depl0 = ta1 .'DEPLACEMENTS'.( dim0 - 1 ) ; 
  114. sigm0 = ta1 .'CONTRAINTES'.( dim0 - 1 ) ;  
  115. temp0 = ta1 .'TEMPS'.( dim0 - 1 ) ; 
  116. var0 = ta1 .'VARIABLES_INTERNES'.( dim0 - 1 ) ; 
  117. eps0 = 'EPSI' mod0 depl0 ; 
  118. *  
  119. defo0 = 'DEFO' mail1 depl0 0. ;
  120. defo1 = 'DEFO' mail1 depl0 10. 'ROUG' ;
  121. *titr ' deformee finale' ;
  122. *trace ( defo0 'ET' defo1 ) ;
  123. *titr ' Deplacement suivant X' ;
  124. *trace ( 'EXCO' depl0 'UX' ) mail1 ;
  125. *titr ' Deplacement suivant Y' ;
  126. *trace ( 'EXCO' depl0 'UY' ) mail1 ;
  127. *titr ' Deplacement suivant Z' ;
  128. *trace ( 'EXCO' depl0 'UZ' ) mail1 ;
  129. sigm0 = 'CHAN' 'CHPO' sigm0 mod0 ;
  130. *titr ' Contrainte suivant X' ;
  131. *trace ( 'EXCO' sigm0 'SMXX' ) mail1 ;
  132. *titr ' Contrainte suivant Y' ;
  133. *trace ( 'EXCO' sigm0 'SMYY' ) mail1 ;
  134. *titr ' Contrainte suivant Z' ;
  135. *trace ( 'EXCO' sigm0 'SMZZ' ) mail1 ;
  136. *titr ' Cisaillement X Y' ;
  137. *trace ( 'EXCO' sigm0 'SMYZ' ) mail1 ;
  138. eps0 = 'CHAN' 'CHPO' eps0 mod0 ;
  139. *titr ' Deformation suivant X' ;
  140. *trace ( 'EXCO' eps0 'EPXX' ) mail1 ;
  141. *titr ' Deformation suivant Y' ;
  142. *trace ( 'EXCO' eps0 'EPYY' ) mail1 ;
  143. *titr ' Deformation suivant Z' ;
  144. *trace ( 'EXCO' eps0 'EPZZ' ) mail1 ;
  145. *titr ' Deformation X Y' ;
  146. *trace ( 'EXCO' eps0 'GAYZ' ) mail1 ;
  147. * 
  148. depl1 = 'EXTR' ( 'EXCO' eps0 'EPZZ' 'SCAL' ) 'SCAL' p10 ;
  149. *  
  150. opti echo 0 dime 3 elem tri3 ;  
  151. *  
  152. *  Maillage  
  153. *  
  154. l1 = 1. ; 
  155. l2 = 1. ; 
  156. p1 = 0. 0. 0. ; 
  157. p4 = 0. l2 0. ; 
  158. p10 = 0. 0. l1 ;
  159. p40 = 0. l2 l1 ; 
  160. li1 = d p1 1 p4 ; 
  161. li2 = d p4 1 p40 ; 
  162. li3 = d p40 1 p10 ;  
  163. li4 = d p10 1 p1 ; 
  164. mail1 = 'DALL' li1 li2 li3 li4 ; 
  165. *  
  166. *  Modele 
  167. *  
  168. mod0 = 'MODE' mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 
  169.        'FLUAGE' 'NORTON' 'DKT' ; 
  170. mat0 = 'MATE' mod0 'YOUN' youn0 'NU' nu0 'AF1' AF1 
  171.        'AF2' AF2 'AF3' AF3 'SMAX' SMAX ; 
  172. car0 = 'CARA' mod0 'EPAI' ep0 ; 
  173. *  
  174. *  Conditions aux limites 
  175. *  
  176. cl1 = 'BLOQ' li1 'UZ' ; 
  177. cl2 = 'BLOQ' ( p1 'ET' p10 ) 'UY' ; 
  178. cl3 = 'BLOQ' mail1 'UX' ; 
  179. cl0 = cl1 'ET' cl2 'ET' cl3 ; 
  180. *  
  181. *  Chargement 
  182. *  
  183. dep0 = 'FORCE' li3 'FX' 0. 'FY' 0. 'FZ' depmax ;  
  184. ev0 = 'EVOL' 'MANU' temps ( 'PROG' 0. 1.E-2 1000. )
  185. y ( 'PROG' 0. 1. 1. ) ; 
  186. cha0 = 'CHAR' 'MECA' dep0 ev0 ; 
  187. *    
  188. *  Appel à PASAPAS 
  189. *  
  190. ta1 = 'TABLE' ; 
  191. ta1 .'MODELE' = mod0 ; 
  192. ta1 .'CARACTERISTIQUES' = mat0 'ET' car0 ; 
  193. ta1 .'CHARGEMENT' = cha0 ; 
  194. ta1 .'BLOCAGES_MECANIQUES' = cl0 ;
  195. ta1 .'TEMPS_CALCULES' = pro0 ;  
  196. *
  197. TMASAU=table;
  198. ta1 . 'MES_SAUVEGARDES'=TMASAU;
  199. TMASAU .'DEFTO'=VRAI;
  200. TMASAU .'DEFIN'=VRAI;  
  201. PASAPAS ta1 ;  
  202. * 
  203. *  Post traitement  
  204. *  
  205. dim0 = 'DIME' ta1 .'TEMPS' ; 
  206. * 
  207. *  Résultats calculés  
  208. * 
  209. def0 = ta1 .'DEFORMATIONS_INELASTIQUES'.( dim0 - 1 ) ;  
  210. depl0 = ta1 .'DEPLACEMENTS'.( dim0 - 1 ) ; 
  211. sigm0 = ta1 .'CONTRAINTES'.( dim0 - 1 ) ;  
  212. temp0 = ta1 .'TEMPS'.( dim0 - 1 ) ; 
  213. var0 = ta1 .'VARIABLES_INTERNES'.( dim0 - 1 ) ; 
  214. eps0 = 'EPSI' mod0 depl0 car0 ; 
  215. sigm0 = 'RTENS' sigm0 mod0 ( 0 1 0 ) ;
  216. eps0 = 'RTENS' eps0 mod0 ( 0 1 0 ) ; 
  217. *  
  218. defo0 = 'DEFO' mail1 depl0 0. ;
  219. defo1 = 'DEFO' mail1 depl0 10. 'ROUG' ;
  220. *titr ' deformee finale' ;
  221. *trace ( defo0 'ET' defo1 ) ;
  222. *titr ' Deplacement suivant X' ;
  223. *trace ( 'EXCO' depl0 'UX' ) mail1 ;
  224. *titr ' Deplacement suivant Y' ;
  225. *trace ( 'EXCO' depl0 'UY' ) mail1 ;
  226. *titr ' Deplacement suivant Z' ;
  227. *trace ( 'EXCO' depl0 'UZ' ) mail1 ;
  228. sigm0 = 'CHAN' 'CHPO' sigm0 mod0 ;
  229. *titr ' Contrainte suivant Y' ;
  230. *trace ( 'EXCO' sigm0 'N11' ) mail1 ;
  231. *titr ' Contrainte suivant Z' ;
  232. *trace ( 'EXCO' sigm0 'N22' ) mail1 ;
  233. *titr ' Cisaillement Y Z' ;
  234. *trace ( 'EXCO' sigm0 'N12' ) mail1 ;
  235. eps0 = 'CHAN' 'CHPO' eps0 mod0 ;
  236. *titr ' Deformation suivant Y' ;
  237. *trace ( 'EXCO' eps0 'EPSS' ) mail1 ;
  238. *titr ' Deformation suivant Z' ;
  239. *trace ( 'EXCO' eps0 'EPTT' ) mail1 ;
  240. *titr ' Deformation Y Z' ;
  241. *trace ( 'EXCO' eps0 'GAST' ) mail1 ;
  242. *
  243. depl2 = 'EXTR' ( 'EXCO' eps0 'EPTT' 'SCAL' ) 'SCAL' p10 ;
  244. * 
  245. opti echo 0 dime 3 elem qua8 ;  
  246. *  
  247. *  Maillage  
  248. *  
  249. l1 = 1. ; 
  250. l2 = 1. ; 
  251. p1 = 0. 0. 0. ; 
  252. p4 = 0. l2 0. ; 
  253. p10 = 0. 0. l1 ;
  254. p40 = 0. l2 l1 ; 
  255. p50 = 0. ( l2 / 2. ) l1 ; 
  256. li1 = d p1 1 p4 ; 
  257. li2 = d p4 1 p40 ; 
  258. li3 = d p40 1 p10 ;  
  259. li4 = d p10 1 p1 ; 
  260. mail1 = 'DALL' li1 li2 li3 li4 ; 
  261. *  
  262. *  Modele 
  263. *  
  264. mod0 = 'MODE' mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 
  265.        'FLUAGE' 'NORTON' 'COQ8' ; 
  266. mat0 = 'MATE' mod0 'YOUN' youn0 'NU' nu0 'AF1' AF1 
  267.        'AF2' AF2 'AF3' AF3 'SMAX' SMAX ; 
  268. car0 = 'CARA' mod0 'EPAI' ep0 ; 
  269. *  
  270. *  Conditions aux limites 
  271. *  
  272. p51 = mail1 'POIN' 'PROC' p50 ; 
  273. cl1 = 'BLOQ' li1 'UZ' ; 
  274. cl2 = 'BLOQ' ( p1 'ET' p10 ) 'UY' ; 
  275. cl3 = 'BLOQ' mail1 'UX' ; 
  276. cl4 = 'RELA' 'ENSE' 'UZ' li3 ;  
  277. cl5 = 'RELA' 'ENSE' 'UY' li4 ; 
  278. cl6 = 'RELA' 'ENSE' 'UY' li2 ; 
  279. cl0 = cl1 'ET' cl2 'ET' cl3 'ET' cl4 
  280.       'ET' cl5 'ET' cl6 ; 
  281. *  
  282. *  Chargement 
  283. *  
  284. li5 = ( p10 'ET' p40 'ET' p51 ) ; 
  285. dep1 = 'FORCE' li5 'FX' 0. 'FY' 0. 'FZ' depmax ;  
  286. ev0 = 'EVOL' 'MANU' temps ( 'PROG' 0. 1.E-2 1000. )
  287. y ( 'PROG' 0. 1. 1. ) ; 
  288. cha0 = 'CHAR' 'MECA' dep1 ev0 ; 
  289. *    
  290. *  Appel à PASAPAS 
  291. *  
  292. ta1 = 'TABLE' ; 
  293. ta1 .'MODELE' = mod0 ; 
  294. ta1 .'CARACTERISTIQUES' = mat0 'ET' car0 ; 
  295. ta1 .'CHARGEMENT' = cha0 ; 
  296. ta1 .'BLOCAGES_MECANIQUES' = cl0 ;
  297. ta1 .'TEMPS_CALCULES' = pro0 ;  
  298. *
  299. TMASAU=table;
  300. ta1 . 'MES_SAUVEGARDES'=TMASAU;
  301. TMASAU .'DEFTO'=VRAI;
  302. TMASAU .'DEFIN'=VRAI;  
  303. PASAPAS ta1 ;  
  304. * 
  305. *  Post traitement  
  306. *  
  307. dim0 = 'DIME' ta1 .'TEMPS' ; 
  308. * 
  309. *  Résultats calculés  
  310. * 
  311. def0 = ta1 .'DEFORMATIONS_INELASTIQUES'.( dim0 - 1 ) ;  
  312. depl0 = ta1 .'DEPLACEMENTS'.( dim0 - 1 ) ; 
  313. sigm0 = ta1 .'CONTRAINTES'.( dim0 - 1 ) ;  
  314. temp0 = ta1 .'TEMPS'.( dim0 - 1 ) ; 
  315. var0 = ta1 .'VARIABLES_INTERNES'.( dim0 - 1 ) ; 
  316. eps0 = 'EPSI' mod0 depl0 car0 ; 
  317. sigm0 = 'RTENS' sigm0 mod0 ( 0 1 0 ) car0 ;
  318. eps0 = 'RTENS' eps0 mod0 ( 0 1 0 ) car0 ; 
  319. *  
  320. defo0 = 'DEFO' mail1 depl0 0. ;
  321. defo1 = 'DEFO' mail1 depl0 10. 'ROUG' ;
  322. *titr ' deformee finale' ;
  323. *trace ( defo0 'ET' defo1 ) ;
  324. *titr ' Deplacement suivant X' ;
  325. *trace ( 'EXCO' depl0 'UX' ) mail1 ;
  326. *titr ' Deplacement suivant Y' ;
  327. *trace ( 'EXCO' depl0 'UY' ) mail1 ;
  328. *titr ' Deplacement suivant Z' ;
  329. *trace ( 'EXCO' depl0 'UZ' ) mail1 ;
  330. sigm0 = 'CHAN' 'CHPO' sigm0 mod0 ;
  331. *titr ' Contrainte suivant Y' ;
  332. *trace ( 'EXCO' sigm0 'SMSS' ) mail1 ;
  333. *titr ' Contrainte suivant Z' ;
  334. *trace ( 'EXCO' sigm0 'SMTT' ) mail1 ;
  335. *titr ' Cisaillement Y Z' ;
  336. *trace ( 'EXCO' sigm0 'SMST' ) mail1 ;
  337. eps0 = 'CHAN' 'CHPO' eps0 mod0 ;
  338. *titr ' Deformation suivant Y' ;
  339. *trace ( 'EXCO' eps0 'EPSS' ) mail1 ;
  340. *titr ' Deformation suivant Z' ;
  341. *trace ( 'EXCO' eps0 'EPTT' ) mail1 ;
  342. *titr ' Deformation Y Z' ;
  343. *trace ( 'EXCO' eps0 'GAST' ) mail1 ;
  344. *
  345. depl3 = 'EXTR' ( 'EXCO' eps0 'EPTT' 'SCAL' ) 'SCAL' p10 ;
  346. * 
  347. 'MESS' ' Déformation suivant Z, calcul massif        =' depl1 ; 
  348. 'MESS' ' Déformation suivant Z, calcul coque mince   =' depl2 ; 
  349. 'MESS' ' Déformation suivant Z, calcul coque épaisse =' depl3 ; 
  350. *  
  351. 'SI' ( depl1 '>' 0. ) ; 
  352.   err1 = ( 'ABS' ( depl1 - depl2 ) ) / depl1 ; 
  353.   err1 = err1 * 100. ; 
  354.   err2 = ( 'ABS' ( depl1 - depl3 ) ) / depl1 ; 
  355.   err2 = err2 * 100. ; 
  356. 'SINON' ; 
  357.   err1 = 0. ; 
  358.   err2 = 0. ; 
  359. 'FINSI' ; 
  360. 'MESS' ' Erreurs relatives obtenues, en %' err1 err2 ; 
  361. *  
  362. 'SI' ( ( err1 '>' .1 ) 'OU' ( err2 '>' .3 ) ) ;
  363.   'ERRE' 5 ; 
  364. 'SINON' ; 
  365.   'ERRE' 0 ; 
  366. 'FINSI' ; 
  367. *  
  368. 'FIN' ;    
  369.  
  370.  
  371.  
  372.  
  373.  
  374.  
  375.  

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