* ************************************************************ * * * Maillage du Four HEC avec navette * * * ************************************************************ * symétrie plan * on place de la poudre en bas * * x est l'axe du four * * résistance de contact entre poudre et navette * mais noeuds communs gaz/métal et gaz/poudre * OPTION DIME 3 ELEM CUB8 PRI6; OPTIO SAUV 'Four_HEC.mgibi'; * OPTIO IMPI 1; DENSITE 0.001 ; * oeil = 25. -50. 50. ; oeilx = -50. 0. 0. ; oeily = 0 -50. 0. ; oeilz = 0. 0. 50. ; * x00 = 0. 0. 0. ; x01 = 0. 1. 0. ; x02 = 0. 0. 1. ; * parametres du maillage * toutes les dimensions en m * * Epaisseur du métal de la navette, partout * valeur navette n°4 Ep_nav = 0.0015 ; * * Longueur intérieure de la navette * valeur navette n°4 * on suppose la navette bien centrée Lo_nav = 0.150 ; Lo_poud = Lo_nav - (2*Ep_nav); * * Diamètre interne tube métallique Dia_tubm = 0.071 ; * * Epaisseur du tube métallique Ep_tub = 0.002 ; * * attention le tube n'est pas monté de manière symétrique dans le four * Longueur du tube métallique L_tub = 1.600 ; L_t1 = 0.29785 ; L_t2 = 0.31215 ; * * Longueur extérieure du capot métallique L_capot = 0.990; ep_ref = 0.009 ; * * Rayon exterieur du capot métallique R_capot = 0.264 ; * * Longueur de l'isolant L_isol = L_capot -(2 * ep_ref); * * Rayon externe de la couche d'isolant R_isol = R_capot - ep_ref; * * Longueur totale des résistors * A noter que l'on a 3 résistors bout à bout L_resi = 0.710 ; * * Epaisseur du résistor Ep_resi = 0.01 ; * * Diamètre ouverture de la navette pour laisser passer le gaz * valeur arbitraire Dia_navo = 0.030 ; * * Hauteur du lit de poudre * selon ET, entre 30mm et 7mm H_poud = 0.010 ; * * calcul de grandeur dérivées Di_poud = Dia_tubm - (2. * Ep_nav); Ri_poud = Di_poud / 2. ; Ri_tubm = Dia_tubm / 2. ; Ria_navo = Dia_navo / 2 ; * * paramétrisation du nombre d'éléments nel_tole = 1; nel_poud = 15; * * dessin des génératrices lg0 = (((-1./2.)*L_capot) - L_t1) 0. 0.; l_tub1 = L_t1 ; lg1 = lg0 PLUS (l_tub1 0. 0.); l_tub2 = ep_ref ; lg2 = lg1 PLUS (l_tub2 0. 0.); l_tub3 = (L_isol- L_resi)/ 2. ; lg3 = lg2 PLUS (l_tub3 0. 0.); l_tub4 = L_resi / 3. ; lg4 = lg3 PLUS (l_tub4 0. 0.); lg41 = ((-1./2.)* Lo_nav) 0. 0.; lg42 = lg41 PLUS (Ep_nav 0. 0.); lg43 = ((0.5* Lo_nav) - Ep_nav) 0. 0.; lg44 = lg43 PLUS (Ep_nav 0. 0.); lg5 = lg4 PLUS (l_tub4 0. 0.); lg6 = lg5 PLUS (l_tub4 0. 0.); lg7 = (L_isol / 2.) 0. 0.; lg8 = lg7 PLUS (l_tub2 0. 0.); l_tub5 = L_t2; lg9 = lg8 PLUS (l_tub5 0. 0.); * lgen_n1 = (D nel_tole lg41 lg42); lgen_n2 = (D nel_poud lg42 lg43); lgen_n3 = (D nel_tole lg43 lg44); lgen_nav = lgen_n1 ET lgen_n2 ET lgen_n3; * lgen_r1 = (D 6 lg3 lg4); lgen_r21 = (D 2 lg4 lg41); lgen_r22 = (D 2 lg44 lg5); lgen_r2 = lgen_r21 ET lgen_nav ET lgen_r22 ; lgen_r3 = (D 6 lg5 lg6); lgen_res = lgen_r1 ET lgen_r2 ET lgen_r3 ; ELIM 1.E-9 lgen_res; * lgen_is1 = (D 6 lg2 lg3); lgen_is2 = (D 6 lg6 lg7); lgen_iso = lgen_is1 ET lgen_r1 ET lgen_r2 ET lgen_r3 ET lgen_is2; ELIM 1.E-9 lgen_iso; * lgen_cp1 = (D 4 lg1 lg2) ; lgen_cp2 = (D 4 lg7 lg8) ; lgen_cpm = lgen_cp1 ET lgen_is1 ET lgen_r1 ET lgen_r2 ET lgen_r3 ET lgen_is2 ET lgen_cp2 ; ELIM 1.E-9 lgen_cpm; * lgen_t1 = (D 4 lg0 lg1); lgen_t2 = (D 4 lg8 lg9); lgen_tf1 = lgen_t1 ET (D 1 lg1 lg2); lgen_tf2 = lgen_t2 ET (D 1 lg7 lg8); * lgen_tub = lgen_tf1 ET lgen_is1 ET lgen_r1 ET lgen_r2 ET lgen_r3 ET lgen_is2 ET lgen_tf2; ELIM 1.E-9 lgen_tub; * lgen_gzg = lgen_tf1 ET lgen_is1 ET lgen_r1 ET lgen_r21; ELIM 1.E-9 lgen_gzg; lgen_gzd = lgen_r22 ET lgen_r3 ET lgen_is2 ET lgen_tf2; ELIM 1.E-9 lgen_gzd; * lgen_res lgen_r1 lgen_r2 lgen_r3 = lgen_res lgen_r1 lgen_r2 lgen_r3 PLUS (0. 0. (Ri_tubm + Ep_tub)); lgen_iso = lgen_iso PLUS (0. 0. (Ri_tubm + Ep_tub + Ep_resi)); lgen_is1 = lgen_is1 PLUS (0. 0. (Ri_tubm + Ep_tub + Ep_resi)); lgen_cpm = lgen_cpm PLUS (0. 0. R_capot); lgen_tub = lgen_tub PLUS (0. 0. Ri_tubm); lgen_t1 = lgen_t1 PLUS (0. 0. (Ri_tubm + ep_tub)); lgen_t2 = lgen_t2 PLUS (0. 0. (Ri_tubm + ep_tub)); lgen_gzg = lgen_gzg PLUS (0. 0. Ri_tubm); lgen_gzd = lgen_gzd PLUS (0. 0. Ri_tubm); * * dessin de la poudre yb0 = Ri_poud*(SIN (ACOS ((Ri_poud-H_poud)/Ri_poud))); a00 = ((-1.)*(Lo_poud/2)) 0. 0. ; p01 = ((-1.)*(Lo_poud/2)) 0. (H_poud-Ri_poud) ; p02 = ((-1.)*(Lo_poud/2)) yb0 (H_poud-Ri_poud) ; p03 = ((-1.)*(Lo_poud/2)) 0. ((-1.)*Ri_poud) ; l_ver = D 4 p01 p03 ; l_arc = CERCLE 7 p03 a00 p02 ; l_hor = D 6 p01 p02 ; ELIM 1.E-9 (l_ver ET l_arc ET l_hor); * surface laterale de la poudre, à droite s_pou_d = SURF PLANE (l_ver ET l_arc ET l_hor); poudre = VOLU TRAN s_pou_d nel_poud (Lo_poud 0. 0.); * surface ronde en contact avec la navette s_pou_h = TRAN l_arc nel_poud (Lo_poud 0. 0.); ELIM 1.E-9 s_pou_h poudre; * surface laterale de la poudre, à gauche s_pou_g = FACE 2 poudre; * pour compatibilité avec maillage four R4 s_pouX = s_pou_g; s_pouX2 = s_pou_d; * surface du haut s_pou_r = TRAN l_hor nel_poud (Lo_poud 0. 0.); s_pou_r = ORIENTER s_pou_r (0. 0. 1.) ; ELIM 1.E-9 s_pou_r poudre; * surface en symétrie s_pou_s = TRAN l_ver nel_poud (Lo_poud 0. 0.); ELIM 1.E-9 s_pou_s poudre; * poudre = poudre COUL VERT ; * * dessin de la navette * partie cylindrique en contact avec la poudre yb1 = Ri_tubm*(SIN (ACOS ((Ri_poud-H_poud)/Ri_tubm))); a00 = ((-1.)*(Lo_poud/2)) 0. 0. ; p01 = ((-1.)*(Lo_poud/2)) 0. ((-1.)*Ri_poud) ; p02 = ((-1.)*(Lo_poud/2)) 0. ((-1.)*Ri_tubm) ; p03 = ((-1.)*(Lo_poud/2)) yb1 (H_poud-Ri_poud) ; p04 = ((-1.)*(Lo_poud/2)) yb0 (H_poud-Ri_poud) ; l_ver = D nel_tole p01 p02; l_arc1 = CERCLE 5 p02 a00 p03; l_hor = D nel_tole p03 p04; l_arc2 = CERCLE 7 p04 a00 p01; ELIM 1.E-9 (l_ver ET l_arc1 ET l_hor ET l_arc2); s_nav1 = SURF PLANE (l_ver ET l_arc1 ET l_hor ET l_arc2); * nav1 = VOLU TRAN s_nav1 nel_poud (Lo_poud 0. 0.); * surface du haut pour contact avec poudre s_fon_ph = TRAN l_arc2 nel_poud (Lo_poud 0. 0.); ELIM 1.E-9 s_fon_ph nav1; * * fond troué droit et gauche, partie basse tubulaire nav1g = VOLU TRAN s_nav1 nel_tole (((-1.)*Ep_nav) 0. 0.); nav1d = nav1g PLUS ((Lo_poud + Ep_nav) 0. 0.); * surfaces latérales, pour échange avec gaz et rayonnement s_nav1g = FACE 2 nav1g; s_nav1d = FACE 1 nav1d; s_nav1g = ORIENTER s_nav1g (-1. 0. 0.); s_nav1d = ORIENTER s_nav1d (1. 0. 0.); * * partie cylindrique haute de la navette a00 = ((-1.)*(Lo_poud/2)) 0. 0. ; p01 = ((-1.)*(Lo_poud/2)) yb0 (H_poud-Ri_poud) ; p02 = ((-1.)*(Lo_poud/2)) Ri_poud 0. ; p03 = ((-1.)*(Lo_poud/2)) 0. Ri_poud ; p04 = ((-1.)*(Lo_poud/2)) 0. Ri_tubm ; p05 = ((-1.)*(Lo_poud/2)) Ri_tubm 0. ; p06 = ((-1.)*(Lo_poud/2)) yb1 (H_poud-Ri_poud) ; l_arc1 = CERCLE 5 p01 a00 p02 ; l_arc2 = CERCLE 5 p02 a00 p03 ; l_ver = D nel_tole p03 p04; l_arc3 = CERCLE 5 p04 a00 p05 ; l_arc4 = CERCLE 5 p05 a00 p06 ; l_hor = D nel_tole p06 p01; ELIM 1.E-9 (l_arc1 ET l_arc2 ET l_ver ET l_arc3 ET l_arc4 ET l_hor ); s_nav2 = SURF PLANE (l_arc1 ET l_arc2 ET l_ver ET l_arc3 ET l_arc4 ET l_hor ); * nav2 = VOLU TRAN s_nav2 nel_poud (Lo_poud 0. 0.); * surface intérieure qui échange avec le gaz, et qui rayonne s_nav2m = TRAN (l_arc1 ET l_arc2) nel_poud (Lo_poud 0. 0.); s_nav2m = ORIENTER s_nav2m (0. -1. 0.); * * fond lateral droit et gauche de la navette en contact avec poudre s_nav3 = s_pou_d PLUS (0. 0. 0.); nav3g = VOLU TRAN s_nav3 nel_tole (((-1.)*Ep_nav) 0. 0.); nav3d = VOLU TRAN (s_nav3 PLUS (Lo_poud 0. 0.)) nel_tole (Ep_nav 0. 0.); * surfaces interieures qui échangent avec la poudre s_nav3pg = nav3g FACE 1; s_nav3pd = nav3d FACE 1; * surfaces exterieures qui échangent avec le gaz, et qui rayonnent s_nav3g = nav3g FACE 2; s_nav3d = nav3d FACE 2; s_nav3g = ORIENTER s_nav3g (-1. 0. 0.); s_nav3d = ORIENTER s_nav3d (1. 0. 0.); * * fond lateral droit et gauche en contact avec cylindre nav4g = VOLU TRAN s_nav2 nel_tole (((-1.)*Ep_nav) 0. 0.); nav4d = VOLU TRAN (s_nav2 PLUS (Lo_poud 0. 0.)) nel_tole (Ep_nav 0. 0.); * surfaces exterieures qui échangent avec le gaz, et qui rayonnent s_nav4g = nav4g FACE 2; s_nav4d = nav4d FACE 2; s_nav4g = ORIENTER s_nav4g (-1. 0. 0.); s_nav4d = ORIENTER s_nav4d (1. 0. 0.); * * partie generale des fonds latéraux yb0 = Ri_poud*(SIN (ACOS ((Ri_poud-H_poud)/Ri_poud))); a00 = ((-1.)*(Lo_poud/2)) 0. 0. ; p01 = ((-1.)*(Lo_poud/2)) yb0 (H_poud-Ri_poud) ; p02 = ((-1.)*(Lo_poud/2)) 0. (H_poud-Ri_poud) ; p03= ((-1.)*(Lo_poud/2)) 0. ((-1.)*Ria_navo) ; p04= ((-1.)*(Lo_poud/2)) Ria_navo 0. ; p05= ((-1.)*(Lo_poud/2)) 0. Ria_navo ; p06= ((-1.)*(Lo_poud/2)) 0. Ri_poud ; p07= ((-1.)*(Lo_poud/2)) Ri_poud 0. ; l_hor = D 6 p01 p02 ; l_ver1 = D 4 p02 p03 ; l_arc1 = CERCLE 4 p03 a00 p04; l_arc2 = CERCLE 4 p04 a00 p05; l_ver2 = D 4 p05 p06; l_arc3 = CERCLE 5 p06 a00 p07; l_arc4 = CERCLE 5 p07 a00 p01; ELIM 1.E-9 (l_hor ET l_ver1 ET l_arc1 ET l_arc2 ET l_ver2 ET l_arc3 ET l_arc4); s_nav5 = (l_hor ET l_ver1 ET l_arc1 ET l_arc2 ET l_ver2 ET l_arc3 ET l_arc4) SURF PLANE; * nav5g = VOLU TRAN s_nav5 nel_tole (((-1.)*Ep_nav) 0. 0.); nav5d = VOLU TRAN (s_nav5 PLUS (Lo_poud 0. 0.)) nel_tole (Ep_nav 0. 0.); * surfaces interieures qui échangent avec le gaz, et qui rayonnent s_nav5ig = nav5g FACE 1; s_nav5id = nav5d FACE 1; s_nav5ig = ORIENTER s_nav5ig (-1. 0. 0.); s_nav5id = ORIENTER s_nav5id (1. 0. 0.); * surfaces dans les ouvertures qui échangent avec le gaz, et qui rayonnent s_nav5og = TRAN (l_arc1 ET l_arc2) nel_tole (((-1.)*Ep_nav) 0. 0.); ELIM 1.E-9 s_nav5og nav5g; s_nav5od = TRAN ((l_arc1 ET l_arc2) PLUS (Lo_poud 0. 0.)) nel_tole (Ep_nav 0. 0.); ELIM 1.E-9 s_nav5od nav5d; s_nav5og = ORIENTER s_nav5og (0. -1. 0.); s_nav5od = ORIENTER s_nav5od (0. -1. 0.); * surfaces exterieures qui échangent avec le gaz, et qui rayonnent s_nav5g = nav5g FACE 2; s_nav5d = nav5d FACE 2; s_nav5g = ORIENTER s_nav5g (-1. 0. 0.); s_nav5d = ORIENTER s_nav5d (1. 0. 0.); * navette = nav1 ET nav1g ET nav1d ET nav2 ET nav3g ET nav3d ET nav4g ET nav4d ET nav5g ET nav5d ; ELIM 1.E-9 navette; navette = navette COUL BLEU ; * * dessin de l'espace occupé par le gaz * gaz dans la navette, hors noyau de ventilation s_gaz1 lg_hor lg_arc3 lg_arc4 lg_ver1 lg_ver2 = s_nav5 l_hor l_arc3 l_arc4 l_ver1 l_ver2 PLUS (0. 0. 0.); gaz1 = VOLU TRAN s_gaz1 nel_poud (Lo_poud 0. 0.) ; * surfaces pour échange convectif s_gng = gaz1 FACE 1; s_gnd = gaz1 FACE 2; s_gni = gaz1 FACE 3; s_gn_h = TRAN lg_hor nel_poud (Lo_poud 0. 0.); s_gn_l = TRAN (lg_arc3 ET lg_arc4) nel_poud (Lo_poud 0. 0.); ELIM 1.E-9 gaz1 (s_gn_h ET s_gn_l); * surface en symétrie s_gn_s = TRAN (lg_ver1 ET lg_ver2) nel_poud (Lo_poud 0. 0.); ELIM 1.E-9 s_gn_s gaz1; * * gaz dans la navette, dans le noyau de ventilation a00 = ((-1.)*(Lo_poud/2)) 0. 0. ; p01= ((-1.)*(Lo_poud/2)) 0. ((-1.)*Ria_navo) ; p02= ((-1.)*(Lo_poud/2)) Ria_navo 0. ; p03= ((-1.)*(Lo_poud/2)) 0. Ria_navo ; l_arc = (CERCLE 4 p01 a00 p02)ET (CERCLE 4 p02 a00 p03); l_ver = D 6 p03 p01 ; s_gaz2 = SURF PLANE (l_arc ET l_ver); gaz2 = VOLU TRAN s_gaz2 nel_poud (Lo_poud 0. 0.) ; * * gaz dans les ouvertures de la navette gaz3g = VOLU TRAN s_gaz2 nel_tole (((-1.)*Ep_nav) 0. 0.); gaz3d = VOLU TRAN (s_gaz2 PLUS (Lo_poud 0. 0.)) nel_tole (Ep_nav 0. 0.); * surfaces pour echange avec navette s_gaz3g = TRAN l_arc nel_tole (((-1.)*Ep_nav) 0. 0.); ELIM 1.E-9 s_gaz3g gaz3g ; s_gaz3d = TRAN (l_arc PLUS (Lo_poud 0. 0.)) nel_tole (Ep_nav 0. 0.); ELIM 1.E-9 s_gaz3d gaz3d ; * gaz = gaz1 ET gaz2 ET gaz3g ET gaz3d ; ELIM 1.E-9 gaz ; gaz = gaz COUL GRIS ; * * dessin du tube métallique a00 = (((-1./2.)*L_capot) - L_t1) 0. 0. ; p01 = (((-1./2.)*L_capot) - L_t1) 0. ((-1.)*Ri_tubm) ; p02 = (((-1./2.)*L_capot) - L_t1) 0. ((-1.)*(Ri_tubm + Ep_tub)); p03 = (((-1./2.)*L_capot) - L_t1) (Ri_tubm + Ep_tub) 0. ; p04 = (((-1./2.)*L_capot) - L_t1) 0. (Ri_tubm + Ep_tub); p05 = (((-1./2.)*L_capot) - L_t1) 0. Ri_tubm ; p06 = (((-1./2.)*L_capot) - L_t1) Ri_tubm 0. ; p07 = (((-1./2.)*L_capot) - L_t1) 0. 0. ; m01 = (Lo_nav / 2.) 0. ((-1.)*Ri_tubm) ; m02 = (Lo_nav / 2.) 0. ((-1.)*(Ri_tubm + Ep_tub)); m03 = (Lo_nav / 2.) (Ri_tubm + Ep_tub) 0. ; l_ver1 = D nel_tole p01 p02 ; l_arc1 = CERCLE 5 p02 a00 p03 ; l_arc2 = CERCLE 5 p03 a00 p04 ; l_ver2 = D nel_tole p04 p05 ; l_arc3 = CERCLE 5 p05 a00 p06; l_arc4 = CERCLE 5 p06 a00 p01; l_ver3 = (D 4 p05 p07) ET (D 4 p07 p01); s_tub = SURF PLANE (l_ver1 ET l_arc1 ET l_arc2 ET l_ver2 ET l_arc3 ET l_arc4); * tube = s_tub VOLU 'GENE' lgen_tub; tube = tube COUL BLEU ; * * création de la surface extérieure du tube refroidie par air stub_a1 = (l_arc1 ET l_arc2) GENE lgen_tf1 ; stub_a2 = ((l_arc1 ET l_arc2)PLUS ((L_t1 + ep_ref + L_isol) 0. 0.)) GENE lgen_tf2 ; stub_air = stub_a1 ET stub_a2 ; ELIM 1.E-9 stub_air tube ; * * creation des surfaces internes du tube, rayonnantes et en contact avec gaz larcr = (l_arc3 ET l_arc4) PLUS (0. 0. 0.); s_tug1 = larcr GENE lgen_gzg; s_tug2 = (larcr PROJ 'CYLI' (-1. 0. 0.) 'PLAN' m01 m02 m03) GENE lgen_gzd; s_tug = s_tug1 ET s_tug2 ; s_tug = s_tug COUL BLEU; s_tug = s_tug ORIENTER (0. -1. 0.); * rajout des bouts du tube, partie cylindrique s_tub = SURF PLANE (l_ver1 ET l_arc1 ET l_arc2 ET l_ver2 ET l_arc3 ET l_arc4); bt1_t_g = VOLU TRAN s_tub 1 (((-1.)*Ep_tub) 0. 0.); sbr_frg = TRAN (l_arc1 ET l_arc2) nel_tole (((-1.)*Ep_tub) 0. 0.); bt1_t_d = bt1_t_g PLUS ((L_tub + Ep_tub) 0. 0.); sbr_frd = sbr_frg PLUS ((L_tub + Ep_tub) 0. 0.); tube = tube ET bt1_t_d ET bt1_t_g ; ELIM 1.E-9 tube ; tube = tube COUL BLEU ; * * création des surfaces des brides refroidies ELIM 1.E-9 (sbr_frg ET sbr_frd) tube ; * * rajout de la partie courante des bouts du tube s_bt_1 = ((s_gaz2 MOINS (Ep_nav 0. 0.)) ET s_nav1g ET s_nav3g ET s_nav4g ET s_nav5g) PLUS (0. 0. 0.); s_bt_1 = s_bt_1 PROJ 'CYLI' (-1. 0. 0.) 'PLAN' p01 p02 p03 ; bt2_t_d = VOLU TRAN s_bt_1 nel_tole (((-1.)*Ep_tub) 0. 0.); bt2_t_g = bt2_t_d PLUS ((L_tub + Ep_tub) 0. 0.); s_bt_2 = s_bt_1 PLUS ((L_tub) 0. 0.); ELIM 1.E-9 s_bt_2 bt2_t_g; s_bt_1 = ORIENTER s_bt_1 (1. 0. 0.); s_bt_2 = ORIENTER s_bt_2 (-1. 0. 0.); tube = tube ET bt2_t_d ET bt2_t_g ; ELIM 1.E-9 tube ; tube = tube COUL BLEU ; * * dessin du gaz à l'intérieur du tube s_gz_tub = s_nav1g ET s_nav3g ET s_nav4g ET s_nav5g ; ELIM 1.E-9 s_gz_tub ; s_gz_tub = s_gz_tub PROJ 'CYLI' (-1. 0. 0.) 'PLAN' p01 p02 p03 ; atm_tub1 = s_gz_tub VOLU 'GENE' lgen_gzg ; s_gz_tub = s_gz_tub PROJ 'CYLI' (1. 0. 0.) 'PLAN' m01 m02 m03 ; atm_tub2 = s_gz_tub VOLU 'GENE' lgen_gzd ; * s_gz_noy = s_gaz2 MOINS (Ep_nav 0. 0.); s_gz_noy = s_gz_noy PROJ 'CYLI' (-1. 0. 0.) 'PLAN' p01 p02 p03 ; atm_noy1 = s_gz_noy VOLU 'GENE' lgen_gzg ; s_gz_noy = s_gz_noy PROJ 'CYLI' (1. 0. 0.) 'PLAN' m01 m02 m03 ; atm_noy2 = s_gz_noy VOLU 'GENE' lgen_gzd ; gaz = gaz ET atm_tub1 ET atm_tub2 ET atm_noy1 ET atm_noy2 ; ELIM 1.E-9 gaz ; gaz = gaz COUL GRIS ; atm_noy = gaz2 ET gaz3g ET gaz3d ET atm_noy1 ET atm_noy2; * * surfaces pour échange avec la navette s_gz_gm = (s_nav1g ET s_nav3g ET s_nav4g ET s_nav5g) PLUS (0. 0. 0.); s_gz_dm = (s_nav1g ET s_nav3g ET s_nav4g ET s_nav5g) PLUS ((Lo_poud + Ep_nav + Ep_nav) 0. 0.); ELIM 1.E-9 (s_gz_gm ET s_gz_dm) gaz; * * surfaces planes en extrémité pour injecter et faire sortir le gaz procédé * attention dans cette version * on considère que le gaz ne circule que dans le noyau de ventilation s_gz_in = s_gaz2 PROJ 'CYLI' (-1. 0. 0.) 'PLAN' p01 p02 p03 ; s_gz_out = s_gz_in PLUS (L_tub 0. 0.); ELIM 1.E-9 (s_gz_in ET s_gz_out) gaz ; * * dessin du tube resistor a00 = ((-1./2.)*L_resi) 0. 0. ; p01 = ((-1./2.)*L_resi) 0. ((-1.)*(Ri_tubm + Ep_tub)) ; p02 = ((-1./2.)*L_resi) 0. ((-1.)*(Ri_tubm + Ep_tub + Ep_resi)); p03 = ((-1./2.)*L_resi) (Ri_tubm + Ep_tub + Ep_resi) 0. ; p04 = ((-1./2.)*L_resi) 0. (Ri_tubm + Ep_tub + Ep_resi) ; p05 = ((-1./2.)*L_resi) 0. (Ri_tubm + Ep_tub) ; p06 = ((-1./2.)*L_resi) (Ri_tubm + Ep_tub) 0. ; l_ver1 = D 3 p01 p02 ; l_arc1 = CERCLE 5 p02 a00 p03 ; l_arc2 = CERCLE 5 p03 a00 p04 ; l_ver2 = D 3 p04 p05 ; l_arc3 = CERCLE 5 p05 a00 p06; l_arc4 = CERCLE 5 p06 a00 p01; s_resi = SURF PLANE (l_ver1 ET l_arc1 ET l_arc2 ET l_ver2 ET l_arc3 ET l_arc4); tub_res1 = s_resi VOLU 'GENE' lgen_r1; tub_res2 = (s_resi PLUS ((L_resi / 3.) 0. 0.)) VOLU 'GENE' lgen_r2; tub_res3 = (s_resi PLUS ((2. * L_resi / 3.) 0. 0.)) VOLU 'GENE' lgen_r3; tub_resi = tub_res1 ET tub_res2 ET tub_res3; ELIM 1.E-9 tub_resi ; tub_resi = tub_resi COUL ROUGE ; * * dessin de l'isolant * partie basse de l'isolant, en contact avec tube, à coté des résistors a00 = ((-1./2.)*L_isol) 0. 0. ; p01 = ((-1./2.)*L_isol) 0. ((-1.)*(Ri_tubm + Ep_tub)) ; p02 = ((-1./2.)*L_isol) 0. ((-1.)*(Ri_tubm + Ep_tub + Ep_resi)) ; p03 = ((-1./2.)*L_isol) (Ri_tubm + Ep_tub + Ep_resi) 0. ; p04 = ((-1./2.)*L_isol) 0. (Ri_tubm + Ep_tub + Ep_resi) ; p05 = ((-1./2.)*L_isol) 0. (Ri_tubm + Ep_tub) ; p06 = ((-1./2.)*L_isol) (Ri_tubm + Ep_tub) 0. ; l_ver1 = D 3 p01 p02 ; l_arc1 = CERCLE 5 p02 a00 p03 ; l_arc2 = CERCLE 5 p03 a00 p04 ; l_ver2 = D 3 p04 p05 ; l_arc3 = CERCLE 5 p05 a00 p06; l_arc4 = CERCLE 5 p06 a00 p01; s_isolb = SURF PLANE (l_ver1 ET l_arc1 ET l_arc2 ET l_ver2 ET l_arc3 ET l_arc4); l_isol1 = (L_isol- L_resi)/ 2. ; isolb1 = s_isolb VOLU 'GENE' lgen_is1 ; isolb11 = isolb1 SYME 'PLAN' x00 x01 x02 ; isolant = isolb1 ET isolb11 ; * partie haute de l'isolant, qui englobe les résistors p01 = ((-1./2.)*L_isol) 0. ((-1.)*(Ri_tubm + Ep_tub + Ep_resi)) ; p02 = ((-1./2.)*L_isol) 0. ((-1.)* R_isol) ; p03 = ((-1./2.)*L_isol) R_isol 0. ; p04 = ((-1./2.)*L_isol) 0. R_isol ; p05 = ((-1./2.)*L_isol) 0. (Ri_tubm + Ep_tub + Ep_resi) ; p06 = ((-1./2.)*L_isol) (Ri_tubm + Ep_tub + Ep_resi) 0. ; l_ver1 = D 7 p01 p02 ; l_arc1 = CERCLE 5 p02 a00 p03 ; l_arc2 = CERCLE 5 p03 a00 p04 ; l_ver2 = D 7 p04 p05 ; l_arc3 = CERCLE 5 p05 a00 p06; l_arc4 = CERCLE 5 p06 a00 p01; s_isol = SURF PLANE (l_ver1 ET l_arc1 ET l_arc2 ET l_ver2 ET l_arc3 ET l_arc4); isolant = isolant ET (s_isol VOLU 'GENE' lgen_iso); isolant = isolant COUL JAUNE; ELIM 1.E-9 isolant ; * * dessin du capot métallique * partie basse du capot, à proximité du tube métallique a00 = ((-1./2.)*L_capot) 0. 0. ; p01 = ((-1./2.)*L_capot) 0. ((-1.)*(Ri_tubm + Ep_tub + Ep_resi)) ; p02 = ((-1./2.)*L_capot) 0. ((-1.)* R_isol) ; p03 = ((-1./2.)*L_capot) R_isol 0. ; p04 = ((-1./2.)*L_capot) 0. R_isol ; p05 = ((-1./2.)*L_capot) 0. (Ri_tubm + Ep_tub + Ep_resi) ; p06 = ((-1./2.)*L_capot) (Ri_tubm + Ep_tub + Ep_resi) 0. ; l_ver1 = D 7 p01 p02 ; l_arc1 = CERCLE 5 p02 a00 p03 ; l_arc2 = CERCLE 5 p03 a00 p04 ; l_ver2 = D 7 p04 p05 ; l_arc3 = CERCLE 5 p05 a00 p06; l_arc4 = CERCLE 5 p06 a00 p01; s_capotb = SURF PLANE (l_ver1 ET l_arc1 ET l_arc2 ET l_ver2 ET l_arc3 ET l_arc4); ep_ref = (L_capot-L_isol)/ 2. ; capotb1 = VOLU TRAN s_capotb 4 (ep_ref 0. 0.); capotb11 = capotb1 SYME 'PLAN' x00 x01 x02 ; * partie haute du capot p01 = ((-1./2.)*L_capot) 0. ((-1.)* R_capot) ; p02 = ((-1./2.)*L_capot) 0. ((-1.)* R_isol) ; p03 = ((-1./2.)*L_capot) R_isol 0. ; p04 = ((-1./2.)*L_capot) 0. R_isol ; p05 = ((-1./2.)*L_capot) 0. R_capot ; p06 = ((-1./2.)*L_capot) R_capot 0. ; l_ver1 = D 4 p01 p02 ; l_arc1 = CERCLE 5 p02 a00 p03 ; l_arc2 = CERCLE 5 p03 a00 p04 ; l_ver2 = D 4 p04 p05 ; l_arc3 = CERCLE 5 p05 a00 p06; l_arc4 = CERCLE 5 p06 a00 p01; s_capot = SURF PLANE (l_ver1 ET l_arc1 ET l_arc2 ET l_ver2 ET l_arc3 ET l_arc4); capotg = s_capot VOLU 'GENE' lgen_cpm ; sf_cap = (l_arc3 ET l_arc4) GENE lgen_cpm; capotm = capotg ET capotb1 ET capotb11; ELIM 1.E-9 capotm ; capotm = capotm COUL BLEU; * * tout la surface extérieure du capot metallique est refroidie * on a déjà incorporé le refroidissement des brides, dans le refroidissement par air s_froid = s_capotb ET (s_capotb SYME 'PLAN' x00 x01 x02) ET sf_cap ; ELIM 1.E-9 s_froid capotm ; capotb1 = capotb1 COUL ROUGE ; capotg = capotg COUL BLEU ; *TRAC oeil cach (capotb1 ET capotg); * tout = poudre ET navette ET gaz ET tube ET tub_resi ET isolant ET capotm ; TITR 'Four DTEC/SDTC avec poudre contenue dans une navette'; TRAC CACH oeil FACE tout; *optio donn 5; * * ecriture des liaisons * ecriture de la liaison entre poudre et navette * attention il faut que ces surfaces soient bien attachées à leurs solides LIST (NBNO s_pou_d);LIST (NBNO s_nav3pg); lia_pm1 = LIAISON 1.E-9 s_pou_d s_nav3pg ; LIST (NBNO s_pou_g);LIST (NBNO s_nav3pd); lia_pm2 = LIAISON 1.E-9 s_pou_g s_nav3pd ; LIST (NBNO s_pou_h);LIST (NBNO s_fon_ph); lia_pm3 = LIAISON 1.E-9 s_pou_h s_fon_ph ; lia_pmet = lia_pm1 ET lia_pm2 ET lia_pm3 ; * * ecriture de la liaison entre tube et navette ELIM 1.E-9 navette tube; * * ecriture de la liaison entre le tube et le reste du four ELIM 1.E-9 (tube ET tub_resi ET isolant ET capotm); * * ecriture de la liaison entre poudre et gaz * on suppose une liaison parfaite entre les 2 ELIM 1.E-9 poudre gaz; * * ecriture de la liaison entre gaz et navette ELIM 1.E-9 gaz navette; * * ecriture de la liaison entre gaz et tube ELIM 1.E-9 gaz tube; * * definition des surfaces rayonnantes internes s_met_r = s_bt_1 ET s_bt_2 ET s_tug ET s_nav2m ET s_nav5ig ET s_nav5id ET s_nav5og ET s_nav5od ET (s_nav1g ET s_nav3g ET s_nav4g ET s_nav5g) ET (s_nav1d ET s_nav3d ET s_nav4d ET s_nav5d); ELIM 1.E-9 s_met_r ; ELIM 1.E-9 s_met_r tout; *TRAC CACH oeil s_met_r ; * ELIM 1.E-9 s_pou_r poudre; ELIM 1.E-9 s_met_r s_pou_r; cavite = s_met_r ET s_pou_r; * * emplacement des thermocouples de régulation Th_Regul = TABLE ; Th_Regul . point = TABLE ; Th_Regul . point . 1 = tub_res1 POINT 'PROC' ((-1*L_resi/3) 0. (Ri_tubm + Ep_tub + (0.5 * Ep_resi))) ; Th_Regul . point . 2 = tub_res2 POINT 'PROC' (0. 0. (Ri_tubm + Ep_tub + (0.5 * Ep_resi))) ; Th_Regul . point . 3 = tub_res3 POINT 'PROC' ((L_resi/3) 0. (Ri_tubm + Ep_tub + (0.5 * Ep_resi))) ; * * emplacement des thermocouples extérieurs Th_Ext = TABLE ; Th_Ext . point = TABLE ; Th_Ext . point . 1 = tube POINT 'PROC' (-0.692 0. (Ri_tubm + Ep_tub)) ; Th_Ext . point . 2 = tube POINT 'PROC' (0.648 0. (Ri_tubm + Ep_tub)) ; Th_Ext . point . 3 = capotm POINT 'PROC' (0. 0. R_capot) ; * * Vérification surface de rayonnement *m_r_fer = MODE s_met_r RAYONNEMENT; *em_fer = MATE m_r_fer 'EMIS' V_em_fer; *m_r_poud = MODE s_pou_r RAYONNEMENT; *em_poud = MATE m_r_poud 'EMIS' V_em_pou; *mrt = m_r_fer ET m_r_poud; *TRACE CACH (s_met_r ET s_pou_r); *em_tout = em_fer ET em_poud; *fft = FFOR mrt 'SYME' (0. 0. 0.) (1. 0. 1.) (-1. 0. 1.) ; *chamr = RAYE mrt fft em_tout ; *nbpsup = NBNO cavite; * LIST (NBNO tout); LIST (NBEL tout); * MESURE poudre ; TRAC CACH navette; SAUV ; FIN ; *