Télécharger FissVoil.dgibi

Retour à la liste

Numérotation des lignes :

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

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