Télécharger syco_3D_defpla.dgibi

Retour à la liste

Numérotation des lignes :

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

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