Télécharger test3_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 contraint courant par rapport à la surface
  5. * de capacité
  6. *
  7. * Test pour l'option pouteau long
  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 'P_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. * Test2 N = Nybl+ et My = Mybl+ et Mz = 0
  215. ******************************************************************
  216. *
  217. TB1.'NN' = NBIL_YP;
  218. TB1.'MY' = MBIL_YP;
  219. VAL3 = G_ULTIFR TB1;
  220. *
  221. * Test3 N = Nzbl+ et M0 = 0 et Mz = Mzbl+
  222. ******************************************************************
  223. *
  224. TB1.'NN' = NBIL_ZP;
  225. TB1.'MZ' = MBIL_ZP;
  226. TB1.'MY' = 0.0;
  227. VAL4 = G_ULTIFR TB1;
  228. *
  229. * Test2 N = Nybl- et My = Mybl- et Mz = 0
  230. ******************************************************************
  231. *
  232. TB1.'NN' = NBIL_YN;
  233. TB1.'MY' = MBIL_YN;
  234. TB1.'MZ' = 0.0;
  235. VAL5 = G_ULTIFR TB1;
  236. *
  237. * Test3 N = Nzbl- et M0 = 0 et Mz = Mzbl-
  238. ******************************************************************
  239. *
  240. TB1.'NN' = NBIL_ZN;
  241. TB1.'MZ' = MBIL_ZN;
  242. TB1.'MY' = 0.0;
  243. VAL6 = G_ULTIFR TB1;
  244. *
  245. * Test3 N = Nzbl- et M0 = Mybl- et Mz = Mzbl-
  246. ******************************************************************
  247. *
  248. TB1.'NN' = NBIL_ZN;
  249. TB1.'MZ' = MBIL_ZN;
  250. TB1.'MY' = MBIL_YN;
  251. VAL7 = G_ULTIFR TB1;
  252. * Test
  253. *
  254. SI ((ABS(VAL1)) > 1.E-8);
  255. ERRE 5;
  256. FINSI;
  257. *
  258. SI ((ABS(VAL2)) > 1.E-8);
  259. ERRE 5;
  260. FINSI;
  261. *
  262. SI ((ABS(VAL3)) > 1.E-8);
  263. ERRE 5;
  264. FINSI;
  265. *
  266. SI ((ABS(VAL4)) > 1.E-8);
  267. ERRE 5;
  268. FINSI;
  269. *
  270. SI ((ABS(VAL5)) > 1.E-8);
  271. ERRE 5;
  272. FINSI;
  273. *
  274. SI ((ABS(VAL6)) > 1.E-8);
  275. ERRE 5;
  276. FINSI;
  277. *
  278. SI (((VAL7)) > 1.E-8);
  279. ERRE 5;
  280. FINSI;
  281. *
  282. FIN;
  283.  
  284.  
  285.  

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