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