Télécharger la plaquette Cast3M


fr en

Télécharger lispel.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  lispel.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. OPTION ECHO 0;
  5. GRAPH = 'N';
  6. SAUT PAGE;
  7. *_________________________________________________________
  8. *
  9. *  PLAQUE AVEC FISSURE SEMI-ELLIPTIQUE DEBOUCHANTE
  10. *             a/c = a/t = 0.2 
  11. *
  12. *  Comparaison des valeurs du facteur d'intensite de contrainte
  13. *  avec la solution de Raju-Newman tirée de RCC-MR A16 (dmt94-043) 
  14. *
  15. *  Dimensions plaque : largeur=longueur=160mm , epai 20mm
  16. *  Dimensions fissure : a=4mm  c=20mm
  17. *  Chargement : contrainte de membrane uniforme SIG=1MPa
  18. *_________________________________________________________
  19. *
  20. * Materiau:
  21. xYoung = 210000.;       
  22. *____________________________________________
  23. *
  24. *  MAILLAGE DE LA PLAQUE
  25. *____________________________________________
  26. *
  27. OPTION DIME 3 ELEM qua4;
  28.  
  29. * Geometrie:
  30. larg1 = 80.;long1 = 80.;
  31. ww  = 20.;flength = 20.0; fdepth  =  4.0;
  32.  
  33. * Nombre d'elements: 
  34. *nb1 = 10;nb2 = 15;nb3 = 16;
  35. nb = 15 ; nb1 = enti (nb / 2.) ;
  36.  
  37. dens (flength / (10. * nb1));
  38. p6  = 0.    flength 0.;
  39.  
  40. dens (flength / nb1);
  41. p0  = 0.    0.   0.;
  42.  
  43. dens (flength / (0.75 * nb1));
  44. p4  = (flength / 2)  0.    0.;
  45. p8  = (flength / 2) flength 0.;
  46.  
  47. dens (flength / (0.2 * nb1));
  48. P5  = long1        0.    0.;
  49. p7  = 0.           larg1   0.;
  50. p9  = long1        larg1   0.;
  51.  
  52. d1 = DROI (-1 * nb1) p0 p4;
  53. d2 = DROI (-1 * nb) p6 p0;
  54. d3 = DROI (-1 * nb1) p8 p6;
  55. d4 = DROI (-1 * nb) p4 p8;
  56. d5 = DROI (-1 * nb) p4 p5;
  57. d6 = DROI (-1 * nb) p8 p9;
  58. d7 = DROI (-1 * nb) p6 p7;
  59. d8 = DROI (-1 * nb) p5 p9;
  60. d9 = DROI (-1 * nb1) p9 p7;
  61.  
  62. su1 = DALL d1 d2 d3 d4 PLAN;
  63. su2 = DALL d5 d8 (INVE d6) (INVE d4) PLAN;
  64. su3 = DALL d6 d9 (INVE d7) (INVE d3) PLAN;
  65. su6 = su1 ET su2 ET su3;
  66.  
  67. d8a d2a d7a su8 = d8 d2 d7 su6 SYME 
  68.                 PLAN (0. 0. 0.) (0. 1. 0.) (0. 0. 1.);
  69. p06 = P6 SYME PLAN (0. 0. 0.) (0. 1. 0.) (0. 0. 1.);
  70.  
  71. su7 = su6 ET su8;
  72.  
  73. d7ap = chan poi1 d7a;ep06 = manu poi1 p06;
  74. ELIM .0001 d7ap ep06;
  75. ELIM .0001 su7 (d7ap diff ep06);
  76.  
  77. su9 = ORIENT su7 DIRECTION (0. 0. 10.);      
  78.  
  79. *_____________________________________________________________________
  80. *
  81. *--- modele et caracteristiques du materiau, coque
  82. mod1 =MODE su9 MECANIQUE ELASTIQUE coq4;
  83. mat1 =MATE mod1 YOUN xYoung NU .3 EPAI ww;
  84.  
  85. *_____________________________________________________________________
  86. *
  87. *               MAILLAGE DE LA FISSURE
  88. *_____________________________________________________________________
  89. *
  90. *--- modele et caracteristiques du materiau, fissure
  91.  
  92. d2    = INVE d2; d2a   = INVE d2a;
  93.  
  94. ii = 0;REPETER matFiss (NBEL d2);ii = ii+1;
  95.  
  96.   pp1a = d2  POIN ii;
  97.   pp2a = d2  POIN (ii+1);
  98.   pp1b = d2a POIN ii;
  99.   pp2b = d2a POIN (ii+1);
  100.   elemA = d2  ELEM COMP pp1a pp2a;
  101.   elemB = d2a ELEM COMP pp1b pp2b;
  102.   fiss0 = RACC 0.0001 (INVE elemA) elemB;
  103.   xx1 = COOR 2 pp1a;
  104.   xx2 = COOR 2 pp2a;
  105.   yy1 = fdepth * ( (1. - ((xx1 / flength)**2.))**.5 );
  106.   yy2 = fdepth * ( (1. - ((xx2 / flength)**2.))**.5 );
  107.   aai = ( yy1 + yy2 ) / 2.;
  108.  
  109.   mod0 = MODE fiss0 MECANIQUE ELASTIQUE LISP;
  110.   mat0 = MATE mod0 YOUN xYoung NU .3  EPAI ww FISS aai
  111.          VX 0.  VY 0.  VZ -1.;
  112. *
  113.   SI (ii EGA 1);
  114.      mod2  = mod0;
  115.      mat2  = mat0;
  116.      fiss1 = fiss0;
  117.    SINON;
  118.      mod2  = mod2 ET mod0;
  119.      mat2  = mat2 ET mat0;
  120.      fiss1 = fiss1 ET fiss0;
  121.   FINSI;
  122. *
  123. FIN matFiss;
  124. *
  125.  
  126. *  Assemblage des 2 maillages
  127. surf1 = su9 ET fiss1;
  128. SI ( NEG GRAPH 'N' );trac surf1 (0. 0. 1.E6); FINSI;
  129.  
  130. *_____________________________________________________________________
  131. *
  132. *--- conditions aux limites
  133. *
  134. enc1 = SYMT DEPL ROTA (0. 0. 0.) (1. 0. 0.) (0. 0. 1.) surf1 0.0001;
  135. enc2 = BLOQ d8a UX UZ RX ;
  136. encT = enc1 ET enc2;
  137. *_____________________________________________________________________
  138. *
  139. mod = mod1 et mod2 ;
  140. mat = mat1 et mat2 ;
  141. *__________________________________________________________
  142. *
  143. *--- chargement : contrainte de membrane en traction
  144. *
  145. SIG = 1.;
  146. effortN = larg1 * ww * SIG;
  147. ff = ((nb + 1.)/ nb) * effortN;
  148. f1 = ff 0. 0.;
  149. f2 = (ff / (-2. * (nb + 1.))) 0. 0.;
  150.  
  151. char1 = (FORC f1 d8) et (FORC f2 P5) et (FORC f2 P9);
  152.  
  153. *__________________________________________________________
  154. *
  155. * --- calcul rigidite
  156. *
  157. rig0 = (RELA 1. UX p06 - 1. UX p6) ET (RELA 1. UY p06 - 1. UY p6)
  158.         ET (RELA 1. UZ p06 - 1. UZ p6) ET (RELA 1. RX p06 - 1. RX p6)
  159.         ET (RELA 1. RY p06 - 1. RY p6) ET (RELA 1. RZ p06 - 1. RZ p6);
  160. *
  161. rig1 = (RIGI mod1 mat1) ET rig0;
  162. rig2 = (RIGI mod2 mat2);
  163. *__________________________________________
  164. *
  165. * --- Resolution:
  166. *
  167. * Force:
  168. rig = rig1 et rig2 et encT;
  169. re1  = RESO rig char1;
  170.  
  171. *_____________________________________________________________________
  172. *
  173. * ---  FACTEUR D'INTENSITE DE CONTRAINTE
  174. *  détermination du facteur F = K / (SIG * (pi * a / Q)^0.5) 
  175. *
  176.  
  177. * ---  valeurs calculées
  178. sigf = SIGM mod2 mat2 re1 ;
  179.  
  180. pKcal = prog (nbel fiss1) * 0.;py1 = prog (nbel fiss1) * 0.;
  181. i=0; repe bbb (nbel fiss1) ; i = i + 1;
  182.   elemi = elem fiss1 i;
  183.   sigi = redu sigf elemi;
  184.   vki = extr sigi 'KI' 1 1 2;remp pKcal i vki;
  185.   yi = (coor 2 (bary elemi));remp py1 i yi;
  186. fin bbb;
  187.  
  188. * ---  valeurs de Raju-Newman
  189. *
  190.  
  191. a = fdepth ; t = ww ; c = flength ;
  192. b = 4 * c;ac = a/c; at = a/t;
  193. Q  = (1.0) + (1.464*((ac)**1.65));
  194. K0 = SIG * ((pi * a / Q) ** 0.5);
  195. pFrn = prog 1.178 1.171 1.151 1.117 1.07 1.01
  196.             0.938 0.851 0.749 0.64 0.587;
  197. pKrn = pFrn * K0;
  198. pphin2 = prog 1. pas -0.1 0.;
  199. py2= (cos (pphin2 * 90.)) * c;
  200.  
  201. * ---  COMPARAISON
  202. *  écart maxi en excluant les 3 elts près de la surface libre
  203.  
  204. nn=nbel fiss1; pecart = prog (nn - 3) * 0.;
  205.  
  206. i=0 ; repe bbb (nn - 3); i = i + 1;
  207. vKcal = extr pKcal i;
  208. vy1 = extr py1 i;
  209. vKrn = ipol vy1 py2 pKrn;
  210. vecart= (vKcal - vKrn) / vKrn; vecart = abs vecart;
  211. remp pecart i vecart;
  212. fin bbb;
  213.  
  214. ERR = maxi pecart;
  215. SI (ERR <  3.4);ERRE 0;SINO;ERRE 5;FINSI;
  216.  
  217. SI ( NEG GRAPH 'N' );
  218. t=tabl;
  219. t.1='MARQ CARR ';
  220. t.2='MARQ ETOI ';
  221. t.'TITRE'=tabl;
  222. t.'TITRE'. 1 = mot 'RAJU-NEWMANN';
  223. t.'TITRE'. 2 = mot 'LINESPRING';
  224. evcal = evol manu 'y (mm)' py1 'K (MPam^0.5)' pKcal;
  225. evrn = evol manu 'y (mm)' py2 'K (MPam^0.5)' pKrn;
  227. 'K = f (Y) - plaque en tract. - defaut semi-ellip. a/c=a/t=0.2';
  228. dess (evrn et evcal) lege t logo xbor 0. 20. ybor 0. 5.;
  229. FINSI;
  230. fin;
  231.  
  232.  
  233.  
  234.  
  235.  
  236.  
  237.  
  238.  

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