Télécharger soudage6.dgibi

Retour à la liste

* fichier soudage6.dgibi
* section : thermique conduction convection
*----------------------------------------------------------------------*
*                      S O U D A G E 6 . D G I B I                     *
*----------------------------------------------------------------------*
* Objet :
* -------
*
*    Exemple de simulation thermique du soudage d'un raidisseur sur une
* plaque avec apport de matiere (4 passes). Les materiaux sont en acier
* inoxydable (316 L).
 
*    On modelise la conduction, la convection avec le milieu ambiant
* et la source de chaleur selon une distribution gaussienne.
*
*    L'apport de matiere est defini a l'aide d'un chargement de nom MODE
* qui decrit l'evolution du modele au cours du temps. Un chargement de
* nom MATE definit celle de ses caracteristiques.
*
*    La procedure SOUDAGE permet de definir la sequence de soudage et
* d'indexer les elements du maillage des cordons de soudure en fonction
* de l'evolution de l'outil le long de la trajectoire de soudage.
*
* Description :
* -------------
*
* Type de calcul : Thermique Transitoire
* Mode de calcul : 3D
* Type d'element : CUB8, PRY6
* Chargement     : Source de chaleur, Convection, Apport de matiere
*
*----------------------------------------------------------------------*
*
opti dime 3 elem cub8 mode trid ;
 
* Pour affichage traces, mettre IG1 a VRAI :
IG1      = faux ;
 
* Pour calcul complet, mettre icomplet a VRAI :
icomplet = faux ;
 
*opti trac psc ;
ig1      = ig1 ou (ega (vale trac) 'PSC') ;
 
*----------------------------------------------------------------------*
*                          Donnees du probleme                         *
*----------------------------------------------------------------------*
*
*----------------------- Parametres geometriques ----------------------*
 
* Geometrie pieces a souder :
* la1    : largeur plaque (m)
* hr1    : hauteur du raidisseur (m)
* lw1    : longueur a souder (m)
* ep1    : epaisseur plaque (m)
* er1    : epaisseur raidisseur (m)
la1      = 60.e-3 ;
hr1      = 40.e-3 ;
lw1      = 30.e-3 ;
ep1      = 30.e-3 ;
er1      = 10.e-3 ;
de1      =  1.5e-3 ;
de2      =  5.e-3 ;
de3      = 15.e-3 ;
 
* Geometrie cordons :
* lc1    : largeur cordons (m)
lc1      = 6.e-3 ;
* Pour la geometrie modelisee, section de cordon identique :
Sc1      = 0.5 * lc1 * lc1 ;
* Positions cordons :
Pc1      = (0.5*er1*(1 0 0)) plus ((lc1 lc1 0.) * 0.5) ;
Pc2      = Pc1 plus (lc1 0. 0.) ;
Pc3      = Pc1 plus ((lc1 lc1 0.) * 0.5) ;
Pc4      = Pc1 plus (0. lc1 0.) ;
 
*--------------------- Donnees probleme thermique ---------------------*
 
* CI / CL :
* Tini1  : temperature initiale :
* Te1    : temperature milieu ambiant :
* h1     : coefficient d'echange convectif (W/m2/K)
* em1    : emissivite 
* tps0   : temps de maintien avant debut deplacement torche (s)
Tini1    = 20. ;
Te1      = 20. ;
h1       = 20. ;
em1      = 0.5 ;
tps0     =  1. ;
 
* Parametres soudage :
* Us1    : tension de soudage (V)
* Is1    : Intensite de soudage (A)
* Vs1    : Vitesse de soudage (m/s)
* Vd1    : Vitesse de deplacement (m/s)
* Vf1    : vitesse de debit de fil (m/s)
* Df1    : diametre du fil (m) : pas utilisee dans le calcul
* Xpas1  : pas de discretisation de l'apport de matiere (m)
Us1      = 12. ;
Is1      = 400. ;
Vs1      = 0.3 / 60. ;
Vd1      = 0.1 * Vs1 ;
Df1      = 1.2e-3 ;
Xpas1    = 5.e-3 ;
 
