Диссертация (Щёточные уплотнения в роторных системах авиационных двигателей), страница 77
Описание файла
Файл "Диссертация" внутри архива находится в папке "Щёточные уплотнения в роторных системах авиационных двигателей". PDF-файл из архива "Щёточные уплотнения в роторных системах авиационных двигателей", который расположен в категории "". Всё это находится в предмете "технические науки" из Аспирантура и докторантура, которые можно найти в файловом архиве МАИ. Не смотря на прямую связь этого архива с МАИ, его также можно найти и в других разделах. , а ещё этот архив представляет собой докторскую диссертацию, поэтому ещё представлен в разделе всех диссертаций на соискание учёной степени доктора технических наук.
Просмотр PDF-файла онлайн
Текст 77 страницы из PDF
9 5 } ] [ expr {$Rr+3 . 7 0 } ] $Rslappend RadiusN [ llength $Radius ]makeauxcurve $Radius $strAux $ s t r C r $N $centreAx# AUX22set N [ expr {$N+1}]set centreAx0lappend centreAx 0 0475Файлы для создания геометрии уплотнительного узла SSB-1 в ICEM CFD430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489set Radius$Rrlappend Radius $Rs $Rdlappend RadiusN [ llength $Radius ]makeauxcurve $Radius $strAux $ s t r C r $N $centreAx# AUX23set N [ expr {$N+1}]set centreAx−$Ldlappend centreAx −$Ldset Radius$Rrlappend Radius $Rdlappend RadiusN [ llength $Radius ]makeauxcurve $Radius $strAux $ s t r C r $N $centreAx###################################################### SHIFTING ROTOR TO ECCENTRIC POSITION#####################################################i f { $ECCENTRICITY } { \for { set i 1} { $ i <$ r o t o r P t +1} { incr i } { \set pt $ s t r R o t o r P t ; \ic_move_geometry p o i n t names [ append pt $ i ] t r a n s l a t e $ v e c t o r E c c \}; \for { set i 1} { $ i <$ r o t o r P t } { incr i } { \set c r $ s t r R o t o r C r ; \set s f $ s t r R o t o r S f ; \ic_move_geometry c u r v e names [ append c r $ i ] t r a n s l a t e $ v e c t o r E c c ; \ic_move_geometry s u r f a c e names [ append s f $ i ] t r a n s l a t e $ v e c t o r E c c \}; \for { set i 1} { $ i <=$N} { incr i } { \i f { $ i <10} { \set c r "AUX0" ; \} else { \set c r "AUX" ; \}ic_move_geometry c u r v e names [ append c r $ i "CR1" ] t r a n s l a t e $ v e c t o r E c c \}; \set c r "AUX07" ; \ic_move_geometry c u r v e names [ append c r "CR2" ] t r a n s l a t e $ v e c t o r E c c ; \set c r "AUX12" ; \ic_move_geometry c u r v e names [ append c r "CR2" ] t r a n s l a t e $ v e c t o r E c c ; \}###################################################### CREATION OF INLET AND OUTLET SURFACES###################################################### Inlet surfaceset INLETSF [ ic_geo_create_surface_from_curves INLET INLETSF1e−6 \{AUX01CR1 AUX01CR2} ]# Outlet surfaceset OUTLETSF [ ic_geo_create_surface_from_curves OUTLET OUTLETSF 1e−6 \[ lappend a r e s "AUX${N}CR1" "AUX${N}CR2" ] ]# B r i s t l e pack s u r f a c e sset BRUSHSF1 [ ic_geo_create_surface_from_curves BRUSH1 BRUSHSF1 1e−6 \{AUX16CR1 AUX16CR2} ]set BRUSHSF2 [ ic_geo_create_surface_from_curves BRUSH2 BRUSHSF2 1e−6 \{AUX17CR1 AUX17CR2} ]###################################################### CREATION OF MATERIAL POINTS476Файлы для создания геометрии уплотнительного узла SSB-1 в ICEM CFD490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521#####################################################ic_geo_create_volume $mPtSEAL_UP_B SEAL_UP_B SEAL_UP_Bic_geo_set_part body SEAL_UP_B SEAL_UP_B 0ic_geo_create_volume $mPtSEAL_DOWN_B SEAL_DOWN_B SEAL_DOWN_Bic_geo_set_part body SEAL_DOWN_B SEAL_DOWN_B 0ic_geo_create_volume $mPtBRUSH_B BRUSH_B BRUSH_Bic_geo_set_part body BRUSH_B BRUSH_B 0###################################################### CREATION OF POINTS ON AUXILIARY CURVES# FOR ASSOCIATION WITH BLOCK MODEL VERTICES#####################################################set s t r P t "PT"fo r { set k 1} {$k<=$N} { incr k} { \i f {$k <10} { \set strAux "AUX0" ; \} else { \set strAux "AUX" ; \}for { set i 1} { $ i <=[ lindex $RadiusN [ expr {$k−1 } ] ] } { incr i } { \fo r { set j 1} { $j <=4} { incr j } { \set c r $strAux ; set c r 1 $ s t r C r ; set pt1 $ s t r P t ; \set s t r 1 $strAux ; set s t r 2 $strAux ; \set s t r 3 [ append s t r 1 $k $ c r 1 $ i ] ; \lappend s t r 3 [ expr {0 .
1 2 5 * ( 2 *$j−1 ) } ] ; \i c _ p o i n t crv_par [ append c r $k ] [ append s t r 2 $k $pt1 $ i $ j ] \$str3 \} \} \}477Приложение ГПакет POLLS для автоматизированногорасчёта щёточного уплотнения в OpenFOAMГ.1. Функция генерации параметров моделипористой средыЛистинг Г.1. pollsPorous.m1function p o l l s P o r o u s ( bb , d , N, phi , f i l e _ o u t )234567%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% EXPRESSIONS%p h i = p h i * pi / 1 8 0 ;89e p s i l o n = 1 − ( ( pi * d^2 * N) / ( 4 * bb * cos ( p h i ) ) ) ;101112C = ( 1 / e p s i l o n ^3) * ( ( 1 − e p s i l o n ) / d ) ^ 2 ;D = ( 1 / e p s i l o n ^3) * ( ( 1 − e p s i l o n ) / d ) ;131415tau_Wu = 1/2 * ( 1 + 0 . 5 * ( 1 − e p s i l o n ) ^ 0 .
5 + ( 1 − e p s i l o n ) ^ 0 . 5 * . . .( ( 1 / ( 1 − e p s i l o n ) ^ 0 . 5 − 1)^2 + 0 . 2 5 ) ^ 0 . 5 / ( 1 − ( 1 − e p s i l o n ) ^ 0 . 5 ) ) ;1617beta_Wu = 1 / ( 1 − ( 1 − e p s i l o n ) ^ 0 . 5 ) ;181920212223242526a_n = 72 * tau_Wu * C;a_s = a_n / 6 0 ;b_n = 3 * tau_Wu / 4 * ( 3 / 2 + 1/beta_Wu^4 − 5 / ( 2 * beta_Wu^ 2 ) ) * D;b_s = 0 . 0 ;%% END: EXPRESSIONS%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%272829303132%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% WRITE%f i d = fopen ( f i l e _ o u t , ’w ’ ) ;33343536373839404142fprintf ( fidfprintf ( fidfprintf ( fidfprintf ( fidfprintf ( fidfprintf ( fidfprintf ( fidfprintf ( fidfprintf ( fid,,,,,,,,,’ FoamFile \n ’ ) ;’ {\n ’ ) ;’version’format’class’location’object’ }\n ’ ) ;’ \n ’ ) ;2.0;\n ’ );a s c i i ;\n ’ );dictionary ;\n ’ );" constant ";\n ’ ) ;porosityProperties ;\n ’ );43478Пакет POLLS для автоматизированного расчёта щёточного уплотнения в OpenFOAM4445464748495051525354555657585960616263646566676869f p r i n t f ( f i d , ’ p o r o s i t y 1 \n ’ ) ;f p r i n t f ( f i d , ’ {\n ’ ) ;fprintf ( fid , ’typeDarcyForchheimer ; \ n ’ ) ;fprintf ( fid , ’activeyes ; \ n ’ ) ;fprintf ( fid , ’cellZoneB_BRUSH; \ n ’ ) ;f p r i n t f ( f i d , ’ \n ’ ) ;fprintf ( fid , ’D a r c y F o r c h h e i m e r C o e f f s \n ’ ) ;fprintf ( fid , ’{\n ’ ) ;fprintf ( fid , ’dd [ 0 −2 0 0 0 0 0 ] (%10.6 e %10.6 e %10.6 e ) ; \ n ’ , .
. .a_n , a_n , a_s ) ;fprintf ( fid , ’ff [ 0 −1 0 0 0 0 0 ] (%10.6 e %10.6 e %10.6 e ) ; \ n ’ , . . .b_n , b_n , b_s ) ;f p r i n t f ( f i d , ’ \n ’ ) ;fprintf ( fid , ’c o o r d i n a t e S y s t e m \n ’ ) ;fprintf ( fid , ’{\n ’ ) ;fprintf ( fid , ’typecartesian ;\n ’ );fprintf ( fid , ’o r i g i n (0 0 0 ) ; \ n ’ ) ;fprintf ( fid , ’c o o r d i n a t e R o t a t i o n \n ’ ) ;fprintf ( fid , ’{\n ’ ) ;fprintf ( fid , ’typeaxesRotation ;\ n ’ ) ;fprintf ( fid , ’e1 ( 0 %10.6 f %10.6 f ) ; \ n ’ , sin ( p h i ) , cos ( p h i ) ) ;fprintf ( fid , ’e2 ( 1 0 0 ) ; \ n ’ ) ;fprintf ( fid , ’}\n ’ ) ;fprintf ( fid , ’}\n ’ ) ;fprintf ( fid , ’}\n ’ ) ;f p r i n t f ( f i d , ’ }\n ’ ) ;707172737475fclose ( f i d ) ;%% END: WRITE%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%7677end479Листинг Г.2.
pollsGeom.m1function pollsGeom ( )2345678%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% INPUT%f i l e _ o u t = ’ blockMeshDict ’ ;% Output f i l e w i t h geometryf i l e _ o u t 2 = ’ p o r o s i t y P r o p e r t i e s ’ ; % Output f i l e w i t h p o r ou s medium model9101112Drotor = 3 9 1 . 4 4 7 e −3;Hfence = 1 .
5 e −3;Lb = 9 . 2 5 e −3;% Rotor d i a m e t e r (m)% Backing p l a t e r a d i a l c l e a r a n c e (m)% B r i s t l e l e n g t h (m)Bb = 1 . 5 9 1 0 e −3;Hb = 0 . 1 4 2 2 e −3;% A c t u a l b r i s t l e pack t h i c k n e s s (m)% A c t u a l b r i s t l e pack r a d i a l c l e a r a n c e (m)H s t a t o r = 3 e −3;Lt = 5 3 . 2 e −3;Bbr = 3 e −3;Bbp = 1 . 5 2 e −3;Bfp = 1 . 5 2 e −3;%%%%%angle = 1 ;% 3D s e a l segment w i d t h ( d e g r e e )131448015161718192021S e a l c a v i t y h e i g h t (m)S e a l l e n g t h (m)Brush s e a l t h i c k n e s s (m)Backing p l a t e t h i c k n e s s (m)Front p l a t e t h i c k n e s s (m)22232425262728293031d = 70 e −6;% Wire d i a m e t e r (m)N = 240 e3 ;% Packing d e n s i t y ( b r i s t l e s /m c i r c .
)phi = 45;% Lay a n g l e ( d e g r e e )%% END: INPUT%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%3233343536%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% POROSITY PROPERTIES%Пакет POLLS для автоматизированного расчёта щёточного уплотнения в OpenFOAMГ.2. Скрипт генерации геометрии для blockMesh38394041p o l l s P o r o u s (Bb , d , N, phi , f i l e _ o u t 2 ) ;%% END: POROSITY PROPERTIES%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%42434445464748495051525354554815657585960616263%%%%%%%%%%%%%%%%%%%%%37−−−−31−−−−25−−−−19−−−−13|||||| 24 | 19 | 14 | 9 ||||||36−−−−30−−−−24−−−−18−−−−12|||||| 23 | 18 | 13 | 8 ||||||45−−−−−−−−−−4135−−−−29−−−−23−−−−17−−−−117−−−−−−−−−−−−−−−−3||||||||||29|| 22 | 17 | 12 | 7 ||2||||||||||44−−−−−−−−−−40−−−−−34−−−−28−−−−22−−−−16−−−−10−−−−−−−−6−−−−−−−−−−−−−−−−2||||||||||28| 26 | 21 | 16 | 11 | 6 |4|1||||||||||43−−−−−−−−−−39−−−−−33−−−−27−−−−21−−−−15−−−−9−−−−−−−−−5−−−−−−−−−−−−−−−−1||||||||||27| 25 | 20 | 15 | 10 | 5 |3|0||||||||||42−−−−−−−−−−38−−−−−32−−−−26−−−−20−−−−14−−−−8−−−−−−−−−4−−−−−−−−−−−−−−−−0646566676869%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% PARAMETERS%Npt = 4 6 ; pt = zeros ( Npt , 3 ) ; e d g e s = zeros ( Npt , 3 ) ;70717273747576D s t a t o r = Drotor + H s t a t o r * 2 ; Dfence = Drotor + Hfence * 2 ;%% END: PARAMETERS%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%Пакет POLLS для автоматизированного расчёта щёточного уплотнения в OpenFOAM3778798081%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% X SECTIONS%Nxi = 9 ; x i = zeros ( Nxi , 1 ) ;82838485868788x i ( 1 ) = Lt ; x i ( 2 ) = Lt − 30 e−3 + Bfp ; x i ( 3 ) = Lt − 30 e −3; x i ( 4 ) = x i ( 3 ) − ( Bbr − Bb) / 2 ; x i ( 5 ) = x i ( 3 ) − Bbr + Bb ;x i ( 6 ) = x i ( 3 ) − Bbr + Bb / 2 ; x i ( 7 ) = x i ( 3 ) − Bbr ; x i ( 8 ) = x i ( 7 ) − Bbp ; x i ( 9 ) = 0 ;%% END: X SECTIONS%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%899091929394%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Y SECTIONS%Nyi = 6 ; y i = zeros ( Nyi , 1 ) ;9548296979899100101102103104105106y i ( 1 ) = Drotor / 2 ;i f (Hb > 0 .