Télécharger sissib_cov.dgibi

Retour à la liste

Numérotation des lignes :

  1. OPTI DIME 3 MODE TRID ELEM QUA4;
  2. *
  3. **********************************************************
  4. * Test sur la procedure SISSIB, suite à la modification
  5. * pour le calcul des efforts maximaux et des covariances sur
  6. * base modale des efforts generalisée pour les coques du modele
  7. * a trois couches de MARTI
  8. *
  9. * On considere une plaque de 15m de larguer et 1m de profondeur
  10. * soumise à trois mouvements en direction x, y et z (reppresentés
  11. * par les spectres de reponse)
  12. * La plaque est encastrée sur un coté et libre sur les autres.
  13. * Les termes de masse sont presents seulement sur un coté.
  14. * Le comportment dynamique est reppresenté par deux modes de
  15. * vibration.
  16. *
  17. * Develloppé par Alberto FRAU /DEN/DANS/DM2S/SEMT/EMSI
  18. * et Nicolas ILE /DEN/DANS/DM2S/SEMT/EMSI
  19. *
  20. **********************************************************
  21. *
  22. * option graphiques
  23. GRAP1 = 0;
  24. *
  25. * Proprietés geometriques
  26. H1 = 0.1;
  27. ER1 = 0.025;
  28. ER2 = 0.025;
  29. COT1 = 1.43;
  30. *
  31. * Definition du chargement sismique x, y et z
  32. LIS_FREQ = PROG 0.1 0.25 2.5 8.0 30.0 40.0 50.0 400.0 ;
  33. * liste des acc?l?rations pour seisme suivant x et y [m.s-2]
  34. LIS_ACCX = PROG 0.0175 0.0175 0.184 0.208 0.1 0.1 0.1 0.1 ;
  35. * liste des acc?l?rations pour s?isme suivant z [m.s-2]
  36. LIS_ACCZ = (LIS_ACCX)*0.6666667 ;
  37. *
  38. LIS_ACCX = (LIS_ACCX * 9.810)*3.0 ;
  39. LIS_ACCZ = (LIS_ACCZ * 9.810)*3.0 ;
  40. *
  41. SEISME_X = EVOL MANU '[Hz]' LIS_FREQ '[m.s-2]' LIS_ACCX ;
  42. SEISME_Y = EVOL MANU '[Hz]' LIS_FREQ '[m.s-2]' LIS_ACCX ;
  43. SEISME_Z = EVOL ROUG MANU '[Hz]' LIS_FREQ '[m.s-2]' LIS_ACCZ ;
  44. *
  45. * -------------------------------------------------------------
  46. P1 = 0. 0. 0.;
  47. P2 = 15. 0. 0.;
  48. P3 = 15. 1. 0.;
  49. P4 = 0. 1. 0.;
  50. *
  51. L1 = D 15 P1 P2;
  52. L2 = D 3 P2 P3;
  53. L3 = D 15 P3 P4;
  54. L3 = D 15 P3 P4;
  55. L4 = D 3 P4 P1;
  56. *
  57. PT1 = 3. 0. 0.;
  58. PT2 = 3. 1. 0.;
  59. LT1 = D 1 PT1 PT2;
  60. *
  61. * definition maillage
  62. SUR1 = DALL L1 L2 L3 L4;
  63. ELIM (SUR1 ET LT1) 0.0001;
  64. *
  65. * modele et materiau
  66. MOD1 = MODE SUR1 MECANIQUE ELASTIQUE ISOTROPE COQ4;
  67. MAT1 = MATE MOD1 YOUNG 30000.E6 NU 0.2 RHO 0. EPAI H1;
  68. *
  69. * matrices de masse et rigidite
  70. RIG1 = RIGI MOD1 MAT1;
  71. MAS1 = (MASS 'UX' L2 80000.);
  72. *
  73. * blocages
  74. BL1 = (BLOQUER 'DEPLA' 'ROTA' L4);
  75. * (BLOQUER 'UY' L3) ET (BLOQUER 'UY' L1);
  76. *
  77. * definition de la base modale
  78. OPTI EPSI LINEAIRE;
  79. TB1 = VIBR INTERVALLE 0. 100. BASSE 2
  80. (RIG1 ET BL1) (MAS1) IMPR ;
  81. *
  82. * calcul contraintes modales
  83. TB1 = SIGSOL MOD1 MAT1 TB1 ;
  84. SIG_M1 = (TB1.MODES. 1 .CONTRAINTE_MODALE);
  85. SIG_M2 = (TB1.MODES. 2 .CONTRAINTE_MODALE);
  86. SIG_M1 = RTENS SIG_M1 MOD1 MAT1 (1.0 0. 0.) (0. 1. 0.);
  87. SIG_M2 = RTENS SIG_M2 MOD1 MAT1 (1.0 0. 0.) (0. 1. 0.);
  88. SIG_M1 = 'CHAN' 'GRAVITE' SIG_M1 MOD1;
  89. SIG_M2 = 'CHAN' 'GRAVITE' SIG_M2 MOD1;
  90. SIG_M1 = REDU SIG_M1 (SUR1 ELEM 1);
  91. SIG_M2 = REDU SIG_M2 (SUR1 ELEM 1);
  92. *
  93. * extraction des efforts globaux pour l'element 1
  94. * mode 1
  95. N11_M1 = EXTR SIG_M1 'N11' 1 1 1;
  96. N22_M1 = EXTR SIG_M1 'N22' 1 1 1;
  97. N12_M1 = EXTR SIG_M1 'N12' 1 1 1;
  98. M11_M1 = EXTR SIG_M1 'M11' 1 1 1;
  99. M22_M1 = EXTR SIG_M1 'M22' 1 1 1;
  100. M12_M1 = EXTR SIG_M1 'M12' 1 1 1;
  101. V1_M1 = EXTR SIG_M1 'V1' 1 1 1;
  102. V2_M1 = EXTR SIG_M1 'V2' 1 1 1;
  103. *
  104. * extraction des efforts globaux pour l'element 1
  105. * mode 2
  106. N11_M2 = EXTR SIG_M2 'N11' 1 1 1;
  107. N22_M2 = EXTR SIG_M2 'N22' 1 1 1;
  108. N12_M2 = EXTR SIG_M2 'N12' 1 1 1;
  109. M11_M2 = EXTR SIG_M2 'M11' 1 1 1;
  110. M22_M2 = EXTR SIG_M2 'M22' 1 1 1;
  111. M12_M2 = EXTR SIG_M2 'M12' 1 1 1;
  112. V1_M2 = EXTR SIG_M2 'V1' 1 1 1;
  113. V2_M2 = EXTR SIG_M2 'V2' 1 1 1;
  114. *
  115. * Calcul des reactions modales
  116. TB1 = REAC BL1 TB1 ;
  117. *
  118. SI (GRAP1 EGA 1);
  119. NN1 = (DIME TB1.MODES) - 2;
  120. I = 1;
  121. REPE IND1 (NN1);
  122. CH1 = TB1.MODES. I .DEFORMEE_MODALE;
  123. FREQ1 = TB1.MODES. I .FREQUENCE;
  124. TIT1 = CHAINE 'Mode ' I ' - Frequence ' FREQ1;
  125. DEF0 = DEFO CH1 SUR1 0. ROUG;
  126. DEF1 = DEFO CH1 SUR1 VERT;
  127. TRAC (0. 0. 1.E6) (DEF1 ET DEF0) TITR TIT1;
  128. I = I + 1;
  129. FIN IND1;
  130. FINSI;
  131. *
  132. * Repose base modale
  133. *
  134. NB_MODE = (DIME TB1.MODES) - 2;
  135. * liste des amortissements
  136. LIS_AMOR = PROG 7.e-2 ;
  137. BAS_AMOR = PROG NB_MODE * 7.e-2 ;
  138. *
  139. TAB1 = 'TABLE' ;
  140. TAB3 = 'TABLE' 'EXCITATION' ;
  141. TAB1 . 'STRUCTURE' = TB1 ;
  142. *
  143. TAB1 . 'AMORTISSEMENT' = BAS_AMOR ;
  144. TAB1 . 'EXCITATION' = TAB3 ;
  145. *TAB1 . 'TRONCATURE' = VRAI;
  146. * Excitation suivant x
  147. TAB3 . 1 = 'TABLE' ;
  148. TAB3 . 1 . 'DIRECTION' = 'X' ;
  149. TAB3 . 1 . 'SPECTRE' = SEISME_X ;
  150. TAB3 . 1 . 'AMORTISSEMENT' = LIS_AMOR ;
  151. *
  152. * Excitation suivant y
  153. TAB3 . 2 = 'TABLE' ;
  154. TAB3 . 2 . 'DIRECTION' = 'Y' ;
  155. TAB3 . 2 . 'SPECTRE' = (SEISME_Y) ;
  156. TAB3 . 2 . 'AMORTISSEMENT' = LIS_AMOR ;
  157. *
  158. * Excitation suivant z
  159. TAB3 . 3 = 'TABLE' ;
  160. TAB3 . 3 . 'DIRECTION' = 'Z' ;
  161. TAB3 . 3 . 'SPECTRE' = (SEISME_Z);
  162. TAB3 . 3 . 'AMORTISSEMENT' = LIS_AMOR ;
  163. *
  164. * recombinaisson quadratique complete
  165. TAB1 . 'RECOMBINAISON_MODES' = 'CQC';
  166. * Recombinaison des directions de seisme
  167. TAB1 . 'RECOMBINAISON_DIRECTIONS' = 'QUADRATIQUE';
  168. TAB1 . 'IMPRESSION' = FAUX ;
  169. *
  170. TAB4 = 'TABLE' 'SORTIES' ;
  171. TAB1 . 'SORTIES' = TAB4 ;
  172. *
  173. * domaine sur lequel porte les resultats
  174. TAB4 . 'DOMAINE' = MOD1;
  175. TAB4 . 'DEPLACEMENTS' = VRAI ;
  176. TAB4 . 'REACTIONS' = VRAI ;
  177. TAB4 . 'CONTRAINTES' = VRAI ;
  178. TAB4 . 'ACCELERATIONS' = VRAI;
  179. *
  180. TAB5 = TABLE;
  181. TAB5. 1 = TABLE;
  182. TAB5. 1 .'MAIL' = SUR1;
  183. TAB5. 1 .'MODE' = MOD1;
  184. TAB5. 1 .'MATE' = MAT1;
  185. TAB5. 1 .'EPAISSEUR' = H1;
  186. TAB5. 1 .'ENROBAGE_E' = ER1;
  187. TAB5. 1 .'ENROBAGE_I' = ER2;
  188. TAB5. 1 .'COT_THETA' = COT1;
  189. TAB5. 1 .'DIR1' = (1.0 0. 0.);
  190. TAB5. 1 .'DIR2' = (0.0 1. 0.);
  191. *
  192. * Calcul SISSIB
  193. TAB2 = SISSIB TAB1 TAB5;
  194. *
  195. * Determination analique des efforts selon
  196. * le modele de Marti pour l'element 1 pour
  197. * chaque champ des efforts modaux
  198. N11_M1E = N11_M1/2;
  199. N22_M1E = N22_M1/2;
  200. N12_M1E = N12_M1/2;
  201. N11_M1I = N11_M1/2;
  202. N22_M1I = N22_M1/2;
  203. N12_M1I = N12_M1/2;
  204. N11_M2E = N11_M2/2;
  205. N22_M2E = N22_M2/2;
  206. N12_M2E = N12_M2/2;
  207. N11_M2I = N11_M2/2;
  208. N22_M2I = N22_M2/2;
  209. N12_M2I = N12_M2/2;
  210. *
  211. * On extrait les efforts selon MARTI deduits par
  212. * la procedure pour chaque champ des efforts modaux - Mode 1
  213. N11_M1ET = EXTR
  214. (REDU (TB1.MODES. 1 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  215. 'N11E' 1 1 1;
  216. N22_M1ET = EXTR
  217. (REDU (TB1.MODES. 1 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  218. 'N22E' 1 1 1;
  219. N12_M1ET = EXTR
  220. (REDU (TB1.MODES. 1 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  221. 'N12E' 1 1 1;
  222. N11_M1IT = EXTR
  223. (REDU (TB1.MODES. 1 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  224. 'N11I' 1 1 1;
  225. N22_M1IT = EXTR
  226. (REDU (TB1.MODES. 1 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  227. 'N22I' 1 1 1;
  228. N12_M1IT = EXTR
  229. (REDU (TB1.MODES. 1 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  230. 'N12I' 1 1 1;
  231. * On extrait les efforts selon MARTI deduits par
  232. * la procedure pour chaque champ des efforts modaux - Mode 2
  233. N11_M2ET = EXTR
  234. (REDU (TB1.MODES. 2 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  235. 'N11E' 1 1 1;
  236. N22_M2ET = EXTR
  237. (REDU (TB1.MODES. 2 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  238. 'N22E' 1 1 1;
  239. N12_M2ET = EXTR
  240. (REDU (TB1.MODES. 2 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  241. 'N12E' 1 1 1;
  242. N11_M2IT = EXTR
  243. (REDU (TB1.MODES. 2 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  244. 'N11I' 1 1 1;
  245. N22_M2IT = EXTR
  246. (REDU (TB1.MODES. 2 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  247. 'N22I' 1 1 1;
  248. N12_M2IT = EXTR
  249. (REDU (TB1.MODES. 2 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  250. 'N12I' 1 1 1;
  251. *
  252. * Control
  253. SI ((ABS(N11_M1ET - N11_M1E)) > 1.E-12);
  254. ERRE 5;
  255. FINSI;
  256. SI ((ABS(N22_M1ET - N22_M1E)) > 1.E-12);
  257. ERRE 5;
  258. FINSI;
  259. SI ((ABS(N12_M1ET - N12_M1E)) > 1.E-12);
  260. ERRE 5;
  261. FINSI;
  262. SI ((ABS(N11_M1IT - N11_M1I)) > 1.E-12);
  263. ERRE 5;
  264. FINSI;
  265. SI ((ABS(N22_M1IT - N22_M1I)) > 1.E-12);
  266. ERRE 5;
  267. FINSI;
  268. SI ((ABS(N12_M1IT - N12_M1I)) > 1.E-12);
  269. ERRE 5;
  270. FINSI;
  271. *
  272. SI ((ABS(N11_M2ET - N11_M2E)) > 1.E-12);
  273. ERRE 5;
  274. FINSI;
  275. SI ((ABS(N22_M2ET - N22_M2E)) > 1.E-12);
  276. ERRE 5;
  277. FINSI;
  278. SI ((ABS(N12_M2ET - N12_M2E)) > 1.E-12);
  279. ERRE 5;
  280. FINSI;
  281. SI ((ABS(N11_M2IT - N11_M2I)) > 1.E-12);
  282. ERRE 5;
  283. FINSI;
  284. SI ((ABS(N22_M2IT - N22_M2I)) > 1.E-12);
  285. ERRE 5;
  286. FINSI;
  287. SI ((ABS(N12_M2IT - N12_M2I)) > 1.E-12);
  288. ERRE 5;
  289. FINSI;
  290. *
  291. * Determination de la Pseudo acceleration pour les Modes 1 et 2
  292. * pour les direction x, y et z
  293. PSE_A1X = TAB2.'REPRISE'.'SPECTRE'. 1 . 1;
  294. PSE_A1Y = TAB2.'REPRISE'.'SPECTRE'. 1 . 2;
  295. PSE_A1Z = TAB2.'REPRISE'.'SPECTRE'. 1 . 3;
  296. PSE_A2X = TAB2.'REPRISE'.'SPECTRE'. 2 . 1;
  297. PSE_A2Y = TAB2.'REPRISE'.'SPECTRE'. 2 . 2;
  298. PSE_A2Z = TAB2.'REPRISE'.'SPECTRE'. 2 . 3;
  299. *
  300. * Determination de l'amplification pour les Modes 1 et 2
  301. * pour les direction x, y et z
  302. AMP_1X = (TB1.MODES. 1 .FACT_PART_NORM. 1)*(PSE_A1X);
  303. AMP_1Y = (TB1.MODES. 1 .FACT_PART_NORM. 2)*(PSE_A1Y);
  304. AMP_1Z = (TB1.MODES. 1 .FACT_PART_NORM. 3)*(PSE_A1Z);
  305. AMP_2X = (TB1.MODES. 2 .FACT_PART_NORM. 1)*(PSE_A2X);
  306. AMP_2Y = (TB1.MODES. 2 .FACT_PART_NORM. 2)*(PSE_A2Y);
  307. AMP_2Z = (TB1.MODES. 2 .FACT_PART_NORM. 3)*(PSE_A2Z);
  308. *
  309. * Efforts projetés amplifiés (mode 1 ; directions x, y, z)
  310. N11_M1EX = N11_M1E*AMP_1X;
  311. N22_M1EX = N22_M1E*AMP_1X;
  312. N12_M1EX = N12_M1E*AMP_1X;
  313. N11_M1IX = N11_M1I*AMP_1X;
  314. N22_M1IX = N22_M1I*AMP_1X;
  315. N12_M1IX = N12_M1I*AMP_1X;
  316. N11_M1EY = N11_M1E*AMP_1Y;
  317. N22_M1EY = N22_M1E*AMP_1Y;
  318. N12_M1EY = N12_M1E*AMP_1Y;
  319. N11_M1IY = N11_M1I*AMP_1Y;
  320. N22_M1IY = N22_M1I*AMP_1Y;
  321. N12_M1IY = N12_M1I*AMP_1Y;
  322. N11_M1EZ = N11_M1E*AMP_1Z;
  323. N22_M1EZ = N22_M1E*AMP_1Z;
  324. N12_M1EZ = N12_M1E*AMP_1Z;
  325. N11_M1IZ = N11_M1I*AMP_1Z;
  326. N22_M1IZ = N22_M1I*AMP_1Z;
  327. N12_M1IZ = N12_M1I*AMP_1Z;
  328. * Efforts projetés amplifiés (mode 2 ; directions x, y, z)
  329. N11_M2EX = N11_M2E*AMP_2X;
  330. N22_M2EX = N22_M2E*AMP_2X;
  331. N12_M2EX = N12_M2E*AMP_2X;
  332. N11_M2IX = N11_M2I*AMP_2X;
  333. N22_M2IX = N22_M2I*AMP_2X;
  334. N12_M2IX = N12_M2I*AMP_2X;
  335. N11_M2EY = N11_M2E*AMP_2Y;
  336. N22_M2EY = N22_M2E*AMP_2Y;
  337. N12_M2EY = N12_M2E*AMP_2Y;
  338. N11_M2IY = N11_M2I*AMP_2Y;
  339. N22_M2IY = N22_M2I*AMP_2Y;
  340. N12_M2IY = N12_M2I*AMP_2Y;
  341. N11_M2EZ = N11_M2E*AMP_2Z;
  342. N22_M2EZ = N22_M2E*AMP_2Z;
  343. N12_M2EZ = N12_M2E*AMP_2Z;
  344. N11_M2IZ = N11_M2I*AMP_2Z;
  345. N22_M2IZ = N22_M2I*AMP_2Z;
  346. N12_M2IZ = N12_M2I*AMP_2Z;
  347. *
  348. * determination du terme N11N11 couche externe
  349. C1111E_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1EX)*(N11_M1EX)));
  350. C1111E_X = C1111E_X +
  351. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2EX)*(N11_M1EX)));
  352. C1111E_X = C1111E_X +
  353. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1EX)*(N11_M2EX)));
  354. C1111E_X = C1111E_X +
  355. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2EX)*(N11_M2EX)));
  356. C1111E_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1EY)*(N11_M1EY)));
  357. C1111E_Y = C1111E_Y +
  358. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2EY)*(N11_M1EY)));
  359. C1111E_Y = C1111E_Y +
  360. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1EY)*(N11_M2EY)));
  361. C1111E_Y = C1111E_Y +
  362. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2EY)*(N11_M2EY)));
  363. C1111E_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1EZ)*(N11_M1EZ)));
  364. C1111E_Z = C1111E_Z +
  365. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2EZ)*(N11_M1EZ)));
  366. C1111E_Z = C1111E_Z +
  367. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1EZ)*(N11_M2EZ)));
  368. C1111E_Z = C1111E_Z +
  369. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2EZ)*(N11_M2EZ)));
  370. C1111E = C1111E_X + C1111E_Y + C1111E_Z;
  371. C1111EA = EXTR
  372. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  373. 'C11E' 1 1 1;
  374. * determination du terme N22N22 couche externe
  375. C2222E_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1EX)*(N22_M1EX)));
  376. C2222E_X = C2222E_X +
  377. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2EX)*(N22_M1EX)));
  378. C2222E_X = C2222E_X +
  379. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1EX)*(N22_M2EX)));
  380. C2222E_X = C2222E_X +
  381. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2EX)*(N22_M2EX)));
  382. C2222E_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1EY)*(N22_M1EY)));
  383. C2222E_Y = C2222E_Y +
  384. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2EY)*(N22_M1EY)));
  385. C2222E_Y = C2222E_Y +
  386. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1EY)*(N22_M2EY)));
  387. C2222E_Y = C2222E_Y +
  388. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2EY)*(N22_M2EY)));
  389. C2222E_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1EZ)*(N22_M1EZ)));
  390. C2222E_Z = C2222E_Z +
  391. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2EZ)*(N22_M1EZ)));
  392. C2222E_Z = C2222E_Z +
  393. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1EZ)*(N22_M2EZ)));
  394. C2222E_Z = C2222E_Z +
  395. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2EZ)*(N22_M2EZ)));
  396. C2222E = C2222E_X + C2222E_Y + C2222E_Z;
  397. C2222EA = EXTR
  398. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  399. 'C22E' 1 1 1;
  400. * determination du terme N12N12 couche externe
  401. C1212E_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N12_M1EX)*(N12_M1EX)));
  402. C1212E_X = C1212E_X +
  403. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M2EX)*(N12_M1EX)));
  404. C1212E_X = C1212E_X +
  405. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M1EX)*(N12_M2EX)));
  406. C1212E_X = C1212E_X +
  407. ((TAB2.REPRISE.PIJ. 2 . 2)*((N12_M2EX)*(N12_M2EX)));
  408. C1212E_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N12_M1EY)*(N12_M1EY)));
  409. C1212E_Y = C1212E_Y +
  410. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M2EY)*(N12_M1EY)));
  411. C1212E_Y = C1212E_Y +
  412. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M1EY)*(N12_M2EY)));
  413. C1212E_Y = C1212E_Y +
  414. ((TAB2.REPRISE.PIJ. 2 . 2)*((N12_M2EY)*(N12_M2EY)));
  415. C1212E_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N12_M1EZ)*(N12_M1EZ)));
  416. C1212E_Z = C1212E_Z +
  417. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M2EZ)*(N12_M1EZ)));
  418. C1212E_Z = C1212E_Z +
  419. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M1EZ)*(N12_M2EZ)));
  420. C1212E_Z = C1212E_Z +
  421. ((TAB2.REPRISE.PIJ. 2 . 2)*((N12_M2EZ)*(N12_M2EZ)));
  422. C1212E = C1212E_X + C1212E_Y + C1212E_Z;
  423. C1212EA = EXTR
  424. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  425. 'C33E' 1 1 1;
  426. * determination du terme N11N22 couche externe
  427. C1122E_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1EX)*(N22_M1EX)));
  428. C1122E_X = C1122E_X +
  429. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2EX)*(N22_M1EX)));
  430. C1122E_X = C1122E_X +
  431. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1EX)*(N22_M2EX)));
  432. C1122E_X = C1122E_X +
  433. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2EX)*(N22_M2EX)));
  434. C1122E_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1EY)*(N22_M1EY)));
  435. C1122E_Y = C1122E_Y +
  436. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2EY)*(N22_M1EY)));
  437. C1122E_Y = C1122E_Y +
  438. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1EY)*(N22_M2EY)));
  439. C1122E_Y = C1122E_Y +
  440. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2EY)*(N22_M2EY)));
  441. C1122E_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1EZ)*(N22_M1EZ)));
  442. C1122E_Z = C1122E_Z +
  443. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2EZ)*(N22_M1EZ)));
  444. C1122E_Z = C1122E_Z +
  445. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1EZ)*(N22_M2EZ)));
  446. C1122E_Z = C1122E_Z +
  447. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2EZ)*(N22_M2EZ)));
  448. C1122E = C1122E_X + C1122E_Y + C1122E_Z;
  449. C1122EA = EXTR
  450. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  451. 'C12E' 1 1 1;
  452. * determination du terme N11N12 couche externe
  453. C1112E_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1EX)*(N12_M1EX)));
  454. C1112E_X = C1112E_X +
  455. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2EX)*(N12_M1EX)));
  456. C1112E_X = C1112E_X +
  457. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1EX)*(N12_M2EX)));
  458. C1112E_X = C1112E_X +
  459. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2EX)*(N12_M2EX)));
  460. C1112E_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1EY)*(N12_M1EY)));
  461. C1112E_Y = C1112E_Y +
  462. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2EY)*(N12_M1EY)));
  463. C1112E_Y = C1112E_Y +
  464. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1EY)*(N12_M2EY)));
  465. C1112E_Y = C1112E_Y +
  466. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2EY)*(N12_M2EY)));
  467. C1112E_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1EZ)*(N12_M1EZ)));
  468. C1112E_Z = C1112E_Z +
  469. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2EZ)*(N12_M1EZ)));
  470. C1112E_Z = C1112E_Z +
  471. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1EZ)*(N12_M2EZ)));
  472. C1112E_Z = C1112E_Z +
  473. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2EZ)*(N12_M2EZ)));
  474. C1112E = C1112E_X + C1112E_Y + C1112E_Z;
  475. C1112EA = EXTR
  476. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  477. 'C13E' 1 1 1;
  478. * determination du terme N22N12 couche externe
  479. C2212E_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1EX)*(N12_M1EX)));
  480. C2212E_X = C2212E_X +
  481. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2EX)*(N12_M1EX)));
  482. C2212E_X = C2212E_X +
  483. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1EX)*(N12_M2EX)));
  484. C2212E_X = C2212E_X +
  485. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2EX)*(N12_M2EX)));
  486. C2212E_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1EY)*(N12_M1EY)));
  487. C2212E_Y = C2212E_Y +
  488. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2EY)*(N12_M1EY)));
  489. C2212E_Y = C2212E_Y +
  490. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1EY)*(N12_M2EY)));
  491. C2212E_Y = C2212E_Y +
  492. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2EY)*(N12_M2EY)));
  493. C2212E_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1EZ)*(N12_M1EZ)));
  494. C2212E_Z = C2212E_Z +
  495. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2EZ)*(N12_M1EZ)));
  496. C2212E_Z = C2212E_Z +
  497. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1EZ)*(N12_M2EZ)));
  498. C2212E_Z = C2212E_Z +
  499. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2EZ)*(N12_M2EZ)));
  500. C2212E = C2212E_X + C2212E_Y + C2212E_Z;
  501. C2212EA = EXTR
  502. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  503. 'C23E' 1 1 1;
  504. * determination du terme N11N11 couche interne
  505. C1111I_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1IX)*(N11_M1IX)));
  506. C1111I_X = C1111I_X +
  507. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2IX)*(N11_M1IX)));
  508. C1111I_X = C1111I_X +
  509. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1IX)*(N11_M2IX)));
  510. C1111I_X = C1111I_X +
  511. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2IX)*(N11_M2IX)));
  512. C1111I_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1IY)*(N11_M1IY)));
  513. C1111I_Y = C1111I_Y +
  514. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2IY)*(N11_M1IY)));
  515. C1111I_Y = C1111I_Y +
  516. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1IY)*(N11_M2IY)));
  517. C1111I_Y = C1111I_Y +
  518. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2IY)*(N11_M2IY)));
  519. C1111I_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1IZ)*(N11_M1IZ)));
  520. C1111I_Z = C1111I_Z +
  521. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2IZ)*(N11_M1IZ)));
  522. C1111I_Z = C1111I_Z +
  523. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1IZ)*(N11_M2IZ)));
  524. C1111I_Z = C1111I_Z +
  525. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2IZ)*(N11_M2IZ)));
  526. C1111I = C1111I_X + C1111I_Y + C1111I_Z;
  527. C1111IA = EXTR
  528. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  529. 'C11I' 1 1 1;
  530. * determination du terme N22N22 couche interne
  531. C2222I_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1IX)*(N22_M1IX)));
  532. C2222I_X = C2222I_X +
  533. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2IX)*(N22_M1IX)));
  534. C2222I_X = C2222I_X +
  535. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1IX)*(N22_M2IX)));
  536. C2222I_X = C2222I_X +
  537. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2IX)*(N22_M2IX)));
  538. C2222I_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1IY)*(N22_M1IY)));
  539. C2222I_Y = C2222I_Y +
  540. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2IY)*(N22_M1IY)));
  541. C2222I_Y = C2222I_Y +
  542. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1IY)*(N22_M2IY)));
  543. C2222I_Y = C2222I_Y +
  544. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2IY)*(N22_M2IY)));
  545. C2222I_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1IZ)*(N22_M1IZ)));
  546. C2222I_Z = C2222I_Z +
  547. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2IZ)*(N22_M1IZ)));
  548. C2222I_Z = C2222I_Z +
  549. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1IZ)*(N22_M2IZ)));
  550. C2222I_Z = C2222I_Z +
  551. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2IZ)*(N22_M2IZ)));
  552. C2222I = C2222I_X + C2222I_Y + C2222I_Z;
  553. C2222IA = EXTR
  554. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  555. 'C22I' 1 1 1;
  556. * determination du terme N12N12 couche interne
  557. C1212I_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N12_M1IX)*(N12_M1IX)));
  558. C1212I_X = C1212I_X +
  559. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M2IX)*(N12_M1IX)));
  560. C1212I_X = C1212I_X +
  561. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M1IX)*(N12_M2IX)));
  562. C1212I_X = C1212I_X +
  563. ((TAB2.REPRISE.PIJ. 2 . 2)*((N12_M2IX)*(N12_M2IX)));
  564. C1212I_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N12_M1IY)*(N12_M1IY)));
  565. C1212I_Y = C1212I_Y +
  566. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M2IY)*(N12_M1IY)));
  567. C1212I_Y = C1212I_Y +
  568. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M1IY)*(N12_M2IY)));
  569. C1212I_Y = C1212I_Y +
  570. ((TAB2.REPRISE.PIJ. 2 . 2)*((N12_M2IY)*(N12_M2IY)));
  571. C1212I_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N12_M1IZ)*(N12_M1IZ)));
  572. C1212I_Z = C1212I_Z +
  573. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M2IZ)*(N12_M1IZ)));
  574. C1212I_Z = C1212I_Z +
  575. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M1IZ)*(N12_M2IZ)));
  576. C1212I_Z = C1212I_Z +
  577. ((TAB2.REPRISE.PIJ. 2 . 2)*((N12_M2IZ)*(N12_M2IZ)));
  578. C1212I = C1212I_X + C1212I_Y + C1212I_Z;
  579. C1212IA = EXTR
  580. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  581. 'C33I' 1 1 1;
  582. * determination du terme N11N22 couche interne
  583. C1122I_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1IX)*(N22_M1IX)));
  584. C1122I_X = C1122I_X +
  585. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2IX)*(N22_M1IX)));
  586. C1122I_X = C1122I_X +
  587. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1IX)*(N22_M2IX)));
  588. C1122I_X = C1122I_X +
  589. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2IX)*(N22_M2IX)));
  590. C1122I_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1IY)*(N22_M1IY)));
  591. C1122I_Y = C1122I_Y +
  592. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2IY)*(N22_M1IY)));
  593. C1122I_Y = C1122I_Y +
  594. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1IY)*(N22_M2IY)));
  595. C1122I_Y = C1122I_Y +
  596. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2IY)*(N22_M2IY)));
  597. C1122I_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1IZ)*(N22_M1IZ)));
  598. C1122I_Z = C1122I_Z +
  599. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2IZ)*(N22_M1IZ)));
  600. C1122I_Z = C1122I_Z +
  601. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1IZ)*(N22_M2IZ)));
  602. C1122I_Z = C1122I_Z +
  603. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2IZ)*(N22_M2IZ)));
  604. C1122E = C1122I_X + C1122I_Y + C1122I_Z;
  605. C1122EA = EXTR
  606. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  607. 'C12I' 1 1 1;
  608. * determination du terme N11N12 couche interne
  609. C1112I_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1IX)*(N12_M1IX)));
  610. C1112I_X = C1112I_X +
  611. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2IX)*(N12_M1IX)));
  612. C1112I_X = C1112I_X +
  613. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1IX)*(N12_M2IX)));
  614. C1112I_X = C1112I_X +
  615. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2IX)*(N12_M2IX)));
  616. C1112I_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1IY)*(N12_M1IY)));
  617. C1112I_Y = C1112I_Y +
  618. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2IY)*(N12_M1IY)));
  619. C1112I_Y = C1112I_Y +
  620. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1IY)*(N12_M2IY)));
  621. C1112I_Y = C1112I_Y +
  622. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2IY)*(N12_M2IY)));
  623. C1112I_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1IZ)*(N12_M1IZ)));
  624. C1112I_Z = C1112I_Z +
  625. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2IZ)*(N12_M1IZ)));
  626. C1112I_Z = C1112I_Z +
  627. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1IZ)*(N12_M2IZ)));
  628. C1112I_Z = C1112I_Z +
  629. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2IZ)*(N12_M2IZ)));
  630. C1112I = C1112I_X + C1112I_Y + C1112I_Z;
  631. C1112IA = EXTR
  632. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  633. 'C13I' 1 1 1;
  634. * determination du terme N22N12 couche interne
  635. C2212I_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1IX)*(N12_M1IX)));
  636. C2212I_X = C2212I_X +
  637. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2IX)*(N12_M1IX)));
  638. C2212I_X = C2212I_X +
  639. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1IX)*(N12_M2IX)));
  640. C2212I_X = C2212I_X +
  641. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2IX)*(N12_M2IX)));
  642. C2212I_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1IY)*(N12_M1IY)));
  643. C2212I_Y = C2212I_Y +
  644. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2IY)*(N12_M1IY)));
  645. C2212I_Y = C2212I_Y +
  646. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1IY)*(N12_M2IY)));
  647. C2212I_Y = C2212I_Y +
  648. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2IY)*(N12_M2IY)));
  649. C2212I_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1IZ)*(N12_M1IZ)));
  650. C2212I_Z = C2212I_Z +
  651. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2IZ)*(N12_M1IZ)));
  652. C2212I_Z = C2212I_Z +
  653. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1IZ)*(N12_M2IZ)));
  654. C2212I_Z = C2212I_Z +
  655. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2IZ)*(N12_M2IZ)));
  656. C2212I = C2212I_X + C2212I_Y + C2212I_Z;
  657. C2212IA = EXTR
  658. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  659. 'C23I' 1 1 1;
  660. *
  661. * Control
  662. SI ((ABS(C1111E - C1111EA)) > 1.E-12);
  663. ERRE 5;
  664. FINSI;
  665. SI ((ABS(C2222E - C2222EA)) > 1.E-12);
  666. ERRE 5;
  667. FINSI;
  668. SI ((ABS(C1212E - C1212EA)) > 1.E-12);
  669. ERRE 5;
  670. FINSI;
  671. SI ((ABS(C1112E - C1112EA)) > 1.E-12);
  672. ERRE 5;
  673. FINSI;
  674. SI ((ABS(C2212E - C2212EA)) > 1.E-12);
  675. ERRE 5;
  676. FINSI;
  677. SI ((ABS(C1111I - C1111IA)) > 1.E-12);
  678. ERRE 5;
  679. FINSI;
  680. SI ((ABS(C2222I - C2222IA)) > 1.E-12);
  681. ERRE 5;
  682. FINSI;
  683. SI ((ABS(C1212I - C1212IA)) > 1.E-12);
  684. ERRE 5;
  685. FINSI;
  686. SI ((ABS(C1112I - C1112IA)) > 1.E-12);
  687. ERRE 5;
  688. FINSI;
  689. SI ((ABS(C2212I - C2212IA)) > 1.E-12);
  690. ERRE 5;
  691. FINSI;
  692. *
  693. FIN;
  694.  
  695.  
  696.  
  697.  
  698.  

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