Télécharger la plaquette Cast3M


fr en

Télécharger xfem3d_01.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : xfem3d_01.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *******************************************************
  5. *opti echo 0;
  6.  
  7. MESS '' ;
  8. MESS '  CALCUL ELASTIQUE DE PROPAGATION STATIQUE DE FISSURE ';
  9. MESS '                    en 3D XFEM ' ;
  10. MESS '' ;
  11.  
  12.  
  13. *******************************************************
  14. *** Options  
  15.  
  16. *** options de calcul 
  17. opti dime 3 elem CUB8 mode trid;
  18. *** options de trace
  19. graph = faux ;
  20. * graph = vrai ;
  21. opti 'TRAC' PSC ;
  22. * opti 'TRAC' X ;
  23.  
  24.  
  25. *******************************************************
  26. *** Maillage : cube sain 1 x 1 x 1
  27.  
  28. * finesse de la discretisation
  29. * n1 = 36;
  30. * n1 = 24;
  31. * n1 = 12;
  32. n1 = 8;
  33.  
  34. * geoemtrie
  35. L1 = 1. ;
  36. dx1  = (L1) / (flot n1) ;
  37. isov1 = dx1 * (prog -6. PAS 0.5 6.);
  38.  
  39. * face avant (a)
  40. DENS (dx1) ;
  41. pa1 = (-0.5*L1) 0. (-0.5*L1) ;
  42. pa2 = ( 0.5*L1) 0. (-0.5*L1) ;
  43. pa3 = ( 0.5*L1) 0. ( 0.5*L1) ;
  44. pa4 = (-0.5*L1) 0. ( 0.5*L1) ;
  45. dra1 = pa1 droi n1 pa2;
  46.  
  47. * repartition 0.25 // 0.5 // 0.25 
  48. vz1 = 0. 0. (0.25*L1);
  49. sa1 = dra1 TRAN (n1/4) vz1;
  50. dra10 = sa1 cote 3;
  51. sa2 = dra10 TRAN (n1/2) (2.*vz1);
  52. dra20 = sa2 cote 3;
  53. sa3 = dra20 TRAN (n1/4) vz1;
  54. dra3 = sa3 cote 3;
  55.  
  56. * volume = cube
  57. ex = 1. 0. 0.;
  58. ey = 0. 1. 0.;
  59. ez = 0. 0. 1.;
  60. ny1 = n1;
  61.  
  62. vol1 = sa2 VOLU 'TRAN' ny1 (L1*ey);
  63. vol2 = sa1 VOLU 'TRAN' ny1 (L1*ey);
  64. vol3 = sa3 VOLU 'TRAN' ny1 (L1*ey);
  65. dx1elim = 1.E-5*dx1;
  66. ELIM vol1 vol2 dx1elim;
  67. ELIM vol1 vol3 dx1elim;
  68. vol23 = vol2 et vol3;
  69. vol0 = vol1 et vol23;
  70.  
  71. * qq recup 
  72. pa0 = poin vol1 'PROCH' (0. 0. 0.);
  73. pa0x = poin vol1 'PROCH' (0. L1 0.);
  74. pa0z = (poin vol1 'PROCH' ((-0.5*L1) L1 0.))
  75.    et  (poin vol1 'PROCH' ((0.5*L1) L1 0.));
  76. sb1 = (dra3 TRAN ny1 (L1*ey)) coul 'BLEU';
  77. ELIM sb1 vol3 (1.E-4 * dx1);
  78. sc1 = (dra1 TRAN ny1 (L1*ey)) coul 'OCEA';
  79. ELIM sc1 vol2 (1.E-4 * dx1);
  80.  
  81. con0 = aret vol0;
  82. env0 = ENVE vol0;
  83.  
  84. eye1 = 100. -50. 30. ;
  85. eye2 = 0. -0. 100. ;
  86. si(graph);
  87.   trac eye1 (vol0 et sb1 et sc1) 'CACH';
  88. fins;
  89. mess (nbno vol23) (nbno vol1) (nbno vol0); 
  90. mess (nbel vol23) (nbel vol1) (nbel vol0);
  91.  
  92.  
  93. **********************************************************
  94. * tables et indice de propagation
  95. it = 0;
  96. pcrack1 = tabl ;
  97. crack1 = tabl;
  98. che1X = tabl;
  99. Rel1X = tabl;
  100.  
  101.  
  102. ***********************************************************
  103. *** Maillage de la fissure
  104.  
  105. *** fissure semi elliptique aplatie
  106.  
  107.   a0 = 0.50;  
  108.   rapab0 = 0.2;
  109.   b0 = a0 * rapab0;
  110.   bbord = 0.20;
  111. * finesse du maillage de la fissure
  112.   da0 = dx1 / 4.; 
  113.   na0 = enti (a0 / da0);
  114.   nb0 = enti (b0 / da0);
  115.   &ELEM = VALE 'ELEM';
  116.   opti elem TRI3;
  117. *  da000 = rapab0 * da0;
  118.   da000 =  da0;
  119. * points de base
  120.   DENS (da000);
  121.   pf1 =      a0   0.     0. ;
  122.   pf2 =      a0   bbord  0. ;
  123.   pf4 = (-1.*a0)  bbord  0. ;
  124.   pf5 = (-1.*a0)  0.     0. ;
  125.   DENS (da0);
  126.   pf0 = 0. 0. 0.;
  127.   pf3 = 0. bbord 0.;
  128. * lignes de base
  129.   lf1 = (pf5 droi pf0) droi pf1;
  130.   lf2 = (pf4 droi pf3) droi pf2;
  131. * profil elliptique
  132.   un2 = manu lf2 'CHPO' 'SCAL' 1.;
  133.   x2  = coor lf2 1;
  134.   vy2 = b0 * ((un2 - ((x2/a0)**2))**0.5);
  135.   vy2 = vy2 nomc 'UY';
  136.   lfron1 = (lf2 plus vy2) coul 'ORAN';
  137. * surface reglee 
  138.   sfron1 = (regl lfron1 lf1) coul 'ROUG';
  139.  
  140. * tres important 
  142.  
  143. si(graph);
  144.   trac eye2 sfron1;
  145.   trac eye1 (con0 et sfron1 ) 'CACH';
  146.   trac eye1 (con0 et lfron1) 'CACH';
  147. fins;
  148.  
  149. pcrack1.it = lfron1;
  150. crack1.it =  sfron1 ;
  151.  
  152. psi1 phi1 = PSIP vol1 crack1.it  'DEUX'  pcrack1.it;
  153.  
  154. si(graph);
  156.   trac phi1 vol0 isov1 'COUP' pf0 (0. 0. 1.) pf3;
  157.   trac psi1 vol0 isov1 'COUP' pf0 pf1 pf3;
  158. fins;
  159.  
  160.  
  161. *opti donn 5 ;
  162. **********************************************************
  163. ***   MODELE & MATERIAU   ***
  164. **********************************************************
  165.  
  166. * coef elastique
  167.   Ey0 = 2.E11;
  168.   nu0 = 0.3 ;
  169.   rho0= 7800. ;
  170.   mu0 = (0.5 * Ey0) / (1. + nu0) ;
  171.   Estar0 = Ey0 / (1. - (nu0**2));
  172.  
  173. * Loi de Paris  (on convertit C : mm -> m)
  174.   mParis = 3.;
  175.   CParis = 1.E-9 * 1.E-3;
  176.  
  177.  
  178. **********************************************************
  179. *** Modele & materiau  ***
  180. *
  181. * zone de propagation (X-FEM)
  182.   mod1 = MODE vol1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 'XC8R';
  183.   mat1 = MATE mod1 'YOUN' Ey0 'NU' nu0 'RHO' rho0  ;
  184.   che1X.it = TRIE mod1 psi1 phi1 ;
  185.  
  186. * constructionsion des blocages des ddl X-fem non actifs dans 
  187. * les éléments de transition.
  188.   Rel1x . it = RELA mod1;
  189. *
  190. * elements standards
  191.   mod23 = MODE vol23 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 'CUB8';
  192.   mat23 = MATE mod23 'YOUN' Ey0 'NU' nu0 'RHO' rho0 ;
  193. *
  194.   mod1tot = mod1 et mod23;
  195.   mat1tot = mat1 et mat23 ;
  196.  
  197.  
  198.  
  199. **********************************************************
  200. ***   CL et DEPLACEMENTS IMPOSES   ***
  201. **********************************************************
  202. syy0 = -100.E6;
  203.  
  204.  
  205. *** MODE I ***
  206. * CL mode I
  207.   cl1 = bloq pa0 'UX' 'UY' 'UZ'; 
  208.   cl2 = bloq pa0x 'UZ';
  209.   cl3 = bloq pa0z 'UX' 'UZ';
  210.   cl1tot =  (cl1 ET cl2 ET cl3) ;
  211. * chargement mode I
  212.   for1 = PRES 'MASS' mod23 sb1 syy0;
  213.   for2 = PRES 'MASS' mod23 sc1 syy0;
  214.   for1tot = for1 et for2;
  215.  
  216. si(graph);
  217.   trac (vect for1tot  'FORC' 'BLEU') 
  218.       (con0 et pcrack1.it) 'CACH';
  219. fins;
  220.  
  221.  
  222. *** matrice elastique et assemblage
  223. Rig1tot = RIGI mod1tot mat1tot;
  224. K1tot = Rig1tot et cl1tot et (Rel1x . it);
  225.  
  226.  
  227. *******************************************************
  228. * INITIALISATIONS POUR LA BOUCLE TEMPORELLE
  229. *******************************************************
  230. logda = FAUX;
  231. u1 = tabl;
  232. u1 . 0 = RESO K1tot (0.*for1tot);
  233.  
  234. TKchpo1 = tabl;  TK1ev = tabl; TK2ev = tabl; TK3ev = tabl;
  235. *******************************************************
  236. * BOUCLE TEMPORELLE POUR FAIRE PROPAGER LA FISSURE
  237. *******************************************************
  238. *nt1=n1;
  239. nt1=2;
  240.  
  241.  REPE BT1 nt1;
  242.  
  243.   it = it + 1;
  245.   mess '***********************' it '*************************';
  246.   TITR (chai 'iteration ' it);
  247.  
  248. *******************************************************
  249. *         PROPAGATION ÉLÉMENTAIRE DE FISSURE         
  250. *******************************************************
  251. si(logda);
  252.  
  253. *** Angle de propagation ***
  254.   un1 = manu 'CHPO' pcrack1 . (it - 1)  'SCAL' 1.;
  255. * angle
  256.   sgnK2 = (masq K2G 'EGSUPE' 0.) - (masq K2G 'INFERIEUR' 0.);
  257. * modif pour etre tranquille
  258.   masK2 = masq (abs K2G) 'EGINFE' (1.E-4*(abs K1G));
  259.   K2G = (masK2) + ((un1 - masK2) * K2G);
  260.     tet0c = K1G / K2G;
  261.     tet8c = (tet0c**2) + (8. * un1);
  262.     tet1c = tet0c - ( sgnK2 * (tet8c**0.5) );
  263.     tet1c = 2. * (atg (0.25 * tet1c));
  264.     tet1c = (0. * masK2 ) + ((un1 - masK2) * tet1c);
  265. * on recupere la direction du repere local de la fissure (CHPOINT)
  266.   ETIP1 = KTAB . 'UTILTET1' . 'DIRECTION1';
  267.   ETIP2 = KTAB . 'UTILTET1' . 'DIRECTION2';
  268. * on fait tourner de l'angle critique 
  269.   DA1VEC = ((cos tet1c) * ETIP1) + ((sin tet1c) * ETIP2);
  270. *trac (vect DA1VEC) (pcrack1 . it et con0);
  271.  
  272. *** vitesse de propagation ***  
  273.   K1eff = ( (K1G**2) + (K2G**2) )**0.5;
  274.   K1eff = ( (K1G**2) + (K2G**2) + ((Estar0/mu0)*(K3G**2)) )**0.5;
  275. * Passage en MPa.m^0.5
  276.   K1eff = K1eff * 1.E-6;
  277. * quel est le dn correspondant a da1min ?
  278.   K1min= mini K1eff 'ABS';
  279.   dadn = CParis * (K1min**mParis);
  280.   dn = da1min / dadn;
  281.   dadn = (CParis * dn) * (K1eff**mParis) ;
  282.   DA1VEC = dadn * DA1VEC;
  283.  
  284. *** on cree un maillage par translation ****
  285.   trac (vect DA1VEC 1) (pcrack1 . (it - 1) et con0); 
  286.   pcrk1 = pcrack1 . (it - 1) PLUS DA1VEC;  
  287.  
  288. *** quelques chpo utiles ***
  289. * on calcul la distance signee a la paroi 
  290.   dis0 = PSIP pcrack1 . (it - 1) env0;
  291.   dis1 = PSIP pcrk1              env0;
  292. * on cree un chpoint permettant de savoir si un point est a l interieur 
  293.   unvol1 = manu 'CHML' vol0 'SCAL' 1.;
  294.   uncrk1 = PROI pcrk1 unvol1;
  295.  
  296. *** Etape 1 : traitement des segments "a cheval" sur le bord  ***
  297. * boucle sur les segments de pcrk1 
  298.   dismall1 = 1.E-6 * da0;
  299.   ncrk1 = nbel pcrk1;
  300.   icrk1 = 0;
  301.   repe Bcrk1 ncrk1;
  302.     icrk1 = icrk1 + 1;
  303. *   recup des segments et des points    
  304.     pcrk1i = ELEM pcrk1              icrk1;
  305.     pt1a = poin pcrk1i 'INITIAL';    pt1b = poin pcrk1i 'FINAL';
  306. *   test pour savoir si a cheval
  307.     un1a = extr uncrk1 pt1a 'SCAL';  un1b = extr uncrk1 pt1b 'SCAL';
  308.     si((abs(un1a - un1b)) >eg 0.5);
  309.     phi1a = EXTR dis1 pt1a 'PHI';    phi1b = EXTR dis1 pt1b 'PHI';
  310.     si(un1a < 0.5);
  311. *   il faut ramener le point pt1a    
  312.        v1a = phi1a / (phi1a - phi1b);
  313.        v1a = v1a * (pt1b moin pt1a);
  314.        DEPL pt1a 'PLUS' v1a;
  315.     sino;
  316. *   il faut ramener le point pt1b    
  317.        v1b = phi1b / (phi1b - phi1a);
  318.        v1b = v1b * (pt1a moin pt1b);
  319.        DEPL pt1b 'PLUS' v1b;
  320.     fins;
  321.     fins;
  322.   fin Bcrk1;
  323.  
  324. *** Etape 2 : traitement des points encore a l'exterieur *** 
  325. * mise a jour de uncrk1
  326.   uncrk1 = PROI pcrk1 unvol1;
  327.   icrk1 = 0;
  328.   repe Bcrk2 ncrk1;
  329.     icrk1 = icrk1 + 1;
  330. *   recup des segments et des points    
  331.     pcrk1i = ELEM pcrk1              icrk1;
  332.     pt1a = poin pcrk1i 'INITIAL'; 
  333. *   pt1a a l exterieur ?
  334.     un1a = EXTR uncrk1 pt1a 'SCAL';
  335.     si(un1a < 0.5);
  336.       pcrk0i = ELEM pcrack1 . (it - 1) icrk1;
  337.       pt0a = poin pcrk0i 'INITIAL'; 
  338. *     distance signée associée
  339.       phi0a = EXTR dis0 pt0a 'PHI';
  340.       phi1a = EXTR dis1 pt1a 'PHI';
  341. *     il faut ramener le point pt1a    
  342.       v1a = phi1a / (phi0a - phi1a);
  343.       v1a = v1a * (pt1a moin pt0a);
  344.       DEPL pt1a 'PLUS' v1a;
  345.     fins;
  346.   fin Bcrk2;
  347. * sans oublier le dernier point du dernier segment
  348.     pt1a = poin pcrk1i 'FINAL'; 
  349. *   pt1a a l exterieur ?
  350.     un1a = EXTR uncrk1 pt1a 'SCAL';
  351.     si(un1a < 0.5);
  352. *     recup des segments et des points du front precedent
  353.       pcrk0i = ELEM pcrack1 . (it - 1) ncrk1;
  354.       pt0a = poin pcrk0i 'FINAL'; 
  355. *     distance signée associée
  356.       phi0a = EXTR dis0 pt0a 'PHI';
  357.       phi1a = EXTR dis1 pt1a 'PHI';
  358. *     il faut ramener le point pt1a    
  359.       v1a = phi1a / (phi0a - phi1a);
  360.       v1a = v1a * (pt1a moin pt0a);
  361.       DEPL pt1a 'PLUS' v1a;
  362.     fins;
  363.  
  364. *** si tous les points du front sont sortis 
  365. *   c est qu on a tout cassé !!! => on quitte
  366.    meshun1 = extr uncrk1 'MAILLAGE';
  367.    si(ega (TYPE meshun1) 'ENTIER');   quit BT1;   fins;
  368.  
  369. *** Etape 3 : points sur le bord qui sont entrés a linterieur *** 
  370. * mise a jour de dis1
  371.   dis1 = PSIP pcrk1 env0;
  372. * 1er point du 1er segment  
  373.     pcrk0i = ELEM pcrack1 . (it - 1) 1;
  374.     pt0a = poin pcrk0i 'INITIAL';
  375.     phi0a = EXTR dis0 pt0a 'PHI';
  376.     si((abs phi0a) <eg dismall1);
  377.       pcrk1i = ELEM pcrk1 1;
  378.       pt1a = poin pcrk1i 'INITIAL';
  379.       phi1a = EXTR dis1 pt1a 'PHI';
  380.       si((abs phi1a) > dismall1);     
  381.         pt0b = poin pcrk0i 'FINAL';
  382.         phi0b= EXTR dis0 pt0a 'PHI';
  383.         pt1b = poin pcrk1i 'FINAL';
  384.         phi1b= EXTR dis1 pt1b 'PHI';
  385.         v1a = phi1a / (phi1b - phi1a) ;
  386.         v1a = v1a * (pt1a moin pt1b);
  387.         DEPL pt1a 'PLUS' v1a;
  388.       fins;
  389.     fins;
  390. * dernier point du dernier segment  
  391.     pcrk0i = ELEM pcrack1 . (it - 1) ncrk1;
  392.     pt0b = poin pcrk0i 'FINAL';
  393.     phi0b= EXTR dis0 pt0a 'PHI';
  394.     si((abs phi0b) <eg dismall1);
  395.       pcrk1i = ELEM pcrk1 ncrk1;
  396.       pt1b = poin pcrk1i 'FINAL';
  397.       phi1b= EXTR dis1 pt1b 'PHI';
  398.       si((abs phi1b) > dismall1);     
  399.         pt0a = poin pcrk0i 'INITIAL';
  400.         phi0a = EXTR dis0 pt0a 'PHI';
  401.         pt1a = poin pcrk1i 'INITIAL';
  402.         phi1a = EXTR dis1 pt1a 'PHI';
  403.         v1b = phi1b / (phi1a - phi1b) ;
  404.         v1b = v1b * (pt1b moin pt1a);
  405.         DEPL pt1b 'PLUS' v1b;
  406.       fins;
  407.     fins;
  408.  
  409. *** on recupere le vrai nouveau front ***
  410.   pcrack1 . it = ELEM pcrk1 'APPUYE' 'STRICT' meshun1;
  411.   si((flot (it / 2)) ega ((flot it) / 2.));
  412.     pcrack1 . it = pcrack1 . it coul 'ROUG';
  413.   sinon;
  414.     pcrack1 . it = pcrack1 . it coul 'ORAN';
  415.   fins; 
  416.  
  417. *** on maille la surface fissuree ***
  418.   &ELEM = VALE 'ELEM';
  419.   opti 'ELEM' TRI3;
  420.   nregl1 = 1 + (enti (da1min / da0));
  421.   sfiss1 = regl pcrk1 nregl1 pcrack1 . (it - 1) ;
  422. *** on retraite la surface fissuree cree au debut ***
  423.   ELIM sfiss1 dismall1;
  424.   sfiss1 = (REGE sfiss1) elem 'TRI3';
  425.   crack1 . it = sfiss1 et crack1 . (it - 1) ; 
  427.  
  428. *** fin de l avancee du maillage du front **********************
  429.  
  430.  
  431. *** Actualisations ***
  432. * ...des level set
  433.   psi1 phi1 = PSIP vol1 crack1.it  'DEUX'  pcrack1.it;
  434.   si(graph);
  435.     trac phi1 vol0 isov1 'COUP' pf0 (0. 0. 1.) pf3;
  436.     trac psi1 vol0 isov1 'COUP' pf0 pf1 pf3;
  437.   fins;
  438.  
  439. * ...du modele, de la rigidite,... (inutile pour le materiau)
  440.   Che1X . it = TRIE mod1tot psi1 phi1;
  441.  
  442. * ... des blocages des ddl X-fem non actifs dans 
  443. * les éléments de transition.
  444.   Rel1x . it = RELA mod1tot;
  445.   Rig1tot = RIGI mod1tot mat1tot;
  446.   K1tot = Rig1tot et cl1tot et (Rel1x . it) ;
  447.  
  448.  
  449. *******************************************************
  450. *           PAS DE PROPA 
  451. *******************************************************
  452. sino;
  453.  
  454.     crack1 . it = crack1 . (it - 1);
  455.     pcrack1 . it = pcrack1 . (it - 1);
  456.  
  457. fins;
  458.  
  459. trac (con0 et crack1 . it) eye2; 
  460.  
  461. ****************************
  462. *** resolution ***
  463.   u1 . it  = RESO K1tot for1tot;
  464.  
  465.   u1phy = XFEM 'RECO' u1 . it mod1tot ;
  466.   sig1  = SIGM mod1tot u1 . it mat1tot  line;
  467.   ucrk1u ucrk1d = XFEM 'FISS' crack1 . it  u1 . it mod1tot ;
  468.   si(graph);
  469. *    amp1 = 25.;
  470.     amp1 = 10.;
  471.     def1 = (DEFO u1phy vol0 amp1) ;
  472.     TRAC sig1 mod1tot def1 eye1;
  473.     def2 =  (DEFO u1phy  con0  amp1)
  474.          et (DEFO ucrk1u crack1 . it amp1)
  475.          et (DEFO ucrk1d crack1 . it amp1);
  476.     TRAC def2 eye1;
  477.     def3 =  def1
  478.          et (DEFO ucrk1u crack1 . it amp1)
  479.          et (DEFO ucrk1d crack1 . it amp1);
  480.     TRAC def3 eye1 'COUP' pf0 (0. 0. 1.) pf3 'CACH';
  481.   finsi;
  482.  
  483. ***   Calcul des FIC par methode integrale G_THETA   ***
  484.   KTAB = tabl ;
  485.   KTAB . 'OBJECTIF' = MOT 'DECOUPLAGE';
  486.   KTAB .'PSI' = psi1;
  487.   KTAB .'PHI' = phi1;
  488.   KTAB .'FRONT_FISSURE' = pcrack1 . it ; 
  489.   KTAB .'MODELE' = mod1tot;
  490.   KTAB .'CARACTERISTIQUES' = mat1tot;
  491.   KTAB .'SOLUTION_RESO' = u1 . it;
  492.   KTAB .'CHARGEMENTS_MECANIQUES' = for1tot;
  493.   KTAB . 'COUCHE' = 1;
  494. *  KTAB . 'COUCHE' = 2;
  495.  
  496.   G_THETA KTAB;
  497.  
  498.   Kchpo1 = KTAB . 'CHPO_RESULTATS' ;
  499.   TKchpo1 . it = Kchpo1;
  500.  
  501.   K1ev = EVOL BLEU 'CHPO' Kchpo1 'K1' pcrack1 . it;
  502.   K2ev = EVOL DEFA 'CHPO' Kchpo1 'K2' pcrack1 . it;
  503.   K3ev = EVOL ROUG 'CHPO' Kchpo1 'K3' pcrack1 . it;
  504.   TK1ev . it = K1ev;
  505.   TK2ev . it = K2ev;
  506.   TK3ev . it = K3ev;
  507.  
  508.   S1pr = EXTR K1ev 'ABSC';
  509.   K1pr = EXTR K1ev 'ORDO';
  510.   K2pr = EXTR K2ev 'ORDO';
  511.   K3pr = EXTR K3ev 'ORDO';
  512.  
  513. *  K1eqpr = ((K1pr**2) + (K2pr**2))**0.5; 
  514.   K1eqpr = ( (K1pr**2) + (K2pr**2) + ((Estar0/mu0)*(K3pr**2)) )**0.5; 
  515.  
  516.   K1eqev = EVOL 'VERT' 'MANU' 's' S1pr 'K1eq' K1eqpr;
  517.   si(it ega 1);
  518.     tdess1 = tabl;
  519.     tdess1 . 1 = 'MOT' 'MARQ PLUS TIRR ';
  520.     tdess1 . 2 = 'MOT' 'MARQ CROI TIRR ';
  521.     tdess1 . 3 = 'MOT' 'TIRM ';
  522.   fins;
  523.   dess (K1ev et K2ev et K3ev et K1eqev) tdess1 'TITR' 'SIF';
  524.  
  525.   dadnev = CParis * (K1eqev**mParis);
  526.   dess dadnev 'TITR' 'DADN';
  527.  
  528.   K1G = EXCO Kchpo1 'K1';
  529.   K2G = EXCO Kchpo1 'K2';
  530.   K3G = EXCO Kchpo1 'K3';
  531.  
  532.   GTAB = tabl ;
  533.   GTAB .'OBJECTIF' = MOT 'J   ';
  534.   GTAB .'PSI' = psi1;
  535.   GTAB .'PHI' = phi1;
  536.   GTAB .'FRONT_FISSURE' = pcrack1 . it ; 
  537.   GTAB .'MODELE' = mod1tot;
  538.   GTAB .'CARACTERISTIQUES' = mat1tot;
  539.   GTAB .'SOLUTION_RESO' = u1 . it;
  540.   GTAB .'CHARGEMENTS_MECANIQUES' = for1tot;
  541.   GTAB .'COUCHE' = 1;
  542.  
  543.   G_THETA GTAB;  
  544.   Jchpo1 = GTAB . CHPO_RESULTATS;
  545.   Jev1 = EVOL ROSE 'CHPO' Jchpo1 'J' pcrack1 . it;
  546.   dess Jev1;
  547.  
  548.  
  549. *** Increment de propagation ***
  550.   logda = VRAI;
  551. * increment manuel
  552.   da1min = 1.4 * dx1;
  553.  
  554.  
  555. * OPTI DONN 5 echo 1 trac X;
  556.  
  557. * ON VERIFIE QUE LA FORMULE D'IRWIN EST VERIFIE A LA TOLERANCE PRES
  558.     GEV1 = ((1.-(NU0**2))/EY0*((K1EV**2) + (K2EV**2))) + ((K3EV**2)/(2.*MU0)) ;
  559.     ERR1 = (INTG ((GEV1 - JEV1)**2))**0.5 ;
  560.     TOL1 = 0.05 ;
  561.     CRI1 = ((INTG (JEV1**2))**0.5) * TOL1 ;
  562.     SAUT 'LIGNE' ;
  563.     MESS 'ERREUR SUR LA FORMULE D''IRWIN :' ERR1 ;
  564.     MESS 'CRITERE' ' :'*31 CRI1 ;
  565.     SAUT 'LIGNE' ;
  566.     SI (ERR1 >EG CRI1) ;
  567.         MESS 'ERREUR : L''ERREUR SUR LA FORMULE D''IRWIN A L''ITERATION' ;
  568.         MESS '        ' &BT1 ' EST SUPERIEURE A' ' ' (ENTI (100.*TOL1)) '%' ;
  569.         ERRE 5 ;
  570.     FINSI ;
  571.  
  572.  FIN BT1;
  573.  
  574. FIN ;
  575.  
  576.  
  577.  
  578.  

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