OPTI DIME 3 ELEM pri6; *===================================================== * ETUDE THERMIQUE DU TUNNEL EST DE L'A86 *===================================================== * * Prarmètres de l'étude * R=5.20; ev=0.42; ed=0.25; L=1; *T=25; * * La densité spécifie la taille des éléments à générer * densite 0.8; * *----------------------------------------------------- * Définition des points *----------------------------------------------------- * O= 0. 0. 0.; * *Points des dalles Sup/Inf * p1= (R). 0. 0.; p5= (R)*(-1). 0. 0.; p6=(((R)*(-1))*((3**0.5)/2)). (((R)*(-1))/2). 0.; p10=((R)*((3**0.5)/2)). (((R)*(-1))/2). 0.; * *Points des voussoirs * p11= (R). 0. 0.; p102= (R). 0. 0.; p12= ((R)/(2**0.5)). (R/(2**0.5)). 0.; p21= ((R)/(2**0.5)). (R/(2**0.5)). 0.; p22= 0. (R). 0.; p31= 0. (R). 0.; p32= (((R)*(-1))/(2**0.5)). ((R)/(2**0.5)). 0.; p41= (((R)*(-1))/(2**0.5)). ((R)/(2**0.5)). 0.; p42= (R)*(-1). 0. 0.; p51= (R)*(-1). 0. 0.; p6v=(((R)*(-1))*((3**0.5)/2)). (((R)*(-1))/2). 0.; p62=(((R)*(-1))/(2**0.5)). (((R)*(-1))/(2**0.5)). 0.; p71=(((R)*(-1))/(2**0.5)). (((R)*(-1))/(2**0.5)). 0.; p72= 0. ((R)*(-1)). 0.; p81= 0. ((R)*(-1)). 0.; p82= ((R)/(2**0.5)). (((R)*(-1))/(2**0.5)). 0.; p91= ((R)/(2**0.5)). (((R)*(-1))/(2**0.5)). 0.; p10v=((R)*((3**0.5)/2)). (((R)*(-1))/2). 0.; * * Creation du terrain * pa=(2*R*(-1)). (2*R). 0.; pb=(2*R). (2*R). 0.; pc=(2*R). (2*R*(-1)). 0.; pd=(2*R*(-1)). (2*R*(-1)). 0.; * *----------------------------------------------------- * Construction de la figure *----------------------------------------------------- * *Dalles * ds= d p5 p1; di= d p6 p10; * *Voussoirs * a1= p11 cerc O p12 ; a2= p21 cerc O p22 ; a3= p31 cerc O p32 ; a4= p41 cerc O p42 ; a5= p51 cerc O p6v ; a6= p6v cerc O p62 ; a7= p71 cerc O p72; a8= p81 cerc O p82 ; a9= p91 cerc O p10v; a10= p10v cerc O p102 ; * * Terrain * da= d 6 pa pb; db= d 6 pb pc; dc= d 6 pc pd; dd= d 6 pd pa; * * On reunit les elements de meme nature * c1 = a1 et a2 et a3 et a4; c2 = a5 et a6 et a7 et a8 et a9 et a10; voutot2 = a1 et a2 et a3 et a4 et a5 et a6 et a7 et a8 et a9 et a10; dalletot = ds et di; elim voutot2; contter = da et db et dc et dd; terrtot= surf (contter et voutot2) plane; v0 = 0. 0. (-30).; vol = terrtot volu tran v0; sd = tran dalletot v0; sc1= tran c1 v0; sc2= tran c2 v0; tunvol = sd et sc1 et sc2; elim tunvol; tronc= vol et (tunvol coul vert); elim tronc; trac tronc; * *----------------------------------------------------- * Conditions Limites *----------------------------------------------------- * cl1 = bloq t c1; rig1 = symt tronc depl rota pa pb pc; rigtot = rig1; * *----------------------------------------------------- * Modèle des éléments *----------------------------------------------------- * motun = mode tunvol mecanique elastique isotrope coq4 cons un; motunth = mode tunvol thermique isotrope cons deux; movol = mode vol mecanique elastique isotrope cons un; movolth = mode vol thermique isotrope cons deux; motuntot = motun et motunth; movoltot = movol et movolth; moth = motuntot et movoltot; * *----------------------------------------------------- * Caractéristiques Matérielles géométriques *----------------------------------------------------- * matun = mate motun alpha 13.6E-6 youn 30.E9 nu 0.2 epai ev; matunth = mate motunth K 2.3 'C' 1000. RHO 2500.; mavol = mate movol alpha 7.6E-6 youn 8.E6 nu 0.45; mavolth = mate movolth K 1.5 'C' 2000. RHO 2000.; matuntot = matun et matunth; mavoltot = mavol et mavolth; math = matuntot et mavoltot; * *===================================================================* * Températures imposées * *===================================================================* * zer1=zero moth temperatures; zer2=chan chpo zer1 moth; zer3=zer2**0; zer4=20. * zer3; *tem0 = manu chpo c1 1 t 20 ; trac (exco zer4 T) tronc; trac (exco zer4 Tsup) tronc; temi = depi cl1 zer4; evo1 = evol manu 'temps' (prog 0 3600 10800) 't' (prog 20 1000 1000); char1 = char temi evo1; * *===================================================================* * Champs relatifs a la procedure PASAPAS * *===================================================================* * tab1=table; tab1.modele=moth; tab1.caracteristiques=math; tab1.blocages_mecaniques=rigtot; tab1.blocages_thermiques = cl1; tab1.temperatures = table; tab1.temperatures . 0 = zer4; tab1. talpha_reference= 20.; tab1.chargement = char1; tab1.temps_calcules= prog 0. pas 500. 3500. ; pasapas tab1;