Numérotation des lignes :

* fichier : pilotage_indirect_2.dgibi*+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++* Problem description: ** In the example below, we study the behaviour of a holed plate under* tensile loading. The mesh is irregular and the load is applied at * the right edge under displacement control. At this edge, only the * direction of the displacement is known a priori. The displacement * intensity is calculated in such a way to respect the path equation * specified by the user. Here, the defaut path equation (CMSI) is * utilised.**+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*'OPTI' 'ECHO' 0;'OPTI' 'DIME' 2 'ELEM' 'TRI3' 'MODE' 'PLAN CONTR';'OPTI' 'TRAC' 'PSC';*+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++* export psc figuresPLT1 = faux;* export output textEXP1 = faux;*+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++** DEFINING THE PARAMETRIC GEOMETRY**+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++pig = 3.141592654;R = 0.5;e  = (pig/80.)*R; e1 = (pig/80.)*R; CC = 1. 0.85; DTETA = ((e/R)/pig)*180.;SINDT = 'SIN'(DTETA);COSDT = 'COS'(DTETA); DTETA1 = ((e1/R)/pig)*180.;SINDT1 = 'SIN'(DTETA1);COSDT1 = 'COS'(DTETA1); dens1 = 1./5.;dens2 = 1./8.;dens3 = 2.*R*(SINDT);dens4 = 2.*R*(SINDT1); PA = 0. 0.;PB = 2. 0.;PC = 2. 2.;PD = 0. 2.; PE = (1.-(R*SINDT)) 0.;PF = (1.+(R*SINDT)) 0.;PG = (1.+(R*SINDT1)) 2.;PH = (1.-(R*SINDT1)) 2.; PI = 2. 0.85;PL = 0. 0.85; PC1 = (1.+(R*SINDT)) (0.35+(R*(1-COSDT)));PC2 = 1.5 0.85; PC3 = (1.+(R*SINDT1)) (1.35-(R*(1-COSDT1)));PC4 = (1.-(R*SINDT1)) (1.35-(R*(1-COSDT1)));PC5 = 0.5 0.85;PC6 = (1.-(R*SINDT)) (0.35+(R*(1-COSDT))); LAE = 'DROI' PA PE 'DINI' dens1 'DFIN' dens3;LC6E = 'DROI' PE PC6 'DINI' dens3 'DFIN' dens3;LC5C6 = 'CERC' PC6 CC PC5 'DINI' dens3 'DFIN' dens2;LC5L = 'DROI' PC5 PL 'DINI' dens2 'DFIN' dens1;LLA = 'DROI' PL PA 'DINI' dens1 'DFIN' dens1;LSW = LAE 'ET' LC6E 'ET' LC5C6 'ET' LC5L 'ET' LLA;SSW = 'SURF' LSW; LEF = 'DROI' 1 PE PF;LC1F = 'DROI' PF PC1 'DINI' dens3 'DFIN' dens3;LC6C1 = 'CERC' PC1 CC PC6 'DINI' dens3 'DFIN' dens3;LS = LEF 'ET' LC1F 'ET' LC6C1 'ET' ('INVE' LC6E);SS = 'SURF' LS; LFB = 'DROI' PF PB 'DINI' dens3 'DFIN' dens1;LBI = 'DROI' PB PI 'DINI' dens1 'DFIN' dens1;LC2I = 'DROI' PI PC2 'DINI' dens1 'DFIN' dens2;LC2C1 = 'CERC' PC2 CC PC1 'DINI' dens2 'DFIN' dens3;LSE = LFB 'ET' LBI 'ET' LC2I 'ET' LC2C1 'ET' ('INVE' LC1F);SSE = 'SURF' LSE; LIC = 'DROI' PI PC 'DINI' dens1 'DFIN' dens1;LCG = 'DROI' PC PG 'DINI' dens1 'DFIN' dens4;LGC3 = 'DROI' PG PC3 'DINI' dens4 'DFIN' dens4;LC3C2 = 'CERC' PC3 CC PC2 'DINI' dens4 'DFIN' dens2;LNE = ('INVE' LC2I) 'ET' LIC 'ET' LCG 'ET' LGC3 'ET' LC3C2;SNE = 'SURF' LNE; LC3C4 = 'CERC' PC3 CC PC4 'DINI' dens4 'DFIN' dens4;LGH = 'DROI' 1 PG PH;LHC4 = 'DROI' PH PC4 'DINI' dens4 'DFIN' dens4;LN = LC3C4 'ET' ('INVE' LGC3) 'ET' LGH 'ET' LHC4;SN = 'SURF' LN;  LC4C5 = 'CERC' PC5 CC PC4 'DINI' dens2 'DFIN' dens4;LHD = 'DROI' PH PD 'DINI'  dens4 'DFIN' dens1;LDL = 'DROI' PD PL 'DINI' dens1 'DFIN' dens1;LNW = ('INVE' LC5L) 'ET' LC4C5 'ET' ('INVE' LHC4) 'ET' LHD 'ET' LDL;SNW = 'SURF' LNW; DAMMPART = SS 'ET' SN;ELASPART = SSW 'ET' SSE 'ET' SNE 'ET' SNW; LAT_E = LDL 'ET' LLA;LAT_D = LIC 'ET' LBI; ESTRUT = (DAMMPART 'ET' ELASPART); 'SI' PLT1;   'TRAC' ESTRUT;'FINSI';*+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++** END OF PARAMETRIC GEOMETRY**+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++** DEFINING MODEL AND MATERIAL**+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++AT=1.0e0;AC=1.2e0;BC=1.5e+03;BT=1.0E4;epsd0=1.0E-4;YG=1.0E+9;BBT=1.06e0;NNU=0.0;RO =2.3E+03; MOD1 = 'MODE' ELASPART 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE';MAT1 = 'MATE' MOD1 'YOUN' YG 'NU' NNU 'RHO' RO; MOD2 = 'MODE' DAMMPART 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'     'ENDOMMAGEMENT' 'MAZARS';MAT2 = 'MATE' MOD2 'YOUN' YG 'NU' NNU 'RHO' RO 'BTRA' BT     'KTR0' epsd0 'ATRA' AT 'ACOM' AC 'BCOM' BC 'BETA' BBT; MODTOT = MOD1 'ET' MOD2;MATTOT = MAT1 'ET' MAT2; *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++** END OF MODEL AND MATERIAL**+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++** DEFINING BOUNDARY CONDITIONS (known a priori)**+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++CL1 = 'BLOQ' 'UX' LAT_E;CL2 = 'BLOQ' 'UY' PA;CLT = CL1 'ET' CL2; LA = ESTRUT 'POIN' 'DROIT' PA PB 0.00001;CL4 = 'BLOQ' 'UY' LA;CLT = CLT 'ET' CL4; *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++** END OF BOUNDARY CONDITIONS**+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++** PROBLEM SOLUTION**+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++* * number of time incrementsntincrem = 100;** control parameterTAU = 1.0E-4;** direction of the applied displacement (u_hat)CL3 = 'BLOQ' 'UX' LAT_D;CLT = CLT 'ET' CL3;uch = 'DEPI' CL3 1.0 ;*lis_incr= 'PROG' 0. 'PAS' 1    (ntincrem. * 1) ;lis_tau = 'PROG' 0. 'PAS' TAU  (ntincrem. * TAU) ;ev1 = 'EVOL' 'MANU' 't' lis_incr 'tau' lis_tau ;*****  preparing PASAPAS input tabletab1 = 'TABLE';tab1 . 'MOVA' = 'D' ;tab1 . 'BLOCAGES_MECANIQUES' = CLT ;tab1 . 'CARACTERISTIQUES' = mattot ;tab1 . 'MODELE' = modtot ;tab1 . 'TEMPS_CALCULES' = lis_incr ; ****  lines to be included for ativating path-following control ******tab1 . 'PILOTAGE_INDIRECT'     = VRAI;tab1 . 'DEPLACEMENTS_PILOTES'  = uch ;tab1 . 'PARAMETRE_DE_PILOTAGE' = ev1 ;tab1 . 'PRECISION' = 1.0e-7; PASAPAS tab1; *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++** END OF SOLUTION**+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++** BEGIN OF POST-PROCESSING**+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*****  procedure @TOTAL adapted'DEBPROC' RESFORCE CH1*'CHPOINT' GEO*'MAILLAGE' COMP1*'MOT' ;  CH2 = 'REDU' CH1 GEO ;  CHTOT = 'RESU' CH2 ;  P1 = ('EXTR' CHTOT MAIL) 'POIN' 1 ;                                 TOT1 = 'EXTR' CHTOT COMP1 P1 ;'FINPROC' TOT1 ;**nn = 'DIME' tab1 . 'TEMPS';*PROF = 'PROG' 0.;PROT = 'PROG' 0.;PROD = 'PROG' 0.;*'REPE' LAB (nn - 1 ) ;  ii = &LAB ;  TT     = tab1 . 'TEMPS' . ii ;  deplac = tab1 . 'DEPLACEMENTS' . ii ;  force_indirect = 'REAC' CL1 deplac  ;  FF = RESFORCE force_indirect LAT_E 'FX' ;  DD = 'EXTR' deplac 'UX' PB ;   PROT = PROT 'ET' ( 'PROG' TT) ;  PROF = PROF 'ET' ( 'PROG' (-1.*FF)) ;  PROD = PROD 'ET' ( 'PROG' DD) ;'FIN' LAB ; 'SI' PLT1;   PROETA = TAB1 . 'COEFFICIENT_DE_PILOTAGE' ;  evo11 = 'EVOL' 'TURQ' 'MANU' 'disp' PROD 'forc' PROF ;   evo12 = 'EVOL' 'ROUG' 'MANU' 'increment' PROT 'eta' PROETA ;   'DESS' evo11 ;  'DESS' evo12 ;'FINSI'; 'SI' EXP1;  @EXCEL1 evo11 'forc_displ_cmsi.txt' ;'FINSI'; *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++** END OF POST-PROCESSING* *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++** BEGIN OF ERROR TESTING* *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ refx = 'PROG' 0.;refx = refx 'ET' ('PROG' 7.25234E-05);refx = refx 'ET' ('PROG' 9.52043E-05);refx = refx 'ET' ('PROG' 8.90173E-05);refx = refx 'ET' ('PROG' 8.23523E-05);refx = refx 'ET' ('PROG' 7.52624E-05);refx = refx 'ET' ('PROG' 6.99765E-05);refx = refx 'ET' ('PROG' 6.64275E-05);refx = refx 'ET' ('PROG' 6.47229E-05);refx = refx 'ET' ('PROG' 6.31895E-05);refx = refx 'ET' ('PROG' 6.25926E-05);refx = refx 'ET' ('PROG' 6.19542E-05);refx = refx 'ET' ('PROG' 6.13881E-05);refx = refx 'ET' ('PROG' 6.11670E-05);refx = refx 'ET' ('PROG' 6.06850E-05);refx = refx 'ET' ('PROG' 6.06745E-05);refx = refx 'ET' ('PROG' 6.05893E-05);refx = refx 'ET' ('PROG' 6.04648E-05);refx = refx 'ET' ('PROG' 6.02701E-05);refx = refx 'ET' ('PROG' 5.98903E-05);refx = refx 'ET' ('PROG' 6.06235E-05);refx = refx 'ET' ('PROG' 6.10159E-05);refx = refx 'ET' ('PROG' 6.09750E-05);refx = refx 'ET' ('PROG' 6.15300E-05);refx = refx 'ET' ('PROG' 6.13401E-05);refx = refx 'ET' ('PROG' 6.16858E-05);refx = refx 'ET' ('PROG' 6.21556E-05);refx = refx 'ET' ('PROG' 6.19956E-05);refx = refx 'ET' ('PROG' 6.23896E-05);refx = refx 'ET' ('PROG' 6.21802E-05);refx = refx 'ET' ('PROG' 6.23712E-05);refx = refx 'ET' ('PROG' 6.20276E-05);refx = refx 'ET' ('PROG' 6.14528E-05);refx = refx 'ET' ('PROG' 6.25645E-05);refx = refx 'ET' ('PROG' 6.24348E-05);refx = refx 'ET' ('PROG' 6.21010E-05);refx = refx 'ET' ('PROG' 6.28804E-05);refx = refx 'ET' ('PROG' 6.19484E-05);refx = refx 'ET' ('PROG' 6.18783E-05);refx = refx 'ET' ('PROG' 6.12946E-05);refx = refx 'ET' ('PROG' 5.94620E-05);refx = refx 'ET' ('PROG' 5.94464E-05);refx = refx 'ET' ('PROG' 5.48346E-05);refx = refx 'ET' ('PROG' 5.00176E-05);refx = refx 'ET' ('PROG' 3.55282E-05);refx = refx 'ET' ('PROG' 3.02823E-05);refx = refx 'ET' ('PROG' 2.82446E-05);refx = refx 'ET' ('PROG' 2.79204E-05);refx = refx 'ET' ('PROG' 2.95119E-05);refx = refx 'ET' ('PROG' 3.21104E-05);refx = refx 'ET' ('PROG' 3.53427E-05);refx = refx 'ET' ('PROG' 3.89299E-05);refx = refx 'ET' ('PROG' 4.26923E-05);refx = refx 'ET' ('PROG' 4.65366E-05);refx = refx 'ET' ('PROG' 5.04177E-05);refx = refx 'ET' ('PROG' 5.43149E-05);refx = refx 'ET' ('PROG' 5.82192E-05);refx = refx 'ET' ('PROG' 6.21265E-05);refx = refx 'ET' ('PROG' 6.60349E-05);refx = refx 'ET' ('PROG' 6.99438E-05);refx = refx 'ET' ('PROG' 7.38529E-05);refx = refx 'ET' ('PROG' 7.77621E-05);refx = refx 'ET' ('PROG' 8.16713E-05);refx = refx 'ET' ('PROG' 8.55805E-05);refx = refx 'ET' ('PROG' 8.94897E-05);refx = refx 'ET' ('PROG' 9.33989E-05);refx = refx 'ET' ('PROG' 9.73081E-05);refx = refx 'ET' ('PROG' 1.01217E-04);refx = refx 'ET' ('PROG' 1.05127E-04);refx = refx 'ET' ('PROG' 1.01457E-04);refx = refx 'ET' ('PROG' 9.27572E-05);refx = refx 'ET' ('PROG' 8.45790E-05);refx = refx 'ET' ('PROG' 7.88846E-05);refx = refx 'ET' ('PROG' 7.51973E-05);refx = refx 'ET' ('PROG' 7.26282E-05);refx = refx 'ET' ('PROG' 7.03561E-05);refx = refx 'ET' ('PROG' 6.88321E-05);refx = refx 'ET' ('PROG' 6.75941E-05);refx = refx 'ET' ('PROG' 6.65391E-05);refx = refx 'ET' ('PROG' 6.60840E-05);refx = refx 'ET' ('PROG' 6.49857E-05);refx = refx 'ET' ('PROG' 6.34483E-05);refx = refx 'ET' ('PROG' 6.19822E-05);refx = refx 'ET' ('PROG' 6.05333E-05);refx = refx 'ET' ('PROG' 5.77906E-05);refx = refx 'ET' ('PROG' 5.31490E-05);refx = refx 'ET' ('PROG' 4.25939E-05);refx = refx 'ET' ('PROG' 3.43240E-05);refx = refx 'ET' ('PROG' 2.93451E-05);refx = refx 'ET' ('PROG' 2.74954E-05);refx = refx 'ET' ('PROG' 2.81269E-05);refx = refx 'ET' ('PROG' 3.02495E-05);refx = refx 'ET' ('PROG' 3.32799E-05);refx = refx 'ET' ('PROG' 3.67737E-05);refx = refx 'ET' ('PROG' 4.04981E-05);refx = refx 'ET' ('PROG' 4.43358E-05);refx = refx 'ET' ('PROG' 4.82257E-05);refx = refx 'ET' ('PROG' 5.21371E-05);refx = refx 'ET' ('PROG' 5.60573E-05);refx = refx 'ET' ('PROG' 5.99810E-05);refx = refx 'ET' ('PROG' 6.39061E-05);**  suite a modification de mazar pour mettre un endommagement residuel nouvelle liste de valeur*REFX = 'PROG' 0.000000    7.25234E-05  9.52044E-05  8.90173E-05  8.23523E-05  7.52625E-05  6.99765E-05  6.64275E-05  6.47228E-05;REFX = REFX 'ET' ('PROG'  6.31895E-05  6.25926E-05  6.19542E-05  6.13881E-05  6.11671E-05  6.06880E-05  6.06790E-05  6.05976E-05);REFX = REFX 'ET' ('PROG'  6.04768E-05  6.02842E-05  5.99090E-05  6.06467E-05  6.10430E-05  6.10057E-05  6.15633E-05  6.13763E-05);REFX = REFX 'ET' ('PROG'  6.17253E-05  6.21982E-05  6.20421E-05  6.24373E-05  6.22293E-05  6.24231E-05  6.20775E-05  6.15096E-05);REFX = REFX 'ET' ('PROG'  6.26263E-05  6.24904E-05  6.21702E-05  6.29478E-05  6.20093E-05  6.19512E-05  6.13524E-05  5.95224E-05);REFX = REFX 'ET' ('PROG'  5.95014E-05  5.48704E-05  5.00570E-05  4.35154E-05  3.60484E-05  3.07921E-05  2.85441E-05  2.81884E-05);REFX = REFX 'ET' ('PROG'  2.95751E-05  3.21451E-05  3.53807E-05  3.89711E-05  4.28356E-05  4.67677E-05  5.06999E-05  5.46320E-05);REFX = REFX 'ET' ('PROG'  5.85641E-05  6.24962E-05  6.64283E-05  7.03604E-05  7.42926E-05  7.82247E-05  8.21568E-05  8.60889E-05);REFX = REFX 'ET' ('PROG'  9.00210E-05  9.39531E-05  9.78852E-05  1.01817E-04  1.05750E-04  1.04886E-04  9.61616E-05  8.74241E-05);REFX = REFX 'ET' ('PROG'  8.11869E-05  7.66412E-05  7.36086E-05  7.09466E-05  6.91937E-05  6.81026E-05  6.69824E-05  6.62705E-05);REFX = REFX 'ET' ('PROG'  6.56305E-05  6.39545E-05  6.23256E-05  6.09259E-05  5.86086E-05  5.49990E-05  4.57171E-05  3.65636E-05);REFX = REFX 'ET' ('PROG'  3.05317E-05  2.78188E-05  2.77164E-05  2.94886E-05  3.23281E-05  3.57223E-05  3.93977E-05  4.33015E-05);REFX = REFX 'ET' ('PROG'  4.72388E-05  5.11762E-05  5.51135E-05  5.90509E-05);  refy = 'PROG' 0.;refy = refy 'ET' ('PROG' 4.45463E+04);refy = refy 'ET' ('PROG' 5.83719E+04);refy = refy 'ET' ('PROG' 5.43956E+04);refy = refy 'ET' ('PROG' 5.00694E+04);refy = refy 'ET' ('PROG' 4.53975E+04);refy = refy 'ET' ('PROG' 4.17647E+04);refy = refy 'ET' ('PROG' 3.92294E+04);refy = refy 'ET' ('PROG' 3.77753E+04);refy = refy 'ET' ('PROG' 3.63999E+04);refy = refy 'ET' ('PROG' 3.56303E+04);refy = refy 'ET' ('PROG' 3.47410E+04);refy = refy 'ET' ('PROG' 3.39867E+04);refy = refy 'ET' ('PROG' 3.33804E+04);refy = refy 'ET' ('PROG' 3.25917E+04);refy = refy 'ET' ('PROG' 3.21308E+04);refy = refy 'ET' ('PROG' 3.15410E+04);refy = refy 'ET' ('PROG' 3.09996E+04);refy = refy 'ET' ('PROG' 3.03773E+04);refy = refy 'ET' ('PROG' 2.96708E+04);refy = refy 'ET' ('PROG' 2.95883E+04);refy = refy 'ET' ('PROG' 2.92486E+04);refy = refy 'ET' ('PROG' 2.87976E+04);refy = refy 'ET' ('PROG' 2.86039E+04);refy = refy 'ET' ('PROG' 2.79892E+04);refy = refy 'ET' ('PROG' 2.77311E+04);refy = refy 'ET' ('PROG' 2.74429E+04);refy = refy 'ET' ('PROG' 2.69136E+04);refy = refy 'ET' ('PROG' 2.66394E+04);refy = refy 'ET' ('PROG' 2.60137E+04);refy = refy 'ET' ('PROG' 2.56496E+04);refy = refy 'ET' ('PROG' 2.49869E+04);refy = refy 'ET' ('PROG' 2.42949E+04);refy = refy 'ET' ('PROG' 2.43173E+04);refy = refy 'ET' ('PROG' 2.37479E+04);refy = refy 'ET' ('PROG' 2.32228E+04);refy = refy 'ET' ('PROG' 2.30970E+04);refy = refy 'ET' ('PROG' 2.22625E+04);refy = refy 'ET' ('PROG' 2.18520E+04);refy = refy 'ET' ('PROG' 2.11716E+04);refy = refy 'ET' ('PROG' 2.00706E+04);refy = refy 'ET' ('PROG' 1.96542E+04);refy = refy 'ET' ('PROG' 1.76046E+04);refy = refy 'ET' ('PROG' 1.56278E+04);refy = refy 'ET' ('PROG' 1.03518E+04);refy = refy 'ET' ('PROG' 8.46019E+03);refy = refy 'ET' ('PROG' 7.54396E+03);refy = refy 'ET' ('PROG' 7.11231E+03);refy = refy 'ET' ('PROG' 7.28604E+03);refy = refy 'ET' ('PROG' 7.77220E+03);refy = refy 'ET' ('PROG' 8.46328E+03);refy = refy 'ET' ('PROG' 9.27464E+03);refy = refy 'ET' ('PROG' 1.01473E+04);refy = refy 'ET' ('PROG' 1.10496E+04);refy = refy 'ET' ('PROG' 1.19656E+04);refy = refy 'ET' ('PROG' 1.28878E+04);refy = refy 'ET' ('PROG' 1.38130E+04);refy = refy 'ET' ('PROG' 1.47395E+04);refy = refy 'ET' ('PROG' 1.56665E+04);refy = refy 'ET' ('PROG' 1.65938E+04);refy = refy 'ET' ('PROG' 1.75211E+04);refy = refy 'ET' ('PROG' 1.84486E+04);refy = refy 'ET' ('PROG' 1.93760E+04);refy = refy 'ET' ('PROG' 2.03034E+04);refy = refy 'ET' ('PROG' 2.12308E+04);refy = refy 'ET' ('PROG' 2.21583E+04);refy = refy 'ET' ('PROG' 2.30857E+04);refy = refy 'ET' ('PROG' 2.40088E+04);refy = refy 'ET' ('PROG' 2.49235E+04);refy = refy 'ET' ('PROG' 2.38237E+04);refy = refy 'ET' ('PROG' 2.15515E+04);refy = refy 'ET' ('PROG' 1.93316E+04);refy = refy 'ET' ('PROG' 1.76276E+04);refy = refy 'ET' ('PROG' 1.64231E+04);refy = refy 'ET' ('PROG' 1.54226E+04);refy = refy 'ET' ('PROG' 1.45048E+04);refy = refy 'ET' ('PROG' 1.37513E+04);refy = refy 'ET' ('PROG' 1.29944E+04);refy = refy 'ET' ('PROG' 1.23544E+04);refy = refy 'ET' ('PROG' 1.17698E+04);refy = refy 'ET' ('PROG' 1.10863E+04);refy = refy 'ET' ('PROG' 1.03002E+04);refy = refy 'ET' ('PROG' 9.53999E+03);refy = refy 'ET' ('PROG' 8.76775E+03);refy = refy 'ET' ('PROG' 7.77889E+03);refy = refy 'ET' ('PROG' 6.50887E+03);refy = refy 'ET' ('PROG' 4.50567E+03);refy = refy 'ET' ('PROG' 2.98532E+03);refy = refy 'ET' ('PROG' 1.85411E+03);refy = refy 'ET' ('PROG' 1.07094E+03);refy = refy 'ET' ('PROG' 5.96771E+02);refy = refy 'ET' ('PROG' 3.10709E+02);refy = refy 'ET' ('PROG' 1.56888E+02);refy = refy 'ET' ('PROG' 7.53307E+01);refy = refy 'ET' ('PROG' 3.32747E+01);refy = refy 'ET' ('PROG' 1.35795E+01);refy = refy 'ET' ('PROG' 5.44144E+00);refy = refy 'ET' ('PROG' 2.16530E+00);refy = refy 'ET' ('PROG' 8.56639E-01);refy = refy 'ET' ('PROG' 3.37225E-01);refy = refy 'ET' ('PROG' 1.32182E-01);REFY = 'PROG'   0.0000       44546.       58372.       54396.       50069.       45398.       41765.       39229.       37775.;REFY = REFY 'ET' ('PROG'     36400.       35630.       34741.       33987.       33381.       32594.       32134.       31547.);REFY = REFY 'ET' ('PROG'     31009.       30388.       29685.       29606.       29270.       28822.       28630.       28019.);REFY = REFY 'ET' ('PROG'     27763.       27478.       26952.       26679.       26056.       25694.       25031.       24343.);REFY = REFY 'ET' ('PROG'     24370.       23799.       23281.       23156.       22321.       21916.       21230.       20130.);REFY = REFY 'ET' ('PROG'     19713.       17655.       15677.       13114.       10574.       8673.6       7680.6       7238.2);REFY = REFY 'ET' ('PROG'     7333.2       7809.2       8504.0       9319.4       10233.       11173.       12112.       13052.);REFY = REFY 'ET' ('PROG'     13991.       14930.       15870.       16809.       17748.       18688.       19627.       20567.);REFY = REFY 'ET' ('PROG'     21506.       22445.       23385.       24318.       25244.       24883.       22607.       20270.);REFY = REFY 'ET' ('PROG'     18464.       17029.       15949.       14884.       14098.       13383.       12723.       12117.);REFY = REFY 'ET' ('PROG'     11474.       10675.       9827.5       9087.1       8125.5       7045.6       5140.8       3478.0);REFY = REFY 'ET' ('PROG'     2212.9       1326.5       748.85       419.76       243.85       152.52       108.16       109.91);REFY = REFY 'ET' ('PROG'     119.91       129.90       139.89       149.89);   tol1 = 1.0e-2;maxres = 0.;'REPE' BTEST (nn - 1 ) ; V1 = 'EXTR' refx (&BTEST + 1); V2 = 'EXTR' prod (&BTEST + 1); res1 = 'ABS' ((V1 - V2)/V1); 'SI' (res1 > maxres); maxres = res1; 'FINSI';  W1 = 'EXTR' refy (&BTEST + 1); W2 = 'EXTR' prof (&BTEST + 1); res2 = 'ABS' ((W1 - W2)/W1); 'SI' (res2 > maxres); maxres = res2; 'FINSI'; 'FIN' BTEST ; err1 = maxres > tol1;  'SI' err1; 'ERRE' 5 ;'FINSI'; *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++** END OF ERROR TESTING* *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++'FIN';

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