* L'intensite de soudage est deduite du debit volumique de matiere et
* de l'energie thermique qu'il faut fournir pour la fondre (on ne tient
* pas compte du rendement).
* A titre indicatif, calcul effectue pour l'estimer (avec : rho = 7.6e3 kg/m3,
* Cp = 500. J/kg/m3, dT = 1500. degC, Qlat = 1.e9 J/m3)
Is0      = Sc1 * Vs1 * (7.6e3 * 500. * 1500. + 1.e9) / Us1 ;
list Is0 ;
 
* La vitesse de fil est deduite de la section des cordons et de la vitesse
* de soudage. Elle n'intervient pas dans le calcul car geometrie cordon fixee.
* A titre indicatif, son calcul et affichage (m/min.) :
Vf1      = Sc1 * Vs1 / (pi * Df1 * Df1 /4.) ;
List (Vf1 * 60.) ;
 
* Parametres modele source Gaussienne 3D :
* Qtot1  : puissance thermique de la source (J)
* Rg1    : "rayon" de la source Gaussienne (m)
* Zg1    : "profondeur" de la source Gaussienne (m)
Qtot1    = Us1 * Is1 ;
Rg1      = 6.0e-3 ;
Zg1      = 4.5e-3 ;
 
*-------------------------- Donnees materiau --------------------------*
* Ref. :
* Camille Cambon, Issam Bendaoud, Sebastien Rouquette, Fabien Soulie.
* “Influence of the first weld bead on strain and stress states in wire+arc additive manufacturing”.
* The 12th International Seminar ”Numerical Analysis of Weldability”,
* Institute for Materials Science, Joining and Forming (IMAT),Sep 2018, Seggau, Austria.
* hal-01954354 https://hal.archives-ouvertes.fr/hal-01954354
*
* rho1   : masse volumique (kg/m3)
* Tfus1  : temperature de changement de phase solide-liquide (degC)
* Qlat1  : chaleur latente de fusion de l'acier (J/m3)
rho1     = 7760. ;
Tfus1    = 1450. ;
Qlat1    =  2.e9 ;
 
* k1     : conductivite thermique (W/m/K) en fonction de la temperature
lt1      = prog 20. 100. pas 100. 1500. ;
lk1      = prog 14. 15.2 16.6 17.9 19. 20.6 21.8 23.1 24.3 26. 27.3 28.6 29.9 45. 60. 60. ;
k1       = evol vert manu 'T' lt1 'K' lk1 ;
si ig1 ;
  dess k1 titr 'Conductivite thermique 316L (W/m/K)' ;
fins ;
 
* Cp1    : capacite thermique massique (J/m3) en fonction de la temperature
* On derive l'enthalpie :
lt1      = prog 0. pas 100. 1500. ;
lh1      = prog 0.983 1.36 1.76 2.19 2.62 3.06 3.51 3.96 4.43 4.92 5.42 5.92 6.44 6.97 9.51 ;
lh2      = prog 1.36 1.76 2.19 2.62 3.06 3.51 3.96 4.43 4.92 5.42 5.92 6.44 6.97 9.51 10.1 ;
dhdt1    = (lh2 - lh1) / 100. ;
lc0      = dhdt1 * 1.e9 / rho1 ;
lt0      = prog 50. pas 100. 1450. ;
cp0      = evol manu 'T' lt0 'C' lc0 ;
idx0     = posi (maxi lc0) dans lc0 ;
lc2      = lc0 enle idx0 ;
lt2      = lt0 enle idx0 ;
cp2      = evol oran manu 'T' lt2 'C' lc2 ;
cp1      = cp2 ;
si ig1 ;
  tleg1    = table ;
  tleg1 . titre = table ;
  tleg1 . 1     = 'Cp derivee publi.' ;
  tleg1 . 2     = 'Cp utilisee' ;
  dess (cp0 et cp1) titr ' Capacite thermique massique 316L (J/m3)' lege tleg1 ;;
fins ;
 
*----------------------------------------------------------------------*
*                               Maillage                               *
*----------------------------------------------------------------------*
 
* Repere :
O1       = 0 0 0 ;
X1       = 1 0 0 ;
Y1       = 0 1 0 ;
Z1       = 0 0 -1 ;
oeil0    = 100 +100 100 ;
 
*-------------------- Maillage Plaque et Raidisseur -------------------*
 
