Télécharger la plaquette Cast3M


fr en

Télécharger sissib_cov.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : sissib_cov.dgibi
  2. OPTI DIME 3 MODE TRID ELEM QUA4;
  3. *
  4. **********************************************************
  5. *  Test sur la procedure SISSIB, suite à la modification  
  6. *  pour le calcul des efforts maximaux et des covariances sur 
  7. *  base modale des efforts generalisés pour les coques du modele 
  8. *  a trois couches de MARTI
  9. *
  10. *  On considere une plaque de 15m de larguer et 1m de profondeur
  11. *  soumise à trois mouvements en direction x, y et z (reppresentés
  12. *    par les spectres de reponse)
  13. *  La plaque est encastrée sur un coté et libre sur les autres. 
  14. *  Les termes de masse sont presents seulement sur un coté.
  15. *  Le comportment dynamique est reppresenté par deux modes de 
  16. *  vibration.
  17. *
  18. *  Develloppé par Alberto FRAU /DEN/DANS/DM2S/SEMT/EMSI
  19. *              et Nicolas ILE  /DEN/DANS/DM2S/SEMT/EMSI
  20. *
  21. **********************************************************
  22. *
  23. * option graphiques
  24. GRAP1 = 0;
  25. *
  26. * Proprietés geometriques
  27. H1 = 0.1;
  28. ER1 = 0.025;
  29. ER2 = 0.025;
  30. COT1 = 1.43;
  31. *
  32. * Definition du chargement sismique x, y et z  
  33. LIS_FREQ = PROG 0.1 0.25 2.5 8.0 30.0 40.0 50.0 400.0 ;
  34. * liste des acc?l?rations pour seisme suivant x et y [m.s-2]
  35. LIS_ACCX = PROG 0.0175 0.0175 0.184 0.208 0.1 0.1 0.1 0.1 ;
  36. * liste des acc?l?rations pour s?isme suivant z [m.s-2]
  37. LIS_ACCZ = (LIS_ACCX)*0.6666667 ;
  38. *
  39. LIS_ACCX = (LIS_ACCX * 9.810)*3.0 ;
  40. LIS_ACCZ = (LIS_ACCZ * 9.810)*3.0 ;
  41. *
  42. SEISME_X = EVOL MANU '[Hz]' LIS_FREQ '[m.s-2]' LIS_ACCX ;
  43. SEISME_Y = EVOL MANU '[Hz]' LIS_FREQ '[m.s-2]' LIS_ACCX ;
  44. SEISME_Z = EVOL ROUG MANU '[Hz]' LIS_FREQ '[m.s-2]' LIS_ACCZ ;
  45. *
  46. * -------------------------------------------------------------
  47. P1 = 0. 0. 0.;
  48. P2 = 15. 0. 0.;
  49. P3 = 15. 1. 0.;
  50. P4 = 0. 1. 0.;
  51. *
  52. L1 = D 15 P1 P2;
  53. L2 = D 3 P2 P3;
  54. L3 = D 15 P3 P4;
  55. L3 = D 15 P3 P4;
  56. L4 = D 3 P4 P1;
  57. *
  58. PT1 = 3. 0. 0.;
  59. PT2 = 3. 1. 0.;
  60. LT1 = D 1 PT1 PT2;
  61. *
  62. * definition maillage
  63. SUR1 = DALL L1 L2 L3 L4;
  64. ELIM (SUR1 ET LT1) 0.0001;
  65. *
  66. * modele et materiau
  67. MOD1 = MODE SUR1 MECANIQUE ELASTIQUE ISOTROPE COQ4;
  68. MAT1 = MATE MOD1 YOUNG 30000.E6 NU 0.2 RHO 0. EPAI H1;
  69. *
  70. * matrices de masse et rigidite
  71. RIG1 = RIGI MOD1 MAT1;
  72. MAS1 = (MASS 'UX' L2 80000.);
  73. *
  74. * blocages
  75. BL1 = (BLOQUER 'DEPLA' 'ROTA' L4);
  76. *      (BLOQUER 'UY' L3) ET (BLOQUER 'UY' L1);
  77. *
  78. * definition de la base modale
  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.  
  700.  
  701.  
  702.  
  703.  

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