Télécharger concyc.dgibi

Retour à la liste

Numérotation des lignes :

  1. *******************************************************
  2. ************************************************************************
  3. ************************************************************************
  4. * Cas test de l'implantation numérique du modele
  5. * RICRAG 3D LOCAL/NON LOCAL
  6. *******************************************************
  7. * Développé par :
  8. * Maxime VASSAUX
  9. * Benjamin RICHARD
  10. *
  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
  18. *******************************************************
  19. * Choix du cas de charge
  20. ncas = 1;
  21. graph=mot 'N';
  22. * Test du fichier compatible avec le non local
  23. * nloc0 = 0; Cas local
  24. * nloc0 = 1; Cas non local
  25.  
  26. nloc0 = 0;
  27. *******************************************************
  28.  
  29. *
  30. *-------------- Options de calcul ---------------------
  31. *
  32. OPTION DIME 3 ELEM CUB8;
  33.  
  34. *
  35. *-------------- Definition de la geometrie ------------
  36. *
  37. dim1 = 1.0;
  38. dim2 = 1.0;
  39. dim3 = 1.0;
  40. P1 = 0. 0. 0.;
  41. P2 = dim1 0. 0.;
  42. P3 = dim1 dim2 0.;
  43. P4 = 0. dim2 0.;
  44. P5 = 0. 0. dim3;
  45. P6 = dim1 0. dim3;
  46. P7 = dim1 dim2 dim3;
  47. P8 = 0. dim2 dim3;
  48.  
  49. nelem1 = 1;
  50. nelem2 = 1;
  51. nelem3 = 1;
  52.  
  53. L1 = P1 DROIT nelem1 P2 ;
  54. L2 = P2 DROIT nelem2 P3 ;
  55. L3 = P3 DROIT nelem1 P4 ;
  56. L4 = P4 DROIT nelem2 P1 ;
  57.  
  58. L5 = P5 DROIT nelem1 P6 ;
  59. L6 = P6 DROIT nelem2 P7 ;
  60. L7 = P7 DROIT nelem1 P8 ;
  61. L8 = P8 DROIT nelem2 P5 ;
  62.  
  63. L9 = P1 DROIT nelem3 P5 ;
  64. L10 = P2 DROIT nelem3 P6 ;
  65. L11 = P3 DROIT nelem3 P7 ;
  66. L12 = P4 DROIT nelem3 P8 ;
  67.  
  68. LTOT1 = L1 ET L2 ET L3 ET L4 ;
  69. SURF1 = SURF LTOT1 PLANE;
  70.  
  71. LTOT2 = L5 ET L6 ET L7 ET L8 ;
  72. SURF2 = SURF LTOT2 PLANE;
  73.  
  74. LTOT3 = L1 ET L5 ET L10 ET L9 ;
  75. SURF3 = SURF LTOT3 PLANE;
  76.  
  77. LTOT4 = L2 ET L11 ET L6 ET L10 ;
  78. SURF4 = SURF LTOT4 PLANE;
  79.  
  80. LTOT5 = L3 ET L12 ET L7 ET L11 ;
  81. SURF5 = SURF LTOT5 PLANE;
  82.  
  83. LTOT6 = L4 ET L9 ET L8 ET L12 ;
  84. SURF6 = SURF LTOT6 PLANE coul 'ROUG';
  85.  
  86. STOT = SURF1 et SURF2 et SURF3 et SURF4 et SURF5 et SURF6;
  87.  
  88. VOLTO = surf1 volu nelem3 tran (P5);
  89.  
  90. VOLTOT = VOLTO;
  91.  
  92. elim stot voltot 0.0001;
  93.  
  94. *
  95. *------- Définition des conditions aux limites --------
  96. *----------- et des déplacements imposés --------------
  97. *
  98. CL1 = BLOQ 'UZ' SURF1;
  99. CL3 = BLOQ 'UY' SURF3;
  100. CL6 = BLOQ 'UX' SURF6;
  101.  
  102. BCL = CL1 ET CL3 ET CL6;
  103. CL2 = BLOQ SURF2 UZ;
  104.  
  105. D1 = DEPI CL2 1;
  106.  
  107. *
  108. *----- Définition du modèle ---------------------------
  109. *
  110. MOD1 = MODE VOLTOT MECANIQUE ELASTIQUE ISOTROPE
  111. ENDOMMAGEMENT CONCYC;
  112.  
  113. *
  114. *----- Paramètres matériaux ---------------------------
  115. *
  116. * Module de Young
  117. youngn = 30.0E9;
  118.  
  119. * Coefficient de Poisson
  120. nun = 0.2;
  121.  
  122. * Choix du modèle d'endommagement pour 'CONCYC' : Ricrag (1) - consolidation modifiee (2) - critere modifie (3)
  123. nendn = 2;
  124.  
  125. * Contrainte limite elastique
  126. sigtn = 3.00E6;
  127.  
  128. * Parametre controlant l'energie de fissuration
  129. atran = 0.004;
  130.  
  131. * RicRag classique *
  132. SI (EGA nendn 1);
  133. btran = 0.0;
  134. FINSI;
  135.  
  136. * RicRag avec consolidation modifiee *
  137. SI (EGA nendn 2);
  138. btran = 4.5;
  139. FINSI;
  140.  
  141. * RicRag avec critere modifie *
  142. SI (EGA nendn 3);
  143. btran = 0.31;
  144. FINSI;
  145.  
  146. * "Vitesse" de refermeture de fissure
  147. qpn = 6.5;
  148.  
  149. * Coefficient de frottement des fissures
  150. cfn = 2.89;
  151.  
  152. MAT1 = MATE MOD1 'YOUN' youngn 'NU ' nun
  153. 'NEND' nendn 'SIGT' sigtn
  154. 'ATRA' atran 'BTRA' btran
  155. 'QP ' qpn 'CF ' cfn ;
  156.  
  157. *
  158. *----- Définition des cas de charge -------------------
  159. *
  160.  
  161. SI (EGA ncas 1) ;
  162. LI1 = PROG 0. 1.;
  163. LI2 = PROG 0. 4.0e-4;
  164. LIS1 = PROG 0. PAS 0.0025 1.;
  165. FINSI;
  166.  
  167. SI (EGA ncas 2) ;
  168. LI1 = PROG 0. 1.;
  169. LI2 = PROG 0. -8.0e-3;
  170. LIS1 = PROG 0. PAS 0.0025 1.;
  171. FINSI;
  172.  
  173. SI (EGA ncas 3) ;
  174. LI1 = PROG 0. 1. 2. 3. 4. 5.;
  175. LI2 = PROG 0. 1.5e-4 -2.0e-4 2.0E-4 -2.0e-4 2.5E-4;
  176. LIS1 = PROG 0. PAS 0.0025 5.;
  177. FINSI;
  178.  
  179. SI (EGA ncas 4) ;
  180. LI1 = PROG 0. 1. 2. 3. 4. 5.;
  181. LI2 = PROG 0. -3.0e-3 -1.5e-4 -5.0E-3 -3.5E-4 -8.0E-3;
  182. LIS1 = PROG 0. PAS 0.0025 5.;
  183. FINSI;
  184.  
  185. SI (EGA ncas 5) ;
  186. LI1 = PROG 0. 1. 2.;
  187. LI2 = PROG 0. 2.0e-4 -4.0e-3;
  188. LIS1 = PROG 0. PAS 0.0025 2;
  189. FINSI;
  190.  
  191. EV = EVOL MANU LI1 LI2 ;
  192. CHA1 = CHAR 'DIMP' D1 EV ;
  193.  
  194. *
  195. *----------- Calcul par l'operateur PASAPAS ------------
  196. *
  197. LC = 1.0e-10;
  198. CO1 = CONNEC mod1 LC NORMAL;
  199.  
  200. TAB1 = TABLE ;
  201. TAB1.'BLOCAGES_MECANIQUES' = BCL ET CL2;
  202. TAB1.'MODELE' = MOD1;
  203. TAB1.'MOVA' = 'D ';
  204. TAB1.'CHARGEMENT' = CHA1;
  205. TAB1.'CARACTERISTIQUES' = MAT1;
  206. TAB1.'TEMPS_CALCULES' = LIS1;
  207. TAB1.'MAXITERATION' = 10;
  208.  
  209. SI (EGA nloc0 1);
  210. TAB1.NON_LOCAL='MOT' 'MOYE' ;
  211. TAB1.LVIAM=MOTS EPTI;
  212. TAB1.LCAR=LC;
  213. TAB1.CONN = CO1;
  214. FINSI;
  215.  
  216. PASAPAS TAB1 ;
  217.  
  218. *
  219. *----------- Courbe effort-deplacement -----------------
  220. *
  221. ev2=@global tab1 CL2 EV fz;
  222. ee = extr ev2 ordo 1;
  223. aa = extr ee ( dime ee);
  224. list aa;
  225. err = abs (((aa) - (5.14929E+05)) / (5.14929E+05)) ;
  226. message ' erreur relative ' err;
  227. si (err > 1.e-3);
  228. erreur (5);
  229. finsi;
  230.  
  231. si ( ega graph 'O');
  232. DESS EV2;
  233. finsi;
  234. list ev2;
  235. fin;
  236.  
  237.  
  238.  
  239.  
  240.  
  241.  
  242.  
  243.  
  244.  
  245.  
  246.  

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