* Plaque : maillage surface dans plan (0,xy)
Pp1      = 0.5*er1*X1 ;
Pp2      = Pp1 plus (3.*lc1*X1) ;
lp1      = ((0 0 0) droi Pp1 dini de1 dfin de1) et (Pp1 droi Pp2 dini de1 dfin de1) ;
sp1      = lp1 tran (-0.4*ep1*Y1) dini de1 dfin de1 ;
sp1      = sp1 coul gris ;
lp2      = sp1 cote 2 ;
lp3      = sp1 cote 3 ;
Pp3      = lp2 poin proc (Pp2 moin Y1) ;
PP4      = Pp3 moin (0.6*ep1*Y1) ;
Pp5      = PP4 moin (0.5*er1*X1) moin (3.*lc1*X1) ;
Pp6      = lp3 poin proc (-1 -1 0) ;
CntSp2   = Pp3 droi Pp4 dini de1 dfin de2 ;
CntSp2   = CntSp2 et (Pp4 droi Pp5 dini de2 dfin de2) ;
CntSp2   = CntSp2 et (Pp5 droi Pp6 dini de2 dfin de1) ;
CntSp2   = CntSp2 et (inve lp3) ;
Sp2      = surf CntSp2 plan ;
Sp2      = sp2 coul gris ;
lp22     = CntSp2 elem comp Pp3 Pp4 ;
Pp7      = la1*X1 ;
Pp8      = Pp7 moin (ep1*Y1) ;
Pp9      = Lp22 poin proc (Pp2 moin (ep1*Y1)) ;
CntSp3   = Pp2 droi Pp7 dini de1 dfin de3 ;
CntSp3   = CntSp3 et (Pp7 droi Pp8 dini de3 dfin de3) ;
CntSp3   = CntSp3 et (Pp8 droi Pp9 dini de3 dfin de2) ;
CntSp3   = CntSp3 et (inve (Lp2 et lp22)) ;
Sp3      = surf CntSp3 plan ;
Sp3      = Sp3 coul gris ;
spz0     = sp1 et Sp2 et Sp3 ;
 
* Raidisseur : maillage surface dans plan (0,xy)
* Pr1 = Pp1 : Noeuds l'arete a l'angle de la plaque et du raidisseur commun aux 2 pieces :
Pr2      = Pp1 plus (2.*lc1*Y1) ;
lr1      = (0 0 0) droi Pp1 dini de1 dfin de1 ;
sr1      = lr1 tran (2.*lc1*Y1) dini de1 dfin de1 ;
Pr2      = sr1 poin proc Pr2 ;
sr1      = sr1 coul vert ;
lr3      = sr1 cote 3 ;
sr2      = lr3 tran ((hr1-lc1-lc1)*Y1) dini de1 dfin de2 ;
sr2      = sr2 coul vert ;
srz0     = sr1 et sr2 ;
 
* Extrusion maillage volumique :
Sz0      = Spz0 et Srz0 ;
Vp0      = spz0 volu tran (lw1*Z1) dini de1 dfin de1 ;
Vr0      = srz0 volu tran (lw1*Z1) dini de1 dfin de1 ;
Vs0      = Vp0 et Vr0 ;
 
*------------------------ Maillage des cordons ------------------------*
 
* Maillage surface cordons dans plan (O,x,y) :
opti elem tet4 ;
 
* Cordon 1 :
Pc12     = lp1 poin proc (Pp1 plus (lc1*X1)) ;
Pc13     = (sr1 cote 2) poin proc (Pp1 plus (lc1*Y1)) ;
CntC1    = lp1 elem comp Pp1 Pc12 ;
CntC1    = CntC1 et (Pc12 droi Pc13 dini de1 dfin de1) ;
CntC1    = CntC1 et ((sr1 cote 2) elem comp Pc13 Pp1) ;
SuC1     = surf CntC1 plan ;
Suc1     = Suc1 coul 1 ;
 
* Cordon 2 :
Pc22     = lp1  poin proc (Pp1 plus (2.*lc1*X1)) ;
Pc23     = Pc22 plus ((-0.5*lc1) (0.5*lc1) 0.) ;
Pc24     = Pc23 moin (lc1*X1) ;
Pc24     = Cntc1  poin proc Pc24 ;
CntC2    = lp1   elem comp Pc12 Pc22 ;
CntC2    = CntC2 et (Pc22  droi Pc23 dini de1 dfin de1) ;
CntC2    = CntC2 et (Pc23 droi Pc24 dini de1 dfin de1) ;
CntC2    = CntC2 et (((CntC1 inve) elem comp Pc24 Pc12) ) ;
Suc2     = Cntc2 surf plan ;
Suc2     = Suc2 coul 2 ;
 
