Télécharger la plaquette Cast3M


fr en

Télécharger test2_fun_gultifr.dgibi

Retour à la liste

Numérotation des lignes : 

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

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