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.  
  25. graph = faux ;
  26. *graph = vrai ;
  27.  
  28. *complet = vrai;
  29. complet = faux;
  30.  
  31. opti 'TRAC' X ;
  32. *opti 'TRAC' psc ;
  33.  
  34.  
  35. *******************************************************
  36. *** Maillage : cube sain avec raff
  37.  
  38.  
  39. * Dimensions du maillage
  40. L1 = 0.2 ;
  41. L2 = L1 / 2.;
  42. L4 = L1 / 4.;
  43. a0 = L1 / 10 ;
  44. *a0 = L1 / 7.5 ;
  45.  
  46.  
  47. na1 = 4;
  48.  
  49.  
  50. * paramettre de raffinement
  51.  
  52. dmin= 1.5 * L4;
  53. densmax= L4 ;
  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. dmax2=dmin2 + densmax2;
  62. finsi;
  63.  
  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. p001 = 0. 0. 1.;
  245. p010 = 0. 1. 0.;
  246. trac va01 eye1 'COUP' pf0 pf1 pf3 'CACH'
  247. 'TITR' 'Maillage avant raff';
  248. trac dens1 va01 eye1 'COUP' pf0 pf1 pf3 'TITR' 'Champ de densite';
  249. trac dens1 va01 eye1 'COUP' pf0 p010 p001 'TITR' 'Champ de densite';
  250. trac vb01 eye1 'COUP' pf0 pf1 pf3 'CACH' 'TITR'
  251. 'Maillage apres raff';
  252. trac vb01 eye1 'COUP' pf0 p010 p001 'CACH' 'TITR'
  253. 'Maillage apres raff';
  254. si (complet);
  255. trac dens2 vb01 eye1 'COUP' pf0 pf1 pf3
  256. 'TITR' ' Deuxieme champ de densite';
  257. trac vc01 eye1 'COUP' pf0 pf1 pf3 'CACH' 'TITR'
  258. 'Maillage apres raff';
  259. finsi;
  260. finsi ;
  261.  
  262. **********************************************************************
  263. * séparation en deux maillages (fem et X-fem)
  264. x0 y0 z0 = coor vol0;
  265.  
  266. si (complet);
  267. zz0 = COLI y0 (cos beta0) z0 (-1.* (sin beta0));
  268. yy0 = COLI y0 (sin beta0) z0 (cos beta0);
  269. pxfem_r = ((( x0**2) + (yy0**2))**0.5) POIN 'COMPRIS'
  270. (0.2*a0) (0.6 * dmin) ;
  271.  
  272. pxfem_h = (abs zz0) POIN 'EGINFE' (0.12 * dmin) ;
  273. pxfem = INTE pxfem_r pxfem_h;
  274. sinon;
  275. * pxfem =(((y0**2) + (x0**2.) + (z0**2))**0.5) POIN
  276. * 'EGINFE' (0.6 * dmin) ;
  277. zz0 = COLI y0 (cos beta0) z0 (-1.* (sin beta0));
  278. yy0 = COLI y0 (sin beta0) z0 (cos beta0);
  279. pxfem_r = (((x0**2) + (yy0**2))**0.5) POIN 'COMPRIS'
  280. (0.6*a0) (0.3 * dmin) ;
  281. pxfem_h = (abs zz0) POIN 'EGINFE' (0.12 * dmin) ;
  282. pxfem = INTE pxfem_r pxfem_h;
  283. finsi;
  284. *
  285.  
  286. vol1 = vol0 ELEM 'APPU' 'LARG' pxfem;
  287. *vol1 = vol1 coul 'ROUG';
  288. vol2 = vol0 DIFF vol1;
  289.  
  290. *** construction des level set ***
  291. * psi1 phi1 = PSIP vol1 scrack0 'DEUX' sfron0;
  292. psi1 phi1 = PSIP vol0 scrack0 'DEUX' sfron0;
  293.  
  294. * qq objets sympas pour le post traitement
  295. con0 = aret vol0;
  296. env0 = ENVE vol0;
  297.  
  298. * qq objets utiles pour les CL
  299. x2 y2 z2 = coor vol2;
  300. lignx0 = ((x2**2)+(z2**2)) poin 'EGINFE' dx1elim;
  301. poiy0 = pa00 ;
  302. facz0 = z0 poin 'EGINFE' ((-1. * L1) + dx1elim);
  303. facz2 = z0 poin 'EGSUPE' (L1 - dx1elim);
  304. facz0 = env0 elem 'APPU' 'STRICT' facz0;
  305. facz2 = env0 elem 'APPU' 'STRICT' facz2;
  306.  
  307. * tracés
  308. mess ' RAFF va01 = vol0 ';
  309. mess ' NBNO=' (nbno va01) (nbno vol0);
  310. mess ' NBEL=' (nbel va01) (nbel vol0);
  311.  
  312.  
  313. * opti donn 5 ;
  314. **********************************************************
  315. *** MODELE & MATERIAU ***
  316. **********************************************************
  317.  
  318. * coef elastique
  319. Ey0 = 2.E11;
  320. * nu0 = 0.03;
  321. nu0 = 0.3 ;
  322. rho0= 7800. ;
  323. mu0 = (0.5 * Ey0) / (1. + nu0) ;
  324. Estar0 = Ey0 / (1. - (nu0**2));
  325.  
  326. * Loi de Paris (on convertit C : mm -> m)
  327. * mParis = 3.;
  328. * CParis = 1.E-9 * 1.E-3;
  329.  
  330.  
  331. **********************************************************
  332. *** Modele & materiau ***
  333. *
  334. * zone de propagation (X-FEM)
  335. mod1 = MODE vol1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 'XC8R';
  336. mat1 = MATE mod1 'YOUN' Ey0 'NU' nu0 'RHO' rho0 ;
  337. che1X = TRIE mod1 psi1 phi1 ;
  338.  
  339. * constructionsion des blocages des ddl X-fem non actifs dans
  340. * les éléments de transition et des relations de conformité dues au
  341. *raffinement dans cette zone.
  342. rel1X = RELA mod1;
  343.  
  344. * elements standards
  345. mod2 = MODE vol2 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 'CUB8';
  346. mat2 = MATE mod2 'YOUN' Ey0 'NU' nu0 'RHO' rho0 ;
  347. *
  348. * constructionsion et des relations de conformité dues au
  349. *raffinement dans cette zone.
  350. rel2 = RELA mod2;
  351.  
  352.  
  353. mod1tot = mod1 et mod2 ;
  354. mat1tot = mat1 et mat2 ;
  355. rel1tot = rel1x et rel2;
  356.  
  357.  
  358. **********************************************************
  359. *** CL et DEPLACEMENTS IMPOSES ***
  360. **********************************************************
  361. szz0 = -100.E6;
  362.  
  363.  
  364. *** CL ***
  365. px = L1 0. 0. ;
  366. px = POIN vol2 'PROC' px;
  367. py = 0. L1 0. ;
  368. py = POIN vol2 'PROC' py;
  369. pz = 0. 0. L1 ;
  370. pz = POIN vol2 'PROC' pz;
  371.  
  372. *blocage des mouvements de coprs rigide
  373.  
  374. cl0x = bloq px 'UY' 'UZ';
  375. cl0y = bloq py 'UX' 'UZ';
  376. cl0z = bloq pz 'UX' 'UY';
  377.  
  378. cl1tot = (cl0x et cl0y et cl0z ) ;
  379.  
  380. *** chargement en mode mixte ***
  381.  
  382. for1 = PRES 'MASS' mod1tot facz2 szz0;
  383. for2 = PRES 'MASS' mod1tot facz0 szz0;
  384. for1tot = for1 et for2 ;
  385.  
  386. si(graph);
  387. trac (vect for1tot 'FORC' 'BLEU') ((coul (px et py et pz )
  388. 'VERT') et con0 et scrack0 et sfron0)
  389. 'TITR' 'CHARGEMENT ET CONDITIONS AUX LIMITES';
  390. trac (con0 et scrack0 et ((enve vol1) coul ROUG))
  391. TITR 'Erichissement';
  392. fins;
  393.  
  394.  
  395. *** matrice elastique et assemblage
  396. Rig1tot = RIGI mod1tot mat1tot;
  397. K1tot = Rig1tot et cl1tot et rel1tot;
  398.  
  399.  
  400. *******************************************************
  401. * CALCUL ELASTIQUE
  402. *******************************************************
  403.  
  404. u1 = RESO K1tot for1tot ;
  405.  
  406.  
  407. * ****************************
  408. * *** deformees ***
  409. *
  410. si(graph);
  411. u1phy = XFEM 'RECO' u1 mod1tot ;
  412. sig1 = SIGM mod1tot u1 mat1tot;
  413. ucrk1u ucrk1d = XFEM 'FISS' scrack0 u1 mod1tot ;
  414. amp1 = 200.;
  415. def2 = (DEFO u1phy con0 amp1)
  416. et (DEFO ucrk1u scrack0 amp1 'ORAN')
  417. et (DEFO ucrk1d scrack0 amp1 'ROUG');
  418. TRAC def2 eye3 'CACH' 'TITR' 'Deformee';
  419. TRAC u1 vol0 eye1 'COUP' pf0 pf1 pf3 'TITR' 'Deplacements';
  420. TRAC u1 vol0 eye1 'COUP' pf0 p010 p001 'TITR' 'Deplacements';
  421. TRAC sig1 mod1tot eye1 'COUP' pf0 pf1 pf3 'TITR' 'Contraintes';
  422. TRAC sig1 mod1tot eye1 'COUP' pf0 p010 p001 'TITR' 'Contraintes';
  423. finsi;
  424. * opti donn 5;
  425.  
  426.  
  427. ************************************************************************
  428. * On calcule les facteurs d'intensité de contraintes *
  429. * sur une moitier du front car G_theta ne prévoi pas *
  430. * de faire ce calcul en 3D pour des front déssinant de *
  431. * contours fermé. *
  432. ************************************************************************
  433.  
  434.  
  435.  
  436. xc yc zc = coor scrack0;
  437. xneg = xc poin 'EGINF' 0.0 ;
  438. xpos = xc poin 'EGSUP' 0.0 ;
  439. sfron1 = sfron0 elem 'APPU' 'STRI' xpos COUL VERT;
  440. sfron2 = sfron0 elem 'APPU' 'STRI' xneg;
  441.  
  442. si (graph);
  443. TRAC (SFRON1 et SFRON2) 'TITR' 'Deux demi fronts';
  444. finsi;
  445. *** Calcul des FIC par methode integrale G_THETA ***
  446. KTAB = tabl ;
  447. KTAB . 'OBJECTIF' = MOT 'DECOUPLAGE';
  448. *KTAB . 'OBJECTIF' = MOT 'J';
  449. KTAB .'PSI' = psi1;
  450. KTAB .'PHI' = phi1;
  451. KTAB .'FRONT_FISSURE' = sfron1 ;
  452. KTAB .'MODELE' = mod1tot;
  453. KTAB .'CARACTERISTIQUES' = mat1tot;
  454. KTAB .'SOLUTION_RESO' = u1;
  455. KTAB .'CHARGEMENTS_MECANIQUES' = for1tot;
  456. KTAB . 'POINT_CENTRE' = PA00;
  457. KTAB . 'COUCHE' = 2;
  458.  
  459. G_THETA KTAB;
  460.  
  461. Kchpo1 = KTAB . 'CHPO_RESULTATS' ;
  462.  
  463.  
  464. K1ev = EVOL BLEU 'CHPO' Kchpo1 'K1' sfron1;
  465. K2ev = EVOL DEFA 'CHPO' Kchpo1 'K2' sfron1;
  466. K3ev = EVOL ROUG 'CHPO' Kchpo1 'K3' sfron1;
  467.  
  468. S1pr = EXTR K1ev 'ABSC';
  469. K1pr = EXTR K1ev 'ORDO';
  470. K2pr = EXTR K2ev 'ORDO';
  471. K3pr = EXTR K3ev 'ORDO';
  472.  
  473. ************************************************************************
  474. * Formule analytique des FIC en milieu infinie *
  475. * pour une fissure circulaire inclinée *
  476. * "The stress analysis of crack Handook" H.TADA P.C0.Paris G.R.IRWIN *
  477. ************************************************************************
  478.  
  479. K1A = (2. * (-1.0)* szz0 * ((sin (-1.0 * beta0))**2) *
  480. (( a0 / pi )**0.5))+ ( S1pr * 0.);
  481.  
  482. CA = (4. * (-1.0)* szz0 * (sin (-1.0 * beta0)) * (cos (-1.0 * beta0)) *
  483. (( a0 / pi )**0.5)) / (2. - nu0);
  484.  
  485. K2A= CA * ( cos ((S1pr * 180. ) / (a0 * pi )));
  486.  
  487. K3A= CA * (1. - nu0)* ( sin ((S1pr * 180. ) / (a0 * pi )));
  488.  
  489. K1Aev = EVOL BLEU 'MANU' S1pr K1A;
  490. K2Aev = EVOL DEFA 'MANU' S1pr K2A;
  491. K3Aev = EVOL ROUG 'MANU' S1pr K3A;
  492. si (graph);
  493. ** Calcul des FIC par methode integrale G_THETAsur l autre demi front **
  494. KSTAB = tabl ;
  495. KSTAB . 'OBJECTIF' = MOT 'DECOUPLAGE';
  496. * KSTAB . 'OBJECTIF' = MOT 'J';
  497. KSTAB .'PSI' = psi1;
  498. KSTAB .'PHI' = phi1;
  499. KSTAB .'FRONT_FISSURE' = sfron2 ;
  500. KSTAB .'MODELE' = mod1tot;
  501. KSTAB .'CARACTERISTIQUES' = mat1tot;
  502. KSTAB .'SOLUTION_RESO' = u1;
  503. KSTAB .'CHARGEMENTS_MECANIQUES' = for1tot;
  504. KSTAB . 'POINT_CENTRE' = PA00;
  505. KSTAB . 'COUCHE' = 2;
  506.  
  507. G_THETA KSTAB;
  508.  
  509. KSchpo1 = KSTAB . 'CHPO_RESULTATS' ;
  510.  
  511. KS1ev = EVOL VERT 'CHPO' KSchpo1 'K1' sfron2;
  512. KS2ev = EVOL ORAN 'CHPO' KSchpo1 'K2' sfron2;
  513. KS3ev = EVOL ROSE 'CHPO' KSchpo1 'K3' sfron2;
  514.  
  515. SS1pr = EXTR KS1ev 'ABSC';
  516. KS1pr = EXTR KS1ev 'ORDO';
  517. KS2pr = EXTR KS2ev 'ORDO';
  518. KS3pr = EXTR KS3ev 'ORDO';
  519.  
  520. KS1A = (2. * (-1.0)* szz0 * ((sin (-1.0 * beta0))**2) *
  521. (( a0 / pi )**0.5))+ ( S1pr * 0.);
  522.  
  523. CSA = (4. * (-1.0)* szz0 * (sin (-1.0 * beta0)) *
  524. (cos (-1.0 * beta0)) * (( a0 / pi )**0.5)) / (2. - nu0);
  525.  
  526. KS2A= CSA * ( cos ((S1pr * 180. ) / (a0 * pi ) + 180.));
  527.  
  528. KS3A= CSA * (1. - nu0)* ( sin ((SS1pr * 180. ) / (a0 * pi ) + 180. ));
  529.  
  530. KS1Aev = EVOL VERT 'MANU' SS1pr KS1A;
  531. KS2Aev = EVOL ORAN 'MANU' SS1pr KS2A;
  532. KS3Aev = EVOL ROSE 'MANU' SS1pr KS3A;
  533.  
  534.  
  535. tdess1 = tabl;
  536. tdess1.'TITRE'= TABL;
  537. tdess1.'TITRE' . 1= MOT 'KI sim';
  538. tdess1.'TITRE' . 2= MOT 'KII sim ';
  539. tdess1.'TITRE' . 3= MOT 'KIII sim ';
  540. tdess1.'TITRE' . 4= MOT 'KI ana';
  541. tdess1.'TITRE' . 5= MOT 'KII ana ';
  542. tdess1.'TITRE' . 6= MOT 'KIII ana ';
  543. tdess1 . 4 = 'MOT' ' TIRR ';
  544. tdess1 . 5 = 'MOT' ' TIRR ';
  545. tdess1 . 6 = 'MOT' 'TIRR ';
  546. dess (K1ev et K2ev et K3ev et K1Aev et K2Aev et K3Aev )
  547. 'LEGE' tdess1 'TITR' 'SIF sur un demi front';
  548. dess (KS1ev et KS2ev et KS3ev et KS1Aev et KS2Aev et KS3Aev)
  549. 'LEGE' tdess1 'TITR' 'SIF sur l autre demi front';
  550. fins;
  551. ************************************************************************
  552. * Erreur en norme L2 des FIC le long du front *
  553. * (On s'eloigne des bords du demi front sur lequel on a fait *
  554. * le calcul g-theta pour ignorer les effet de bords dues a G-theta) *
  555. ************************************************************************
  556. si (complet);
  557. dint = 17;
  558. fint = 92;
  559. ip6 = 55;
  560. ip5 = 28;
  561. sinon ;
  562. dint = 6;
  563. fint = 51;
  564. ip6 = 29;
  565. ip5 = 15;
  566. finsi;
  567.  
  568. DiffL2K1 =((K1Aev - K1ev)**2) / ((EXTR K1A ip6)**2);
  569. DiffL2K2 =((K2Aev - K2ev)**2) / ((EXTR K2A 1)**2);
  570. DiffL2K3 =((K3Aev - K3ev)**2) / ((EXTR K3A ip6)**2);
  571.  
  572. ErrL2K1 =((intg DiffL2K1 dint fint)/
  573. ((EXTR S1pr fint ) -(EXTR S1pr dint)))**0.5;
  574. ErrL2K2 =((intg DiffL2K2 dint fint)/
  575. ((EXTR S1pr fint ) - (EXTR S1pr dint)))**0.5;
  576. ErrL2K3 =((intg DiffL2K3 dint fint)/
  577. ((EXTR S1pr fint ) -(EXTR S1pr dint)))**0.5;
  578. ErrL2K1 = extr ErrL2K1 1;
  579. ErrL2K2 = extr ErrL2K2 1;
  580. ErrL2K3 = extr ErrL2K3 1;
  581. list ErrL2K1;
  582. list ErrL2K2;
  583. list ErrL2K3;
  584.  
  585.  
  586. si (graph);
  587. opti donn 5;
  588. finsi;
  589.  
  590. SI ((ErrL2K1 < (0.05)) et (ErrL2K2 < (0.10)) et (ErrL2K3 < 0.15));
  591. 'MESS' ' ----------------------' ;
  592. 'MESS' ' SUCCES DU CAS-TEST !' ;
  593. 'MESS' ' ----------------------' ;
  594. Erre 0;
  595. Sinon;
  596. 'MESS' ' ---------------------' ;
  597. 'MESS' ' ECHEC DU CAS-TEST !' ;
  598. 'MESS' ' ---------------------' ;
  599. Erre 5;
  600. finsi;
  601. si (graph);
  602. ind1K = inde (KTAB . 'CHAMP_THET1');
  603. teta1 = KTAB . 'CHAMP_THET1' . (ind1K . 1);
  604. vtheta1 = vect teta1 'BLEU' ;
  605. trac vtheta1 (vol1 et scrack0 ) 'COUP' pa00 vzL1 pa03 'CACH';
  606. finsi;
  607.  
  608. TEMP 'IMPR' 'MAXI' 'CPU';
  609. FIN ;
  610.  
  611.  
  612.  
  613.  
  614.  
  615.  
  616.  
  617.  
  618.  
  619.  

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