Télécharger raff02.dgibi

Retour à la liste

Numérotation des lignes :

  1. ************************************************************
  2. ** raff02.dgibi
  3. ************************************************************
  4. * TEST de la procedure de raffinement en cour de calcul *
  5. * RAFF - PASAPAS - PROI en 3D avec plasticite *
  6. * un seul raffinements *
  7. * Calcul sur une CT comparaison avec l'essais 520RY-AT *
  8. * realise au DM2S/SEMT/LISN *
  9. * *
  10. *création : gg, le 15.03.2017 *
  11. ************************************************************
  12.  
  13.  
  14. ************************************************************
  15. * *
  16. * MAILLAGE CT25 *
  17. * Metal de base 707 + cordon de soudure *
  18. * *
  19. ************************************************************
  20. ************************************************************
  21.  
  22. TITRE 'EPROUVETTE CT25 3D' ;
  23. OPTI DIME 3 ELEM CUB8 ;
  24. *---------------------------------------------------
  25. *---------------------------------------------------
  26. * Option de trace
  27. *graph = vrai;
  28. graph = faux;
  29. *graph_ps = vrai;
  30. graph_ps = faux;
  31. *---------------------------------------------------
  32.  
  33. *---------------------------------------------------
  34. *---------------------------------------------------
  35. * DONNEES
  36. * unites = m N Pa
  37. *---------------------------------------------------
  38. *---------------------------------------------------
  39.  
  40. dn = 1.0E-7;
  41.  
  42. * Fissure
  43. *---------------------------------------------------
  44. * Taille de maille fissure ;
  45.  
  46. tm =5.0E-5;
  47. * longueur prefissure
  48.  
  49. lpfis = 0.02482;
  50.  
  51.  
  52. * CHOISIR longueur fissure maillage
  53. xfis10 = lpfis;
  54.  
  55.  
  56. * Eprouvettete
  57. *---------------------------------------------------
  58. a=0.02;
  59.  
  60. le = 0.004;
  61.  
  62. ep = 0.025 ;
  63.  
  64. *paramètres de raffinement
  65.  
  66. dmin=0.0030;
  67. densmax=0.008 ;
  68. densmin=densmax/8;
  69. *densmin=densmax/16.;
  70. dmax=dmin + (2*densmax) ;
  71.  
  72. *---------------------------------------------------
  73. * Nb elements
  74. *---------------------------------------------------
  75.  
  76. * epaisseur eprouvette
  77. ne_ep = 3 ;
  78. np_ep =ne_ep + 1;
  79. * maillage fissure
  80. nf1 = 59;
  81.  
  82. nf2 = 2;
  83.  
  84.  
  85.  
  86. * cercle goupille
  87. * nc1 multiple de 4
  88. nc1 = 32 ;
  89. *
  90. * nc2 multiple de 2
  91. nc2 = 2 ;
  92. ne1 = 2 ;
  93. ne2 = 2 ;
  94. ne3 = 4 ;
  95.  
  96.  
  97. *---------------------------------------------------
  98. * Eprouvette
  99. *---------------------------------------------------
  100.  
  101. * rayons
  102. r2 = 0.005 ;
  103.  
  104. r3 = 0.00625 ;
  105.  
  106. * angles
  107. a1 = 15. ;
  108. a2 = 30. ;
  109.  
  110. * ordonnees
  111. y1 = 0.0015 ;
  112. y2 = le/2. ;
  113. y3 = 0.006 ;
  114. y5 = 0.0155 ;
  115. y6 = 0.03 ;
  116. *
  117. yf1 = tm*3*3*2;
  118. list yf1;
  119.  
  120. * abscisses
  121. x1 = -0.0125 ;
  122. x2 = 0.00458 ;
  123. x3 = x2 + r2 ;
  124. x5 = a;
  125. x6 = 0.05;
  126. *
  127. xf1 = x6 - (nf1*3*3*tm);
  128. list xf1;
  129.  
  130. * coordonnees supplementaires
  131. y4 = y3 + ((-1.)*(x1*(tan(a1)))) ;
  132. x4 = x5 + ((-1.)*(y1/(tan(a2)))) ;
  133. x7 = x3 + (y3 - y2);
  134.  
  135.  
  136.  
  137.  
  138. *---------------------------------------------------
  139. *---------------------------------------------------
  140. * MAILLAGE
  141. *
  142. *---------------------------------------------------
  143. *---------------------------------------------------
  144.  
  145. *---------------------------------------------------
  146. * Les points
  147. *---------------------------------------------------
  148.  
  149. * centres des cercles
  150. O1 = 0. y5 0.;
  151. O2 = x2 0. 0.;
  152.  
  153. * contour eprouvette
  154. p1 = x6 0. 0.;
  155. p2 = x6 y6 0.;
  156. p3 = x1 y6 0.;
  157. p4 = x1 y4 0.;
  158. p5 = 0. y3 0.;
  159. p7 = 0. y2 0.;
  160. p8 = x3 y1 0.;
  161. p9 = x4 y1 0.;
  162. p10 = x5 0. 0.;
  163.  
  164. * sous decoupe eprouvette
  165. p11 = x7 y6 0.;
  166. p12 = (p11 COOR 1) (p5 COOR 2) 0.;
  167. p13 = (p11 COOR 1) (p9 COOR 2) 0.;
  168. p14 = x4 y6 0.;
  169. p15 = (p14 COOR 1) (p5 COOR 2) 0.;
  170. p16 = (p1 COOR 1) (p9 COOR 2) 0.;
  171. p17 = (p1 COOR 1) (p15 COOR 2) 0.;
  172. p18 = (r3*((2**(0.5))/2)) (y5+(r3*((2**(0.5))/2))) 0.;
  173. p19 = ((-1.)*(p18 COOR 1)) (p18 COOR 2) 0.;
  174.  
  175. p20 = xfis10 0. 0. ;
  176. p21 = xfis10 y1 0. ;
  177. p22 = xfis10 y3 0. ;
  178. p23 = xfis10 y6 0. ;
  179.  
  180. *--------------------------------------------------
  181. * Les lignes
  182. *---------------------------------------------------
  183. * eprouvette (sans fissure) ----------
  184. * se1
  185. p9p13 = p9 DROI ne2 p13;
  186. p15p12 = P15 DROI ne2 p12;
  187. * se2
  188. p14p11 = p14 DROI ne2 p11;
  189. * se3
  190. p13p12 = p13 DROI ne1 p12;
  191. p12p5 = p12 DROI (1+(nc1/8)) p5;
  192. p13p5 = p13p12 ET p12p5;
  193. ELIM dn p13p5 ; p13p5 = REGE p13p5;
  194. *p8p7 = CERC ((-1)*(ne1+1+(nc1/8))) p8 O2 p7 'DINI' 0.1 'DFIN' 2.;
  195. p8p7 = CERC ((-1)*(ne1+1+(nc1/8))) p8 O2 p7 'DINI' 0.0001 'DFIN' 0.002;
  196. p13p8 = p13 DROI nc2 p8;
  197. p7p5 = p7 DROI (nc2) p5;
  198. ELIM dn p7p5 ; p7p5 = REGE p7p5;
  199. *se4
  200. p18p19 = CERC (nc1/4) p18 O1 p19;
  201. p11p3 = p11 DROI (nc1/4) p3;
  202. c11 = CERC ((3*(nc1/4))/2) p19 O1 (O1 PLUS (0. ((-1)*r3) 0.));
  203. c12 = CERC ((3*(nc1/4))/2) (O1 PLUS (0. ((-1)*r3) 0.)) O1 p18;
  204. c1 = c11 ET c12;
  205. ELIM dn c1; c1 = REGE c1;
  206. p3p4 = p3 DROI (nc1/4) p4;
  207. p4p5 = p4 DROI ((-1)+(nc1/8)) p5;
  208. p5p12 = p5 DROI ((+1)+(nc1/8)) p12;
  209. p12p11 = p12 DROI (nc1/4) p11;
  210. c1bis = p3p4 ET p4p5 ET p5p12 ET p12p11;
  211.  
  212.  
  213. *--------------------------------------------------
  214. * Les surfaces
  215. *---------------------------------------------------
  216.  
  217. * eprouvette (sans fissure) *-----------------------
  218. se1 = p9p13 REGL ne1 p15p12;
  219. se2 = p15p12 REGL (nc1/4) p14p11;
  220. se3 = DALL (INVE p13p5) p13p8 p8p7 p7p5;
  221. se4a = p18p19 REGL 4 p11p3;
  222. se4b = c1 REGL 4 c1bis;
  223. se4 = se4a ET se4b;
  224. ELIM dn se4; se4 = REGE se4;
  225.  
  226. * eprouvette avec fissure *-----------------------
  227.  
  228. p15p22 = droi 3 p15 p22;
  229. p14p23 = droi 3 p14 p23;
  230. se5 = regl 8 p15p22 p14p23;
  231. elim dn se2 se5;
  232.  
  233. p22p17 = droi 9 p22 p17;
  234. p23p2 = droi 9 p23 p2;
  235. se6 = regl 8 p22p17 p23p2;
  236. elim dn se5 se6;
  237.  
  238. p9p21 = droi 3 p9 p21;
  239. se7 = regl 1 p9p21 p15p22;
  240. elim dn se1 se7;
  241. elim dn se6 se7;
  242.  
  243. p21p16 = droi 9 p21 p16;
  244. se8 = regl 1 p21p16 p22p17;
  245. elim dn se6 se8;
  246. elim dn se7 se8;
  247.  
  248. p10p20 = droi 3 p10 p20;
  249. se9 = regl 1 p10p20 p9p21;
  250. elim dn se9 se7;
  251.  
  252. p20p1 = droi 9 p20 p1;
  253. se10 = regl 1 p20p1 p21p16;
  254. elim dn se10 se9;
  255. elim dn se10 se8;
  256.  
  257. si (graph) ;
  258. trac (se1 et se2 et se3 et se4 et se5 et se6 et se7 et se8 et
  259. se9 et se10);
  260. finsi ;
  261. * eprouvette totale *--------------------------------
  262. seprou = se1 et se2 et se3 et se4 et se5 et se6 et se7 et se8 et
  263. se9 et se10;
  264. ELIM dn seprou;
  265. seprou = REGE seprou;
  266.  
  267.  
  268. * goupille *-----------------------------------------
  269. O1p18 = O1 DROI 1 p18;
  270. sgoup= O1p18 ROTA (nc1/4) (90.) O1 (O1 PLUS (0. 0. 1.));
  271. ELIM dn sgoup; sgoup = REGE sgoup;
  272. * ELIM dn (vgoupet seprou);
  273. *ELIM 0.001 (vgoupet seprou);
  274. ELIM 0.000001 (sgoup et seprou);
  275. *ELIM 0.001 seprou ;
  276. ELIM 0.000001 seprou ;
  277.  
  278. *--------------------------------------------------
  279. * Les volumes
  280. *---------------------------------------------------
  281. veprou= seprou VOLU ne_ep 'TRAN' (0. 0. (ep/2.));
  282. vcub8 = veprou;
  283.  
  284. vgoup = sgoup VOLU ne_ep 'TRAN' (0. 0. (ep/2.));
  285.  
  286. ELIM 0.000001 (vgoup et veprou );
  287.  
  288. *--------------------------------------------------
  289. * Boite pour la visualisation
  290. *--------------------------------------------------
  291. p24 = (1.2*xfis10) 0. 0. ;
  292. p9p24 = droi 1 p9 p24;
  293. sboit = p9p24 TRAN 1 (0. 0. (ep/2.));
  294.  
  295.  
  296. *----------------------------------------------------
  297. * Geometries pour CL en 3D et fissure
  298. *----------------------------------------------------
  299. xeprou yeprou zeprou = coor veprou ;
  300. pmil = yeprou poin 'INFERIEUR' (1.E-4*tm);
  301. lmil = (enve veprou) elem 'APPUYE' pmil ;
  302. fronfis = DROI ne_ep p20 (xfis10 0. (ep/2.));
  303. psym1 = xeprou poin 'SUPERIEUR' xfis10;
  304. pfis1 = xeprou poin 'INFERIEUR' xfis10;
  305. ssym1 = lmil elem 'APPUYE' 'LARGEMENT' psym1 ;
  306. sfis10 = lmil elem 'APPUYE' 'LARGEMENT' pfis1 ;
  307. surfz0= seprou;
  308.  
  309. xgoup ygoup zgoup = coor vgoup;
  310. lchar = ygoup poin 'INFERIEUR' (y5+(1.E-4*tm));
  311.  
  312. si (graph) ;
  313.  
  314. trac (vcub8 et (vgoup coul bleu) et (ssym1 coul roug)
  315. et (sfis10 coul vert)) 'CACH';
  316. OPTI TRAC PSC ;
  317. * OPTI FTRA 'proj_images/mail_n.ps';
  318. trac (vcub8 et (vgoup coul bleu) et (ssym1 coul roug)
  319. et (sfis10 coul vert)) 'CACH';
  320. finsi ;
  321.  
  322. *----------------------------------------------------
  323. * Modèle + Matériau
  324. *----------------------------------------------------
  325.  
  326. *Modèles éprouvette
  327. MODEPR = MODE VEPROU MECANIQUE ELASTIQUE PLASTIQUE EPSI LINEAIRE ;
  328. MODC8 = MODE VCUB8 MECANIQUE ELASTIQUE PLASTIQUE EPSI LINEAIRE ;
  329.  
  330. *Caractéristique materiaux
  331.  
  332. E1=2.09086E+11;
  333. nu1=0.3;
  334.  
  335. *Courbe de traction matériau éprouvette
  336. EPS1 = PROG 0. 2.50007E-3 5.39261E-3 1.40558E-2 2.28078E-2 3.17654E-2
  337. 4.1161E-2 5.14294E-2 6.33414E-2 7.82031E-2 9.81451E-2 0.12653 0.1685
  338. 0.23177 0.32756 0.47195 0.68754 1.0056 1.4685 2.1336 3.0765 4.3966
  339. 6.223 8.7212 12.102 16.633 22.647 30.56 40.883 54.247 71.419 93.330
  340. 121.10 ;
  341.  
  342. SG1 = PROG 0. 522.73 532.73 562.73 592.73 622.73 652.73 682.73 712.73
  343. 742.73 772.73 802.73 832.73 862.73 892.73 922.73 952.73 982.73 1012.7
  344. 1042.7 1072.7 1102.7 1132.7 1162.7 1192.7 1222.7 1252.7 1282.7 1312.7
  345. 1342.7 1372.7 1402.7 1432.7 ;
  346. SG1=SG1*1.0E+6;
  347. EVTRAC = EVOL 'MANU' EPS1 SG1 ;
  348.  
  349. si (graph);
  350. dess EVTRAC 'TITR' ' courbe traction eprouvette';
  351. finsi ;
  352.  
  353. *Matériaux éprouvette
  354. MATEPR = MATE MODEPR YOUN E1 NU nu1 TRAC EVTRAC ;
  355. MATC8 = MATE MODC8 YOUN E1 NU nu1 TRAC EVTRAC ;
  356.  
  357. *Modèle goupille
  358. MODGOUP = MODE VGOUP MECANIQUE ELASTIQUE EPSI LINEAIRE ;
  359.  
  360. *Matériau goupille
  361. MATGOUP = MATE MODGOUP YOUN E1 NU nu1 ;
  362.  
  363. MATold=(MATC8 et MATGOUP);
  364. MODold=(MODC8 et MODGOUP);
  365.  
  366. *----------------------------------------------------
  367. * CL + Chargement
  368. *----------------------------------------------------
  369.  
  370. KEPR = RIGI MODEPR MATEPR ;
  371. KRELA = EXTR KEPR 'RIGI' 'MULT';
  372. KRELAold=KRELA;
  373.  
  374. cl1 = bloq ssym1 uy ;
  375. cl2 = bloq lchar uy ;
  376. cl3 = bloq (surfz0 et sgoup) uz;
  377. cl4 = bloq lchar ux ;
  378. CL = CL1 ET CL2 et cl3 et CL4;
  379. CLodl=CL;
  380.  
  381. uchar=0.0016;
  382. f1 = depi cl2 uchar/2;
  383.  
  384.  
  385. LI1 = PROG 0. 1. ;
  386. LI2 = PROG 0. 1. ;
  387. EV1 = EVOL MANU T LI1 F(T) LI2 ;
  388.  
  389. CHA1 = CHAR 'DIMP' F1 EV1 ;
  390. list CHA1;
  391. *----------------------------------------------------
  392. * Résolution de tout le problème
  393. * sur le maillage grossier
  394. *----------------------------------------------------
  395. f4 = depi cl2 uchar;
  396. LI4 = PROG 0.0 1.0 ;
  397. EV4 = EVOL MANU T LI1 F(T) LI4 ;
  398. CHA4 = CHAR 'DIMP' F4 EV4 ;
  399.  
  400. TABPAP4=0.;
  401. TABPAP4= TABL;
  402. TABPAP4 . 'TEMPS_CALCULES' = PROG 0. pas 0.05 1. ;
  403. TABPAP4 . 'CARACTERISTIQUES' = (MATC8 et matgoup) ;
  404. TABPAP4 . 'CHARGEMENT' = CHA4 ;
  405. TABPAP4 . 'BLOCAGES_MECANIQUES' = (CL ET KRELA) ;
  406. TABPAP4 . 'MODELE' = (MODC8 et modgoup) ;
  407. TABPAP4 . 'LAGRANGIEN' = mot 'REACTUALISE' ;
  408. PASAPAS TABPAP4 ;
  409.  
  410. eglo4=PROG;
  411. I=0;
  412. repe bcl4 20 ;
  413. I=I+1;
  414.  
  415. sig1 = TABPAP4 . 'CONTRAINTES' . i;
  416. eps1 = EPSI (MODC8 et modgoup) TABPAP4 . 'DEPLACEMENTS'. i;
  417. eloc1 = ENER (MODC8 et modgoup) eps1 sig1;
  418. eglo1 = INTG eloc1 (MODC8 et modgoup) ;
  419. eglo4 = eglo4 et eglo1 ;
  420.  
  421.  
  422.  
  423. fin bcl4;
  424.  
  425. *DESS ev1;
  426. *----------------------------------------------------
  427. * Premier pas de résolution résolution
  428. *----------------------------------------------------
  429. temps=TABL;
  430.  
  431. *Résolution PASAPAS
  432. TABPAP = TABLE ;
  433. TABPAP . 'TEMPS_CALCULES' = PROG 0. pas 0.1 1. ;
  434. TABPAP . 'CARACTERISTIQUES' = (MATC8 et matgoup) ;
  435. TABPAP . 'CHARGEMENT' = CHA1 ;
  436. TABPAP . 'BLOCAGES_MECANIQUES' = CL ET KRELA ;
  437. TABPAP . 'MODELE' = (MODC8 et modgoup) ;
  438. TABPAP . 'LAGRANGIEN' = mot 'REACTUALISE' ;
  439. PASAPAS TABPAP ;
  440. TABRAF=TABL;
  441. TABRAF.1=TABPAP;
  442.  
  443.  
  444. ouv2 = PROG 0.0;
  445. for2 = PROG 0.0;
  446. eglo2= PROG 0.0;
  447. KI2 = PROG 0.0;
  448. Sig_m2 = PROG 0.0;
  449. KIamorc = 1213.;
  450.  
  451. repe bcl1 10 ;
  452.  
  453. ouv2 = ouv2 ET (2*(EXTR (TABPAP . 'DEPLACEMENTS' . &bcl1) 'UY' P7));
  454.  
  455. for1=0 ;
  456. REPE bcl11 4;
  457. for1 = for1 + ((EXTR (TABPAP. 'REACTIONS' . &bcl1 )'FY'
  458. (POIN lchar &bcl11))*2);
  459. fin bcl11;
  460. for2 = for2 ET for1;
  461. sig1 = TABPAP . 'CONTRAINTES' . &bcl1;
  462. eps1 = EPSI (MODC8 et modgoup) TABPAP . 'DEPLACEMENTS'. &bcl1;
  463. eloc1 = ENER (MODC8 et modgoup) eps1 sig1;
  464. eglo1 = INTG eloc1 (MODC8 et modgoup) ;
  465. eglo2 = eglo2 et eglo1 ;
  466.  
  467.  
  468. fin bcl1;
  469. *-----------------------------------------------------------------------
  470. * - raffiner le maillage
  471. * - receer le mod et le mat
  472. * - recreer les cl
  473. * - incrémenter le chargement
  474. * - projeter les variables issu
  475. * du dernier pas de pasapas
  476. * - réaliser la suite du calcul
  477. * pasapas
  478. *-----------------------------------------------------------------------
  479.  
  480.  
  481. *-----------------------------------------------------------------------
  482. * Raffinement autour de la pointe de fissure
  483. *-----------------------------------------------------------------------
  484.  
  485. xeprou yeprou zeprou= coor veprou ;
  486.  
  487. DF1= (((yeprou**2) + ((xfis10 -xeprou)**2.))**0.5);
  488. *DF1= (((yeprou**2) + ((BORN ((xfis10 - 0.0032)-xeprou)
  489. *'SCAL' 'MINIMUM' 0.0)**2.))**0.5);
  490.  
  491. aff = (((densmax - densmin)*(DF1 - dmax))/(dmax - dmin)) + densmax ;
  492. dens1 = born aff 'SCAL' 'COMPRIS' (densmin*1.000001)
  493. (densmax*1.000001);
  494.  
  495. si (graph) ;
  496. trac dens1 veprou ;
  497. finsi ;
  498.  
  499. veprou = raff veprou dens1 ;
  500. * LIST (NBEL SEPROU);
  501.  
  502.  
  503. VOLURL = ELEM VEPROU SURE ;
  504. VOLD= (VCUB8 et VGOUP);
  505. VOLDC8=VCUB8;
  506.  
  507. VCUB8 = VEPROU DIFF VOLURL ;
  508.  
  509. ELIM 0.000001 (Vgoup et Veprou);
  510. ELIM 0.000001 Veprou ;
  511. ELIM 0.000001 VCUB8 ;
  512. *----------------------------------------------------
  513. * Geometries pour CL en 3D et fissure
  514. *----------------------------------------------------
  515.  
  516. xeprou yeprou zeprou = coor veprou ;
  517. pmil = yeprou poin 'INFERIEUR' (1.E-4*tm);
  518. lmil = (enve veprou) elem 'APPUYE' pmil ;
  519. fronfis = DROI ne_ep (xfis10 0. 0.) (xfis10 0. (ep/2.));
  520. psym1 = xeprou poin 'SUPERIEUR' xfis10;
  521. pfis1 = xeprou poin 'INFERIEUR' xfis10;
  522. ssym1 = lmil elem 'APPUYE' 'LARGEMENT' psym1 ;
  523. sfis10 = lmil elem 'APPUYE' 'LARGEMENT' pfis1 ;
  524. poiz0 = zeprou poin 'INFERIEUR' (1.E-4*tm);
  525. surfz0 = (enve veprou) elem 'APPUYE' poiz0 ;
  526.  
  527. si (graph) ;
  528.  
  529. trac (vcub8 et (vgoup coul bleu) et (ssym1 coul roug)
  530. et (sfis10 coul vert)) 'CACH';
  531.  
  532. finsi ;
  533.  
  534. *----------------------------------------------------
  535. * Modèle + Matériau
  536. *----------------------------------------------------
  537.  
  538. MODOLD = MODC8;
  539. MODEPR = MODE VEPROU MECANIQUE ELASTIQUE PLASTIQUE EPSI LINEAIRE ;
  540. MODC8 = MODE VCUB8 MECANIQUE ELASTIQUE PLASTIQUE EPSI LINEAIRE ;
  541. MATEPR = MATE MODEPR YOUN E1 NU nu1 TRAC EVTRAC ;
  542. MATC8 = MATE MODC8 YOUN E1 NU nu1 TRAC EVTRAC ;
  543. *----------------------------------------------------
  544. * CL + Chargement
  545. *----------------------------------------------------
  546.  
  547. KEPR = RIGI MODEPR MATEPR ;
  548. KRELA = EXTR KEPR 'RIGI' 'MULT';
  549.  
  550. cl1 = bloq ssym1 uy ;
  551.  
  552. cl2 = bloq lchar uy ;
  553.  
  554. cl3 = bloq (surfz0 et sgoup) uz;
  555. cl4 = bloq lchar ux ;
  556. CL = CL1 ET CL2 et cl3 et CL4;
  557. CLodl=CL;
  558.  
  559. f1 = depi cl2 (uchar/2);
  560. LI2 = PROG (1.0) (2.0) ;
  561. EV1 = EVOL MANU T LI1 F(T) LI2 ;
  562. CHA1 = CHAR 'DIMP' F1 EV1 ;
  563. list CHA1;
  564.  
  565.  
  566. *-----------------------------------------------------------------------
  567. * Projection des champs issues de l'iteration précédente
  568. *-----------------------------------------------------------------------
  569.  
  570. * DEPLACEMENT
  571. uold=tabpap.continuation.deplacements;
  572. ugoup= redu uold vgoup;
  573. * passage de u en chamelem
  574. uproj1 = chan 'CHAM' (redu voldc8 uold) voldc8;
  575. * projection
  576. uproj2 = PROI uproj1 VCUB8;
  577. * passage de uproj en chpoin
  578. uproj3 = chan 'CHPO' uproj2 modC8;
  579.  
  580. uproj= uproj3 et ugoup ;
  581.  
  582. si (graph) ;
  583. trac uold vold TITR 'u old';
  584. trac uproj VCUB8 TITR 'u proj';
  585. finsi;
  586.  
  587. *CONTRAINTES
  588.  
  589. sigold=tabpap.continuation.contraintes;
  590. siggoup= redu sigold vgoup;
  591.  
  592. * passage point de gausses -> noeuds du chaml
  593. sigproj1 = CHAN 'NOEUD' MODOLD (redu voldc8 sigold) ;
  594. * projection
  595. sigproj2= PROI modC8 sigproj1 'STRESSES';
  596. * passage noeuds -> point de gausses du chaml
  597. * sigproj3 = CHAN 'STRESSES' MODC8 sigproj2 ;
  598. * sigproj4 = CHAN 'TYPE' sigproj3 'CONTRAINTES';
  599. sigproj4 = CHAN 'TYPE' sigproj2 'CONTRAINTES';
  600. sigproj = sigproj4 et siggoup ;
  601.  
  602.  
  603.  
  604. si (graph) ;
  605. trac sigold modold TITR 'sig old';
  606. trac sigproj modC8 TITR 'sig proj';
  607.  
  608. finsi;
  609.  
  610. *VARIABLES_INTERNES
  611.  
  612. pold=tabpap.continuation.variables_internes;
  613. pgoup= redu pold vgoup;
  614.  
  615. * passage point de gausses -> noeuds du chaml
  616. pproj1 = CHAN 'NOEUD' MODOLD (redu voldc8 pold) ;
  617. * projection
  618. pproj2= PROI modC8 pproj1 'STRESSES';
  619. pproj3= born pproj2 'EPSE' 'MINIMUM' 0.0;
  620. * passage noeuds -> point de gausses du chaml
  621. * pproj3 = CHAN 'STRESSES' MODC8 pproj2 ;
  622. pproj= pgoup et pproj3;
  623.  
  624. si (graph) ;
  625. trac pold modold TITR 'p old';
  626. trac pproj modC8 TITR 'p proj';
  627. finsi;
  628.  
  629. * calcul de l'énergie de déformation
  630.  
  631. epsold = EPSI (MODOLD et modgoup) uold;
  632. elocold = ENER (MODOLD et modgoup) epsold sigold;
  633. egloold = INTG elocold (MODOLD et modgoup) ;
  634.  
  635. epsproj = EPSI (modC8 et modgoup ) uproj;
  636. elocproj = ENER (modC8 et modgoup ) epsproj sigproj;
  637. egloproj = INTG elocproj (MODC8 et modgoup) ;
  638.  
  639. step1=(0.0 + 5.0E06)/22 ;
  640. liii1= PROG (0.0) PAS step1 5.0E06;
  641. step1=(0.0 + 3.50E08)/22 ;
  642. liii1= PROG (0.0) PAS step1 3.50E08;
  643.  
  644. si (graph) ;
  645. trac liii1 elocold modold TITR 'energie champ avant proj';
  646. trac liii1 elocproj modC8 TITR 'energie champprojete';
  647.  
  648. trac elocold modold TITR 'energie champ avant proj';
  649. trac elocproj modC8 TITR 'energie champprojete';
  650. finsi;
  651.  
  652. *Résolution PASAPAS
  653. TABPAP=0.;
  654. TABPAP= TABL;
  655. TABPAP . 'TEMPS_CALCULES' = PROG 0. pas 0.1 1. ;
  656. TABPAP . 'CARACTERISTIQUES' = (MATC8 et matgoup) ;
  657. TABPAP . 'CHARGEMENT' = CHA1 ;
  658. TABPAP . 'BLOCAGES_MECANIQUES' = CL ET KRELA ;
  659. TABPAP . 'MODELE' = (MODC8 et modgoup) ;
  660. TABPAP . 'LAGRANGIEN' = mot 'REACTUALISE' ;
  661. TABPAP . 'DEPLACEMENTS' = TABL;
  662. TABPAP . 'DEPLACEMENTS' . 0 = uproj;
  663. TABPAP . 'CONTRAINTES' = TABL;
  664. TABPAP . 'CONTRAINTES' . 0 = sigproj;
  665. TABPAP . 'VARIABLES_INTERNES' = TABL;
  666. TABPAP . 'VARIABLES_INTERNES' . 0 = pproj;
  667. PASAPAS TABPAP ;
  668. TABRAF.2=TABPAP;
  669.  
  670. * opti donn 5;
  671. J=0;
  672. REPE bcl2 10 ;
  673. J=J+1;
  674.  
  675. ouv1 =(2*(EXTR (TABPAP . 'DEPLACEMENTS' . j) 'UY' P7));
  676.  
  677.  
  678. ouv2 = ouv2 ET ouv1;
  679. K=0;
  680. for1=0;
  681. REPE bc21 np_ep;
  682. K=K+1;
  683. for1 = for1 + ((EXTR (TABPAP. 'REACTIONS' . j ) 'FY'
  684. (POIN lchar K))*2);
  685. fin bc21;
  686. for2 = for2 ET for1;
  687. sig1 = TABPAP . 'CONTRAINTES' . j;
  688. eps1 = EPSI (MODC8 et modgoup) TABPAP . 'DEPLACEMENTS'. j;
  689. eloc1 = ENER (MODC8 et modgoup) eps1 sig1;
  690. eglo1 = INTG eloc1 (MODC8 et modgoup) ;
  691. eglo2 = eglo2 et eglo1 ;
  692.  
  693.  
  694. fin bcl2;
  695.  
  696. *----------------------------------------------------------------------*
  697. * Resolution en partant directement du dernier maillage *
  698. *----------------------------------------------------------------------*
  699.  
  700. * chargement
  701. f1 = depi cl2 uchar;
  702. LI2 = PROG 0.0 1.0 ;
  703. EV1 = EVOL MANU T LI1 F(T) LI2 ;
  704. CHA1 = CHAR 'DIMP' F1 EV1 ;
  705.  
  706. TABPAP3=0.;
  707. TABPAP3= TABL;
  708. TABPAP3 . 'TEMPS_CALCULES' = PROG 0. pas 0.05 1. ;
  709. TABPAP3 . 'CARACTERISTIQUES' = (MATC8 et matgoup) ;
  710. TABPAP3 . 'CHARGEMENT' = CHA1 ;
  711. TABPAP3 . 'BLOCAGES_MECANIQUES' = CL ET KRELA ;
  712. TABPAP3 . 'MODELE' = (MODC8 et modgoup) ;
  713. TABPAP3 . 'LAGRANGIEN' = mot 'REACTUALISE' ;
  714. PASAPAS TABPAP3 ;
  715. *temps.(I+2)=DATE;
  716.  
  717. ouv3= prog;
  718. for3= prog;
  719. eglo3=prog;
  720. ouv4= prog;
  721. for4= prog;
  722.  
  723. J=0;
  724. REPE bcl3 20 ;
  725. J=J+1;
  726.  
  727. ouv1 =(2*(EXTR (TABPAP3 . 'DEPLACEMENTS' . j) 'UY' P7));
  728. ouv3 = ouv3 ET ouv1;
  729.  
  730. ouv4 = ouv4 et (2*(EXTR (TABPAP4 . 'DEPLACEMENTS' . j) 'UY' P7));
  731. K=0;
  732. for1=0;
  733. REPE bcl31 np_ep;
  734. K=K+1;
  735. for1 = for1 + ((EXTR (TABPAP3. 'REACTIONS' . j )'FY'
  736. (POIN lchar K))*2);
  737. fin bcl31;
  738. for3 = for3 ET for1;
  739. K=0;
  740. for1=0;
  741. REPE bcl32 np_ep;
  742. K=K+1;
  743. for1 = for1 + ((EXTR (TABPAP4. 'REACTIONS' . j )'FY'
  744. (POIN lchar K))*2);
  745. fin bcl32;
  746. for4 = for4 ET for1;
  747. sig1 = TABPAP3 . 'CONTRAINTES' . j;
  748. eps1 = EPSI (MODC8 et modgoup) TABPAP3 . 'DEPLACEMENTS'. j;
  749. eloc1 = ENER (MODC8 et modgoup) eps1 sig1;
  750. eglo1 = INTG eloc1 (MODC8 et modgoup) ;
  751. eglo3 = eglo3 et eglo1 ;
  752.  
  753. fin bcl3;
  754.  
  755. TAB1= TABL;
  756. TAB1 . 'TITRE'= 'TABLE';
  757. TAB1 . 'TITRE' . 1 = 'MOT' 'Raff en cour de calcul';
  758. TAB1 . 'TITRE' . 2 = 'MOT' 'Maillage fin' ;
  759. TAB1 . 'TITRE' . 3 = 'MOT' 'Maillage grossier' ;
  760. TAB1 . 'TITRE' . 4 = 'MOT' 'Courbe experimentale' ;
  761. *TAB1 . 1 = 'MOT' 'MARQ PLUS ';
  762. TAB1 . 2 = 'MOT' 'TIRR';
  763. TAB1 . 3 = 'MOT' 'TIRR';
  764.  
  765. * Données expérimentales pour comparaison :
  766. * (essais 5FTRA20RY-AT au LISN)
  767.  
  768.  
  769.  
  770. OUVEXP1 = PROG 0 2.41661E-03 5.30094E-03 7.63577E-03 1.00449E-02
  771. 1.27110E-02 1.56070E-02 1.86920E-02 2.16473E-02 2.45408E-02
  772. 2.77426E-02 3.07212E-02 3.34408E-02 3.67156E-02 4.03118E-02
  773. 4.37164E-02 4.62874E-02 5.01488E-02 5.31909E-02 5.61370E-02
  774. 5.99129E-02 6.34043E-02 6.65528E-02 6.99978E-02 7.32631E-02
  775. 7.70875E-02 8.08755E-02 8.46865E-02 8.85926E-02 9.18141E-02
  776. 9.52742E-02 9.88186E-02 0.10224 0.10657 0.10996 0.11396 0.11699 ;
  777. OUVEXP2 = PROG
  778. 0.12151 0.12498 0.12893 0.13283 0.13624 0.14064 0.14422 0.14830 0.15246
  779. 0.15596 0.16034 0.16425 0.16871 0.17227 0.17591 0.17957 0.18386
  780. 0.18799 0.19211 0.19610 0.20000 0.20378 0.20782 0.21256 0.21643
  781. 0.22014 0.22449 0.22856 0.23254 0.23670 0.24096 0.24470 0.24898
  782. 0.25300 0.25716 0.26109 0.26519 0.26925 0.27342 0.27816 0.28178
  783. 0.28576 0.28991 0.29410 0.29791 0.30235 0.30668 0.31095 0.31533
  784. 0.31948 0.32415 0.32819 0.33253;
  785. OUVEXP3 = PROG 0.33713
  786. 0.34150 0.34597 0.35026 0.35450 0.35883 0.36311 0.36743 0.37209 0.37653
  787. 0.38110 0.38537 0.38967 0.39411 0.39873 0.40309 0.40792 0.41242 0.41704
  788. 0.42137 0.42601 0.43042 0.43556 0.44009 0.44495 0.44955 0.45372 0.45806
  789. 0.46273 0.46761 0.47230 0.47712 0.48179 0.48671 0.49134 0.49650 0.50160
  790. 0.50623 0.51082 0.51634 0.52069 ;
  791. ouvexp = OUVEXP1 et OUVEXP2 ET OUVEXP3 ;
  792. * convertion en m
  793. ouvexp = 0.001*ouvexp;
  794.  
  795.  
  796. forexp1 = PROG
  797. 456.08621 852.54 1164.9 1548.8 1858.2 2286.3 2581.0
  798. 2930.3 3287.3 3573.5 3998.0 4388.6 4773.6 5174.7 5624.6 5985.0
  799. 6383.5 6811.0 7200.5 7646.5 8078.4 8530.0 8972.0 9405.0 9875.8
  800. 10339. 10802. 11255. 11715. 12197. 12659. 13113. 13599. 14104.
  801. 14587. 15079. 15533. 15990. 16466. 16991. 17503. 17984. 18501.;
  802. forexp2 = PROG
  803. 18979. 19474. 19970. 20461. 20989. 21498. 22031. 22501. 23011.
  804. 23529. 24052. 24538. 25023. 25511. 26017. 26507. 27024. 27557.
  805. 28015. 28540. 29037. 29512. 30038. 30523. 31008. 31491. 32003.
  806. 32556. 33043. 33557. 34057. 34541. 35032. 35524. 35980. 36502.
  807. 36983. 37509. 38014. 38534. 39051. 39536. 40020. 40547. 41007.
  808. 41510. 41980. 42478. 42946. 43439. 43923.;
  809. forexp3 = PROG
  810. 44400. 44871. 45352. 45827. 46339. 46847. 47294. 47771. 48213.
  811. 48679. 49142. 49596. 50044. 50490. 50935. 51389. 51841. 52291.
  812. 52750. 53187. 53645. 54089. 54507. 54902. 55338. 55790. 56257.
  813. 56719. 57154. 57562. 57987. 58401. 58821. 59262. 59687. 60060.
  814. 60441.;
  815. forexp= forexp1 et forexp2 et forexp3;
  816.  
  817.  
  818. EVFOUV2 = EVOL 'MANU' ouv2 for2 ;
  819. EVFOUV3 = EVOL 'ROUG' 'MANU' ouv3 for3 ;
  820. EVFOUV4 = EVOL 'BLEU' 'MANU' ouv4 for4 ;
  821. EVFOUVE = EVOL 'VERT' 'MANU' ouvexp forexp ;
  822.  
  823. EVEOUV2 = EVOL 'MANU' ouv2 eglo2 ;
  824. EVEOUV3 = EVOL 'ROUG' 'MANU' ouv3 eglo3 ;
  825. EVEOUV4 = EVOL 'BLEU' 'MANU' ouv4 eglo4 ;
  826.  
  827. res2 = TABRAF. 2 .'DEPLACEMENTS'. 10;
  828. ressig2 = TABRAF. 2 .'CONTRAINTES'. 10;
  829. resp2 = TABRAF. 2 .'VARIABLES_INTERNES'. 10;
  830. res3 = TABPAP3.'DEPLACEMENTS'. 20;
  831. ressig3 = TABPAP3 .'CONTRAINTES'. 20;
  832. resp3 = TABPAP3.'VARIABLES_INTERNES'. 20;
  833.  
  834. * deformees finales
  835.  
  836. def0 = DEFO 0.0 res2 veprou 'GRIS' ;
  837. def2 = DEFO 1 res2 veprou ;
  838. def3 = DEFO 1 res3 veprou 'ROUGE';
  839.  
  840.  
  841. si (graph);
  842.  
  843.  
  844. DESS (EVFOUV2 et EVFOUV3 et EVFOUV4 et EVFOUVE) 'TITR'
  845. 'Courbe Force Ouverture' 'LEGE' TAB1;
  846.  
  847. DESS (EVEOUV2 et EVEOUV3 et EVEOUV4) 'TITR'
  848. 'Courbe Energie totale Ouverture' 'LEGE' TAB1;
  849. TRAC (def0 et def2 et def3) 'CACH';
  850. finsi ;
  851.  
  852. * Champ différence
  853. udiffx = (EXCO res2 'UX') - (EXCO res3 'UX') ;
  854. udiffy = (EXCO res2 'UY') - (EXCO res3 'UY') ;
  855. udiffz = (EXCO res2 'UZ') - (EXCO res3 'UZ') ;
  856.  
  857. udiff = ( (udiffx**2)+(udiffy**2)+(udiffz**2))**0.5;
  858. umax = maxi ((((EXCO res3 'UX')**2)+
  859. ((EXCO res2 'UY')**2)+((EXCO res3 'UZ')**2))**0.5);
  860.  
  861. uerr = udiff /umax;
  862.  
  863. si (graph);
  864. trac udiff veprou TITR 'Ecart sur le deplacement';
  865. trac uerr veprou TITR 'Ecart relatif sur le deplacement';
  866. finsi;
  867.  
  868. perr = (resp2 - resp3) / (maxi resp3);
  869.  
  870.  
  871.  
  872. eps2 = EPSI (MODC8 et modgoup) res2;
  873. eloc2 = ENER (MODC8 et modgoup) eps2 ressig2;
  874. eps3 = EPSI (MODC8 et modgoup) res3;
  875. eloc3 = ENER (MODC8 et modgoup) eps3 ressig3;
  876.  
  877. elocerr = (ABS (eloc2 - eloc3)) / (maxi eloc3);
  878.  
  879. si (graph_ps);
  880. OPTI 'TRAC' 'PSC' 'EPTR' 7 ;
  881. OPTI 'POTR' 'HELVETICA_16';
  882.  
  883. * OPTI FTRA 'proj_images/fouv.ps';
  884. DESS (EVFOUV2 et EVFOUV3 et EVFOUV4 et EVFOUVE) 'TITR'
  885. 'Courbe Force Ouverture' 'LEGE' TAB1;
  886.  
  887. * OPTI FTRA 'proj_images/ufin.ps';
  888. trac (exco res2 'UY') veprou TITR 'Deplacement final';
  889.  
  890. * OPTI FTRA 'proj_images/ufinz.ps';
  891. trac (exco res2 'UY') veprou TITR 'Deplacement final' boit sboit;
  892.  
  893. * OPTI FTRA 'proj_images/sfin.ps';
  894. trac ressig2 modC8 TITR 'Contraintes finales';
  895.  
  896. * OPTI FTRA 'proj_images/sfinz.ps';
  897. trac ressig2 modC8 TITR 'Contraintes finales' boit sboit;
  898.  
  899. * OPTI FTRA 'proj_images/pfin.ps';
  900. trac resp2 modC8 TITR 'Deformation plastique finale';
  901.  
  902. * OPTI FTRA 'proj_images/pfinz.ps';
  903. trac resp2 modC8 TITR 'Deformation plastique finale' boit sboit;
  904.  
  905. * OPTI FTRA 'proj_images/eouv.ps';
  906. DESS (EVEOUV2 et EVEOUV3 et EVEOUV4) 'TITR'
  907. 'Courbe Energie totale Ouverture' 'LEGE' TAB1;
  908.  
  909. * OPTI 'TRAC' 'PSC' 'EPTR' 1 ;
  910. OPTI FTRA 'proj_images/udiff.ps';
  911. trac udiff veprou TITR 'Ecart sur le deplacement ';
  912.  
  913. * OPTI FTRA 'proj_images/uerr.ps';
  914. trac uerr veprou TITR 'Ecart relatif sur le deplacement';
  915.  
  916. * OPTI FTRA 'proj_images/perr.ps';
  917. trac perr modC8 TITR 'Ecart relatif sur la deformation
  918. plastique';
  919.  
  920. * OPTI FTRA 'proj_images/elocerr.ps';
  921. trac elocerr modC8 TITR 'Ecart relatif sur l energie de
  922. deformation locale';
  923.  
  924. * OPTI FTRA 'proj_images/uold.ps';
  925. trac (exco uold 'UY') vold TITR 'deplacement UY avent proj ';
  926.  
  927. * OPTI FTRA 'proj_images/uproj.ps';
  928. trac (exco uproj 'UY') VCUB8 TITR 'deplacement UY apres proj';
  929.  
  930. * OPTI FTRA 'proj_images/sigold.ps';
  931. trac sigold modold TITR 'sigma avent proj';
  932.  
  933. * OPTI FTRA 'proj_images/sigproj.ps';
  934. trac sigproj modC8 TITR 'sigma apres proj';
  935.  
  936. * OPTI FTRA 'proj_images/pold.ps';
  937. trac pold modold TITR 'variable interne avent proj';
  938.  
  939. * OPTI FTRA 'proj_images/pproj.ps';
  940. trac pproj modC8 TITR 'variable interne apres proj';
  941.  
  942. step1=(0.0 + 3.50E08)/22 ;
  943. liii1= PROG (0.0) PAS step1 3.50E08;
  944.  
  945. * OPTI FTRA 'proj_images/eold1.ps';
  946. trac liii1 elocold (modold )
  947. TITR 'energie champ avant proj';
  948.  
  949. * OPTI FTRA 'proj_images/eproj1.ps';
  950. trac liii1 elocproj (modC8 )TITR 'energie champ projete';
  951.  
  952. * OPTI FTRA 'proj_images/eold2.ps';
  953. trac elocold (modold ) TITR 'energie champ avant proj';
  954.  
  955. * OPTI FTRA 'proj_images/eproj2.ps';
  956. trac elocproj (modC8 )TITR 'energie champ projete';
  957.  
  958. OPTI 'TRAC' 'PSC' 'EPTR' 2 ;
  959.  
  960. * OPTI FTRA 'proj_images/mail_raff.ps';
  961. trac (vcub8 et (vgoup coul bleu) et (ssym1 coul roug)
  962. et (sfis10 coul vert)) 'CACH';
  963.  
  964. * OPTI FTRA 'proj_images/mail_raffz.ps';
  965. trac (vcub8 et (vgoup coul bleu) et (ssym1 coul roug)
  966. et (sfis10 coul vert)) 'CACH' boit sboit;
  967.  
  968.  
  969. finsi;
  970.  
  971.  
  972. *----------------------------------------------------------------------*
  973. * Ecart experience simulation sur le maillage raffiné *
  974. *----------------------------------------------------------------------*
  975. OPTI SAUV 'RAFF02.SAUV';
  976. nrjsim = extr (intg EVFOUV3 1 3) 1;
  977. ouvf = EXTR ouv3 3;
  978.  
  979. repe bcl4 (dime ouvexp);
  980. si ((EXTR ouvexp &bcl4) > ouvf) ;
  981. indf = &bcl4;
  982. quit bcl4;
  983. finsi;
  984. fin bcl4;
  985.  
  986.  
  987. nrjexp = extr (intg EVFOUVE 1 indf) 1;
  988. errexp = (ABS (nrjsim - nrjexp))/nrjexp;
  989.  
  990. *----------------------------------------------------------------------*
  991. * Erreur de projection *
  992. *----------------------------------------------------------------------*
  993.  
  994. nrjdir = extr (intg EVFOUV3 11 20) 1 ;
  995. nrjproj = extr (intg EVFOUV2 11 20) 1;
  996. errproj = (ABS (nrjdir - nrjproj))/nrjdir;
  997. list errexp;
  998. list errproj;
  999.  
  1000. SI ((errproj < (1.0E-2)) et (errexp < (1.0E-1)));
  1001. 'MESS' ' ----------------------' ;
  1002. 'MESS' ' SUCCES DU CAS-TEST !' ;
  1003. 'MESS' ' ----------------------' ;
  1004. Erre 0;
  1005. Sinon;
  1006. 'MESS' ' ---------------------' ;
  1007. 'MESS' ' ECHEC DU CAS-TEST !' ;
  1008. 'MESS' ' ---------------------' ;
  1009. Erre 5;
  1010. finsi;
  1011.  
  1012. FIN ;
  1013.  
  1014.  
  1015.  

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