Télécharger la plaquette Cast3M


fr en

Télécharger norton_cis2.dgibi

Retour à la liste

Numérotation des lignes : 

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

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