Télécharger nlsb_operateur.dgibi

Retour à la liste

Numérotation des lignes :

  1. ***********************************************************************
  2. *_______________________________________________________________________
  3. *
  4. * Simulation COMPACT TENSION TEST
  5. *_______________________________________________________________________
  6. *-----------------------------------------------------------------------
  7. * ^ ^ ^ ^ ^ ^
  8. * | | | | | |
  9. * ________________
  10. * | |
  11. * | |
  12. * | |
  13. * | |
  14. * | |
  15. * | |
  16. * >|________________|
  17. * ^^^^^^^^^
  18. *-----------------------------------------------------------------------
  19. *
  20. * Analyse de la description du champs nonlocal en pointe d'entaille
  21. * C. Giry, F. Dufour, J. Mazars
  22. * -----------------------------
  23. ************************************************************************
  24. graph = 'N';
  25.  
  26. opti dime 2 elem qua4 mode plan cont;
  27.  
  28. * Taille element
  29. * --------------
  30. telem = 0.02;
  31.  
  32.  
  33. dens telem;
  34.  
  35. *-----------------------------------------------------------------------
  36. * Maillage
  37. *-----------------------------------------------------------------------
  38.  
  39.  
  40. *----- Largeur de l'éprouvette
  41. b = 1.;
  42.  
  43. *------Hauteur de l'éprouvette
  44. h = 1.;
  45.  
  46. p1 = 0. 0.;
  47. p1bis = (b/2.) 0.;
  48. p2 = b 0.;
  49. p3 = b h;
  50. p3bis = (b/2.) h;
  51. p4 = 0. h;
  52. d1 = d p1 p1bis ;
  53. d1bis = d p1bis p2 ;
  54. d2 = d p2 p3 ;
  55. d2bis = d p3 p3bis ;
  56. d3 = d p3bis p4 ;
  57. d4 = d p4 p1 ;
  58. s1 = (d1 et d1bis et d2 et d2bis et d3 et d4);
  59. s1 = surf s1 ;
  60.  
  61. *-----------------------------------------------------------------------
  62. * Modèle et matériau
  63. *-----------------------------------------------------------------------
  64.  
  65. MOD1 = MODE S1 MECANIQUE ELASTIQUE ISOTROPE ;
  66.  
  67. *
  68. * Paramètres
  69. * ----------
  70. * ft : limite de traction
  71. * younc : module d'Young
  72. * nuc : coefficient de Poisson
  73. * y0 : seuil en deformation de premiere fissuration
  74. * ll : longueur interne
  75.  
  76. ft = 3.E+6;
  77. younc = 3.D+9 ;
  78. nuc = 0.2;
  79. y0 = ft/younc;
  80. ll = 0.3 ;
  81.  
  82.  
  83. MAT1 = MATE mod1 'YOUN' younc 'NU' nuc 'RHO' 2.5e+3;
  84.  
  85. *-----------------------------------------------------------------------
  86. * Connectivites
  87. *-----------------------------------------------------------------------
  88. conn1 = CONN mod1 ll 'NORMAL' 'INTERIEUR' ;
  89. conn2 = CONN mod1 ll 'DROITE' P1bis P2 'CENTRE' ;
  90. connt = conn1 et conn2;
  91.  
  92. *-----------------------------------------------------------------------
  93. * Conditions limites
  94. *-----------------------------------------------------------------------
  95. cl1 = bloq (d1bis) UY;
  96. cl2 = bloq p4 UX;
  97. CL3 = RELA 'ENSE' UY (d2bis et d3);
  98. CLTOT = CL1 ET CL2 ET CL3;
  99. forc1 = FORC (d2bis et d3) (0. 1.) ;
  100.  
  101. *-----------------------------------------------------------------------
  102. * Rigidité
  103. *-----------------------------------------------------------------------
  104. rbet = RIGI mod1 mat1 ;
  105. rigt = rbet ET CLTOT ;
  106.  
  107. *-----------------------------------------------------------------------
  108. * Analyse elastique
  109. *-----------------------------------------------------------------------
  110. depl1 = RESO rigt forc1 ;
  111.  
  112. cont1 = SIGM mod1 mat1 depl1;
  113. conp1 = PRIN mod1 cont1;
  114.  
  115. contmax = MAXI (EXCO conp1 'SI11');
  116. f2 = ft/contmax;
  117. forc1 = FORC (d2bis et d3) (0. f2) ;
  118. depl1 = RESO rigt forc1 ;
  119.  
  120. *-----------------------------------------------------------------------
  121. * Analyse du champ nonlocal
  122. *-----------------------------------------------------------------------
  123. defo1 = EPSI mod1 depl1;
  124. cont1 = SIGM mod1 mat1 depl1;
  125. defp1 = PRIN mod1 defo1;
  126. conp1 = PRIN mod1 cont1;
  127.  
  128.  
  129. * Deformation equivalente Mazars
  130. * ------------------------------
  131. epeq11 = EXCO defp1 'SI11';
  132. epeq12 = MASQ epeq11 'SUPERIEUR' 0.;
  133. epeq1 = NOMC 'EPTI' (epeq11 * epeq12);
  134. epeq21 = EXCO defp1 'SI22';
  135. epeq22 = MASQ epeq21 'SUPERIEUR' 0.;
  136. epeq2 = NOMC 'EPTI' (epeq21 * epeq22);
  137. epeq31 = EXCO defp1 'SI33';
  138. epeq32 = MASQ epeq31 'SUPERIEUR' 0.;
  139. epeq3 = NOMC 'EPTI' (epeq31 * epeq32);
  140.  
  141. epeqm = ((epeq1)**2) + ((epeq2)**2) + ((epeq3)**2);
  142. epeqm = (epeqm)**(0.5);
  143.  
  144. * Calcul de la longueur par element et limite de traction
  145. * -------------------------------------------------------
  146. CHAMLO = MANU 'CHML' MOD1 'LONG' 1.0 'STRESSES';
  147. CHAML1 = INTG CHAMLO MOD1 ELEM 'LONG';
  148. CHAML2 = CHAML1**(1./2.);
  149. CHAML3 = NOMC 'LONG' (CHAN 'STRESSES' MOD1 CHAML2);
  150. CHAMFT = MANU 'CHML' MOD1 'FT' ft 'STRESSES';
  151.  
  152. * Calcul du champ nonlocal de deformation equivalente
  153. * ---------------------------------------------------
  154. lmot1 = MOTS 'EPTI' 'SI11' 'SI22' 'SI33' 'COX1' 'COY1'
  155. 'COZ1' 'COX2' 'COY2' 'COZ2' 'COX3' 'COY3'
  156. 'COZ3' 'FT' 'LONG';
  157. epeqmsb = epeqm + conp1 + chamft + chaml3;
  158. epeqnlm = 'NLOC' 'SB' epeqmsb connt lmot1;
  159. Pmax1 = POIN (EXCO epeqnlm 'EPTI') 'MAXI';
  160. XX1 YY1 = COOR PMAX1;
  161. lerr = ABS((MAXI XX1) - (0.5*b));
  162.  
  163. list lerr;
  164.  
  165. si (ega graph 'O');
  166. TRAC epeqnlm mod1;
  167. finsi;
  168.  
  169. * Test sur la position du maximum du champ nonlocal par rapport a la
  170. * pointe de l'entaille
  171. * --------------------
  172.  
  173. si (lerr > telem);
  174. erre 5;
  175. sinon;
  176. erre 0;
  177. finsi;
  178.  
  179.  
  180. fin;
  181.  
  182.  
  183.  
  184.  
  185.  
  186.  
  187.  
  188.  

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