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és 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 TBAS;
  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 .'TYPE' = CHAINE 'VOILE';
  183. TAB5. 1 .'MAIL' = SUR1;
  184. TAB5. 1 .'MODE' = MOD1;
  185. TAB5. 1 .'MATE' = MAT1;
  186. TAB5. 1 .'EPAISSEUR' = H1;
  187. TAB5. 1 .'ENROBAGE_E' = ER1;
  188. TAB5. 1 .'ENROBAGE_I' = ER2;
  189. TAB5. 1 .'COT_THETA' = COT1;
  190. TAB5. 1 .'DIR1' = (1.0 0. 0.);
  191. TAB5. 1 .'DIR2' = (0.0 1. 0.);
  192. *
  193. * Calcul SISSIB
  194. TAB2 = SISSIB TAB1 TAB5;
  195. *
  196. * Determination analique des efforts selon
  197. * le modele de Marti pour l'element 1 pour
  198. * chaque champ des efforts modaux
  199. N11_M1E = N11_M1/2;
  200. N22_M1E = N22_M1/2;
  201. N12_M1E = N12_M1/2;
  202. N11_M1I = N11_M1/2;
  203. N22_M1I = N22_M1/2;
  204. N12_M1I = N12_M1/2;
  205. N11_M2E = N11_M2/2;
  206. N22_M2E = N22_M2/2;
  207. N12_M2E = N12_M2/2;
  208. N11_M2I = N11_M2/2;
  209. N22_M2I = N22_M2/2;
  210. N12_M2I = N12_M2/2;
  211. *
  212. * On extrait les efforts selon MARTI deduits par
  213. * la procedure pour chaque champ des efforts modaux - Mode 1
  214. N11_M1ET = EXTR
  215. (REDU (TB1.MODES. 1 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  216. 'N11E' 1 1 1;
  217. N22_M1ET = EXTR
  218. (REDU (TB1.MODES. 1 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  219. 'N22E' 1 1 1;
  220. N12_M1ET = EXTR
  221. (REDU (TB1.MODES. 1 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  222. 'N12E' 1 1 1;
  223. N11_M1IT = EXTR
  224. (REDU (TB1.MODES. 1 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  225. 'N11I' 1 1 1;
  226. N22_M1IT = EXTR
  227. (REDU (TB1.MODES. 1 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  228. 'N22I' 1 1 1;
  229. N12_M1IT = EXTR
  230. (REDU (TB1.MODES. 1 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  231. 'N12I' 1 1 1;
  232. * On extrait les efforts selon MARTI deduits par
  233. * la procedure pour chaque champ des efforts modaux - Mode 2
  234. N11_M2ET = EXTR
  235. (REDU (TB1.MODES. 2 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  236. 'N11E' 1 1 1;
  237. N22_M2ET = EXTR
  238. (REDU (TB1.MODES. 2 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  239. 'N22E' 1 1 1;
  240. N12_M2ET = EXTR
  241. (REDU (TB1.MODES. 2 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  242. 'N12E' 1 1 1;
  243. N11_M2IT = EXTR
  244. (REDU (TB1.MODES. 2 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  245. 'N11I' 1 1 1;
  246. N22_M2IT = EXTR
  247. (REDU (TB1.MODES. 2 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  248. 'N22I' 1 1 1;
  249. N12_M2IT = EXTR
  250. (REDU (TB1.MODES. 2 .CONTRAINTE_MEMBRANE.VOI1) (SUR1 ELEM 1))
  251. 'N12I' 1 1 1;
  252. *
  253. * Control
  254. SI ((ABS(N11_M1ET - N11_M1E)) > 1.E-12);
  255. ERRE 5;
  256. FINSI;
  257. SI ((ABS(N22_M1ET - N22_M1E)) > 1.E-12);
  258. ERRE 5;
  259. FINSI;
  260. SI ((ABS(N12_M1ET - N12_M1E)) > 1.E-12);
  261. ERRE 5;
  262. FINSI;
  263. SI ((ABS(N11_M1IT - N11_M1I)) > 1.E-12);
  264. ERRE 5;
  265. FINSI;
  266. SI ((ABS(N22_M1IT - N22_M1I)) > 1.E-12);
  267. ERRE 5;
  268. FINSI;
  269. SI ((ABS(N12_M1IT - N12_M1I)) > 1.E-12);
  270. ERRE 5;
  271. FINSI;
  272. *
  273. SI ((ABS(N11_M2ET - N11_M2E)) > 1.E-12);
  274. ERRE 5;
  275. FINSI;
  276. SI ((ABS(N22_M2ET - N22_M2E)) > 1.E-12);
  277. ERRE 5;
  278. FINSI;
  279. SI ((ABS(N12_M2ET - N12_M2E)) > 1.E-12);
  280. ERRE 5;
  281. FINSI;
  282. SI ((ABS(N11_M2IT - N11_M2I)) > 1.E-12);
  283. ERRE 5;
  284. FINSI;
  285. SI ((ABS(N22_M2IT - N22_M2I)) > 1.E-12);
  286. ERRE 5;
  287. FINSI;
  288. SI ((ABS(N12_M2IT - N12_M2I)) > 1.E-12);
  289. ERRE 5;
  290. FINSI;
  291. *
  292. * Determination de la Pseudo acceleration pour les Modes 1 et 2
  293. * pour les direction x, y et z
  294. PSE_A1X = TAB2.'REPRISE'.'SPECTRE'. 1 . 1;
  295. PSE_A1Y = TAB2.'REPRISE'.'SPECTRE'. 1 . 2;
  296. PSE_A1Z = TAB2.'REPRISE'.'SPECTRE'. 1 . 3;
  297. PSE_A2X = TAB2.'REPRISE'.'SPECTRE'. 2 . 1;
  298. PSE_A2Y = TAB2.'REPRISE'.'SPECTRE'. 2 . 2;
  299. PSE_A2Z = TAB2.'REPRISE'.'SPECTRE'. 2 . 3;
  300. *
  301. * Determination de l'amplification pour les Modes 1 et 2
  302. * pour les direction x, y et z
  303. AMP_1X = (TB1.MODES. 1 .FACT_PART_NORM. 1)*(PSE_A1X);
  304. AMP_1Y = (TB1.MODES. 1 .FACT_PART_NORM. 2)*(PSE_A1Y);
  305. AMP_1Z = (TB1.MODES. 1 .FACT_PART_NORM. 3)*(PSE_A1Z);
  306. AMP_2X = (TB1.MODES. 2 .FACT_PART_NORM. 1)*(PSE_A2X);
  307. AMP_2Y = (TB1.MODES. 2 .FACT_PART_NORM. 2)*(PSE_A2Y);
  308. AMP_2Z = (TB1.MODES. 2 .FACT_PART_NORM. 3)*(PSE_A2Z);
  309. *
  310. * Efforts projetés amplifiés (mode 1 ; directions x, y, z)
  311. N11_M1EX = N11_M1E*AMP_1X;
  312. N22_M1EX = N22_M1E*AMP_1X;
  313. N12_M1EX = N12_M1E*AMP_1X;
  314. N11_M1IX = N11_M1I*AMP_1X;
  315. N22_M1IX = N22_M1I*AMP_1X;
  316. N12_M1IX = N12_M1I*AMP_1X;
  317. N11_M1EY = N11_M1E*AMP_1Y;
  318. N22_M1EY = N22_M1E*AMP_1Y;
  319. N12_M1EY = N12_M1E*AMP_1Y;
  320. N11_M1IY = N11_M1I*AMP_1Y;
  321. N22_M1IY = N22_M1I*AMP_1Y;
  322. N12_M1IY = N12_M1I*AMP_1Y;
  323. N11_M1EZ = N11_M1E*AMP_1Z;
  324. N22_M1EZ = N22_M1E*AMP_1Z;
  325. N12_M1EZ = N12_M1E*AMP_1Z;
  326. N11_M1IZ = N11_M1I*AMP_1Z;
  327. N22_M1IZ = N22_M1I*AMP_1Z;
  328. N12_M1IZ = N12_M1I*AMP_1Z;
  329. * Efforts projetés amplifiés (mode 2 ; directions x, y, z)
  330. N11_M2EX = N11_M2E*AMP_2X;
  331. N22_M2EX = N22_M2E*AMP_2X;
  332. N12_M2EX = N12_M2E*AMP_2X;
  333. N11_M2IX = N11_M2I*AMP_2X;
  334. N22_M2IX = N22_M2I*AMP_2X;
  335. N12_M2IX = N12_M2I*AMP_2X;
  336. N11_M2EY = N11_M2E*AMP_2Y;
  337. N22_M2EY = N22_M2E*AMP_2Y;
  338. N12_M2EY = N12_M2E*AMP_2Y;
  339. N11_M2IY = N11_M2I*AMP_2Y;
  340. N22_M2IY = N22_M2I*AMP_2Y;
  341. N12_M2IY = N12_M2I*AMP_2Y;
  342. N11_M2EZ = N11_M2E*AMP_2Z;
  343. N22_M2EZ = N22_M2E*AMP_2Z;
  344. N12_M2EZ = N12_M2E*AMP_2Z;
  345. N11_M2IZ = N11_M2I*AMP_2Z;
  346. N22_M2IZ = N22_M2I*AMP_2Z;
  347. N12_M2IZ = N12_M2I*AMP_2Z;
  348. *
  349. * determination du terme N11N11 couche externe
  350. C1111E_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1EX)*(N11_M1EX)));
  351. C1111E_X = C1111E_X +
  352. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2EX)*(N11_M1EX)));
  353. C1111E_X = C1111E_X +
  354. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1EX)*(N11_M2EX)));
  355. C1111E_X = C1111E_X +
  356. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2EX)*(N11_M2EX)));
  357. C1111E_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1EY)*(N11_M1EY)));
  358. C1111E_Y = C1111E_Y +
  359. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2EY)*(N11_M1EY)));
  360. C1111E_Y = C1111E_Y +
  361. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1EY)*(N11_M2EY)));
  362. C1111E_Y = C1111E_Y +
  363. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2EY)*(N11_M2EY)));
  364. C1111E_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1EZ)*(N11_M1EZ)));
  365. C1111E_Z = C1111E_Z +
  366. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2EZ)*(N11_M1EZ)));
  367. C1111E_Z = C1111E_Z +
  368. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1EZ)*(N11_M2EZ)));
  369. C1111E_Z = C1111E_Z +
  370. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2EZ)*(N11_M2EZ)));
  371. C1111E = C1111E_X + C1111E_Y + C1111E_Z;
  372. C1111EA = EXTR
  373. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  374. 'C11E' 1 1 1;
  375. * determination du terme N22N22 couche externe
  376. C2222E_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1EX)*(N22_M1EX)));
  377. C2222E_X = C2222E_X +
  378. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2EX)*(N22_M1EX)));
  379. C2222E_X = C2222E_X +
  380. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1EX)*(N22_M2EX)));
  381. C2222E_X = C2222E_X +
  382. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2EX)*(N22_M2EX)));
  383. C2222E_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1EY)*(N22_M1EY)));
  384. C2222E_Y = C2222E_Y +
  385. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2EY)*(N22_M1EY)));
  386. C2222E_Y = C2222E_Y +
  387. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1EY)*(N22_M2EY)));
  388. C2222E_Y = C2222E_Y +
  389. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2EY)*(N22_M2EY)));
  390. C2222E_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1EZ)*(N22_M1EZ)));
  391. C2222E_Z = C2222E_Z +
  392. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2EZ)*(N22_M1EZ)));
  393. C2222E_Z = C2222E_Z +
  394. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1EZ)*(N22_M2EZ)));
  395. C2222E_Z = C2222E_Z +
  396. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2EZ)*(N22_M2EZ)));
  397. C2222E = C2222E_X + C2222E_Y + C2222E_Z;
  398. C2222EA = EXTR
  399. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  400. 'C22E' 1 1 1;
  401. * determination du terme N12N12 couche externe
  402. C1212E_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N12_M1EX)*(N12_M1EX)));
  403. C1212E_X = C1212E_X +
  404. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M2EX)*(N12_M1EX)));
  405. C1212E_X = C1212E_X +
  406. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M1EX)*(N12_M2EX)));
  407. C1212E_X = C1212E_X +
  408. ((TAB2.REPRISE.PIJ. 2 . 2)*((N12_M2EX)*(N12_M2EX)));
  409. C1212E_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N12_M1EY)*(N12_M1EY)));
  410. C1212E_Y = C1212E_Y +
  411. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M2EY)*(N12_M1EY)));
  412. C1212E_Y = C1212E_Y +
  413. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M1EY)*(N12_M2EY)));
  414. C1212E_Y = C1212E_Y +
  415. ((TAB2.REPRISE.PIJ. 2 . 2)*((N12_M2EY)*(N12_M2EY)));
  416. C1212E_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N12_M1EZ)*(N12_M1EZ)));
  417. C1212E_Z = C1212E_Z +
  418. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M2EZ)*(N12_M1EZ)));
  419. C1212E_Z = C1212E_Z +
  420. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M1EZ)*(N12_M2EZ)));
  421. C1212E_Z = C1212E_Z +
  422. ((TAB2.REPRISE.PIJ. 2 . 2)*((N12_M2EZ)*(N12_M2EZ)));
  423. C1212E = C1212E_X + C1212E_Y + C1212E_Z;
  424. C1212EA = EXTR
  425. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  426. 'C33E' 1 1 1;
  427. * determination du terme N11N22 couche externe
  428. C1122E_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1EX)*(N22_M1EX)));
  429. C1122E_X = C1122E_X +
  430. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2EX)*(N22_M1EX)));
  431. C1122E_X = C1122E_X +
  432. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1EX)*(N22_M2EX)));
  433. C1122E_X = C1122E_X +
  434. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2EX)*(N22_M2EX)));
  435. C1122E_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1EY)*(N22_M1EY)));
  436. C1122E_Y = C1122E_Y +
  437. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2EY)*(N22_M1EY)));
  438. C1122E_Y = C1122E_Y +
  439. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1EY)*(N22_M2EY)));
  440. C1122E_Y = C1122E_Y +
  441. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2EY)*(N22_M2EY)));
  442. C1122E_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1EZ)*(N22_M1EZ)));
  443. C1122E_Z = C1122E_Z +
  444. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2EZ)*(N22_M1EZ)));
  445. C1122E_Z = C1122E_Z +
  446. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1EZ)*(N22_M2EZ)));
  447. C1122E_Z = C1122E_Z +
  448. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2EZ)*(N22_M2EZ)));
  449. C1122E = C1122E_X + C1122E_Y + C1122E_Z;
  450. C1122EA = EXTR
  451. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  452. 'C12E' 1 1 1;
  453. * determination du terme N11N12 couche externe
  454. C1112E_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1EX)*(N12_M1EX)));
  455. C1112E_X = C1112E_X +
  456. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2EX)*(N12_M1EX)));
  457. C1112E_X = C1112E_X +
  458. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1EX)*(N12_M2EX)));
  459. C1112E_X = C1112E_X +
  460. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2EX)*(N12_M2EX)));
  461. C1112E_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1EY)*(N12_M1EY)));
  462. C1112E_Y = C1112E_Y +
  463. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2EY)*(N12_M1EY)));
  464. C1112E_Y = C1112E_Y +
  465. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1EY)*(N12_M2EY)));
  466. C1112E_Y = C1112E_Y +
  467. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2EY)*(N12_M2EY)));
  468. C1112E_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1EZ)*(N12_M1EZ)));
  469. C1112E_Z = C1112E_Z +
  470. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2EZ)*(N12_M1EZ)));
  471. C1112E_Z = C1112E_Z +
  472. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1EZ)*(N12_M2EZ)));
  473. C1112E_Z = C1112E_Z +
  474. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2EZ)*(N12_M2EZ)));
  475. C1112E = C1112E_X + C1112E_Y + C1112E_Z;
  476. C1112EA = EXTR
  477. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  478. 'C13E' 1 1 1;
  479. * determination du terme N22N12 couche externe
  480. C2212E_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1EX)*(N12_M1EX)));
  481. C2212E_X = C2212E_X +
  482. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2EX)*(N12_M1EX)));
  483. C2212E_X = C2212E_X +
  484. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1EX)*(N12_M2EX)));
  485. C2212E_X = C2212E_X +
  486. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2EX)*(N12_M2EX)));
  487. C2212E_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1EY)*(N12_M1EY)));
  488. C2212E_Y = C2212E_Y +
  489. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2EY)*(N12_M1EY)));
  490. C2212E_Y = C2212E_Y +
  491. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1EY)*(N12_M2EY)));
  492. C2212E_Y = C2212E_Y +
  493. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2EY)*(N12_M2EY)));
  494. C2212E_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1EZ)*(N12_M1EZ)));
  495. C2212E_Z = C2212E_Z +
  496. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2EZ)*(N12_M1EZ)));
  497. C2212E_Z = C2212E_Z +
  498. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1EZ)*(N12_M2EZ)));
  499. C2212E_Z = C2212E_Z +
  500. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2EZ)*(N12_M2EZ)));
  501. C2212E = C2212E_X + C2212E_Y + C2212E_Z;
  502. C2212EA = EXTR
  503. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  504. 'C23E' 1 1 1;
  505. * determination du terme N11N11 couche interne
  506. C1111I_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1IX)*(N11_M1IX)));
  507. C1111I_X = C1111I_X +
  508. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2IX)*(N11_M1IX)));
  509. C1111I_X = C1111I_X +
  510. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1IX)*(N11_M2IX)));
  511. C1111I_X = C1111I_X +
  512. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2IX)*(N11_M2IX)));
  513. C1111I_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1IY)*(N11_M1IY)));
  514. C1111I_Y = C1111I_Y +
  515. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2IY)*(N11_M1IY)));
  516. C1111I_Y = C1111I_Y +
  517. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1IY)*(N11_M2IY)));
  518. C1111I_Y = C1111I_Y +
  519. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2IY)*(N11_M2IY)));
  520. C1111I_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1IZ)*(N11_M1IZ)));
  521. C1111I_Z = C1111I_Z +
  522. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2IZ)*(N11_M1IZ)));
  523. C1111I_Z = C1111I_Z +
  524. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1IZ)*(N11_M2IZ)));
  525. C1111I_Z = C1111I_Z +
  526. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2IZ)*(N11_M2IZ)));
  527. C1111I = C1111I_X + C1111I_Y + C1111I_Z;
  528. C1111IA = EXTR
  529. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  530. 'C11I' 1 1 1;
  531. * determination du terme N22N22 couche interne
  532. C2222I_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1IX)*(N22_M1IX)));
  533. C2222I_X = C2222I_X +
  534. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2IX)*(N22_M1IX)));
  535. C2222I_X = C2222I_X +
  536. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1IX)*(N22_M2IX)));
  537. C2222I_X = C2222I_X +
  538. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2IX)*(N22_M2IX)));
  539. C2222I_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1IY)*(N22_M1IY)));
  540. C2222I_Y = C2222I_Y +
  541. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2IY)*(N22_M1IY)));
  542. C2222I_Y = C2222I_Y +
  543. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1IY)*(N22_M2IY)));
  544. C2222I_Y = C2222I_Y +
  545. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2IY)*(N22_M2IY)));
  546. C2222I_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1IZ)*(N22_M1IZ)));
  547. C2222I_Z = C2222I_Z +
  548. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2IZ)*(N22_M1IZ)));
  549. C2222I_Z = C2222I_Z +
  550. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1IZ)*(N22_M2IZ)));
  551. C2222I_Z = C2222I_Z +
  552. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2IZ)*(N22_M2IZ)));
  553. C2222I = C2222I_X + C2222I_Y + C2222I_Z;
  554. C2222IA = EXTR
  555. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  556. 'C22I' 1 1 1;
  557. * determination du terme N12N12 couche interne
  558. C1212I_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N12_M1IX)*(N12_M1IX)));
  559. C1212I_X = C1212I_X +
  560. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M2IX)*(N12_M1IX)));
  561. C1212I_X = C1212I_X +
  562. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M1IX)*(N12_M2IX)));
  563. C1212I_X = C1212I_X +
  564. ((TAB2.REPRISE.PIJ. 2 . 2)*((N12_M2IX)*(N12_M2IX)));
  565. C1212I_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N12_M1IY)*(N12_M1IY)));
  566. C1212I_Y = C1212I_Y +
  567. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M2IY)*(N12_M1IY)));
  568. C1212I_Y = C1212I_Y +
  569. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M1IY)*(N12_M2IY)));
  570. C1212I_Y = C1212I_Y +
  571. ((TAB2.REPRISE.PIJ. 2 . 2)*((N12_M2IY)*(N12_M2IY)));
  572. C1212I_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N12_M1IZ)*(N12_M1IZ)));
  573. C1212I_Z = C1212I_Z +
  574. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M2IZ)*(N12_M1IZ)));
  575. C1212I_Z = C1212I_Z +
  576. ((TAB2.REPRISE.PIJ. 2 . 1)*((N12_M1IZ)*(N12_M2IZ)));
  577. C1212I_Z = C1212I_Z +
  578. ((TAB2.REPRISE.PIJ. 2 . 2)*((N12_M2IZ)*(N12_M2IZ)));
  579. C1212I = C1212I_X + C1212I_Y + C1212I_Z;
  580. C1212IA = EXTR
  581. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  582. 'C33I' 1 1 1;
  583. * determination du terme N11N22 couche interne
  584. C1122I_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1IX)*(N22_M1IX)));
  585. C1122I_X = C1122I_X +
  586. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2IX)*(N22_M1IX)));
  587. C1122I_X = C1122I_X +
  588. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1IX)*(N22_M2IX)));
  589. C1122I_X = C1122I_X +
  590. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2IX)*(N22_M2IX)));
  591. C1122I_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1IY)*(N22_M1IY)));
  592. C1122I_Y = C1122I_Y +
  593. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2IY)*(N22_M1IY)));
  594. C1122I_Y = C1122I_Y +
  595. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1IY)*(N22_M2IY)));
  596. C1122I_Y = C1122I_Y +
  597. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2IY)*(N22_M2IY)));
  598. C1122I_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1IZ)*(N22_M1IZ)));
  599. C1122I_Z = C1122I_Z +
  600. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2IZ)*(N22_M1IZ)));
  601. C1122I_Z = C1122I_Z +
  602. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1IZ)*(N22_M2IZ)));
  603. C1122I_Z = C1122I_Z +
  604. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2IZ)*(N22_M2IZ)));
  605. C1122E = C1122I_X + C1122I_Y + C1122I_Z;
  606. C1122EA = EXTR
  607. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  608. 'C12I' 1 1 1;
  609. * determination du terme N11N12 couche interne
  610. C1112I_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1IX)*(N12_M1IX)));
  611. C1112I_X = C1112I_X +
  612. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2IX)*(N12_M1IX)));
  613. C1112I_X = C1112I_X +
  614. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1IX)*(N12_M2IX)));
  615. C1112I_X = C1112I_X +
  616. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2IX)*(N12_M2IX)));
  617. C1112I_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1IY)*(N12_M1IY)));
  618. C1112I_Y = C1112I_Y +
  619. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2IY)*(N12_M1IY)));
  620. C1112I_Y = C1112I_Y +
  621. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1IY)*(N12_M2IY)));
  622. C1112I_Y = C1112I_Y +
  623. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2IY)*(N12_M2IY)));
  624. C1112I_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N11_M1IZ)*(N12_M1IZ)));
  625. C1112I_Z = C1112I_Z +
  626. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M2IZ)*(N12_M1IZ)));
  627. C1112I_Z = C1112I_Z +
  628. ((TAB2.REPRISE.PIJ. 2 . 1)*((N11_M1IZ)*(N12_M2IZ)));
  629. C1112I_Z = C1112I_Z +
  630. ((TAB2.REPRISE.PIJ. 2 . 2)*((N11_M2IZ)*(N12_M2IZ)));
  631. C1112I = C1112I_X + C1112I_Y + C1112I_Z;
  632. C1112IA = EXTR
  633. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  634. 'C13I' 1 1 1;
  635. * determination du terme N22N12 couche interne
  636. C2212I_X = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1IX)*(N12_M1IX)));
  637. C2212I_X = C2212I_X +
  638. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2IX)*(N12_M1IX)));
  639. C2212I_X = C2212I_X +
  640. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1IX)*(N12_M2IX)));
  641. C2212I_X = C2212I_X +
  642. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2IX)*(N12_M2IX)));
  643. C2212I_Y = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1IY)*(N12_M1IY)));
  644. C2212I_Y = C2212I_Y +
  645. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2IY)*(N12_M1IY)));
  646. C2212I_Y = C2212I_Y +
  647. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1IY)*(N12_M2IY)));
  648. C2212I_Y = C2212I_Y +
  649. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2IY)*(N12_M2IY)));
  650. C2212I_Z = ((TAB2.REPRISE.PIJ. 1 . 1)*((N22_M1IZ)*(N12_M1IZ)));
  651. C2212I_Z = C2212I_Z +
  652. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M2IZ)*(N12_M1IZ)));
  653. C2212I_Z = C2212I_Z +
  654. ((TAB2.REPRISE.PIJ. 2 . 1)*((N22_M1IZ)*(N12_M2IZ)));
  655. C2212I_Z = C2212I_Z +
  656. ((TAB2.REPRISE.PIJ. 2 . 2)*((N22_M2IZ)*(N12_M2IZ)));
  657. C2212I = C2212I_X + C2212I_Y + C2212I_Z;
  658. C2212IA = EXTR
  659. (REDU (TAB2.REPONSE_TOTALE.COVARIANCE.VOI1) (SUR1 ELEM 1))
  660. 'C23I' 1 1 1;
  661. *
  662. * Control
  663. SI ((ABS(C1111E - C1111EA)) > 1.E-12);
  664. ERRE 5;
  665. FINSI;
  666. SI ((ABS(C2222E - C2222EA)) > 1.E-12);
  667. ERRE 5;
  668. FINSI;
  669. SI ((ABS(C1212E - C1212EA)) > 1.E-12);
  670. ERRE 5;
  671. FINSI;
  672. SI ((ABS(C1112E - C1112EA)) > 1.E-12);
  673. ERRE 5;
  674. FINSI;
  675. SI ((ABS(C2212E - C2212EA)) > 1.E-12);
  676. ERRE 5;
  677. FINSI;
  678. SI ((ABS(C1111I - C1111IA)) > 1.E-12);
  679. ERRE 5;
  680. FINSI;
  681. SI ((ABS(C2222I - C2222IA)) > 1.E-12);
  682. ERRE 5;
  683. FINSI;
  684. SI ((ABS(C1212I - C1212IA)) > 1.E-12);
  685. ERRE 5;
  686. FINSI;
  687. SI ((ABS(C1112I - C1112IA)) > 1.E-12);
  688. ERRE 5;
  689. FINSI;
  690. SI ((ABS(C2212I - C2212IA)) > 1.E-12);
  691. ERRE 5;
  692. FINSI;
  693. *
  694. FIN;
  695.  
  696.  
  697.  
  698.  
  699.  

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