Télécharger syco_3D_contpla.dgibi

Retour à la liste

Numérotation des lignes :

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

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