Télécharger la plaquette Cast3M


fr en

Télécharger FissVoil.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : FissVoil.dgibi
  2. **********************************
  3. *Section: Mecanique Nonlineaire  *
  4. *Section: Mecanique Transitoire  *
  5. **********************************
  6. graph = faux ;
  7. *
  8. * CAS TEST DU 13/11/14     PROVENANCE : TEST                                    
  9. SAUT PAGE;                                                                      
  10. *                                                                               
  11. *******************************************************        
  12. *                                                     *            
  13. *  TEST PROCEDURE OUVCOR POUR UN PANEL EN CISSAILEMENT*            
  14. *  4. * 1. * 0.2                                      *            
  15. *                                                     *             
  16. *                                                     *            
  17. *  Le maillage est en 3D a l'aide d'elements          *
  18. *  massif CUB8. L'acier est maille a l'aide d'elements*
  19. *  barre.                                             *
  20. *                                                     *
  21. *******************************************************            
  22.  
  23.  
  24. *******************************************************
  25. * MAILLAGE
  26. *******************************************************
  28. opti elem cub8 ;
  29.  
  30.  
  31. dfi = 0.15 ;
  32.  
  33. P1 = 0. 0. 0 ;
  34. P2 = 4. 0. 0 ;
  35.  
  36. DT = DROI P1 P2 dini dfi dfin dfi ;
  37.  
  38. SM = DT TRANS (0 1. 0) dini dfi dfin dfi ;
  39. DT = SM COTE 3 ;
  40.  
  41. VT = VOLU SM TRANS (0 0 0.2) dini dfi dfin dfi ;
  42.  
  43. elim 1.e-4 VT ;
  44. *** 
  45. *** definition des renforcements dans la pertie centrale
  46. ** barre verticale
  47. PA1 = 2. 0.05 0.05 ;
  48. PA2 = PA1 PLUS (0 (1. - 0.05) 0) ;
  49.  
  50. DV = DROI PA1 PA2 dini dfi dfin dfi ;
  51.  
  52. ** barres horizontales
  53.  
  54. P1 = 0.05 0.5 0.055 ;
  55. P2 = P1 PLUS ((4. - 0.05) 0 0) ;
  56.  
  57. DH = DROI P1 P2 dini dfi dfin dfi ;
  58.  
  59. STRUCT = vt et dv et dh ;
  60.  
  61. MODB = MODE VT MECANIQUE ELASTIQUE ISOTROPE 
  62.                  ENDOMMAGEMENT DAMAGE_TC ;
  63.  
  64. * aciers
  65.  
  66. MODA = MODE (DH ET DV) MECANIQUE ELASTIQUE ISOTROPE 
  67.                 PLASTIQUE ISOTROPE BARR ;
  68.  
  69. *
  70. MODTOT = MODB et MODA;
  71.  
  72. * propriétés mécaniques
  73.  
  74.  yg = 22.e9 ;
  75.  nub = 0.19 ;
  76.  rhom = 2410 ;
  77.  Gvalm = 300. ;
  78.  ftulm = 3.3e6 ;
  79.  redcm = 1.7e6 ;
  80.  fc01m = -25.e6 ;
  81.  rt45m = 1.18 ;
  82.  fcu1m = -42.5e6 ;
  83.  extum = -0.015;
  84.  strpm = -22.e6 ;
  85.  extpm = -0.001 ;
  86.  ext1m = -0.006 ;
  87.  str1m = -35.e6 ;
  88.  ext2m = -0.008 ;
  89.  str2m = -22.e6 ;
  90. ncrim = 1 ;
  91.  
  92. jaco1 = jaco modb ;
  93. jaco2 = chan 'RIGIDITE' modb jaco1 ;
  94. hlenm = jaco2**(1./3.) ; 
  95.  
  96. matb = mate modb  YOUN  yg NU nub RHO rhom HLEN hlenm
  97.                    GVAL Gvalm FTUL ftulm REDC redcm FC01 fc01m
  98.                    RT45 rt45m FCU1 fcu1m STRU extum EXTP extpm
  99.                    STRP strpm EXT1 ext1m STR1 str1m EXT2 ext2m
  100.                    STR2 str2m NCRI ncrim ;
  101.  
  102. * aciers renforcement 
  103.  
  104. EA = 189274e6 ;
  105. NUA = 0.3 ;
  106.  
  107. fe=554.e6;
  108. ETa=3245.e6;
  109.  
  110. epsel  = fe/Ea;
  111. epsfin = 2.7e-2 .;
  112. leps   = prog 0. epsel epsfin;
  113. sigfin = fe + (eta*(epsfin - epsel));
  114. lsig   = prog 0. fe sigfin;
  115. evsig  = evol manu leps lsig;
  116. lsig   =  lsig enle 1 ;
  117. leps   =  (leps enle 1) - (lsig / ea) ;
  118. evecr  =  evol vert manu eps leps sig lsig ;
  119. si graph ;
  120.   dess (evsig et evecr) titr 'Courbes de traction et d ecrouissage (vert)' ;
  121. fins ;
  122.  
  123.  
  124.  
  125. MATA = MATE MODA YOUN EA NU NUA RHO 7850 
  126.                    ECRO evecr SECT 7.85e-5  ;
  127.  
  128.  
  129. *
  130. MATOT = MATB et MATA;
  131.  
  132. *******************************************************
  133. * CONDITIONS AUX LIMITES
  134. *******************************************************
  135. p1 = vt poin proc (0 0 0);
  136. p2 = vt poin proc (1 0 0);
  137. p3 = vt poin proc (1 0 1);
  138. sb = vt poin plan p1 p2 p3 1e-4;
  139.  
  140. p1 = vt poin proc (p1 plus (0 1 0));
  141. p2 = vt poin proc (p2 plus (0 1 0));
  142. p3 = vt poin proc (p3 plus (0 1 0));
  143. sh = vt poin plan p1 p2 p3 1e-4;
  144.  
  145. pp1 = vt poin proc (0 0 0);
  146. pp2 = vt poin proc (1 0 0);
  147. pp3 = vt poin proc (0 4 0);
  148. ss = vt poin plan pp1 pp2 pp3 1e-4;
  149.  
  150. cl1 = BLOQ sh ux;
  151. cl2 = BLOQ sb ux uy;
  152. cl3 = BLOQ sh uy ;
  153. cl4 = BLOQ ss uz ;
  154. cl5 = RELA ACCRO (DV ET DH) VT  1.e-4 ;
  155.  
  156.  
  157. cltot = cl1 et cl2 et cl3 et cl4 et cl5;   
  158.  
  159.  
  160. * chargement en déplacement 
  161.  
  162. FO1 = DEPI cl1 (1) ;
  163.  
  164. list1 = prog 0 1;
  165. list2 = prog 0 1.5e-3 ;
  166.  
  167. evol1 = evol manu list1 list2 ;
  168.  
  169. char1= CHAR MECA FO1 evol1;
  170.  
  171. *******************************************************
  172. * RESOLUTION : PASAPAS
  173. *******************************************************
  174.  
  175. listt = prog 0 pas 0.05 1 ;
  176.  
  177. tab1 = TABLE;
  178.  
  179. tab1 . MODELE = modtot ;
  180.  
  181. tab1 . CARACTERISTIQUES = matot ;
  182.  
  183. tab1 . CHARGEMENT = char1 ;
  184.  
  185. tab1 . BLOCAGES_MECANIQUES = cltot ;
  186.  
  187. tab1 . TEMPS_CALCULES = listt ;
  188.  
  189. tab1 . MOVA = 'RIEN' ;
  190.  
  191.  
  192. PASAPAS tab1 ;
  193.  
  194. *******************************************************
  195. * POST-TRAITEMENT AVEC OUVCOR
  196. *******************************************************
  197.  
  198. * initou
  199. tab1.geo = vt;
  200. tab1.poi = vt poin proc (0. 0. 0.);
  201. tab1.hor = vrai;
  202. tab1.pla = 'XY';
  203. tab1.pas = 20;
  204. tab1.lh = 2.;
  205. tab1.lb = 0.;
  206. tab1.lg = 0.;
  207. tab1.ld = 4.;
  208. tab1.crito = 1.e-9;
  209. tab1.critp = 1.e-5;
  210.  
  211. initou tab1;
  212.  
  213.  
  214. * zonfis
  215. tab1.droi = faux;
  216. hau = 2.;
  217. bas = 0.;
  218. nha = 4;
  219. nba = 2;
  220. pha = 1;
  221. pba = 1;
  222. gau = -0.05;
  223. dro = 0.05;
  224. alpha = 0.;
  225. zonfis tab1 pha pba nha nba hau bas gau dro alpha;
  226.  
  227.  
  228. * postou
  229. dist = 0.03;
  230. tab1.critt = 0.5;
  231. *postou tab1 dist;
  232.  
  233. fin;
  234.  
  235.  
  236.  
  237.  
  238.  
  239.  
  240.  
  241.  
  242.  
  243.  

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