Télécharger test1_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 court
  8. *
  9. * b=0.3
  10. * h=0.4
  11. * c=0.025
  12. * armature long 4f12
  13. * cadre phi 8
  14. * s = 0.95 m
  15. * ly=1.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. * donnée d'entrée
  25. B_Y1 = 0.3;
  26. B_Z1 = 0.4;
  27. LTOT1Y = 1.0;
  28. LTOT1Z = 1.5;
  29. S_CAD1 = 0.95;
  30. PHI_LON1 = 12.0;
  31. ENR1 = 0.025;
  32. A_CADRE1 = (2.0)*((PI)*((0.004)**(2.0)));
  33. A_CADRE2 = (2.0)*((PI)*((0.004)**(2.0)));
  34. *
  35. * calcul des armatures longitudinales
  36. ALL1 = (4.0)*((PI)*(((PHI_LON1*PHI_LON1)/(1.E6))/(4.0)));
  37. * om_sy et om_sz
  38. OM_SY1 = ((ALL1)*(500.E6))/(((B_Y1)*(B_Z1))*(25.E6));
  39. OM_SZ1 = ((ALL1)*(500.E6))/(((B_Y1)*(B_Z1))*(25.E6));
  40. * om_wy et om_wz
  41. OM_WZ1 = ((A_CADRE1)*(500.E6))/(((B_Y1)*(S_CAD1))*(25.E6));
  42. OM_WY1 = ((A_CADRE2)*(500.E6))/(((B_Z1)*(S_CAD1))*(25.E6));
  43. * lamy ety lamz
  44. LAMB_Y1 = ((LTOT1Z)/(B_Y1));
  45. LAMB_Z1 = ((LTOT1Y)/(B_Z1));
  46. * chiy et chiz
  47. CHI_Z1 = ((B_Z1 - ((2.0)*(ENR1)))/(B_Z1));
  48. CHI_Y1 = ((B_Y1 - ((2.0)*(ENR1)))/(B_Y1));
  49. * test pour savoir si on est dans le cas du poteau court ou long
  50. VAL_TT1 = LAMB_Y1 <EG (OM_SY1/OM_WY1);
  51. VAL_TT2 = LAMB_Z1 <EG (OM_SZ1/OM_WZ1);
  52. * calcul de Vymax et Vzmax
  53. TAN_Y1 = (((LAMB_Y1**2.0) + 1.0)**(0.5)) - LAMB_Y1;
  54. TAN_Z1 = (((LAMB_Z1**2.0) + 1.0)**(0.5)) - LAMB_Z1;
  55. VY_MAX1 = OM_WY1*CHI_Y1 + ((0.5 - OM_WY1)*(TAN_Y1));
  56. VZ_MAX1 = OM_WZ1*CHI_Z1 + ((0.5 - OM_WZ1)*(TAN_Z1));
  57. VY_MAX1 = ((VY_MAX1)*(((B_Y1)*(B_Z1))*(25.e6)));
  58. VZ_MAX1 = ((VZ_MAX1)*(((B_Y1)*(B_Z1))*(25.e6)));
  59. *
  60. * calcul des differentes valeur de N pour les trois zone dans le
  61. * diagramme d'interaction
  62. N_TRAC1 = ((OM_SY1)*(((B_Y1)*(B_Z1))*(25.e6)));
  63. N_COMP1 = ((1.0 + OM_SY1)*(((B_Y1)*(B_Z1))*(-25.e6)));
  64. N_LIMY1 = 0.5 + OM_SY1 - ((OM_WY1)*(LAMB_Y1 + 1.0 - CHI_Y1));
  65. N_LIMY2 = 0.5 - OM_SY1 + ((OM_WY1)*(LAMB_Y1 - 1.0 + CHI_Y1));
  66. N_LIMZ1 = 0.5 + OM_SZ1 - ((OM_WZ1)*(LAMB_Z1 + 1.0 - CHI_Z1));
  67. N_LIMZ2 = 0.5 - OM_SZ1 + ((OM_WZ1)*(LAMB_Z1 - 1.0 + CHI_Z1));
  68. N_LIMY1 = ((N_LIMY1)*(((B_Y1)*(B_Z1))*(-25.e6)));
  69. N_LIMY2 = ((N_LIMY2)*(((B_Y1)*(B_Z1))*(-25.e6)));
  70. N_LIMZ1 = ((N_LIMZ1)*(((B_Y1)*(B_Z1))*(-25.e6)));
  71. N_LIMZ2 = ((N_LIMZ2)*(((B_Y1)*(B_Z1))*(-25.e6)));
  72. *
  73. * Test1 (Zone 2) Vy= Vymax N=(N1y+N2y)/2
  74. ******************************************************************
  75. TB1 = TABLE;
  76. TB1.'TYPE' = CHAINE 'P_COURT';
  77. TB1.'NN' = (N_LIMY1 + N_LIMY2)/(2.0);
  78. TB1.'VY' = VY_MAX1;
  79. TB1.'VZ' = 0.0;
  80. TB1.'MT' = 0.0;
  81. TB1.'MY' = 0.0;
  82. TB1.'MZ' = 0.0;
  83. TB1.'BY' = B_Y1;
  84. TB1.'BZ' = B_Z1;
  85. TB1.'FCD' = 25.e6;
  86. TB1.'FSD' = 500.e6;
  87. TB1.'WSY' = OM_SY1;
  88. TB1.'WSZ' = OM_SZ1;
  89. TB1.'WWY' = OM_WY1;
  90. TB1.'WWZ' = OM_WZ1;
  91. TB1.'LY' = LAMB_Y1;
  92. TB1.'LZ' = LAMB_Z1;
  93. TB1.'XIY' = CHI_Y1;
  94. TB1.'XIZ' = CHI_Z1;
  95. *
  96. * Test2 (Zone 2) Vy= Vymax N=(N1y) et N=(N2y)
  97. ******************************************************************
  98. *
  99. VAL1Y = G_ULTIFR TB1;
  100. TB1.'NN' = (N_LIMY1);
  101. TB1.'VY' = VY_MAX1;
  102. VAL2Y = G_ULTIFR TB1;
  103. TB1.'NN' = (N_LIMY2);
  104. TB1.'VY' = VY_MAX1;
  105. VAL3Y = G_ULTIFR TB1;
  106. *
  107. * Test3 (Zone 2) Vz= Vzmax N=(N1z+N2z)/2.0
  108. ******************************************************************
  109. *
  110. TB1.'NN' = (N_LIMZ1 + N_LIMZ2)/(2.0);
  111. TB1.'VY' = 0.0;
  112. TB1.'VZ' = VZ_MAX1;
  113. *
  114. VAL1Z = G_ULTIFR TB1;
  115. *
  116. * Test4 (Zone 2) Vz= Vzmax N=(N1z) et N=N2z
  117. TB1.'NN' = (N_LIMZ1);
  118. VAL2Z = G_ULTIFR TB1;
  119. *
  120. TB1.'NN' = (N_LIMZ2);
  121. VAL3Z = G_ULTIFR TB1;
  122. *
  123. * Test5 N=Ntrac
  124. ******************************************************************
  125. *
  126. TB1.'NN' = N_TRAC1;
  127. TB1.'VY' = 0.0;
  128. TB1.'VZ' = 0.0;
  129. VAL4 = G_ULTIFR TB1;
  130. *
  131. * Test6 N=Ncomp
  132. ******************************************************************
  133. *
  134. TB1.'NN' = N_COMP1;
  135. TB1.'VY' = 0.0;
  136. TB1.'VZ' = 0.0;
  137. VAL5 = G_ULTIFR TB1;
  138. *
  139. *
  140. * controle
  141. SI ((ABS(VAL1Y)) > 1.E-8);
  142. ERRE 5;
  143. FINSI;
  144. *
  145. SI ((ABS(VAL2Y)) > 1.E-8);
  146. ERRE 5;
  147. FINSI;
  148. *
  149. SI ((ABS(VAL3Y)) > 1.E-8);
  150. ERRE 5;
  151. FINSI;
  152. *
  153. *
  154. *
  155. *
  156. SI ((ABS(VAL1Z)) > 1.E-8);
  157. ERRE 5;
  158. FINSI;
  159. *
  160. SI ((ABS(VAL2Z)) > 1.E-8);
  161. ERRE 5;
  162. FINSI;
  163. *
  164. SI ((ABS(VAL3Z)) > 1.E-8);
  165. ERRE 5;
  166. FINSI;
  167. *
  168. *
  169. *
  170. *
  171. SI ((ABS(VAL4)) > 1.E-8);
  172. ERRE 5;
  173. FINSI;
  174. SI ((ABS(VAL5)) > 1.E-8);
  175. ERRE 5;
  176. FINSI;
  177. *
  178. *
  179. *
  180. *
  181. FIN;
  182.  
  183.  
  184.  
  185.  
  186.  
  187.  

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