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

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