Télécharger la plaquette Cast3M


fr en

Télécharger fluendo3d_def_rag_thcm.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : fluendo3d_rag_def_thcm.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. * test de la formulation RAG(AAR) et RSI(DEF) du modele fluendo3d
  6. * ----------------------------------------------------------------------
  7. *
  8. * Alain Sellier, Stephane Multon, Pierre Morenon
  9. * mercredi 4 janvier 2023
  10. *
  11. *-----------------------------------------------------------------------
  12. * Exemple de calcul d un element en beton soumis a la RAG et a la RSI
  13. * en condition thermo-hydro-chimiques variables et gonflement libre en
  14. * un premier temps, et soumis a une contrainte axiale de 10 MPa par la
  15. * suite
  16. *-----------------------------------------------------------------------
  17. graph1=FAUX;
  18. opti dime 2 elem qua4 mode axis;
  19. *geometrie
  20. p1=0. 0.;
  21. p2=0.065 0.;
  22. p3=0.065 0.24;
  23. p4=0. 0.24;
  24. n1=1;
  25. d12=d n1 p1 p2;
  26. d23=d n1 p2 p3;
  27. d34=d n1 p3 p4;
  28. d41=d n1 p4 p1;
  29. surf1=daller d12  d23  d34  d41;
  30. *trac surf1;
  31.  
  32. ****************************modele et materiau *************************
  33.  
  34. mod1=mode surf1 mecanique elastique 
  35. viscoplastique FLUENDO3D;
  36.  
  37. *****beton *************************************************************
  38. *rc et E du beton
  39. rc180=72.;
  40. ec180=44000.;
  41.  
  42. *caracteristiques  beton deduite avec regles EC2
  43. rt180=4.4;
  44. moyrt00=4.4;
  45. sigrt00=0.10*moyrt00;
  46. rt0 = moyrt00;
  47. E0=ec180;
  48. ep00=moyrt00/E0;
  49. ep0=(E0**-1)*rt0;
  50.  
  51. *caracteristiques des aciers
  52. Eacier=1.82e5;
  53. fyacier=430.; 
  54. *contrainte equivalente de cisaillement pour le calcul de ls (MC2010)
  55. taueq0=1.8*moyrt00;
  56.  
  57. *evolution CNAD-NA
  58. evcna1=evol manu 'NA' (prog 0. 1.) 'NALD' (prog 0. 1.);
  59.  
  60. *evolution VW-W
  61. evcw1=evol manu 'W' (prog 0. 1.) 'VW' (prog 0. 1.);
  62.  
  63. *finite element size for Hillerborgh method
  64. tail1=tail direction mod1; 
  65. ray1=coor 1 surf1;
  66. ray2=0.25*(chan 'CHAM' ray1 mod1 'RIGIDITE'  'CARACTERISTIQUES');
  67.  
  68.  
  69. *elasticity, plasticity and damage
  70. matd1=MATE mod1 YOUN e0 NU 0.2 RHO 2400. ALPH 1.0e-5 NREN 0 DIM3 0.1
  71.  YORF e0 NURF 0.2 HYDR 1. HREF 1. HYDS 0.2 RT rt0 EPT 1.1e-4 GFT 1.e-4
  72.  RC rc180 EPC 2.e-3 DELT 1.  BETA 0.15 EKDC 2.0e-3 REF 5. GFR 1.0e-4 
  73.  TSTH 45. DT80 0.15 TREF 20. TALP 20.;
  74.  
  75. *water
  76. matgw1=MATE mod1 TTKW 40. VW evcw1 PORO 0.12 BSHR 0.25 MSHR 44. 
  77.  MVGN 0.49  SKDW 17. DCDW 1. CSHR 0.50 HSHR (3.0*e0) TTRW 20. 
  78.  KWRT 0.25 KWRC 0.1 ; 
  79.  
  80. *visco elasticity and visco plasticity (creep)
  81. matf1=MATE mod1 TTRF 20. TAUK 2. YKSY 3.5 TAUM 2. EKFL 3.8e-4
  82. XFLU 2. NRJM 45000. DFMX 0. MDTT 15000. TDTT 1.0 WDTT 0.02 PDTT 1.0 ;
  83.  
  84. *rag
  85. matgr1=MATE mod1  VRAG 1.0e-2  TRAG 120. NRJR 40000. TTRG 38. DCDG 0.15
  86. SRSR 0.5 VVRG 1.0e-3 KRAG 2000. CRAG 1. HRAG 1000. EKDG 3.0e-3;
  87.  
  88. *def
  89. matgd1=MATE mod1  VDEF 6.0e-2  SSAD 0.6 NALD evcna1 TPRD 30. TTKD 80.
  90. NAKD 0.28 NRJP 40000. TTRP 20. SRSD 0.95  NABD 0.92   TDID (65./24.) 
  91. NRJD 70000. TTKF 70. TFID (30./24.) NRJF 180000.  EXND 0.18 EXMD 3. 
  92. VVDF 6.0e-3 KDEF 2000. CDEF 1. HDEF 1000. EKDS 3.0e-4 DCDS 0.015;
  93.  
  94. *Weibull scale effect : note cf below to prepar Vmax field
  95. matw1=MATE mod1 VREF 2.1e-4 VMAX 2.2e-3 CVRT 0.08;
  96.  
  97. *material chracteristics assemblage
  98. mat1=matd1 ET matf1  ET  matgw1 et matgr1 ET matgd1 ET matw1
  99.  ET tail1;
  100.  
  101. *******************conditions aux limites en pression ******************
  102. cl1=bloq p1 UR UZ;
  103. cl2=bloq d12 UZ;
  104. cl3=bloq d41 UR;
  105. sig2=press mass mod1 d34 1.;
  106. *chargement vertical
  107. evsig2=evol manu 'temps'
  108.             (prog 0.  350. 351. 800.)
  109.             'coeff'
  110.             (prog 1.0e-2 1.0e-2  10. 10.) coul bleu;
  111. si (graph1) ;
  112.     dess evsig2 titr 'contrainte verticale';
  113. finsi ; 
  114. charsig2=char 'MECA' sig2 evsig2; 
  115. sig3=press mass mod1 d23 1.;
  116. *chargement radial
  117. evsig3=evol manu 'temps'
  118.             (prog 0.  350. 351. 800.)
  119.             'coeff'
  120.             (prog 1.0e-2 1.0e-2 1.0e-2 1.0e-2) coul bleu;
  121. si (graph1) ;
  122.     dess evsig3 titr 'contrainte radiale'; 
  123. finsi ;
  124.  
  125. charsig3=char 'MECA' sig3 evsig3; 
  126.  
  127. *******************histoire de temperature *****************************
  128. chth1=manu chpo  surf1 'T' 1.;
  129. evth1=evol manu 'temps' (prog 0. 1. 20. 21. 48. 49. 800.) 'coef' 
  130.                         (prog 20. 65. 65. 20. 20. 38. 38.);
  131. si(graph1) ;
  132.      dess evth1 titr 'temperature(C)-temps(j)';
  133. finsi ;
  134. chath1=char 'T' chth1 evth1;
  135.  
  136. ******************histoire de alcalins *********************************
  137. chna1=manu chpo  surf1 'NA' 1.;
  138. evna1=evol manu 'temps' (prog 0.  800.) 'coef' 
  139.                         (prog 0.92 0.0);
  140. si(graph1) ;
  141.   dess evna1 titr 'Na(mol/l)-temps(j)';
  142. finsi ;
  143. chana1=char 'NA' chna1 evna1;
  144.  
  145. *****************histoire hydrique *************************************
  146. chw1=manu chpo  surf1 'W' 0.12;
  147. ltw1=prog 0. pas 10. 800.;
  148. sr0=0.85;
  149. sr1=1.0;
  150. tsech1=50.;
  151. coeft1=-1.*(tsech1**-1);
  152. lcw1=((sr0-sr1)*(exp((coeft1)*ltw1)))+(sr1*(exp(0.*ltw1)));
  153. evw1=evol manu 'temps' ltw1 'coef'  lcw1;
  154. si (graph1);
  155.   dess evw1 titr 'Evolution Sr(t)';
  156. finsi;
  157. chaw1=char 'W' chw1 evw1;
  158.  
  159. ************************************************************************
  160. tab1=table;
  161. ltc1=prog 0. pas 0.5 1. pas 1. 30. pas 2. 60. pas 5. 350. pas 0.1 352. 
  162.  pas 2. 450. pas 5.  500. ;  
  163. tab1.temps_calcules=ltc1;
  164. tab1.caracteristiques=mat1;
  165. tab1.modele=mod1;
  166. tab1.blocages_mecaniques=cl1 et cl2 et cl3;
  167. tab1.mova='VAFT';
  168. tab1.chargement=charsig2  et chana1 et chath1 et chaw1 et charsig3;
  169. tab1.precision=1.0e-5;
  170. tab1.processeurs=mot 'COMPORTEMENT';
  171. pasapas tab1;
  172.  
  173.  
  174. ************************************************************************
  175. *exploitation des resultats
  176. ************************************************************************
  177.  
  178. *-----------------------------------------------------------------------
  179.  
  180. *courbe de deplacement au centre de l eprouvette
  181. *et variables internes
  182. n1=dime (tab1.temps);
  183. i1=0;
  184. dp3=prog;
  185. lt1=prog;
  186. epr1=prog;
  187. epz1=prog;
  188. repeter bloc0 n1;
  189.  z1=extr tab1.deplacements.i1 UZ p1;
  190.  depi1=tab1.deplacements.i1;
  191.  epsi1=epsi depi1 mod1;
  192.  sig1=tab1.contraintes.i1;
  193.  epzz1=extr ( exco epsi1 'EPZZ') 'EPZZ' 1 1 1;
  194.  szz1=extr ( exco sig1 'SMZZ') 'SMZZ' 1 1 1;
  195.  eprr1=extr(exco epsi1 'EPRR') 'EPRR' 1 1 1;
  196.  srr1=extr ( exco sig1 'SMRR') 'SMRR' 1 1 1;
  197.  epr1=epr1 et (prog (eprr1));
  198.  epz1=epz1 et (prog (epzz1));
  199.  lt1=lt1 et (prog (tab1.temps.i1));
  200.  i1=i1+1;
  201. fin bloc0;
  202.  
  203. *courbes deformations temps theorique
  204. evezth1=(evol manu 'TEMPS' lt1 'EPZZ' (1.e6*epz1)) coul roug;
  205. everth1=(evol manu 'TEMPS' lt1 'EPRR' (1.e6*epr1)) coul vert;
  206. si (graph1);
  207.     dess (evezth1 et everth1) titr 'deplacements axial et radial';
  208. finsi ;
  209.  
  210.  
  211. *-----------------------------------------------------------------------
  212. *test de conformite
  213. test1=VRAI;
  214. si (TEST1) ;
  215.     n1=184; fref1 = 2.8984E-03;
  216.     n2=184; fref2 = 6.6484E-03;
  217.     fcal1='EXTR' epz1 (n1);
  218.     fcal2='EXTR' epr1 (n2);
  219.     err1 = abs ((fref1 - fcal1)/fref1) ; 
  220.     err2 = abs ((fref2 - fcal2)/fref2) ;
  221.     si ((err1 < 1.0e-3) et (err2 < 1.0e-3)) ;
  222.         erre 0 ;
  223.     sinon ;
  224.         erre 5 ;
  225.     finsi ;
  226. finsi;
  227.  
  228.  
  229. si (graph1) ;
  230.  
  231. *-----------------------------------------------------------------------
  232. * evolution de la chimie
  233.  
  234. lt0=prog;
  235.  
  236. tabn2=table;
  237. tabn2.1='AFT1';
  238. tabn2.2='AFM1';
  239. tabn2.3='AFT2';
  240. tabn2.4='VAFT';
  241. tabn2.5='ATIL';
  242. tabn2.6='STIL';
  243.  
  244. tabl2=table;
  245. tabl2.1=prog;
  246. tabl2.2=prog;
  247. tabl2.3=prog;
  248. tabl2.4=prog;
  249. tabl2.5=prog;
  250. tabl2.6=prog;
  251.  
  252. tabe5=table;
  253. tabe5.1=bleu;
  254. tabe5.2=roug;
  255. tabe5.3=rose;
  256. tabe5.4=ROUG;
  257. tabe5.5=CYAN;
  258. tabe5.6=VERT;
  259. tabe6=table;
  260. tabe6.TITRE=table;
  261. tabe6.TITRE.1='AFT1';
  262. tabe6.TITRE.2='AFM1';
  263. tabe6.TITRE.3='AFT2';
  264. tabe6.TITRE.4='VAFT';
  265. tabe6.TITRE.5='ATIL';
  266. tabe6.TITRE.6='STIL';
  267. tabe6.1= 'MARQ CROI';
  268. tabe6.2= 'MARQ ROND';
  269. tabe6.3= 'MARQ CARR';
  270. tabe6.4= 'MARQ CROI';
  271. tabe6.5= 'MARQ PLUS';
  272. tabe6.6= 'MARQ CARR';
  273.  
  274. dim0=dime tab1.temps;
  275.  
  276. repe bouc1 (dim0-1);
  277.  vari1=tab1.variables_internes.(&bouc1);
  278.  repe bouc2 (dime(tabn2));
  279.   t3=exco (tabn2.(&bouc2)) vari1;
  280.   t2=extr t3 (tabn2.(&bouc2))  1 1 1;
  281.   tabl2.(&bouc2)=tabl2.(&bouc2) et (prog t2);
  282.  fin bouc2;
  283.  t0=tab1.temps.(&bouc1);
  284.  lt0=lt0 et (prog t0);
  285. fin bouc1;
  286.  
  287. tabe3=table;
  288. tabe5.'TITRE'=table;
  289. x1=(log(10))**-1;
  290. lnt0=x1*(log(lt0));
  291.  
  292. repe bouc3 (dime(tabn2));
  293.  tabe5.'TITRE'.(&bouc3)=tabn2.(&bouc3); 
  294.  si((&bouc3) ega 1);
  295.   evol1=evol manu 'temps' lt0 
  296.   tabn2.(&bouc3) tabl2.(&bouc3);
  297.   evol2=evol manu 'temps' lnt0 
  298.   tabn2.(&bouc3) tabl2.(&bouc3);
  299.  sinon;
  300.   evol1=evol1 et(evol manu 'temps' lt0 
  301.   tabn2.(&bouc3) tabl2.(&bouc3) coul tabe5.&bouc3);
  302.   evol2=evol2 et (evol manu 'temps' lnt0 
  303.   tabn2.(&bouc3) tabl2.(&bouc3) coul tabe5.&bouc3);
  304.  finsi;
  305. fin bouc3;
  306. si (graph1);
  307. dess evol1 'LEGE' tabe6;
  308. finsi; 
  309.  
  310. *-----------------------------------------------------------------------
  311.  
  312. lrag1=prog;
  313. ldef1=prog;
  314. ltemp1=prog;
  315. repe bouc2 (dim0-1);
  316.  t3=tab1.variables_internes.(&bouc2);
  317.  aar1=extr t3 'AGEL' 1 1 1 ;
  318.  lrag1=lrag1 et (prog aar1);
  319.  ade1=extr t3 'ADEF' 1 1 1 ;
  320.  ldef1=ldef1 et (prog ade1);
  321.  ltemp1=ltemp1 et (prog (tab1.temps.(&bouc2)));
  322. fin bouc2;
  323. evar1=(evol manu 'temps' ltemp1 'aar' lrag1) coul bleu;
  324. evde1=(evol manu 'temps' ltemp1 'def' ldef1) coul roug;
  325. tabe6=table;
  326. tabe6.TITRE=table;
  327. tabe6.TITRE.1='AAR';
  328. tabe6.TITRE.2='DEF';
  329. dess (evar1 et evde1) 'LEGE' tabe6;
  330.  
  331. *-----------------------------------------------------------------------
  332. lrag1=prog;
  333. ldef1=prog;
  334. lwat1=prog;
  335. ltemp1=prog;
  336. repe bouc2 (dim0-1);
  337.  t3=tab1.variables_internes.(&bouc2);
  338.  aar1=extr t3 'PGEL' 1 1 1 ;
  339.  lrag1=lrag1 et (prog aar1);
  340.  ade1=extr t3 'PAFT' 1 1 1 ;
  341.  ldef1=ldef1 et (prog ade1);
  342.  awa1=extr t3 'PWAT' 1 1 1 ;
  343.  lwat1=lwat1 et (prog awa1);
  344.  ltemp1=ltemp1 et (prog (tab1.temps.(&bouc2)));
  345. fin bouc2;
  346. evar1=(evol manu 'temps' ltemp1 'PGEL' lrag1) coul rose;
  347. evde1=(evol manu 'temps' ltemp1 'PAFT' ldef1) coul roug;
  348. evwe1=(evol manu 'temps' ltemp1 'PWAT' lwat1) coul cyan;
  349. tabe6=table;
  350. tabe6.TITRE=table;
  351. tabe6.TITRE.1='PGEL';
  352. tabe6.TITRE.2='PAFT';
  353. tabe6.TITRE.3='PWAT';
  354. dess (evar1 et evde1 et evwe1) 'LEGE' tabe6;
  355.  
  356.  
  357. *points experimentaux RAG libre
  358. ltpsex1=prog 7. 28.   56.  91. 182.  273. 365.;
  359. lepzex1=1.e4*(prog 0. 0.015 0.05 0.1 0.175 0.21 0.22);
  360. leprex1=1.e4*(prog 0. 0.05  0.1  0.175 0.24 0.25 0.27);
  361.  
  362. *correction pt exp partie elastique
  363. defelz1=(0./37200.)*1.0e6;
  364. defelr1=(-0.2*defelz1);
  365. nls1=dime lepzex1;
  366. lsepex1=prog nls1 * defelz1 ;
  367. lseper1=prog nls1 * defelr1 ;
  368. lepzex1=lepzex1+lsepex1;
  369. leprex1=leprex1+lseper1; 
  370. evezex1=evol manu 'temps' ltpsex1 'epzz exp' lepzex1; 
  371. everex1=evol manu 'temps' ltpsex1 'eprr exp' leprex1;
  372.  
  373. tab3=table;
  374. tab3.1 = 'MARQ CROI NOLI' ;
  375. tab3.3 = 'TIRR';
  376. tab3.4 = 'MARQ  LOSA NOLI';
  377. tab3.'TITRE' = 'TABLE' ;
  378. tab3.'TITRE'. 1 = MOT 'epzz exp' ;
  379. tab3.'TITRE'. 2 = MOT 'epzz th' ;
  380. tab3.'TITRE'. 3 = MOT 'eprr th' ;
  381. tab3.'TITRE'. 4 = MOT 'eprr exp';
  382. dess (evezex1 et evezth1 et everth1 et everex1) lege tab3; 
  383.  
  384. OPTI 'SORT' 'Def_et_rag_axial_00_szz_ezz.xls';
  385. SORT 'EXCE' evezth1 evezex1 'SEPA' 'PVIR';
  386. OPTI 'SORT' 'Def_et_rag_axial_00_srr_err.xls';
  387. SORT 'EXCE' everth1 everex1 'SEPA' 'PVIR';
  388.  
  389. *déformations sur l'éprouvette
  390.  
  391. def0=defo tab1.deplacements.(n1-1) surf1 0.;
  392. def1=defo tab1.deplacements.(n1-1) surf1  roug;
  393. *epsi1=epsi def1 mod1;
  394. trac (def1 et def0);
  395. *------------------------------------------------------------------------
  396. *evol de quelques variables internes
  397. lt0=prog;
  398.  
  399. tabn2=table;
  400. tabn2.1='EVT1';
  401. tabn2.2='EVT2';
  402. tabn2.3='EVT3';
  403. tabn2.4='EVT4';
  404. tabn2.5='EVT5';
  405. tabn2.6='EVT6';
  406.  
  407. tabl2=table;
  408. tabl2.1=prog;
  409. tabl2.2=prog;
  410. tabl2.3=prog;
  411. tabl2.4=prog;
  412. tabl2.5=prog;
  413. tabl2.6=prog;
  414.  
  415. tabe5=table;
  416. tabe5.1=bleu;
  417. tabe5.2=roug;
  418. tabe5.3=rose;
  419. tabe5.4=ROUG;
  420. tabe5.5=CYAN;
  421. tabe5.6=VERT;
  422.  
  423. tabe6=table;
  424. tabe6.TITRE=table;
  425. tabe6.TITRE.1='EVT1';
  426. tabe6.TITRE.2='EVT2';
  427. tabe6.TITRE.3='EVT3';
  428. tabe6.TITRE.4='EVT4';
  429. tabe6.TITRE.5='EVT5';
  430. tabe6.TITRE.6='EVT6';
  431. tabe6.1= 'MARQ CROI';
  432. tabe6.2= 'MARQ ROND';
  433. tabe6.3= 'MARQ CARR';
  434. tabe6.4= 'MARQ CROI';
  435. tabe6.5= 'MARQ PLUS';
  436. tabe6.6= 'MARQ CARR';
  437.  
  438. dim0=dime tab1.temps;
  439.  
  440. repe bouc1 (dim0-1);
  441.  vari1=tab1.variables_internes.(&bouc1);
  442.  repe bouc2 (dime(tabn2));
  443.   t3=exco (tabn2.(&bouc2)) vari1;
  444.   t2=extr t3 (tabn2.(&bouc2))  1 1 1;
  445.   tabl2.(&bouc2)=tabl2.(&bouc2) et (prog t2);
  446.  fin bouc2;
  447.  t0=tab1.temps.(&bouc1);
  448.  lt0=lt0 et (prog t0);
  449. fin bouc1;
  450.  
  451. tabe3=table;
  452. tabe5.'TITRE'=table;
  453. x1=(log(10))**-1;
  454. lnt0=x1*(log(lt0));
  455.  
  456. repe bouc3 (dime(tabn2));
  457.  tabe5.'TITRE'.(&bouc3)=tabn2.(&bouc3); 
  458.  si((&bouc3) ega 1);
  459.   evol1=evol manu 'temps' lt0 
  460.   tabn2.(&bouc3) tabl2.(&bouc3);
  461.   evol2=evol manu 'temps' lnt0 
  462.   tabn2.(&bouc3) tabl2.(&bouc3);
  463.  sinon;
  464.   evol1=evol1 et(evol manu 'temps' lt0 
  465.   tabn2.(&bouc3) tabl2.(&bouc3) coul tabe5.&bouc3);
  466.   evol2=evol2 et (evol manu 'temps' lnt0 
  467.   tabn2.(&bouc3) tabl2.(&bouc3) coul tabe5.&bouc3);
  468.  finsi;
  469. fin bouc3;
  470.  
  471. dess evol1 'LEGE' tabe6;
  472.  
  473. *-----------------------------------------------------------------------
  474.  
  475. *------------------------------------------------------------------------
  476. *evol de quelques variables internes
  477. lt0=prog;
  478.  
  479. tabn2=table;
  480. tabn2.1='EPG1';
  481. tabn2.2='EPG2';
  482. tabn2.3='EPG3';
  483. tabn2.4='EPG4';
  484. tabn2.5='EPG5';
  485. tabn2.6='EPG6';
  486.  
  487. tabl2=table;
  488. tabl2.1=prog;
  489. tabl2.2=prog;
  490. tabl2.3=prog;
  491. tabl2.4=prog;
  492. tabl2.5=prog;
  493. tabl2.6=prog;
  494.  
  495. tabe5=table;
  496. tabe5.1=bleu;
  497. tabe5.2=roug;
  498. tabe5.3=rose;
  499. tabe5.4=ROUG;
  500. tabe5.5=CYAN;
  501. tabe5.6=VERT;
  502.  
  503. tabe6=table;
  504. tabe6.TITRE=table;
  505. tabe6.TITRE.1='EPG1';
  506. tabe6.TITRE.2='EPG2';
  507. tabe6.TITRE.3='EPG3';
  508. tabe6.TITRE.4='EPG4';
  509. tabe6.TITRE.5='EPG5';
  510. tabe6.TITRE.6='EPG6';
  511. tabe6.1= 'MARQ CROI';
  512. tabe6.2= 'MARQ ROND';
  513. tabe6.3= 'MARQ CARR';
  514. tabe6.4= 'MARQ CROI';
  515. tabe6.5= 'MARQ PLUS';
  516. tabe6.6= 'MARQ CARR';
  517.  
  518. dim0=dime tab1.temps;
  519.  
  520. repe bouc1 (dim0-1);
  521.  vari1=tab1.variables_internes.(&bouc1);
  522.  repe bouc2 (dime(tabn2));
  523.   t3=exco (tabn2.(&bouc2)) vari1;
  524.   t2=extr t3 (tabn2.(&bouc2))  1 1 1;
  525.   tabl2.(&bouc2)=tabl2.(&bouc2) et (prog t2);
  526.  fin bouc2;
  527.  t0=tab1.temps.(&bouc1);
  528.  lt0=lt0 et (prog t0);
  529. fin bouc1;
  530.  
  531. tabe3=table;
  532. tabe5.'TITRE'=table;
  533. x1=(log(10))**-1;
  534. lnt0=x1*(log(lt0));
  535.  
  536. repe bouc3 (dime(tabn2));
  537.  tabe5.'TITRE'.(&bouc3)=tabn2.(&bouc3); 
  538.  si((&bouc3) ega 1);
  539.   evol1=evol manu 'temps' lt0 
  540.   tabn2.(&bouc3) tabl2.(&bouc3);
  541.   evol2=evol manu 'temps' lnt0 
  542.   tabn2.(&bouc3) tabl2.(&bouc3);
  543.  sinon;
  544.   evol1=evol1 et(evol manu 'temps' lt0 
  545.   tabn2.(&bouc3) tabl2.(&bouc3) coul tabe5.&bouc3);
  546.   evol2=evol2 et (evol manu 'temps' lnt0 
  547.   tabn2.(&bouc3) tabl2.(&bouc3) coul tabe5.&bouc3);
  548.  finsi;
  549. fin bouc3;
  550.  
  551. dess evol1 'LEGE' tabe6;
  552.  
  553. *-----------------------------------------------------------------------
  554. *-----------------------------------------------------------------------
  555. lt0=prog;
  556.  
  557. tabn2=table;
  558. tabn2.1='EPS1';
  559. tabn2.2='EPS2';
  560. tabn2.3='EPS3';
  561. tabn2.4='EPS4';
  562. tabn2.5='EPS5';
  563. tabn2.6='EPS6';
  564.  
  565. tabl2=table;
  566. tabl2.1=prog;
  567. tabl2.2=prog;
  568. tabl2.3=prog;
  569. tabl2.4=prog;
  570. tabl2.5=prog;
  571. tabl2.6=prog;
  572.  
  573. tabe5=table;
  574. tabe5.1=bleu;
  575. tabe5.2=roug;
  576. tabe5.3=rose;
  577. tabe5.4=ROUG;
  578. tabe5.5=CYAN;
  579. tabe5.6=VERT;
  580. tabe6=table;
  581. tabe6.TITRE=table;
  582. tabe6.TITRE.1='EPS1';
  583. tabe6.TITRE.2='EPS2';
  584. tabe6.TITRE.3='EPS3';
  585. tabe6.TITRE.4='EPS4';
  586. tabe6.TITRE.5='EPS5';
  587. tabe6.TITRE.6='EPS6';
  588. tabe6.1= 'MARQ CROI';
  589. tabe6.2= 'MARQ ROND';
  590. tabe6.3= 'MARQ CARR';
  591. tabe6.4= 'MARQ CROI';
  592. tabe6.5= 'MARQ PLUS';
  593. tabe6.6= 'MARQ CARR';
  594.  
  595. dim0=dime tab1.temps;
  596.  
  597. repe bouc1 (dim0-1);
  598.  vari1=tab1.variables_internes.(&bouc1);
  599.  repe bouc2 (dime(tabn2));
  600.   t3=exco (tabn2.(&bouc2)) vari1;
  601.   t2=extr t3 (tabn2.(&bouc2))  1 1 1;
  602.   tabl2.(&bouc2)=tabl2.(&bouc2) et (prog t2);
  603.  fin bouc2;
  604.  t0=tab1.temps.(&bouc1);
  605.  lt0=lt0 et (prog t0);
  606. fin bouc1;
  607.  
  608. tabe3=table;
  609. tabe5.'TITRE'=table;
  610. x1=(log(10))**-1;
  611. lnt0=x1*(log(lt0));
  612.  
  613. repe bouc3 (dime(tabn2));
  614.  tabe5.'TITRE'.(&bouc3)=tabn2.(&bouc3); 
  615.  si((&bouc3) ega 1);
  616.   evol1=evol manu 'temps' lt0 
  617.   tabn2.(&bouc3) tabl2.(&bouc3);
  618.   evol2=evol manu 'temps' lnt0 
  619.   tabn2.(&bouc3) tabl2.(&bouc3);
  620.  sinon;
  621.   evol1=evol1 et(evol manu 'temps' lt0 
  622.   tabn2.(&bouc3) tabl2.(&bouc3) coul tabe5.&bouc3);
  623.   evol2=evol2 et (evol manu 'temps' lnt0 
  624.   tabn2.(&bouc3) tabl2.(&bouc3) coul tabe5.&bouc3);
  625.  finsi;
  626. fin bouc3;
  627.  
  628. dess (evol1 ) 'LEGE' tabe6;
  629.  
  630. *-----------------------------------------------------------------------
  631. lt0=prog;
  632.  
  633. tabn2=table;
  634. tabn2.1='EPT1';
  635. tabn2.2='EPT2';
  636. tabn2.3='EPT3';
  637. tabn2.4='EPT4';
  638. tabn2.5='EPT5';
  639. tabn2.6='EPT6';
  640.  
  641. tabl2=table;
  642. tabl2.1=prog;
  643. tabl2.2=prog;
  644. tabl2.3=prog;
  645. tabl2.4=prog;
  646. tabl2.5=prog;
  647. tabl2.6=prog;
  648.  
  649. tabe5=table;
  650. tabe5.1=bleu;
  651. tabe5.2=roug;
  652. tabe5.3=rose;
  653. tabe5.4=ROUG;
  654. tabe5.5=CYAN;
  655. tabe5.6=VERT;
  656. tabe6=table;
  657. tabe6.TITRE=table;
  658. tabe6.TITRE.1='EPT1';
  659. tabe6.TITRE.2='EPT2';
  660. tabe6.TITRE.3='EPT3';
  661. tabe6.TITRE.4='EPT4';
  662. tabe6.TITRE.5='EPT5';
  663. tabe6.TITRE.6='EPT6';
  664. tabe6.1= 'MARQ CROI';
  665. tabe6.2= 'MARQ ROND';
  666. tabe6.3= 'MARQ CARR';
  667. tabe6.4= 'MARQ CROI';
  668. tabe6.5= 'MARQ PLUS';
  669. tabe6.6= 'MARQ CARR';
  670.  
  671. dim0=dime tab1.temps;
  672.  
  673. repe bouc1 (dim0-1);
  674.  vari1=tab1.variables_internes.(&bouc1);
  675.  repe bouc2 (dime(tabn2));
  676.   t3=exco (tabn2.(&bouc2)) vari1;
  677.   t2=extr t3 (tabn2.(&bouc2))  1 1 1;
  678.   tabl2.(&bouc2)=tabl2.(&bouc2) et (prog t2);
  679.  fin bouc2;
  680.  t0=tab1.temps.(&bouc1);
  681.  lt0=lt0 et (prog t0);
  682. fin bouc1;
  683.  
  684. tabe3=table;
  685. tabe5.'TITRE'=table;
  686. x1=(log(10))**-1;
  687. lnt0=x1*(log(lt0));
  688.  
  689. repe bouc3 (dime(tabn2));
  690.  tabe5.'TITRE'.(&bouc3)=tabn2.(&bouc3); 
  691.  si((&bouc3) ega 1);
  692.   evol1=evol manu 'temps' lt0 
  693.   tabn2.(&bouc3) tabl2.(&bouc3);
  694.   evol2=evol manu 'temps' lnt0 
  695.   tabn2.(&bouc3) tabl2.(&bouc3);
  696.  sinon;
  697.   evol1=evol1 et(evol manu 'temps' lt0 
  698.   tabn2.(&bouc3) tabl2.(&bouc3) coul tabe5.&bouc3);
  699.   evol2=evol2 et (evol manu 'temps' lnt0 
  700.   tabn2.(&bouc3) tabl2.(&bouc3) coul tabe5.&bouc3);
  701.  finsi;
  702. fin bouc3;
  703.  
  704. dess evol1 'LEGE' tabe6;
  705. *-----------------------------------------------------------------------
  706. *-----------------------------------------------------------------------
  707. lt0=prog;
  708.  
  709. tabn2=table;
  710. tabn2.1='EPC1';
  711. tabn2.2='EPC2';
  712. tabn2.3='EPC3';
  713. tabn2.4='EPC4';
  714. tabn2.5='EPC5';
  715. tabn2.6='EPC6';
  716.  
  717. tabl2=table;
  718. tabl2.1=prog;
  719. tabl2.2=prog;
  720. tabl2.3=prog;
  721. tabl2.4=prog;
  722. tabl2.5=prog;
  723. tabl2.6=prog;
  724.  
  725. tabe5=table;
  726. tabe5.1=bleu;
  727. tabe5.2=roug;
  728. tabe5.3=rose;
  729. tabe5.4=ROUG;
  730. tabe5.5=CYAN;
  731. tabe5.6=VERT;
  732. tabe6=table;
  733. tabe6.TITRE=table;
  734. tabe6.TITRE.1='EPC1';
  735. tabe6.TITRE.2='EPC2';
  736. tabe6.TITRE.3='EPC3';
  737. tabe6.TITRE.4='EPC4';
  738. tabe6.TITRE.5='EPC5';
  739. tabe6.TITRE.6='EPC6';
  740. tabe6.1= 'MARQ CROI';
  741. tabe6.2= 'MARQ ROND';
  742. tabe6.3= 'MARQ CARR';
  743. tabe6.4= 'MARQ CROI';
  744. tabe6.5= 'MARQ PLUS';
  745. tabe6.6= 'MARQ CARR';
  746.  
  747. dim0=dime tab1.temps;
  748.  
  749. repe bouc1 (dim0-1);
  750.  vari1=tab1.variables_internes.(&bouc1);
  751.  repe bouc2 (dime(tabn2));
  752.   t3=exco (tabn2.(&bouc2)) vari1;
  753.   t2=extr t3 (tabn2.(&bouc2))  1 1 1;
  754.   tabl2.(&bouc2)=tabl2.(&bouc2) et (prog t2);
  755.  fin bouc2;
  756.  t0=tab1.temps.(&bouc1);
  757.  lt0=lt0 et (prog t0);
  758. fin bouc1;
  759.  
  760. tabe3=table;
  761. tabe5.'TITRE'=table;
  762. x1=(log(10))**-1;
  763. lnt0=x1*(log(lt0));
  764.  
  765. repe bouc3 (dime(tabn2));
  766.  tabe5.'TITRE'.(&bouc3)=tabn2.(&bouc3); 
  767.  si((&bouc3) ega 1);
  768.   evol1=evol manu 'temps' lt0 
  769.   tabn2.(&bouc3) tabl2.(&bouc3);
  770.   evol2=evol manu 'temps' lnt0 
  771.   tabn2.(&bouc3) tabl2.(&bouc3);
  772.  sinon;
  773.   evol1=evol1 et(evol manu 'temps' lt0 
  774.   tabn2.(&bouc3) tabl2.(&bouc3) coul tabe5.&bouc3);
  775.   evol2=evol2 et (evol manu 'temps' lnt0 
  776.   tabn2.(&bouc3) tabl2.(&bouc3) coul tabe5.&bouc3);
  777.  finsi;
  778. fin bouc3;
  779.  
  780. lrag1=prog;
  781. ldef1=prog;
  782. ltemp1=prog;
  783. repe bouc2 (dim0-1);
  784.  t3=tab1.variables_internes.(&bouc2);
  785.  aar1=extr t3 'AGEL' 1 1 1 ;
  786.  lrag1=lrag1 et (prog aar1);
  787.  ade1=extr t3 'ADEF' 1 1 1 ;
  788.  ldef1=ldef1 et (prog ade1);
  789.  ltemp1=ltemp1 et (prog (tab1.temps.(&bouc2)));
  790. fin bouc2;
  791. evar1=(evol manu 'temps' ltemp1 'aar' lrag1) coul bleu;
  792. evde1=(evol manu 'temps' ltemp1 'def' ldef1) coul roug;
  793. dess (evar1 et evde1);
  794.  
  795. dess evol1 'LEGE' tabe6;
  796. *-----------------------------------------------------------------------
  797. *-----------------------------------------------------------------------
  798. lt0=prog;
  799.  
  800. tabn2=table;
  801. tabn2.1='EPE1';
  802. tabn2.2='EPE2';
  803. tabn2.3='EPE3';
  804. tabn2.4='EPE4';
  805. tabn2.5='EPE5';
  806. tabn2.6='EPE6';
  807.  
  808. tabl2=table;
  809. tabl2.1=prog;
  810. tabl2.2=prog;
  811. tabl2.3=prog;
  812. tabl2.4=prog;
  813. tabl2.5=prog;
  814. tabl2.6=prog;
  815.  
  816. tabe5=table;
  817. tabe5.1=bleu;
  818. tabe5.2=roug;
  819. tabe5.3=rose;
  820. tabe5.4=ROUG;
  821. tabe5.5=CYAN;
  822. tabe5.6=VERT;
  823. tabe6=table;
  824. tabe6.TITRE=table;
  825. tabe6.TITRE.1='EPE1';
  826. tabe6.TITRE.2='EPE2';
  827. tabe6.TITRE.3='EPE3';
  828. tabe6.TITRE.4='EPE4';
  829. tabe6.TITRE.5='EPE5';
  830. tabe6.TITRE.6='EPE6';
  831. tabe6.1= 'MARQ CROI';
  832. tabe6.2= 'MARQ ROND';
  833. tabe6.3= 'MARQ CARR';
  834. tabe6.4= 'MARQ CROI';
  835. tabe6.5= 'MARQ PLUS';
  836. tabe6.6= 'MARQ CARR';
  837.  
  838. dim0=dime tab1.temps;
  839.  
  840. repe bouc1 (dim0-1);
  841.  vari1=tab1.variables_internes.(&bouc1);
  842.  repe bouc2 (dime(tabn2));
  843.   t3=exco (tabn2.(&bouc2)) vari1;
  844.   t2=extr t3 (tabn2.(&bouc2))  1 1 1;
  845.   tabl2.(&bouc2)=tabl2.(&bouc2) et (prog t2);
  846.  fin bouc2;
  847.  t0=tab1.temps.(&bouc1);
  848.  lt0=lt0 et (prog t0);
  849. fin bouc1;
  850.  
  851. tabe3=table;
  852. tabe5.'TITRE'=table;
  853. x1=(log(10))**-1;
  854. lnt0=x1*(log(lt0));
  855.  
  856. repe bouc3 (dime(tabn2));
  857.  tabe5.'TITRE'.(&bouc3)=tabn2.(&bouc3); 
  858.  si((&bouc3) ega 1);
  859.   evol1=evol manu 'temps' lt0 
  860.   tabn2.(&bouc3) tabl2.(&bouc3);
  861.   evol2=evol manu 'temps' lnt0 
  862.   tabn2.(&bouc3) tabl2.(&bouc3);
  863.  sinon;
  864.   evol1=evol1 et(evol manu 'temps' lt0 
  865.   tabn2.(&bouc3) tabl2.(&bouc3) coul tabe5.&bouc3);
  866.   evol2=evol2 et (evol manu 'temps' lnt0 
  867.   tabn2.(&bouc3) tabl2.(&bouc3) coul tabe5.&bouc3);
  868.  finsi;
  869. fin bouc3;
  870. dess (evol1 ) 'LEGE' tabe6;
  871.  
  872. lt0=prog;
  873.  
  874. tabn2=table;
  875. tabn2.1='EPK1';
  876. tabn2.2='EPK2';
  877. tabn2.3='EPK3';
  878. tabn2.4='EPK4';
  879. tabn2.5='EPK5';
  880. tabn2.6='EPK6';
  881.  
  882. tabl2=table;
  883. tabl2.1=prog;
  884. tabl2.2=prog;
  885. tabl2.3=prog;
  886. tabl2.4=prog;
  887. tabl2.5=prog;
  888. tabl2.6=prog;
  889.  
  890. tabe5=table;
  891. tabe5.1=bleu;
  892. tabe5.2=roug;
  893. tabe5.3=rose;
  894. tabe5.4=ROUG;
  895. tabe5.5=CYAN;
  896. tabe5.6=VERT;
  897. tabe6=table;
  898. tabe6.TITRE=table;
  899. tabe6.TITRE.1='EPK1';
  900. tabe6.TITRE.2='EPK2';
  901. tabe6.TITRE.3='EPK3';
  902. tabe6.TITRE.4='EPK4';
  903. tabe6.TITRE.5='EPK5';
  904. tabe6.TITRE.6='EPK6';
  905. tabe6.1= 'MARQ CROI';
  906. tabe6.2= 'MARQ ROND';
  907. tabe6.3= 'MARQ CARR';
  908. tabe6.4= 'MARQ CROI';
  909. tabe6.5= 'MARQ PLUS';
  910. tabe6.6= 'MARQ CARR';
  911.  
  912. dim0=dime tab1.temps;
  913.  
  914. repe bouc1 (dim0-1);
  915.  vari1=tab1.variables_internes.(&bouc1);
  916.  repe bouc2 (dime(tabn2));
  917.   t3=exco (tabn2.(&bouc2)) vari1;
  918.   t2=extr t3 (tabn2.(&bouc2))  1 1 1;
  919.   tabl2.(&bouc2)=tabl2.(&bouc2) et (prog t2);
  920.  fin bouc2;
  921.  t0=tab1.temps.(&bouc1);
  922.  lt0=lt0 et (prog t0);
  923. fin bouc1;
  924.  
  925. tabe3=table;
  926. tabe5.'TITRE'=table;
  927. x1=(log(10))**-1;
  928. lnt0=x1*(log(lt0));
  929.  
  930. repe bouc3 (dime(tabn2));
  931.  tabe5.'TITRE'.(&bouc3)=tabn2.(&bouc3); 
  932.  si((&bouc3) ega 1);
  933.   evol1=evol manu 'temps' lt0 
  934.   tabn2.(&bouc3) tabl2.(&bouc3);
  935.   evol2=evol manu 'temps' lnt0 
  936.   tabn2.(&bouc3) tabl2.(&bouc3);
  937.  sinon;
  938.   evol1=evol1 et(evol manu 'temps' lt0 
  939.   tabn2.(&bouc3) tabl2.(&bouc3) coul tabe5.&bouc3);
  940.   evol2=evol2 et (evol manu 'temps' lnt0 
  941.   tabn2.(&bouc3) tabl2.(&bouc3) coul tabe5.&bouc3);
  942.  finsi;
  943. fin bouc3;
  944.  
  945. dess (evol1 ) 'LEGE' tabe6;
  946.  
  947. *-----------------------------------------------------------------------
  948. lt0=prog;
  949.  
  950. tabn2=table;
  951. tabn2.1='EPM1';
  952. tabn2.2='EPM2';
  953. tabn2.3='EPM3';
  954. tabn2.4='EPM4';
  955. tabn2.5='EPM5';
  956. tabn2.6='EPM6';
  957.  
  958. tabl2=table;
  959. tabl2.1=prog;
  960. tabl2.2=prog;
  961. tabl2.3=prog;
  962. tabl2.4=prog;
  963. tabl2.5=prog;
  964. tabl2.6=prog;
  965.  
  966. tabe5=table;
  967. tabe5.1=bleu;
  968. tabe5.2=roug;
  969. tabe5.3=rose;
  970. tabe5.4=ROUG;
  971. tabe5.5=CYAN;
  972. tabe5.6=VERT;
  973. tabe6=table;
  974. tabe6.TITRE=table;
  975. tabe6.TITRE.1='EPM1';
  976. tabe6.TITRE.2='EPM2';
  977. tabe6.TITRE.3='EPM3';
  978. tabe6.TITRE.4='EPM4';
  979. tabe6.TITRE.5='EPM5';
  980. tabe6.TITRE.6='EPM6';
  981. tabe6.1= 'MARQ CROI';
  982. tabe6.2= 'MARQ ROND';
  983. tabe6.3= 'MARQ CARR';
  984. tabe6.4= 'MARQ CROI';
  985. tabe6.5= 'MARQ PLUS';
  986. tabe6.6= 'MARQ CARR';
  987.  
  988. dim0=dime tab1.temps;
  989.  
  990. repe bouc1 (dim0-1);
  991.  vari1=tab1.variables_internes.(&bouc1);
  992.  repe bouc2 (dime(tabn2));
  993.   t3=exco (tabn2.(&bouc2)) vari1;
  994.   t2=extr t3 (tabn2.(&bouc2))  1 1 1;
  995.   tabl2.(&bouc2)=tabl2.(&bouc2) et (prog t2);
  996.  fin bouc2;
  997.  t0=tab1.temps.(&bouc1);
  998.  lt0=lt0 et (prog t0);
  999. fin bouc1;
  1000.  
  1001. tabe3=table;
  1002. tabe5.'TITRE'=table;
  1003. x1=(log(10))**-1;
  1004. lnt0=x1*(log(lt0));
  1005.  
  1006. repe bouc3 (dime(tabn2));
  1007.  tabe5.'TITRE'.(&bouc3)=tabn2.(&bouc3); 
  1008.  si((&bouc3) ega 1);
  1009.   evol1=evol manu 'temps' lt0 
  1010.   tabn2.(&bouc3) tabl2.(&bouc3);
  1011.   evol2=evol manu 'temps' lnt0 
  1012.   tabn2.(&bouc3) tabl2.(&bouc3);
  1013.  sinon;
  1014.   evol1=evol1 et(evol manu 'temps' lt0 
  1015.   tabn2.(&bouc3) tabl2.(&bouc3) coul tabe5.&bouc3);
  1016.   evol2=evol2 et (evol manu 'temps' lnt0 
  1017.   tabn2.(&bouc3) tabl2.(&bouc3) coul tabe5.&bouc3);
  1018.  finsi;
  1019. fin bouc3;
  1020. dess (evol1 ) 'LEGE' tabe6;
  1021.  
  1022. dim0=dime tab1.temps;
  1023. trac (exco dtra tab1.variables_internes.(dim0-1)) mod1;
  1024.  
  1025. dim0=dime tab1.temps;
  1026. trac (exco dcom tab1.variables_internes.(dim0-1)) mod1;
  1027.  
  1028.  
  1029. *trac d'une variable interne
  1030. nom0='WPL0';
  1031. dim0=dime tab1.temps;
  1032. list1=prog;
  1033. list2=prog;
  1034. repe bouc4 dim0;
  1035.  vi0=tab1.variables_internes.(&bouc4-1);
  1036.  xi0=exco nom0 vi0;
  1037.  yi0=extr xi0 nom0 1 1 1;
  1038.  t0=tab1.temps.(&bouc4-1);
  1039.  list1=list1 et (prog t0);
  1040.  list2=list2 et (prog yi0);
  1041. fin bouc4;
  1042. evolpt3=evol manu 'temps' list1 nom0 list2;
  1043. dess (evolpt3 );
  1044.  
  1045. *trac d'une variable interne
  1046. nom0='DTRA';
  1047. dim0=dime tab1.temps;
  1048. list1=prog;
  1049. list2=prog;
  1050. repe bouc4 dim0;
  1051.  vi0=tab1.variables_internes.(&bouc4-1);
  1052.  xi0=exco nom0 vi0;
  1053.  yi0=extr xi0 nom0 1 1 1;
  1054.  t0=tab1.temps.(&bouc4-1);
  1055.  list1=list1 et (prog t0);
  1056.  list2=list2 et (prog yi0);
  1057. fin bouc4;
  1058. evolpt3=evol manu 'temps' list1 nom0 list2;
  1059. dess (evolpt3 );
  1060.  
  1061. *trac d'une variable interne
  1062. nom0='DTPP';
  1063. dim0=dime tab1.temps;
  1064. list1=prog;
  1065. list2=prog;
  1066. repe bouc4 dim0;
  1067.  vi0=tab1.variables_internes.(&bouc4-1);
  1068.  xi0=exco nom0 vi0;
  1069.  yi0=extr xi0 nom0 1 1 1;
  1070.  t0=tab1.temps.(&bouc4-1);
  1071.  list1=list1 et (prog t0);
  1072.  list2=list2 et (prog yi0);
  1073. fin bouc4;
  1074. evolpt3=evol manu 'temps' list1 nom0 list2;
  1075. dess (evolpt3 );
  1076.  
  1077. *trac d'une variable interne
  1078. nom0='DTM0';
  1079. dim0=dime tab1.temps;
  1080. list1=prog;
  1081. list2=prog;
  1082. repe bouc4 dim0;
  1083.  vi0=tab1.variables_internes.(&bouc4-1);
  1084.  xi0=exco nom0 vi0;
  1085.  yi0=extr xi0 nom0 1 1 1;
  1086.  t0=tab1.temps.(&bouc4-1);
  1087.  list1=list1 et (prog t0);
  1088.  list2=list2 et (prog yi0);
  1089. fin bouc4;
  1090. evolpt3=evol manu 'temps' list1 nom0 list2;
  1091. dess (evolpt3 );
  1092.  
  1093. *trac d'une variable interne
  1094. nom0='DTW0';
  1095. dim0=dime tab1.temps;
  1096. list1=prog;
  1097. list2=prog;
  1098. repe bouc4 dim0;
  1099.  vi0=tab1.variables_internes.(&bouc4-1);
  1100.  xi0=exco nom0 vi0;
  1101.  yi0=extr xi0 nom0 1 1 1;
  1102.  t0=tab1.temps.(&bouc4-1);
  1103.  list1=list1 et (prog t0);
  1104.  list2=list2 et (prog yi0);
  1105. fin bouc4;
  1106. evolpt3=evol manu 'temps' list1 nom0 list2;
  1107. dess (evolpt3 );
  1108.  
  1109. *trac d'une variable interne
  1110. nom0='DTG0';
  1111. dim0=dime tab1.temps;
  1112. list1=prog;
  1113. list2=prog;
  1114. repe bouc4 dim0;
  1115.  vi0=tab1.variables_internes.(&bouc4-1);
  1116.  xi0=exco nom0 vi0;
  1117.  yi0=extr xi0 nom0 1 1 1;
  1118.  t0=tab1.temps.(&bouc4-1);
  1119.  list1=list1 et (prog t0);
  1120.  list2=list2 et (prog yi0);
  1121. fin bouc4;
  1122. evolpt3=evol manu 'temps' list1 nom0 list2;
  1123. dess (evolpt3 );
  1124.  
  1125. *trac d'une variable interne
  1126. nom0='DTS0';
  1127. dim0=dime tab1.temps;
  1128. list1=prog;
  1129. list2=prog;
  1130. repe bouc4 dim0;
  1131.  vi0=tab1.variables_internes.(&bouc4-1);
  1132.  xi0=exco nom0 vi0;
  1133.  yi0=extr xi0 nom0 1 1 1;
  1134.  t0=tab1.temps.(&bouc4-1);
  1135.  list1=list1 et (prog t0);
  1136.  list2=list2 et (prog yi0);
  1137. fin bouc4;
  1138. evolpt3=evol manu 'temps' list1 nom0 list2;
  1139. dess (evolpt3 );
  1140.  
  1141. *trac d'une variable interne
  1142. nom0='WPL0';
  1143. dim0=dime tab1.temps;
  1144. list1=prog;
  1145. list2=prog;
  1146. repe bouc4 dim0;
  1147.  vi0=tab1.variables_internes.(&bouc4-1);
  1148.  xi0=exco nom0 vi0;
  1149.  yi0=extr xi0 nom0 1 1 1;
  1150.  t0=tab1.temps.(&bouc4-1);
  1151.  list1=list1 et (prog t0);
  1152.  list2=list2 et (prog yi0);
  1153. fin bouc4;
  1154. evolpt3=evol manu 'temps' list1 nom0 list2;
  1155. dess (evolpt3 );
  1156.  
  1157. *trac d'une variable interne
  1158. nom0='DCOM';
  1159. dim0=dime tab1.temps;
  1160. list1=prog;
  1161. list2=prog;
  1162. repe bouc4 dim0;
  1163.  vi0=tab1.variables_internes.(&bouc4-1);
  1164.  xi0=exco nom0 vi0;
  1165.  yi0=extr xi0 nom0 1 1 1;
  1166.  t0=tab1.temps.(&bouc4-1);
  1167.  list1=list1 et (prog t0);
  1168.  list2=list2 et (prog yi0);
  1169. fin bouc4;
  1170. evolpt3=evol manu 'temps' list1 nom0 list2;
  1171. dess (evolpt3 );
  1172.  
  1173. *trac d'une variable interne
  1174. nom0='DCM0';
  1175. dim0=dime tab1.temps;
  1176. list1=prog;
  1177. list2=prog;
  1178. repe bouc4 dim0;
  1179.  vi0=tab1.variables_internes.(&bouc4-1);
  1180.  xi0=exco nom0 vi0;
  1181.  yi0=extr xi0 nom0 1 1 1;
  1182.  t0=tab1.temps.(&bouc4-1);
  1183.  list1=list1 et (prog t0);
  1184.  list2=list2 et (prog yi0);
  1185. fin bouc4;
  1186. evolpt3=evol manu 'temps' list1 nom0 list2;
  1187. dess (evolpt3 );
  1188.  
  1189. *trac d'une variable interne
  1190. nom0='DCW0';
  1191. dim0=dime tab1.temps;
  1192. list1=prog;
  1193. list2=prog;
  1194. repe bouc4 dim0;
  1195.  vi0=tab1.variables_internes.(&bouc4-1);
  1196.  xi0=exco nom0 vi0;
  1197.  yi0=extr xi0 nom0 1 1 1;
  1198.  t0=tab1.temps.(&bouc4-1);
  1199.  list1=list1 et (prog t0);
  1200.  list2=list2 et (prog yi0);
  1201. fin bouc4;
  1202. evolpt3=evol manu 'temps' list1 nom0 list2;
  1203. dess (evolpt3 );
  1204.  
  1205. *trac d'une variable interne
  1206. nom0='DCG0';
  1207. dim0=dime tab1.temps;
  1208. list1=prog;
  1209. list2=prog;
  1210. repe bouc4 dim0;
  1211.  vi0=tab1.variables_internes.(&bouc4-1);
  1212.  xi0=exco nom0 vi0;
  1213.  yi0=extr xi0 nom0 1 1 1;
  1214.  t0=tab1.temps.(&bouc4-1);
  1215.  list1=list1 et (prog t0);
  1216.  list2=list2 et (prog yi0);
  1217. fin bouc4;
  1218. evolpt3=evol manu 'temps' list1 nom0 list2;
  1219. dess (evolpt3 );
  1220.  
  1221. *trac d'une variable interne
  1222. nom0='DCS0';
  1223. dim0=dime tab1.temps;
  1224. list1=prog;
  1225. list2=prog;
  1226. repe bouc4 dim0;
  1227.  vi0=tab1.variables_internes.(&bouc4-1);
  1228.  xi0=exco nom0 vi0;
  1229.  yi0=extr xi0 nom0 1 1 1;
  1230.  t0=tab1.temps.(&bouc4-1);
  1231.  list1=list1 et (prog t0);
  1232.  list2=list2 et (prog yi0);
  1233. fin bouc4;
  1234. evolpt3=evol manu 'temps' list1 nom0 list2;
  1235. dess (evolpt3 );
  1236.  
  1237. *trac d'une variable interne
  1238. nom0='DFLU';
  1239. dim0=dime tab1.temps;
  1240. list1=prog;
  1241. list2=prog;
  1242. repe bouc4 dim0;
  1243.  vi0=tab1.variables_internes.(&bouc4-1);
  1244.  xi0=exco nom0 vi0;
  1245.  yi0=extr xi0 nom0 1 1 1;
  1246.  t0=tab1.temps.(&bouc4-1);
  1247.  list1=list1 et (prog t0);
  1248.  list2=list2 et (prog yi0);
  1249. fin bouc4;
  1250. evolpt3=evol manu 'temps' list1 nom0 list2;
  1251. dess (evolpt3 );
  1252.  
  1253. *trac d'une variable interne
  1254. nom0='DTHE';
  1255. dim0=dime tab1.temps;
  1256. list1=prog;
  1257. list2=prog;
  1258. repe bouc4 dim0;
  1259.  vi0=tab1.variables_internes.(&bouc4-1);
  1260.  xi0=exco nom0 vi0;
  1261.  yi0=extr xi0 nom0 1 1 1;
  1262.  t0=tab1.temps.(&bouc4-1);
  1263.  list1=list1 et (prog t0);
  1264.  list2=list2 et (prog yi0);
  1265. fin bouc4;
  1266. evolpt3=evol manu 'temps' list1 nom0 list2;
  1267. dess (evolpt3 );
  1268.  
  1269. *trac d'une variable interne
  1270. nom0='BAFT';
  1271. dim0=dime tab1.temps;
  1272. list1=prog;
  1273. list2=prog;
  1274. repe bouc4 dim0;
  1275.  vi0=tab1.variables_internes.(&bouc4-1);
  1276.  xi0=exco nom0 vi0;
  1277.  yi0=extr xi0 nom0 1 1 1;
  1278.  t0=tab1.temps.(&bouc4-1);
  1279.  list1=list1 et (prog t0);
  1280.  list2=list2 et (prog yi0);
  1281. fin bouc4;
  1282. evolpt3=evol manu 'temps' list1 nom0 list2;
  1283. dess (evolpt3 );
  1284.  
  1285. *trac d'une variable interne
  1286. nom0='PAFT';
  1287. dim0=dime tab1.temps;
  1288. list1=prog;
  1289. list2=prog;
  1290. repe bouc4 dim0;
  1291.  vi0=tab1.variables_internes.(&bouc4-1);
  1292.  xi0=exco nom0 vi0;
  1293.  yi0=extr xi0 nom0 1 1 1;
  1294.  t0=tab1.temps.(&bouc4-1);
  1295.  list1=list1 et (prog t0);
  1296.  list2=list2 et (prog yi0);
  1297. fin bouc4;
  1298. evolpt3=evol manu 'temps' list1 nom0 list2;
  1299. dess (evolpt3 );
  1300.  
  1301. *trac d'une variable interne
  1302. nom0='BGEL';
  1303. dim0=dime tab1.temps;
  1304. list1=prog;
  1305. list2=prog;
  1306. repe bouc4 dim0;
  1307.  vi0=tab1.variables_internes.(&bouc4-1);
  1308.  xi0=exco nom0 vi0;
  1309.  yi0=extr xi0 nom0 1 1 1;
  1310.  t0=tab1.temps.(&bouc4-1);
  1311.  list1=list1 et (prog t0);
  1312.  list2=list2 et (prog yi0);
  1313. fin bouc4;
  1314. evolpt3=evol manu 'temps' list1 nom0 list2;
  1315. dess (evolpt3 );
  1316.  
  1317. *trac d'une variable interne
  1318. nom0='PGEL';
  1319. dim0=dime tab1.temps;
  1320. list1=prog;
  1321. list2=prog;
  1322. repe bouc4 dim0;
  1323.  vi0=tab1.variables_internes.(&bouc4-1);
  1324.  xi0=exco nom0 vi0;
  1325.  yi0=extr xi0 nom0 1 1 1;
  1326.  t0=tab1.temps.(&bouc4-1);
  1327.  list1=list1 et (prog t0);
  1328.  list2=list2 et (prog yi0);
  1329. fin bouc4;
  1330. evolpt3=evol manu 'temps' list1 nom0 list2;
  1331. dess (evolpt3 );
  1332.  
  1333. *-----------------------------------------------------------------------
  1334. *evol de quelques variables internes
  1335. lt0=prog;
  1336. tabn2=table;
  1337.  
  1338. tabn2.1='WPL1';
  1339. tabn2.2='WPL2';
  1340. tabn2.3='WPL3';
  1341. tabn2.4='WPL4';
  1342. tabn2.5='WPL5';
  1343. tabn2.6='WPL6';
  1344.  
  1345. tabl2=table;
  1346. tabl2.1=prog;
  1347. tabl2.2=prog;
  1348. tabl2.3=prog;
  1349. tabl2.4=prog;
  1350. tabl2.5=prog;
  1351. tabl2.6=prog;
  1352.  
  1353. tabe5=table;
  1354. tabe5.1=bleu;
  1355. tabe5.2=roug;
  1356. tabe5.3=rose;
  1357. tabe5.4=ROUG;
  1358. tabe5.5=CYAN;
  1359. tabe5.6=VERT;
  1360.  
  1361. tabe6=table;
  1362. tabe6.TITRE=table;
  1363. tabe6.TITRE.1='WPL1';
  1364. tabe6.TITRE.2='WPL2';
  1365. tabe6.TITRE.3='WPL3';
  1366. tabe6.TITRE.4='WPL4';
  1367. tabe6.TITRE.5='WPL5';
  1368. tabe6.TITRE.6='WPL6';
  1369. tabe6.1= 'MARQ CROI';
  1370. tabe6.2= 'MARQ PLUS';
  1371. tabe6.3= 'MARQ CARR';
  1372. tabe6.4= 'MARQ CROI';
  1373. tabe6.5= 'MARQ PLUS';
  1374. tabe6.6= 'MARQ CARR';
  1375.  
  1376. dim0=dime tab1.temps;
  1377.  
  1378. repe bouc1 (dim0-1);
  1379.  vari1=tab1.variables_internes.(&bouc1);
  1380.  repe bouc2 (dime(tabn2));
  1381.   t3=exco (tabn2.(&bouc2)) vari1;
  1382.   t2=extr t3 (tabn2.(&bouc2))  1 1 1;
  1383.   tabl2.(&bouc2)=tabl2.(&bouc2) et (prog t2);
  1384.  fin bouc2;
  1385.  t0=tab1.temps.(&bouc1);
  1386.  lt0=lt0 et (prog t0);
  1387. fin bouc1;
  1388.  
  1389. tabe3=table;
  1390. tabe5.'TITRE'=table;
  1391. x1=(log(10))**-1;
  1392. lnt0=x1*(log(lt0));
  1393.  
  1394. repe bouc3 (dime tabn2);
  1395.  tabe5.'TITRE'.(&bouc3)=tabn2.(&bouc3); 
  1396.  si((&bouc3) ega 1);
  1397.   evol1=evol manu 'temps' lt0 
  1398.   tabn2.(&bouc3) tabl2.(&bouc3);
  1399.   evol2=evol manu 'temps' lnt0 
  1400.   tabn2.(&bouc3) tabl2.(&bouc3);
  1401.  sinon;
  1402.   evol1=evol1 et(evol manu 'temps' lt0 
  1403.   tabn2.(&bouc3) tabl2.(&bouc3) coul tabe5.&bouc3);
  1404.   evol2=evol2 et (evol manu 'temps' lnt0 
  1405.   tabn2.(&bouc3) tabl2.(&bouc3) coul tabe5.&bouc3);
  1406.  finsi;
  1407. fin bouc3;
  1408.  
  1409. dess evol1 'LEGE' tabe6;
  1410. finsi;
  1411. fin;
  1412.  
  1413.  
  1414.  
  1415.  
  1416.  
  1417.  
  1418.  
  1419.  
  1420.  
  1421.  
  1422.  

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