Télécharger la plaquette Cast3M


fr en

Télécharger mrcframe_test.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : mrcframe_test.dgibi
  2. *
  3. **********************************************************
  4. *  Test sur la procedure MRCFRAME, fonction pour determiner 
  5. *  la marge de securité d'un element type POUT et TIMO soumis à un 
  6. *  chargement sismique. Les effort sont calculées selon la methode
  7. *  spectrale avec et sans la prise en compte des termes de covariances
  8. *
  9. *
  10. *  Section
  11. *
  12. *  b=0.3
  13. *  h=0.4
  14. *  c=0.025
  15. *  armature long 4f12
  16. *  cadre phi 8
  17. *  At = f8 mm
  18. *
  19. *  4 Cast test:
  20. *    - V_TEST1 = 1 Poteau court         ly=1.0 lz=1.0 s = 0.95 m
  21. *    - V_TEST1 = 2 Poutre courte        ly=1.0 lz=1.0 s = 0.95 m
  22. *    - V_TEST1 = 3 Poteau long          ly=2.0 lz=1.5 s = 0.1 m
  23. *    - V_TEST1 = 4 Poutre longe         ly=2.0 lz=1.5 s = 0.1 m
  24. *
  25. *  Pour V_TEST1 = 1, 3 type de chargement:
  26. *     - Cha1: N0 = -1600.0 kN
  27. *             Nmax = 10 kN
  28. *             Vymax = 200 kNm
  29. *             Vzmax = 10 kNm
  30. *     - Cha2: N0 = -1600.0 kN
  31. *             Nmax = 10 kN
  32. *             Vymax = 10 kNm
  33. *             Vzmax = 200 kNm
  34. *     - Cha3: N0 = -1600.0 kN
  35. *             Nmax = 700 kN
  36. *             Vymax = 10 kNm
  37. *             Vzmax = 10 kNm    
  38. *  Pour V_TEST1 = 2, 3 type de chargement:   
  39. *     - Cha1: N0 = -1600.0 kN
  40. *             Nmax = 10 kN
  41. *             Vymax = 10 kNm
  42. *             Vzmax = 200 kNm
  43. *     - Cha3: N0 = -1600.0 kN
  44. *             Nmax = 700 kN
  45. *             Vymax = 10 kNm
  46. *             Vzmax = 10 kNm       
  47. *  Pour V_TEST1 = 3, 3 type de chargement:
  48. *     - Cha1: N0 = -1600.0 kN
  49. *             Nmax = 1 kN
  50. *             Mymax = 10 kNm
  51. *             Mzmax = 1 kNm
  52. *     - Cha2: N0 = -1600.0 kN
  53. *             Nmax = 1 kN
  54. *             Mymax = 1 kNm
  55. *             Mzmax = 10 kNm
  56. *     - Cha3: N0 = -1600.0 kN
  57. *             Nmax = 700 kN
  58. *             Mymax = 1 kNm
  59. *             Mzmax = 1 kNm  
  60. *  Pour V_TEST1 = 4, 2 type de chargement:
  61. *     - Cha1: N0 = -1600.0 kN
  62. *             Nmax = 1 kN
  63. *             Mymax = 10 kNm
  64. *             Mzmax = 1 kNm
  65. *     - Cha3: N0 = -1600.0 kN
  66. *             Nmax = 700 kN
  67. *             Mymax = 1 kNm
  68. *             Mzmax = 1 kNm  
  69. *
  70. *
  71. *
  72. *  Develloppé par Alberto FRAU /DEN/DANS/DM2S/SEMT/EMSI
  73. *              et Nicolas ILE  /DEN/DANS/DM2S/SEMT/EMSI
  74. *
  75. **********************************************************
  76. *
  77. OPTI DIME 3 MODE TRID ELEM SEG2;
  78. OPTI TRAC PSC;
  79. *
  80. * dimensions de la section 
  81. B_Y1 = 0.3;
  82. B_Z1 = 0.4;
  83. * barres longitudianles
  84. PHI_LON1 = 12.0;
  85. * enrobage
  86. ENR1 = 0.025;
  87. * Armature transversales
  88. A_CADRE1 = (2.0)*((PI)*((0.004)**(2.0)));
  89. A_CADRE2 = (2.0)*((PI)*((0.004)**(2.0)));
  90. * barres longitudianles
  91. PHI_LON1 = PROG 12.0 12.0 12.0 12.0;
  92. * Position y barres
  93. PHI_Y = PROG ((-1.0)*((B_Y1 - (2.0*ENR1))/(2.0))) 
  94.              ((+1.0)*((B_Y1 - (2.0*ENR1))/(2.0)))
  95.              ((+1.0)*((B_Y1 - (2.0*ENR1))/(2.0)))
  96.              ((-1.0)*((B_Y1 - (2.0*ENR1))/(2.0)));
  97. * Position z barres
  98. PHI_Z = PROG ((-1.0)*((B_Z1 - (2.0*ENR1))/(2.0))) 
  99.              ((-1.0)*((B_Z1 - (2.0*ENR1))/(2.0)))
  100.              ((+1.0)*((B_Z1 - (2.0*ENR1))/(2.0)))
  101.              ((+1.0)*((B_Z1 - (2.0*ENR1))/(2.0)));
  102. *
  103. *
  104. V_TEST1 = 1;
  105. *V_TEST1 = 2;
  106. *V_TEST1 = 3;
  107. *V_TEST1 = 4;
  108. *
  109. *
  110. SI (V_TEST1 EGA 1);
  111.   V_CHA1 = 1;
  112. *  V_CHA1 = 2;
  113. *  V_CHA1 = 3;
  114. FINSI;
  115. SI (V_TEST1 EGA 2);
  116.   V_CHA1 = 1;
  117. *  V_CHA1 = 3;
  118. FINSI;
  119. SI (V_TEST1 EGA 3);
  120.   V_CHA1 = 1;
  121. *  V_CHA1 = 2;
  122. *  V_CHA1 = 3;
  123. FINSI;
  124. SI (V_TEST1 EGA 4);
  125.   V_CHA1 = 1;
  126. *  V_CHA1 = 3;
  127. FINSI;
  128. *
  129. *
  130. *
  131. P1 = 0. 0. 0.;
  132. P2 = 1. 0. 0.;
  133. L1 = D 1 P1 P2;
  134. *
  135. *
  136. *------------ Cas Pouteau et Poutre courte ----------------------
  137. SI ((V_TEST1 EGA 1) OU  (V_TEST1 EGA 2));
  138. *
  139. * longuer ly lz
  140.   LTOT1Y = 1.0;
  141.   LTOT1Z = 1.5;
  142. * espacement cadres
  143.   S_CAD1 = 0.95;
  144. * calcul aire armature longitudinales
  145.   ALL1 = 0.0;
  146.   I = 1;
  147.   NN1 = DIME PHI_LON1;
  148. *
  149.   REPE IND1 (NN1);
  150.     PHI_I1 = EXTR PHI_LON1 I;
  151.     AA1 = (PI)*(((PHI_I1)*(PHI_I1))/(4.0));
  152.     AA1 = AA1/(1.E6);
  153.     ALL1 = ALL1 + AA1;
  154.     I = I + 1;
  155.   FIN IND1;
  156. * calcul omega_s et omega_w
  157.   OM_SY1 = ((ALL1)*(500.E6))/(((B_Y1)*(B_Z1))*(25.E6));
  158.   OM_SZ1 = ((ALL1)*(500.E6))/(((B_Y1)*(B_Z1))*(25.E6));
  159.   OM_WZ1 = ((A_CADRE1)*(500.E6))/(((B_Y1)*(S_CAD1))*(25.E6));
  160.   OM_WY1 = ((A_CADRE2)*(500.E6))/(((B_Z1)*(S_CAD1))*(25.E6));
  161.   LAMB_Y1 = ((LTOT1Z)/(B_Y1));
  162.   LAMB_Z1 = ((LTOT1Y)/(B_Z1));
  163.   CHI_Z1 = ((B_Z1 - ((2.0)*(ENR1)))/(B_Z1));
  164.   CHI_Y1 = ((B_Y1 - ((2.0)*(ENR1)))/(B_Y1));
  165. * test 
  166.   VAL_TT1 = LAMB_Y1 <EG (OM_SY1/OM_WY1);
  167.   VAL_TT2 = LAMB_Z1 <EG (OM_SZ1/OM_WZ1);
  168. * vymax et vzmax
  169.   TAN_Y1 = (((LAMB_Y1**2.0) + 1.0)**(0.5)) - LAMB_Y1;
  170.   TAN_Z1 = (((LAMB_Z1**2.0) + 1.0)**(0.5)) - LAMB_Z1;
  171.   VY_MAX1 = OM_WY1*CHI_Y1 + ((0.5 - OM_WY1)*(TAN_Y1));
  172.   VZ_MAX1 = OM_WZ1*CHI_Z1 + ((0.5 - OM_WZ1)*(TAN_Z1));
  173.   VY_MAX1 = ((VY_MAX1)*(((B_Y1)*(B_Z1))*(25.e6)));
  174.   VZ_MAX1 = ((VZ_MAX1)*(((B_Y1)*(B_Z1))*(25.e6)));
  175. * n traction et n compression plus nlim
  176.   N_TRAC1 = ((OM_SY1)*(((B_Y1)*(B_Z1))*(25.e6)));
  177.   N_COMP1 = ((1.0 + OM_SY1)*(((B_Y1)*(B_Z1))*(-25.e6)));
  178.   N_LIMY1 = 0.5 + OM_SY1 - ((OM_WY1)*(LAMB_Y1 + 1.0 - CHI_Y1));
  179.   N_LIMY2 = 0.5 - OM_SY1 + ((OM_WY1)*(LAMB_Y1 - 1.0 + CHI_Y1));
  180.   N_LIMZ1 = 0.5 + OM_SZ1 - ((OM_WZ1)*(LAMB_Z1 + 1.0 - CHI_Z1));
  181.   N_LIMZ2 = 0.5 - OM_SZ1 + ((OM_WZ1)*(LAMB_Z1 - 1.0 + CHI_Z1));
  182.   N_LIMY1 = ((N_LIMY1)*(((B_Y1)*(B_Z1))*(-25.e6)));
  183.   N_LIMY2 = ((N_LIMY2)*(((B_Y1)*(B_Z1))*(-25.e6)));
  184.   N_LIMZ1 = ((N_LIMZ1)*(((B_Y1)*(B_Z1))*(-25.e6)));
  185.   N_LIMZ2 = ((N_LIMZ2)*(((B_Y1)*(B_Z1))*(-25.e6)));
  186. *
  187.   VY_MAX1 = VY_MAX1/1000.;
  188.   VZ_MAX1 = VZ_MAX1/1000.;
  189.   N_TRAC1 = N_TRAC1/1000.;
  190.   N_COMP1 = N_COMP1/1000.;
  191.   N_LIMY1 = N_LIMY1/1000.;
  192.   N_LIMY2 = N_LIMY2/1000.;
  193.   N_LIMZ1 = N_LIMZ1/1000.;
  194.   N_LIMZ2 = N_LIMZ2/1000.;
  195. *
  196.   SI (V_TEST1 EGA 1);
  197.     TOL_1 = 1.E-8;
  198.     TOL_2 = 1.E-6;
  199.   FINSI;
  200.   SI (V_TEST1 EGA 1);
  201.     TOL_1 = 2.E-8;
  202.     TOL_2 = 1.E-6;
  203.   FINSI;
  204. *
  205. FINSI;
  206. *
  207. *
  208. *
  209. *
  210. *------------ Cas Pouteau et Poutre long ----------------------
  211. SI ((V_TEST1 EGA 3) OU  (V_TEST1 EGA 4));
  212. *
  213. * longuer ly lz
  214.   LTOT1Y = 2.0;
  215.   LTOT1Z = 1.5;
  216. * espacement cadres
  217.   S_CAD1 = 0.1;
  218. * calcul aire armature longitudinales
  219.   ALL1 = 0.0;
  220.   I = 1;
  221.   NN1 = DIME PHI_LON1;
  222. *
  223.   REPE IND1 (NN1);
  224.     PHI_I1 = EXTR PHI_LON1 I;
  225.     AA1 = (PI)*(((PHI_I1)*(PHI_I1))/(4.0));
  226.     AA1 = AA1/(1.E6);
  227.     ALL1 = ALL1 + AA1;
  228.     I = I + 1;
  229.   FIN IND1;
  230. * calcul omega_s et omega_w
  231.   OM_SY1 = ((ALL1)*(500.E6))/(((B_Y1)*(B_Z1))*(25.E6));
  232.   OM_SZ1 = ((ALL1)*(500.E6))/(((B_Y1)*(B_Z1))*(25.E6));
  233.   OM_WZ1 = ((A_CADRE1)*(500.E6))/(((B_Y1)*(S_CAD1))*(25.E6));
  234.   OM_WY1 = ((A_CADRE2)*(500.E6))/(((B_Z1)*(S_CAD1))*(25.E6));
  235.   LAMB_Y1 = ((LTOT1Z)/(B_Y1));
  236.   LAMB_Z1 = ((LTOT1Y)/(B_Z1));
  237.   CHI_Z1 = ((B_Z1 - ((2.0)*(ENR1)))/(B_Z1));
  238.   CHI_Y1 = ((B_Y1 - ((2.0)*(ENR1)))/(B_Y1));
  239. * test 
  240.   VAL_TT1 = LAMB_Y1 <EG (OM_SY1/OM_WY1);
  241.   VAL_TT2 = LAMB_Z1 <EG (OM_SZ1/OM_WZ1);
  242. * vymax et vzmax
  243.   TAN_Y1 = (((LAMB_Y1**2.0) + 1.0)**(0.5)) - LAMB_Y1;
  244.   TAN_Z1 = (((LAMB_Z1**2.0) + 1.0)**(0.5)) - LAMB_Z1;
  245.   VY_MAX1 = OM_WY1*CHI_Y1 + ((0.5 - OM_WY1)*(TAN_Y1));
  246.   VZ_MAX1 = OM_WZ1*CHI_Z1 + ((0.5 - OM_WZ1)*(TAN_Z1));
  247.   VY_MAX1 = ((VY_MAX1)*(((B_Y1)*(B_Z1))*(25.e6)));
  248.   VZ_MAX1 = ((VZ_MAX1)*(((B_Y1)*(B_Z1))*(25.e6)));
  249. * n traction et n compression plus nlim
  250.   N_TRAC1 = ((OM_SY1)*(((B_Y1)*(B_Z1))*(25.e6)));
  251.   N_COMP1 = ((1.0 + OM_SY1)*(((B_Y1)*(B_Z1))*(-25.e6)));
  252.   N_LIMY1 = 0.5 + OM_SY1 - ((OM_WY1)*(LAMB_Y1 + 1.0 - CHI_Y1));
  253.   N_LIMY2 = 0.5 - OM_SY1 + ((OM_WY1)*(LAMB_Y1 - 1.0 + CHI_Y1));
  254.   N_LIMZ1 = 0.5 + OM_SZ1 - ((OM_WZ1)*(LAMB_Z1 + 1.0 - CHI_Z1));
  255.   N_LIMZ2 = 0.5 - OM_SZ1 + ((OM_WZ1)*(LAMB_Z1 - 1.0 + CHI_Z1));
  256.   N_LIMY1 = ((N_LIMY1)*(((B_Y1)*(B_Z1))*(-25.e6)));
  257.   N_LIMY2 = ((N_LIMY2)*(((B_Y1)*(B_Z1))*(-25.e6)));
  258.   N_LIMZ1 = ((N_LIMZ1)*(((B_Y1)*(B_Z1))*(-25.e6)));
  259.   N_LIMZ2 = ((N_LIMZ2)*(((B_Y1)*(B_Z1))*(-25.e6)));
  260. *
  261.   VY_MAX1 = VY_MAX1/1000.;
  262.   VZ_MAX1 = VZ_MAX1/1000.;
  263.   N_TRAC1 = N_TRAC1/1000.;
  264.   N_COMP1 = N_COMP1/1000.;
  265.   N_LIMY1 = N_LIMY1/1000.;
  266.   N_LIMY2 = N_LIMY2/1000.;
  267.   N_LIMZ1 = N_LIMZ1/1000.;
  268.   N_LIMZ2 = N_LIMZ2/1000.;
  269. *N traction
  270.   N_TRAC1 =((ALL1)*(500.E6));
  271. *N compression
  272.   N_COMP1 = (((B_Y1)*(B_Z1))*(25.E6))+((ALL1)*(500.E6));
  273.   N_COMP1 = ((-1.0)*(N_COMP1));
  274. *
  275. *DEFORMATION DU BETON
  276.   EPSI_BET = 0.0035;
  277. *DEFORMATION DE L'ACIER
  278.   EPSI_ACI = 0.045;
  279. *
  280. *HAUTEUR NETTE
  281.   D_YP = B_Z1-ENR1;
  282.   D_ZP = B_Y1-ENR1;
  283.   D_YN = (B_Z1-ENR1);
  284.   D_ZN = (B_Y1-ENR1);
  285. *
  286. *position de L'AXE NEUTRE
  287.   AXEN_YP = ((EPSI_BET)/(EPSI_BET+EPSI_ACI))*D_YP;
  288.   AXEN_ZP = ((EPSI_BET)/(EPSI_BET+EPSI_ACI))*D_ZP;
  289.   AXEN_YN = ((EPSI_BET)/(EPSI_BET+EPSI_ACI))*D_YN;
  290.   AXEN_ZN = ((EPSI_BET)/(EPSI_BET+EPSI_ACI))*D_ZN;
  291. *
  292.   FF_YP = ((EPSI_BET)/(AXEN_YP));
  293.   FF_ZP = ((EPSI_BET)/(AXEN_ZP));
  294.   FF_YN = (-1.0)*((EPSI_BET)/(AXEN_YN));
  295.   FF_ZN = (-1.0)*((EPSI_BET)/(AXEN_ZN));
  296.   EPS0_YP = (FF_YP)*((B_Z1/2.0) - AXEN_YP);
  297.   EPS0_ZP = (FF_ZP)*((B_Y1/2.0) - AXEN_ZP);
  298.   EPS0_YN = (FF_YN)*(AXEN_YN - (B_Z1/2.0));
  299.   EPS0_ZN = (FF_ZN)*(AXEN_ZN - (B_Y1/2.0));
  300. *
  301. *N BALANCED
  302.   NBIL_YP = ((-0.8)*(25.E6))*((B_Y1)*(AXEN_YP));
  303.   NBIL_YN = ((-0.8)*(25.E6))*((B_Y1)*(AXEN_YN));
  304.   NBIL_ZP = ((-0.8)*(25.E6))*((B_Z1)*(AXEN_ZP));
  305.   NBIL_ZN = ((-0.8)*(25.E6))*((B_Z1)*(AXEN_ZN));
  306. *M BALANCED
  307.   MBIL_YP = ((1.0)*(NBIL_YP))*((B_Z1/2.0) - ((0.4)*(AXEN_YP)));
  308.   MBIL_ZP = ((1.0)*(NBIL_ZP))*((B_Y1/2.0) - ((0.4)*(AXEN_ZP)));
  309.   MBIL_YN = ((-1.0)*(NBIL_YN))*((B_Z1/2.0) - ((0.4)*(AXEN_YN)));
  310.   MBIL_ZN = ((-1.0)*(NBIL_ZN))*((B_Y1/2.0) - ((0.4)*(AXEN_ZN)));
  311. *
  312.   I = 1;
  313.   REPE IND1 (NN1);
  314. * Calcolo di epsi
  315.     EPSI_YP = EPS0_YP - ((FF_YP)*(EXTR PHI_Z I));
  316.     EPSI_ZP = EPS0_ZP - ((FF_ZP)*(EXTR PHI_Y I));
  317.     EPSI_YN = EPS0_YN - ((FF_YN)*(EXTR PHI_Z I));
  318.     EPSI_ZN = EPS0_ZN - ((FF_ZN)*(EXTR PHI_Y I));
  319. *calcolo di sigma
  320.     SIGI_YP = (210000.0E6)*(EPSI_YP);
  321.     SIGI_YN = (210000.0E6)*(EPSI_YN);
  322.     SIGI_ZP = (210000.0E6)*(EPSI_ZP);
  323.     SIGI_ZN = (210000.0E6)*(EPSI_ZN);
  324.     SI ((ABS(SIGI_YP)) > 500.E6);
  325.       SIGI_YP = (SIGI_YP/(ABS(SIGI_YP)))*(500.E6);
  326.     FINSI;
  327.     SI ((ABS(SIGI_YN)) > 500.E6);
  328.       SIGI_YN = (SIGI_YN/(ABS(SIGI_YN)))*(500.E6);
  329.     FINSI;
  330.     SI ((ABS(SIGI_ZP)) > 500.E6);
  331.       SIGI_ZP = (SIGI_ZP/(ABS(SIGI_ZP)))*(500.E6);
  332.     FINSI;
  333.     SI ((ABS(SIGI_ZN)) > 500.E6);
  334.       SIGI_ZN = (SIGI_ZN/(ABS(SIGI_ZN)))*(500.E6);
  335.     FINSI;
  336. *
  337.     PHI_I1 = EXTR PHI_LON1 I;
  338.     AA1 = (PI)*(((PHI_I1)*(PHI_I1))/(4.0));
  339.     AA1 = AA1/(1.E6);
  340. *
  341.     NBIL_YP = NBIL_YP + ((AA1)*(SIGI_YP));
  342.     NBIL_YN = NBIL_YN + ((AA1)*(SIGI_YN));
  343.     NBIL_ZP = NBIL_ZP + ((AA1)*(SIGI_ZP));
  344.     NBIL_ZN = NBIL_ZN + ((AA1)*(SIGI_ZN));
  345. *
  346.     MBIL_YP = MBIL_YP + (((AA1)*(SIGI_YP))*(EXTR PHI_Z I));
  347.     MBIL_YN = MBIL_YN + (((AA1)*(SIGI_YN))*(EXTR PHI_Z I));
  348.     MBIL_ZP = MBIL_ZP + (((AA1)*(SIGI_ZP))*(EXTR PHI_Y I));
  349.     MBIL_ZN = MBIL_ZN + (((AA1)*(SIGI_ZN))*(EXTR PHI_Y I));
  350. *
  351.   I = I + 1;
  352.   FIN IND1;
  353. *
  354.   MBIL_YP = MBIL_YP/1000.;
  355.   MBIL_YN = MBIL_YN/1000.;
  356.   MBIL_ZP = MBIL_ZP/1000.;
  357.   MBIL_ZN = MBIL_ZN/1000.;
  358.   NBIL_YP = NBIL_YP/1000.;
  359.   NBIL_YN = NBIL_YN/1000.;
  360.   NBIL_ZP = NBIL_ZP/1000.;
  361.   NBIL_ZN = NBIL_ZN/1000.;
  362.   N_TRAC1 = N_TRAC1/1000.;
  363.   N_COMP1 = N_COMP1/1000.;
  364. *
  365.   SI (V_TEST1 EGA 3);
  366.     TOL_1 = 1.E-8;
  367.     TOL_2 = 1.E-6;
  368.   FINSI;
  369.   SI (V_TEST1 EGA 4);
  370.     TOL_1 = 1.E-6;
  371.     TOL_2 = 1.E-6;
  372.   FINSI;
  373. *
  374. FINSI;
  375. *
  376. *
  377. *
  378. *
  379. * -------------------- Cas Pouteau courte ----------------------
  380. SI (V_TEST1 EGA 1);
  381. * Caracteristiques
  382.   CH_CAR1 = MANU CHML L1 'B_Y' B_Y1  'B_Z'  B_Z1  'LIBY' LTOT1Y
  383.                          'LIBZ' LTOT1Z 
  384.                         'SCAD' S_CAD1 'ENRB' ENR1  'PFER' PHI_LON1
  385.                         'YFER' PHI_Y   'ZFER' PHI_Z 'ASWY' A_CADRE2
  386.                         'ASWZ' A_CADRE1    'YACI' 210000.0E6  
  387.                         'EPSB' 0.0035 'EPSA' 0.045 'TYFR' 'POTEAU';
  388. *
  389. *
  390. * Chargement CAS1
  391.   SI (V_CHA1 EGA 1);
  392.     CH_SEI1 = MANU CHML L1 
  393.         'C11' (10.E3*10.E3) 'C22' (10.E3*10.E3) 'C33' (200.E3*200.E3)
  394.               'C44' (00.0*00.0) 'C55' (00.0*00.0) 'C66' (00.0*00.0)
  395.               'C12' (0.E3*0.E3) 'C13' (00.0*00.0) 'C14' (00.0*00.0)
  396.               'C15' (00.0*00.0) 'C16' (00.0*00.0);
  397.     CH_SEI1 = CH_SEI1 ET (MANU CHML L1 
  398.               'C23' (00.0*00.0) 'C24' (00.0*00.0) 'C25' (00.0*00.0)
  399.               'C26' (00.0*00.0) 'C34' (00.0*00.0) 'C35' (00.0*00.0)
  400.               'C36' (00.0*00.0) 'C45' (00.0*00.0) 'C46' (00.0*00.0)
  401.                          'C56' (00.0*00.0)) ;
  402.     CH_STA1 = MANU CHML L1 'EFFX' -1600.E3 'EFFY' 0.0 'EFFZ' 0.0
  403.                            'MOMX' 0.0 'MOMY' 0.0 'MOMZ' 0.0;
  404.   FINSI;
  405. * Chargement CAS2
  406.   SI (V_CHA1 EGA 2);
  407.     CH_SEI1 = MANU CHML L1 
  408.         'C11' (10.E3*10.E3) 'C22' (200.E3*200.E3) 'C33' (10.E3*10.E3)
  409.               'C44' (00.0*00.0) 'C55' (00.0*00.0) 'C66' (00.0*00.0)
  410.               'C12' (0.E3*0.E3) 'C13' (00.0*00.0) 'C14' (00.0*00.0)
  411.               'C15' (00.0*00.0) 'C16' (00.0*00.0);
  412.     CH_SEI1 = CH_SEI1 ET (MANU CHML L1 
  413.               'C23' (00.0*00.0) 'C24' (00.0*00.0) 'C25' (00.0*00.0)
  414.               'C26' (00.0*00.0) 'C34' (00.0*00.0) 'C35' (00.0*00.0)
  415.               'C36' (00.0*00.0) 'C45' (00.0*00.0) 'C46' (00.0*00.0)
  416.                          'C56' (00.0*00.0)) ;
  417.     CH_STA1 = MANU CHML L1 'EFFX' -1600.E3 'EFFY' 0.0 'EFFZ' 0.0
  418.                            'MOMX' 0.0 'MOMY' 0.0 'MOMZ' 0.0;
  419.   FINSI;
  420. * Chargement CAS3
  421.   SI (V_CHA1 EGA 3);
  422.     CH_SEI1 = MANU CHML L1 
  423.         'C11' (700.E3*700.E3) 'C22' (10.E3*50.E3) 'C33' (10.E3*50.E3)
  424.               'C44' (00.0*00.0) 'C55' (00.0*00.0) 'C66' (00.0*00.0)
  425.               'C12' (0.E3*0.E3) 'C13' (00.0*00.0) 'C14' (00.0*00.0)
  426.               'C15' (00.0*00.0) 'C16' (00.0*00.0);
  427.     CH_SEI1 = CH_SEI1 ET (MANU CHML L1 
  428.               'C23' (00.0*00.0) 'C24' (00.0*00.0) 'C25' (00.0*00.0)
  429.               'C26' (00.0*00.0) 'C34' (00.0*00.0) 'C35' (00.0*00.0)
  430.               'C36' (00.0*00.0) 'C45' (00.0*00.0) 'C46' (00.0*00.0)
  431.                        'C56' (00.0*00.0)) ;
  432.     CH_STA1 = MANU CHML L1 'EFFX' -1600.E3 'EFFY' 0.0 'EFFZ' 0.0
  433.                            'MOMX' 0.0 'MOMY' 0.0 'MOMZ' 0.0;
  434.   FINSI;
  435. FINSI;
  436. *
  437. *
  438. * -------------------- Cas Poutre courte ----------------------
  439. SI (V_TEST1 EGA 2);
  440. * Caracteristiques
  441.   CH_CAR1 = MANU CHML L1 'B_Y' B_Y1  'B_Z'  B_Z1  'LIBY' LTOT1Y
  442.                          'LIBZ' LTOT1Z 
  443.                        'SCAD' S_CAD1 'ENRB' ENR1  'PFER' PHI_LON1
  444.                        'YFER' PHI_Y   'ZFER' PHI_Z 'ASWY' A_CADRE2
  445.                        'ASWZ' A_CADRE1    'YACI' 210000.0E6  
  446.                        'EPSB' 0.0035 'EPSA' 0.045 'TYFR' 'POUTRE';
  447. *
  448. * Chargement CAS1
  449.   SI (V_CHA1 EGA 1);
  450.     CH_SEI1 = MANU CHML L1 
  451.         'C11' (10.E3*10.E3) 'C22' (10.E3*10.E3) 'C33' (200.E3*200.E3)
  452.               'C44' (00.0*00.0) 'C55' (00.0*00.0) 'C66' (00.0*00.0)
  453.               'C12' (0.E3*0.E3) 'C13' (00.0*00.0) 'C14' (00.0*00.0)
  454.               'C15' (00.0*00.0) 'C16' (00.0*00.0);
  455.     CH_SEI1 = CH_SEI1 ET (MANU CHML L1 
  456.               'C23' (00.0*00.0) 'C24' (00.0*00.0) 'C25' (00.0*00.0)
  457.               'C26' (00.0*00.0) 'C34' (00.0*00.0) 'C35' (00.0*00.0)
  458.               'C36' (00.0*00.0) 'C45' (00.0*00.0) 'C46' (00.0*00.0)
  459.                          'C56' (00.0*00.0)) ;
  460.     CH_STA1 = MANU CHML L1 'EFFX' -1600.E3 'EFFY' 0.0 'EFFZ' 0.0
  461.                            'MOMX' 0.0 'MOMY' 0.0 'MOMZ' 0.0;
  462.   FINSI;
  463. * Chargement CAS3
  464.   SI (V_CHA1 EGA 3);
  465.     CH_SEI1 = MANU CHML L1 
  466.         'C11' (700.E3*700.E3) 'C22' (10.E3*50.E3) 'C33' (10.E3*50.E3)
  467.               'C44' (00.0*00.0) 'C55' (00.0*00.0) 'C66' (00.0*00.0)
  468.               'C12' (0.E3*0.E3) 'C13' (00.0*00.0) 'C14' (00.0*00.0)
  469.               'C15' (00.0*00.0) 'C16' (00.0*00.0);
  470.     CH_SEI1 = CH_SEI1 ET (MANU CHML L1 
  471.               'C23' (00.0*00.0) 'C24' (00.0*00.0) 'C25' (00.0*00.0)
  472.               'C26' (00.0*00.0) 'C34' (00.0*00.0) 'C35' (00.0*00.0)
  473.               'C36' (00.0*00.0) 'C45' (00.0*00.0) 'C46' (00.0*00.0)
  474.                        'C56' (00.0*00.0)) ;
  475.     CH_STA1 = MANU CHML L1 'EFFX' -1600.E3 'EFFY' 0.0 'EFFZ' 0.0
  476.                            'MOMX' 0.0 'MOMY' 0.0 'MOMZ' 0.0;
  477.   FINSI;
  478. FINSI;
  479. *
  480. *
  481. *
  482. * -------------------- Cas Pouteau long ----------------------
  483. SI (V_TEST1 EGA 3);
  484. * Caracteristiques
  485.   CH_CAR1 = MANU CHML L1 'B_Y' B_Y1  'B_Z'  B_Z1  'LIBY' LTOT1Y
  486.                          'LIBZ' LTOT1Z 
  487.                          'SCAD' S_CAD1 'ENRB' ENR1  'PFER' PHI_LON1
  488.                          'YFER' PHI_Y   'ZFER' PHI_Z 'ASWY' A_CADRE2
  489.                          'ASWZ' A_CADRE1    'YACI' 210000.0E6  
  490.                          'EPSB' 0.0035 'EPSA' 0.045 'TYFR' 'POTEAU';
  491. *
  492. *
  493. * Chargement CAS1
  494.   SI (V_CHA1 EGA 1);
  495.     CH_SEI1 = MANU CHML L1 
  496.         'C11' (1.E3*1.E3) 'C22' (0.E3*0.E3) 'C33' (0.E3*0.E3)
  497.               'C44' (00.0*00.0) 'C55' (10.E3*10.E3) 'C66' (1.E3*1.E3)
  498.               'C12' (0.E3*0.E3) 'C13' (00.0*00.0) 'C14' (00.0*00.0)
  499.               'C15' (00.0*00.0) 'C16' (00.0*00.0);
  500.     CH_SEI1 = CH_SEI1 ET (MANU CHML L1 
  501.               'C23' (00.0*00.0) 'C24' (00.0*00.0) 'C25' (00.0*00.0)
  502.               'C26' (00.0*00.0) 'C34' (00.0*00.0) 'C35' (00.0*00.0)
  503.               'C36' (00.0*00.0) 'C45' (00.0*00.0) 'C46' (00.0*00.0)
  504.                        'C56' (00.0*00.0)) ;
  505.     CH_STA1 = MANU CHML L1 'EFFX' (NBIL_YP*1.E3) 'EFFY' 0.0 'EFFZ' 0.0
  506.                            'MOMX' 0.0 'MOMY' 0.0 'MOMZ' 0.0;
  507.   FINSI;
  508. * Chargement CAS2
  509.   SI (V_CHA1 EGA 2);
  510.     CH_SEI1 = MANU CHML L1 
  511.           'C11' (1.E3*1.E3) 'C22' (0.E3*0.E3) 'C33' (0.E3*0.E3)
  512.           'C44' (00.0*00.0) 'C55' (1.E3*1.E3) 'C66' (10.E3*10.E3)
  513.           'C12' (0.E3*0.E3) 'C13' (00.0*00.0) 'C14' (00.0*00.0)
  514.           'C15' (00.0*00.0) 'C16' (00.0*00.0);
  515.      CH_SEI1 = CH_SEI1 ET (MANU CHML L1 
  516.           'C23' (00.0*00.0) 'C24' (00.0*00.0) 'C25' (00.0*00.0)
  517.           'C26' (00.0*00.0) 'C34' (00.0*00.0) 'C35' (00.0*00.0)
  518.           'C36' (00.0*00.0) 'C45' (00.0*00.0) 'C46' (00.0*00.0)
  519.                        'C56' (00.0*00.0)) ;
  520.      CH_STA1 = MANU CHML L1 'EFFX' (NBIL_ZP*1.E3) 'EFFY' 0.0 'EFFZ' 0.0
  521.                             'MOMX' 0.0 'MOMY' 0.0 'MOMZ' 0.0;
  522.   FINSI;
  523. * Chargement CAS3
  524.   SI (V_CHA1 EGA 3);
  525.      CH_SEI1 = MANU CHML L1 
  526.           'C11' (700.E3*700.E3) 'C22' (0.E3*0.E3) 'C33' (0.E3*0.E3)
  527.           'C44' (00.0*00.0) 'C55' (1.E3*1.E3) 'C66' (1.E3*1.E3)
  528.           'C12' (0.E3*0.E3) 'C13' (00.0*00.0) 'C14' (00.0*00.0)
  529.           'C15' (00.0*00.0) 'C16' (00.0*00.0);
  530.      CH_SEI1 = CH_SEI1 ET (MANU CHML L1 
  531.           'C23' (00.0*00.0) 'C24' (00.0*00.0) 'C25' (00.0*00.0)
  532.           'C26' (00.0*00.0) 'C34' (00.0*00.0) 'C35' (00.0*00.0)
  533.           'C36' (00.0*00.0) 'C45' (00.0*00.0) 'C46' (00.0*00.0)
  534.                        'C56' (00.0*00.0)) ;
  535.      CH_STA1 = MANU CHML L1 'EFFX' (-1600*1.E3) 'EFFY' 0.0 'EFFZ' 0.0
  536.                             'MOMX' 0.0 'MOMY' 0.0 'MOMZ' 0.0;
  537.   FINSI;
  538. FINSI;
  539. *
  540. *
  541. *
  542. * -------------------- Cas Poutre long ----------------------
  543. SI (V_TEST1 EGA 4);
  544. * Caracteristiques
  545.   CH_CAR1 = MANU CHML L1 'B_Y' B_Y1  'B_Z'  B_Z1  'LIBY' LTOT1Y
  546.                          'LIBZ' LTOT1Z 
  547.                          'SCAD' S_CAD1 'ENRB' ENR1  'PFER' PHI_LON1
  548.                          'YFER' PHI_Y   'ZFER' PHI_Z 'ASWY' A_CADRE2
  549.                          'ASWZ' A_CADRE1    'YACI' 210000.0E6  
  550.                          'EPSB' 0.0035 'EPSA' 0.045 'TYFR' 'POUTRE';
  551. *
  552. *
  553. * Chargement CAS1
  554.   SI (V_CHA1 EGA 1);
  555.     CH_SEI1 = MANU CHML L1 
  556.           'C11' (1.E3*1.E3) 'C22' (0.E3*0.E3) 'C33' (0.E3*0.E3)
  557.           'C44' (00.0*00.0) 'C55' (10.E3*10.E3) 'C66' (1.E3*1.E3)
  558.           'C12' (0.E3*0.E3) 'C13' (00.0*00.0) 'C14' (00.0*00.0)
  559.           'C15' (00.0*00.0) 'C16' (00.0*00.0);
  560.     CH_SEI1 = CH_SEI1 ET (MANU CHML L1 
  561.           'C23' (00.0*00.0) 'C24' (00.0*00.0) 'C25' (00.0*00.0)
  562.           'C26' (00.0*00.0) 'C34' (00.0*00.0) 'C35' (00.0*00.0)
  563.           'C36' (00.0*00.0) 'C45' (00.0*00.0) 'C46' (00.0*00.0)
  564.                        'C56' (00.0*00.0)) ;
  565.     CH_STA1 = MANU CHML L1 'EFFX' (NBIL_YP*1.E3) 'EFFY' 0.0 'EFFZ' 0.0
  566.                            'MOMX' 0.0 'MOMY' 0.0 'MOMZ' 0.0;
  567.   FINSI;
  568. * Chargement CAS3
  569.   SI (V_CHA1 EGA 3);
  570.     CH_SEI1 = MANU CHML L1 
  571.           'C11' (700.E3*700.E3) 'C22' (0.E3*0.E3) 'C33' (0.E3*0.E3)
  572.           'C44' (00.0*00.0) 'C55' (1.E3*1.E3) 'C66' (1.E3*1.E3)
  573.           'C12' (0.E3*0.E3) 'C13' (00.0*00.0) 'C14' (00.0*00.0)
  574.           'C15' (00.0*00.0) 'C16' (00.0*00.0);
  575.     CH_SEI1 = CH_SEI1 ET (MANU CHML L1 
  576.           'C23' (00.0*00.0) 'C24' (00.0*00.0) 'C25' (00.0*00.0)
  577.           'C26' (00.0*00.0) 'C34' (00.0*00.0) 'C35' (00.0*00.0)
  578.           'C36' (00.0*00.0) 'C45' (00.0*00.0) 'C46' (00.0*00.0)
  579.                        'C56' (00.0*00.0)) ;
  580.     CH_STA1 = MANU CHML L1 'EFFX' (-1600*1.E3) 'EFFY' 0.0 'EFFZ' 0.0
  581.                            'MOMX' 0.0 'MOMY' 0.0 'MOMZ' 0.0;
  582.   FINSI;
  583. FINSI;
  584. *
  585. *
  586. *
  587. *
  588. TAB_INP                       = TABLE;
  589. TAB_INP.'MAILLAGE'            = L1;
  590. TAB_INP.'EFFORT_SEISME'       = CH_SEI1;
  591. TAB_INP.'EFFORT_STATIQUE'     = CH_STA1;
  592. TAB_INP.'CARACTERISTIQUES'    = CH_CAR1;
  593. TAB_INP.'FC_BET'              = 30.E6;
  594. TAB_INP.'GAM_C'               = 1.20;
  595. TAB_INP.'ALP_C'               = 1.00;
  596. TAB_INP.'FS_ACI'              = 500.E6;
  597. TAB_INP.'GAM_S'               = 1.0;
  598. *
  599. LLL1 = LECT 1;
  600.  
  601. TAB22 = MRCFRAME 'ELLIPSOIDE' TAB_INP 2.E-8 LLL1 ;
  602. TAB33 = MRCFRAME 'RECTANGLE' TAB_INP 1.E-6 LLL1 ;
  603. *
  604. *
  605. LAMB1 = ('EXTR' TAB22.CH_LAMBDA 'LAMB' 1 1 1);
  606. LAMB2 = ('EXTR' TAB33.CH_LAMBDA 'LAMB' 1 1 1);
  607. *
  608. *-------------------------------------------------------------------
  609. SI (V_TEST1 EGA 1);
  610.   SI (V_CHA1 EGA 1);
  611.     LAMB3 = (VZ_MAX1/200.);
  612.   FINSI;
  613.   SI (V_CHA1 EGA 2);
  614.     LAMB3 = (VY_MAX1/200.);
  615.   FINSI;
  616.   SI (V_CHA1 EGA 3);
  617.     LAMB3 = ((-1600.0 - N_COMP1)/(700.));
  618.   FINSI;
  619. FINSI;
  620. *
  621. *-------------------------------------------------------------------
  622. SI (V_TEST1 EGA 2);
  623.   SI (V_CHA1 EGA 1);
  624.     LAMB3 = (VZ_MAX1/200.);
  625.   FINSI;
  626.   SI (V_CHA1 EGA 3);
  627.     LAMB3 = ((-1600.0 - N_COMP1)/(700.));
  628.   FINSI;
  629. FINSI;
  630. *
  631. *-------------------------------------------------------------------
  632. SI (V_TEST1 EGA 3);
  633.   SI (V_CHA1 EGA 1);
  634.     LAMB3 = (MBIL_YN/10.);
  635.   FINSI;
  636.   SI (V_CHA1 EGA 2);
  637.     LAMB3 = (MBIL_ZN/10.);
  638.   FINSI;
  639.   SI (V_CHA1 EGA 3);
  640.     LAMB3 = ((-1600.0 - N_COMP1)/(700.));
  641.   FINSI;
  642. FINSI;
  643. *
  644. *-------------------------------------------------------------------
  645. SI (V_TEST1 EGA 4);
  646.   SI (V_CHA1 EGA 1);
  647.     LAMB3 = (MBIL_YN/10.);
  648.   FINSI;
  649.   SI (V_CHA1 EGA 3);
  650.     LAMB3 = ((-1600.0 - N_COMP1)/(700.));
  651.   FINSI;
  652. FINSI;
  653. *
  654. *
  655. *
  656. ERR_11 = (ABS(LAMB1 - LAMB3))/LAMB3;
  657. SI (ERR_11 > 1.E-4);
  658.   ERRE 5;
  659. FINSI;
  660. SI (LAMB2 > LAMB1);
  661.   ERRE 5;
  662. FINSI;
  663. LIST LAMB1;
  664. LIST LAMB2;
  665. LIST LAMB3;
  666. FIN;
  667.  
  668.  
  669.  
  670.  
  671.  
  672.  
  673.  
  674.  

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