Télécharger mrcshell.dgibi

Retour à la liste

Numérotation des lignes :

  1. OPTI DIME 3 MODE TRID ELEM QUA4;
  2. *
  3. **********************************************************
  4. * Test sur la procedure de MRCSHELL pour calculer
  5. * les marges de securites pour les elements en beton arme
  6. * de type coque.
  7. * On cosidere deux cas:
  8. * Cas1 CAS_1 = 1;
  9. * On considere deux elements sur lesquels on a definit
  10. * les efforts statiques et l'enveloppe des efforts sismiques
  11. * Les proprietes sont les suivantes:
  12. * epaisseur = 0.12
  13. * enrobage_ext = 0.005
  14. * enrobage_int = 0.005
  15. * ferrailage_ext = 0.001
  16. * ferrailage_int = 0.01
  17. * ferrailage_cor = 0.01
  18. * fck = 30 MPa
  19. * fsk = 500 MPa
  20. * gamma_c = 1.5
  21. * gamma_s = 1.15
  22. * nu_c = 1.0
  23. * alpa = 1.0
  24. * cot(theta) = 1.0
  25. * Pour l'element 1 on connais les marges analitiquement.
  26. * lambda*N12 = fc/2 (externe ou interne)
  27. * lambda*T1 = fs*rho_t*cot(theta) (intermediaire)
  28. * avec fc = nu*alpa*fck/gamma_c
  29. * fs = fsk/gamma_s
  30. * Pour l'element 2 l'etat statique des effort est en dehors
  31. * des surfaces ultimes. En consequence:
  32. * lambda = 0
  33. * Cas2 CAS_1 = 2;
  34. * Calcul complet des marges d'une voile qui compose une batiment
  35. * a plusieurs etages. La chene du calcul complet est effectué.
  36. * Les procedures EFFMARTI et SISSIB (avec calcul des covariances)
  37. * sont appelées
  38. *
  39. * Develloppé par Alberto FRAU /DEN/DANS/DM2S/SEMT/EMSI
  40. * et Nicolas ILE /DEN/DANS/DM2S/SEMT/EMSI
  41. *
  42. **********************************************************
  43. *
  44. CAS_1 = 1;
  45. GRAP1 = 'N';
  46. * ---------------------- Cas1 -------------------------------
  47. SI (CAS_1 EGA 1);
  48. *
  49. * Definition des points
  50. P1 = 0. 0. 0.;
  51. P2 = 1. 0. 0.;
  52. P3 = 1. 1. 0.;
  53. P4 = 0. 1. 0.;
  54. * Definition des lignes
  55. L1 = D 1 P1 P2;
  56. L2 = D 1 P2 P3;
  57. L3 = D 1 P3 P4;
  58. L4 = D 1 P4 P1;
  59. * Definition du maillage
  60. ELE1 = DALL L1 L2 L3 L4;
  61. ELE2 = ELE1 PLUS (1. 0. 0.);
  62. ELE_TOT = ELE1 ET ELE2;
  63. ELIM 0.001 ELE_TOT;
  64. * Modele
  65. MOD1 = MODE ELE1 MECANIQUE ELASTIQUE ISOTROPE COQ4;
  66. MOD2 = MODE ELE2 MECANIQUE ELASTIQUE ISOTROPE COQ4;
  67. MOD_TOT = MOD1 ET MOD2;
  68. * Champ par element covariances
  69. * [ N11e² N11eN22e N11eN12e [1.e8 0.e8 0.e8
  70. * N11eN22e N22e² N11eN12e = 0.e8 1.e8 0.e8
  71. * N11eN22e N11eN12e N12e² ] 0.e8 0.e8 25.e8]
  72. * [ N11i² N11iN22i N11iN12i [1.e8 0.e8 0.e8
  73. * N11iN22i N22i² N11iN12i = 0.e8 1.e8 0.e8
  74. * N11iN22i N11eN12i N12i² ] 0.e8 0.e8 4.e8]
  75. * [ V1² V1V2 [4.e10 0.e8
  76. * V1V2 V2²]= 0.e8 1.0e8]
  77. *
  78. CH1 = (MANU CHML MOD1
  79. 'C11E' ((1.0E4)**(2.0)) 'C22E' ((1.0E4)**(2.0)) 'C33E' ((5.0E4)**(2.0))
  80. 'C12E' ((0.0E6)**(2.0)) 'C13E' ((0.0E6)**(2.0)) 'C23E' ((0.0E6)**(2.0))
  81. 'C11I' ((1.0E4)**(2.0)) 'C22I' ((1.0E4)**(2.0)) 'C33I' ((2.0E4)**(2.0))
  82. 'C12I' ((0.0E6)**(2.0)) 'C13I' ((0.0E6)**(2.0)) 'C23I' ((0.0E6)**(2.0))
  83. 'CV11' ((2.0E5)**(2.0)) 'CV22' ((1.0E4)**(2.0))
  84. 'CV12' ((0.0E6)**(2.0)));
  85. CH1 = CH1 ET
  86. (MANU CHML MOD2
  87. 'C11E' ((1.0E4)**(2.0)) 'C22E' ((1.0E4)**(2.0)) 'C33E' ((5.0E4)**(2.0))
  88. 'C12E' ((0.0E6)**(2.0)) 'C13E' ((0.0E6)**(2.0)) 'C23E' ((0.0E6)**(2.0))
  89. 'C11I' ((1.0E4)**(2.0)) 'C22I' ((1.0E4)**(2.0)) 'C33I' ((4.0E4)**(2.0))
  90. 'C12I' ((0.0E6)**(2.0)) 'C13I' ((0.0E6)**(2.0)) 'C23I' ((0.0E6)**(2.0))
  91. 'CV11' ((2.0E5)**(2.0)) 'CV22' ((1.0E4)**(2.0))
  92. 'CV12' ((0.0E6)**(2.0)));
  93. *
  94. * Etat Statique
  95. * N11e = -10.e4
  96. * N22e = -10.e4
  97. * N12e = -0.e4
  98. * N11i = -10.e4
  99. * N22i = -10.e4
  100. * N12i = -0.e4
  101. * V1 = 0.0
  102. * V2 = 0.0
  103. CH2 = (MANU CHML MOD1
  104. 'N11E' ((-10.0E4)) 'N22E' ((-10.0E4))
  105. 'N12E' ((0.0E4)) 'N11I' ((-10.0E4))
  106. 'N22I' ((-10.0E4)) 'N12I' ((0.0E4))
  107. 'M11T' ((0.0E6)) 'M22T' ((0.0E6))
  108. 'M12T' ((0.0E6)) 'V1T' ((0.0E6))
  109. 'V2T' ((0.0E6)) 'VR' ((0.1E6)));
  110. CH2 = CH2 ET (MANU CHML MOD2
  111. 'N11E' ((-100.0E4)) 'N22E' ((-100.0E4))
  112. 'N12E' ((0.0E4)) 'N11I' ((-100.0E4))
  113. 'N22I' ((-100.0E4)) 'N12I' ((0.0E4))
  114. 'M11T' ((0.0E6)) 'M22T' ((0.0E6))
  115. 'M12T' ((0.0E6)) 'V1T' ((100.0E6))
  116. 'V2T' ((100.0E6)) 'VR' ((0.1E6)));
  117. *
  118. * Definition des caracteristiques rho1e = 0.001
  119. * rho2e = 0.01
  120. * enro_e = 0.005
  121. * rho1i = 0.001
  122. * rho2i = 0.001
  123. * enro_i = 0.005
  124. CAR_EXT1 = MANU CHML MOD_TOT 'RHO1' 0.001 'RHO2' 0.01 'ENRO' 0.005;
  125. CAR_INT1 = MANU CHML MOD_TOT 'RHO1' 0.001 'RHO2' 0.01 'ENRO' 0.005;
  126. CAR_COR1 = MANU CHML MOD_TOT 'H' 0.12 'RHOT' 0.01 'COTH' 1.0;
  127. *
  128. * DEfinition de la table d'entree
  129. TAB_MAR = TABLE;
  130. TAB_MAR.'MAILLAGE' = ELE_TOT;
  131. TAB_MAR.'EFFORT_SEISME' = CH1;
  132. TAB_MAR.'EFFORT_STATIQUE' = CH2;
  133. TAB_MAR.'CARACTERISTIQUES_EXTERNE' = CAR_EXT1;
  134. TAB_MAR.'CARACTERISTIQUES_INTERNE' = CAR_INT1;
  135. TAB_MAR.'CARACTERISTIQUES_CORE' = CAR_COR1;
  136. TAB_MAR.'FC_BET' = 30.E6;
  137. TAB_MAR.'GAM_C' = 1.5;
  138. TAB_MAR.'ALP_C' = 1.0;
  139. TAB_MAR.'NU_C' = 1.0;
  140. TAB_MAR.'FS_ACI' = 500.0E6;
  141. TAB_MAR.'GAM_S' = 1.15;
  142. *
  143. * Calcul des marges
  144. VAL_TOL1 = 0.000001;
  145. ELE_LI1 = LECT 1;
  146. OPTI EPSILON LINEAIRE;
  147. T_LAM_EL = MRCSHELL 'ELLIPSOIDE' TAB_MAR VAL_TOL1 ELE_LI1;
  148. T_LAM_RE = MRCSHELL 'RECTANGLE' TAB_MAR VAL_TOL1 ELE_LI1;
  149. *
  150. * extrapolation des resultats
  151. LAM1_E = EXTR (T_LAM_EL.CH_LAMBDA_E) 'LAME' 1 1 1;
  152. LAM1_I = EXTR (T_LAM_EL.CH_LAMBDA_I) 'LAMI' 1 1 1;
  153. LAM1_C = EXTR (T_LAM_EL.CH_LAMBDA_C) 'LAMC' 1 1 1;
  154. LAM2_E = EXTR (T_LAM_RE.CH_LAMBDA_E) 'LAME' 1 1 1;
  155. LAM2_I = EXTR (T_LAM_RE.CH_LAMBDA_I) 'LAMI' 1 1 1;
  156. LAM2_C = EXTR (T_LAM_RE.CH_LAMBDA_C) 'LAMC' 1 1 1;
  157. *
  158. LIM_E = T_LAM_EL.GRAPHIQUES. 1 .LIMITE_E;
  159. LIM_I = T_LAM_EL.GRAPHIQUES. 1 .LIMITE_I;
  160. LIM_C = T_LAM_EL.GRAPHIQUES. 1 .LIMITE_C;
  161. LIM_E = LIM_E COUL ROUG;
  162. LIM_I = LIM_I COUL ROUG;
  163. LIM_C = LIM_C COUL ROUG;
  164. *
  165. ELP_E = T_LAM_EL.GRAPHIQUES. 1 .ELLIPSOIDE_E;
  166. ELP_I = T_LAM_EL.GRAPHIQUES. 1 .ELLIPSOIDE_I;
  167. ELP_C = T_LAM_EL.GRAPHIQUES. 1 .ELLIPSOIDE_C;
  168. ELP_E = ELP_E COUL BLEU;
  169. ELP_I = ELP_I COUL BLEU;
  170. ELP_C = ELP_C COUL BLEU;
  171. *
  172. REC_E = T_LAM_EL.GRAPHIQUES. 1 .RECTANGLE_E;
  173. REC_I = T_LAM_EL.GRAPHIQUES. 1 .RECTANGLE_I;
  174. REC_C = T_LAM_EL.GRAPHIQUES. 1 .RECTANGLE_C;
  175. REC_E = REC_E COUL BLEU;
  176. REC_I = REC_I COUL BLEU;
  177. REC_C = REC_C COUL BLEU;
  178. *
  179. ELP_E_A = T_LAM_EL.GRAPHIQUES. 1 .ELLIPSOIDE_E_AUG;
  180. ELP_I_A = T_LAM_EL.GRAPHIQUES. 1 .ELLIPSOIDE_I_AUG;
  181. ELP_C_A = T_LAM_EL.GRAPHIQUES. 1 .ELLIPSOIDE_C_AUG;
  182. ELP_E_A = ELP_E_A COUL VERT;
  183. ELP_I_A = ELP_I_A COUL VERT;
  184. ELP_C_A = ELP_C_A COUL VERT;
  185. *
  186. REC_E_A = T_LAM_RE.GRAPHIQUES. 1 .RECTANGLE_E_AUG;
  187. REC_I_A = T_LAM_RE.GRAPHIQUES. 1 .RECTANGLE_I_AUG;
  188. REC_C_A = T_LAM_RE.GRAPHIQUES. 1 .RECTANGLE_C_AUG;
  189. REC_E_A = REC_E_A COUL VERT;
  190. REC_I_A = REC_I_A COUL VERT;
  191. REC_C_A = REC_C_A COUL VERT;
  192. *
  193. * Valeurs analytiques
  194. FC = 30.E6/1.5;
  195. LAM_E_T = (FC/2)/(((5.E4)/(0.01)));
  196. LAM_I_T = (FC/2)/(((2.E4)/(0.01)));
  197. FS = (500.E6/1.15)*(0.01);
  198. LAM_C_T = (FS)/(((2.E5)/(0.1)));
  199. *
  200. * Control
  201. SI ((ABS(LAM1_E - LAM_E_T)) > (1.E-8));
  202. ERRE 5;
  203. FINSI;
  204. SI ((ABS(LAM1_I - LAM_I_T)) > (1.E-8));
  205. ERRE 5;
  206. FINSI;
  207. SI ((ABS(LAM1_C - LAM_C_T)) > (1.E-8));
  208. ERRE 5;
  209. FINSI;
  210. *
  211. SI ((LAM1_E) < (LAM2_E));
  212. ERRE 5;
  213. FINSI;
  214. SI ((LAM1_I) < (LAM2_I));
  215. ERRE 5;
  216. FINSI;
  217. SI ((LAM1_C) < (LAM2_C));
  218. ERRE 5;
  219. FINSI;
  220. *
  221. * Trace surface
  222. SI ('EGA' GRAP1 'O');
  223. TRAC (1.E6 0. 0.) (LIM_E ET ELP_E ET REC_E ET ELP_E_A ET REC_E_A);
  224. TRAC (0. 1.E6 0.) (LIM_E ET ELP_E ET REC_E ET ELP_E_A ET REC_E_A);
  225. TRAC (1.E6 0. 0.) (LIM_I ET ELP_I ET REC_I ET ELP_I_A ET REC_I_A);
  226. TRAC (0. 1.E6 0.) (LIM_I ET ELP_I ET REC_I ET ELP_I_A ET REC_I_A);
  227. TRAC (0. 0. 1.E6) (LIM_C ET ELP_C ET REC_C ET ELP_C_A ET REC_C_A);
  228.  
  229. FINSI;
  230. *
  231. FIN;
  232. FINSI;
  233. SI (CAS_1 EGA 2);
  234.  
  235. *Les longueurs sont en cm
  236. *Les Forces sont en daN
  237. *Les masses sont en kg * 10-3
  238. *Les masses volumiques kg/m3 * 10-9
  239. *Les accelerations cm/sec2
  240. ******************************************************
  241. *
  242. OPTI DIME 3 ELEM QUA4;
  243. OPTI EPSI LINEAIRE ;
  244. *OPTI ECHO 0;
  245. *
  246. *
  247. *POINTS DE VUE
  248. OEIL = 3000.0 -20000.0 6000.0 ;
  249. ******************************************************
  250. * CONTRUCTION MAILLAGE
  251. ******************************************************
  252. *points � la base
  253. 0P0 = 0.0 0.0 0.0 ;
  254. 0P1 = 0P0 PLUS (700.0 0.0 0.0) ;
  255. 0P2 = 0P1 PLUS (700.0 0.0 0.0) ;
  256. 0P3 = 0P2 PLUS (0.0 700.0 0.0) ;
  257. 0P4 = 0P1 PLUS (0.0 700.0 0.0) ;
  258. 0P5 = 0P0 PLUS (0.0 700.0 0.0) ;
  259. *points haut du RDC
  260. *
  261. 1P0 = 0P0 PLUS (0.0 0.0 400.0) ;
  262. 1P1 = 0P1 PLUS (0.0 0.0 400.0) ;
  263. 1P2 = 0P2 PLUS (0.0 0.0 400.0) ;
  264. 1P3 = 0P3 PLUS (0.0 0.0 400.0) ;
  265. 1P4 = 0P4 PLUS (0.0 0.0 400.0) ;
  266. 1P5 = 0P5 PLUS (0.0 0.0 400.0) ;
  267. *
  268. *points haut du 1er etage
  269. 2P0 = 1P0 PLUS (0.0 0.0 400.0) ;
  270. 2P1 = 1P1 PLUS (0.0 0.0 400.0) ;
  271. 2P2 = 1P2 PLUS (0.0 0.0 400.0) ;
  272. 2P3 = 1P3 PLUS (0.0 0.0 400.0) ;
  273. 2P4 = 1P4 PLUS (0.0 0.0 400.0) ;
  274. 2P5 = 1P5 PLUS (0.0 0.0 400.0) ;
  275. *
  276. *points haut du 1er etage
  277. 3P0 = 2P0 PLUS (0.0 0.0 400.0) ;
  278. 3P1 = 2P1 PLUS (0.0 0.0 400.0) ;
  279. 3P2 = 2P2 PLUS (0.0 0.0 400.0) ;
  280. 3P3 = 2P3 PLUS (0.0 0.0 400.0) ;
  281. 3P4 = 2P4 PLUS (0.0 0.0 400.0) ;
  282. 3P5 = 2P5 PLUS (0.0 0.0 400.0) ;
  283. *
  284. *CONSTRUCTION DES SUPPORTS
  285. *GEOMETRIQUES POUR LES ELEMENTS
  286. *
  287. *poteaux RDC
  288. 1POT1 = (DROI 7 0P0 1P0) ;
  289. 1POT2 = (DROI 7 0P1 1P1) ;
  290. 1POT3 = (DROI 7 0P5 1P5) ;
  291. *
  292. *poteaux 1er etage
  293. 2POT1 = (DROI 7 1P0 2P0) ;
  294. 2POT2 = (DROI 7 1P1 2P1) ;
  295. 2POT3 = (DROI 7 1P5 2P5) ;
  296. *
  297. *poteaux 2eme etage
  298. 3POT1 = (DROI 7 2P0 3P0) ;
  299. 3POT2 = (DROI 7 2P1 3P1) ;
  300. 3POT3 = (DROI 7 2P5 3P5) ;
  301. *
  302. * Maillage Total tous les poteaux
  303. *
  304. SPOT = 1POT1 ET 1POT2 ET 1POT3
  305. ET 2POT1 ET 2POT2 ET 2POT3
  306. ET 3POT1 ET 3POT2 ET 3POT3 ;
  307. ELIM 0.001 SPOT ;
  308. *
  309. *TRAC OEIL CACHE SPOT ; FIN ;
  310. *poutres haut du RDC longitudinales (suivant X)
  311. 1POUL1 = (DROI 10 1P0 1P1) ;
  312. 1POUL2 = (DROI 10 1P1 1P2) ;
  313. 1POUL3 = (DROI 10 1P5 1P4) ;
  314. *
  315. *poutres haut du 1er etage longitudinales (suivant X)
  316. 2POUL1 = (DROI 10 2P0 2P1) ;
  317. 2POUL2 = (DROI 10 2P1 2P2) ;
  318. 2POUL3 = (DROI 10 2P5 2P4) ;
  319. *
  320. *poutres haut du 2eme etage longitudinales (suivant X)
  321. 3POUL1 = (DROI 10 3P0 3P1) ;
  322. 3POUL2 = (DROI 10 3P1 3P2) ;
  323. 3POUL3 = (DROI 10 3P5 3P4) ;
  324. *
  325. *Maillage pour tous les poutres longitudinales
  326. *
  327. SPOUTL = 1POUL1 ET 1POUL2 ET 1POUL3 ET
  328. 2POUL1 ET 2POUL2 ET 2POUL3 ET 3POUL1 ET 3POUL2
  329. ET 3POUL3 ;
  330. SPOUTL = SPOUTL COUL VERTE ;
  331. ELIM 0.001 SPOUTL ;
  332. *
  333. *poutres haut du RDC transversales (suivant Y)
  334. 1POUT1 = (DROI 10 1P0 1P5) ;
  335. 1POUT2 = (DROI 10 1P1 1P4) ;
  336. *
  337. *poutres haut du 1er etage transversales (suivant Y)
  338. 2POUT1 = (DROI 10 2P0 2P5) ;
  339. 2POUT2 = (DROI 10 2P1 2P4) ;
  340. *
  341. *poutres haut du 2eme etage transversales (suivant Y)
  342. 3POUT1 = (DROI 10 3P0 3P5) ;
  343. 3POUT2 = (DROI 10 3P1 3P4) ;
  344. *
  345. *
  346. * Maillage total toutes les poutres longitudinales
  347. *
  348. SPOUTT = 1POUT1 ET 1POUT2 ET 2POUT1 ET 2POUT2 ET
  349. 3POUT1 ET 3POUT2 ;
  350. SPOUTT = SPOUTT COUL JAUNE ;
  351. ELIM 0.001 SPOUTT ;
  352. *
  353. * Maillage voile suivant x
  354. *
  355. LV1 = DROI 10 0P2 0P3 ;
  356. LV2 = DROI 21 0P3 3P3 ;
  357. LV3 = DROI 10 3P3 3P2 ;
  358. LV4 = DROI 21 3P2 0P2 ;
  359. SVOILEX = DALLE LV1 LV2 LV3 LV4 'PLANE' ;
  360. ELIM 0.001 SVOILEX ;
  361. *
  362. * Maillage voile suivant y
  363. *
  364. AV1 = DROI 10 0P3 0P4 ;
  365. AV2 = DROI 21 0P4 3P4 ;
  366. AV3 = DROI 10 3P4 3P3 ;
  367. AV4 = DROI 21 3P3 0P3 ;
  368. SVOILEY = DALLE AV1 AV2 AV3 AV4 'PLANE' ;
  369. ELIM 0.001 SVOILEY ;
  370. *TRAC OEIL CACHE SVOILEY ;
  371. *
  372. SVOILEY = ORIENTER SVOILEY 'POINT' (1.0 0.0 0.0) ;
  373. SVOILEY = VERSENS SVOILEY ;
  374. *
  375. SVOILEX = ORIENTER SVOILEX 'POINT' (0.0 1.0 0.0) ;
  376. SVOILEX = VERSENS SVOILEX ;
  377. *
  378. *TRAC OEIL CACHE SVOILEY ; FIN ;
  379. *
  380. * Maillage Dalles haut du rdc et haut du 1er etage
  381. *
  382. LDA1 = DROI 20 1P0 1P2 ;
  383. LDA2 = DROI 10 1P2 1P3 ;
  384. LDA3 = DROI 20 1P3 1P5 ;
  385. LDA4 = DROI 10 1P5 1P0 ;
  386. SDAL1 = DALLE LDA1 LDA2 LDA3 LDA4 'PLANE' ;
  387. ELIM 0.001 SDAL1 ;
  388. *
  389. HDA1 = DROI 20 2P0 2P2 ;
  390. HDA2 = DROI 10 2P2 2P3 ;
  391. HDA3 = DROI 20 2P3 2P5 ;
  392. HDA4 = DROI 10 2P5 2P0 ;
  393. SDAL2 = DALLE HDA1 HDA2 HDA3 HDA4 'PLANE' ;
  394. ELIM 0.001 SDAL2 ;
  395. *
  396. SDAL3 = SDAL2 PLUS (0.0 0.0 400.0) ;
  397. *
  398. *MAILLAGE TOTAL DALLES
  399. *
  400. SDAL = SDAL1 ET SDAL2 ET SDAL3 ;
  401. *
  402. *MAILLAGE TOTAL DE LA STRUCTURE
  403. *
  404. STOT = SPOT ET SPOUTL ET SPOUTT ET SVOILEX ET
  405. SVOILEY ET SDAL ;
  406. ELIM 0.001 STOT ;
  407. SI ('EGA' GRAP1 'O');
  408. TRAC OEIL CACHE STOT ;
  409. TRAC OEIL CACHE SVOILEY ;
  410. FINSI;
  411. *
  412. ******************************************************
  413. * CONTRUCTION MODELES
  414. ******************************************************
  415. *
  416. *vecteur orientant l'axe local y des poteaux
  417. PORIEC = 1.0 0.0 0.0 ;
  418. *vecteur orientant l'axe local y des poutres longitudinales
  419. PORIEL = 0.0 1.0 0.0 ;
  420. *vecteur orientant l'axe local y des poutres transversales
  421. PORIET = 1.0 0.0 0.0 ;
  422. *
  423. *poteaux
  424. MODPOT = MODEL SPOT MECANIQUE ELASTIQUE POUT ;
  425. *poutres longitudinales
  426. MODPOUTL = MODEL SPOUTL MECANIQUE ELASTIQUE POUT ;
  427. *poutres transversales
  428. MODPOUTT = MODEL SPOUTT MECANIQUE ELASTIQUE POUT ;
  429. *voile X
  430. MODVOIX = MODEL SVOILEX MECANIQUE ELASTIQUE COQ4 ;
  431. *voile Y
  432. MODVOIY = MODEL SVOILEY MECANIQUE ELASTIQUE COQ4 ;
  433. *les deux dalles
  434. MODAL1 = MODEL SDAL1 MECANIQUE ELASTIQUE COQ4 ;
  435. MODAL2 = MODEL SDAL2 MECANIQUE ELASTIQUE COQ4 ;
  436. MODAL3 = MODEL SDAL3 MECANIQUE ELASTIQUE COQ4 ;
  437. MODAL = MODAL1 ET MODAL2 ET MODAL3 ;
  438. *
  439. MODLTOT = MODPOT ET MODPOUTL ET MODPOUTT ET MODVOIX
  440. ET MODVOIY ET MODAL;
  441. *
  442. ******************************************************
  443. * DEFINITION MATERIAUX
  444. ******************************************************
  445. *
  446. *module du b�ton
  447. EBET = 330000.0 ;
  448. *coeficient de poisson
  449. cpoi = 0.20 ;
  450. * densit� b�ton pour les poutres, les poteaux et les voiles
  451. RHOBET = 2500.0E-09 ;
  452. *
  453. *densit� b�ton dalle + charges d'exploitation dalles
  454. *(inclues dans la densit� des dalles)
  455. RHODAL = 2500.0E-09*1.67 ;
  456. *
  457. *caract�ristiques des coques et des poutres
  458. *epaisseur voiles
  459. epv = 25.0 ;
  460. *epaisseur dalles
  461. epa = 30.0 ;
  462. *excentrement des voiles et des dalles par rapport
  463. *au plan moyen
  464. ecpm = 0.0 ;
  465. *
  466. *sections poteaux
  467. bp = 50.0 ;
  468. hp = 50.0 ;
  469. *sections poutres
  470. bpout = 40.0 ;
  471. hpout = 80.0 ;
  472. *caracteristiques inertielles poteaux
  473. IIY = (bp*hp*hp*hp)/(12.0) ;
  474. IIZ = (hp*bp*bp*bp)/(12.0) ;
  475. IIX = (0.141)*(bp*hp*hp*hp) ;
  476. *caracteristiques inertielles poutres
  477. IIYP = (bpout*hpout*hpout*hpout)/(12.0) ;
  478. IIZP = (bpout*hpout*hpout*hpout)/(12.0) ;
  479. IIXP = (0.196)*(hpout*bpout*bpout*bpout) ;
  480. *section poteaux
  481. APOT = bp*hp ;
  482. *sections poutres
  483. ABEAM = bpout*hpout ;
  484. *
  485. *poteaux
  486. MATPOT = MATER MODPOT 'YOUN' EBET 'NU' cpoi 'SECT' APOT 'VECT'
  487. PORIEC 'INRY' IIY 'INRZ' IIZ 'TORS' IIX 'RHO' RHOBET ;
  488. *poutres longitudinales
  489. MATPOUTL = MATER MODPOUTL 'YOUN' EBET 'NU' cpoi 'SECT' ABEAM
  490. 'VECT' PORIEL 'INRY' IIYP 'INRZ' IIZP 'TORS' IIXP
  491. 'RHO' RHOBET ;
  492. *poutres transversales
  493. MATPOUTT = MATER MODPOUTT 'YOUN' EBET 'NU' cpoi 'SECT' ABEAM
  494. 'VECT' PORIET 'INRY' IIYP 'INRZ' IIZP 'TORS' IIXP
  495. 'RHO' RHOBET ;
  496. *voile X
  497. MATVOIX = MATER MODVOIX 'YOUN' EBET 'NU' cpoi 'EPAI' epv
  498. 'RHO' RHOBET 'EXCE' ecpm ;
  499. *voile Y
  500. MATVOIY = MATER MODVOIY 'YOUN' EBET 'NU' cpoi 'EPAI' epv
  501. 'RHO' RHOBET 'EXCE' ecpm ;
  502. *dalles
  503. MATDAL1 = MATER MODAL1 'YOUN' EBET 'NU' cpoi 'EPAI' epa
  504. 'RHO' RHODAL 'EXCE' ecpm ;
  505. *
  506. MATDAL2 = MATER MODAL2 'YOUN' EBET 'NU' cpoi 'EPAI' epa
  507. 'RHO' RHODAL 'EXCE' ecpm ;
  508. *
  509. MATDAL3 = MATER MODAL3 'YOUN' EBET 'NU' cpoi 'EPAI' epa
  510. 'RHO' RHODAL 'EXCE' ecpm ;
  511. *
  512. MATDAL = MATDAL1 ET MATDAL2 ET MATDAL3 ;
  513. *
  514. MATOT = MATPOT ET MATPOUTL ET MATPOUTT ET MATVOIX ET
  515. MATVOIY ET MATDAL ;
  516. *
  517. ******************************************************
  518. * BLOCAGES
  519. ******************************************************
  520. *
  521. CL1 = BLOQ 'DEPLA' 'ROTA' 0P0 ;
  522. CL2 = BLOQ 'DEPLA' 'ROTA' 0P1 ;
  523. CL3 = BLOQ 'DEPLA' 'ROTA' 0P5 ;
  524. CL4 = BLOQ 'DEPLA' 'ROTA' LV1 ;
  525. CL5 = BLOQ 'DEPLA' 'ROTA' AV1 ;
  526. *
  527. CLT = CL1 ET CL2 ET CL3 ET CL4 ET CL5 ;
  528. *
  529. ******************************************************
  530. * DEFINITION MATRICES RIGIDITE ET MASSE
  531. ******************************************************
  532. *
  533. RIGSTR = RIGI MODLTOT MATOT ;
  534. *rigidite totale
  535. RIGT = RIGSTR ET CLT ;
  536. MASTOT = MASS MODLTOT MATOT ;
  537. *
  538. ******************************************************
  539. * DEFINITION FORCE STATIQUE
  540. ******************************************************
  541. *
  542. CPSTAT = MANU CHPO STOT 1 UZ -981.0 ;
  543. FSTAT = MASTOT*CPSTAT ;
  544. *
  545. *
  546. ******************************************************
  547. * DEFINITION FORCE STATIQUE
  548. ******************************************************
  549. *
  550. DEPL = RESO RIGT FSTAT ;
  551. *
  552. * Contraintes statiques
  553. CS = SIGMA MODLTOT MATOT DEPL ;
  554. *
  555. *TRAC OEIL CACHE SVOILEX ; FIN ;
  556. *contraintes voile VY
  557. CSVY = REDU CS SVOILEX ;
  558. CSVY = 1.35*CSVY ;
  559. *liste CSVY ;
  560. CSVY = CHAN 'GRAVITE' MODVOIX CSVY ;
  561. *liste CSVY ; FIN ;
  562. *
  563. *********************************************
  564. *DEFINITION DU BRAS DE LEVIER
  565. *******************************
  566. couv = 2.5 ;
  567. haut = 25.0 ;
  568. d = haut-(2*couv) ;
  569. *liste d ; fin ;
  570. ttv = 0.70 ;
  571. *
  572. ******************************************************
  573. * EFFORT STATIQUE SELON LE MODELE DE MARTI
  574. ******************************************************
  575. *
  576. VECT_11 = 0. 0. 1.;
  577. VECT_22 = 0. 1. 0.;
  578. *
  579. EFF0 = EFFMARTI CSVY MODVOIX MATVOIX
  580. VECT_11 VECT_22 25. 2.5 2.5 0.7;
  581. *
  582. ******************************************************
  583. * ANALYSE MODALE
  584. ******************************************************
  585. *
  586. *
  587. * masse modales dans les 3 directions pour selection des modes utiles
  588. MMT_1 = 0. ;
  589. MMT_2 = 0. ;
  590. MMT_3 = 0. ;
  591. *
  592. * masse totale de l'ouvrage complet
  593. CHP3 = MANU CHPO STOT UX 1. UY 1. UZ 1. ;
  594. CHP3 = MASTOT * CHP3 ;
  595. *
  596. MT_1 = @TOTAL STOT CHP3 FX ;
  597. MT_2 = @TOTAL STOT CHP3 FY ;
  598. MT_3 = @TOTAL STOT CHP3 FZ ;
  599. *
  600. *
  601. *plage de fr�quences d'analyse
  602. *si l'on veut on peut aller bien sur plus loin
  603. *on peut par exemple mettre FREQSUP = 50, etc.
  604. FREQINF = 0.05 ;
  605. FREQSUP = 33.0 ;
  606. *
  607. * Analyse modale
  608. TB1 = VIBRATION 'INTERVALLE' FREQINF FREQSUP
  609. RIGT MASTOT 'MULT' TBAS ;
  610. *
  611. * nombre de modes dans la table TB1 issue de VIBR
  612. NB_MODE = (DIME (TB1 . 'MODES')) - 2 ;
  613. *
  614. * Calcul des contraintes modales
  615. TB1 = SIGSOL MODLTOT MATOT TB1 ;
  616. *
  617. *Calcul des reactions modales
  618. TB1 = REAC CLT TB1 ;
  619. *
  620. ******************************************************
  621. * CHARGEMENT SISMIQUE
  622. ******************************************************
  623. *
  624. *** liste des fr�quences [Hz]
  625. LIS_FREQ = PROG 0.1 0.25 2.5 8.0 30.0 40.0 50.0 400.0 ;
  626. *** liste des acc�l�rations pour seisme suivant x et y [m.s-2]
  627. LIS_ACCX = PROG 0.0175 0.0175 0.184 0.208 0.1 0.1 0.1 0.1 ;
  628. *** liste des acc�l�rations pour s�isme suivant z [m.s-2]
  629. LIS_ACCZ = (LIS_ACCX)*0.6666667 ;
  630. ***************************************************************
  631. LIS_ACCX = (LIS_ACCX * 981.0)*3.0 ; ;
  632. LIS_ACCZ = (LIS_ACCZ * 981.0)*3.0 ; ;
  633. *
  634. SEISME_X = EVOL MANU '[Hz]' LIS_FREQ '[cm.s-2]' LIS_ACCX ;
  635. SEISME_Y = EVOL MANU '[Hz]' LIS_FREQ '[cm.s-2]' LIS_ACCX ;
  636. SEISME_Z = EVOL ROUG MANU '[Hz]' LIS_FREQ '[cm.s-2]' LIS_ACCZ ;
  637. SI ('EGA' GRAP1 'O');
  638. DESS (SEISME_X ET SEISME_Z) ;
  639. DESS (SEISME_Y ET SEISME_Z) ;
  640. FINSI;
  641. *
  642. ******************************************************
  643. * REPONSE BASE MODALE
  644. ******************************************************
  645. *
  646. *ATENTION: ICI ON SPECIFIE L'AMORTISSEMENT
  647. *DES SPECTRES ET DES MODES
  648. *
  649. * liste des amortissements
  650. LIS_AMOR = PROG 7.e-2 ;
  651. BAS_AMOR = PROG NB_MODE * 7.e-2 ;
  652. *
  653. TAB1 = 'TABLE' ;
  654. TAB3 = 'TABLE' 'EXCITATION' ;
  655. TAB1 . 'STRUCTURE' = TB1 ;
  656. *
  657. TAB1 . 'AMORTISSEMENT' = BAS_AMOR ;
  658. TAB1 . 'EXCITATION' = TAB3 ;
  659. *TAB1 . 'TRONCATURE' = VRAI;
  660. * Excitation suivant x (axe longitudinal de l'ouvrage)
  661. TAB3 . 1 = 'TABLE' ;
  662. TAB3 . 1 . 'DIRECTION' = 'X' ;
  663. TAB3 . 1 . 'SPECTRE' = SEISME_X ;
  664. TAB3 . 1 . 'AMORTISSEMENT' = LIS_AMOR ;
  665. *
  666. * Excitation suivant y (axe transversal de l'ouvrage)
  667. TAB3 . 2 = 'TABLE' ;
  668. TAB3 . 2 . 'DIRECTION' = 'Y' ;
  669. TAB3 . 2 . 'SPECTRE' = SEISME_Y ;
  670. TAB3 . 2 . 'AMORTISSEMENT' = LIS_AMOR ;
  671. *
  672. * Excitation suivant z (axe vertical de l'ouvrage)
  673. TAB3 . 3 = 'TABLE' ;
  674. TAB3 . 3 . 'DIRECTION' = 'Z' ;
  675. TAB3 . 3 . 'SPECTRE' = SEISME_Z ;
  676. TAB3 . 3 . 'AMORTISSEMENT' = LIS_AMOR ;
  677. *
  678. * recombinaisson quadratique compl�te
  679. TAB1 . 'RECOMBINAISON_MODES' = 'CQC';
  680. * Recombinaison des directions de seisme
  681. TAB1 . 'RECOMBINAISON_DIRECTIONS' = 'QUADRATIQUE';
  682. TAB1 . 'IMPRESSION' = FAUX ;
  683. *
  684. TAB4 = 'TABLE' 'SORTIES' ;
  685. TAB1 . 'SORTIES' = TAB4 ;
  686. *
  687. * domaine sur lequel porte les resultats
  688. TAB4 . 'DOMAINE' = MODLTOT;
  689. TAB4 . 'DEPLACEMENTS' = VRAI ;
  690. TAB4 . 'REACTIONS' = VRAI ;
  691. TAB4 . 'CONTRAINTES' = VRAI ;
  692. TAB4 . 'ACCELERATIONS' = VRAI;
  693. *
  694. TAB5 = TABLE;
  695. TAB5. 1 = TABLE;
  696. TAB5. 1 .'MAIL' = SVOILEX;
  697. TAB5. 1 .'MODE' = MODVOIX;
  698. TAB5. 1 .'MATE' = MATVOIX;
  699. TAB5. 1 .'EPAISSEUR' = 25.;
  700. TAB5. 1 .'ENROBAGE_E' = 2.5;
  701. TAB5. 1 .'ENROBAGE_I' = 2.5;
  702. TAB5. 1 .'COT_THETA' = 0.7;
  703. TAB5. 1 .'DIR1' = VECT_11;
  704. TAB5. 1 .'DIR2' = VECT_22;
  705. *
  706. * Calcul SISSIB
  707. TAB2 = SISSIB TAB1 TAB5;
  708. **************************************************************
  709. *
  710. *
  711. SI ('EGA' GRAP1 'O');
  712. TRAC 'CACH' (EXCO EFF0 'N11E') MODVOIX 'TITR' 'N11E - Statique';
  713. TRAC 'CACH' (EXCO EFF0 'N22E') MODVOIX 'TITR' 'N22E - Statique';
  714. TRAC 'CACH' (EXCO EFF0 'N12E') MODVOIX 'TITR' 'N12E - Statique';
  715. TRAC 'CACH' (EXCO EFF0 'N11I') MODVOIX 'TITR' 'N11I - Statique';
  716. TRAC 'CACH' (EXCO EFF0 'N22I') MODVOIX 'TITR' 'N22I - Statique';
  717. TRAC 'CACH' (EXCO EFF0 'N12I') MODVOIX 'TITR' 'N12I - Statique';
  718. *
  719. CHEL1 = TAB2.'REPONSE_TOTALE'.'COVARIANCE'.'VOI1';
  720. TRAC 'CACH' ((EXCO CHEL1 'C11E')**(0.5)) MODVOIX 'TITR' 'N11E - Seisme';
  721. TRAC 'CACH' ((EXCO CHEL1 'C22E')**(0.5)) MODVOIX 'TITR' 'N22E - Seisme';
  722. TRAC 'CACH' ((EXCO CHEL1 'C33E')**(0.5)) MODVOIX 'TITR' 'N12E - Seisme';
  723. TRAC 'CACH' ((EXCO CHEL1 'C11I')**(0.5)) MODVOIX 'TITR' 'N11I - Seisme';
  724. TRAC 'CACH' ((EXCO CHEL1 'C22I')**(0.5)) MODVOIX 'TITR' 'N22I - Seisme';
  725. TRAC 'CACH' ((EXCO CHEL1 'C33I')**(0.5)) MODVOIX 'TITR' 'N12I - Seisme';
  726. OUBL CHEL1;
  727. FINSI;
  728. *
  729. ******************************************************
  730. * CALCUL DES MARGES
  731. ******************************************************
  732. *
  733. * On selection plusieurs zones de la voile SVOILEX
  734. CH1 CH2 CH3 = 'COOR' (SVOILEX);
  735. ELE1 = CH1 'POIN' 'COMPRIS' -2000. 2000.;
  736. CH2 = 'REDU' CH2 ELE1 ;
  737. ELE1 = CH2 'POIN' 'COMPRIS' -1000.001 1000.01;
  738. CH3 = 'REDU' CH3 ELE1 ;
  739. ELE1 = CH3 'POIN' 'COMPRIS' -0.001 400.01;
  740. CH3 = 'REDU' CH3 ELE1 ;
  741. SVOIX1 = (SVOILEX) 'ELEM' 'APPUYE'
  742. 'STRICTEMENT' ELE1;
  743. *
  744. CH1 CH2 CH3 = 'COOR' (SVOILEX);
  745. ELE1 = CH1 'POIN' 'COMPRIS' -2000. 2000.;
  746. CH2 = 'REDU' CH2 ELE1 ;
  747. ELE1 = CH2 'POIN' 'COMPRIS' -0.001 700.01;
  748. CH3 = 'REDU' CH3 ELE1 ;
  749. ELE1 = CH3 'POIN' 'COMPRIS' 399.98 800.01;
  750. CH3 = 'REDU' CH3 ELE1 ;
  751. SVOIX2 = (SVOILEX) 'ELEM' 'APPUYE'
  752. 'STRICTEMENT' ELE1;
  753. *
  754. CH1 CH2 CH3 = 'COOR' (SVOILEX);
  755. ELE1 = CH1 'POIN' 'COMPRIS' -2000. 2000.;
  756. CH2 = 'REDU' CH2 ELE1 ;
  757. ELE1 = CH2 'POIN' 'COMPRIS' -0.001 700.01;
  758. CH3 = 'REDU' CH3 ELE1 ;
  759. ELE1 = CH3 'POIN' 'COMPRIS' 799.98 1200.01;
  760. CH3 = 'REDU' CH3 ELE1 ;
  761. SVOIX3 = (SVOILEX) 'ELEM' 'APPUYE'
  762. 'STRICTEMENT' ELE1;
  763. *
  764. MODVOIX1 = MODE SVOIX1 MECANIQUE ELASTIQUE ISOTROPE COQ4;
  765. MODVOIX2 = MODE SVOIX2 MECANIQUE ELASTIQUE ISOTROPE COQ4;
  766. MODVOIX3 = MODE SVOIX3 MECANIQUE ELASTIQUE ISOTROPE COQ4;
  767. *
  768. * definition des caracteristiques
  769. CAR_EXT1 = (MANU CHML MODVOIX1 'RHO1' 0.045 'RHO2' 0.02 'ENRO' 2.5)
  770. ET (MANU CHML MODVOIX2 'RHO1' 0.035 'RHO2' 0.02 'ENRO' 2.5)
  771. ET (MANU CHML MODVOIX3 'RHO1' 0.02 'RHO2' 0.015 'ENRO' 2.5);
  772. CAR_INT1 = (MANU CHML MODVOIX1 'RHO1' 0.045 'RHO2' 0.02 'ENRO' 2.5)
  773. ET (MANU CHML MODVOIX2 'RHO1' 0.035 'RHO2' 0.02 'ENRO' 2.5)
  774. ET (MANU CHML MODVOIX3 'RHO1' 0.02 'RHO2' 0.015 'ENRO' 2.5);
  775. CAR_COR1 = (MANU CHML MODVOIX1 'H' 25. 'RHOT' 0.001 'COTH' 2.5)
  776. ET (MANU CHML MODVOIX2 'H' 25. 'RHOT' 0.001 'COTH' 2.5)
  777. ET (MANU CHML MODVOIX3 'H' 25. 'RHOT' 0.001 'COTH' 2.5);
  778. *
  779. TAB_MAR = TABLE;
  780. *
  781. FBET1 = 300.;
  782. FACI1 = 4000.;
  783. GAM_C1 = 1.0;
  784. ALP_C1 = 1.0;
  785. NU_C1 = 1.0;
  786. GAM_A1 = 1.0;
  787. * definition table d'entres
  788. TAB_MAR.'MAILLAGE' = SVOILEX;
  789. TAB_MAR.'EFFORT_SEISME' = (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1);
  790. TAB_MAR.'EFFORT_STATIQUE' = EFF0;
  791. TAB_MAR.'CARACTERISTIQUES_EXTERNE' = CAR_EXT1;
  792. TAB_MAR.'CARACTERISTIQUES_INTERNE' = CAR_INT1;
  793. TAB_MAR.'CARACTERISTIQUES_CORE' = CAR_COR1;
  794. TAB_MAR.'FC_BET' = FBET1;
  795. TAB_MAR.'GAM_C' = GAM_C1;
  796. TAB_MAR.'ALP_C' = ALP_C1;
  797. TAB_MAR.'NU_C' = NU_C1;
  798. TAB_MAR.'FS_ACI' = FACI1;
  799. TAB_MAR.'GAM_S' = GAM_A1;
  800. *
  801. * liste des elements
  802. ELE_LI1 = LECT 1 10 20 30 40 50 60 70 80 90 100;
  803. *ELE_LI1 = LECT 169;
  804. VAL_TOL1 = 0.04;
  805. * verification elipsoide
  806. T_LAM_EL = MRCSHELL 'ELLIPSOIDE' TAB_MAR VAL_TOL1 ELE_LI1;
  807. * verification rectangle
  808. T_LAM_RE = MRCSHELL 'RECTANGLE' TAB_MAR VAL_TOL1 ELE_LI1;
  809. *
  810. ******************************************************
  811. * Dessin
  812. ******************************************************
  813. * extraction des sorties du calcul
  814. LAM_RE_E = T_LAM_RE.'CH_LAMBDA_E';
  815. LAM_RE_I = T_LAM_RE.'CH_LAMBDA_I';
  816. LAM_RE_C = T_LAM_RE.'CH_LAMBDA_C';
  817. CAR_RE_E = T_LAM_RE.'CARTE_E';
  818. CAR_RE_I = T_LAM_RE.'CARTE_I';
  819. CAR_RE_C = T_LAM_RE.'CARTE_C';
  820. *
  821. LAM_EL_E = T_LAM_EL.'CH_LAMBDA_E';
  822. LAM_EL_I = T_LAM_EL.'CH_LAMBDA_I';
  823. LAM_EL_C = T_LAM_EL.'CH_LAMBDA_C';
  824. CAR_EL_E = T_LAM_EL.'CARTE_E';
  825. CAR_EL_I = T_LAM_EL.'CARTE_I';
  826. CAR_EL_C = T_LAM_EL.'CARTE_C';
  827. *
  828. * trace
  829. SI ('EGA' GRAP1 'O');
  830. V_MAX1 = MAXI LAM_RE_E;
  831. V_MIN1 = MINI LAM_RE_E;
  832. L_ISO1 = PROG V_MIN1 PAS ((V_MAX1 - V_MIN1)/(20.)) V_MAX1;
  833. TITRE 'outer layer: rectangle margin' ;
  834. TRAC 'CACH' (1.E9 0. 0.) (LAM_RE_E) MODVOIX L_ISO1;
  835. V_MAX1 = MAXI LAM_RE_I;
  836. V_MIN1 = MINI LAM_RE_I;
  837. L_ISO1 = PROG V_MIN1 PAS ((V_MAX1 - V_MIN1)/(20.)) V_MAX1;
  838. TITRE ' inner layer: rectangle margin' ;
  839. TRAC 'CACH' (1.E9 0. 0.) (LAM_RE_I) MODVOIX L_ISO1;
  840.  
  841. V_MAX1 = MAXI LAM_RE_C;
  842. V_MIN1 = MINI LAM_RE_C;
  843. L_ISO1 = PROG V_MIN1 PAS ((V_MAX1 - V_MIN1)/(20.)) V_MAX1;
  844. TITRE 'core layer: ellipsoid margin lamda' ;
  845. TRAC 'CACH' (1.E9 0. 0.) (LAM_RE_C) MODVOIX L_ISO1;
  846.  
  847. V_MAX1 = MAXI LAM_EL_E;
  848. V_MIN1 = MINI LAM_EL_E;
  849. L_ISO1 = PROG V_MIN1 PAS ((V_MAX1 - V_MIN1)/(20.)) V_MAX1;
  850. TITRE 'outer layer: ellipsoid margin lamda' ;
  851. TRAC 'CACH' (1.E9 0. 0.) (LAM_EL_E) MODVOIX L_ISO1;
  852. V_MAX1 = MAXI LAM_EL_I;
  853. V_MIN1 = MINI LAM_EL_I;
  854. L_ISO1 = PROG V_MIN1 PAS ((V_MAX1 - V_MIN1)/(20.)) V_MAX1;
  855. TITRE 'inner layer: ellipsoid margin lamda' ;
  856. TRAC 'CACH' (1.E9 0. 0.) (LAM_EL_I) MODVOIX L_ISO1;
  857.  
  858. V_MAX1 = MAXI LAM_EL_C;
  859. V_MIN1 = MINI LAM_EL_C;
  860. L_ISO1 = PROG V_MIN1 PAS ((V_MAX1 - V_MIN1)/(20.)) V_MAX1;
  861. TITRE 'core layer: ellipsoid margin lamda' ;
  862. TRAC 'CACH' (1.E9 0. 0.) (LAM_EL_C) MODVOIX L_ISO1;
  863.  
  864. *
  865. TRAC 'CACH' 'FACE' (1.E9 0. 0.)
  866. ((CAR_RE_E) ET (CAR_EL_E PLUS (0. 800. 0.)))
  867. 'TITR' 'outer layer: left RECTANGLE right ELLIPSOID
  868. RED margin < 1, GREEN margin > 1';
  869. TRAC 'CACH' 'FACE' (1.E9 0. 0.)
  870. ((CAR_RE_I) ET (CAR_EL_I PLUS (0. 800. 0.)))
  871. 'TITR' 'inner layer: left RECTANGLE right ELLIPSOID
  872. RED margin < 1, GREEN margin > 1';
  873. TRAC 'CACH' 'FACE' (1.E9 0. 0.)
  874. ((CAR_RE_C) ET (CAR_EL_C PLUS (0. 800. 0.)))
  875. 'TITR' 'core layer: left RECTANGLE right ELLIPSOID
  876. RED margin < 1, GREEN margin > 1';
  877. FINSI;
  878. FIN;
  879. FINSI;
  880.  
  881.  
  882.  
  883.  
  884.  

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