Télécharger la plaquette Cast3M


fr en

Télécharger syco_3D_defpla.dgibi

Retour à la liste

Numérotation des lignes : 

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

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