Télécharger idvar5.eso

Retour à la liste

Numérotation des lignes :

idvar5
  1. C IDVAR5 SOURCE BR232186 16/12/05 21:39:48 9238
  2.  
  3. *_______________________________________________________________________
  4. * RECHERCHE DES NOMS DE COMPOSANTES DES VARIABLES INTERNES
  5. * --------------------------------------------------------
  6. * CAS DE L'ENDOMMAGEMENT
  7. * ---------------------
  8. *
  9. * ENTREES :
  10. * ---------
  11. * MATEPL NUMERO DU MATERIAU ENDOMAGEABLE
  12. * MFR NUMERO DE FORMULATION
  13. * IFOUR VALEUR DE IFOU DANS CCOPTIO
  14. * IFLAG 0=TOUTES LES VARIABLES
  15. * 1=FONCTIONS SEUIL
  16. *
  17. * SORTIES :
  18. * ---------
  19. * IPCOMP POINTEUR SUR LA LISTE DES NOMS DE COMPOSANTES
  20. * OBLIGATOIRES ET FACULTATIVES
  21. * NBROBL NOMBRE DE COMPOSANTES OBLIGATOIRES
  22. * NBRFAC NOMBRE DE COMPOSANTES FACTULTATIVES
  23. *
  24. * PP 16/8/91 ajout var.int pour l'endommagement
  25. *_______________________________________________________________________
  26.  
  27. SUBROUTINE IDVAR5 (MATEPL,MFR,IFOUR,IFLAG,IPCOMP,NBROBL,NBRFAC)
  28.  
  29. IMPLICIT INTEGER(I-N)
  30. IMPLICIT REAL*8 (A-H,O-Z)
  31.  
  32. SEGMENT NOMID
  33. CHARACTER*8 LESOBL(NBROBL),LESFAC(NBRFAC)
  34. ENDSEGMENT
  35.  
  36. NBROBL=0
  37. NBRFAC=0
  38. NOMID=0
  39. C PRINt *,'ENTREE IDVAR5 MATEPL, MFR =',MATEPL,MFR
  40.  
  41. C 1 - Cas du modele SCALAIRE
  42. C ============================
  43. IF (MATEPL.EQ.30.OR.MATEPL.EQ.118) THEN
  44. IF (IFLAG.EQ.0) THEN
  45. NBROBL=2
  46. SEGINI,NOMID
  47. LESOBL(1)='EPTI'
  48. LESOBL(2)='D'
  49. ELSE
  50. NBROBL=1
  51. SEGINI,NOMID
  52. LESOBL(1)='EPTI'
  53. ENDIF
  54.  
  55. C 2 - Cas du modele UNILATERAL simple ou avec viscosite
  56. C =======================================================
  57. ELSE IF (MATEPL.EQ.31.OR.MATEPL.EQ. 97) THEN
  58. IF (IFLAG.EQ.0) THEN
  59. NBROBL=4
  60. SEGINI,NOMID
  61. LESOBL(1)='YLI1'
  62. LESOBL(2)='YLI2'
  63. LESOBL(3)='DOM1'
  64. LESOBL(4)='DOM2'
  65. ELSE
  66. NBROBL=2
  67. SEGINI,NOMID
  68. LESOBL(1)='YLI1'
  69. LESOBL(2)='YLI2'
  70. ENDIF
  71.  
  72. C 3 - Cas du modele ROTATING_CRACK
  73. C ==================================
  74. ELSE IF (MATEPL.EQ.37) THEN
  75. NBROBL=2
  76. SEGINI,NOMID
  77. LESOBL(1)='EPSM'
  78. LESOBL(2)='EPSN'
  79.  
  80. C 4 - Cas du modele SIC_SIC
  81. C Modele disponible seulement en MASSIF 3D
  82. C ==============================================
  83. ELSE IF (MATEPL.EQ.88) THEN
  84. IF (MFR.EQ.1) THEN
  85. NBROBL=10
  86. SEGINI,NOMID
  87. LESOBL( 1)='EPSE'
  88. LESOBL( 2)='DOM1'
  89. LESOBL( 3)='DOM2'
  90. LESOBL( 4)='DOM3'
  91. LESOBL( 5)='DF11'
  92. LESOBL( 6)='DF22'
  93. LESOBL( 7)='DF33'
  94. LESOBL( 8)='DF23'
  95. LESOBL( 9)='DF13'
  96. LESOBL(10)='DF12'
  97. ELSE
  98. CALL ERREUR (19)
  99. ENDIF
  100.  
  101. C 5 - Cas du modele MICROPLANS ORTHOTROPE
  102. C =========================================
  103. ELSE IF (MATEPL.EQ.96) THEN
  104. NBROBL=14
  105. SEGINI,NOMID
  106. LESOBL( 1)='EPSE'
  107. LESOBL( 2)='R '
  108. C= Tenseur d'endommagement
  109. LESOBL( 3)='D11 '
  110. LESOBL( 4)='D22 '
  111. LESOBL( 5)='D33 '
  112. LESOBL( 6)='D12 '
  113. LESOBL( 7)='D13 '
  114. LESOBL( 8)='D23 '
  115. C= Tenseur des contraintes anelastiques dues a l'endommagement
  116. LESOBL( 9)='SA11'
  117. LESOBL(10)='SA22'
  118. LESOBL(11)='SA33'
  119. LESOBL(12)='SA12'
  120. LESOBL(13)='SA13'
  121. LESOBL(14)='SA23'
  122.  
  123. C 5 - Cas du modele MICROPLANS ISOTROPE
  124. C =======================================
  125. ELSE IF (MATEPL.EQ.98) THEN
  126. NBROBL=9
  127. SEGINI,NOMID
  128. LESOBL(1)='EPSE'
  129. LESOBL(2)='R '
  130. C= Scalaire d'endommagement
  131. LESOBL(3)='DOM'
  132. C= Tenseur des contraintes anelastiques dues a l'endommagement
  133. LESOBL(4)='SA11'
  134. LESOBL(5)='SA22'
  135. LESOBL(6)='SA33'
  136. LESOBL(7)='SA12'
  137. LESOBL(8)='SA13'
  138. LESOBL(9)='SA23'
  139.  
  140. * RDj
  141. *
  142. ELSEIF(MATEPL.EQ.124)THEN
  143. *
  144. * CAS DU MODELE SCALAIRE AVEC DEFORMATION ANELASTIQUE
  145. *
  146. NBROBL=24
  147. SEGINI,NOMID
  148. LESOBL(1) ='ECM'
  149. LESOBL(2) ='D'
  150. LESOBL(3) ='DELT'
  151. LESOBL(4) ='EIC1'
  152. LESOBL(5) ='EIC2'
  153. LESOBL(6) ='EIC3'
  154. LESOBL(7) ='EIC4'
  155. LESOBL(8) ='EIC5'
  156. LESOBL(9) ='EIC6'
  157. LESOBL(10) ='EIT1'
  158. LESOBL(11) ='EIT2'
  159. LESOBL(12) ='EIT3'
  160. LESOBL(13) ='EIT4'
  161. LESOBL(14) ='EIT5'
  162. LESOBL(15) ='EIT6'
  163. LESOBL(16) ='EPTI'
  164. LESOBL(17) ='VER1'
  165. LESOBL(18) ='VER2'
  166. LESOBL(19) ='VER3'
  167. LESOBL(20) ='VER4'
  168. LESOBL(21) ='VER5'
  169. LESOBL(22) ='VER6'
  170. LESOBL(23) ='ENDP'
  171. LESOBL(24) ='ENDM'
  172. * RDj
  173. *
  174. * CCd
  175. *
  176. ELSEIF(MATEPL.EQ.125)THEN
  177. *
  178. * MODEL DAMAGE_TC - R. Faria & J. Oliver
  179. *
  180. C NBROBL=12
  181. NBROBL=6
  182. SEGINI NOMID
  183. LESOBL(1)='THRT'
  184. LESOBL(2)='EQUT'
  185. LESOBL(3)='DAMT'
  186. LESOBL(4)='THRC'
  187. LESOBL(5)='EQUC'
  188. LESOBL(6)='DAMC'
  189. *
  190. C LESOBL(7) ='SE11'
  191. C LESOBL(8) ='SE22'
  192. C LESOBL(9) ='SE33'
  193. C LESOBL(10)='SE12'
  194. C LESOBL(11)='SE13'
  195. C LESOBL(12)='SE23'
  196. *
  197.  
  198. C 11 - Cas du modele SICf_SIC Scalaire ONERA
  199. C SICSCAL
  200. C Modele disponible seulement en MASSIF 3D
  201. C ==============================================
  202. ELSE IF (MATEPL.EQ.131) THEN
  203. IF (MFR.EQ.1) THEN
  204. NBROBL=19
  205. SEGINI,NOMID
  206. * Somme des defo inel
  207. LESOBL( 1)='EPSE'
  208. * Dommage dans les différentes directions
  209. LESOBL( 2)='DOM1'
  210. LESOBL( 3)='DOM2'
  211. LESOBL( 4)='DOM3'
  212. * Déformation repère orthotrope ordre ONERA
  213. LESOBL( 5)='DF11'
  214. LESOBL( 6)='DF22'
  215. LESOBL( 7)='DF33'
  216. LESOBL( 8)='DF23'
  217. LESOBL( 9)='DF13'
  218. LESOBL(10)='DF12'
  219. * Indice d'activation des fissures
  220. LESOBL(11)='NUA1'
  221. LESOBL(12)='NUA2'
  222. LESOBL(13)='NUA3'
  223. * Déformations inélastiques repère orthotrope ordre ONERA
  224. LESOBL(14)='EPI1'
  225. LESOBL(15)='EPI2'
  226. LESOBL(16)='EPI3'
  227. LESOBL(17)='EPI4'
  228. LESOBL(18)='EPI5'
  229. LESOBL(19)='EPI6'
  230. ELSE
  231. CALL ERREUR (19)
  232. ENDIF
  233. C 12 - Cas du modele SICf_SIC pseudo Tensoriel ONERA
  234. C SICTENS
  235. C Modele disponible seulement en MASSIF 3D
  236. C ==============================================
  237. ELSE IF (MATEPL.EQ.132) THEN
  238. IF (MFR.EQ.1) THEN
  239. NBROBL=23
  240. SEGINI,NOMID
  241. * somme defo inel
  242. LESOBL( 1)='EPSE'
  243. * 5 variables de dommage
  244. LESOBL( 2)='DOM1'
  245. LESOBL( 3)='DOM2'
  246. LESOBL( 4)='DOM3'
  247. LESOBL( 5)='DOM4'
  248. LESOBL( 6)='DOM5'
  249. * Défo dans le repère ortho ordre ONERA
  250. LESOBL( 7)='DF11'
  251. LESOBL( 8)='DF22'
  252. LESOBL( 9)='DF33'
  253. LESOBL( 10)='DF23'
  254. LESOBL( 11)='DF13'
  255. LESOBL( 12)='DF12'
  256. * Indices d'activation
  257. LESOBL(13)='NUA1'
  258. LESOBL(14)='NUA2'
  259. LESOBL(15)='NUA3'
  260. LESOBL(16)='NUA4'
  261. LESOBL(17)='NUA5'
  262. * Déformations inélastique repère ortho ordre ONERA
  263. LESOBL(18)='EPI1'
  264. LESOBL(19)='EPI2'
  265. LESOBL(20)='EPI3'
  266. LESOBL(21)='EPI4'
  267. LESOBL(22)='EPI5'
  268. LESOBL(23)='EPI6'
  269. ELSE
  270. CALL ERREUR (19)
  271. ENDIF
  272. *
  273. ELSEIF(MATEPL.EQ.134)THEN
  274. *
  275. * 13. CAS DU MODELE DESMORAT
  276. *
  277. NBROBL=23
  278. SEGINI,NOMID
  279. *
  280. * flag pour le non local
  281. LESOBL(1) ='EPTI'
  282. *
  283. * tenseur d'endommagement
  284. *
  285. LESOBL(2) ='D11 '
  286. LESOBL(3) ='D22 '
  287. LESOBL(4) ='D33 '
  288. LESOBL(5) ='D12 '
  289. LESOBL(6) ='D13 '
  290. LESOBL(7) ='D23 '
  291. *
  292. * tenseur de la base d'endommagement
  293. *
  294. LESOBL(8) ='RT11'
  295. LESOBL(9) ='RT12'
  296. LESOBL(10) ='RT13'
  297. LESOBL(11) ='RT21'
  298. LESOBL(12) ='RT22'
  299. LESOBL(13) ='RT23'
  300. LESOBL(14) ='RT31'
  301. LESOBL(15) ='RT32'
  302. LESOBL(16) ='RT33'
  303. *
  304. * flag pour l'endommagement maximum
  305. LESOBL(17) ='FAIL'
  306. *
  307. * tenseur des deformations
  308. *
  309. LESOBL(18) ='EP11'
  310. LESOBL(19) ='EP22'
  311. LESOBL(20) ='EP33'
  312. LESOBL(21) ='EP12'
  313. LESOBL(22) ='EP13'
  314. LESOBL(23) ='EP23'
  315. ELSEIF(MATEPL.EQ.141)THEN
  316. *
  317. C 14 - Cas du modele 'FATSIN'
  318. C ==============================================
  319. IF (IFLAG.EQ.0) THEN
  320. NBROBL=3
  321. SEGINI,NOMID
  322. LESOBL(1)='EPTI'
  323. LESOBL(2) ='D'
  324. LESOBL(3) ='EPT0'
  325. ELSE
  326. NBROBL=1
  327. SEGINI,NOMID
  328. LESOBL(1)='EPTI'
  329. ENDIF
  330. ELSEIF(MATEPL.EQ.144)THEN
  331. *
  332. C 15 - Cas du modele 'RICRAG'
  333. C ==============================================
  334. NBROBL=26
  335. SEGINI,NOMID
  336. *
  337. * POUR LE NON LOCAL
  338. *
  339. LESOBL(1) ='EPTI'
  340. *
  341. * TENSEUR DE GLISSEMENT
  342. *
  343. LESOBL(2) ='EPI1'
  344. LESOBL(3) ='EPI2'
  345. LESOBL(4) ='EPI3'
  346. LESOBL(5) ='EPI4'
  347. LESOBL(6) ='EPI5'
  348. LESOBL(7) ='EPI6'
  349. *
  350. * TENSEUR D ECROUISSAGE CINEMATIQUE
  351. *
  352. LESOBL(8) ='ABC1'
  353. LESOBL(9) ='ABC2'
  354. LESOBL(10) ='ABC3'
  355. LESOBL(11) ='ABC4'
  356. LESOBL(12) ='ABC5'
  357. LESOBL(13) ='ABC6'
  358. *
  359. * QUELQUES VARIABLES INTERNES SCALAIRES
  360. *
  361. LESOBL(14) ='D '
  362. LESOBL(15) ='DOND'
  363. LESOBL(16) ='DONI'
  364. LESOBL(17) ='FFFZ'
  365. LESOBL(18) ='TYPC'
  366. *
  367. * TENSEUR DE DEFORMATION TOTALES
  368. *
  369. LESOBL(19) ='EPA1'
  370. LESOBL(20) ='EPA2'
  371. LESOBL(21) ='EPA3'
  372. LESOBL(22) ='EPA4'
  373. LESOBL(23) ='EPA5'
  374. LESOBL(24) ='EPA6'
  375. *
  376. * FLAGS POUR ENDOMMAGEMENT ECCESSIF ET ECROUISSAGE
  377. *
  378. LESOBL(25) ='FZF1'
  379. LESOBL(26) ='FZF2'
  380. C +BR
  381. ELSEIF(MATEPL.EQ.157)THEN
  382. *
  383. C 16 - Cas du modele 'GLRC_DM'
  384. C ==============================================
  385. NBROBL=12
  386. SEGINI,NOMID
  387. LESOBL( 1) ='VA01'
  388. LESOBL( 2) ='VA02'
  389. LESOBL( 3) ='VA03'
  390. LESOBL( 4) ='VA04'
  391. LESOBL( 5) ='VA05'
  392. LESOBL( 6) ='VA06'
  393. LESOBL( 7) ='VA07'
  394. LESOBL( 8) ='VA08'
  395. LESOBL( 9) ='XD1 '
  396. LESOBL(10) ='XD2 '
  397. LESOBL(11) ='XDM '
  398. LESOBL(12) ='XDF '
  399.  
  400. ELSEIF(MATEPL.EQ.158)THEN
  401. *
  402. C 17 - Cas du modele 'RICBET'
  403. C ==============================================
  404. NBROBL=42
  405. SEGINI,NOMID
  406.  
  407. LESOBL( 1)='EPTI'
  408.  
  409. LESOBL( 2)='EPI1'
  410. LESOBL( 3)='EPI2'
  411. LESOBL( 4)='EPI3'
  412. LESOBL( 5)='EPI4'
  413. LESOBL( 6)='EPI5'
  414. LESOBL( 7)='EPI6'
  415.  
  416. LESOBL( 8)='ABC1'
  417. LESOBL( 9)='ABC2'
  418. LESOBL(10)='ABC3'
  419. LESOBL(11)='ABC4'
  420. LESOBL(12)='ABC5'
  421. LESOBL(13)='ABC6'
  422.  
  423. LESOBL(14)='D '
  424.  
  425. LESOBL(15)='FZFT'
  426.  
  427. LESOBL(16)='EPA1'
  428. LESOBL(17)='EPA2'
  429. LESOBL(18)='EPA3'
  430. LESOBL(19)='EPA4'
  431. LESOBL(20)='EPA5'
  432. LESOBL(21)='EPA6'
  433.  
  434. LESOBL(22)='SIC1'
  435. LESOBL(23)='SIC2'
  436. LESOBL(24)='SIC3'
  437. LESOBL(25)='SIC4'
  438. LESOBL(26)='SIC5'
  439. LESOBL(27)='SIC6'
  440.  
  441. LESOBL(28)='DEFF'
  442.  
  443. LESOBL(29)='SIP1'
  444. LESOBL(30)='SIP2'
  445. LESOBL(31)='SIP3'
  446. LESOBL(32)='SIP4'
  447. LESOBL(33)='SIP5'
  448. LESOBL(34)='SIP6'
  449.  
  450. LESOBL(35)='EPP1'
  451. LESOBL(36)='EPP2'
  452. LESOBL(37)='EPP3'
  453. LESOBL(38)='EPP4'
  454. LESOBL(39)='EPP5'
  455. LESOBL(40)='EPP6'
  456.  
  457. LESOBL(41)='RECO'
  458.  
  459. LESOBL(42)='DEPP'
  460.  
  461. ELSEIF(MATEPL.EQ.159)THEN
  462. *
  463. C 18 - Cas du modele 'RICCOQ'
  464. C ==============================================
  465. NBROBL=18
  466. SEGINI,NOMID
  467.  
  468. LESOBL(1) ='EPM1'
  469. LESOBL(2) ='EPM2'
  470. LESOBL(3) ='EPM3'
  471.  
  472. LESOBL(4) ='EPF1'
  473. LESOBL(5) ='EPF2'
  474. LESOBL(6) ='EPF3'
  475.  
  476. LESOBL(7) ='DT '
  477. LESOBL(8) ='ZT '
  478. LESOBL(9) ='DC '
  479. LESOBL(10)='ZC '
  480.  
  481. LESOBL(11)='EPP1'
  482. LESOBL(12)='EPP2'
  483. LESOBL(13)='EPP3'
  484.  
  485. LESOBL(14)='XEC1'
  486. LESOBL(15)='XEC2'
  487. LESOBL(16)='XEC3'
  488.  
  489. LESOBL(17)='XYDT'
  490. LESOBL(18)='XYDC'
  491.  
  492. ELSEIF(MATEPL.EQ.163)THEN
  493. *
  494. C 19 - Cas du modele 'CONCYC'
  495. C ==============================================
  496. NBROBL=35
  497. SEGINI,NOMID
  498. LESOBL( 1) ='EP11'
  499. LESOBL( 2) ='EP22'
  500. LESOBL( 3) ='EP33'
  501. LESOBL( 4) ='EP12'
  502. LESOBL( 5) ='EP13'
  503. LESOBL( 6) ='EP23'
  504.  
  505. LESOBL( 7) ='SF11'
  506. LESOBL( 8) ='SF22'
  507. LESOBL( 9) ='SF33'
  508. LESOBL(10) ='SF12'
  509. LESOBL(11) ='SF13'
  510. LESOBL(12) ='SF23'
  511.  
  512. LESOBL(13) ='QF11'
  513. LESOBL(14) ='QF22'
  514. LESOBL(15) ='QF33'
  515. LESOBL(16) ='QF12'
  516. LESOBL(17) ='QF13'
  517. LESOBL(18) ='QF23'
  518.  
  519. LESOBL(19) ='D '
  520. LESOBL(20) ='ZT '
  521. LESOBL(21) ='PHI '
  522.  
  523. LESOBL(22) ='EM11'
  524. LESOBL(23) ='EM22'
  525. LESOBL(24) ='EM33'
  526. LESOBL(25) ='EM12'
  527. LESOBL(26) ='EM13'
  528. LESOBL(27) ='EM23'
  529.  
  530. LESOBL(28) ='DP11'
  531. LESOBL(29) ='DP22'
  532. LESOBL(30) ='DP33'
  533. LESOBL(31) ='DP12'
  534. LESOBL(32) ='DP13'
  535. LESOBL(33) ='DP23'
  536.  
  537. LESOBL(34) ='EDIS'
  538.  
  539. LESOBL(35) ='EPTI'
  540.  
  541. C -BR
  542. ELSEIF(MATEPL.EQ.177)THEN
  543.  
  544. C 20 - Cas du modele 'EFEM'
  545. C ==============================================
  546. NBROBL=59
  547. SEGINI,NOMID
  548.  
  549. LESOBL(1) ='RHO1'
  550. LESOBL(2) ='RHO2'
  551. LESOBL(3) ='RHO4'
  552. LESOBL(4) ='RHO7'
  553.  
  554. LESOBL(5) ='EPS1'
  555. LESOBL(6) ='EPS2'
  556. LESOBL(7) ='EPS3'
  557. LESOBL(8) ='EPS4'
  558. LESOBL(9) ='EPS5'
  559. LESOBL(10) ='EPS6'
  560.  
  561. LESOBL(11) ='T1 '
  562. LESOBL(12) ='T2 '
  563.  
  564. LESOBL(13) ='OUV1'
  565. LESOBL(14) ='OUV2'
  566.  
  567. LESOBL(15) ='Z1 '
  568. LESOBL(16) ='Z2 '
  569. LESOBL(17) ='Z4 '
  570. LESOBL(18) ='Z7 '
  571.  
  572. LESOBL(19) ='EPSZ'
  573.  
  574. LESOBL(20) ='MNOU'
  575.  
  576. LESOBL(21) ='ZIN1'
  577. LESOBL(22) ='ZIN2'
  578. LESOBL(23) ='ZIN4'
  579. LESOBL(24) ='ZIN7'
  580.  
  581. LESOBL(25) ='DGL1'
  582. LESOBL(26) ='DGL2'
  583. LESOBL(27) ='DGL3'
  584.  
  585. LESOBL(28) ='SGL1'
  586. LESOBL(29) ='SGL2'
  587. LESOBL(30) ='SGL3'
  588.  
  589. LESOBL(31) ='NFLA'
  590. LESOBL(32) ='NIT '
  591.  
  592. LESOBL(33) ='XNX1'
  593. LESOBL(34) ='YN '
  594.  
  595. LESOBL(35) ='CPLA'
  596. LESOBL(36) ='Z_U '
  597. LESOBL(37) ='UMAX'
  598. LESOBL(38) ='OUVN'
  599. LESOBL(39) ='OUVT'
  600. LESOBL(40) ='RESI'
  601. LESOBL(41) ='XNT '
  602. LESOBL(42) ='YNT '
  603. LESOBL(43) ='TUN '
  604. LESOBL(44) ='TUT '
  605. LESOBL(45) ='XNOR'
  606. LESOBL(46) ='YNOR'
  607. LESOBL(47) ='XPSI'
  608. LESOBL(48) ='VPT1'
  609. LESOBL(49) ='VPT2'
  610. LESOBL(50) ='XNVP'
  611. C-------POUR LE NON LOCAL
  612.  
  613. LESOBL(51) ='ENL1'
  614. LESOBL(52) ='ENL2'
  615. LESOBL(53) ='ENL3'
  616. LESOBL(54) ='SXX'
  617. LESOBL(55) ='SYY'
  618. LESOBL(56) ='SZZ'
  619. LESOBL(57) ='SXY'
  620. LESOBL(58) ='SLM'
  621. LESOBL(59) ='DIFF'
  622.  
  623.  
  624. ENDIF
  625. IPCOMP=NOMID
  626. SEGDES,NOMID
  627.  
  628.  
  629. RETURN
  630. END
  631.  
  632.  
  633.  
  634.  
  635.  
  636.  
  637.  
  638.  
  639.  
  640.  
  641.  

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