Télécharger la plaquette Cast3M


fr en

Télécharger lispnl.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  lispnl.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. OPTION ECHO 0;
  5. GRAPH = 'N';
  6. SAUT PAGE;
  7. *                                                                               
  8. *   PLAQUE FISSUREE SOLLICITEE EN TRACTION PURE
  9. *        PLASTICITE PARFAITE                                       
  10. *                                            
  11. opti dime 3 elem qua4;
  12. a =  25.;b = 51.5;c=b - a;h=525 / 2.;r = a / b;
  13. E = 207000.;Eprim=E  / (1. - (0.3 ** 2));vsigf=471.;
  14. p1 =   0.001  0.  0. ;p2 =   0.001 50.  0. ;
  15. p3 =   0.  0. 0.;p4 =   0. 50. 0.;
  16. d1 = p1 d 1 p2;d2 = p3 d 1 p4 ;
  17. cad   = d1  trans 5 ( h 0. 0.) ;
  18. lchar = cad cote 3 ;pchar = lchar poin init;
  19. cag   = d2  trans 5 ((h * -1) 0. 0.) ;
  20. lbloc = cag cote 3 ;
  21. tole  = cad et cag ;
  22. *
  23. * Approximation de la courbe de traction parfaite
  24. * une loi de Ramberg-Osgood avec alpha=1 et n=20
  25. *  Comparaison des courbes
  26. n=20;sig0=410.;
  27. SI ( NEG GRAPH 'N' );
  28. Psig=prog 0. pas 10. sig0 pas 1. 510.;
  29. Peps=prog (dime Psig) * 0.;
  30. evcal=evol manu (prog 0. (vsigf / E) 0.3) (prog 0. vsigf vsigf);
  31. i=0;repeter bbb (dime Psig);i=i + 1;
  32. sig=extr Psig i;
  33. eps=(sig0 / E)  * ((sig / sig0) + ((sig / sig0) ** n));
  34. remp Peps i eps;
  35. fin bbb;
  36. evshih=evol manu 'eps' Peps 'sig (MPa)' Psig;
  37. t=tabl;
  38. t.1='MARQ CARR ';
  39. t.2='MARQ ETOI ';
  40. t.'TITRE'=tabl;
  41. t.'TITRE'. 1 = mot 'PARF';
  42. t.'TITRE'. 2 = mot 'N=20';
  43. titr 'Courbes de traction plastique parfaite et ecrouissage n=20';
  44. dess (evcal et evshih) t logo;
  45. FINSI;
  46. *
  47. *
  48. mod1 = mode tole MECANIQUE ELASTIQUE PLASTIQUE PARFAIT coq4;
  49. mat1 = mate mod1 youn E nu 0.3 epai b sigy sig0;
  50. fis =    racc ( d1) ( d2) 0.001 ;
  51. mod2 = mode fis MECANIQUE ELASTIQUE PLASTIQUE LINESPRING LISP;
  52. mat2 = mate mod2  youn E  nu 0.3  epai b  fiss a
  53.        vx 0. vy 0. vz 1. SIGF vsigf;
  54. mo = mod1 et mod2 ;
  55. ma = mat1 et mat2 ;
  56. total = tole et fis ;
  57. *
  58. blo1 = bloq lbloc UX UZ RX RY RZ ;
  59. blo2 =bloq P1 UY;
  60. blo4 = bloq lchar UX;
  61. charx = depi blo4 1.;
  62. blo = blo1 et blo2 et blo4;
  63. rig = (RIGI mo ma) et blo ;
  64. *
  65. EVT=EVOL MANU ABS (PROG 0. 12.) ORD (PROG 0. 1.2);
  66. PGT=PROG 1. PAS 1 8. pas 0.2 9.8;
  67. FDT=CHAR charx EVT 'MECA';
  68. *
  69. TENTR = TABLE;
  70. TENTR.CHARGEMENT  = FDT;
  71. TENTR.BLOCAGES_MECANIQUES=blo;
  72. TENTR.MODELE=mo;
  73. TENTR.CARACTERISTIQUES=ma;
  74. TENTR.TEMPS_CALCULES=PGT;  
  75. PASAPAS TENTR;
  76. *
  77. tdep=TENTR.'DEPLACEMENTS';
  78. tsig=TENTR.'CONTRAINTES';
  79. tvar=TENTR.'VARIABLES_INTERNES';
  80. *
  81. F=1.12 - (0.23 * r) + (10.6 * (r ** 2)) - (21.7 * (r ** 3))
  82. + (30.4 * (r ** 4));coef=F * ((pi * a) ** 0.5);
  83. H1=0.0391;
  84. nu1=((1. + ((a / c) ** 2)) ** 0.5) - (a / c);
  85. P0=(2. / (3. ** 0.5)) * 1.26 * nu1 * c * sig0;
  86. *
  87. NPAS=dime PGT;MPAS=NPAS + 1; 
  88. PU=PROG MPAS * 0.;
  89. PJEcal=PROG MPAS * 0.;PJEshih=PROG MPAS * 0.;
  90. PJPcal=PROG MPAS * 0.;PJPshih=PROG MPAS * 0.;
  91. PJcal=PROG MPAS * 0.;PJshih=PROG MPAS * 0.;
  92. *
  93. repeter blo NPAS;i=&blo;  
  94.  VU=extr (tdep.i) 'UX' pchar;
  95.  sig=tsig.i;var=tvar.i; 
  96.  RSIG=REDU SIG fis;
  97.  P=extr RSIG 'NZZ' 1 1 1;
  98.  RVAR=REDU VAR fis;
  99.  VK1cal=EXTR RSIG KI 1 1 1;  VJEcal=(VK1cal ** 2) / Eprim;
  100.  VJPcal=EXTR RVAR VHJP 1 1 1;
  101.   VJcal=VJEcal + VJPcal;
  102.  VK1shih=coef * (P / b);VJEshih=(VK1shih ** 2) / Eprim;
  103.  VJPshih=(sig0 ** 2) / E * c * a / b * h1 * ((P / P0) ** n);
  104.   VJshih=VJEshih + VJPshih;
  105.  REMP PU (i + 1) VU;
  106.  REMP PJEcal (i + 1) VJEcal;REMP PJEshih (i + 1) VJEshih;
  107.  REMP PJPcal (i + 1) VJPcal;REMP PJPshih (i + 1) VJPshih;
  108.  REMP PJcal (i + 1) VJcal;REMP PJshih (i + 1) VJshih;
  109. SI (i ega NPAS);
  110. ERRelas=abs ((VJEcal - VJEshih) / VJEshih);
  111. ERRplas=abs ((VJPcal - VJPshih) / VJPshih);
  112. SI (ERRelas <  3.5E-2);ERRE 0;SINO;ERRE 5;FINSI;
  113. SI (ERRplas <  8.E-2);ERRE 0;SINO;ERRE 5;FINSI;
  114. FINSI; 
  115. fin blo;
  116. *
  117. SI ( NEG GRAPH 'N' );
  118. EVEcal=EVOL VERT MANU 'U (mm)' PU 'J (KJ/M^2)' PJEcal;
  119. EVEshih=EVOL ROUG MANU 'U (mm)' PU 'J (KJ/M^2)' PJEshih;
  120. EVPcal=EVOL VERT MANU 'U (mm)' PU 'J (KJ/M^2)' PJPcal;
  121. EVPshih=EVOL ROUG MANU 'U (mm)' PU 'J (KJ/M^2)' PJPshih;
  122. EVcal=EVOL VERT MANU 'U (mm)' PU 'J (KJ/M^2)' PJcal;
  123. EVshih=EVOL ROUG MANU 'U (mm)' PU 'J (KJ/M^2)' PJshih;
  124. dess (evEcal et evEshih) t logo;
  125. dess (evPcal et evPshih) t logo;
  126. dess (evcal et evshih) t logo;
  127. FINSI;
  128. FIN;
  129.  
  130.  
  131.  
  132.  
  133.  
  134.  
  135.  
  136.  

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