* Cordon 3 :
Pc33     = Pc23 plus ((-0.5*lc1) (0.5*lc1) 0.) ;
CntC3    = (inve CntC2) elem comp Pc24 Pc23 ;
CntC3    = CntC3 et (Pc23 droi Pc33 dini de1 dfin de1) ;
CntC3    = CntC3 et (Pc33 droi Pc13 dini de1 dfin de1) ;
CntC3    = CntC3 et ((inve CntC1) elem comp Pc13 Pc24) ;
SuC3     = surf CntC3 plan ;
SuC3     = SuC3 coul 3 ;
 
* Cordon 4 :
CntC4    = (inve CntC3) elem comp Pc13 Pc33 ;
CntC4    = CntC4 et (Pc33 droi Pr2  dini de1 dfin de1) ;
CntC4    = CntC4 et ((sr1 cote 2) elem comp Pr2  Pc13) ;
SuC4     = CntC4 surf plan ;
SuC4     = SuC4 coul 4 ;
 
* Volume cordons :
Suc0     = Suc1 et Suc2 et Suc3 et Suc4 ;
*trac (0 0 1000) Suc0 ;
Vc0      = Suc0 volu tran (lw1*Z1) dini de1 dfin de1 ;
 
* Maillage volumique de chaque cordon :
Vc1      = Vc0 elem coul 1 ;
Vc2      = Vc0 elem coul 2 ;
Vc3      = Vc0 elem coul 3 ;
Vc4      = Vc0 elem coul 4 ;
 
* Fusion noeuds communs cordons/plaque et cordons/raidisseur mais pas plaque/raidisseur
elim ((enve Vc0) et (enve Vp0)) 1.e-5 ;
elim ((enve Vc0) et (enve Vr0)) 1.e-5 ;
 
* Maillage total :
V0       = Vc0 et Vs0 ;
 
* Affichages :
si IG1 ;
  mot1     = chai ' Maillage plaque et raidisseur :' (nbno Vs0) ' noeuds,' (nbel Vs0) ' elem.' ;
  trac oeil0 face Vs0 titr mot1 ;
  mbox2    = boite (Sp1 et Sr1 et Vc0) ;
  mot1     = chai ' Maillage plaque, raidisseur et cordons :' (nbno V0) ' noeuds,' (nbel V0) ' elem.' ;
  trac oeil0 face V0 titr mot1 ;
  trac (0 1 10) face V0 boit mbox2 titr mot1 ;
fins ;
 
*----------------------------------------------------------------------*
*                          Sequence de soudage                         *
*----------------------------------------------------------------------*
 
* Initialisation table de soudage :
tsoud1                          = table ;
tsoud1 . vitesse_de_soudage     = Vs1 ;
tsoud1 . vitesse_de_deplacement = Vd1 ;
tsoud1 . puissance_de_soudage   = Qtot1 ;
tsoud1 . diametre_de_fil        = Df1 ;
tsoud1 . vitesse_de_fil         = Vf1 ;
tsoud1 . point_de_depart        = Pc1 ;
 
* Specification de la sequence :
* Delai de 0,5 s avant de faire la 1ere passe :
soudage tsoud1 point tps0 ;
* 1ere passe :
soudage tsoud1 passe droi (lw1*Z1) ;
* 2e passe :
soudage tsoud1 depla droi Pc2 abso ;
soudage tsoud1 point 600. puis 0. debi 0. ;
soudage tsoud1 passe droi (lw1*Z1) ;
* 3e passe :
soudage tsoud1 depla droi Pc3 abso ;
soudage tsoud1 point 600. puis 0. debi 0. ;
soudage tsoud1 passe droi (lw1*Z1) ;
* 4e passe :
soudage tsoud1 depla droi Pc4 abso ;
soudage tsoud1 point 600. puis 0. debi 0. ;
soudage tsoud1 passe droi (lw1*Z1) ;
* Refroidissement 
soudage tsoud1 point 600. puis 0. debi 0. ;
 
* Sequencage maillage des cordons le long de la trajectoire :
tab2     = soudage tsoud1 'MAIL' vc0 'PAS' Xpas1 'TEMP' ;
 
