Télécharger devlfa.eso

Retour à la liste

Numérotation des lignes :

  1. C DEVLFA SOURCE BP208322 18/12/20 21:15:35 10048
  2. SUBROUTINE DEVLFA(Q1,Q2,FTOTA,NA1,IPALA,IPLIA,XPALA,XVALA,
  3. & NLIAA,PDT,T,NPAS,IND,FINERT,IVINIT,FTEST,FTOTA0)
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. *--------------------------------------------------------------------*
  7. * *
  8. * Operateur DYNE : algorithme de Fu - de Vogelaere *
  9. * ________________________________________________ *
  10. * *
  11. * Calcul des forces de choc base A. *
  12. * *
  13. * Parametres: *
  14. * *
  15. c c * e Q1(.,.) Vecteur des deplacements generalises. *
  16. * e Q2(.,.) Vecteur des vitesses generalisees. *
  17. * es FTOTA Forces exterieures totalisees sur la base A. *
  18. * e NA1 Nombre total d'inconnues en base A. *
  19. * e IPALA Renseigne sur le type de la liaison. *
  20. * e IPLIA Tableau contenant les numeros "DYNE" de la liaison. *
  21. * e XPALA Tableau contenant les parametres de la liaison. *
  22. * es XVALA Tableau contenant les variables internes des liaisons *
  23. * e NLIAA Nombre de liaisons sur la base A. *
  24. * e PDT pas de temps *
  25. * e T temps *
  26. * e NPAS Numero du pas de temps *
  27. * e IND Indice du demi-pas de temps *
  28. * = 2 si 1er demi-pas = 1 si 2eme demi-pas *
  29. * es FINERT Forces d'inertie *
  30. * e IVINIT =1 si vitesses initiales, =0 sinon *
  31. * *
  32. * Auteur, date de creation: *
  33. * *
  34. * Lionel VIVAN, le 20 aout 1989. *
  35. * E. de LANGRE, 15 fevrier 91 correction jeu negatif *
  36. * Lionel VIVAN, 12 mars 1991 ajout de la liaison FLUIDE *
  37. * Denis ROBERT, 30 avril 1992 ajout liaison POLYNOMIALE *
  38. * *
  39. * remarque: *
  40. * ========= *
  41. * Si jeu negatif (cas particulier de la base A ou il n'y a pas de *
  42. * normale), on renverse les variables avec XNORM. *
  43. * *
  44. *--------------------------------------------------------------------*
  45. *
  46. INTEGER IPALA(NLIAA,*),IPLIA(NLIAA,*)
  47. REAL*8 XPALA(NLIAA,*),Q1(NA1,*),Q2(NA1,*),FTOTA(NA1,*)
  48. REAL*8 XVALA(NLIAA,4,*),FINERT(NA1,*)
  49. PARAMETER (XZERO = 0.D0)
  50. REAL*8 FTest(nA1,4)
  51. REAL*8 ftOTa0(NA1,4)
  52.  
  53. XFIN = 0.D0
  54. PDTS2 = 0.5D0 * PDT
  55.  
  56. *--------------------------------------------------------------------*
  57. * BOUCLE SUR LES LIAISONS
  58. *--------------------------------------------------------------------*
  59.  
  60. DO 10 I = 1,NLIAA
  61.  
  62. ITYP = IPALA(I,1)
  63. icond= IPALA(I,2)
  64. iannul= 0
  65.  
  66. IF (ICOND .NE. 1 ) GOTO 199
  67. * >>> CAS DES LIAISONS CONDITIONNELLES <<<
  68.  
  69. * --> boucle sur j
  70. DO 101 j = 4,20
  71.  
  72. jliai = ipala(i,j)
  73. jpliai = abs ( jliai)
  74. if ( jliai . EQ. 0 ) goto 101
  75.  
  76. jtyp = ipala(jpliai,1)
  77. do 102 ik = 1,nA1
  78. do 102 jk = 1,4
  79. ftest(ik,jk) = 0d0
  80. ftota0 (ik,jk) = ftota(ik,jk)
  81. 102 continue
  82.  
  83. * 1> CALCUL DES LIAISONS POUR TEST
  84. *
  85. * ------ choc elementaire POINT_PLAN sans amortissement
  86. *
  87. IF (JTYP.EQ.1) THEN
  88. XRAID = XPALA(jpliai,1)
  89. XJEU = XPALA(jpliai,2)
  90. XNORM = 1.
  91. IF (XJEU.LT.0D0) THEN
  92. XNORM = -1.
  93. XJEU = -1*XJEU
  94. ENDIF
  95. INA1 = IPLIA(jpliai,1)
  96. XDEP = XNORM*Q1(INA1,IND)
  97. CALL DYCHEL(XDEP,XRAID,XJEU, XFL,iannul)
  98. XVALA(jpliai,IND,1) = XNORM*XFL
  99. XVALA(jpliai,IND,4) = XNORM*XDEP
  100. FTest(INA1,IND) =Ftest(INA1,IND) + XNORM*XFL
  101. *
  102. * ------ choc elementaire POINT_PLAN avec amortissement
  103. *
  104. ELSE IF (JTYP.EQ.2) THEN
  105. XRAID = XPALA(jpliai,1)
  106. XJEU = XPALA(jpliai,2)
  107. XAMO = XPALA(jpliai,3)
  108. XNORM = 1.
  109. IPERM = 0
  110. IF (XJEU.LT.XZERO) THEN
  111. XNORM = -1.
  112. XJEU = -1*XJEU
  113. ENDIF
  114. INA1 = IPLIA(jpliai,1)
  115. XDEP = XNORM * Q1(INA1,IND)
  116. IF ((NPAS.EQ.1).AND.(IND.EQ.3)) THEN
  117. XVIT = XNORM * Q2(INA1,IND)
  118. ELSE
  119. IND2 = IND + 1
  120. XDEPM1 = XNORM * Q1(INA1,IND2)
  121. XVIT = (XDEP - XDEPM1) / PDTS2
  122. ENDIF
  123. XVALA(jpliai,IND,3) = XNORM*XVIT
  124. CALL DYCHAM(XDEP,XVIT,XRAID,XJEU,XAMO,XFL,IPERM,iannul)
  125. *+*
  126. IF (XDEP.GT.0.D0 .AND. XFL.GT.0.D0) XFL = 0.D0
  127. IF (XDEP.LT.0.D0 .AND. XFL.LT.0.D0) XFL = 0.D0
  128. *+*
  129. XVALA(jpliai,IND,1) = XNORM*XFL
  130. XVALA(jpliai,IND,4) = XNORM*XDEP
  131. FTest(INA1,IND) =Ftest(INA1,IND) + XNORM*XFL
  132. *
  133. * ------ choc elementaire POINT_PLAN_FLUIDE
  134. *
  135. ELSE IF (JTYP.EQ.3) THEN
  136. XINER = XPALA(jpliai,1)
  137. XCONV = XPALA(jpliai,2)
  138. XVISC = XPALA(jpliai,3)
  139. XPCEL = XPALA(jpliai,4)
  140. XPCRA = XPALA(jpliai,5)
  141. XJEU = XPALA(jpliai,6)
  142. INA1 = IPLIA(jpliai,1)
  143. XDEP = Q1(INA1,IND)
  144. IF ((NPAS.EQ.1).AND.(IND.EQ.3)) THEN
  145. XVIT = Q2(INA1,IND)
  146. ELSE
  147. IND2 = IND + 1
  148. XDEPM1 = Q1(INA1,IND2)
  149. XVIT = (XDEP - XDEPM1) / PDTS2
  150. ENDIF
  151. IF (XJEU.GT.0D0) THEN
  152. XDH= XJEU - XDEP
  153. XNORM = 1.
  154. ELSE
  155. XDH= XDEP - XJEU
  156. XNORM = -1.
  157. ENDIF
  158. * Calcul de la masse ajoutee
  159. XXIN = -XINER / XDH
  160. FINERT(INA1,IND) = FINERT(INA1,IND) + XXIN
  161. * Calcul de la force de convection
  162. CALL DYFCON(XDH,XDEP,XVIT,XJEU,XCONV,XFCO,iannul)
  163. * Calcul de la force de viscosite
  164. CALL DYFVIS(XDH,XDEP,XVIT,XJEU,XVISC,XFVI,iannul)
  165. * Calcul de la force de perte de charge
  166. CALL DYFPDC(XDH,XDEP,XVIT,XJEU,XPCEL,XPCRA,XFPE,iannul)
  167. XFL = (XFCO* XNORM ) + XFVI + XFPE
  168. XVALA(jpliai,IND,1) = XDEP
  169. XVALA(jpliai,IND,2) = XVIT
  170. XVALA(jpliai,IND,3) = XXIN
  171. XVALA(jpliai,IND,4) = XFCO*XNORM
  172. XVALA(jpliai,IND,5) = XFVI
  173. XVALA(jpliai,IND,6) = XFPE
  174. FTest(INA1,IND) =Ftest(iNA1,IND) + XFL
  175. *
  176. * ------ force elementaire de COUPLAGE EN VITESSE
  177. *
  178. ELSE IF (JTYP.EQ.4) THEN
  179. INA1 = IPLIA(jpliai,1)
  180. INA2 = IPLIA(jpliai,2)
  181. XDEP = Q1(INA2,IND)
  182. IF ((NPAS.EQ.1).AND.(IND.EQ.3)) THEN
  183. XVIT = Q2(INA2,IND)
  184. ELSE
  185. IND2 = IND + 1
  186. XDEPM1 = Q1(INA2,IND2)
  187. XVIT = (XDEP - XDEPM1) / PDTS2
  188. ENDIF
  189. XCPLGE = XPALA(jpliai,1)
  190. CALL DYCPLV(XVIT,XCPLGE,XFL,iannul)
  191. XVALA(jpliai,IND,3) = XVIT
  192. XVALA(jpliai,IND,1) = XFL
  193. XVALA(jpliai,IND,4) = XDEP
  194. FTest(INA1,IND) =Ftest(INA1,IND) + XFL
  195. *
  196. * ------ force elementaire de COUPLAGE EN DEPLACEMENT
  197. *
  198. ELSE IF (JTYP.EQ.5) THEN
  199. INA1 = IPLIA(jpliai,1)
  200. INA2 = IPLIA(jpliai,2)
  201. XDEP = Q1(INA2,IND)
  202. XCPLGE = XPALA(jpliai,1)
  203. jfonct = ipala(jpliai,3)
  204. cbp calcul eventuel d'un produit de convolution :
  205. c XDEP=\int_0^T h(\tau)*Qj(t-\tau) d\tau
  206. if(jfonct.eq.100) then
  207. IP1=IPALA(jpliai,4)
  208. c t_{n+1} ou t_{n}
  209. if (IND.eq.1.or.IND.eq.3) then
  210. IP2=IPALA(jpliai,5)
  211. c t_{n+1/2} ou t_{n-1/2}
  212. else
  213. IP2=IPALA(jpliai,6)
  214. endif
  215. CALL DYCPL1(IP1,IP2,XDEP,NPAS,PDT,XCONV)
  216. XFL=XCPLGE*XCONV
  217. c WRITE(*,*) 't=',(NPAS*PDT),'x=',XDEP,'F1=',XFL
  218. elseif(jfonct.eq.101) then
  219. cbp calcul d'un produit de convolution GRANGER_PAIDOUSSIS
  220. IP1=IPALA(jpliai,4)
  221. IP2=IPALA(jpliai,5)
  222. IP3=IPALA(jpliai,6)
  223. VSD=XPALA(jpliai,2)
  224. XA0=XPALA(jpliai,3)
  225. CALL DYCPL2(IP1,IP2,IP3,VSD,XA0,XDEP,NPAS,PDTS2,XCONV)
  226. XFL=XCPLGE*XCONV
  227. c WRITE(*,*) '>DYCPL2:t=',T,' x=',XDEP,' F1=',XFL
  228. else
  229. c produit XFL=XCPLGE*(XDEP**EXPOSANT) ou XCPLGE*XDEP
  230. c CALL DYCPLD(XDEP,XCPLGE,XFL,iannul)
  231. cbp on mutliplie eventuellement par une fonction temporelle
  232. if(jfonct.ge.1.and.jfonct.le.2) then
  233. CALL DYCPLD(XDEP,XCPLGE,XFL,iannul)
  234. XFREQ = XPALA(jpliai,2)
  235. XTIME = DBLE(NPAS-1)*PDT
  236. if(IND.EQ.2) XTIME=XTIME-PDTS2
  237. if(jfonct.eq.1) XFONCT = COS(XFREQ*XTIME)
  238. if(jfonct.eq.2) XFONCT = SIN(XFREQ*XTIME)
  239. XFL = XFL * XFONCT
  240. else
  241. if(iannul.eq.0) then
  242. Xexpo = XPALA(jpliai,2)
  243. XFL = XCPLGE * (XDEP**Xexpo)
  244. else
  245. XFL=0.D0
  246. endif
  247. endif
  248. endif
  249. XVALA(jpliai,IND,1) = XFL
  250. XVALA(jpliai,IND,4) = XDEP
  251. FTest(INA1,IND) =Ftest(INA1,IND) + XFL
  252. *
  253. * ------ force elementaire de type POLYNOMIALE
  254. *
  255. ELSE IF (JTYP.EQ.6) THEN
  256. * nombre de modes "origine"
  257. NMOD = IPALA(I,2)
  258. CALL DYPOL1(Q1,Q2,NA1,IPLIA,XPALA,XVALA,NLIAA,IND,PDT,
  259. & NPAS,jpliai,NMOD,FTest,IVINIT,iannul)
  260. *
  261. * ------ choc elementaire ...
  262. *
  263. * ELSE IF (JTYP.EQ. ) THEN
  264. * .......
  265. * .......
  266. *
  267. ENDIF
  268.  
  269. * 2> TESTS POUR VOIR SI ON ANNULERA
  270. xff = 0d0
  271. do 104 ik = 1,na1
  272. do 105 jk = 1,4
  273. xff = xff + ( ftest(ik,jk) ** 2)
  274. 105 continue
  275. 104 continue
  276. xff = xff ** .5
  277. if ( ((xff .le. 1e-20 ) .and. ( jliai .gt. 0) )
  278. & .OR. ((xff .gt. 1e-20 ) .and. ( jliai .lt. 0) ) )
  279. & then
  280. iannul = 1
  281. endif
  282.  
  283. 101 CONTINUE
  284. * --> fin de boucle sur j
  285.  
  286. 199 CONTINUE
  287. * >>> FIN CAS DES LIAISONS CONDITIONNELLES <<<
  288. *
  289.  
  290. * ------ choc elementaire POINT_PLAN sans amortissement
  291. *
  292. IF (ITYP.EQ.1) THEN
  293. XRAID = XPALA(I,1)
  294. XJEU = XPALA(I,2)
  295. XNORM = 1.
  296. IF (XJEU.LT.0D0) THEN
  297. XNORM = -1.
  298. XJEU = -1*XJEU
  299. ENDIF
  300. INA1 = IPLIA(I,1)
  301. XDEP = XNORM*Q1(INA1,IND)
  302. CALL DYCHEL(XDEP,XRAID,XJEU, XFL,iannul)
  303. XVALA(I,IND,1) = XNORM*XFL
  304. XVALA(I,IND,4) = XNORM*XDEP
  305. FTOTA(INA1,IND) = FTOTA(INA1,IND) + XNORM*XFL
  306. *
  307. * ------ choc elementaire POINT_PLAN avec amortissement
  308. *
  309. ELSE IF (ITYP.EQ.2) THEN
  310. XRAID = XPALA(I,1)
  311. XJEU = XPALA(I,2)
  312. XAMO = XPALA(I,3)
  313. XNORM = 1.
  314. IPERM = 0
  315. IF (XJEU.LT.XZERO) THEN
  316. XNORM = -1.
  317. XJEU = -1*XJEU
  318. ENDIF
  319. INA1 = IPLIA(I,1)
  320. XDEP = XNORM * Q1(INA1,IND)
  321. IF ((NPAS.EQ.1).AND.(IND.EQ.3)) THEN
  322. XVIT = XNORM * Q2(INA1,IND)
  323. ELSE
  324. IND2 = IND + 1
  325. XDEPM1 = XNORM * Q1(INA1,IND2)
  326. XVIT = (XDEP - XDEPM1) / PDTS2
  327. ENDIF
  328. XVALA(I,IND,3) = XNORM*XVIT
  329. CALL DYCHAM(XDEP,XVIT,XRAID,XJEU,XAMO,XFL,IPERM,iannul)
  330. *+*
  331. IF (XDEP.GT.0.D0 .AND. XFL.GT.0.D0) XFL = 0.D0
  332. IF (XDEP.LT.0.D0 .AND. XFL.LT.0.D0) XFL = 0.D0
  333. *+*
  334. XVALA(I,IND,1) = XNORM*XFL
  335. XVALA(I,IND,4) = XNORM*XDEP
  336. FTOTA(INA1,IND) = FTOTA(INA1,IND) + XNORM*XFL
  337. *
  338. * ------ choc elementaire POINT_PLAN_FLUIDE
  339. *
  340. ELSE IF (ITYP.EQ.3) THEN
  341. XINER = XPALA(I,1)
  342. XCONV = XPALA(I,2)
  343. XVISC = XPALA(I,3)
  344. XPCEL = XPALA(I,4)
  345. XPCRA = XPALA(I,5)
  346. XJEU = XPALA(I,6)
  347. INA1 = IPLIA(I,1)
  348. XDEP = Q1(INA1,IND)
  349. IF ((NPAS.EQ.1).AND.(IND.EQ.3)) THEN
  350. XVIT = Q2(INA1,IND)
  351. ELSE
  352. IND2 = IND + 1
  353. XDEPM1 = Q1(INA1,IND2)
  354. XVIT = (XDEP - XDEPM1) / PDTS2
  355. ENDIF
  356. IF (XJEU.GT.0D0) THEN
  357. XDH= XJEU - XDEP
  358. XNORM = 1.
  359. ELSE
  360. XDH= XDEP - XJEU
  361. XNORM = -1.
  362. ENDIF
  363. * Calcul de la masse ajoutee
  364. XXIN = -XINER / XDH
  365. FINERT(INA1,IND) = FINERT(INA1,IND) + XXIN
  366. * Calcul de la force de convection
  367. CALL DYFCON(XDH,XDEP,XVIT,XJEU,XCONV,XFCO,iannul)
  368. * Calcul de la force de viscosite
  369. CALL DYFVIS(XDH,XDEP,XVIT,XJEU,XVISC,XFVI,iannul)
  370. * Calcul de la force de perte de charge
  371. CALL DYFPDC(XDH,XDEP,XVIT,XJEU,XPCEL,XPCRA,XFPE,
  372. & iannul)
  373. XFL = (XFCO* XNORM ) + XFVI + XFPE
  374. XVALA(I,IND,1) = XDEP
  375. XVALA(I,IND,2) = XVIT
  376. XVALA(I,IND,3) = XXIN
  377. XVALA(I,IND,4) = XFCO*XNORM
  378. XVALA(I,IND,5) = XFVI
  379. XVALA(I,IND,6) = XFPE
  380. FTOTA(INA1,IND) = FTOTA(INA1,IND) + XFL
  381. *
  382. * ------ force elementaire de COUPLAGE EN VITESSE
  383. *
  384. ELSE IF (ITYP.EQ.4) THEN
  385. INA1 = IPLIA(I,1)
  386. INA2 = IPLIA(I,2)
  387. XDEP = Q1(INA2,IND)
  388. IF ((NPAS.EQ.1).AND.(IND.EQ.3)) THEN
  389. XVIT = Q2(INA2,IND)
  390. ELSE
  391. IND2 = IND + 1
  392. XDEPM1 = Q1(INA2,IND2)
  393. XVIT = (XDEP - XDEPM1) / PDTS2
  394. ENDIF
  395. XCPLGE = XPALA(I,1)
  396. CALL DYCPLV(XVIT,XCPLGE,XFL,iannul)
  397. c WRITE(*,*) 't=',(NPAS*PDT),'dx/dt=',XVIT,'F2=',XFL
  398. XVALA(I,IND,3) = XVIT
  399. XVALA(I,IND,1) = XFL
  400. XVALA(I,IND,4) = XDEP
  401. FTOTA(INA1,IND) = FTOTA(INA1,IND) + XFL
  402. *
  403. * ------ force elementaire de COUPLAGE EN DEPLACEMENT
  404. *
  405. ELSE IF (ITYP.EQ.5) THEN
  406. INA1 = IPLIA(I,1)
  407. INA2 = IPLIA(I,2)
  408. XDEP = Q1(INA2,IND)
  409. XCPLGE = XPALA(I,1)
  410. jfonct = ipala(I,3)
  411. cbp calcul eventuel d'un produit de convolution :
  412. c XDEP=\int_0^T h(\tau)*Qj(t-\tau) d\tau
  413. if(jfonct.eq.100) then
  414. IP1=IPALA(I,4)
  415. c t_{n+1} ou t_{n}
  416. if (IND.eq.1.or.IND.eq.3) then
  417. IP2=IPALA(I,5)
  418. c t_{n+1/2} ou t_{n-1/2}
  419. else
  420. IP2=IPALA(I,6)
  421. endif
  422. CALL DYCPL1(IP1,IP2,XDEP,NPAS,PDT,XCONV)
  423. XFL=XCPLGE*XCONV
  424. c WRITE(*,*) 't=',(NPAS*PDT),'x=',XDEP,'F1=',XFL
  425. elseif(jfonct.eq.101) then
  426. cbp calcul d'un produit de convolution GRANGER_PAIDOUSSIS
  427. IP1=IPALA(I,4)
  428. IP2=IPALA(I,5)
  429. IP3=IPALA(I,6)
  430. VSD=XPALA(I,2)
  431. XA0=XPALA(I,3)
  432. CALL DYCPL2(IP1,IP2,IP3,VSD,XA0,XDEP,NPAS,PDTS2,XCONV)
  433. XFL=XCPLGE*XCONV
  434. c WRITE(*,*) '>DYCPL2:t=',T,' x=',XDEP,' F1=',XFL
  435. else
  436. c produit XFL=XCPLGE*XDEP
  437. c CALL DYCPLD(XDEP,XCPLGE,XFL,iannul)
  438. cbp on mutliplie eventuellement par une fonction temporelle
  439. if(jfonct.ge.1.and.jfonct.le.2) then
  440. CALL DYCPLD(XDEP,XCPLGE,XFL,iannul)
  441. XFREQ = XPALA(I,2)
  442. XTIME = DBLE(NPAS-1)*PDT
  443. if(IND.EQ.2) XTIME=XTIME-PDTS2
  444. if(jfonct.eq.1) XFONCT = COS(XFREQ*XTIME)
  445. if(jfonct.eq.2) XFONCT = SIN(XFREQ*XTIME)
  446. XFL = XFL * XFONCT
  447. else
  448. if(iannul.eq.0) then
  449. Xexpo = XPALA(I,2)
  450. if(Xexpo.eq.1.D0) then
  451. XFL = XCPLGE * XDEP
  452. else
  453. XFL = XCPLGE * (XDEP**Xexpo)
  454. endif
  455. else
  456. XFL=0.D0
  457. endif
  458. endif
  459. endif
  460. XVALA(I,IND,1) = XFL
  461. XVALA(I,IND,4) = XDEP
  462. FTOTA(INA1,IND) = FTOTA(INA1,IND) + XFL
  463. *
  464. * ------ force elementaire de type POLYNOMIALE
  465. *
  466. ELSE IF (ITYP.EQ.6) THEN
  467. * nombre de modes "origine"
  468. NMOD = IPALA(I,2)
  469. CALL DYPOL1(Q1,Q2,NA1,IPLIA,XPALA,XVALA,NLIAA,IND,PDT,
  470. & NPAS,I,NMOD,FTOTA,IVINIT,iannul)
  471. *
  472. * ------ choc elementaire ...
  473. *
  474. * ELSE IF (ITYP.EQ. ) THEN
  475. * .......
  476. * .......
  477. *
  478. ENDIF
  479.  
  480. *
  481. * la suite n'est plus utile car on passe iannul aux s_p de calcul des
  482. * forces de liaisons.
  483.  
  484. * si la liaison etait annulee on l'annule
  485. * if ( ( icond.eq. 1 ) .and. ( iannul.eq.1)) then
  486. * on annulle l'increment de ftotb
  487. * do 112 ik = 1,na1
  488. * do 113 jk = 1,4
  489. * ftota (ik,jk) = ftota0(ik,jk)
  490. * 113 continue
  491. * 112 continue
  492.  
  493. * end if
  494.  
  495. 10 CONTINUE
  496. *--------------------------------------------------------------------*
  497. * Fin de Boucle sur le nombre de liaisons
  498. *--------------------------------------------------------------------*
  499. *
  500. END
  501.  
  502.  
  503.  
  504.  
  505.  
  506.  
  507.  
  508.  
  509.  
  510.  
  511.  
  512.  
  513.  

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