Télécharger la plaquette Cast3M


fr en

Télécharger test4_fun_gultifr.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : test4_fun_gultifr.dgibi
  2. *
  3. **********************************************************
  4. *  Test sur la procedure G_ULTIFR (fonction pour determiner 
  5. *  la position de l'etat de contrainte courant par rapport à la surface 
  6. *  de capacité
  7. *
  8. *  Test pour l'option poutre longue
  9. *
  10. *  b=0.3
  11. *  h=0.4
  12. *  c=0.025
  13. *  armature long 1f20 + 1f30 + 1f12 + 1f12
  14. *  cadre phi 8
  15. *  s = 0.1 m
  16. *  ly=2.0
  17. *  lz=1.5
  18. *  
  19. *
  20. *  Develloppé par Alberto FRAU /DEN/DANS/DM2S/SEMT/EMSI
  21. *              et Nicolas ILE  /DEN/DANS/DM2S/SEMT/EMSI
  22. *
  23. **********************************************************
  24. *
  25. *dimensions de la section
  26. *Y horizontal
  27. *Z vertical
  28. B_Y1 = 0.3;
  29. B_Z1 = 0.4;
  30. *
  31. *Ltot = longueur libre d'inflexion
  32. LTOT1Y = 2.0;
  33. LTOT1Z = 1.5;
  34. *distance entre les cadres
  35. S_CAD1 = 0.1;
  36. *enrobage
  37. ENR1 = 0.025;
  38. *
  39. *diametre ferraillage longitudinal 
  40. PHI_LON1 = PROG 20.0 30.0 12.0 12.0;
  41. PHI_Y = PROG ((-1.0)*((B_Y1 - (2.0*ENR1))/(2.0))) 
  42.              ((+1.0)*((B_Y1 - (2.0*ENR1))/(2.0)))
  43.              ((+1.0)*((B_Y1 - (2.0*ENR1))/(2.0)))
  44.              ((-1.0)*((B_Y1 - (2.0*ENR1))/(2.0)));
  45. PHI_Z = PROG ((-1.0)*((B_Z1 - (2.0*ENR1))/(2.0))) 
  46.              ((-1.0)*((B_Z1 - (2.0*ENR1))/(2.0)))
  47.              ((+1.0)*((B_Z1 - (2.0*ENR1))/(2.0)))
  48.              ((+1.0)*((B_Z1 - (2.0*ENR1))/(2.0)));
  49. *
  50.  
  51. * A tot des cadres selon Y et selon Z
  52. A_CADRE1 = (2.0)*((PI)*((0.004)**(2.0)));
  53. A_CADRE2 = (2.0)*((PI)*((0.004)**(2.0)));
  54. *
  55. *A longitudinale
  56. ALL1 = 0.0;
  57. I = 1;
  58. NN1 = DIME PHI_LON1;
  59. *
  60. REPE IND1 (NN1);
  61.   PHI_I1 = EXTR PHI_LON1 I;
  62.   AA1 = (PI)*(((PHI_I1)*(PHI_I1))/(4.0));
  63.   AA1 = AA1/(1.E6);
  64.   ALL1 = ALL1 + AA1;
  65.   I = I + 1;
  66. FIN IND1;
  67. *
  68. * omsy1 et omsz
  69. OM_SY1 = ((ALL1)*(500.E6))/(((B_Y1)*(B_Z1))*(25.E6));
  70. OM_SZ1 = ((ALL1)*(500.E6))/(((B_Y1)*(B_Z1))*(25.E6));
  71. *
  72. * omwy1 et omwz
  73. OM_WZ1 = ((A_CADRE1)*(500.E6))/(((B_Y1)*(S_CAD1))*(25.E6));
  74. OM_WY1 = ((A_CADRE2)*(500.E6))/(((B_Z1)*(S_CAD1))*(25.E6));
  75. *
  76. * lamy1 et lamy2
  77. LAMB_Y1 = ((LTOT1Z)/(B_Y1));
  78. LAMB_Z1 = ((LTOT1Y)/(B_Z1));
  79. *
  80. * test court ou long
  81. VAL_TT1 = LAMB_Y1 <EG (OM_SY1/OM_WY1);
  82. VAL_TT2 = LAMB_Z1 <EG (OM_SZ1/OM_WZ1);
  83. *
  84. *si la relation est satisfaite --> error
  85. SI ((VAL_TT1) ET (VAL_TT2));
  86.   ERRE 5;
  87. FINSI;
  88. *
  89. *N traction
  90. N_TRAC1 =((ALL1)*(500.E6));
  91. *N compression
  92. N_COMP1 = (((B_Y1)*(B_Z1))*(25.E6))+((ALL1)*(500.E6));
  93. N_COMP1 = ((-1.0)*(N_COMP1));
  94. *
  95. *DEFORMATION DU BETON
  96. EPSI_BET = 0.0035;
  97. *DEFORMATION DE L'ACIER
  98. EPSI_ACI = 0.045;
  99. *
  100. *HAUTEUR NETTE
  101. D_YP = B_Z1-ENR1;
  102. D_ZP = B_Y1-ENR1;
  103. D_YN = (B_Z1-ENR1);
  104. D_ZN = (B_Y1-ENR1);
  105. *
  106. *position de L'AXE NEUTRE
  107. AXEN_YP = ((EPSI_BET)/(EPSI_BET+EPSI_ACI))*D_YP;
  108. AXEN_ZP = ((EPSI_BET)/(EPSI_BET+EPSI_ACI))*D_ZP;
  109. AXEN_YN = ((EPSI_BET)/(EPSI_BET+EPSI_ACI))*D_YN;
  110. AXEN_ZN = ((EPSI_BET)/(EPSI_BET+EPSI_ACI))*D_ZN;
  111. *
  112. FF_YP = ((EPSI_BET)/(AXEN_YP));
  113. FF_ZP = ((EPSI_BET)/(AXEN_ZP));
  114. FF_YN = (-1.0)*((EPSI_BET)/(AXEN_YN));
  115. FF_ZN = (-1.0)*((EPSI_BET)/(AXEN_ZN));
  116. EPS0_YP = (FF_YP)*((B_Z1/2.0) - AXEN_YP);
  117. EPS0_ZP = (FF_ZP)*((B_Y1/2.0) - AXEN_ZP);
  118. EPS0_YN = (FF_YN)*(AXEN_YN - (B_Z1/2.0));
  119. EPS0_ZN = (FF_ZN)*(AXEN_ZN - (B_Y1/2.0));
  120. *
  121. *N BALANCED
  122. NBIL_YP = ((-0.8)*(25.E6))*((B_Y1)*(AXEN_YP));
  123. NBIL_YN = ((-0.8)*(25.E6))*((B_Y1)*(AXEN_YN));
  124. NBIL_ZP = ((-0.8)*(25.E6))*((B_Z1)*(AXEN_ZP));
  125. NBIL_ZN = ((-0.8)*(25.E6))*((B_Z1)*(AXEN_ZN));
  126. *M BALANCED
  127. MBIL_YP = ((1.0)*(NBIL_YP))*((B_Z1/2.0) - ((0.4)*(AXEN_YP)));
  128. MBIL_ZP = ((1.0)*(NBIL_ZP))*((B_Y1/2.0) - ((0.4)*(AXEN_ZP)));
  129. MBIL_YN = ((-1.0)*(NBIL_YN))*((B_Z1/2.0) - ((0.4)*(AXEN_YN)));
  130. MBIL_ZN = ((-1.0)*(NBIL_ZN))*((B_Y1/2.0) - ((0.4)*(AXEN_ZN)));
  131. *
  132. I = 1;
  133.  
  134. REPE IND1 (NN1);
  135. * Calcolo di epsi
  136.   EPSI_YP = EPS0_YP - ((FF_YP)*(EXTR PHI_Z I));
  137.   EPSI_ZP = EPS0_ZP - ((FF_ZP)*(EXTR PHI_Y I));
  138.   EPSI_YN = EPS0_YN - ((FF_YN)*(EXTR PHI_Z I));
  139.   EPSI_ZN = EPS0_ZN - ((FF_ZN)*(EXTR PHI_Y I));
  140. *calcolo di sigma
  141.   SIGI_YP = (210000.0E6)*(EPSI_YP);
  142.   SIGI_YN = (210000.0E6)*(EPSI_YN);
  143.   SIGI_ZP = (210000.0E6)*(EPSI_ZP);
  144.   SIGI_ZN = (210000.0E6)*(EPSI_ZN);
  145.   SI ((ABS(SIGI_YP)) > 500.E6);
  146.     SIGI_YP = (SIGI_YP/(ABS(SIGI_YP)))*(500.E6);
  147.   FINSI;
  148.   SI ((ABS(SIGI_YN)) > 500.E6);
  149.     SIGI_YN = (SIGI_YN/(ABS(SIGI_YN)))*(500.E6);
  150.   FINSI;
  151.   SI ((ABS(SIGI_ZP)) > 500.E6);
  152.     SIGI_ZP = (SIGI_ZP/(ABS(SIGI_ZP)))*(500.E6);
  153.   FINSI;
  154.   SI ((ABS(SIGI_ZN)) > 500.E6);
  155.     SIGI_ZN = (SIGI_ZN/(ABS(SIGI_ZN)))*(500.E6);
  156.   FINSI;
  157. *
  158.   PHI_I1 = EXTR PHI_LON1 I;
  159.   AA1 = (PI)*(((PHI_I1)*(PHI_I1))/(4.0));
  160.   AA1 = AA1/(1.E6);
  161. *
  162.   NBIL_YP = NBIL_YP + ((AA1)*(SIGI_YP));
  163.   NBIL_YN = NBIL_YN + ((AA1)*(SIGI_YN));
  164.   NBIL_ZP = NBIL_ZP + ((AA1)*(SIGI_ZP));
  165.   NBIL_ZN = NBIL_ZN + ((AA1)*(SIGI_ZN));
  166. *
  167.   MBIL_YP = MBIL_YP + (((AA1)*(SIGI_YP))*(EXTR PHI_Z I));
  168.   MBIL_YN = MBIL_YN + (((AA1)*(SIGI_YN))*(EXTR PHI_Z I));
  169.   MBIL_ZP = MBIL_ZP + (((AA1)*(SIGI_ZP))*(EXTR PHI_Y I));
  170.   MBIL_ZN = MBIL_ZN + (((AA1)*(SIGI_ZN))*(EXTR PHI_Y I));
  171. *
  172. I = I + 1;
  173. FIN IND1;
  174. *
  175. * Test1 N = Ncomp
  176. ******************************************************************
  177. TB1                     = TABLE;
  178. TB1.'TYPE'              = CHAINE 'PT_LONG';
  179. TB1.'NN'                = N_COMP1;
  180. TB1.'VY'                = 0.0;
  181. TB1.'VZ'                = 0.0;
  182. TB1.'MT'                = 0.0;
  183. TB1.'MY'                = 0.0;
  184. TB1.'MZ'                = 0.0;
  185. TB1.'BY'                = B_Y1;
  186. TB1.'BZ'                = B_Z1;
  187. TB1.'FCD'               = 25.e6;
  188. TB1.'FSD'               = 500.e6;
  189. TB1.'WSY'               = OM_SY1;
  190. TB1.'WSZ'               = OM_SZ1;
  191. TB1.'WWY'               = OM_WY1;
  192. TB1.'WWZ'               = OM_WZ1;
  193. TB1.'LY'                = LAMB_Y1;
  194. TB1.'LZ'                = LAMB_Z1;
  195. TB1.'XIY'               = CHI_Y1;
  196. TB1.'XIZ'               = CHI_Z1;
  197. TB1.'PHI_LON'           = PHI_LON1;
  198. TB1.'PHI_Y'             = PHI_Y;
  199. TB1.'PHI_Z'             = PHI_Z;
  200. TB1.'EPSI_BET'          = EPSI_BET;
  201. TB1.'EPSI_ACI'          = EPSI_ACI;
  202. TB1.'YOUNG_ACIER'       = 210000.E6;
  203. *
  204. ******************************************************************
  205. *
  206. VAL1 = G_ULTIFR TB1;
  207. *
  208. * Test2 N = Ntrac
  209. ******************************************************************
  210. *
  211. TB1.'NN'                = N_TRAC1;
  212. TB1.'VY'                = 0.0;
  213. VAL2 = G_ULTIFR TB1;
  214. *
  215. * Test3 N = Nybl+ et My = Mybl+
  216. ******************************************************************
  217. *
  218. TB1.'NN'                = NBIL_YP;
  219. TB1.'MY'                = MBIL_YP;
  220. VAL3 = G_ULTIFR TB1;
  221. *
  222. * Test4 N = Nybl- et My = Mybl-
  223. ******************************************************************
  224. *
  225. TB1.'NN'                = NBIL_YN;
  226. TB1.'MY'                = MBIL_YN;
  227. TB1.'MZ'                = 0.0;
  228. VAL4 = G_ULTIFR TB1;
  229. *  Test 
  230. *
  231. SI ((ABS(VAL1)) > 1.E-8);
  232.   ERRE 5;
  233. FINSI;
  234. *
  235. SI ((ABS(VAL2)) > 1.E-8);
  236.   ERRE 5;
  237. FINSI;
  238. *
  239. SI ((ABS(VAL3)) > 1.E-8);
  240.   ERRE 5;
  241. FINSI;
  242. *
  243. SI ((ABS(VAL4)) > 1.E-8);
  244.   ERRE 5;
  245. FINSI;
  246. *
  247. FIN;
  248.  
  249.  
  250.  
  251.  
  252.  

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