Télécharger la plaquette Cast3M


fr en

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-9 ;
  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-9 ;
  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-9 ;
  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-9 ;
  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.  
  443.  
  444.  

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