Télécharger pilotage_indirect_2.dgibi

Retour à la liste

Numérotation des lignes :

  1. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  2. * Problem description:
  3. *
  4. * In the example below, we study the behaviour of a holed plate under
  5. * tensile loading. The mesh is irregular and the load is applied at
  6. * the right edge under displacement control. At this edge, only the
  7. * direction of the displacement is known a priori. The displacement
  8. * intensity is calculated in such a way to respect the path equation
  9. * specified by the user. Here, the defaut path equation (CMSI) is
  10. * utilised.
  11. *
  12. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  13. 'OPTI' 'ECHO' 0;
  14. 'OPTI' 'DIME' 2 'ELEM' 'TRI3' 'MODE' 'PLAN CONTR';
  15. 'OPTI' 'TRAC' 'PSC';
  16. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  17. * export psc figures
  18. PLT1 = VRAI;
  19. * export output text
  20. EXP1 = FAUX;
  21. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  22. *
  23. * DEFINING THE PARAMETRIC GEOMETRY
  24. *
  25. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  26. pig = 3.141592654;
  27. R = 0.5;
  28. e = (pig/80.)*R;
  29. e1 = (pig/80.)*R;
  30. CC = 1. 0.85;
  31.  
  32. DTETA = ((e/R)/pig)*180.;
  33. SINDT = 'SIN'(DTETA);
  34. COSDT = 'COS'(DTETA);
  35.  
  36. DTETA1 = ((e1/R)/pig)*180.;
  37. SINDT1 = 'SIN'(DTETA1);
  38. COSDT1 = 'COS'(DTETA1);
  39.  
  40. dens1 = 1./5.;
  41. dens2 = 1./8.;
  42. dens3 = 2.*R*(SINDT);
  43. dens4 = 2.*R*(SINDT1);
  44.  
  45. PA = 0. 0.;
  46. PB = 2. 0.;
  47. PC = 2. 2.;
  48. PD = 0. 2.;
  49.  
  50. PE = (1.-(R*SINDT)) 0.;
  51. PF = (1.+(R*SINDT)) 0.;
  52. PG = (1.+(R*SINDT1)) 2.;
  53. PH = (1.-(R*SINDT1)) 2.;
  54.  
  55. PI = 2. 0.85;
  56. PL = 0. 0.85;
  57.  
  58. PC1 = (1.+(R*SINDT)) (0.35+(R*(1-COSDT)));
  59. PC2 = 1.5 0.85;
  60. PC3 = (1.+(R*SINDT1)) (1.35-(R*(1-COSDT1)));
  61. PC4 = (1.-(R*SINDT1)) (1.35-(R*(1-COSDT1)));
  62. PC5 = 0.5 0.85;
  63. PC6 = (1.-(R*SINDT)) (0.35+(R*(1-COSDT)));
  64.  
  65. LAE = 'DROI' PA PE 'DINI' dens1 'DFIN' dens3;
  66. LC6E = 'DROI' PE PC6 'DINI' dens3 'DFIN' dens3;
  67. LC5C6 = 'CERC' PC6 CC PC5 'DINI' dens3 'DFIN' dens2;
  68. LC5L = 'DROI' PC5 PL 'DINI' dens2 'DFIN' dens1;
  69. LLA = 'DROI' PL PA 'DINI' dens1 'DFIN' dens1;
  70. LSW = LAE 'ET' LC6E 'ET' LC5C6 'ET' LC5L 'ET' LLA;
  71. SSW = 'SURF' LSW;
  72.  
  73. LEF = 'DROI' 1 PE PF;
  74. LC1F = 'DROI' PF PC1 'DINI' dens3 'DFIN' dens3;
  75. LC6C1 = 'CERC' PC1 CC PC6 'DINI' dens3 'DFIN' dens3;
  76. LS = LEF 'ET' LC1F 'ET' LC6C1 'ET' ('INVE' LC6E);
  77. SS = 'SURF' LS;
  78.  
  79. LFB = 'DROI' PF PB 'DINI' dens3 'DFIN' dens1;
  80. LBI = 'DROI' PB PI 'DINI' dens1 'DFIN' dens1;
  81. LC2I = 'DROI' PI PC2 'DINI' dens1 'DFIN' dens2;
  82. LC2C1 = 'CERC' PC2 CC PC1 'DINI' dens2 'DFIN' dens3;
  83. LSE = LFB 'ET' LBI 'ET' LC2I 'ET' LC2C1 'ET' ('INVE' LC1F);
  84. SSE = 'SURF' LSE;
  85.  
  86. LIC = 'DROI' PI PC 'DINI' dens1 'DFIN' dens1;
  87. LCG = 'DROI' PC PG 'DINI' dens1 'DFIN' dens4;
  88. LGC3 = 'DROI' PG PC3 'DINI' dens4 'DFIN' dens4;
  89. LC3C2 = 'CERC' PC3 CC PC2 'DINI' dens4 'DFIN' dens2;
  90. LNE = ('INVE' LC2I) 'ET' LIC 'ET' LCG 'ET' LGC3 'ET' LC3C2;
  91. SNE = 'SURF' LNE;
  92.  
  93. LC3C4 = 'CERC' PC3 CC PC4 'DINI' dens4 'DFIN' dens4;
  94. LGH = 'DROI' 1 PG PH;
  95. LHC4 = 'DROI' PH PC4 'DINI' dens4 'DFIN' dens4;
  96. LN = LC3C4 'ET' ('INVE' LGC3) 'ET' LGH 'ET' LHC4;
  97. SN = 'SURF' LN;
  98.  
  99. LC4C5 = 'CERC' PC5 CC PC4 'DINI' dens2 'DFIN' dens4;
  100. LHD = 'DROI' PH PD 'DINI' dens4 'DFIN' dens1;
  101. LDL = 'DROI' PD PL 'DINI' dens1 'DFIN' dens1;
  102. LNW = ('INVE' LC5L) 'ET' LC4C5 'ET' ('INVE' LHC4) 'ET' LHD 'ET' LDL;
  103. SNW = 'SURF' LNW;
  104.  
  105. DAMMPART = SS 'ET' SN;
  106. ELASPART = SSW 'ET' SSE 'ET' SNE 'ET' SNW;
  107.  
  108. LAT_E = LDL 'ET' LLA;
  109. LAT_D = LIC 'ET' LBI;
  110.  
  111. ESTRUT = (DAMMPART 'ET' ELASPART);
  112.  
  113. 'SI' PLT1;
  114. 'TRAC' ESTRUT;
  115. 'FINSI';
  116. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  117. *
  118. * END OF PARAMETRIC GEOMETRY
  119. *
  120. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  121.  
  122. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  123. *
  124. * DEFINING MODEL AND MATERIAL
  125. *
  126. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  127. AT=1.0e0;
  128. AC=1.2e0;
  129. BC=1.5e+03;
  130. BT=1.0E4;
  131. epsd0=1.0E-4;
  132. YG=1.0E+9;
  133. BBT=1.06e0;
  134. NNU=0.0;
  135. RO =2.3E+03;
  136.  
  137. MOD1 = 'MODE' ELASPART 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE';
  138. MAT1 = 'MATE' MOD1 'YOUN' YG 'NU' NNU 'RHO' RO;
  139.  
  140. MOD2 = 'MODE' DAMMPART 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  141. 'ENDOMMAGEMENT' 'MAZARS';
  142. MAT2 = 'MATE' MOD2 'YOUN' YG 'NU' NNU 'RHO' RO 'BTRA' BT
  143. 'KTR0' epsd0 'ATRA' AT 'ACOM' AC 'BCOM' BC 'BETA' BBT;
  144.  
  145. MODTOT = MOD1 'ET' MOD2;
  146. MATTOT = MAT1 'ET' MAT2;
  147.  
  148. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  149. *
  150. * END OF MODEL AND MATERIAL
  151. *
  152. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  153.  
  154. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  155. *
  156. * DEFINING BOUNDARY CONDITIONS (known a priori)
  157. *
  158. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  159. CL1 = 'BLOQ' 'UX' LAT_E;
  160. CL2 = 'BLOQ' 'UY' PA;
  161. CLT = CL1 'ET' CL2;
  162.  
  163. LA = ESTRUT 'POIN' 'DROIT' PA PB 0.00001;
  164. CL4 = 'BLOQ' 'UY' LA;
  165. CLT = CLT 'ET' CL4;
  166.  
  167. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  168. *
  169. * END OF BOUNDARY CONDITIONS
  170. *
  171. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  172.  
  173. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  174. *
  175. * PROBLEM SOLUTION
  176. *
  177. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  178. *
  179. * number of time increments
  180. ntincrem = 100;
  181. *
  182. * control parameter
  183. TAU = 1.0E-4;
  184. *
  185. * direction of the applied displacement (u_hat)
  186. CL3 = 'BLOQ' 'UX' LAT_D;
  187. CLT = CLT 'ET' CL3;
  188. uch = 'DEPI' CL3 1.0 ;
  189. *
  190. lis_incr= 'PROG' 0. 'PAS' 1 (ntincrem. * 1) ;
  191. lis_tau = 'PROG' 0. 'PAS' TAU (ntincrem. * TAU) ;
  192. ev1 = 'EVOL' 'MANU' 't' lis_incr 'tau' lis_tau ;
  193. *
  194. **** preparing PASAPAS input table
  195. tab1 = 'TABLE';
  196. tab1 . 'MOVA' = 'D' ;
  197. tab1 . 'BLOCAGES_MECANIQUES' = CLT ;
  198. tab1 . 'CARACTERISTIQUES' = mattot ;
  199. tab1 . 'MODELE' = modtot ;
  200. tab1 . 'TEMPS_CALCULES' = lis_incr ;
  201.  
  202. **** lines to be included for ativating path-following control ******
  203. tab1 . 'PILOTAGE_INDIRECT' = VRAI;
  204. tab1 . 'DEPLACEMENTS_PILOTES' = uch ;
  205. tab1 . 'PARAMETRE_DE_PILOTAGE' = ev1 ;
  206.  
  207. PASAPAS tab1;
  208.  
  209. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  210. *
  211. * END OF SOLUTION
  212. *
  213. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  214.  
  215. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  216. *
  217. * BEGIN OF POST-PROCESSING
  218. *
  219. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  220. *
  221. **** procedure @TOTAL adapted
  222. 'DEBPROC' RESFORCE CH1*'CHPOINT' GEO*'MAILLAGE' COMP1*'MOT' ;
  223. CH2 = 'REDU' CH1 GEO ;
  224. CHTOT = 'RESU' CH2 ;
  225. P1 = ('EXTR' CHTOT MAIL) 'POIN' 1 ;
  226. TOT1 = 'EXTR' CHTOT COMP1 P1 ;
  227. 'FINPROC' TOT1 ;
  228. *
  229. *
  230. nn = 'DIME' tab1 . 'TEMPS';
  231. *
  232. PROF = 'PROG' 0.;
  233. PROT = 'PROG' 0.;
  234. PROD = 'PROG' 0.;
  235. *
  236. 'REPE' LAB (nn - 1 ) ;
  237. ii = &LAB ;
  238. TT = tab1 . 'TEMPS' . ii ;
  239. deplac = tab1 . 'DEPLACEMENTS' . ii ;
  240. force_indirect = 'REAC' CL1 deplac ;
  241. FF = RESFORCE force_indirect LAT_E 'FX' ;
  242. DD = 'EXTR' deplac 'UX' PB ;
  243.  
  244. PROT = PROT 'ET' ( 'PROG' TT) ;
  245. PROF = PROF 'ET' ( 'PROG' (-1.*FF)) ;
  246. PROD = PROD 'ET' ( 'PROG' DD) ;
  247. 'FIN' LAB ;
  248.  
  249. 'SI' PLT1;
  250. PROETA = TAB1 . 'COEFFICIENT_DE_PILOTAGE' ;
  251. evo11 = 'EVOL' 'TURQ' 'MANU' 'disp' PROD 'forc' PROF ;
  252. evo12 = 'EVOL' 'ROUG' 'MANU' 'increment' PROT 'eta' PROETA ;
  253. 'DESS' evo11 ;
  254. 'DESS' evo12 ;
  255. 'FINSI';
  256.  
  257. 'SI' EXP1;
  258. @EXCEL1 evo11 'forc_displ_cmsi.txt' ;
  259. 'FINSI';
  260.  
  261. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  262. *
  263. * END OF POST-PROCESSING
  264. *
  265. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  266.  
  267. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  268. *
  269. * BEGIN OF ERROR TESTING
  270. *
  271. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  272.  
  273. refx = 'PROG' 0.;
  274. refx = refx 'ET' ('PROG' 7.25234E-05);
  275. refx = refx 'ET' ('PROG' 9.52176E-05);
  276. refx = refx 'ET' ('PROG' 8.90416E-05);
  277. refx = refx 'ET' ('PROG' 8.23551E-05);
  278. refx = refx 'ET' ('PROG' 7.54260E-05);
  279. refx = refx 'ET' ('PROG' 6.99694E-05);
  280. refx = refx 'ET' ('PROG' 6.64493E-05);
  281. refx = refx 'ET' ('PROG' 6.47405E-05);
  282. refx = refx 'ET' ('PROG' 6.31880E-05);
  283. refx = refx 'ET' ('PROG' 6.25887E-05);
  284. refx = refx 'ET' ('PROG' 6.19514E-05);
  285. refx = refx 'ET' ('PROG' 6.14393E-05);
  286. refx = refx 'ET' ('PROG' 6.11918E-05);
  287. refx = refx 'ET' ('PROG' 6.08784E-05);
  288. refx = refx 'ET' ('PROG' 6.06724E-05);
  289. refx = refx 'ET' ('PROG' 6.05939E-05);
  290. refx = refx 'ET' ('PROG' 6.05941E-05);
  291. refx = refx 'ET' ('PROG' 6.03417E-05);
  292. refx = refx 'ET' ('PROG' 5.98811E-05);
  293. refx = refx 'ET' ('PROG' 6.05896E-05);
  294. refx = refx 'ET' ('PROG' 6.10311E-05);
  295. refx = refx 'ET' ('PROG' 6.09359E-05);
  296. refx = refx 'ET' ('PROG' 6.15372E-05);
  297. refx = refx 'ET' ('PROG' 6.13111E-05);
  298. refx = refx 'ET' ('PROG' 6.16883E-05);
  299. refx = refx 'ET' ('PROG' 6.23597E-05);
  300. refx = refx 'ET' ('PROG' 6.14553E-05);
  301. refx = refx 'ET' ('PROG' 6.23802E-05);
  302. refx = refx 'ET' ('PROG' 6.21707E-05);
  303. refx = refx 'ET' ('PROG' 6.23885E-05);
  304. refx = refx 'ET' ('PROG' 6.21127E-05);
  305. refx = refx 'ET' ('PROG' 6.11701E-05);
  306. refx = refx 'ET' ('PROG' 6.24796E-05);
  307. refx = refx 'ET' ('PROG' 6.25073E-05);
  308. refx = refx 'ET' ('PROG' 6.18036E-05);
  309. refx = refx 'ET' ('PROG' 6.28600E-05);
  310. refx = refx 'ET' ('PROG' 6.19973E-05);
  311. refx = refx 'ET' ('PROG' 6.18662E-05);
  312. refx = refx 'ET' ('PROG' 6.14553E-05);
  313. refx = refx 'ET' ('PROG' 5.95327E-05);
  314. refx = refx 'ET' ('PROG' 5.95444E-05);
  315. refx = refx 'ET' ('PROG' 5.50465E-05);
  316. refx = refx 'ET' ('PROG' 5.01617E-05);
  317. refx = refx 'ET' ('PROG' 4.69227E-05);
  318. refx = refx 'ET' ('PROG' 3.96286E-05);
  319. refx = refx 'ET' ('PROG' 3.28881E-05);
  320. refx = refx 'ET' ('PROG' 2.90179E-05);
  321. refx = refx 'ET' ('PROG' 2.78755E-05);
  322. refx = refx 'ET' ('PROG' 2.84600E-05);
  323. refx = refx 'ET' ('PROG' 3.05303E-05);
  324. refx = refx 'ET' ('PROG' 3.34501E-05);
  325. refx = refx 'ET' ('PROG' 3.68675E-05);
  326. refx = refx 'ET' ('PROG' 4.05489E-05);
  327. refx = refx 'ET' ('PROG' 4.44735E-05);
  328. refx = refx 'ET' ('PROG' 4.82174E-05);
  329. refx = refx 'ET' ('PROG' 5.21294E-05);
  330. refx = refx 'ET' ('PROG' 5.60083E-05);
  331. refx = refx 'ET' ('PROG' 5.99179E-05);
  332. refx = refx 'ET' ('PROG' 6.38235E-05);
  333. refx = refx 'ET' ('PROG' 6.77313E-05);
  334. refx = refx 'ET' ('PROG' 7.16400E-05);
  335. refx = refx 'ET' ('PROG' 7.55490E-05);
  336. refx = refx 'ET' ('PROG' 7.94581E-05);
  337. refx = refx 'ET' ('PROG' 8.33673E-05);
  338. refx = refx 'ET' ('PROG' 8.72765E-05);
  339. refx = refx 'ET' ('PROG' 9.11857E-05);
  340. refx = refx 'ET' ('PROG' 9.50949E-05);
  341. refx = refx 'ET' ('PROG' 9.90041E-05);
  342. refx = refx 'ET' ('PROG' 1.02913E-04);
  343. refx = refx 'ET' ('PROG' 1.06823E-04);
  344. refx = refx 'ET' ('PROG' 1.10732E-04);
  345. refx = refx 'ET' ('PROG' 1.03558E-04);
  346. refx = refx 'ET' ('PROG' 9.48544E-05);
  347. refx = refx 'ET' ('PROG' 8.60156E-05);
  348. refx = refx 'ET' ('PROG' 7.97706E-05);
  349. refx = refx 'ET' ('PROG' 7.60670E-05);
  350. refx = refx 'ET' ('PROG' 7.32897E-05);
  351. refx = refx 'ET' ('PROG' 7.09948E-05);
  352. refx = refx 'ET' ('PROG' 6.90251E-05);
  353. refx = refx 'ET' ('PROG' 6.78843E-05);
  354. refx = refx 'ET' ('PROG' 6.67594E-05);
  355. refx = refx 'ET' ('PROG' 6.62557E-05);
  356. refx = refx 'ET' ('PROG' 6.55069E-05);
  357. refx = refx 'ET' ('PROG' 6.37642E-05);
  358. refx = refx 'ET' ('PROG' 6.21030E-05);
  359. refx = refx 'ET' ('PROG' 6.08528E-05);
  360. refx = refx 'ET' ('PROG' 5.83282E-05);
  361. refx = refx 'ET' ('PROG' 5.44940E-05);
  362. refx = refx 'ET' ('PROG' 4.48749E-05);
  363. refx = refx 'ET' ('PROG' 3.59752E-05);
  364. refx = refx 'ET' ('PROG' 3.02039E-05);
  365. refx = refx 'ET' ('PROG' 2.76942E-05);
  366. refx = refx 'ET' ('PROG' 2.77927E-05);
  367. refx = refx 'ET' ('PROG' 2.96723E-05);
  368. refx = refx 'ET' ('PROG' 3.25695E-05);
  369. refx = refx 'ET' ('PROG' 3.59909E-05);
  370. refx = refx 'ET' ('PROG' 3.96795E-05);
  371. refx = refx 'ET' ('PROG' 4.34990E-05);
  372. refx = refx 'ET' ('PROG' 4.73811E-05);
  373. refx = refx 'ET' ('PROG' 5.12893E-05);
  374.  
  375. refy = 'PROG' 0.;
  376. refy = refy 'ET' ('PROG' 4.45463E+04);
  377. refy = refy 'ET' ('PROG' 5.83801E+04);
  378. refy = refy 'ET' ('PROG' 5.44108E+04);
  379. refy = refy 'ET' ('PROG' 5.00711E+04);
  380. refy = refy 'ET' ('PROG' 4.54985E+04);
  381. refy = refy 'ET' ('PROG' 4.17602E+04);
  382. refy = refy 'ET' ('PROG' 3.92427E+04);
  383. refy = refy 'ET' ('PROG' 3.77858E+04);
  384. refy = refy 'ET' ('PROG' 3.63997E+04);
  385. refy = refy 'ET' ('PROG' 3.56287E+04);
  386. refy = refy 'ET' ('PROG' 3.47398E+04);
  387. refy = refy 'ET' ('PROG' 3.40165E+04);
  388. refy = refy 'ET' ('PROG' 3.33956E+04);
  389. refy = refy 'ET' ('PROG' 3.27009E+04);
  390. refy = refy 'ET' ('PROG' 3.21314E+04);
  391. refy = refy 'ET' ('PROG' 3.15453E+04);
  392. refy = refy 'ET' ('PROG' 3.10697E+04);
  393. refy = refy 'ET' ('PROG' 3.04185E+04);
  394. refy = refy 'ET' ('PROG' 2.96710E+04);
  395. refy = refy 'ET' ('PROG' 2.95760E+04);
  396. refy = refy 'ET' ('PROG' 2.92618E+04);
  397. refy = refy 'ET' ('PROG' 2.87840E+04);
  398. refy = refy 'ET' ('PROG' 2.86147E+04);
  399. refy = refy 'ET' ('PROG' 2.79788E+04);
  400. refy = refy 'ET' ('PROG' 2.77375E+04);
  401. refy = refy 'ET' ('PROG' 2.75490E+04);
  402. refy = refy 'ET' ('PROG' 2.66834E+04);
  403. refy = refy 'ET' ('PROG' 2.66546E+04);
  404. refy = refy 'ET' ('PROG' 2.60293E+04);
  405. refy = refy 'ET' ('PROG' 2.56770E+04);
  406. refy = refy 'ET' ('PROG' 2.50479E+04);
  407. refy = refy 'ET' ('PROG' 2.41931E+04);
  408. refy = refy 'ET' ('PROG' 2.43071E+04);
  409. refy = refy 'ET' ('PROG' 2.38041E+04);
  410. refy = refy 'ET' ('PROG' 2.31315E+04);
  411. refy = refy 'ET' ('PROG' 2.31210E+04);
  412. refy = refy 'ET' ('PROG' 2.22992E+04);
  413. refy = refy 'ET' ('PROG' 2.18663E+04);
  414. refy = refy 'ET' ('PROG' 2.12580E+04);
  415. refy = refy 'ET' ('PROG' 2.01194E+04);
  416. refy = refy 'ET' ('PROG' 1.97150E+04);
  417. refy = refy 'ET' ('PROG' 1.76897E+04);
  418. refy = refy 'ET' ('PROG' 1.56893E+04);
  419. refy = refy 'ET' ('PROG' 1.42315E+04);
  420. refy = refy 'ET' ('PROG' 1.17564E+04);
  421. refy = refy 'ET' ('PROG' 9.42547E+03);
  422. refy = refy 'ET' ('PROG' 7.94773E+03);
  423. refy = refy 'ET' ('PROG' 7.28952E+03);
  424. refy = refy 'ET' ('PROG' 7.14698E+03);
  425. refy = refy 'ET' ('PROG' 7.46348E+03);
  426. refy = refy 'ET' ('PROG' 8.05049E+03);
  427. refy = refy 'ET' ('PROG' 8.80375E+03);
  428. refy = refy 'ET' ('PROG' 9.64809E+03);
  429. refy = refy 'ET' ('PROG' 1.05819E+04);
  430. refy = refy 'ET' ('PROG' 1.14458E+04);
  431. refy = refy 'ET' ('PROG' 1.23745E+04);
  432. refy = refy 'ET' ('PROG' 1.32890E+04);
  433. refy = refy 'ET' ('PROG' 1.42166E+04);
  434. refy = refy 'ET' ('PROG' 1.51423E+04);
  435. refy = refy 'ET' ('PROG' 1.60691E+04);
  436. refy = refy 'ET' ('PROG' 1.69962E+04);
  437. refy = refy 'ET' ('PROG' 1.79235E+04);
  438. refy = refy 'ET' ('PROG' 1.88509E+04);
  439. refy = refy 'ET' ('PROG' 1.97784E+04);
  440. refy = refy 'ET' ('PROG' 2.07058E+04);
  441. refy = refy 'ET' ('PROG' 2.16332E+04);
  442. refy = refy 'ET' ('PROG' 2.25607E+04);
  443. refy = refy 'ET' ('PROG' 2.34881E+04);
  444. refy = refy 'ET' ('PROG' 2.44128E+04);
  445. refy = refy 'ET' ('PROG' 2.53300E+04);
  446. refy = refy 'ET' ('PROG' 2.62386E+04);
  447. refy = refy 'ET' ('PROG' 2.43671E+04);
  448. refy = refy 'ET' ('PROG' 2.21051E+04);
  449. refy = refy 'ET' ('PROG' 1.97443E+04);
  450. refy = refy 'ET' ('PROG' 1.79317E+04);
  451. refy = refy 'ET' ('PROG' 1.67139E+04);
  452. refy = refy 'ET' ('PROG' 1.56873E+04);
  453. refy = refy 'ET' ('PROG' 1.47359E+04);
  454. refy = refy 'ET' ('PROG' 1.38902E+04);
  455. refy = refy 'ET' ('PROG' 1.31620E+04);
  456. refy = refy 'ET' ('PROG' 1.25045E+04);
  457. refy = refy 'ET' ('PROG' 1.19343E+04);
  458. refy = refy 'ET' ('PROG' 1.12926E+04);
  459. refy = refy 'ET' ('PROG' 1.04793E+04);
  460. refy = refy 'ET' ('PROG' 9.64775E+03);
  461. refy = refy 'ET' ('PROG' 8.93150E+03);
  462. refy = refy 'ET' ('PROG' 7.93536E+03);
  463. refy = refy 'ET' ('PROG' 6.81617E+03);
  464. refy = refy 'ET' ('PROG' 4.90968E+03);
  465. refy = refy 'ET' ('PROG' 3.30525E+03);
  466. refy = refy 'ET' ('PROG' 2.08046E+03);
  467. refy = refy 'ET' ('PROG' 1.21911E+03);
  468. refy = refy 'ET' ('PROG' 6.69900E+02);
  469. refy = refy 'ET' ('PROG' 3.53096E+02);
  470. refy = refy 'ET' ('PROG' 1.81002E+02);
  471. refy = refy 'ET' ('PROG' 8.84480E+01);
  472. refy = refy 'ET' ('PROG' 4.01091E+01);
  473. refy = refy 'ET' ('PROG' 1.65356E+01);
  474. refy = refy 'ET' ('PROG' 6.63704E+00);
  475. refy = refy 'ET' ('PROG' 2.64475E+00);
  476.  
  477. tol1 = 1.0e-3;
  478. maxres = 0.;
  479. 'REPE' BTEST (nn - 1 ) ;
  480. V1 = 'EXTR' refx (&BTEST + 1);
  481. V2 = 'EXTR' prod (&BTEST + 1);
  482. res1 = 'ABS' ((V1 - V2)/V1);
  483. 'SI' (res1 > maxres); maxres = res1; 'FINSI';
  484.  
  485. W1 = 'EXTR' refy (&BTEST + 1);
  486. W2 = 'EXTR' prof (&BTEST + 1);
  487. res2 = 'ABS' ((W1 - W2)/W1);
  488. 'SI' (res2 > maxres); maxres = res2; 'FINSI';
  489.  
  490. 'FIN' BTEST ;
  491.  
  492. err1 = maxres > tol1;
  493.  
  494. 'SI' err1;
  495. 'ERRE' 5 ;
  496. 'FINSI';
  497.  
  498. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  499. *
  500. * END OF ERROR TESTING
  501. *
  502. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  503. 'FIN';
  504.  
  505.  
  506.  
  507.  

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