ttps1    = tab2.evolution_maillage.temps ;
 
si ig1 ;
  titr ' Trajectoire soudage des 4 passes (vert : depl. outil)' ;
  trac (1 2 10) (tsoud1.trajectoire et (aret vc1) et (aret vc2 coul defa) et (aret vc3 coul defa) et (aret vc4 coul defa)) ;
  dess tsoud1.evolution_deplacement titr 'Evolution temporelle du deplacement sur la trajectoire' ;
  dess tsoud1.evolution_puissance titr 'Evolution tempeorelle de la puissance thermique' ;
*
  opti oeil (1. 2. 3.) ;
  waam tab2 visu cach ;
*dess (evol manu tab2.temps_calcules (prog 1. pas 1. (dime tab2.temps_calcules))) titr 'Temps calcules' ;
fins ;
 
*----------------------------------------------------------------------*
*                            Modele Physique                           *
*----------------------------------------------------------------------*
 
* Accrochage DDL thermique bord raidisseur sur plaque :
clt1     = rela accro (enve Vr0) vp0 (mots T) ;
 
* sinte1 : surface interieure : a retirer de la surface qui convecte 
* partie du raidisseur en vis-a-vis de la plaque (surface y = ymin.)
ptry0    = (vr0 coor 2) poin mini ;
sury0    = (vr0 enve) elem appu stri ptry0 ;
sinte1   = sury0 ;
* partie de la plaque en vis-a-vis du raidisseur (maillage non-conforme)
pclt1    = clt1 extr mail nomu ;
supi1    = (enve Vp0) elem appu stri (pclt1 et sury0) ;
sinte1   = sinte1 et supi1 ;
* plan de symetrie en x=0
ptx0     = (V0 coor 1) poin infe 1.e-6 ;
sux0     = (enve V0) elem appu stri ptx0 ;
sinte1   = sinte1 et sux0 ;
 
si ig1 ;
  trac (-5 1 3) cach ((enve V0) diff sinte1) boit (boite vc0) titr 'Surface de convection' ;
fins ;
 
*---------------------- Modele / Carateristiques ----------------------*
 
* Conduction / Convection :
mod1     = mode V0 thermique ;
mat1     = mate mod1 rho rho1 k k1 'C' cp1 'TINI' Tini1 ;
 
* Source thermique :
evqt1    = tsoud1.evolution_puissance ;
evxs1    = tsoud1.evolution_deplacement ;
chxs1    = tsoud1.trajectoire coor curv ;
cgxs1    = char traj chxs1 evxs1 ;
mod3     = mode V0 thermique source gaussienne elliptique ;
mat3     = mate mod3 qtot evqt1 orig cgxs1 rgau Rg1 dire (1. 1. 0.) zgau Zg1 ;
 
* Tables d'evolution des modeles et caracteristiques :
ttps1    = tab2.evolution_maillage.temps ;
tmai1    = tab2.evolution_maillage.maillage ;
nb1      = dime ttps1 ;
tmod1    = table ;
tmod2    = table ;
tmod3    = table ;
tmod4    = table ;
tmat1    = table ;
tmat2    = table ;
tmat3    = table ;
tmat4    = table ;
repe b1 nb1 ;
  geoi1    = tmai1 . (&b1 - 1) ;
  geoi1    = geoi1 et Vp0 et Vr0 ;
  tmod1 . (&b1 - 1) = redu mod1 geoi1 ;
  tmat1 . (&b1 - 1) = redu mat1 (tmod1 . (&b1 - 1)) ;
  sconv1   = (enve geoi1) diff sinte1 ;
*trac nclk cach sconv1 boit (boite vc0) ;
  tmod2 . (&b1 - 1) = mode sconv1 thermique convection ;
  tmat2 . (&b1 - 1) = mate (tmod2 . (&b1 - 1)) 'H' h1 'TC' Te1 ;
  tmod3 . (&b1 - 1) = redu mod3 geoi1 ;
  tmat3 . (&b1 - 1) = redu mat3 (tmod3 . (&b1 - 1)) ;
  tmod4 . (&b1 - 1) = mode sconv1 thermique rayonnement infini ;
  tmat4 . (&b1 - 1) = mate (tmod4 . (&b1 - 1)) 'EMIS' em1 'T_IN' Te1 ;
