Télécharger syco_3D_contpla.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : syco_3D_contpla.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. graph = faux ;
  5. *******************************************************
  6. **-------------------------------------------------------------------
  7. MESS '' ;
  8. MESS ' COMPARAISON 3D et 2D contraintes planes ';
  9. MESS ' pour lois SYCO1 et SYCO2 ';
  10. MESS ' CALCUL DYNAMIQUE COMPORTEMENT VISCOPLASTIQUE SYCO1 ';
  11. MESS ' en 3D ' ;
  12. MESS ' avec 1 face bloquee pour comparer a 2D cont plane' ;
  13. MESS ' Unites: N Pa m kg s' ;
  14. MESS '' ;
  15. **--------------------------------------------------------------------
  16. *******************************************************
  17. *** Options de calcul
  18. opti dime 3 elem CUB8 mode trid;
  19.  
  20. *******************************************************
  21. *** maillage
  22. pa1 = 0. 0. 0. ;
  23. pa2 = 0.0001 0. 0. ;
  24. pa3 = 0.0001 5. 0. ;
  25. pa4 = 0. 5. 0. ;
  26. li12 = pa1 droi 1 pa2;
  27. li23 = pa2 droi 1 pa3;
  28. li34 = pa3 droi 1 pa4;
  29. li41 = pa4 droi 1 pa1;
  30. su1 = dall li12 li23 li34 li41;
  31. vol1 = su1 volu tran 1 (0. 0. 3.);
  32. elim 0.00001 vol1;
  33. *
  34. pa5 = vol1 poin proch (0. 0. 3.);
  35. pa6 = vol1 poin proch (3. 0. 3.);
  36. pa7 = vol1 poin proch (3. 5. 3.);
  37. pa8 = vol1 poin proch (0. 5. 3.);
  38.  
  39. **********************************************************
  40. *** MODELE & MATERIAU ***
  41. **********************************************************
  42. *
  43. ep1 = prog 0. 0.003151475 0.072983512 0.147777886 ;
  44. si1 = (prog 0. 666. 851.6513994 951.65139) * 1.e6;
  45. lcxx = 'EVOL' 'MANU' 'Epsilon' ep1 'Sigma' si1 ;
  46. ab1 = extr lcxx 'ABSC';
  47. ord1 = extr lcxx 'ORDO';
  48. eey1 = (extr ord1 2)/(extr ab1 2);
  49. * Conversion courbe traction > ecrouissage
  50. lsm1 = si1 enle 1 ;
  51. lep1 = (ep1 enle 1) - (lsm1 / eey1) ;
  52. evec = evol vert manu eps lep1 sig lsm1 ;
  53. si graph ;
  54. dess (lcxx et evec) titr ' Courbes de traction et d ecrouissage (vert)' ;
  55. fins ;
  56. *
  57. nu0 = 0.3 ;
  58. rho0= 7800.;
  59. **********************************************************
  60. *** Modele & materiau ***
  61. *
  62. * elements standards
  63. mod1 = 'MODELISER' vol1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  64. 'VISCOPLASTIQUE' 'SYCO1' 'CUB8';
  65. mat1 = 'MATERIAU' mod1'YOUN' eey1 'NU' nu0 'RHO' rho0
  66. 'ECRO' evec 'PSYC' 12. 'DSYC' 1.d8 ;
  67. *
  68.  
  69. **********************************************************
  70. *** CL et DEPLACEMENTS IMPOSES ***
  71. **********************************************************
  72.  
  73. cl1 = bloq (pa1 et pa4 et pa8 et pa5) 'UX';
  74. cl2 = bloq (pa1 et pa2 et pa5 et pa6) 'UY';
  75. cl4 = bloq (pa1 et pa2 et pa3 et pa4) 'UZ';
  76. cl5 = bloq (pa3 et pa4 et pa7 et pa8) 'UY';
  77. cltot = (cl1 'ET' cl2 'ET' cl4 'ET' cl5) ;
  78. *
  79. *... chargement en deplacement impose
  80. umax = 1.;
  81. uimp = depi cl5 umax;
  82.  
  83. *
  84. xx1 = prog 0. 1.;
  85. yy1 = prog 0. 10.;
  86. evochar = evol manu 'temps' xx1 'charge' yy1;
  87. char1 = char meca uimp evochar;
  88. *
  89.  
  90. *=======================================================*
  91. * *
  92. * RESOLUTION *
  93. * *
  94. *=======================================================*
  95.  
  96. dt1 = 2.e-4;
  97. xx0 = prog 0. pas dt1 'NPAS' 30;
  98. *
  99. tab1 = table;
  100. tab1.dynamique = vrai;
  101. tab1.modele = mod1;
  102. tab1.caracteristiques = mat1;
  103. tab1.blocages_mecaniques = cltot;
  104. tab1.chargement = char1;
  105. tab1.temps_calcules = xx0;
  106. tab1.precision = 1.d-10;
  107. tab1.convergence_forcee = faux;
  108. pasapas tab1;
  109. *
  110. *post traitement
  111. ut1 = tab1.deplacements. 30;
  112. var1 = tab1.variables_internes. 30;
  113. *
  114. val1 = extr var1 'EPSE' 1 1 3;
  115. val2 = extr var1 'VP' 1 1 3;
  116. val3 = extr ut1 'UZ' pa8;
  117. **********************************************************
  118. **-------------------------------------------------------------------
  119. MESS '' ;
  120. MESS ' CALCUL DYNAMIQUE COMPORTEMENT VISCOPLASTIQUE SYCO1 ';
  121. MESS ' en 2D cont plane ' ;
  122. MESS ' Unites: N Pa m kg s' ;
  123. MESS '' ;
  124. **--------------------------------------------------------------------
  125. *** Options de calcul
  126. opti dime 2 elem QUA4 mode plan cont;
  127.  
  128. *******************************************************
  129. *** maillage
  130. pa1 = 0. 0.;
  131. pa2 = 5. 0.;
  132. pa3 = 5. 3.;
  133. pa4 = 0. 3.;
  134. li12 = pa1 droi 1 pa2;
  135. li23 = pa2 droi 1 pa3;
  136. li34 = pa3 droi 1 pa4;
  137. li41 = pa4 droi 1 pa1;
  138. su1 = dall li12 li23 li34 li41;
  139. elim 0.001 su1;
  140. *
  141.  
  142. **********************************************************
  143. *** MODELE & MATERIAU ***
  144. **********************************************************
  145. **********************************************************
  146. *** Modele & materiau ***
  147. *
  148. * elements standards
  149. mod1 = 'MODELISER' su1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  150. 'VISCOPLASTIQUE' 'SYCO1' 'QUA4';
  151. mat1 = 'MATERIAU' mod1 'YOUN' eey1 'NU' nu0 'RHO' rho0
  152. 'ECRO' evec 'PSYC' 12. 'DSYC' 1.d8 ;
  153. *
  154.  
  155. **********************************************************
  156. *** CL et DEPLACEMENTS IMPOSES ***
  157. **********************************************************
  158.  
  159. cl1 = bloq (pa1 et pa4) 'UX';
  160. cl2 = bloq (pa1 et pa2) 'UY';
  161. cl3 = bloq (pa2 et pa3) 'UX';
  162. cltot = (cl1 'ET' cl2 'ET' cl3 ) ;
  163. *
  164. *... chargement en deplacement impose
  165. umax = 1.;
  166. uimp = depi cl3 umax;
  167.  
  168. *
  169. xx1 = prog 0. 1.;
  170. yy1 = prog 0. 10.;
  171. evochar = evol manu 'temps' xx1 'charge' yy1;
  172. char1 = char meca uimp evochar;
  173. *
  174. *=======================================================*
  175. * *
  176. * RESOLUTION *
  177. * *
  178. *=======================================================*
  179.  
  180. dt1 = 2.e-4;
  181. *
  182. tab2 = table;
  183. tab2.dynamique = vrai;
  184. tab2.modele = mod1;
  185. tab2.caracteristiques = mat1;
  186. tab2.blocages_mecaniques = cltot;
  187. tab2.chargement = char1;
  188. tab2.temps_calcules = xx0;
  189. tab2.precision = 1.d-10;
  190. tab2.convergence_forcee = faux;
  191. pasapas tab2;
  192. *
  193. *post traitement
  194. ut1 = tab2.deplacements. 30;
  195. var1 = tab2.variables_internes. 30;
  196. *
  197. xal1 = extr var1 'EPSE' 1 1 3;
  198. xal2 = extr var1 'VP' 1 1 3;
  199. xal3 = extr ut1 'UY' pa3;
  200. *
  201. aa1 = ABS(val1 - xal1);
  202. aa2 = ABS(val2 - xal2);
  203. aa3 = ABS(val3 - xal3);
  204. *list aa1;
  205. *list aa2;
  206. *list aa3;
  207. *opti donn 5;
  208. *
  209. si (aa1 >EG 5.e-5);
  210. erre 5;
  211. finsi;
  212. si (aa2 >EG 5.e-3);
  213. erre 5;
  214. finsi;
  215. si (aa3 >EG 5.e-4);
  216. erre 5;
  217. finsi;
  218.  
  219. *******************************************************
  220. **-------------------------------------------------------------------
  221. MESS '' ;
  222. MESS ' CALCUL DYNAMIQUE COMPORTEMENT VISCOPLASTIQUE SYCO2 ';
  223. MESS ' en 3D ' ;
  224. MESS ' avec 1 face bloquee pour comparer a 2D cont plane' ;
  225. MESS ' Unites: N Pa m kg s' ;
  226. MESS '' ;
  227. **--------------------------------------------------------------------
  228. *******************************************************
  229. *** Options de calcul
  230. opti dime 3 elem CUB8 mode trid;
  231.  
  232. *******************************************************
  233. *** maillage
  234. pa1 = 0. 0. 0. ;
  235. pa2 = 0.01 0. 0. ;
  236. pa3 = 0.01 5. 0. ;
  237. pa4 = 0. 5. 0. ;
  238. li12 = pa1 droi 1 pa2;
  239. li23 = pa2 droi 1 pa3;
  240. li34 = pa3 droi 1 pa4;
  241. li41 = pa4 droi 1 pa1;
  242. su1 = dall li12 li23 li34 li41;
  243. vol1 = su1 volu tran 1 (0. 0. 3.);
  244. elim 0.001 vol1;
  245. *
  246. pa5 = vol1 poin proch (0. 0. 3.);
  247. pa6 = vol1 poin proch (3. 0. 3.);
  248. pa7 = vol1 poin proch (3. 5. 3.);
  249. pa8 = vol1 poin proch (0. 5. 3.);
  250.  
  251. **********************************************************
  252. *** MODELE & MATERIAU ***
  253. **********************************************************
  254. *
  255. *** Modele & materiau ***
  256. *
  257. TEMP01 = -125.;
  258. PSY9 = 4.25 - (0.015*TEMP01);
  259. ASY9 = 0.07592 - (0.07182*TEMP01/273.);
  260. BSY9 = 0.06352 - (0.23970*TEMP01/273.);
  261. CSY9 = 0.02528 - (0.04124*TEMP01/273.);
  262. *message 'P A B C' psy9 asy9 bsy9 csy9;
  263. * elements standards
  264. mod1 = 'MODELISER' vol1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  265. 'VISCOPLASTIQUE' 'SYCO2' 'CUB8';
  266. mat1 = 'MATERIAU' mod1'YOUN' eey1 'NU' nu0 'RHO' rho0
  267. 'ECRO' evec 'PSYC' psy9 'ASYC' asy9 'BSYC' bsy9 'CSYC' csy9;
  268. *
  269. **********************************************************
  270. *** CL et DEPLACEMENTS IMPOSES ***
  271. **********************************************************
  272.  
  273. cl1 = bloq (pa1 et pa4 et pa8 et pa5) 'UX';
  274. cl2 = bloq (pa1 et pa2 et pa5 et pa6) 'UY';
  275. cl4 = bloq (pa1 et pa2 et pa3 et pa4) 'UZ';
  276. cl5 = bloq (pa3 et pa4 et pa7 et pa8) 'UY';
  277. cltot = (cl1 'ET' cl2 'ET' cl4 'ET' cl5) ;
  278. *
  279. *... chargement en deplacement impose
  280. umax = 1.;
  281. uimp = depi cl5 umax;
  282.  
  283. *
  284. xx1 = prog 0. 1.;
  285. yy1 = prog 0. 10.;
  286. evochar = evol manu 'temps' xx1 'charge' yy1;
  287. char1 = char meca uimp evochar;
  288. *
  289. *
  290. *=======================================================*
  291. * *
  292. * RESOLUTION *
  293. * *
  294. *=======================================================*
  295.  
  296. dt1 = 2.e-4;
  297. *
  298. tab1 = table;
  299. tab1.dynamique = vrai;
  300. tab1.modele = mod1;
  301. tab1.caracteristiques = mat1;
  302. tab1.blocages_mecaniques = cltot;
  303. tab1.chargement = char1;
  304. tab1.temps_calcules = xx0;
  305. tab1.precision = 1.d-10;
  306. tab1.convergence_forcee = faux;
  307. pasapas tab1;
  308. *
  309. *post traitement
  310.  
  311. sig1 = tab1.contraintes. 30;
  312. ut1 = tab1.deplacements. 30;
  313. vt1 = tab1.vitesses. 30;
  314. tt1 = tab1.temps. 30;
  315. var1 = tab1.variables_internes. 30;
  316. *
  317. *
  318. val1 = extr var1 'EPSE' 1 1 3;
  319. val2 = extr var1 'VP' 1 1 3;
  320. val3 = extr ut1 'UZ' pa7;
  321.  
  322. *******************************************************
  323. **-------------------------------------------------------------------
  324. MESS '' ;
  325. MESS ' CALCUL DYNAMIQUE COMPORTEMENT VISCOPLASTIQUE SYCO2 ';
  326. MESS ' en 2D cont plane ' ;
  327. MESS ' Unites: N Pa m kg s' ;
  328. MESS '' ;
  329. **--------------------------------------------------------------------
  330. *** Options de calcul
  331. opti dime 2 elem QUA4 mode plan cont;
  332. opti 'TRAC' X ;
  333.  
  334. *******************************************************
  335. *** maillage
  336. pa1 = 0. 0.;
  337. pa2 = 5. 0.;
  338. pa3 = 5. 3.;
  339. pa4 = 0. 3.;
  340. li12 = pa1 droi 1 pa2;
  341. li23 = pa2 droi 1 pa3;
  342. li34 = pa3 droi 1 pa4;
  343. li41 = pa4 droi 1 pa1;
  344. su1 = dall li12 li23 li34 li41;
  345. elim 0.001 su1;
  346. *
  347.  
  348. **********************************************************
  349. *** MODELE & MATERIAU ***
  350. **********************************************************
  351. *** Modele & materiau ***
  352. *
  353. TEMP01 = -125.;
  354. PSY9 = 4.25 - (0.015*TEMP01);
  355. ASY9 = 0.07592 - (0.07182*TEMP01/273.);
  356. BSY9 = 0.06352 - (0.23970*TEMP01/273.);
  357. CSY9 = 0.02528 - (0.04124*TEMP01/273.);
  358. *message 'P A B C' psy9 asy9 bsy9 csy9;
  359. * elements standards
  360. mod1 = 'MODELISER' su1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  361. 'VISCOPLASTIQUE' 'SYCO2' 'QUA4';
  362. mat1 = 'MATERIAU' mod1'YOUN' eey1 'NU' nu0 'RHO' rho0
  363. 'ECRO' evec 'PSYC' psy9 'ASYC' asy9 'BSYC' bsy9 'CSYC' csy9;
  364. *
  365.  
  366. **********************************************************
  367. *** CL et DEPLACEMENTS IMPOSES ***
  368. **********************************************************
  369.  
  370. cl1 = bloq (pa1 et pa4) 'UX';
  371. cl2 = bloq (pa1 et pa2) 'UY';
  372. cl3 = bloq (pa2 et pa3) 'UX';
  373. cltot = (cl1 'ET' cl2 'ET' cl3 ) ;
  374. *
  375. *... chargement en deplacement impose
  376. umax = 1.;
  377. uimp = depi cl3 umax;
  378.  
  379. *
  380. xx1 = prog 0. 1.;
  381. yy1 = prog 0. 10.;
  382. evochar = evol manu 'temps' xx1 'charge' yy1;
  383. char1 = char meca uimp evochar;
  384. *
  385. *=======================================================*
  386. * *
  387. * RESOLUTION *
  388. * *
  389. *=======================================================*
  390.  
  391. dt1 = 2.e-4;
  392. *
  393. tab2 = table;
  394. tab2.dynamique = vrai;
  395. tab2.modele = mod1;
  396. tab2.caracteristiques = mat1;
  397. tab2.blocages_mecaniques = cltot;
  398. tab2.chargement = char1;
  399. tab2.temps_calcules = xx0;
  400. tab2.precison = 1.d-10;
  401. tab2.convergence_forcee = faux;
  402. pasapas tab2;
  403. *
  404. *post traitement
  405. sig1 = tab2.contraintes. 30;
  406. ut1 = tab2.deplacements. 30;
  407. vt1 = tab2.vitesses. 30;
  408. tt1 = tab2.temps. 30;
  409. var1 = tab2.variables_internes. 30;
  410. *
  411. *
  412. xal1 = extr var1 'EPSE' 1 1 3;
  413. xal2 = extr var1 'VP' 1 1 3;
  414. xal3 = extr ut1 'UY' pa3;
  415.  
  416. aa1 = ABS(val1 - xal1);
  417. aa2 = ABS(val2 - xal2);
  418. aa3 = ABS(val3 - xal3);
  419. *list aa1;
  420. *list aa2;
  421. *list aa3;
  422. *opti donn 5;
  423. *
  424. si (aa1 >EG 5.e-5);
  425. erre 5;
  426. finsi;
  427. si (aa2 >EG 5.e-3);
  428. erre 5;
  429. finsi;
  430. si (aa3 >EG 5.e-4);
  431. erre 5;
  432. finsi;
  433. fin;
  434.  
  435.  
  436.  
  437.  
  438.  
  439.  
  440.  
  441.  
  442.  

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