Télécharger la plaquette Cast3M


fr en

Télécharger ricbet_uni_2.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  ricbet_uni_2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ******************************************************
  5. * DESCRIPTION
  6. ******************************************************
  7. * Local test - Multifiber analysis
  8. *
  9. * Loading path considered:
  10. *
  11. * - 1 : Monotonic tension test
  12. * - 2 : Monotonic compression test
  13. * - 3 : Cyclic tension test
  14. * - 4 : Cyclic compression test
  15. * - 5 : Cyclic tension/compression test
  16. * - 6 : Tension/compression/tension test
  17. *
  18. *
  19. ******************************************************
  20. * AUTHOR
  21. ******************************************************
  22. * Developped by : Benjamin RICHARD 
  23. *                 CEA-DEN/DANS/DM2S/SEMT/EMSI
  24. *                 Benjamin.Richard@cea.fr
  25. ******************************************************
  26. graph = 'N';
  27. ******************************************************
  28. * Type of loading 
  29. ncas  = 1; 
  30. ****************************************************** 
  31.  
  32. ******************************************************
  33. * OPTIONS
  34. ******************************************************
  35. OPTION DIME 2 ELEM QUA4;
  36.  
  37. * Number of element------------------------------------
  38. nelem = 10;
  39.  
  40. * Concrete parameters----------------------------------
  41.  
  42. * Young
  43. ESSS= 28000E6;
  44.  
  45. * Poisson ratio
  46. nun = 0.20;
  47.  
  48. * Tension strength
  49. C1  = 2.1E+6;
  50.  
  51. * Brittleness in tension
  52. C2  = 9.0E-3;
  53.  
  54. * Kinematic hardening coefficient 1
  55. C3  = 3.0E+9;
  56.  
  57. * Kinematic hardening coefficient 2
  58. C4  = 8.0e-6;
  59.  
  60. * Closure stress
  61. C5  = -3.0e6;
  62.  
  63. * Resistance en compession
  64. C6  = 10.0e6;
  65.  
  66. * XXX (unused)
  67. C7  = 1.0;
  68.  
  69. * XXX (unused)
  70. C8  = 1.0;
  71.  
  72. * Yield surface 1
  73. C9  = 4.0e10;
  74.  
  75. * Yield surface 2
  76. C10  = 600.0;
  77.  
  78. * Parameter that can be implemented in ESOPE 
  79. * (not subjected to high variations)
  80. C11  = -6.0e6;
  81.  
  82. ******************************************************
  83. * GEOMETRICAL PARAMETERS
  84. ******************************************************
  85. SACIER1 = 0. 0.;
  86. SACIER2 = 1. 0.;
  87. SACIER3 = 1. 1.;
  88. SACIER4 = 0. 1.;
  89.  
  90. LL1 = D 1 SACIER1 SACIER2;
  91. LL2 = D 1 SACIER2 SACIER3;
  92. LL3 = D 1 SACIER3 SACIER4;
  93. LL4 = D 1 SACIER4 SACIER1;
  94.  
  95. SACIER  = DALL LL1 LL2 LL3 LL4;
  96.  
  97. ******************************************************
  98. * MODEL
  99. ******************************************************
  100. OPTI DIME 3 ELEM SEG2;
  101.  
  102. *mo1 = mode sacier mecanique elastique quas;
  103.  
  104. *ma1 = MATE mo1 YOUN ESSS NU nun;
  105.  
  106. mo1 = mode sacier mecanique elastique plastique RICBET_UNI 
  107.       quas;
  108.  
  109. ma1 = MATE mo1 YOUN ESSS NU nun
  110. 'FT  ' C1  'ALDI' C2  'GAM1' C3
  111. 'A1  ' C4  'SIGF' C5  'FC  ' C6
  112. 'AF  ' C7  'AG  ' C8  'AC  ' C9
  113. 'BC  ' C10 'SIGU' C11 'HYST'  2;
  114.  
  115. car1 = CARA mo1 ALPY .66 ALPZ .66;
  116.  
  117. modstot = mo1;
  118. matstot = ma1 et car1;
  119.  
  120. ******************************************************
  121. * BEAM ELEMENT
  122. ******************************************************
  123. OPTI ELEM SEG2;
  124.  
  125. P1 = 0. 0. 0.;
  126. P2 = 0. 0. 1.;
  127.  
  128. Lpoutre = d nelem P1 P2;
  129.  
  130. modpout = MODE Lpoutre MECANIQUE ELASTIQUE SECTION 
  131.                PLASTIQUE SECTION TIMO;
  132. matpout = MATE modpout MODS modsTOT MATS matsTOT 
  133.                VECT (0 1 0);
  134.  
  135. ******************************************************
  136. * BOUNDARY CONDITIONS
  137. ******************************************************
  138. con1 = BLOQ DEPL ROTA         P1;
  139. con3 = BLOQ UZ                P2;
  140. con2 = BLOQ UX UY RX RY RZ P2 ;
  141.  
  142. ******************************************************
  143. * LOADING CASES
  144. ******************************************************
  145. fo1  = depi con3 1.0;
  146.  
  147. SI (EGA ncas 1) ; 
  148. LI1  = PROG 0. 1.; 
  149. LI2  = PROG 0. 1.0e-3; 
  150. LIS1 = PROG 0. PAS 0.02 1.; 
  151. FINSI; 
  152.  
  153. SI (EGA ncas 2) ; 
  154. LI1  = PROG 0.  1.; 
  155. LI2  = PROG 0. -5.0e-3;
  156. LIS1 = PROG 0. PAS 0.02 1; 
  157. FINSI; 
  158.  
  159. SI (EGA ncas 3) ; 
  160. LI1  = PROG 0. 1. 2. 3. 4. 5.; 
  161. LI2  = PROG 0. 1.5e-4 9.5e-6 2.0E-4 3.0E-5 10E-4;
  162. LIS1 = PROG 0. PAS 0.02 5.; 
  163. FINSI; 
  164.  
  165. SI (EGA ncas 4) ; 
  166. LI1  = PROG 0. 1. 2. 3. 4. 5.; 
  167. LI2  = PROG 0. -3.0e-3 1.5e-4 -5.0E-3 -3.5E-4 -8.0E-3; 
  168. LIS1 = PROG 0. PAS 0.02 5.; 
  169. FINSI; 
  170.  
  171. SI (EGA ncas 5) ; 
  172. LI1  = PROG 0. 1. 2. 3. 4. 5.; 
  173. LI2  = PROG 0. 3.e-4 -3.0e-4 5.0E-4 -1.5E-4 1.0E-3; 
  174. LIS1 = PROG 0. PAS 0.02 5; 
  175. FINSI;
  176.  
  177. SI (EGA ncas 6) ; 
  178. LI1  = PROG 0. 1. 2. 3. 4. 5.; 
  179. LI2  = PROG 0. 3.e-4 -1.e-5 5.0E-4 -5E-3 -3.0E-3; 
  180. LIS1 = PROG 0. PAS 0.02 5; 
  181. FINSI; 
  182.  
  183. ev1  = evol manu li1 li2;
  184. cha1 = char dimp fo1 ev1;
  185.  
  186. ******************************************************
  187. * SOLVING PROCEDURE
  188. ******************************************************
  189. TAB1 = TABLE;
  190. TAB1 . BLOCAGES_MECANIQUES = con1 et con3 et con2;
  191. TAB1 . MODELE              = modpout;
  192. TAB1 . CARACTERISTIQUES    = matpout;
  193. TAB1 . CHARGEMENT          = cha1;
  194. TAB1 . TEMPS_CALCULES      = lis1;
  195. TAB1 . MAXDEFOR            = 1.E-5;
  196. TAB1 . MOVA                = RIEN;
  197.  
  198. PASAPAS TAB1;
  199.  
  200. ******************************************************
  201. * POST PROCESSING
  202. ******************************************************
  203. dep0 = TAB1.'DEPLACEMENTS';
  204. tim0 = tab1. TEMPS;
  205. sig0 = TAB1.'CONTRAINTES';
  206. var0 = TAB1.'VARIABLES_INTERNES';
  207. ind0 = INDE dep0;
  208. time = PROG (tim0 . 0);
  209.  
  210. uzp0  = PROG 0.;
  211. fzp0  = PROG 0.;
  212.  
  213. ep1 = prog 0.;
  214. si1 = prog 0.;
  215. nn = (('DIME' ind0) - 1);
  216.  
  217. j = 1;
  218. REPE lab1 (('DIME' ind0) - 1);
  219.   j = j + 1;
  220.   indj  = ind0.j;
  221.   depj  = dep0.indj;
  222.   sigj  = sig0.indj;
  223.   epsj  = epsi depj modpout matpout;
  224.   vecj  = bsigma sigj modpout matpout;
  225.   time  = time et (prog tim0. &lab1);
  226.   uzp0  = uzp0 ET (PROG (EXTR depj 'UZ' P2));
  227.   fzp0  = fzp0 ET (PROG (EXTR vecj 'FZ' P2));
  228. FIN lab1;
  229.  
  230. si (ega graph 'O');
  231. evuzfz = evol manu uzp0 (fzp0);
  232. list evuzfz;
  233. dess evuzfz;
  234. finsi;
  235.  
  236. aa = (EXTR vecj 'FZ' P2);
  237.  
  238. message 'erreur relative';
  239. err = (aa - 3.46710E+05) / 3.46710E+05;
  240. list err;
  241.  
  242. si (> err 1E-3);
  243. erreur(5);
  244. finsi;
  245.  
  246. fin;
  247.  
  248.  
  249.  
  250.  
  251.  
  252.  
  253.  
  254.  
  255.  
  256.  
  257.  
  258.  
  259.  
  260.  
  261.  
  262.  
  263.  
  264.  
  265.  
  266.  
  267.  
  268.  
  269.  

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