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