Télécharger mrcframe_test.dgibi

Retour à la liste

Numérotation des lignes :

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

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