* fichier : thme2.dgibi
*----------------------------------------------------------------
*
* Réduction d'un jeu sous l'action d'une sollicitation
* thermique en régime transitoire
*
*----------------------------------------------------------------
*
* Calcul thermo-mécanique ( mécanique et thermique linéaire )
* Utilisation de la procédure PASAPAS
* Utilisation de la procédure utilisateur perso1.
*
*
* | x1 = 0.1 m x2 = 0.1002 m |
* | T1 = 500°C T2 = ? |
* |_______________________________ ____________________________|
* | | | |
* | T_initial = 20°C | | T_initial = 20°C |
* |_______________________________| |____________________________|
* | x0 = 0 m jeu x3 = 0.2002 m |
* | T0 = 20°C T3 = 20°C |
* | Encastrement Encastrement |
*
*
* Le jeu sera considéré comblé quand la distance entre les
* barreaux est inférieure à 0.1mm
*
* on étudie la diffusion thermique dans les barreaux et leurs
* dilatations.
*
* caracteristiques des barreaux :
*
* l = 0.1m
* K = 14.6
* C = 460
* RHO = 7800
* ALPHA = 1E-5
* YOUNG = 2.E11
* NU = 0.00001
*
* Mailles de 0.01 m de cote
*
*----------------------------------------------------------------
*
OPTION ECHO 0;
OPTION DIME 2 ELEM QUA8 MODE PLAN DEFO;
*
*-------------------- Creation du maillage ----------------------
*
DENS 1;
P1 = 0. 0. ;
P2 = 0.1 0. ;
P3 = 0.1 0.01 ;
P4 = 0. 0.01 ;
P5 = 0.1002 0.;
P6 = 0.2002 0.;
P7 = 0.2002 0.01;
P8 = 0.1002 0.01;
P9 = 0.12 0.;
D1 = DROI 10 P1 P2 ;
D2 = DROI 1 P2 P3 ;
D3 = DROI 10 P3 P4 ;
D4 = DROI 1 P4 P1 ;
D5 = DROI 10 P5 P6 ;
D6 = DROI 1 P6 P7 ;
D7 = DROI 10 P7 P8 ;
D8 = DROI 1 P8 P5 ;
L1 = D1 ET D2 ET D3 ET D4;
L2 = D5 ET D6 ET D7 ET D8;
SU1 = SURF L1 PLANE;
SU2 = SURF L2 PLANE;
P9 = SU2 POINT PROCHE P9;
*
*-------------------- Définition du modèle ----------------------
*
MOD1 = MODL (SU1 ET SU2) THERMIQUE ISOTROPE CONS UN;
MOD3 = MODL (SU1 ET SU2) MECANIQUE ELASTIQUE ISOTROPE CONS DEUX ;
MOD = MOD1 et MOD3;
*
*------------------- Définition des matériaux -------------------
*
MAT1 = MATR MOD1 K 14.6 'C' 460. RHO 7800. ;
MAT3 = MATR MOD3 ALPHA 1E-5 YOUN 2E+11 NU 0.00001;
MAT = MAT1 ET MAT3;
*
*------------------- blocages mécaniques ------------------------
*
CL1 = BLOQ UY (SU1 ET SU2);
CL2 = BLOQ UX D4;
CL3 = BLOQ UX D6;
CL_ME = CL1 ET CL2 ET CL3;
*
*------------------- blocages thermiques ------------------------
*
CL4 = BLOQ T D4;
CL5 = BLOQ T D2;
CL6 = BLOQ T D6;
TEM1 = DEPI (CL4 ET CL6) 20.;
TEM2 = DEPI CL5 500.;
CL_TH = CL4 ET CL5 ET CL6;
*
*------------------- températures imposées ----------------------
*
TEM1 = DEPI (CL4 ET CL6) 20.;
TEM2 = DEPI CL5 500.;
TIM = TEM1 ET TEM2;
EV1 = EVOL MANU TEMPS (PROG 0. 200.) y (prog 1. 1.);
CHA = CHAR 'TIMP' TIM EV1;
*
*------------------- température initiale -----------------------
*
TEM0 = MANU CHPO (SU1 ET SU2) 1 T 20.;
*
*--------------- creation de la table pour pasapas ------------
*
ta = table;
ta.modele = mod;
ta.caracteristiques = mat;
ta.blocages_mecaniques = cl_me;
ta.blocages_thermiques = cl_th;
ta.temperatures = table;
ta.temperatures . 0 = tem0;
ta.chargement = cha;
ta.procedure_thermique = DUPONT;
ta.procedure_perso1 = vrai;
ta.ok = faux;
ta.temps_calcules = (prog 0. pas 5. 80.) et (prog 81. pas 1. 89.)
et (prog 90. pas 5. 150.);
*
*------------ procédure utilisateur perso1 --------------------
*
debproc perso1 TTU*TABLE;
dep = ttu.'ESTIMATION'.'DEPLACEMENTS';
WTAB=TTU.'WTABLE';
x = extr dep UX P2;
list x;
si ((WTAB.'TEMPS0' > 84.5) et ( WTAB.'TEMPS0' < 89.5));
list ttu.'ESTIMATION'.'TEMPERATURES';
finsi;
si (('>' x 0.0001) et (non ttu.ok));
ttu.loup = WTAB.'TEMPS0';
mess ' *************************';
mess ' mise a jour du au contact pour le temps : ' ttu.loup ;
mess ' *************************';
CL7 = bloq T D8;
tim = depi cl7 500.;
WTAB.'BLOCAGES_THERMIQUES' = (WTAB.'BLOCAGES_THERMIQUES' ) et cl7;
EV1 = EVOL MANU TEMPS (PROG 0. 200.) y (prog 1. 1.);
CHA7 = CHAR 'TIMP' TIM EV1;
WTAB.'CHARGEMENT_I' = WTAB.'CHARGEMENT_I' et cha7;
WTAB.'CHARGEMENT' = WTAB.'CHARGEMENT' et cha7;
ttu.ok = vrai;
finsi;
finproc;
*
*------------ Appel à la procédure PASAPAS --------------------
*
pasapas ta;
*
*----------- Depouillement ------------------------------------
*
a = (14.6 / 7800.)/460.;
tpr1 = prog;
tpr3 = prog;
tep_th = prog;
x = 0.12;
er1 = 0.;
nb = 0;
idim = (dime (ta.temps)) - 1;
*
repeter bou idim;
ite = &bou;
te = ta.temps.ite;
tem = ta.temperatures.ite;
tep1 = extr tem T p9;
tpr1 = tpr1 et (prog tep1);
tpr3 = tpr3 et (prog te);
si ('>' te (ta.loup)) ;
co1 = erf((x - .1002)/(2 * (( a * (te - (ta.loup))**0.5))));
co2 = erf(0.1/(2.*(( a * (te - (ta.loup)))**0.5)));
sol_te = 500. - ((500. - 20.)*co1/co2);
er1 = er1 + (abs ((sol_te - tep1)/sol_te) * 100.);
nb = nb + 1;
sinon;
sol_te = 20.;
finsi;
tep_th = tep_th et (prog sol_te);
fin bou;
*
er1 = er1/nb;
mess 'erreur sur la température :'er1'%';
*
*lvar = TEXT 'EVOLUTION DE LA TEMPERATURE AU POINT 0.12 m';
*titre lvar;
*ev_te = evol manu temps tpr3 temperature tep_th;
*ev_tet = evol manu temps tpr3 temperature tpr1;
*dess (ev_tet et ev_te) ;
*
si ( er1 > 5. ) ;
mess 'Résultat incorrect' ;
erreur 5 ;
finsi ;
*
fin;
© Cast3M 2003 - All rights reserved.
Disclaimer