Télécharger test4_fun_gultifr.dgibi

Retour à la liste

Numérotation des lignes :

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

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