Télécharger formation_debutant_3_mecanique.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : formation_debutant_3_mecanique.dgibi
  2. *
  3. ************************************************************************
  4. * FORMATION DEBUTER AVEC CAST3M - CALCULS THERMO-MECANIQUES *
  5. * *
  6. * Modelisation du comportement thermo-mecanique d'une plaque perforee *
  7. * *
  8. * Ce fichier est la partie 3 sur 3 de la formation debutant a Cast3M *
  9. * *
  10. * Partie 3 : calculs mecaniques *
  11. * restitution des objets *
  12. * elasticite lineaire, plasticite parfaite *
  13. * chargement thermique *
  14. * caracteristiques variables *
  15. * creation de procedure Gibiane *
  16. * *
  17. ************************************************************************
  18.  
  19.  
  20.  
  21.  
  22. * RESTITUTION DES OBJETS CREES
  23. * ----------------------------
  24. OPTI 'REST' 'formation_debutant_2_thermique.sauv' ;
  25.  
  26.  
  27. * OBJET LOGIQUE POUR INDIQUER SI L'ON FAIT LES TRACES OU NON
  28. * ----------------------------------------------------------
  29. ITRAC = FAUX ;
  30.  
  31.  
  32.  
  33.  
  34.  
  35.  
  36. ************************************************************************
  37. * CHAPITRE 5 : MECANIQUE ELASTIQUE LINEAIRE *
  38. ************************************************************************
  39.  
  40.  
  41. * HYPOTHESE SUR LES CONTRAINTES PLANES
  42. * ------------------------------------
  43. OPTI 'MODE' 'PLAN' 'CONT' ;
  44.  
  45.  
  46. * MODELE MECANIQUE A MATERIAU CONSTANT
  47. * ------------------------------------
  48. MOM1 = MODE SU 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' ;
  49. MAM1 = MATE MOM1 'YOUN' YOUNGMAT 'NU' NUMAT 'ALPH' ALPHAMAT ;
  50.  
  51.  
  52. * CALCUL DE LA MATRICE DE RAIDEUR
  53. * -------------------------------
  54. RI = RIGI MOM1 MAM1 ;
  55.  
  56.  
  57. * CONDITIONS AUX LIMITES: DEPLACEMENT IMPOSE
  58. * ------------------------------------------
  59. BLMX = BLOQ PA 'UX' ;
  60. BLMY = BLOQ LBAS 'UY' ;
  61. RITOT = RI ET BLMX ET BLMY ;
  62.  
  63.  
  64. * CONDITIONS AUX LIMITES: TRACTION IMPOSEE
  65. * ----------------------------------------
  66. TR = PRES 'MASS' MOM1 LHAUT (-100.E6) ;
  67.  
  68.  
  69. * CALCUL DU CHAMP DE DEPLACEMENT PAR APPEL AU SOLVEUR
  70. * ---------------------------------------------------
  71. U5 = RESO RITOT TR ;
  72.  
  73.  
  74. * POST TRAITEMENT
  75. * ---------------
  76. * DEFORMATIONS
  77. DEF = EPSI U5 MOM1 ;
  78. * CONTRAINTES (PAR DEUX OPERATEURS => RESULTATS EQUIVALENTS)
  79. SIG = ELAS DEF MOM1 MAM1 ;
  80. SIG = SIGM U5 MOM1 MAM1 ;
  81. * DEFORMEES
  82. DEF_5 = DEFO SU U5 150. 'ROUG' ;
  83. DEF_INI = DEFO SU U5 0. ;
  84. SI ITRAC ;
  85. TRAC (DEF_INI ET DEF_5) 'TITR'
  86. '[5] Deformee elastique sous traction uniforme' ;
  87. FINSI ;
  88. DEF_5C = DEFO CSU U5 150. 'ROUG' ;
  89. DEF_INIC = DEFO CSU U5 0. ;
  90. SI ITRAC ;
  91. TRAC (DEF_INIC ET DEF_5C) 'TITR'
  92. '[5] Deformee du contour sous traction uniforme' ;
  93. FINSI ;
  94.  
  95. SI ITRAC ;
  96. TRAC SIG MOM1 CSU 15 'TITR' '[5] Contraintes, traction uniforme' ;
  97. FINSI ;
  98. DEF_5B = DEFO SU U5 150. ;
  99. SI ITRAC ;
  100. TRAC SIG MOM1 DEF_5B CSU 15 'TITR'
  101. '[5] Contraintes, traction uniforme, maillage deforme' ;
  102. FINSI ;
  103. * EVOLUTION DE LA CONCENTRATION DE CONTRAINTE DU TROU VERS LE BORD
  104. SIGB = CHAN 'CHPO' SIG MOM1 ;
  105. LBASG = LBAS ELEM 'COMP' PB PA ;
  106. EVSIG = EVOL 'CHPO' SIGB 'SMYY' LBASG ;
  107. EVK = EVSIG / 100.E6 ;
  108. SI ITRAC ;
  109. DESS EVK 'TITR' '[5] Concentration de contrainte le long de LBAS' ;
  110. FINSI ;
  111.  
  112. * EFFORTS DE REACTION AUX APPUIS
  113. REAC1 = REAC U5 (BLMX ET BLMY) ;
  114. VREAC = VECT REAC1 'FORC' 'ROUG' ;
  115. VFIMP = VECT TR 'FORC' 'VERT' ;
  116. SI ITRAC ;
  117. TRAC (VFIMP ET VREAC) CSU 'TITR'
  118. '[5] Efforts imposes(V) et reactions aux appuis(R)' ;
  119. FINSI ;
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128. ************************************************************************
  129. * CHAPITRE 6 : MECANIQUE ELASTIQUE LINEAIRE *
  130. * CHARGEMENT THERMIQUE *
  131. ************************************************************************
  132.  
  133.  
  134. * DEFORMATIONS THERMIQUES PURES AVEC LE DERNIER CHAMP DE TEMPERATURE
  135. * DU CALCUL THERMIQUE AVEC CONVECTION + RAYONNEMENT
  136. DELTA_TE = (TAB1 . 'TEMPERATURES' . (N1 - 1)) - T_INI ;
  137. EPT = EPTH DELTA_TE MOM1 MAM1 ;
  138.  
  139.  
  140. * PSEUDO CONTRAINTES CORRESPONDANT A CES DEFORMATIONS THERMIQUES
  141. * --------------------------------------------------------------
  142. SIT = ELAS EPT MOM1 MAM1 ;
  143.  
  144.  
  145. * FORCES NODALES POUR CETTE DEFORMATION THERMIQUE
  146. * -----------------------------------------------
  147. FFT = BSIG SIT MOM1 MAM1 ;
  148.  
  149.  
  150. * DEPLACEMENTS, PAR APPEL AU SOLVEUR EN SUPERPOSANT LES FORCES
  151. * PUREMENT MECANIQUES ET LES PSEUDO FORCES THERMIQUES
  152. * ------------------------------------------------------------
  153. U6 = RESO RITOT (TR ET FFT) ;
  154.  
  155.  
  156. * DEFORMEE THERMOMECANIQUE
  157. * ------------------------
  158. DEF_6 = DEFO SU U6 150. ;
  159. DEF_6C = DEFO CSU U6 150. 'ORAN' ;
  160. SI ITRAC ;
  161. TRAC (DEF_INIC ET DEF_5C ET DEF_6C) 'TITR'
  162. '[6] Deformees, traction seule(R), traction + temperature(O)' ;
  163. FINSI ;
  164.  
  165.  
  166. * DEFORMATIONS TOTALES
  167. * --------------------
  168. EP = EPSI U6 MOM1 ;
  169.  
  170.  
  171. * DEFORMATIONS ELASTIQUES
  172. * -----------------------
  173. EPE = EP - EPT ;
  174.  
  175.  
  176. * LES CONTRAINTES SONT CALCULEES A PARTIR DES DEFORMATIONS ELASTIQUES
  177. * -------------------------------------------------------------------
  178. SIGT = ELAS EPE MOM1 MAM1 ;
  179. SI ITRAC ;
  180. TRAC SIGT MOM1 MAM1 DEF_6 CSU 15 'TITR'
  181. '[6] Contraintes, traction + temperature' ;
  182. FINSI ;
  183.  
  184.  
  185.  
  186.  
  187.  
  188.  
  189.  
  190. ************************************************************************
  191. * CHAPITRE 7 : MECANIQUE ELASTIQUE LINEAIRE *
  192. * CHARGEMENT THERMIQUE *
  193. * MATERIAU VARIABLE DANS L'ESPACE *
  194. ************************************************************************
  195.  
  196.  
  197. * ON REND LE COEFFICIENT DE DILATATION LINEAIRE ALPHA
  198. * VARIABLE EN FONCTION DE LA COORDONNEE X
  199. * ALPHA(X) = ALPHA0 * (1 + BETA*((X-XMOYEN)/(XMAX-XMIN))**3)
  200.  
  201. * RECUPERATION DE LA COORDONNEE X ET FORMATION D'UN
  202. * CHAMP PAR POINTS (CHPOINT) PORTANT CETTE INFORMATION
  203. * ----------------------------------------------------
  204. XX = COOR SU 1 ;
  205.  
  206.  
  207. * VALEURS EXTREMES ET MOYENNE DE LA COORDONNEE X DU MAILLAGE
  208. * ----------------------------------------------------------
  209. XMAX = MAXI XX ;
  210. XMIN = MINI XX ;
  211. XMOY = 0.5 * (XMAX + XMIN) ;
  212.  
  213.  
  214. * CHAMP PAR POINT SCALAIRE UNITAIRE
  215. * ---------------------------------
  216. XXU = MANU 'CHPO' 1 SU 'SCAL' 1. ;
  217.  
  218.  
  219. * CALCUL DU CHAMP PAR POINTS SCALAIRE DE ALPHA VARIABLE
  220. * -----------------------------------------------------
  221. BETA = 7. ;
  222. CHP_ALPH = ALPHAMAT *
  223. (1. + (BETA * (((XX - XMOY) / (XMAX - XMIN)) ** 3))) ;
  224. SI ITRAC ;
  225. TRAC CHP_ALPH SU 'TITR'
  226. '[7] Coefficient de dilatation thermique lineaire' ;
  227. FINSI ;
  228.  
  229.  
  230. * CONVERSION DU CHPOINT -> EN MCHAML
  231. * ----------------------------------
  232. CHM_ALPH = CHAN 'CHAM' CHP_ALPH MOM1 ;
  233.  
  234.  
  235. * MISE A JOUR DU MCHAML DU MATERIAU AVEC ALPHA VARIABLE
  236. * -----------------------------------------------------
  237. MAM1B = MATE MOM1 'YOUN' YOUNGMAT 'NU' NUMAT 'ALPH' CHM_ALPH ;
  238.  
  239.  
  240. * MISE A JOUR DE LA MATRICE DE RAIDEUR EN FONCTION DU NOUVEAU MCHAML
  241. * DE MATERIAU (MEME SI LA RAIDEUR EST INDEPENDANTE DE ALPHA !)
  242. * ------------------------------------------------------------------
  243. RI = RIGI MOM1 MAM1B ;
  244. RITOT = RI ET BLMX ET BLMY ;
  245.  
  246.  
  247. * MISE A JOUR DES DEFORMATIONS THERMIQUES PURES AVEC ALPHA VARIABLE
  248. * -----------------------------------------------------------------
  249. EPT = EPTH DELTA_TE MOM1 MAM1B ;
  250.  
  251.  
  252. * MISE A JOUR DES PSEUDO CONTRAINTES THERMIQUES AVEC ALPHA VARIABLE
  253. * -----------------------------------------------------------------
  254. SIT = ELAS EPT MOM1 MAM1B ;
  255.  
  256.  
  257. * MISE A JOUR DES FORCES NODALES POUR CETTE DEFORMATION THERMIQUE
  258. * ---------------------------------------------------------------
  259. FFT = BSIG SIT MOM1 MAM1B ;
  260.  
  261.  
  262. * DEPLACEMENTS AVEC ALPHA VARIABLE
  263. * --------------------------------
  264. U7 = RESO RITOT (TR ET FFT) ;
  265.  
  266.  
  267. * DEFORMEE AVEC ALPHA VARIABLE
  268. * ----------------------------
  269. DEF_7 = DEFO SU U7 150. ;
  270. DEF_7C = DEFO CSU U7 150. 'ROUG' ;
  271. SI ITRAC ;
  272. TRAC (DEF_INIC ET DEF_6C ET DEF_7C) 'TITR'
  273. '[7] Deformees, theromeca(O), + alpha var.(R)' ;
  274. FINSI ;
  275.  
  276.  
  277. * DEFORMATIONS TOTALES AVEC ALPHA VARIABLE
  278. * ----------------------------------------
  279. EP = EPSI U7 MOM1 ;
  280.  
  281.  
  282. * DEFORMATION ELASTIQUES AVEC ALPHA VARIABLE
  283. * ------------------------------------------
  284. EPE = EP - EPT ;
  285.  
  286.  
  287. * CONTRAINTES AVEC ALPHA VARIABLE
  288. * -------------------------------
  289. SIGT = ELAS EPE MOM1 MAM1B ;
  290. SI ITRAC ;
  291. TRAC SIGT MOM1 MAM1B DEF_7 CSU 15 'TITR'
  292. '[7] Contraintes, thermomeca, alpha var.' ;
  293. FINSI ;
  294.  
  295.  
  296.  
  297.  
  298.  
  299.  
  300.  
  301.  
  302. ************************************************************************
  303. * CHAPITRE 8 : MECANIQUE ELASTIQUE LINEAIRE ET PLASTIQUE PARFAIT *
  304. * THERMIQUE TRANSITOIRE PRECEDENT (COND. + CONV. + RAYO.) *
  305. * MATERIAU VARIABLE DANS L'ESPACE *
  306. ************************************************************************
  307.  
  308.  
  309. * MISE A JOUR DU MODELE MECANIQUE
  310. * -------------------------------
  311. MOM2 = MODE SU 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  312. 'PLASTIQUE' 'PARFAIT' ;
  313.  
  314.  
  315. * MISE A JOUR DU MATERIAU (AJOUT DE LA LIMITE D'ELASTICITE 'SIGY')
  316. * ----------------------------------------------------------------
  317. MAM2 = MATE MOM2 'YOUN' YOUNGMAT 'NU' NUMAT 'ALPH' CHM_ALPH
  318. 'SIGY' SIGYMAT ;
  319.  
  320.  
  321. * DEFINITION DU CHARGEMENT MECANIQUE INCREMENTAL
  322. * ----------------------------------------------
  323. * NOTA : EN PLASTICITE ON NE PEUT GENERALEMENT PAS APPLIQUER
  324. * BRUTALEMENT UN CHARGEMENT FINAL POUR LEQUEL IL Y AURAIT
  325. * DEPASSEMENT EXCESSIF DE LA LIMITE D'ELASTICITE DES LA
  326. * PREMIERE ITERATION DU CALCUL NON LINEAIRE ITERATIF
  327.  
  328. * EVOLUTION TEMPORELLE DU CHARGEMENT MECANIQUE
  329. * --------------------------------------------
  330. EVTM = EVOL 'MANU' 'Temps' (PROG 0. (0.2 * TPSFIN) TPSFIN)
  331. 'Coef' (PROG 0. 0. 1.) ;
  332. CHM = CHAR 'MECA' EVTM TR ;
  333.  
  334.  
  335. * PREPARATION DE LA TABLE POUR PASAPAS
  336. * ------------------------------------
  337. TAB2 = TABL 'PASAPAS' ;
  338. TAB2 . 'MODELE' = MOT ET MOC ET MOR ET MOM2 ;
  339. TAB2 . 'CARACTERISTIQUES' = MAT ET MAC ET MAR ET MAM2 ;
  340. TAB2 . 'BLOCAGES_MECANIQUES' = BLMX ET BLMY ;
  341. TAB2 . 'BLOCAGES_THERMIQUES' = BLT ;
  342. TAB2 . 'CHARGEMENT' = CHT ET CHACONV ET CHARAYE ET CHM;
  343. TAB2 . 'TEMPS_CALCULES' = PROG 0. 'PAS' (0.02 * TPSFIN) TPSFIN ;
  344. TAB2 . 'TEMPERATURES' = TABL ;
  345. TAB2 . 'TEMPERATURES' . 0 = T_INI ;
  346. TAB2 . 'CELSIUS' = VRAI ;
  347.  
  348.  
  349. * APPEL A PASAPAS
  350. * ---------------
  351. PASAPAS TAB2 ;
  352.  
  353.  
  354. * POST TRAITEMENT
  355. * ---------------
  356. * COMPARAISON DE LA DEFORMEE FINALE AVEC LES CALCULS PRECEDENTS
  357. N2 = DIME (TAB2 . 'DEPLACEMENTS') ;
  358. U8 = TAB2 . 'DEPLACEMENTS' . (N2 - 1) ;
  359. DEF_8 = DEFO SU U8 150. ;
  360. DEF_8C = DEFO CSU U8 150. 'VIOL' ;
  361. SI ITRAC ;
  362. TRAC (DEF_INIC ET DEF_6C ET DEF_7C ET DEF_8C) 'TITR'
  363. '[8] Deformees, thermomeca(O), + alpha var.(R), + plas.(V)' ;
  364. FINSI ;
  365.  
  366. * ON FAIT UNE PROCEDURE GENERIQUE QUI TRACE N'IMPORTE QUEL CHAMP PAR
  367. * ELEMENT DE LA TABLE DE PASAPAS A TOUS LES TEMPS SAUVES SUR LA DEFORMEE
  368.  
  369. DEBP @PTT T1*'TABLE' MOT1*'MOT' MOT2/'MOT' LIS1/'LISTREEL' ;
  370. NPAS = DIME (T1 . 'TEMPS') ;
  371. * MODELE MECANIQUE, CARACTERISTIQUES ET MAILLAGE
  372. MO = EXTR (T1 . 'MODELE') 'FORM' 'MECANIQUE' ;
  373. MA = REDU (T1 . 'CARACTERISTIQUES') MO ;
  374. MAIL1 = EXTR MO 'MAIL' ;
  375. CONT1 = CONT MAIL1 ;
  376. * BOUCLE SUR LES PAS DE TEMPS
  377. REPE B1 NPAS ;
  378. * VALEUR DU PSEUDO TEMPS
  379. TPS1 = T1 . 'TEMPS' . (&B1 - 1) ;
  380. * DEFORMEE AU TEMPS I
  381. DEPL1 = T1 . 'DEPLACEMENTS' . (&B1 - 1) ;
  382. DEF1 = DEFO MAIL1 DEPL1 150. ;
  383. * CHAMP A TRACER
  384. CHAM1 = T1 . MOT1 . (&B1 - 1) ;
  385. * TITRE AUTOMATIQUE
  386. TIT1 = CHAI MOT1 ;
  387. SI (NEG (TYPE MOT2) 'ANNULE') ;
  388. CHAM1 = EXCO CHAM1 MOT2 ;
  389. TIT1 = CHAI TIT1 ', composante ' MOT2 ;
  390. FINSI ;
  391. TIT1 = CHAI TIT1 ', au temps ' TPS1 ;
  392. * TRACE DU CHAMP SUR LA DEFORMEE
  393. SI (EGA (TYPE LIS1) 'ANNULE') ;
  394. TRAC CHAM1 MO MA DEF1 CONT1 15 'TITR' TIT1 ;
  395. SINON ;
  396. TRAC CHAM1 MO MA DEF1 CONT1 'TITR' TIT1 LIS1 ;
  397. FINSI ;
  398. FIN B1 ;
  399. FINP ;
  400.  
  401. * CONTRAINTES ET DEFORMATIONS PLASTIQUES CUMULEES
  402. SI ITRAC ;
  403. @PTT TAB2 'CONTRAINTES' (PROG 0. 'PAS' 10.E6 160.E6) ;
  404. FINSI ;
  405. VI = TAB2 . 'VARIABLES_INTERNES' . (N2 - 1) ;
  406. EQ_MAX = MAXI (EXCO 'EPSE' VI) ;
  407. SI ITRAC ;
  408. @PTT TAB2 'VARIABLES_INTERNES' 'EPSE'
  409. (PROG 0. 'PAS' (EQ_MAX / 15.) EQ_MAX) ;
  410. FINSI ;
  411.  
  412.  
  413.  
  414.  
  415.  
  416.  
  417.  
  418.  
  419.  
  420. ************************************************************************
  421. * CHAPITRE 9 : MECANIQUE ELASTIQUE LINEAIRE ET PLASTIQUE PARFAIT *
  422. * CHARGEMENT THERMIQUE *
  423. * MATERIAU VARIABLE DANS L'ESPACE *
  424. * MATERIAU DEPENDANT DE LA TEMPERATURE *
  425. ************************************************************************
  426.  
  427.  
  428. * VARIATION DE LA LIMITE ELASTIQUE EN FONCTION DE T
  429. * -------------------------------------------------
  430. EVSIGY = EVOL 'MANU' 'T' (PROG 50. 200.) 'SIGY' (PROG 250.E6 120.E6) ;
  431.  
  432.  
  433. * MISE A JOUR DU MATERIAU
  434. * -----------------------
  435. MAM2 = MATE MOM2 'YOUN' YOUNGMAT 'NU' NUMAT 'ALPH' CHM_ALPH
  436. 'SIGY' EVSIGY ;
  437.  
  438.  
  439. * MISE A JOUR DES DONNEES D'ENTREE DANS LA TABLE DE PASAPAS
  440. * ----------------------------------------------------------
  441. TAB2 = TABL 'PASAPAS' ;
  442. TAB2 . 'MODELE' = MOT ET MOC ET MOR ET MOM2 ;
  443. TAB2 . 'CARACTERISTIQUES' = MAT ET MAC ET MAR ET MAM2 ;
  444. TAB2 . 'BLOCAGES_MECANIQUES' = BLMX ET BLMY ;
  445. TAB2 . 'BLOCAGES_THERMIQUES' = BLT ;
  446. TAB2 . 'CHARGEMENT' = CHT ET CHACONV ET CHARAYE ET CHM;
  447. TAB2 . 'TEMPS_CALCULES' = PROG 0. 'PAS' (0.02 * TPSFIN) TPSFIN ;
  448. TAB2 . 'TEMPERATURES' = TABL ;
  449. TAB2 . 'TEMPERATURES' . 0 = T_INI ;
  450. TAB2 . 'CELSIUS' = VRAI ;
  451.  
  452.  
  453. * APPEL A PASAPAS
  454. * ---------------
  455. PASAPAS TAB2 ;
  456.  
  457.  
  458. * POST TRAITEMENT
  459. * ---------------
  460. * COMPARAISON DE LA DEFORMEE FINALE AVEC LES CALCULS PRECEDENTS
  461. N2 = DIME (TAB2 . 'DEPLACEMENTS') ;
  462. U9 = TAB2 . 'DEPLACEMENTS' . (N2 - 1) ;
  463. DEF_9 = DEFO SU U9 150. ;
  464. DEF_9C = DEFO CSU U9 150. 'BLEU' ;
  465. MOT1 = CHAI '[9] Deformees, thermomeca(O), + alpha var.(R), '
  466. '+ plas.(V), + Sy var.(B)' ;
  467. SI ITRAC ;
  468. TRAC (DEF_INIC ET DEF_6C ET DEF_7C ET DEF_8C ET DEF_9C) 'TITR' MOT1 ;
  469. FINSI ;
  470.  
  471. * EFFORTS DE REACTION AUX APPUIS
  472. REAC9 = TAB2 . 'REACTIONS' . (N2 - 1) ;
  473. VREAC = VECT REAC9 'FORC' 'ROUG' ;
  474. SI ITRAC ;
  475. TRAC (VFIMP ET VREAC) CSU 'TITR'
  476. '[9] Efforts imposes(V) et reactions aux appuis(R)' ;
  477. FINSI ;
  478.  
  479. * CONTRAINTES ET DEFORMATIONS PLASTIQUES CUMULEES
  480. SI ITRAC ;
  481. @PTT TAB2 'CONTRAINTES' (PROG 0. 'PAS' 10.E6 160.E6) ;
  482. FINSI ;
  483. VI = TAB2 . 'VARIABLES_INTERNES' . (N2 - 1) ;
  484. EQ_MAX = MAXI (EXCO 'EPSE' VI) ;
  485. SI ITRAC ;
  486. @PTT TAB2 'VARIABLES_INTERNES' 'EPSE'
  487. (PROG 0. 'PAS' (EQ_MAX / 15.) EQ_MAX) ;
  488. FINSI ;
  489.  
  490.  
  491.  
  492.  
  493.  
  494.  
  495. FIN ;
  496.  
  497.  
  498.  
  499.  
  500.  
  501.  
  502.  

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