Télécharger la plaquette Cast3M


fr en

Télécharger inclusion3d_thm.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : inclusion3d_thm.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. * test de la formulation du mdele inclusion3d
  6. * ----------------------------------------------------------------------
  7. *
  8. * Alain Sellier, Elsa Anglade, Aurelie Papon,
  9. * Clement Lacombe, Thierry Vidal
  10. * mercredi 4 janvier 2023
  11. *
  12. *-----------------------------------------------------------------------
  13. * Exemple de calcul d un element soumis a sechage suivi d un essai
  14. * de traction
  15. *-----------------------------------------------------------------------
  16. graph1=FAUX;
  17.  
  18. *maillage d'1/8 de cube (pour avoir les memes conditions limites que dans fluendo3d
  19. OPTI DIME 3 ELEM CUB8;
  20. *
  21. * Repere
  22. *
  23. * Z
  24. * | Y
  25. * |/
  26. * O----X
  27. *
  28. * Points du maillage-----------------------------------------------------------
  29. *
  30. *   8----7
  31. *  /|   /|
  32. * 5----6 |
  33. * | 4--|-3
  34. * |/   |/
  35. * 1----2
  36. *
  37. *Dimensions du cube (HAUTEUR,LARGEUR,PROFONDEUR)
  38. HAUT=1.;
  39. LARG=1.;
  40. PROF=1.;
  41. *
  42. *Nb d'éléments par arrete du cube (HAUTEUR,LARGEUR,PROFONDEUR)
  43. NELH=1;
  44. NELL=1;
  45. NELP=1;
  46. *
  47. *---Creation des points
  48. P1 = POIN 0. 0. 0.;
  49. P2 = POIN LARG 0. 0.;
  50. P3 = POIN LARG PROF 0.;
  51. P4 = POIN 0. PROF 0.;
  52. P5 = POIN 0. 0. HAUT;
  53. P6 = POIN LARG 0. HAUT;
  54. P7 = POIN LARG PROF HAUT;
  55. P8 = POIN 0. PROF HAUT;
  56. *---Creation de la face inferieure du maillage (SINF)
  57. D12=DROI NELL P1 P2;
  58. D23=DROI NELP P2 P3;
  59. D34=DROI NELL P3 P4;
  60. D41=DROI NELP P4 P1;
  61. SINF=DALL d12 d23 d34 d41;
  62. *---Creation du cube par translation
  63. CUB1 = (SINF VOLU NELH TRAN P5);
  64. ELIM (CUB1 ET SINF) 0.0001;
  65. *---Creation et denomination des droites sur Z
  66. D15=DROI NELH P1 P5;
  67. D26=DROI NELH P2 P6;
  68. D37=DROI NELH P3 P7;
  69. D48=DROI NELH P4 P8;
  70. ELIM (CUB1 ET D15 ET D26 ET D37 ET D48) 0.0001;
  71. *---Creation et denomination des droites sur la surface superieure
  72. D56=DROI NELH P5 P6;
  73. D67=DROI NELH P6 P7;
  74. D78=DROI NELH P7 P8;
  75. D85=DROI NELH P8 P5;
  76. *---Denomination des surfaces
  77. SAVA=D12 TRAN NELH P5 ;
  78. SGAU=D41 TRAN NELH P5 ;
  79. SSUP=(SINF PLUS P5) ;
  80. SARR=(SAVA PLUS P4) ;
  81. SDRO=(SGAU PLUS P2) ;
  82. ELIM (CUB1 ET SAVA ET SGAU ET SSUP ET SARR ET SDRO) 0.0001 ;
  83. ELIM (CUB1 ET D56 ET D67 ET D78 ET D85) 0.0001;
  84. *---Affichage du maillage
  85. si (graph1);
  86.  TRAC CUB1 'QUAL' ;
  87. finsi;
  88. *
  89. *
  90. *Paramètres du modèle--------------------------------------------------------
  91. *
  92. *
  93. *--Inclusion
  94. **général
  95. frac1=0.6;
  96. youn1=100000.;
  97. nu1=0.3;
  98. alpi1=5.0e-6;
  99. **retrait
  100. srwi1=1.;
  101. **plasticité
  102. rti1=10.;
  103. rri1=10.;
  104. **réaction chimique
  105. pchi1=0.;
  106. tchi1=2.;
  107. eaci1=46000.;
  108. wsci1=1;
  109. *
  110. *--Matrice
  111. **général
  112. frac0=1-frac1;
  113. youn0=35.;
  114. nu0=0.3;
  115. alpm0=1.0e-6;
  116. **retrait
  117. evol1=evol manu 'EAU' (prog 0. 1.) 'SWP0' (prog 0. 1.);
  118. srwm0=evol1;
  119. porm0=0.48;
  120. mvgm0=0.13;
  121. nvgm0=1.-(1./1.685);
  122. **plasticité
  123. rtm0=0.052;
  124. rrm0=0.052;
  125. **réaction chimique
  126. pchm0=0.;
  127. tchm0=2.;
  128. eacm0=46000.;
  129. wscm0=1;
  130. **fluage
  131. tflm0=0.;
  132. kflm0=0.;
  133. eafm0=17000.;
  134. *
  135. *
  136. *--Matériau homogène
  137. youh0=1./(frac1/youn1+frac0/youn0);
  138. nuh0=0.3;
  139. rho0=2400.;
  140. alph0=1.0e-5;
  141. trf0=20.;
  142. nin0=1.;
  143. *
  144. *Modele et Materiau-----------------------------------------------------
  145. mod1=mode CUB1 mecanique elastique viscoplastique inclusion3d; 
  146. *
  147. *--elasticity, plasticity and damage
  148. math1=MATE mod1 YOUN youh0 NU nuh0 RHO rho0 ALPH alph0 TTRF trf0 
  149. TREF 20. TALP 20. NINC nin0 RT 0. RR 0. DIM3 0.; 
  150.  
  151. *
  152. *--Matrice
  153. *
  154. matm1=MATE mod1 FRA0 frac0  YOU0 youn0  NUP0 nu0 ALP0 alpm0 
  155.   RTI0 rtm0 RTP0 rtm0 RFI0 rtm0 RFP0 rtm0 POR0 porm0 SWP0 srwm0 MVG0 mvgm0 NVG0 nvgm0 
  156.   FLM0 0.  FLK0 0.25 TFL0 0. EAF0 17000. VCH0 0.  ACS0 0. CPH0 (1./4.9)
  157.   SRS0 0.9  TCH0 1. EAC0 46000. COH0 10. DLT0 0.5 BTA0 0.2 KCH0 0.;  
  158.  
  159. *--Inclusion
  160. mati1=MATE mod1  'FRA1' frac1 'YOU1' youn1 'NUP1' nu1 'ALP1' alpi1  
  161.   'RTI1' rtm0 'RFI1' rtm0 'RTP1' rti1 'RFP1' rti1 'POR1' 0.001 'SWP1' srwi1 'MVG1' 0. 'NVG1' 0.5 'CPH1' 0.
  162.   'FLM1' 1.  'FLK1' 0.1 'TFL1' 0. 'EAF1' 17000. 'VCH1' 0.  'ACS1' 0. 
  163.   'SRS1' 0.9  'TCH1' 1. 'EAC1' 46000. 'COH1' 20. 'DLT1' 0.5 'BTA1' 0.2 KCH1 0.; 
  164. *
  165. *--assemblage
  166. mat1=math1 et matm1 ET mati1 ;
  167. *
  168. *
  169. *
  170. *Conditions aux limites-------------------------------------------------------------------------
  171. *
  172. *--Blocages
  173. cl1=bloq SGAU UX;
  174. cl2=bloq SAVA UY;
  175. cl3=bloq SINF UZ;
  176. *
  177. *--Déplacements égaux sur une surface
  178. SDRO2=SDRO MOIN P7;
  179. SARR2=SARR MOIN P7;
  180. SSUP2=SSUP MOIN P7;
  181. rig1=rela 1. 'UX' SDRO2 - 1. 'UX' P7;
  182. rig2=rela 1. 'UY' SARR2 - 1. 'UY' P7;
  183. rig3=rela 1. 'UZ' SSUP2 - 1. 'UZ' P7;
  184. *
  185. *--assemblage des conditions limites
  186. cl0=cl1 et cl2 et cl3 et rig1 et rig2 et rig3;
  187. *
  188. *
  189. *
  190. *Chargements-----------------------------------------------------------------------------------
  191. *
  192. *--Chargement mécanique
  193. FX1= -0.0001 0. 0.;
  194. FY1= 0. -0.0001 0.;
  195. FZ1= 0. 0. -0.0001;
  196. TFIN=1000.;
  197. * 
  198. FOR1=Forc FX1 SDRO;
  199. FOR2=Forc FY1 SARR;
  200. FOR3=forc FZ1 SSUP;
  201. * 
  202.  
  203. evolf1=evol manu 'TEMPS' (prog 0. TFIN) 'COEFF' (prog 0. 1.);
  204. evolf2=evol manu 'TEMPS' (prog 0. TFIN) 'COEFF' (prog 0. 1.);
  205. char1=char 'MECA' FOR1 evolf1;
  206. char2=char 'MECA' FOR2 evolf2;
  207. char3=char 'MECA' FOR3 evolf2;
  208. si (graph1) ;
  209.  dess evolf1 titr 'chargement mecanique';
  210. finsi ;
  211. *
  212. *vec1=(vect FOR1 1. 'FORC') coul roug;
  213. *vec2=(vect FOR2 1. 'FORC') coul bleu;
  214. *vec3=(vect FOR3 1. 'FORC') coul vert;
  215. *trac (vec1 et vec2  et vec3) (CUB1);
  216. *
  217. charm1=char1 et char2 et char3;
  218. *
  219. *
  220. *--Chargement thermique
  221. temp1=manu chpo CUB1 'T' 1.;
  222. ltth1=prog 0. TFIN;
  223. ltemp1=prog 20. 20.;
  224. evolth1=evol manu 'TEMPS' ltth1 'COEFF' ltemp1 coul roug;
  225. si (graph1) ;
  226.  dess evolth1 titr 'CHARGEMENT THERMIQUE';
  227. finsi ;
  228. chart1=char 'T' temp1 evolth1;
  229. *
  230. *
  231. *-- Chargement hydrique
  232. eau1=manu chpo CUB1 'EAU' 1.;
  233. *evolution  en fonction du temps 
  234. liste1=prog 0.   pas 0.5   TFIN;
  235. liste11=prog 0.  pas 0.5   TFIN;
  236. liste2=(liste11/180.)**5.;
  237. liste5=(liste11/160.)**2.;
  238. liste3=0.02+(0.98/(liste11/190.+liste2+liste5+ 1.));
  239. evole4=evol manu 'temps' liste1 'COEFF' liste3;
  240. si (graph1);
  241.  dess evole4 titr 'sechage';
  242. finsi;
  243. *Chargement hydrique sans retrait 
  244. *evolution  en fonction du temps 
  245. * liste1=prog 0.   pas 0.5 TFIN;
  246. * liste2=(liste1/180.)**5.;
  247. * liste5=(liste1/160.)**2.;
  248. * liste3=0.88+(0.12/(liste1/190.+liste2+liste5+ 1.));
  249. *evole4=evol manu 'temps' liste1 'COEFF' liste3;
  250.  
  251. *dess evole4 titr 'Evolution SRW matrice';
  252. charr1=char 'EAU' eau1 evole4;
  253. *
  254. *
  255. *
  256. *Définitions pour le calcul pas a pas du chargement initial--------------------------------
  257. tab1=table;
  258. *opti sauv 'cas_test_inclusion3d.sauv';
  259. *sauv CUB1;
  260. *sauv mod1;
  261. *sauv p7;
  262. ltc1=prog 0.  pas  1. 10. pas 5. 100. pas 10. 700. ;
  263. tab1.temps_calcules=ltc1;
  264. tab1.caracteristiques=mat1;
  265. tab1.modele=mod1;
  266. tab1.blocages_mecaniques=cl0;
  267. tab1.chargement=charm1 et chart1 et charr1;
  268. tab1.precision=1.0e-5;
  269. tab1.temps_sauves=ltc1;
  270. *tab1.temps_sauvegardes=ltc1;
  271. pasapas tab1;
  272. *
  273.  
  274.  
  275. *essai de traction 
  276. n1=dime tab1.temps;
  277. n1=n1-1;
  278. udroi=extr  tab1.deplacements.n1 'UX' p2;
  279. cl14=bloq SDRO UX;
  280. cl0=cl1 et cl2 et cl3 et  rig2 et rig3 et cl14;
  281. depi4=DEPI cl14 1.;
  282. evol4=evol manu 'TEMPS' (prog 0. 700.  TFIN) 
  283.                 'COEFF' (prog 0. udroi 0.1);
  284. *dess evol4;
  285. char1=char 'MECA' depi4 evol4;
  286. charm1=char1 et char2 et char3;
  287.  
  288. tab1.chargement=charm1 et chart1 et charr1;
  289. tab1.blocages_mecaniques=cl0;
  290. ltc1=prog 700.05  pas 0.05  730.  ;
  291. * ltc1=prog 705. pas 5.  TFIN  ;
  292. tab1.temps_calcules=ltc1;
  293. tab1.temps_sauves=ltc1;
  294. *tab1.temps_sauvegardes=ltc1;
  295. pasapas tab1;
  296.  
  297. ************************************************************************
  298. ************************************************************************
  299. *   exploitation des resultats
  300. ************************************************************************
  301. ************************************************************************
  302.  
  303. n1=dime (tab1.temps);
  304. i1=0;
  305. dp5=prog;
  306. e1s=prog;
  307. e2s=prog;
  308. e2d1=prog;
  309. e1d1=prog;
  310. lt5=prog;
  311. epx1=prog;
  312. epy1=prog;
  313. epz1=prog;
  314. six1=prog;
  315. siy1=prog;
  316. siz1=prog;
  317. eel00=prog;
  318. eelt0=prog;
  319. youn0=youn1;
  320. repeter bloc0 n1;
  321.  z5=extr tab1.deplacements.i1 'UX' p3;
  322.  vari1=tab1.variables_internes.i1;
  323.  depi1=tab1.deplacements.i1;
  324.  epsi1=epsi depi1 mod1;
  325.  sig1=tab1.contraintes.i1;
  326.  epxx1=extr ( exco epsi1 'EPXX') 'EPXX' 1 1 1;
  327. * dc1=extr (exco dc tab1.variables_internes.i1) dc 1 1 1 ;
  328.  dc1=0.;
  329.  youn1=youn0*(1.-dc1);
  330.  sxx1=extr ( exco sig1 'SMXX') 'SMXX' 1 1 1;
  331.  epyy1=extr(exco epsi1 'EPYY') 'EPYY' 1 1 1;
  332.  syy1=extr ( exco sig1 'SMYY') 'SMYY' 1 1 1;
  333.  epzz1=extr(exco epsi1 'EPZZ') 'EPZZ' 1 1 1;
  334.  szz1=extr ( exco sig1 'SMZZ') 'SMZZ' 1 1 1;
  335.  epxe1=0.;
  336.  epye1=0.;
  337.  epze1=0.;
  338.  dp5=dp5 et (prog z5);
  339.  epx1=epx1 et (prog (epxx1-epxe1));
  340.  six1=six1 et (prog sxx1);
  341.  epy1=epy1 et (prog (epyy1-epye1));
  342.  siy1=siy1 et (prog syy1);
  343.  epz1=epz1 et (prog (epzz1-epze1));
  344.  siz1=siz1 et (prog szz1);
  345.  lt5=lt5 et (prog (tab1.temps.i1));
  346.  i1=i1+1;
  347. fin bloc0;
  348.  
  349. edp5=evol manu 'TEMPS' lt5 'DP5' dp5;
  350. *points de controle
  351. n1=90; fref1= -7.95742E-03;
  352. n2=684; fref2= 2.73214E-03;
  353. LISTF1='EXTR' edp5 'ORDO';
  354. fcal1='EXTR' LISTF1 (n1);
  355. fcal2='EXTR' LISTF1 (n2);
  356. err1 = abs ((fref1 - fcal1)/fref1) ; 
  357. err2 = abs ((fref2 - fcal2)/fref2) ;
  358. si ((err1 < 1.e-3) et (err2 < 1.e-3)) ;
  359.     erre 0 ;
  360. sinon ;
  361.     erre 5 ;
  362. finsi ;
  363.  
  364.  
  365. si (graph1);
  366. dess edp5 titr 'retrait suivi de la traction';
  367. epx1=epx1;
  368. eepx1=(evol manu 'TEMPS' lt5 'EPXX' epx1) coul roug;
  369. esix1=evol manu 'TEMPS' lt5 'SMXX' six1;
  370. esiepx=evol manu 'EPXX' epx1 'SMXX' six1;
  371. *dess esiepx;
  372. eepx1=eepx1 coul roug;
  373. epy1=epy1;
  374. eepy1=evol manu 'TEMPS' lt5 'EPYY' epy1;
  375. esiy1=evol manu 'TEMPS' lt5 'SMYY' siy1;
  376. eepy1=eepy1 coul vert;
  377. epz1=epz1;
  378. eepz1=evol manu 'TEMPS' lt5 'EPZZ' epz1;
  379. esiz1=evol manu 'TEMPS' lt5 'SMZZ' siz1;
  380.  
  381. tab2=table;
  382. tab2.1 = 'MARQ CROI ' ;
  383. tab2.'TITRE' = 'TABLE' ;
  384. tab2.'TITRE'. 1 = MOT 'Dep P5' ;
  385. dess (eepx1 ) lege tab2;
  386. tab2.'TITRE'. 1 = MOT 'SMXX' ;
  387. dess (esix1 ) lege tab2;
  388. *opti donn 5;
  389.  
  390. *points experimentaux
  391. lt6= prog 0. 1. 4. 10. 20. 28. 29. 30. 34. 40. 50.;
  392. exx0=prog 0. 50. 130. 200. 250. 280. 220. 205. 200.
  393.  200. 200.;
  394. exx0=-1.*exx0;
  395. evx0=evol manu 'TEMPS' lt6 'EPXX exp' exx0;
  396. evx0=evx0 coul roug;
  397. tab3=table;
  398. tab3.1 = 'MARQ CROI NOLI' ;
  399. tab3.'TITRE' = 'TABLE' ;
  400. tab3.'TITRE'. 1 = MOT 'fluageXX exp' ;
  401. tab3.'TITRE'. 2 = MOT 'EPXX th' ;
  402. tab3.'TITRE'. 3 = MOT 'EPYY th' ;
  403. tab3.'TITRE'. 4=  MOT 'EPZZ th' ;
  404. tab3.'TITRE'. 5=  MOT 'SIXX th' ;
  405. tab3.'TITRE'. 6=  MOT 'SIYY th' ;
  406. tab3.'TITRE'. 7=  MOT 'SIZZ th' ;
  407. *dess (evx0 et eepx1 
  408. *et eepy1 et eepz1 et esix1 et esiy1 et esiz1) lege tab3;
  409. *opti donn 5;
  410.  
  411. *OPTI 'SORT' 'gopala_inclusion_uniax.xls' ;
  412. *SORT 'EXCE' eepx1 evx0 eepy1 eepz1;
  413. *OPTI 'SORT' 'EPXX_60S' ;
  414. *SORT 'EXCE' eepx1  'SEPA' 'VIRG';
  415.  
  416. *opti sort 'gopola_uniax_epxx_expe.res';
  417. *evolepx=evol manu lt6 exx0;
  418. *sort 'EXCE' evolepx ;
  419. **'SEPA' 'ESPA';
  420.  
  421.  
  422.  
  423. ************************************************************************
  424. *hydrique
  425.  
  426. *trace d un sclaire
  427. n1=dime tab1.temps;
  428. ecc1=table;
  429. tab2=table;
  430. tab2.'TITRE'=table;
  431. repe bou2 2;
  432.  lt1=prog;
  433.  ls1=prog;
  434.  nom1=chain 'EPH' (&bou2-1);
  435.  repe bou1 n1;
  436.    vari1=tab1.variables_internes.(&bou1-1);
  437.    temp1=tab1.temps.(&bou1-1);
  438.    ds1=extr (exco vari1 nom1) nom1 1 1 1;
  439.    lt1=lt1 et (prog temp1);
  440.   ls1=ls1 et (prog ds1);
  441.  fin bou1;
  442.  ecc1.&bou2=(evol manu 'temps' lt1 nom1 ls1);
  443.  tab2.'TITRE'. &bou2 = nom1 ;
  444. fin bou2;
  445. tab2.1 = 'MARQ CROI ' ;
  446. tab2.2 = 'MARQ PLUS ' ;
  447. tab2.3 = 'MARQ ETOI ' ;
  448. dess ( ecc1.1 et ecc1.2  ) lege tab2;
  449.  
  450.  
  451.  
  452. * opti donn 5;
  453. *trace d un sclaire
  454. n1=dime tab1.temps;
  455. ecc1=table;
  456. tab2=table;
  457. tab2.'TITRE'=table;
  458. repe bou2 2;
  459.   lt1=prog;
  460.   ls1=prog;
  461.   nom1=chain 'SEW' (&bou2-1);
  462.   repe bou1 n1;
  463.     vari1=tab1.variables_internes.(&bou1-1);
  464.     temp1=tab1.temps.(&bou1-1);
  465.     ds1=extr (exco vari1 nom1) nom1 1 1 1;
  466.     lt1=lt1 et (prog temp1);
  467.     ls1=ls1 et (prog ds1);
  468.   fin bou1;
  469.   ecc1.&bou2=(evol manu 'temps' lt1 nom1 ls1);
  470.   tab2.'TITRE'. &bou2 = nom1 ;
  471. fin bou2;
  472. tab2.1 = 'MARQ CROI ' ;
  473. tab2.2 = 'MARQ PLUS ' ;
  474. tab2.3 = 'MARQ ETOI ' ;
  475. dess ( ecc1.1 et ecc1.2  ) lege tab2;
  476.  
  477.  
  478. ************************************************************************
  479. *contraintes dans les phases
  480. ************************************************************************
  481.  
  482. n1=dime tab1.temps;
  483. ecc1=table;
  484. tab2=table;
  485. tab2.'TITRE'=table;
  486. repe bou2 6;
  487.   lt1=prog;
  488.   ls1=prog;
  489.   nom1=chain 'RS0' (&bou2);
  490.   repe bou1 n1;
  491.     vari1=tab1.variables_internes.(&bou1-1);
  492.     temp1=tab1.temps.(&bou1-1);
  493.     ds1=extr (exco vari1 nom1) nom1 1 1 1;
  494.     lt1=lt1 et (prog temp1);
  495.     ls1=ls1 et (prog ds1);
  496.   fin bou1;
  497.   ecc1.&bou2=(evol manu 'temps' lt1 nom1 ls1);
  498.   tab2.'TITRE'. &bou2 = nom1 ;
  499. fin bou2;
  500. tab2.1 = 'MARQ CROI ' ;
  501. tab2.2 = 'MARQ PLUS ' ;
  502. tab2.3 = 'MARQ ETOI ' ;
  503. tab2.4 = 'MARQ CARR ' ;
  504. tab2.5 = 'MARQ LOSA ' ;
  505. tab2.6 = 'MARQ ETOI NOLI' ;
  506. tab2.7 = 'MARQ CARR NOLI' ;
  507. dess ( ecc1.1 et ecc1.2 et ecc1.3 et ecc1.4 
  508.  et ecc1.5 et ecc1.6 ) lege tab2;
  509.  
  510.  
  511.  
  512. n1=dime tab1.temps;
  513. ecc1=table;
  514. tab2=table;
  515. tab2.'TITRE'=table;
  516. repe bou2 6;
  517.   lt1=prog;
  518.   ls1=prog;
  519.   nom1=chain 'RS1' (&bou2);
  520.   repe bou1 n1;
  521.     vari1=tab1.variables_internes.(&bou1-1);
  522.     temp1=tab1.temps.(&bou1-1);
  523.     ds1=extr (exco vari1 nom1) nom1 1 1 1;
  524.     lt1=lt1 et (prog temp1);
  525.     ls1=ls1 et (prog ds1);
  526.   fin bou1;
  527.   ecc1.&bou2=(evol manu 'temps' lt1 nom1 ls1);
  528.   tab2.'TITRE'. &bou2 = nom1 ;
  529. fin bou2;
  530. tab2.1 = 'MARQ CROI ' ;
  531. tab2.2 = 'MARQ PLUS ' ;
  532. tab2.3 = 'MARQ ETOI ' ;
  533. tab2.4 = 'MARQ CARR ' ;
  534. tab2.5 = 'MARQ LOSA ' ;
  535. tab2.6 = 'MARQ ETOI NOLI' ;
  536. tab2.7 = 'MARQ CARR NOLI' ;
  537. dess ( ecc1.1 et ecc1.2 et ecc1.3 et ecc1.4 
  538.  et ecc1.5 et ecc1.6  ) lege tab2;
  539.  
  540. n1=dime tab1.temps;
  541. ecc1=table;
  542. tab2=table;
  543. tab2.'TITRE'=table;
  544. repe bou2 6;
  545.   lt1=prog;
  546.   ls1=prog;
  547.   nom1=chain 'IS1' (&bou2);
  548.   repe bou1 n1;
  549.     vari1=tab1.variables_internes.(&bou1-1);
  550.     temp1=tab1.temps.(&bou1-1);
  551.     ds1=extr (exco vari1 nom1) nom1 1 1 1;
  552.     lt1=lt1 et (prog temp1);
  553.     ls1=ls1 et (prog ds1);
  554.   fin bou1;
  555.   ecc1.&bou2=(evol manu 'temps' lt1 nom1 ls1);
  556.   tab2.'TITRE'. &bou2 = nom1 ;
  557. fin bou2;
  558. tab2.1 = 'MARQ CROI ' ;
  559. tab2.2 = 'MARQ PLUS ' ;
  560. tab2.3 = 'MARQ ETOI ' ;
  561. tab2.4 = 'MARQ CARR ' ;
  562. tab2.5 = 'MARQ LOSA ' ;
  563. tab2.6 = 'MARQ ETOI NOLI' ;
  564. tab2.7 = 'MARQ CARR NOLI' ;
  565. dess ( ecc1.1 et ecc1.2 et ecc1.3 et ecc1.4 
  566.  et ecc1.5 et ecc1.6  ) lege tab2;
  567.  
  568. n1=dime tab1.temps;
  569. ecc1=table;
  570. tab2=table;
  571. tab2.'TITRE'=table;
  572. repe bou2 6;
  573.   lt1=prog;
  574.   ls1=prog;
  575.   nom1=chain 'OS1' (&bou2);
  576.   repe bou1 n1;
  577.     vari1=tab1.variables_internes.(&bou1-1);
  578.     temp1=tab1.temps.(&bou1-1);
  579.     ds1=extr (exco vari1 nom1) nom1 1 1 1;
  580.     lt1=lt1 et (prog temp1);
  581.     ls1=ls1 et (prog ds1);
  582.   fin bou1;
  583.   ecc1.&bou2=(evol manu 'temps' lt1 nom1 ls1);
  584.   tab2.'TITRE'. &bou2 = nom1 ;
  585. fin bou2;
  586. tab2.1 = 'MARQ CROI ' ;
  587. tab2.2 = 'MARQ PLUS ' ;
  588. tab2.3 = 'MARQ ETOI ' ;
  589. tab2.4 = 'MARQ CARR ' ;
  590. tab2.5 = 'MARQ LOSA ' ;
  591. tab2.6 = 'MARQ ETOI NOLI' ;
  592. tab2.7 = 'MARQ CARR NOLI' ;
  593. dess ( ecc1.1 et ecc1.2 et ecc1.3 et ecc1.4 
  594.  et ecc1.5 et ecc1.6 ) lege tab2;
  595.  
  596. ************************************************************************
  597. * deformations elastiques dans les phases
  598. ************************************************************************
  599.  
  600. n1=dime tab1.temps;
  601. ecc1=table;
  602. tab2=table;
  603. tab2.'TITRE'=table;
  604. repe bou2 6;
  605.   lt1=prog;
  606.   ls1=prog;
  607.   nom1=chain 'RE0' (&bou2);
  608.   repe bou1 n1;
  609.     vari1=tab1.variables_internes.(&bou1-1);
  610.     temp1=tab1.temps.(&bou1-1);
  611.     ds1=extr (exco vari1 nom1) nom1 1 1 1;
  612.     lt1=lt1 et (prog temp1);
  613.     ls1=ls1 et (prog ds1);
  614.   fin bou1;
  615.   ecc1.&bou2=(evol manu 'temps' lt1 nom1 ls1);
  616.   tab2.'TITRE'. &bou2 = nom1 ;
  617. fin bou2;
  618. tab2.1 = 'MARQ CROI ' ;
  619. tab2.2 = 'MARQ PLUS ' ;
  620. tab2.3 = 'MARQ ETOI ' ;
  621. tab2.4 = 'MARQ CARR ' ;
  622. tab2.5 = 'MARQ LOSA ' ;
  623. tab2.6 = 'MARQ ETOI NOLI' ;
  624. tab2.7 = 'MARQ CARR NOLI' ;
  625. dess ( ecc1.1 et ecc1.2 et ecc1.3 et ecc1.4 
  626.  et ecc1.5 et ecc1.6 ) lege tab2; 
  627.  
  628.  
  629.  
  630.  n1=dime tab1.temps;
  631. ecc1=table;
  632. tab2=table;
  633. tab2.'TITRE'=table;
  634. repe bou2 6;
  635.   lt1=prog;
  636.   ls1=prog;
  637.   nom1=chain 'RE1' (&bou2);
  638.   repe bou1 n1;
  639.     vari1=tab1.variables_internes.(&bou1-1);
  640.     temp1=tab1.temps.(&bou1-1);
  641.     ds1=extr (exco vari1 nom1) nom1 1 1 1;
  642.     lt1=lt1 et (prog temp1);
  643.     ls1=ls1 et (prog ds1);
  644.   fin bou1;
  645.   ecc1.&bou2=(evol manu 'temps' lt1 nom1 ls1);
  646.   tab2.'TITRE'. &bou2 = nom1 ;
  647. fin bou2;
  648. tab2.1 = 'MARQ CROI ' ;
  649. tab2.2 = 'MARQ PLUS ' ;
  650. tab2.3 = 'MARQ ETOI ' ;
  651. tab2.4 = 'MARQ CARR ' ;
  652. tab2.5 = 'MARQ LOSA ' ;
  653. tab2.6 = 'MARQ ETOI NOLI' ;
  654. tab2.7 = 'MARQ CARR NOLI' ;
  655. dess ( ecc1.1 et ecc1.2 et ecc1.3 et ecc1.4 
  656.  et ecc1.5 et ecc1.6  ) lege tab2;
  657.  
  658.  
  659. n1=dime tab1.temps;
  660. ecc1=table;
  661. tab2=table;
  662. tab2.'TITRE'=table;
  663. repe bou2 6;
  664.   lt1=prog;
  665.   ls1=prog;
  666.   nom1=chain 'IE1' (&bou2);
  667.   repe bou1 n1;
  668.     vari1=tab1.variables_internes.(&bou1-1);
  669.     temp1=tab1.temps.(&bou1-1);
  670.     ds1=extr (exco vari1 nom1) nom1 1 1 1;
  671.     lt1=lt1 et (prog temp1);
  672.     ls1=ls1 et (prog ds1);
  673.   fin bou1;
  674.   ecc1.&bou2=(evol manu 'temps' lt1 nom1 ls1);
  675.   tab2.'TITRE'. &bou2 = nom1 ;
  676. fin bou2;
  677. tab2.1 = 'MARQ CROI ' ;
  678. tab2.2 = 'MARQ PLUS ' ;
  679. tab2.3 = 'MARQ ETOI ' ;
  680. tab2.4 = 'MARQ CARR ' ;
  681. tab2.5 = 'MARQ LOSA ' ;
  682. tab2.6 = 'MARQ ETOI NOLI' ;
  683. tab2.7 = 'MARQ CARR NOLI' ;
  684. dess ( ecc1.1 et ecc1.2 et ecc1.3 et ecc1.4 
  685.  et ecc1.5 et ecc1.6  ) lege tab2;
  686.  
  687.  
  688.  
  689. n1=dime tab1.temps;
  690. ecc1=table;
  691. tab2=table;
  692. tab2.'TITRE'=table;
  693. repe bou2 6;
  694.   lt1=prog;
  695.   ls1=prog;
  696.   nom1=chain 'OE1' (&bou2);
  697.   repe bou1 n1;
  698.     vari1=tab1.variables_internes.(&bou1-1);
  699.     temp1=tab1.temps.(&bou1-1);
  700.     ds1=extr (exco vari1 nom1) nom1 1 1 1;
  701.     lt1=lt1 et (prog temp1);
  702.     ls1=ls1 et (prog ds1);
  703.   fin bou1;
  704.   ecc1.&bou2=(evol manu 'temps' lt1 nom1 ls1);
  705.   tab2.'TITRE'. &bou2 = nom1 ;
  706. fin bou2;
  707. tab2.1 = 'MARQ CROI ' ;
  708. tab2.2 = 'MARQ PLUS ' ;
  709. tab2.3 = 'MARQ ETOI ' ;
  710. tab2.4 = 'MARQ CARR ' ;
  711. tab2.5 = 'MARQ LOSA ' ;
  712. tab2.6 = 'MARQ ETOI NOLI' ;
  713. tab2.7 = 'MARQ CARR NOLI' ;
  714. dess ( ecc1.1 et ecc1.2 et ecc1.3 et ecc1.4 
  715.  et ecc1.5 et ecc1.6  ) lege tab2;
  716.  
  717. ************************************************************************
  718. * deformations plastiques dans les phases
  719. ************************************************************************
  720.  
  721.  
  722. n1=dime tab1.temps;
  723. ecc1=table;
  724. tab2=table;
  725. tab2.'TITRE'=table;
  726. repe bou2 6;
  727.   lt1=prog;
  728.   ls1=prog;
  729.   nom1=chain 'RT0' (&bou2);
  730.   repe bou1 n1;
  731.     vari1=tab1.variables_internes.(&bou1-1);
  732.     temp1=tab1.temps.(&bou1-1);
  733.     ds1=extr (exco vari1 nom1) nom1 1 1 1;
  734.     lt1=lt1 et (prog temp1);
  735.     ls1=ls1 et (prog ds1);
  736.   fin bou1;
  737.   ecc1.&bou2=(evol manu 'temps' lt1 nom1 ls1);
  738.   tab2.'TITRE'. &bou2 = nom1 ;
  739. fin bou2;
  740. tab2.1 = 'MARQ CROI ' ;
  741. tab2.2 = 'MARQ PLUS ' ;
  742. tab2.3 = 'MARQ ETOI ' ;
  743. tab2.4 = 'MARQ CARR ' ;
  744. tab2.5 = 'MARQ LOSA ' ;
  745. tab2.6 = 'MARQ ETOI NOLI' ;
  746. tab2.7 = 'MARQ CARR NOLI' ;
  747. dess ( ecc1.1 et ecc1.2 et ecc1.3 et ecc1.4 
  748.  et ecc1.5 et ecc1.6  ) lege tab2;
  749.  
  750. n1=dime tab1.temps;
  751. ecc1=table;
  752. tab2=table;
  753. tab2.'TITRE'=table;
  754. repe bou2 6;
  755.   lt1=prog;
  756.   ls1=prog;
  757.   nom1=chain 'RT1' (&bou2);
  758.   repe bou1 n1;
  759.     vari1=tab1.variables_internes.(&bou1-1);
  760.     temp1=tab1.temps.(&bou1-1);
  761.     ds1=extr (exco vari1 nom1) nom1 1 1 1;
  762.     lt1=lt1 et (prog temp1);
  763.     ls1=ls1 et (prog ds1);
  764.   fin bou1;
  765.   ecc1.&bou2=(evol manu 'temps' lt1 nom1 ls1);
  766.   tab2.'TITRE'. &bou2 = nom1 ;
  767. fin bou2;
  768. tab2.1 = 'MARQ CROI ' ;
  769. tab2.2 = 'MARQ PLUS ' ;
  770. tab2.3 = 'MARQ ETOI ' ;
  771. tab2.4 = 'MARQ CARR ' ;
  772. tab2.5 = 'MARQ LOSA ' ;
  773. tab2.6 = 'MARQ ETOI NOLI' ;
  774. tab2.7 = 'MARQ CARR NOLI' ;
  775. dess ( ecc1.1 et ecc1.2 et ecc1.3 et ecc1.4 
  776.  et ecc1.5 et ecc1.6  ) lege tab2;
  777.  
  778. n1=dime tab1.temps;
  779. ecc1=table;
  780. tab2=table;
  781. tab2.'TITRE'=table;
  782. repe bou2 6;
  783.   lt1=prog;
  784.   ls1=prog;
  785.   nom1=chain 'IT1' (&bou2);
  786.   repe bou1 n1;
  787.     vari1=tab1.variables_internes.(&bou1-1);
  788.     temp1=tab1.temps.(&bou1-1);
  789.     ds1=extr (exco vari1 nom1) nom1 1 1 1;
  790.     lt1=lt1 et (prog temp1);
  791.     ls1=ls1 et (prog ds1);
  792.   fin bou1;
  793.   ecc1.&bou2=(evol manu 'temps' lt1 nom1 ls1);
  794.   tab2.'TITRE'. &bou2 = nom1 ;
  795. fin bou2;
  796. tab2.1 = 'MARQ CROI ' ;
  797. tab2.2 = 'MARQ PLUS ' ;
  798. tab2.3 = 'MARQ ETOI ' ;
  799. tab2.4 = 'MARQ CARR ' ;
  800. tab2.5 = 'MARQ LOSA ' ;
  801. tab2.6 = 'MARQ ETOI NOLI' ;
  802. tab2.7 = 'MARQ CARR NOLI' ;
  803. dess ( ecc1.1 et ecc1.2 et ecc1.3 et ecc1.4 
  804.  et ecc1.5 et ecc1.6  ) lege tab2;
  805.  
  806. n1=dime tab1.temps;
  807. ecc1=table;
  808. tab2=table;
  809. tab2.'TITRE'=table;
  810. repe bou2 6;
  811.   lt1=prog;
  812.   ls1=prog;
  813.   nom1=chain 'OT1' (&bou2);
  814.   repe bou1 n1;
  815.     vari1=tab1.variables_internes.(&bou1-1);
  816.     temp1=tab1.temps.(&bou1-1);
  817.     ds1=extr (exco vari1 nom1) nom1 1 1 1;
  818.     lt1=lt1 et (prog temp1);
  819.     ls1=ls1 et (prog ds1);
  820.   fin bou1;
  821.   ecc1.&bou2=(evol manu 'temps' lt1 nom1 ls1);
  822.   tab2.'TITRE'. &bou2 = nom1 ;
  823. fin bou2;
  824. tab2.1 = 'MARQ CROI ' ;
  825. tab2.2 = 'MARQ PLUS ' ;
  826. tab2.3 = 'MARQ ETOI ' ;
  827. tab2.4 = 'MARQ CARR ' ;
  828. tab2.5 = 'MARQ LOSA ' ;
  829. tab2.6 = 'MARQ ETOI NOLI' ;
  830. tab2.7 = 'MARQ CARR NOLI' ;
  831. dess ( ecc1.1 et ecc1.2 et ecc1.3 et ecc1.4 
  832.  et ecc1.5 et ecc1.6  ) lege tab2;
  833.  
  834. ************************************************************************
  835. * chemo plasticite dans les phases
  836. ************************************************************************
  837.  
  838. n1=dime tab1.temps;
  839. ecc1=table;
  840. tab2=table;
  841. tab2.'TITRE'=table;
  842. repe bou2 6;
  843.   lt1=prog;
  844.   ls1=prog;
  845.   nom1=chain 'TS1' (&bou2);
  846.   repe bou1 n1;
  847.     vari1=tab1.variables_internes.(&bou1-1);
  848.     temp1=tab1.temps.(&bou1-1);
  849.     ds1=extr (exco vari1 nom1) nom1 1 1 1;
  850.     lt1=lt1 et (prog temp1);
  851.     ls1=ls1 et (prog ds1);
  852.   fin bou1;
  853.   ecc1.&bou2=(evol manu 'temps' lt1 nom1 ls1);
  854.   tab2.'TITRE'. &bou2 = nom1 ;
  855. fin bou2;
  856. tab2.1 = 'MARQ CROI ' ;
  857. tab2.2 = 'MARQ PLUS ' ;
  858. tab2.3 = 'MARQ ETOI ' ;
  859. tab2.4 = 'MARQ CARR ' ;
  860. tab2.5 = 'MARQ LOSA ' ;
  861. tab2.6 = 'MARQ ETOI NOLI' ;
  862. tab2.7 = 'MARQ CARR NOLI' ;
  863. dess ( ecc1.1 et ecc1.2 et ecc1.3 et ecc1.4 
  864.  et ecc1.5 et ecc1.6  ) lege tab2; 
  865.  
  866. n1=dime tab1.temps;
  867. ecc1=table;
  868. tab2=table;
  869. tab2.'TITRE'=table;
  870. repe bou2 6;
  871.   lt1=prog;
  872.   ls1=prog;
  873.   nom1=chain 'TS0' (&bou2);
  874.   repe bou1 n1;
  875.     vari1=tab1.variables_internes.(&bou1-1);
  876.     temp1=tab1.temps.(&bou1-1);
  877.     ds1=extr (exco vari1 nom1) nom1 1 1 1;
  878.     lt1=lt1 et (prog temp1);
  879.     ls1=ls1 et (prog ds1);
  880.   fin bou1;
  881.   ecc1.&bou2=(evol manu 'temps' lt1 nom1 ls1);
  882.   tab2.'TITRE'. &bou2 = nom1 ;
  883. fin bou2;
  884. tab2.1 = 'MARQ CROI ' ;
  885. tab2.2 = 'MARQ PLUS ' ;
  886. tab2.3 = 'MARQ ETOI ' ;
  887. tab2.4 = 'MARQ CARR ' ;
  888. tab2.5 = 'MARQ LOSA ' ;
  889. tab2.6 = 'MARQ ETOI NOLI' ;
  890. tab2.7 = 'MARQ CARR NOLI' ;
  891. dess ( ecc1.1 et ecc1.2 et ecc1.3 et ecc1.4 
  892.  et ecc1.5 et ecc1.6  ) lege tab2;
  893.  
  894. *-----------------------------------------------------------------------
  895. *evol de quelques variables internes
  896. *-----------------------------------------------------------------------
  897. lt0=prog;
  898. tabn2=table;
  899.  
  900. tabn2.1='RL11';
  901. tabn2.2='RL12';
  902. tabn2.3='RL13';
  903. tabn2.4='RL14';
  904. tabn2.5='RL15';
  905. tabn2.6='RL16';
  906.  
  907. tabl2=table;
  908. tabl2.1=prog;
  909. tabl2.2=prog;
  910. tabl2.3=prog;
  911. tabl2.4=prog;
  912. tabl2.5=prog;
  913. tabl2.6=prog;
  914.  
  915. tabe5=table;
  916. tabe5.1=bleu;
  917. tabe5.2=roug;
  918. tabe5.3=rose;
  919. tabe5.4=ROUG;
  920. tabe5.5=CYAN;
  921. tabe5.6=VERT;
  922.  
  923. tabe6=table;
  924. tabe6.TITRE=table;
  925. tabe6.TITRE.1='RL11';
  926. tabe6.TITRE.2='RL12';
  927. tabe6.TITRE.3='RL13';
  928. tabe6.TITRE.4='RL14';
  929. tabe6.TITRE.5='RL15';
  930. tabe6.TITRE.6='RL16';
  931. tabe6.1= 'MARQ CROI';
  932. tabe6.2= 'MARQ PLUS';
  933. tabe6.3= 'MARQ CARR';
  934. tabe6.4= 'MARQ CROI';
  935. tabe6.5= 'MARQ PLUS';
  936. tabe6.6= 'MARQ CARR';
  937.  
  938. dim0=dime tab1.temps;
  939.  
  940. repe bouc1 (dim0-1);
  941.  vari1=tab1.variables_internes.(&bouc1);
  942.  repe bouc2 (dime(tabn2));
  943.   t3=exco (tabn2.(&bouc2)) vari1;
  944.   t2=extr t3 (tabn2.(&bouc2))  1 1 1;
  945.   tabl2.(&bouc2)=tabl2.(&bouc2) et (prog t2);
  946.  fin bouc2;
  947.  t0=tab1.temps.(&bouc1);
  948.  lt0=lt0 et (prog t0);
  949. fin bouc1;
  950.  
  951. tabe3=table;
  952. tabe5.'TITRE'=table;
  953. x1=(log(10))**-1;
  954. lnt0=x1*(log(lt0));
  955.  
  956. repe bouc3 (dime tabn2);
  957.  tabe5.'TITRE'.(&bouc3)=tabn2.(&bouc3); 
  958.  si((&bouc3) ega 1);
  959.   evol1=evol manu 'temps' lt0 
  960.   tabn2.(&bouc3) tabl2.(&bouc3);
  961.   evol2=evol manu 'temps' lnt0 
  962.   tabn2.(&bouc3) tabl2.(&bouc3);
  963.  sinon;
  964.   evol1=evol1 et(evol manu 'temps' lt0 
  965.   tabn2.(&bouc3) tabl2.(&bouc3) coul tabe5.&bouc3);
  966.   evol2=evol2 et (evol manu 'temps' lnt0 
  967.   tabn2.(&bouc3) tabl2.(&bouc3) coul tabe5.&bouc3);
  968.  finsi;
  969. fin bouc3;
  970.  
  971. dess evol1 'LEGE' tabe6;
  972. finsi;
  973. fin;
  974.  
  975.  
  976.  
  977.  
  978.  
  979.  
  980.  

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