Télécharger ricbet_3d.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : ricbet_3d.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *******************************************************
  5. * Cas test de l'implantation numerique du modele
  6. * RICBET LOCAL/NON LOCAL 3D
  7. *******************************************************
  8. * Développé par :
  9. * Benjamin Richard
  10. * Contact :
  11. * Benjamin.Richard@lmt.ens-cachan.fr
  12. * Institution :
  13. * CEA\DEN\DANS\DM2S\SEMT\EMSI
  14. * LMT\ENS Cachan
  15. *******************************************************
  16. * Les cas de charges sont entrés :
  17. * - 1 : Traction monotone
  18. * - 2 : Compression monotone
  19. * - 3 : Traction cyclique
  20. * - 4 : Compression cyclique
  21. * - 5 : Traction/compression cyclique
  22. * - 6 : Compression/traction
  23. * - 7 : Traction/Compression/Traction 1
  24. * - 8 : Traction/Compression/Traction 2
  25. *******************************************************
  26. graph = 'N';
  27.  
  28. * Choix du cas de charge-
  29. ncas = 1;
  30.  
  31. * Test du fichier compatible avec le non local
  32. * nloc0 = 0; Cas local
  33. * nloc0 = 1; Cas non local
  34.  
  35. nloc0 = 0;
  36.  
  37. * Calcul élastiques (pour référence avant fissuration)
  38. * nelas = 1; cas non lineaire
  39. * nelas = 0; cas lineaire
  40.  
  41. nelas = 1;
  42. *******************************************************
  43.  
  44. *
  45. *-------------- Options de calcul ---------------------
  46. *
  47. OPTION DIME 3 ELEM CUB8;
  48.  
  49. *
  50. *-------------- Definition de la geometrie ------------
  51. *
  52. P1 = 0. 0. 0.;
  53. P2 = 1. 0. 0.;
  54. P3 = 1. 1. 0.;
  55. P4 = 0. 1. 0.;
  56. P5 = 0. 0. 1. ;
  57.  
  58. L1 = P1 DROIT 1 P2 ;
  59. L2 = P2 DROIT 1 P3 ;
  60. L3 = P3 DROIT 1 P4 ;
  61. L4 = P4 DROIT 1 P1 ;
  62.  
  63. LTOT = L1 ET L2 ET L3 ET L4 ;
  64. SURF1 = SURF LTOT PLANE;
  65. VOLTO =VOLU SURF1 1 TRANS P5;
  66. SURF2= FACE VOLTO 2;
  67.  
  68. VOLTOT = VOLTO;
  69.  
  70. *
  71. *------- Définition des conditions aux limites --------
  72. *----------- et des déplacements imposés --------------
  73. *
  74. CL = BLOQ SURF1 UZ;
  75. CLL = BLOQ P1 'DEPL';
  76. CL1 = BLOQ SURF2 Uz;
  77. CL2 = bloq (noeud 3) UY;
  78. cl=cl et cl2;
  79.  
  80. D1 = DEPI CL1 1;
  81.  
  82. *
  83. *----- Définition du modèle ---------------------------
  84. *
  85. SI (EGA nelas 1);
  86. SI (EGA nloc0 1);
  87. MOD1=MODE VOLTOT MECANIQUE ELASTIQUE ISOTROPE ENDOMMAGEMENT RICBET
  88. 'NON_LOCAL' 'MOYE' 'V_MOYENNE' (MOTS 'EPTI') ;
  89. SINON;
  90. MOD1=MODE VOLTOT MECANIQUE ELASTIQUE ISOTROPE ENDOMMAGEMENT RICBET;
  91. FINSI ;
  92. SINON;
  93. MOD1=MODE VOLTOT MECANIQUE ELASTIQUE ISOTROPE;
  94. FINSI;
  95. EBETON = 34100E6;
  96. NUBETON = 0.16;
  97. LC = 0.020;
  98. FTN = 3.80E+6;
  99.  
  100. * Young modulus
  101. youngn = EBETON;
  102.  
  103. * Poisson ratio
  104. nun = NUBETON;
  105.  
  106. * Brittleness in tension
  107. aldin = 6.0e-3;
  108.  
  109. * Hardening modulus 1
  110. gam1n = 5.0e9;
  111.  
  112. * Hardening modulus 2
  113. a1n = 8.0e-6;
  114.  
  115. * Closure stress
  116. srefn = -3e6;
  117.  
  118. * Compression strength (initial)
  119. fcn = 6e6;
  120.  
  121. * Asymptotic stress in compression
  122. sigun = -4e6;
  123.  
  124. * Yield surface 1
  125. afn = 0.7;
  126.  
  127. * Yield surface 2
  128. bfn = 0.3;
  129.  
  130. * Plastic potential 1
  131. agn = 0.6;
  132.  
  133. * Plastic potential 2
  134. bgn = 0.45;
  135.  
  136. * Yield hardening 1
  137. acn = 3.2e10;
  138.  
  139. * Yield hardening 1
  140. bcn = 700.0;
  141.  
  142. SI (EGA nelas 1);
  143. mat1 =MATE mod1 YOUN youngn NU nun
  144. FT ftn GAM1 gam1n
  145. A1 a1n ALDI aldin
  146. SREF srefn AF afn BF bfn
  147. AG agn BG bgn AC acn BC bcn
  148. SIGU sigun FC fcn;
  149. SINON;
  150. mat1 =MATE mod1 YOUN youngn NU nun;
  151. FINSI;
  152.  
  153.  
  154.  
  155.  
  156. *
  157. *----- Définition des cas de charge -------------------
  158. *
  159.  
  160. SI (EGA ncas 1) ;
  161. LI1 = PROG 0. 1. ;
  162. LI2 = PROG 0. 1.0e-3;
  163. LIS1 = PROG 0. PAS 0.01 1.;
  164. FINSI;
  165.  
  166. SI (EGA ncas 2) ;
  167. LI1 = PROG 0. 1.;
  168. LI2 = PROG 0. -1;
  169. LIS1 = PROG 0. PAS 1.E-4 6.0E-3;
  170. FINSI;
  171.  
  172. SI (EGA ncas 3) ;
  173. LI1 = PROG 0. 1. 2. 3. 4. 5.;
  174. LI2 = PROG 0. 1.5e-4 9.5e-6 2.0E-4 3.0E-5 2.5E-4;
  175. LIS1 = PROG 0. PAS 0.02 5.;
  176. FINSI;
  177.  
  178. SI (EGA ncas 4) ;
  179. LI1 = PROG 0. 1. 2. 3. 4. 5.;
  180. LI2 = PROG 0. -3.0e-3 -1.5e-4 -5.0E-3 -3.5E-4 -10.0E-3;
  181. LIS1 = PROG 0. PAS 0.02 5.;
  182. FINSI;
  183.  
  184. SI (EGA ncas 5) ;
  185. LI1 = PROG 0. 1. 2. 3. 4. 5. 6. 7.;
  186. LI2 = PROG 0. 1.3e-4 -3.0e-3 -1.5E-4 -5.0E-3 -3.5E-4 -8.0E-3 1.0E-3;
  187. LIS1 = PROG 0. PAS 0.005 7;
  188. FINSI;
  189.  
  190. SI (EGA ncas 6) ;
  191. LI1 = PROG 0. 1. 2.;
  192. LI2 = PROG 0. -2.0e-3 5.0e-4;
  193. LIS1 = PROG 0. PAS 0.005 2;
  194. FINSI;
  195.  
  196. SI (EGA ncas 7) ;
  197. LI1 = PROG 0. pas 1. 11.;
  198. LI2 = PROG 0. 5.0e-4 -5.0e-4 6.0e-4 -5.0e-4
  199. 7.0e-4 -5.0e-4 8.0e-4 -5.0e-4
  200. 9.0e-4 -5.0e-4 10.0e-4;
  201. LIS1 = PROG 0. PAS 0.01 11;
  202. FINSI;
  203.  
  204. SI (EGA ncas 8) ;
  205. LI1 = PROG 0. 1. 2. 3. 4. 5.;
  206. LI2 = PROG 0. 3.e-4 -1.e-5 5.0E-4 -5E-3 -4.0E-3;
  207. LIS1 = PROG 0. PAS 0.02 2;
  208. FINSI;
  209.  
  210.  
  211. EV = EVOL MANU LI1 LI2 ;
  212. CHA1 = CHAR 'DIMP' D1 EV ;
  213.  
  214. *
  215. *----------- Calcul par l'operateur PASAPAS ------------
  216. *
  217. LC = 1.0e-10;
  218. CO1 = CONNEC mod1 LC NORMAL;
  219.  
  220. TAB1 = TABLE ;
  221. TAB1.'BLOCAGES_MECANIQUES' = CL ET CLL ET CL1;
  222. TAB1.'MODELE' = MOD1;
  223. TAB1.'MOVA' = 'D ';
  224. TAB1.'CHARGEMENT' = CHA1;
  225. TAB1.'TEMPS_CALCULES' = LIS1;
  226.  
  227. SI (EGA nloc0 1);
  228. TAB1.CONN = CO1;
  229. MAT1 = MAT1 ET ( MATE MOD1 'LCAR' LC) ;
  230. FINSI;
  231.  
  232. TAB1.'CARACTERISTIQUES' = MAT1;
  233. PASAPAS TAB1 ;
  234.  
  235. dd=dime tab1.deplacements;
  236.  
  237. *
  238. *----------- Courbe effort-deplacement -----------------
  239. *
  240. si (ega graph 'O');
  241. ev2=@global tab1 CL1 EV fz;
  242. @excel1 ev2 'loi_beton-RL5.txt';
  243. DESS EV2;
  244. finsi;
  245.  
  246. ee=tab1.variables_internes.(dd-1);
  247. aa=extr ee 'D ' 1 1 1;
  248. list aa;
  249. message 'erreur relative';
  250.  
  251. err = (aa - 0.99020) / 0.99020;
  252. list err;
  253.  
  254. si (> err 1.0E-3);
  255. erreur(5);
  256. finsi;
  257.  
  258. fin;
  259.  
  260.  
  261.  
  262.  
  263.  
  264.  
  265.  
  266.  
  267.  
  268.  
  269.  
  270.  
  271.  
  272.  
  273.  
  274.  
  275.  
  276.  
  277.  

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