Télécharger ricjoi_3d.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : ricjoi_3d.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *******************************************************
  5. * Cas test de l'implantation numérique du modele
  6. * RICJOI_3D LOCAL
  7. *******************************************************
  8. * Développé par :
  9. * Benjamin Richard
  10. * Contact :
  11. * Benjamin.Richard@lmt.ens-cachan.fr
  12. *******************************************************
  13. * Les cas de charges sont entrés :
  14. * - 1 : Essai cyclique mono cadrant
  15. * - 2 : Essai cyclique bi cadrant (cisaillement)
  16. * - 3 : Essai effet unilatteral
  17. * - 4 : Essai frottement
  18. * - 5 : Essai variations stiffness
  19. * - 6 : Essai compression normale avec rouille monotone
  20. *******************************************************
  21. * Choix du cas de charge
  22. ncas = 1;
  23. graph = mot 'N';
  24. *
  25. *-------------- Options de calcul ---------------------
  26. *
  27. OPTI DIME 3 ELEM QUA4 MODE TRID ;
  28. * ajout de option epsilon lineaire pour la precision des test!
  29. OPTION epsilon lineaire;
  30. *
  31. *-------------- Definition de la geometrie ------------
  32. *
  33. A1 = 0.00 0.00 0.00 ;
  34. B1 = 1.00 0.00 0.00 ;
  35. C1 = 1.00 1.00 0.00 ;
  36. D1 = 0.00 1.00 0.00 ;
  37.  
  38. O1 = 0.00 0.00 10.00 ;
  39.  
  40. H1 = A1 DROIT 1 B1 ;
  41. H2 = B1 DROIT 1 C1 ;
  42. H3 = C1 DROIT 1 D1 ;
  43. H4 = D1 DROIT 1 A1 ;
  44. *
  45. L1 = H1 ET H2 ET H3 ET H4 ;
  46. S1 = SURF L1 PLAN ;
  47.  
  48. IA1 = 0.00 0.00 0.00 ;
  49. IB1 = 1.00 0.00 0.00 ;
  50. IC1 = 1.00 1.00 0.00 ;
  51. ID1 = 0.00 1.00 0.00 ;
  52.  
  53. IH1 = IA1 DROIT 1 IB1 ;
  54. IH2 = IB1 DROIT 1 IC1 ;
  55. IH3 = IC1 DROIT 1 ID1 ;
  56. IH4 = ID1 DROIT 1 IA1 ;
  57.  
  58. IL1 = IH1 ET IH2 ET IH3 ET IH4 ;
  59. IS1 = SURF IL1 PLAN ;
  60.  
  61. OPTION ELEM LIA4 ;
  62. VOL = LIAISON 0.00001 S1 IS1 ;
  63.  
  64. *
  65. *------- Définition des conditions aux limites --------
  66. *----------- et des déplacements imposés --------------
  67. *
  68. CL11 = BLOQ IA1 UX ;
  69. CL12 = BLOQ IA1 UY ;
  70. CL13 = BLOQ IA1 UZ ;
  71. CL1 = CL11 ET CL12 ET CL13 ;
  72.  
  73. CL21 = BLOQ IB1 UX ;
  74. CL22 = BLOQ IB1 UY ;
  75. CL23 = BLOQ IB1 UZ ;
  76. CL2 = CL21 ET CL22 ET CL23 ;
  77.  
  78. CL31 = BLOQ IC1 UX ;
  79. CL32 = BLOQ IC1 UY ;
  80. CL33 = BLOQ IC1 UZ ;
  81. CL3 = CL31 ET CL32 ET CL33 ;
  82.  
  83. CL41 = BLOQ ID1 UX ;
  84. CL42 = BLOQ ID1 UY ;
  85. CL43 = BLOQ ID1 UZ ;
  86. CL4 = CL41 ET CL42 ET CL43 ;
  87.  
  88. SI (OU (EGA ncas 3) (EGA ncas 6));
  89. CL5 = BLOQ S1 UZ;
  90. SINON;
  91. CL5 = BLOQ S1 UY;
  92. FINSI;
  93.  
  94. CL = CL1 ET CL2 ET CL3 ET CL4 ET CL5;
  95.  
  96. *
  97. *----- Définition du modèle ---------------------------
  98. *
  99. MOD1 = MODELI VOL 'MECANIQUE' ELASTIQUE 'PLASTIQUE' INTJOI JOI4 ;
  100.  
  101. *
  102. *----- Paramètres matériaux ---------------------------
  103. *
  104.  
  105. * Rigidite normale
  106. knn = (36000E+6);
  107.  
  108. * Rigidite tangentielle
  109. ksn = (knn/2.4);
  110.  
  111. * Coefficient de dilatation
  112. alpnn = 1.0E-10;
  113.  
  114. * Fragilite
  115. adn = 3.5E-5;
  116.  
  117. * Seuil initial en energie pour l'endommagement
  118. y0n = 50;
  119.  
  120. * Coefficient de dissymétrie entre modes
  121. alpan = 6;
  122.  
  123. * Coefficient d'écrouissage 1
  124. gainn = 2.0E9;
  125.  
  126. * Coefficient d'écrouissage 2
  127. aainn = 5.0E-7;
  128.  
  129. * Coefficient de Gurson modifié 1
  130. q1con = 3.5;
  131.  
  132. * Coefficient de Gurson modifié 2
  133. q2con = 0.9;
  134.  
  135. * Coefficient de Gurson modifié 3
  136. q3con = 0.1;
  137.  
  138. * Contrainte de raidissement seuil
  139. sycon = -1.0E6;
  140.  
  141. * Coefficient loi de contraction 1
  142. ncoen = 2.0;
  143.  
  144. * Coefficient loi de contraction 2
  145. kcoen = 1.0E+10;
  146.  
  147. * Degré de corrosion (en perte de section)
  148. SI(EGA ncas 6);
  149. tcn = 1.0E-10;
  150. SINON;
  151. tcn = 0;
  152. FINSI;
  153.  
  154. * Indicateur de phase de calcul
  155. gonfn = 0;
  156.  
  157. MA1 = MATE MOD1 KS ksn KN knn ALPN alpnn
  158. AD adn Y0 y0n ALPA alpan
  159. GAIN gainn AAIN aainn
  160. Q1CO q1con Q2CO q2con Q3CO q3con
  161. SYCO sycon NCOE ncoen KCOE kcoen
  162. TC tcn GONF gonfn;
  163.  
  164. *
  165. *----- Définition des cas de charge -------------------
  166. *
  167.  
  168. * Essai cyclique mono cadrant
  169. SI (EGA ncas 1);
  170. LI1 = PROG 0. 1. 2. 3. ;
  171. LI2 = PROG 0. 5.0E-3 0.0E0 10.0E-3 ;
  172. LIS1 = PROG 0.0 PAS 0.01 3.0;
  173. FINSI;
  174.  
  175. * Essai cyclique bi cadrant (cisaillement)
  176. SI (EGA ncas 2);
  177. LI1 = PROG 0. 1. 2. 3. 4. 5. 6. 7. 8. 9;
  178. LI2 = PROG 0. -3.0E-3 5.0E-3 -7.0E-3 9.0E-3 -11.0E-3 13.0E-3
  179. -15.0E-3 17E-3 0.0;
  180. LIS1 = PROG 0.0 PAS 0.01 9.0;
  181. FINSI;
  182.  
  183. * Essai effet unilatteral
  184. SI (EGA ncas 3);
  185. LI1 = PROG 0. 1. 2. 3.;
  186. LI2 = PROG 0. 5.0E-3 -1.0E-3 10E-3;
  187. LIS1 = PROG 0. PAS 0.01 3.0;
  188. FINSI;
  189.  
  190. * Essai frottement
  191. SI (EGA ncas 4);
  192. LI1 = PROG 0. pas 1. 7.;
  193. LI2 = PROG 0. 5.0E-3 0.0E-3 6e-3 0.0e-3 8e-3 0.0e-3 10E-3;
  194. LIS1 = PROG 0. PAS 0.01 7.0;
  195. FINSI;
  196.  
  197. * Essai variations stiffness
  198. SI (EGA ncas 5);
  199. LI1 = PROG 0. pas 1. 1.;
  200. LI2 = PROG 0. 10.0E-3 ;
  201. LIS1 = PROG 0. PAS 0.01 1.0;
  202. FINSI;
  203.  
  204. * Essai compression normale avec rouille monotone
  205. SI (EGA ncas 6);
  206. LI1 = PROG 0. pas 1. 1.;
  207. LI2 = PROG 0. -1.0E-4 ;
  208. LIS1 = PROG 0. PAS 0.01 1;
  209. FINSI;
  210.  
  211. RIG1 = DEPI CL5 1.0;
  212. EV1 = EVOL MANU LI1 LI2 ;
  213. CAR1 = CHAR 'DIMP' RIG1 EV1 ;
  214.  
  215. *
  216. *----------- Calcul par l'operateur PASAPAS ------------
  217. *
  218. TAB = TABLE ;
  219. TAB.'MOVA' = 'DINT';
  220. TAB.'BLOCAGES_MECANIQUES' = CL;
  221. TAB.'CARACTERISTIQUES' = MA1;
  222. TAB.'MODELE' = MOD1;
  223. TAB.'TEMPS_CALCULES' = LIS1;
  224. TAB.'CHARGEMENT' = CAR1;
  225. PASAPAS TAB ;
  226.  
  227. *
  228. *----------- Courbe effort-deplacement -----------------
  229. *
  230. dep0 = TAB.'DEPLACEMENTS';
  231. tim0 = TAB. TEMPS;
  232. sig0 = TAB.'CONTRAINTES';
  233. var0 = TAB.'VARIABLES_INTERNES';
  234. ind0 = INDE dep0;
  235. time = PROG (tim0 . 0);
  236.  
  237. uzp0 = PROG 0.;
  238. fzp0 = PROG 0.;
  239.  
  240. ep1 = prog 0.;
  241. si1 = prog 0.;
  242. nn = (('DIME' ind0) - 1);
  243.  
  244. j = 1;
  245. REPE lab1 (('DIME' ind0) - 1);
  246. j = j + 1;
  247. indj = ind0.j;
  248. depj = dep0.indj;
  249. sigj = sig0.indj;
  250. epsj = epsi depj mod1 ma1;
  251. vecj = bsigma sigj mod1 ma1;
  252. time = time et (prog tim0. &lab1);
  253. SI (OU (EGA ncas 3) (EGA ncas 6));
  254. uzp0 = uzp0 ET (PROG (EXTR depj 'UZ' A1));
  255. fzp0 = fzp0 ET (PROG (EXTR vecj 'FZ' IA1));
  256. SINON;
  257. uzp0 = uzp0 ET (PROG (EXTR depj 'UY' A1));
  258. fzp0 = fzp0 ET (PROG (EXTR vecj 'FY' IA1));
  259. FINSI;
  260. FIN lab1;
  261. aa = extr fzp0 ( dime fzp0);
  262. list aa;
  263. err = abs (aa + 1.85275E+06 ) / 1.85275E+06;
  264. message ' erreur relative ' err;
  265. si (err > 1.e-3);
  266. erreur (5);
  267. finsi;
  268. si ( ega graph 'O');
  269. EV2 = EVOL MANU uzp0 (-4.0*fzp0);
  270. DESS EV2;
  271. finsi;
  272. FIN ;
  273.  
  274.  
  275.  
  276.  
  277.  
  278.  
  279.  
  280.  
  281.  
  282.  
  283.  
  284.  
  285.  
  286.  

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