Télécharger test2_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 poutre courte
  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) Vz= Vzmax N=(N1z+N2z)/2.0
  74. ******************************************************************
  75. TB1 = TABLE;
  76. TB1.'TYPE' = CHAINE 'PT_COURT';
  77. TB1.'NN' = (N_LIMZ1 + N_LIMZ2)/(2.0);
  78. TB1.'VY' = 0.0;
  79. TB1.'VZ' = VZ_MAX1;
  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. VAL1Z = G_ULTIFR TB1;
  97. *
  98. * Test2 (Zone 2) Vz= Vzmax N=(N1z+N2z)/2.0
  99. TB1.'NN' = (N_LIMZ1);
  100. VAL2Z = G_ULTIFR TB1;
  101. *
  102. TB1.'NN' = (N_LIMZ2);
  103. VAL3Z = G_ULTIFR TB1;
  104. *
  105. * Test3 N=Ntrac
  106. ******************************************************************
  107. *
  108. TB1.'NN' = N_TRAC1;
  109. TB1.'VY' = 0.0;
  110. TB1.'VZ' = 0.0;
  111. VAL4 = G_ULTIFR TB1;
  112. *
  113. * Test4 N=Ntrac
  114. ******************************************************************
  115. *
  116. TB1.'NN' = N_COMP1;
  117. TB1.'VY' = 0.0;
  118. TB1.'VZ' = 0.0;
  119. VAL5 = G_ULTIFR TB1;
  120. *
  121. *
  122. *
  123. *
  124. SI ((ABS(VAL1Z)) > 1.E-8);
  125. ERRE 5;
  126. FINSI;
  127. *
  128. SI ((ABS(VAL2Z)) > 1.E-8);
  129. ERRE 5;
  130. FINSI;
  131. *
  132. SI ((ABS(VAL3Z)) > 1.E-8);
  133. ERRE 5;
  134. FINSI;
  135. *
  136. *
  137. *
  138. *
  139. SI ((ABS(VAL4)) > 1.E-8);
  140. ERRE 5;
  141. FINSI;
  142. SI ((ABS(VAL5)) > 1.E-8);
  143. ERRE 5;
  144. FINSI;
  145. *
  146. *
  147. *
  148. *
  149. FIN;
  150.  
  151.  
  152.  
  153.  
  154.  
  155.  

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