Télécharger raff04.dgibi

Retour à la liste

Numérotation des lignes :

  1. ************************************************************
  2. ** raff04.dgibi
  3. ************************************************************
  4. * *
  5. * Calcul elasitique 3D des facteurs d'intencité de *
  6. * contraintes pour une fissure circulaire sous un *
  7. * chargempent de traction à l'infini incliné de 45 degrte *
  8. * par rapport au pla de fissure (mode mixte) *
  9. * *
  10. * Test De RAFF avec X-FEM et GTHETA *
  11. * *
  12. * Comparaison de la simulations avec une solution de *
  13. * reference pour k1, k2 et k3 donne par : *
  14. * [TADA, STRESS ANALYSIS HANDBOOK,1973] *
  15. * *
  16. *création : gg, le 15.03.2017 *
  17. ************************************************************
  18. *** Options
  19.  
  20. *** options de calcul
  21. opti dime 3 elem CUB8 mode trid;
  22.  
  23. *** options de trace
  24. graph = faux ;
  25. *graph = vrai ;
  26.  
  27. *complet = vrai;
  28. complet = faux;
  29.  
  30. opti 'TRAC' X ;
  31. *opti 'TRAC' psc ;
  32.  
  33.  
  34. *******************************************************
  35. *** Maillage : cube sain avec raff
  36.  
  37.  
  38. * Dimensions du maillage
  39. L1 = 0.2 ;
  40. L2 = L1 / 2.;
  41. * L2 = L1 * (3./4.);
  42. L4 = L1 / 4.;
  43. * L4 = L1 / 5.;
  44. a0 = L1 / 10 ;
  45.  
  46. na1 = 4;
  47.  
  48.  
  49. * paramettre de raffinement
  50.  
  51. dmin= 1.5 * L4;
  52. densmax= L4 ;
  53. *densmin=densmax/4;
  54. densmin=densmax/8.;
  55. dmax=dmin + densmax;
  56.  
  57. si (complet) ;
  58. dmin2= 0.7 * dmin;
  59. densmax2= dmin ;
  60. densmin2=densmax/16;
  61. *densmin=densmax/16.;
  62. dmax2=dmin2 + densmax2;
  63. finsi;
  64. * maillage de la fissure
  65. dx1 = a0 / 8;
  66.  
  67. dx1elim = 1.E-4 * dx1;
  68. DENS dx1;
  69.  
  70.  
  71.  
  72.  
  73. * oeuil pour la visualisation
  74. eye1 = -100. -70. 30. ;
  75. eye2 = 0. -0. 100. ;
  76. eye3 = -100. -30. 10.;
  77.  
  78. * boite dallée (a)
  79. vzL1 = (0. 0. L1);
  80. pa00 = 0. 0. 0.;
  81. pa01 = L1 0. 0.;
  82. pa02 = L1 L1 0.;
  83. pa03 = 0. L1 0.;
  84. la01 = pa00 droi na1 pa01;
  85. la02 = pa01 droi na1 pa02;
  86. la03 = pa02 droi na1 pa03;
  87. la04 = pa03 droi na1 pa00;
  88.  
  89. * on cree des surface qu on eliminera apres coup
  90. la06 = la02 plus vzL1;
  91. la07 = la03 plus vzL1;
  92. sa02 = la02 regl na1 la06;
  93. sa03 = la03 regl na1 la07;
  94.  
  95. * creation du volume
  96. sa01 = dall la01 la02 la03 la04 'PLAN';
  97. va01 = sa01 volu na1 'TRAN' vzL1;
  98. * recup
  99. pa04 = va01 POIN 'PROCH' (pa00 plus vzL1);
  100. pa05 = va01 POIN 'PROCH' (pa01 plus vzL1);
  101. pa06 = va01 POIN 'PROCH' (pa02 plus vzL1);
  102. pa07 = va01 POIN 'PROCH' (pa03 plus vzL1);
  103. elim va01 sa02 dx1elim;
  104. elim va01 sa03 dx1elim;
  105. sa04 = va01 face 2;
  106.  
  107.  
  108.  
  109. ***********************************************************
  110. *** Maillage de la fissure
  111. &ELEM = VALE 'ELEM';
  112. opti elem TRI3;
  113.  
  114. *** petite fissure circulaire inclinée
  115. * beta0 = 0.;
  116. * beta0 = 20.;
  117. beta0 = 45.;
  118. * a0 = L1 / 10.; deja defini + haut pour caler la finesse du maillage
  119. rapab0 = 1.;
  120. b0 = a0 * rapab0;
  121. * finesse du maillage de la fissure
  122. da0 = dx1 / 4.;
  123. si (complet);
  124. da0 = da0/2.;
  125. finsi;
  126. * da0 = dx1 / 12.;
  127. na0 = enti (a0 / da0);
  128. nb0 = enti (b0 / da0);
  129.  
  130. * on maille un cercle dans le plan
  131. * cad avec b0=a0 : on ajustera apres si ellipse
  132. * on tournera apres si inclinee
  133. DENS da0;
  134. * on fait un carré
  135. pf0 = 0. 0. 0.;
  136. xf1 = a0 * ( (2**(-1.5)) + 0.5) / 2.;
  137. pf1 = xf1 0. 0.;
  138. pf2 = xf1 xf1 0.;
  139. pf3 = 0. xf1 0.;
  140. lf1 = pf0 droi pf1;
  141. lf2 = pf1 droi pf2;
  142. lf6 = pf2 droi pf1;
  143. lf3 = pf2 droi pf3;
  144. lf5 = pf3 droi pf2;
  145. lf4 = pf3 droi pf0;
  146. sf1 = dall lf1 lf2 lf3 lf4 'PLAN';
  147. sf1 = inve sf1;
  148. nf1 = nbel lf1;
  149. * on fait le disque
  150. pf4 = a0. 0. 0.;
  151. xf5 = a0 / (2**0.5);
  152. pf5 = xf5 xf5 0.;
  153. pf6 = 0. a0. 0.;
  154. cf5 = CERC nf1 pf6 pf0 pf5;
  155. cf6 = CERC nf1 pf5 pf0 pf4;
  156. sf2 = regl (cf5 et cf6) (lf5 et lf6);
  157. * assemblage
  158. sfron0 = (cf5 et cf6) coul 'ROUG';
  159. scrack0 = (sf1 et sf2) coul 'ORAN';
  160. * rotation de beta0
  161. DEPL scrack0 'TOUR' beta0 pf0 pf1;
  162. * tres important
  163.  
  164. scrack0 = ELIM scrack0 dx1elim;
  165. sfron0 = ELIM sfron0 dx1elim;
  166.  
  167.  
  168. ***********************************************************
  169. *symétries
  170. * de la structure
  171. * par % a z=0
  172.  
  173. va01bis = va01 syme 'PLAN' pa00 pa01 pa03;
  174.  
  175. va01 = va01 et va01bis;
  176.  
  177. * par % a y=0
  178. va01bis = va01 syme 'PLAN' pa00 pa01 pa04;
  179.  
  180. va01 = va01 et va01bis;
  181. ELIM va01 dx1elim;
  182.  
  183. * par % a X=0
  184.  
  185. va01bis = va01 syme 'PLAN' pa00 pa03 pa04;
  186.  
  187. va01 = va01 et va01bis;
  188. ELIM va01 dx1elim;
  189.  
  190.  
  191. * de la fissure
  192. sfron0b scrack0b = sfron0 scrack0 syme 'DROIT' pa00 pa01;
  193. scrack0b = INVE scrack0b ;
  194. sfron0b = INVE sfron0b ;
  195.  
  196. scrack0 = scrack0 et scrack0b;
  197. sfron0 = sfron0 et sfron0b;
  198.  
  199. sfron0b scrack0b = sfron0 scrack0 syme 'DROIT' pa00 pf6;
  200. scrack0b = INVE scrack0b ;
  201. sfron0b = INVE sfron0b ;
  202.  
  203. scrack0 = scrack0 et scrack0b;
  204. sfron0 = sfron0 et sfron0b;
  205.  
  206. ELIM scrack0 dx1elim;
  207. ELIM sfron0 dx1elim;
  208.  
  209. **********************************************************
  210. *** Raffinement autour de la fissure
  211.  
  212. xa ya za= coor va01 ;
  213.  
  214. DF1= (((ya**2) + (xa**2.) + (za**2))**0.5);
  215.  
  216. aff = (((densmax - densmin)*(DF1 - dmax))/(dmax - dmin)) + densmax ;
  217. dens1 = born aff 'SCAL' 'COMPRIS' (densmin*1.000001)
  218. (densmax*1.000001);
  219.  
  220.  
  221. vb01 = raff dens1 va01;
  222.  
  223. si (complet);
  224. vb01 = raff dens1 va01;
  225.  
  226. xb yb zb= coor vb01 ;
  227.  
  228. DF2= (((yb**2) + (xb**2.) + (zb**2))**0.5);
  229. aff2 = (((densmax2 - densmin2)*(DF2 - dmax2))/(dmax2 - dmin2))
  230. + densmax2 ;
  231. dens2 = born aff2 'SCAL' 'COMPRIS' (densmin2*1.000001)
  232. (densmax*21.000001);
  233.  
  234.  
  235. vc01 = raff dens2 vb01;
  236.  
  237. vol0 = ELIM vc01 dx1elim;
  238. sinon ;
  239. vol0 = ELIM vb01 dx1elim;
  240. finsi;
  241.  
  242.  
  243. si (graph) ;
  244. trac va01 eye1 'COUP' pf0 pf1 pf3 'CACH'
  245. 'TITR' 'Maillage avant raff';
  246. trac dens1 va01 eye1 'COUP' pf0 pf1 pf3 'TITR' 'Champ de densite';
  247. trac vb01 eye1 'COUP' pf0 pf1 pf3 'CACH' 'TITR'
  248. 'Maillage apres raff';
  249. si (complet);
  250. trac dens2 vb01 eye1 'COUP' pf0 pf1 pf3
  251. 'TITR' ' Deuxieme champ de densite';
  252. trac vc01 eye1 'COUP' pf0 pf1 pf3 'CACH' 'TITR'
  253. 'Maillage apres raff';
  254. finsi;
  255. finsi ;
  256.  
  257. **********************************************************************
  258. * séparation en deux maillages (fem et X-fem)
  259. x0 y0 z0 = coor vol0;
  260.  
  261. si (complet);
  262. zz0 = COLI y0 (cos beta0) z0 (-1.* (sin beta0));
  263. yy0 = COLI y0 (sin beta0) z0 (cos beta0);
  264. pxfem_r = ((( x0**2) + (yy0**2))**0.5) POIN 'COMPRIS'
  265. (0.2*a0) (0.6 * dmin) ;
  266. pxfem_h = (abs zz0) POIN 'EGINFE' (0.12 * dmin) ;
  267. pxfem = INTE pxfem_r pxfem_h;
  268. sinon;
  269. * pxfem =(((y0**2) + (x0**2.) + (z0**2))**0.5) POIN
  270. * 'EGINFE' (0.6 * dmin) ;
  271. zz0 = COLI y0 (cos beta0) z0 (-1.* (sin beta0));
  272. yy0 = COLI y0 (sin beta0) z0 (cos beta0);
  273. pxfem_r = (((x0**2) + (yy0**2))**0.5) POIN 'COMPRIS'
  274. (0.6*a0) (0.3 * dmin) ;
  275. pxfem_h = (abs zz0) POIN 'EGINFE' (0.12 * dmin) ;
  276. pxfem = INTE pxfem_r pxfem_h;
  277. finsi;
  278. *
  279.  
  280. vol1 = vol0 ELEM 'APPU' 'LARG' pxfem;
  281. vol1 = vol1 coul 'ROUG';
  282. vol2 = vol0 DIFF vol1;
  283.  
  284. *** construction des level set ***
  285. * psi1 phi1 = PSIP vol1 scrack0 'DEUX' sfron0;
  286. psi1 phi1 = PSIP vol0 scrack0 'DEUX' sfron0;
  287.  
  288. * qq objets sympas pour le post traitement
  289. con0 = aret vol0;
  290. env0 = ENVE vol0;
  291.  
  292. * qq objets utiles pour les CL
  293. x2 y2 z2 = coor vol2;
  294. lignx0 = ((x2**2)+(z2**2)) poin 'EGINFE' dx1elim;
  295. poiy0 = pa00 ;
  296. facz0 = z0 poin 'EGINFE' ((-1. * L1) + dx1elim);
  297. facz2 = z0 poin 'EGSUPE' (L1 - dx1elim);
  298. facz2 = env0 elem 'APPU' 'STRICT' facz2;
  299.  
  300. * tracés
  301. mess ' RAFF va01 = vol0 ';
  302. mess ' NBNO=' (nbno va01) (nbno vol0);
  303. mess ' NBEL=' (nbel va01) (nbel vol0);
  304.  
  305.  
  306. * opti donn 5 ;
  307. **********************************************************
  308. *** MODELE & MATERIAU ***
  309. **********************************************************
  310.  
  311. * coef elastique
  312. Ey0 = 2.E11;
  313. * nu0 = 0.03;
  314. nu0 = 0.3 ;
  315. rho0= 7800. ;
  316. mu0 = (0.5 * Ey0) / (1. + nu0) ;
  317. Estar0 = Ey0 / (1. - (nu0**2));
  318.  
  319. * Loi de Paris (on convertit C : mm -> m)
  320. * mParis = 3.;
  321. * CParis = 1.E-9 * 1.E-3;
  322.  
  323.  
  324. **********************************************************
  325. *** Modele & materiau ***
  326. *
  327. * zone de propagation (X-FEM)
  328. mod1 = MODE vol1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 'XC8R';
  329. mat1 = MATE mod1 'YOUN' Ey0 'NU' nu0 'RHO' rho0 ;
  330. che1X = TRIE mod1 psi1 phi1 ;
  331.  
  332. * elements standards
  333. mod2 = MODE vol2 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 'CUB8';
  334. mat2 = MATE mod2 'YOUN' Ey0 'NU' nu0 'RHO' rho0 ;
  335. *
  336. mod1tot = mod1 et mod2 ;
  337. mat1tot = mat1 et mat2 ;
  338.  
  339.  
  340.  
  341. **********************************************************
  342. *** CL et DEPLACEMENTS IMPOSES ***
  343. **********************************************************
  344. szz0 = -100.E6;
  345.  
  346.  
  347. *** CL ***
  348.  
  349. *blocage des mouvements de coprs rigide
  350.  
  351. cl0x = bloq lignx0 'UX';
  352. cl0y = bloq poiy0 'UY';
  353. cl0z = bloq facz0 'UZ';
  354.  
  355. cl1tot = (cl0x et cl0y et cl0z ) ;
  356.  
  357. *** chargement en mode mixte ***
  358.  
  359. for1tot = PRES 'MASS' mod1tot facz2 szz0;
  360.  
  361. si(graph);
  362. trac (vect for1tot 'FORC' 'BLEU') ((coul lignx0 'VERT') et
  363. (coul facz0 'TURQ') et con0 et scrack0 et sfron0)
  364. 'TITR' 'CHARGEMENT ET CONDITIONS AUX LIMITES';
  365. trac (con0 et scrack0 et (enve vol1)) TITR 'Erichissement';
  366. fins;
  367.  
  368.  
  369. *** matrice elastique et assemblage
  370. Rig1tot = RIGI mod1tot mat1tot;
  371. K1tot = Rig1tot et cl1tot;
  372.  
  373.  
  374. *******************************************************
  375. * CALCUL ELASTIQUE
  376. *******************************************************
  377.  
  378. u1 = RESO K1tot for1tot ;
  379.  
  380.  
  381. * ****************************
  382. * *** deformees ***
  383. *
  384. si(graph);
  385. u1phy = XFEM 'RECO' u1 mod1tot ;
  386. sig1 = SIGM mod1tot u1 mat1tot;
  387. ucrk1u ucrk1d = XFEM 'FISS' scrack0 u1 mod1tot ;
  388. amp1 = 200.;
  389. def2 = (DEFO u1phy con0 amp1)
  390. et (DEFO ucrk1u scrack0 amp1 'ORAN')
  391. et (DEFO ucrk1d scrack0 amp1 'ROUG');
  392. TRAC def2 eye3 'CACH' 'TITR' 'Deformee';
  393. TRAC u1 vol0 eye1 'COUP' pf0 pf1 pf3 'TITR' 'Deplacements';
  394. TRAC sig1 mod1tot eye1 'COUP' pf0 pf1 pf3 'TITR' 'Contraintes';
  395. finsi;
  396. * opti donn 5;
  397.  
  398.  
  399. ************************************************************************
  400. * On calcule les facteurs d'intensité de contraintes *
  401. * sur une moitier du front car G_theta ne prévoi pas *
  402. * de faire ce calcul en 3D pour des front déssinant de *
  403. * contours fermé. *
  404. ************************************************************************
  405.  
  406. xc yc zc = coor scrack0;
  407. xneg = xc poin 'EGINF' 0.0 ;
  408. xpos = xc poin 'EGSUP' 0.0 ;
  409. sfron1 = sfron0 elem 'APPU' 'STRI' xpos COUL VERT;
  410. sfron2 = sfron0 elem 'APPU' 'STRI' xneg;
  411.  
  412. si (graph);
  413. TRAC (SFRON1 et SFRON2) 'TITR' 'Deux demi fronts';
  414. finsi;
  415. *** Calcul des FIC par methode integrale G_THETA ***
  416. KTAB = tabl ;
  417. KTAB . 'OBJECTIF' = MOT 'DECOUPLAGE';
  418. *KTAB . 'OBJECTIF' = MOT 'J';
  419. KTAB .'PSI' = psi1;
  420. KTAB .'PHI' = phi1;
  421. KTAB .'FRONT_FISSURE' = sfron1 ;
  422. KTAB .'MODELE' = mod1tot;
  423. KTAB .'CARACTERISTIQUES' = mat1tot;
  424. KTAB .'SOLUTION_RESO' = u1;
  425. KTAB .'CHARGEMENTS_MECANIQUES' = for1tot;
  426. KTAB . 'POINT_CENTRE' = PA00;
  427. KTAB . 'COUCHE' = 2;
  428.  
  429. G_THETA KTAB;
  430.  
  431. Kchpo1 = KTAB . 'CHPO_RESULTATS' ;
  432.  
  433.  
  434. K1ev = EVOL BLEU 'CHPO' Kchpo1 'K1' sfron1;
  435. K2ev = EVOL DEFA 'CHPO' Kchpo1 'K2' sfron1;
  436. K3ev = EVOL ROUG 'CHPO' Kchpo1 'K3' sfron1;
  437.  
  438. S1pr = EXTR K1ev 'ABSC';
  439. K1pr = EXTR K1ev 'ORDO';
  440. K2pr = EXTR K2ev 'ORDO';
  441. K3pr = EXTR K3ev 'ORDO';
  442.  
  443. ************************************************************************
  444. * Formule analytique des FIC en milieu infinie *
  445. * pour une fissure circulaire inclinée *
  446. * "The stress analysis of crack Handook" H.TADA P.C0.Paris G.R.IRWIN *
  447. ************************************************************************
  448.  
  449. K1A = (2. * (-1.0)* szz0 * ((sin (-1.0 * beta0))**2) *
  450. (( a0 / pi )**0.5))+ ( S1pr * 0.);
  451.  
  452. CA = (4. * (-1.0)* szz0 * (sin (-1.0 * beta0)) * (cos (-1.0 * beta0)) *
  453. (( a0 / pi )**0.5)) / (2. - nu0);
  454.  
  455. K2A= CA * ( cos ((S1pr * 180. ) / (a0 * pi )));
  456.  
  457. K3A= CA * (1. - nu0)* ( sin ((S1pr * 180. ) / (a0 * pi )));
  458.  
  459. K1Aev = EVOL BLEU 'MANU' S1pr K1A;
  460. K2Aev = EVOL DEFA 'MANU' S1pr K2A;
  461. K3Aev = EVOL ROUG 'MANU' S1pr K3A;
  462. si (graph);
  463. ** Calcul des FIC par methode integrale G_THETAsur l autre demi front **
  464. KSTAB = tabl ;
  465. KSTAB . 'OBJECTIF' = MOT 'DECOUPLAGE';
  466. * KSTAB . 'OBJECTIF' = MOT 'J';
  467. KSTAB .'PSI' = psi1;
  468. KSTAB .'PHI' = phi1;
  469. KSTAB .'FRONT_FISSURE' = sfron2 ;
  470. KSTAB .'MODELE' = mod1tot;
  471. KSTAB .'CARACTERISTIQUES' = mat1tot;
  472. KSTAB .'SOLUTION_RESO' = u1;
  473. KSTAB .'CHARGEMENTS_MECANIQUES' = for1tot;
  474. KSTAB . 'POINT_CENTRE' = PA00;
  475. KSTAB . 'COUCHE' = 2;
  476.  
  477. G_THETA KSTAB;
  478.  
  479. KSchpo1 = KSTAB . 'CHPO_RESULTATS' ;
  480.  
  481. KS1ev = EVOL VERT 'CHPO' KSchpo1 'K1' sfron2;
  482. KS2ev = EVOL ORAN 'CHPO' KSchpo1 'K2' sfron2;
  483. KS3ev = EVOL ROSE 'CHPO' KSchpo1 'K3' sfron2;
  484.  
  485. SS1pr = EXTR KS1ev 'ABSC';
  486. KS1pr = EXTR KS1ev 'ORDO';
  487. KS2pr = EXTR KS2ev 'ORDO';
  488. KS3pr = EXTR KS3ev 'ORDO';
  489.  
  490. KS1A = (2. * (-1.0)* szz0 * ((sin (-1.0 * beta0))**2) *
  491. (( a0 / pi )**0.5))+ ( S1pr * 0.);
  492.  
  493. CSA = (4. * (-1.0)* szz0 * (sin (-1.0 * beta0)) *
  494. (cos (-1.0 * beta0)) * (( a0 / pi )**0.5)) / (2. - nu0);
  495.  
  496. KS2A= CSA * ( cos ((S1pr * 180. ) / (a0 * pi ) + 180.));
  497.  
  498. KS3A= CSA * (1. - nu0)* ( sin ((SS1pr * 180. ) / (a0 * pi ) + 180. ));
  499.  
  500. KS1Aev = EVOL VERT 'MANU' SS1pr KS1A;
  501. KS2Aev = EVOL ORAN 'MANU' SS1pr KS2A;
  502. KS3Aev = EVOL ROSE 'MANU' SS1pr KS3A;
  503.  
  504.  
  505. tdess1 = tabl;
  506. tdess1.'TITRE'= TABL;
  507. tdess1.'TITRE' . 1= MOT 'KI sim';
  508. tdess1.'TITRE' . 2= MOT 'KII sim ';
  509. tdess1.'TITRE' . 3= MOT 'KIII sim ';
  510. tdess1.'TITRE' . 4= MOT 'KI ana';
  511. tdess1.'TITRE' . 5= MOT 'KII ana ';
  512. tdess1.'TITRE' . 6= MOT 'KIII ana ';
  513. tdess1 . 4 = 'MOT' ' TIRR ';
  514. tdess1 . 5 = 'MOT' ' TIRR ';
  515. tdess1 . 6 = 'MOT' 'TIRR ';
  516. dess (K1ev et K2ev et K3ev et K1Aev et K2Aev et K3Aev )
  517. 'LEGE' tdess1 'TITR' 'SIF sur un demi front';
  518. dess (KS1ev et KS2ev et KS3ev et KS1Aev et KS2Aev et KS3Aev)
  519. 'LEGE' tdess1 'TITR' 'SIF sur l autre demi front';
  520. fins;
  521. ************************************************************************
  522. * Erreur en norme L2 des FIC le long du front *
  523. * (On s'eloigne des bords du demi front sur lequel on a fait *
  524. * le calcul g-theta pour ignorer les effet de bords dues a G-theta) *
  525. ************************************************************************
  526. si (complet);
  527. dint = 17;
  528. fint = 92;
  529. ip6 = 55;
  530. ip5 = 28;
  531. sinon ;
  532. dint = 6;
  533. fint = 51;
  534. ip6 = 29;
  535. ip5 = 15;
  536. finsi;
  537.  
  538. DiffL2K1 =((K1Aev - K1ev)**2) / ((EXTR K1A ip6)**2);
  539. DiffL2K2 =((K2Aev - K2ev)**2) / ((EXTR K2A 1)**2);
  540. DiffL2K3 =((K3Aev - K3ev)**2) / ((EXTR K3A ip6)**2);
  541.  
  542. ErrL2K1 =((intg DiffL2K1 dint fint)/
  543. ((EXTR S1pr fint ) -(EXTR S1pr dint)))**0.5;
  544. ErrL2K2 =((intg DiffL2K2 dint fint)/
  545. ((EXTR S1pr fint ) - (EXTR S1pr dint)))**0.5;
  546. ErrL2K3 =((intg DiffL2K3 dint fint)/
  547. ((EXTR S1pr fint ) -(EXTR S1pr dint)))**0.5;
  548. ErrL2K1 = extr ErrL2K1 1;
  549. ErrL2K2 = extr ErrL2K2 1;
  550. ErrL2K3 = extr ErrL2K3 1;
  551. list ErrL2K1;
  552. list ErrL2K2;
  553. list ErrL2K3;
  554.  
  555.  
  556. si (graph);
  557. opti donn 5;
  558. finsi;
  559.  
  560. SI ((ErrL2K1 < (0.05)) et (ErrL2K2 < (0.10)) et (ErrL2K3 < 0.15));
  561. 'MESS' ' ----------------------' ;
  562. 'MESS' ' SUCCES DU CAS-TEST !' ;
  563. 'MESS' ' ----------------------' ;
  564. Erre 0;
  565. Sinon;
  566. 'MESS' ' ---------------------' ;
  567. 'MESS' ' ECHEC DU CAS-TEST !' ;
  568. 'MESS' ' ---------------------' ;
  569. Erre 5;
  570. finsi;
  571.  
  572. *ind1K = inde (KTAB . 'CHAMP_THETA');
  573. *teta1 = KTAB . 'CHAMP_THETA' . (ind1K . 1);
  574. *trac (vect teta1 'BLEU') (vol0 et scrack0);
  575.  
  576.  
  577. TEMP 'IMPR' 'MAXI' 'CPU';
  578.  
  579. FIN ;
  580.  
  581.  
  582.  
  583.  
  584.  
  585.  

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