OPTION DIME 2 ELEM QUA4 AXIS ; **************************MAILLAGE************************* ******Piston Bas******** p1 = 0. 0. ; p2 = 19.06e-03 0. ; p3 = 61.56e-03 0. ; p4 = 100.025e-03 0. ; p5 = 100.025e-03 48.0e-03 ; p6 = 62.0e-03 48.0e-03 ; p7 = 20.05e-03 48.0e-03 ; p8 = 20.05e-03 51.0e-03 ; p9 = 0. 51.0e-03 ; p10 = 0. 48.0e-03 ; p11 = 0. 68.0e-03 ; *p12 = 5.3e-03 68.0e-03 ; *p13 = 6.2e-03 68.0e-03 ; *p14 = 7.5e-03 63.33577e-03 ; *p15 = 7.5e-03 66.4044e-03 ; p16 = 0. 68.0e-03 ; p17 = 5.3e-03 68.0e-03 ; p18 = 6.2e-03 68.0e-03 ; p19 = 6.2e-03 95.0e-03; p20 = 5.3e-03 95.0e-03; p21 = 0. 95.0e-03; p22 = 0. 89.0e-03 ; p23 = 3.0e-03 89.0e-03 ; p24 = 4.0e-03 89.0e-03 ; p25 = 4.0e-03 74.0e-03 ; p26 = 3.0e-03 74.0e-03 ; p27 = 2.5e-03 74.0e-03 ; p28 = 0. 74.0e-03 ; p29 = 0. 77.5e-03 ; p30 = 2.5e-03 77.5e-03 ; p31 = 2.5e-03 79.e-03 ; p32 = 2.2e-03 79.e-03 ; p33 = 0. 79.e-03 ; p34 = 2.2e-03 84.e-03 ; p35 = 2.5e-03 84.0e-03 ; p36 = 0. 84.0e-03 ; p37 = 0. 85.5e-03 ; p38 = 2.5e-03 85.5e-03 ; p39 = 2.5e-03 89.0e-03 ; p40 = 6.2e-03 89.0e-03 ; p41 = 6.2e-03 74.0e-03 ; p86 = 0. 63.33577e-3 ; p84 = 10.0e-3 63.33577e-3; l1 = p1 d 48 P2 ; l2 = p2 d 50 p3 ; l3 = p3 d 20 P4 ; l4 = p4 d 20 p5 ; l5 = p5 d 20 P6 ; l6 = p6 d 50 p7 ; l7 = p7 d 48 P10 ; l8 = p9 d 48 p8 ; l9 = p2 d 20 P7 ; l10 = p1 d 20 p10 ; l11 = p7 d 20 p8 ; l12 = p10 d 20 p9 ; l13 = p11 d 36 p86 ; l14 = p84 d 36 p18 ; l15 = p18 d 48 p11; l16 = p3 d 20 p6 ; l120 = p84 d 100 p8; l126 = p86 d 100 p9; *l127 = p2 d 100 p8; l128 = p84 d 48 p86; S1 = dall l1 l9 l7 l10 plan; S2 = dall l11 l8 l12 l7 plan; S3 = dall l14 l15 l13 l128 plan; S4 = dall l2 l16 l6 l9 plan; S5 = dall l3 l4 l5 l16 plan; S6 = dall l126 l8 l120 l128 plan; Pbas = S1 et S2 et S3 et S4 et S5 et S6; elim 1e-5 Pbas ; *trac Pbas; **********Echantillon et Environnement***** na1 = 150 ; na2 = 50 ; na3 = 35 ; na4 = 165 ; na5 = 29 ; l20 = p11 d 41 P17 ; l21 = p18 d 7 p17 ; l22 = p18 d 24 P41 ; l23 = p19 d 7 p20 ; l24 = p20 d 41 P21 ; l25 = p21 d 29 p22 ; l26 = p24 d 30 P20 ; l27 = p17 d 30 p25 ; l28 = p28 d 29 P11 ; l29 = p27 d 25 p28 ; l30 = p39 d 25 p22 ; l31 = p26 d na1 p23 ; l32 = p31 d na2 p35 ; l33 = p30 d 25 p29 ; l34 = p29 d na3 p28; l35 = p27 d na3 p30 ; l36 = p26 d 5 p27 ; l37 = p23 d 5 p39 ; l38 = p25 d 10 p26 ; l39 = p23 d 10 P24 ; l40 = p22 d na3 p37 ; l41 = p37 d 25 P38 ; l42 = p38 d na3 p39 ; l43 = p38 d 15 P35 ; l44 = p37 d 15 p36 ; l45 = p36 d 22 P34 ; l46 = p36 d na2 p33 ; l47 = p33 d 22 P32 ; l48 = p32 d na2 p34 ; l49 = p34 d 3 p35 ; l50 = p32 d 3 p31 ; l51 = p33 d 15 p29 ; l52 = p30 d 15 p31 ; l53 = p27 d 6 p30 ; l54 = p41 d 60 p40 ; l55 = p40 d 24 p19; l56 = p31 d 16 p20 ; l57 = p25 d 40 p28 ; l58 = p18 d na1 p19 ; l59 = p20 d na1 p17 ; l60 = p24 d na1 p25; l61 = p22 d 40 p24 ; l62 = p39 d na1 p27 ; l63 = p31 d 25 p33 ; l64 = p35 d 25 p36 ; inoxbas = dall l20 l27 l57 l28 plan; inoxhaut = dall l24 l25 l61 l26 plan; inox = inoxbas et inoxhaut ; inox = inox coul vert ; pyro1 = dall l60 l39 l31 l38 plan; pyro2 = dall l26 l60 l27 l59 plan; pyro3 = dall l21 l58 l23 l59 plan; pyrobas = dall l35 l33 l34 l29 plan; pyrohaut = dall l40 l41 l42 l30 plan; pyro = pyro1 et pyro2 et pyro3 et pyrobas et pyrohaut ; pyro = pyro coul bleu; BN1 = dall l51 l33 l52 l63 plan; BN2 = dall l50 l32 l49 l48 plan; BN3 = dall l44 l64 l43 l41 plan; BN = BN1 et BN2 et BN3; BN = BN coul rose; sample = dall l48 l45 l46 l47 plan; sample = sample coul roug; four = dall l31 l37 l62 l36 plan; four = four coul jaun; Cellule = inox et pyro et BN et sample et four ; elim 1e-5 Cellule ; *trac Cellule ; **************Ceinture Belt********* p42 = 6.2e-3 74.0e-3 ; p42 = cellule poin proc p42; p43 = 20.98e-3 56.5e-3 ; p44 = 25.0e-03 56.5e-3 ; p45 = 35.79e-03 56.5e-3 ; p46 = 58.85e-03 56.5e-3 ; p47 = 100.025e-03 56.5e-3 ; p48 = 100.025e-03 106.5e-03 ; p49 = 60.35e-03 106.5e-03 ; p50 = 37.29e-3 106.5e-03 ; p51 = 25.0e-3 106.5e-03 ; p52 = 20.98e-3 106.5e-03 ; p53 = 6.2e-03 89.0e-03 ; p53 = cellule poin proc p53; p54 = 10.5e-03 94.09135e-03 ; p55 = 13.0e-03 97.05145e-03 ; p56 = 10.5e-03 68.9086e-03 ; p57 = 13.0e-3 65.9484e-03 ; P85 = 10.e-3 163.e-3 ; p77 = 13e-3 97.05145e-3 ; p78 = 10.0e-3 99.664265e-3; p84 = 10.0e-3 63.33577e-3; l65 = p56 d 10 P46 ; l66 = p43 d 10 p44 ; l67 = p44 d 8 P45 ; l68 = p45 d 15 p46 ; l69 = p46 d 15 P47 ; l70 = p47 d 84 p48 ; l71 = p48 d 15 P49 ; l72 = p49 d 15 p50 ; l73 = p50 d 8 P51 ; l74 = p51 d 10 p52 ; l75 = p52 d 10 p55 ; l76 = p53 d 84 p42 ; l77 = p44 d 84 p51 ; l78 = p45 d 84 p50 ; l79 = p46 d 84 p49; l80 = p56 d 10 p57 ; l81 = p57 d 10 p43 ; l82 = p55 d 10 p54 ; l83 = p54 d 10 p53 ; l84 = p42 d 41 p57 ; l85 = p77 d 41 p53 ; l86 = p43 d 84 p52; l117 = p52 d 41 p77; l119 = p57 d 41 p43; l121 = p77 d 84 p57 ; c1 = dall l84 l121 l85 l76 plan ; c2 = dall l66 l77 l74 l86 plan ; c3 = dall l67 l78 l73 l77 plan ; c4 = dall l68 l79 l72 l78 plan ; c5 = dall l69 l70 l71 l79 plan ; c6 = dall l121 l119 l86 l117 plan; ceinture = c1 et c2 et c3 et c4 et c5 et c6; elim 1e-5 ceinture; *trac ceinture ; ************* Piston Haut ****** p58 = 0. 163.0e-3; p59 = 19.06e-03 163.0e-3; p60 = 61.56e-03 163.0e-3 ; p61 = 100.025e-03 163.0e-3 ; p62 = 100.025e-03 115.0e-3 ; p63 = 62.0e-03 115.0e-3 ; p64 = 20.05e-03 115.0e-3 ; p65 = 20.05e-03 112.0e-3 ; p66 = 0. 112.0e-3 ; p67 = 0. 115.0e-3 ; *p68 = 0. 95.0e-03 ; p69 = 5.3e-03 95.0e-03 ; *p70 = 6.2e-03 95.0e-03 ; p71 = 7.5e-03 99.664265e-03 ; p72 = 7.5e-03 96.59555e-03 ; p78 = 10.0e-3 99.664265e-3; P85 = 0. 99.664265e-3; nb1 = 48 ; nb2 = 36 ; l87 = p58 d nb1 P59 ; l88 = p59 d 50 p60 ; l89 = p60 d 20 P61 ; l90 = p61 d 20 p62 ; l91 = p62 d 20 P63 ; l92 = p63 d 50 p64 ; l93 = p64 d nb1 P67 ; l94 = p66 d nb1 p65 ; l95 = p59 d 20 P64 ; l96 = p58 d 20 p67 ; l97 = p64 d 20 p65 ; l98 = p67 d 20 p66 ; l99 = p21 d nb2 p85 ; l100 = p78 d nb2 p19 ; *l101 = l24 et l23 ; l101 = p19 d nb1 p21; l102 = p60 d 20 p63 ; l118 = p78 d 100 p65; l123 = p85 d nb1 p78 ; l124 = p66 d 100 p85; *l125 = p59 d 36 p85 ; e1 = dall l87 l95 l93 l96 plan; e2 = dall l97 l94 l98 l93 plan; e3 = dall l100 l101 l99 l123 plan; e4 = dall l88 l102 l92 l95 plan; e5 = dall l89 l90 l91 l102 plan; e6 = dall l123 l118 l94 l124 plan; Phaut = e1 et e2 et e3 et e4 et e5 et e6; elim 1e-5 Phaut; *trac Phaut; ***** Joints Pyro + Téflon (Haut et Bas) ************ *p73 = 6.2e-3 89.0e-3 ; p74 = 10.5e-3 94.09135e-3 ; p74 = ceinture poin proc p74; p75 = 7.5e-3 96.59555e-3 ; p75 = Phaut poin proc p75; *p76 = 6.2e-3 95.0e-3 ; p77 = 13e-3 97.05145e-3 ; p77 = ceinture poin proc p77; p78 = 10.0e-3 99.664265e-3; p78 = Phaut poin proc p78; l103 = p53 d 26 p74; l104 = p74 d 15 p77 ; l105 = p77 d 34 p78 ; l106 = p78 d 24 p75 ; l107 = p75 d 12 p19; l108 = p19 d 33 p53 ; l109 = p74 d 33 p75 ; jh1 = dall l103 l109 l107 l108 plan; jh1 = jh1 coul bleu; jh2 = dall l104 l105 l106 l109 plan; jh2 = jh2 coul turq; Jhaut = jh1 et jh2 ; elim 1e-5 Jhaut; *trac Jhaut; p57 = ceinture poin proc p57; *p79 = 6.2e-3 74.0e-3 ; *p80 = 10.5e-3 68.89086e-3 ; p56 = ceinture poin proc p56; p81 = 7.5e-3 66.4044e-3 ; p81 = Pbas poin proc p81; *p82 = 6.2e-3 68.0e-3 ; *p83 = 13e-3 65.9484e-3 ; p84 = 10.0e-3 63.33577e-3; p84 = Pbas poin proc p84; l110 = p42 d 26 p56; l111 = p56 d 15 p57 ; l112 = p57 d 34 p84 ; l113 = p84 d 24 p81 ; l114 = p81 d 12 p18; l115 = p18 d 33 p42 ; l116 = p56 d 33 p81 ; jb1 = dall l110 l116 l114 l115 plan; jb1 = jb1 coul bleu; jb2 = dall l111 l112 l113 l116 plan; jb2 = jb2 coul turq; Jbas = jb1 et jb2 ; elim 1e-5 Jbas; *trac Jbas; ********** Assemblage Maillages ********** elim 1e-5 Cellule Pbas; elim 1e-5 Cellule Jbas ; elim 1e-5 Cellule Jhaut; elim 1e-5 Cellule ceinture; elim 1e-5 Cellule Phaut; Tot = Cellule et Pbas et Phaut et ceinture et Jbas et Jhaut; elim 1e-5 Tot; Tot = rege Tot; *trac Tot ; ****************** MODELE ET MATERIAU NON-SOURCE (CONSTANT) MOD1 = MODEL Pbas THERMIQUE ISOTROPE; MOD2 = MODEL Phaut THERMIQUE ISOTROPE; MOD3 = MODEL Jb1 THERMIQUE ISOTROPE; MOD4 = MODEL Jb2 THERMIQUE ISOTROPE; MOD5 = MODEL Jh1 THERMIQUE ISOTROPE; MOD6 = MODEL Jh2 THERMIQUE ISOTROPE; MOD7 = MODEL ceinture THERMIQUE ISOTROPE; MOD8 = MODEL BN THERMIQUE ISOTROPE; MOD9 = MODEL four THERMIQUE ISOTROPE; *MOD9 = REDU MOD9 FOUR ; MOD10 = MODEL sample THERMIQUE ISOTROPE; MOD11 = MODEL pyro THERMIQUE ISOTROPE; MOD12 = MODEL inox THERMIQUE ISOTROPE; MAT1 = MATE MOD1 K 110. ; MAT2 = MATE MOD2 K 110. ; MAT3 = MATE MOD3 K 2.65 ; MAT4 = MATE MOD4 K 0.25 ; MAT5 = MATE MOD5 K 2.65 ; MAT6 = MATE MOD6 K 0.25 ; MAT7 = MATE MOD7 K 110. ; MAT8 = MATE MOD8 K 24. ; ***************MISE EN PLACE VARIATION K DANS MATERIAU SOURCE******** *evk = evol MANU 'T' (prog 20 3000) 'K' (prog 44 4) ; evk = evol MANU ' Temperature ' (prog 0. 100. 200. 300. 400. 500. 600. 700. 800. 900. 1000. 1200. 1600. 2000. 3000. ) 'Conductivite Four' (prog 44. 39.5 35.64 31.43 27.66 24.51 22.25 20.08 17.98 15.19 14.37 13.2 11.4 9.9 8.) ; *dess evk ; MAT9 = MATE MOD9 K evk; **** MAT10 = MATE MOD10 K 60. ; MAT11 = MATE MOD11 K 2.65; MAT12 = MATE MOD12 K 26.; ************* TEMPERATURE IMPOSEE (CONDITIONS DE BORD Dirichlet)****** *B1 = BLOQUE T l66; *B2 = BLOQUE T l67; *B3 = BLOQUE T l68; *B4 = BLOQUE T l69; *B5 = BLOQUE T l70; *B6 = BLOQUE T l71; *B7 = BLOQUE T l72; *B8 = BLOQUE T l73; *B9 = BLOQUE T l74; *B10 = BLOQUE T l117; *B11= BLOQUE T l119; B12 = BLOQUE T l112; B13 = BLOQUE T l105; B14 = BLOQUE T l5; B15 = BLOQUE T l6; B16 = BLOQUE T l11; B17= BLOQUE T l120; B18 = BLOQUE T l91; B19 = BLOQUE T l92; B20 = BLOQUE T l97; B21= BLOQUE T l118; *B = B1 et B2 et B3 et B4 et B5 et B6 et B7 et B8 et B9 et B10 et B11 *et B12 et B13 et B14 et B15 et B16 et B17 et B18 et B19 et B20 et B21; Bb = BLOQUE T ceinture; B = Bb et B12 et B13 et B14 et B15 et B16 et B17 et B18 et B19 et B20 et B21; E1 = DEPI B 20 ; ************ FLUX NUL (CONDITIONS DE BORD Neumann)**************** *F1 = FLUX MOD4 0. l112 ; *F2 = FLUX MOD6 0. l105 ; F3 = FLUX MOD1 0. l1 ; F4 = FLUX MOD1 0. l2 ; F5 = FLUX MOD1 0. l3 ; F6 = FLUX MOD1 0. l4 ; *F7 = FLUX MOD1 0. l5 ; *F8 = FLUX MOD1 0. l6 ; *F9 = FLUX MOD1 0. l11 ; *F10 = FLUX MOD1 0. l120 ; F11 = FLUX MOD2 0. l87 ; F12 = FLUX MOD2 0. l88 ; F13 = FLUX MOD2 0. l89 ; F14 = FLUX MOD2 0. l90 ; *F15 = FLUX MOD2 0. l91 ; *F16 = FLUX MOD2 0. l92 ; *F17 = FLUX MOD2 0. l97 ; *F18 = FLUX MOD2 0. l118 ; ********************SOURCE VOLUMIQUE DE CHALEUR******************** **********Calcul Puissance Volumique ******* Données **** *** Dimension Four (m) d_int = 2.5e-3 ; d_ext = 3.e-3 ; h = 15.2e-3 ; *** Puissance électrique (W) P = 565. ; de = d_ext * d_ext ; di = d_int * d_int; df = de - di; V = pi * df * h ; Pv = P / V ; MODZ = REDU MOD9 FOUR ; P = SOURCE MODZ Pv four ; *P = REDU P FOUR; Charge = F3 et F4 et F5 et F6 et F11 et F12 et F13 et F14 et P; *Charge = F1 et F2 et F3 et F4 et F5 et F6 et F7 et F8 et F9 et F10 et *F11 et F12 et F13 et F14 et F15 et F16 et F17 et F18 et P; *Charge = P; ******************Résolution Non-Linéaire**************************** ******** Préparation de table pour THERMIC tabth = table THERMIQUE ; tabth . 'BLOCAGE' = b ; tabth . 'IMPOSE' = E1 ; tabth . 'FLUX' = Charge ; tabth . 'INSTANT(0)' = manu chpo Tot 1 T 20 nature diffus ; tabth . 'TABCOND' = table ; tabth . 'TABCOND' . MOD1 = 110. ; tabth . 'TABCOND' . MOD2 = 110.; tabth . 'TABCOND' . MOD3 = 2.65 ; tabth . 'TABCOND' . MOD4 = 0.25 ; tabth . 'TABCOND' . MOD5 = 2.65 ; tabth . 'TABCOND' . MOD6 = 0.25 ; tabth . 'TABCOND' . MOD7 = 110.; tabth . 'TABCOND' . MOD8 = 24. ; tabth . 'TABCOND' . MOD9 = evk; tabth . 'TABCOND' . MOD10 = 60. ; tabth . 'TABCOND' . MOD11 = 2.65 ; tabth . 'TABCOND' . MOD12 = 26. ; tabth . 'NIVEAU' = 1 ; tabth . 'CRITERE'= 1e-5 ; *** Lancement du calcul ...0 thermic tabth NONLINEAIRE ; ****************************POST TRAITEMENT tresu = tabth . TEMPERATURE ; trac sample tresu; T1 = extr tabth . TEMPERATURE T p36; evol1 = EVOL vert CHPO tresu T l46; evol2 = EVOL rouge CHPO tresu T l48; tab1 = TABLE ; tab1 . 'TITRE' = TABLE ; tab1 . 'TITRE' . 1 = 'Ligne centrale'; tab1 . 'TITRE' . 2 = 'Ligne externe'; dess (evol1 et evol2) MIMA TITR 'Temperature echantillon le long de lignes (surf) axiales' LEGE tab1 ; evol3 = EVOL vert CHPO tresu T l45; evol4 = EVOL rouge CHPO tresu T l47; tab2 = TABLE ; tab2 . 'TITRE' = TABLE ; tab2 . 'TITRE' . 1 = 'Ligne superieure' ; tab2 . 'TITRE' . 2 = 'Ligne inferieure'; dess (evol3 et evol4) MIMA TITR 'Temperature echantillon le long de lignes (surf.) radiales' LEGE tab2 ; *fin;