fin b1 ;
 
* Chargements MODE / MATE :
cgmod1   = char mode ttps1 tmod1 ;
cgmod2   = char mode ttps1 tmod2 ;
cgmod3   = char mode ttps1 tmod3 ;
cgmod4   = char mode ttps1 tmod4 ;
cgmod0   = cgmod1 et cgmod2 et cgmod3 et cgmod4 ;
 
cgmat1   = char mate ttps1 tmat1 ;
cgmat2   = char mate ttps1 tmat2 ;
cgmat3   = char mate ttps1 tmat3 ;
cgmat4   = char mate ttps1 tmat4 ;
cgmat0   = cgmat1 et cgmat2 et cgmat3 et cgmat4 ;
 
*------------------------- Resolution PASAPAS -------------------------*
 
* temps calcules :
ltca0    = tab2.temps_calcules ;
it0      = posi tps0 dans ltca0 ;
ltca0    = (prog 0. pas 0.2 tps0) et (ltca0 extr (lect (it0 + 1) pas 1 (dime ltca0))) ;
si (non icomplet) ;
  ltca1    = ltca0 extr (lect 1 pas 1 20) ;
sino ;
  ltca1    = ltca0 ;
fins ;
list ltca1 ;
 
* temps des sauvegardes :
si icomplet ;
  ltsg1    = ltca1 extr (lect 50 pas 50 ((dime ltca1) / 50 * 50)) ;
  si ((maxi ltsg1) neg (maxi ltca1)) ;
    ltsg1    = ltsg1 et (maxi ltca1) ;
  fins ;
  list ltsg1 ;
fins ;
 
tpas1                     = table ;
tpas1.modele              = cgmod0 ;
tpas1.caracteristiques    = cgmat0 ;
*tpas1.chargement          = cgmat0 ;
tpas1.blocages_thermiques = clt1 ;
tpas1.temps_calcules      = ltca1 ;
tpas1.temps_sauves        = ltca1 ;
tpas1.celsius             = vrai ;
 
si icomplet ;
  tpas1.temps_sauvegardes   = ltsg1 ;
  opti sauv 'thermique_SoudageRaidisseur.sauv' ;
  sauv ;
fins ;
 
temp zero ;
 
pasapas tpas1 ;
 
* Affichage temps de calcul :
duree1   = temp horl ;
duree1   = (duree1 / 1000) ;
dmin1    = (duree1  / 60) ;
dsec1    = duree1  - (60 * dmin1) ;
mot1     = chai '*****  DUREE DU CALCUL (s) :' dmin1 ' min' dsec1 ' s' ;
mess mot1 ;
 
*------------------------ Petit post-traitement -----------------------*
 
si ig1 ;
 
* Animation du champ de temperature :
  if1      = (dime tpas1.temps) - 1 ;
  modf1    = tpas1.modeles.if1 ;
  mailf1   = (extr modf1 mate conduction) extr mail ;
  dep0     = manu chpo mailf1 3 ux 0. uy 0. uz 0. ;
  def1     = vide deforme ;
  liso1    = prog 0. 50. 100. PAS 150. 1450. ;
  vtra1    = vale trac ;
  si (ega vtra1 'PSC') ;
    opti ftra 'Temperature_soudage6_anim.ps' ;
  fins ;
  mbox1    = boite Vc0 ;
  repe b1 if1 ;
    modi1    = tpas1.modeles.(&b1-1) ;
    maili1   = (extr modi1 mate conduction) extr mail ;
    chti1    = tpas1.temperatures.(&b1-1) ;
    defoi1   = defo maili1 dep0 chti1 ;
    def1     = def1 et defoi1 ;
    si (ega vtra1 'PSC') ;
      mot1     = chai format '(F6.1)' 'Temperature (degC) au temps (s) :' (tpas1.temps.(&b1-1)) ;
      trac (1 2 3) chti1 maili1 liso1 titr mot1 boit mbox1 ;
    fins ;
  fin b1 ;
  si (neg vtra1 'PSC') ;
    trac (1 2 3) anim def1 liso1 boit mbox1 ;
  fins ;
 
fins ;
 
si icomplet ; sauv ; fins ;
 
*--------------   F I N    S O U D A G E 6 . D G I B I   --------------*
fin ;
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

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