Télécharger ricrag_2d.dgibi

Retour à la liste

Numérotation des lignes :

  1. *******************************************************
  2. ************************************************************************
  3. ************************************************************************
  4. * Cas test de l'implantation numérique du modele
  5. * RICRAG 2D LOCAL/NON LOCAL
  6. *******************************************************
  7. * Développé par :
  8. * Benjamin Richard
  9. * Contact :
  10. * Benjamin.Richard@lmt.ens-cachan.fr
  11. *******************************************************
  12. * Les cas de charges sont entrés :
  13. * - 1 : Traction monotone
  14. * - 2 : Compression monotone
  15. * - 3 : Traction cyclique
  16. * - 4 : Compression cyclique
  17. * - 5 : Traction/compression cyclique
  18. *******************************************************
  19. * Choix du cas de charge
  20. ncas = 1;
  21.  
  22. * Test du fichier compatible avec le non local
  23. * nloc0 = 0; Cas local
  24. * nloc0 = 1; Cas non local
  25. graph= mot 'N';
  26. nloc0 = 0;
  27. *******************************************************
  28.  
  29. *
  30. *-------------- Options de calcul ---------------------
  31. *
  32. OPTION DIME 2 ELEM QUA4 MODE PLAN CONT;
  33.  
  34. *
  35. *-------------- Definition de la geometrie ------------
  36. *
  37. P1 = 0. 0.;
  38. P2 = 1. 0.;
  39. P3 = 1. 1.;
  40. P4 = 0. 1.;
  41.  
  42. L1 = P1 DROIT 1 P2 ;
  43. L2 = P2 DROIT 1 P3 ;
  44. L3 = P3 DROIT 1 P4 ;
  45. L4 = P4 DROIT 1 P1 ;
  46.  
  47. LTOT = L1 ET L2 ET L3 ET L4 ;
  48.  
  49. S1 = SURF LTOT PLANE;
  50.  
  51. *
  52. *------- Définition des conditions aux limites --------
  53. *----------- et des déplacements imposés --------------
  54. *
  55. CL = BLOQ P1 UX UY;
  56. CLL = BLOQ P2 UY;
  57. CL1 = BLOQ L3 UY;
  58.  
  59. D1 = DEPI CL1 1;
  60.  
  61. *
  62. *----- Définition du modèle ---------------------------
  63. *
  64. MOD1 = MODE S1 MECANIQUE ELASTIQUE ISOTROPE
  65. ENDOMMAGEMENT RICRAG;
  66.  
  67. *
  68. *----- Paramètres matériaux ---------------------------
  69. *
  70.  
  71. * Module d'Young
  72. youngn = 36000E+6;
  73.  
  74. * Coefficient de Poisson
  75. nun = 0.2;
  76.  
  77. * Résistance en traction
  78. ftn = 3.6e6;
  79.  
  80. * Fragilité en traction
  81. aldin = 1.0e-2;
  82.  
  83. * Fragilité en compression
  84. alinn = 5.0e-4;
  85.  
  86. * Module d'écrouissage 1
  87. gam1n = 7.0e9;
  88.  
  89. * Moduke d'écrouissage 2
  90. a1n = 7.0e-7;
  91.  
  92. mat1 =MATE mod1 YOUN youngn NU nun
  93. FT ftn ALIN alinn
  94. GAM1 gam1n A1 a1n
  95. ALDI aldin;
  96.  
  97. *
  98. *----- Définition des cas de charge -------------------
  99. *
  100.  
  101. SI (EGA ncas 1) ;
  102. LI1 = PROG 0. 1.;
  103. LI2 = PROG 0. 4.0e-4;
  104. LIS1 = PROG 0. PAS 0.02 1.;
  105. FINSI;
  106.  
  107. SI (EGA ncas 2) ;
  108. LI1 = PROG 0. 1.;
  109. LI2 = PROG 0. -8.0e-3;
  110. LIS1 = PROG 0. PAS 0.02 1.;
  111. FINSI;
  112.  
  113. SI (EGA ncas 3) ;
  114. LI1 = PROG 0. 1. 2. 3. 4. 5.;
  115. LI2 = PROG 0. 1.5e-4 9.5e-6 2.0E-4 3.0E-5 2.5E-4;
  116. LIS1 = PROG 0. PAS 0.02 5.;
  117. FINSI;
  118.  
  119. SI (EGA ncas 4) ;
  120. LI1 = PROG 0. 1. 2. 3. 4. 5.;
  121. LI2 = PROG 0. -3.0e-3 -1.5e-4 -5.0E-3 -3.5E-4 -8.0E-3;
  122. LIS1 = PROG 0. PAS 0.02 5.;
  123. FINSI;
  124.  
  125. SI (EGA ncas 5) ;
  126. LI1 = PROG 0. 1. 2. 3. 4. 5. 6.;
  127. LI2 = PROG 0. 1.3e-4 -3.0e-3 -1.5E-4 -5.0E-3 -3.5E-4 -8.0E-3;
  128. LIS1 = PROG 0. PAS 0.02 6;
  129. FINSI;
  130.  
  131. EV = EVOL MANU LI1 LI2 ;
  132. CHA1 = CHAR 'DIMP' D1 EV ;
  133.  
  134. *
  135. *----------- Calcul par l'operateur PASAPAS ------------
  136. *
  137. SI (EGA nloc0 1);
  138. LC = 0.000001;
  139. CO1 = CONNEC mod1 LC NORMAL;
  140. FINSI;
  141.  
  142. TAB1 = TABLE ;
  143. TAB1.'BLOCAGES_MECANIQUES' = CL ET CLL ET CL1;
  144. TAB1.'MODELE' = MOD1;
  145. TAB1.'MOVA' = 'D ';
  146. TAB1.'CHARGEMENT' = CHA1;
  147. TAB1.'CARACTERISTIQUES' = MAT1;
  148. TAB1.'TEMPS_CALCULES' = LIS1;
  149.  
  150. SI (EGA nloc0 1);
  151. TAB1.NON_LOCAL='MOT' 'MOYE' ;
  152. TAB1.LVIAM=MOTS EPTI;
  153. TAB1.LCAR=LC;
  154. TAB1.CONN = CO1;
  155. FINSI;
  156.  
  157. PASAPAS TAB1 ;
  158.  
  159. *
  160. *----------- Courbe effort-deplacement -----------------
  161. *
  162. ev2=@global tab1 CL1 EV fy;
  163. si ( ega graph 'O');
  164.  
  165. DESS EV2;
  166. finsi;
  167. ee = extr ev2 ordo 1;
  168. aa = extr ee ( dime ee);
  169. list aa;
  170. err = abs (aa - 2.37395E+06) /2.37395E+06 ;
  171. message ' erreur relative ' err;
  172. si (err > 1.e-3);
  173. erreur (5);
  174. finsi;
  175. fin;
  176.  
  177.  
  178.  
  179.  
  180.  
  181.  
  182.  
  183.  
  184.  
  185.  
  186.  
  187.  
  188.  

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