УДК 62.762 DOI: 10.12737/23161
А.В. Иванов
ВЛИЯНИЕ ФОРМЫ ЗАЗОРА БЕСКОНТАКТНОГО УПЛОТНЕНИЯ ВЫСОКООБОРОТНОЙ ТУРБОМАШИНЫ НА ЕГО ПАРАМЕТРЫ
Проанализировано влияние формы уплотни-тельного зазора на подъемную силу в нем, что особенно важно при определении влияния уплотнений на критические частоты вращения ротора. Приведены эпюры давления жидкости в зазоре уплотне-
ния при отклонениях формы и перекосе осей роторного и статорного элементов уплотнения.
Ключевые слова: бесконтактное уплотнение, форма зазора, подъемная сила, перекос, конусность, эллипсность.
Л.У. Туаиоу
INFLUENCE OF GAP FORM IN NON-CONTACT SEAL OF HIGH-SPEED TURBOMACHINE UPON ITS PARAMETERS
Consolidation is one of the most significant elements in any turbomachine defining both economy and capacity for work of a unit. In high-speed pumps and turbines which belong turbo-pump units (TPU) of liquid-propellant engines (LPE) to, one uses mainly non-contact seals and in the first place groove seals. A radial seal gap its value and a form is a basic factor defining parameters and characteristics of sealing. In the paper there is shown the influence of geometrical characteristics of a non-contact seal gap in a highspeed turbomachine upon its characteristics. The influence of the form of a sealing gap upon carrying capacity in it is analyzed that is particularly significant at the definition of seal influence upon critical frequencies of rotor rotation. The mutual misalignments in rotor and stator elements of a seal result in rotor stability reduc-
Уплотнение - один из наиболее ответственных элементов любой турбома-шины, определяющих как экономичность, так и работоспособность агрегата. В высокооборотных насосах и турбинах, к которым относятся турбонасосные агрегаты (ТНА) жидкостных ракетных двигателей (ЖРД), применяют в основном бесконтактные уплотнения, в первую очередь щелевые. Радиальный уплотнительный зазор, его величина и форма - основной фактор, определяющий параметры и характеристики уплотнения. Факторы, влияющие на зазор в уплотнении, приведены на рис. 1.
Стремление к снижению массы конструкции, что особенно важно для двигателей летательных аппаратов, часто приводит к снижению жесткости деталей и
tion, the occurrence of additional radial loads upon bearings. Deviations in a surface form and deformations of sealing elements result in the pressure redistribution in a seal gap that, in its turn, changes forces character affecting a rotor from the side of a seal and, hence, its dynamic characteristics. The analysis of the influence upon carrying capacity in a seal gap of such factors as ellipse, obliquity, barreling, saddling is carried out. The results shown in the paper allow defining the tolerances influence upon manufacturing and deformations of sealing elements upon dynamic characteristics of a high-speed turbomachine rotor at the stage of designing a "rotor-bearings-seals" system.
Key words: noncontact seal, gap form, carrying capacity, misalignment, obliquity, ellipse.
сборочных единиц, что способствует повышению их деформаций в рабочих условиях. Для насосов и турбин это приводит к ухудшению условий работы уплотнений и опор роторов [1].
Статорные элементы уплотнений роторов часто либо входят в состав корпусов, либо жестко с ними завязаны. В процессе работы агрегата они подвергаются существенному силовому и температурному воздействию. В процессе цикла «запуск - работа на номинальном (стационарном) режиме - останов» чередуются динамические и стационарные воздействия. При стационарном режиме материал детали находится в условиях упругости, прочность обусловлена действующими напряжениями и критериями кратковременного разрушения.
Рис. 1. Факторы, влияющие на радиальный зазор в уплотнении
На переходных режимах, включающих запуск, останов, переход с одного стационарного режима на другой, температурные поля неравномерны, разности температур и температурные напряжения, как правило, достигают максимальных значений. Наличие концентраторов напряжений может привести к появлению пластических деформаций при запусках и остановах. При стационарном режиме при воздействии повышенных температур может наблюдаться ползучесть [2]. Таким образом, прочность статорных элементов уплотнения, как замкнутых на корпусные элементы (щелевые и лабиринтные уплотнения), так и не замкнутых на них (уплотнения с плавающими и полуподвижными кольцами), должна оцениваться с учетом характеристик термоусталости и ползучести. Нестационарные (неравномерные) температурные поля приводят к короблению элементов статора, в первую очередь корпусов. Вследствие разности температур частей корпуса происходит искривление оси корпуса и его выпучивание, что приводит к изменению зазоров в уплотнениях, возможному задеванию и износу элементов уплотнения. При значительных задеваниях (касаниях) в уплотнениях ситуация может усугубляться возможным прогрессирующим тепловым искривлением ротора с весьма серьезными последствиями, особенно для кислородных насосов.
При расчете уплотнения основной задачей является определение расхода рабочего тела через него. В то же время анализ динамики роторных систем требует определения сил, действующих на ротор. В системе «ротор - уплотнения - опоры» существенные трудности связаны с расчетом гидромеханических сил в уплотни-тельном зазоре. Для уплотнений проточной части насосов и турбин ТНА ЖРД характерны перепады давления до 35 МПа и окружные скорости роторного элемента уплотнения до 600 м/с, при этом, в зависимости от его конструкции, уплотнение может вести себя как дополнительная опора жидкостного трения.
Моделирование течения рабочего тела в уплотнительном зазоре должно учитывать две особенности: 1) малая (сотни или десятки микрон) величина радиального уплотнительного зазора; 2) достаточно сложная геометрия радиального зазора. Нахождение распределения давления в зазоре между роторным и статорным элементами уплотнения основано на совместном решении уравнений гидромеханики, отражающих три фундаментальных закона сохранения: импульса, массы и энергии. Основными допущениями, принятыми в расчете поля давления, являются следующие:
- рабочее тело уплотнения - сплошная среда, заполняющая весь зазор;
- рабочее тело уплотнения считаем ньютоновской жидкостью;
- среду считаем изотропной;
- вследствие малой величины уп-лотнительного зазора по сравнению с длиной и радиусом уплотнения не учитываем его кривизну;
- изменение термодинамических параметров по толщине смазочного уплот-нительного зазора считаем несущественным;
- полагаем значительными градиенты скоростей только в направлении нормали к уплотнительным поверхностям, а скорость рабочего тела в этом направлении считаем малой;
- полагаем, что движение ротора относительно статора в осевом направлении не влияет на режим течения;
- пренебрегаем действием сил поверхностного натяжения, инерции и тяжести рабочего тела, а также изменением объема рабочего тела, связанным с изменением температуры;
- скольжение рабочего тела относительно роторной и статорной поверхностей уплотнения и ротора отсутствует (явление адсорбции), а скорость граничных слоев
dV - / 2
Р
рабочего тела равна скорости прилегающих поверхностей.
В основу изучения движения рабочего тела в уплотнительном зазоре положены два фундаментальных уравнения гидромеханики, отражающие законы сохранения импульса и массы [3; 4], - уравнение На-вье-Стокса и уравнение неразрывности:
dt
pF-grad| p + 3\divVj + 2div(\defV), dp + div(pV)
= 0,
где V = (V.,:, vz) - вектор скорости рабочего тела; F - плотность распределения массовых сил; def V - деформация поля
вектора; V - тензор скоростей деформации.
Уравнения справедливы как для ламинарного, так и для турбулентного режима течения. С учетом сделанных допуще-
ний о пренебрежении влиянием массовых и инерционных сил, за исключением сил турбулентной вязкости, а также помня о малости градиентов скоростей по осям л и г и незначительности скорости Vу вследствие малой толщины слоя, уравнение На-вье-Стокса в проекциях на оси координат запишем в следующем виде:
dp dx
д
í дк >
|1-
дУ | дУ
d(pvXvy)
дУ
dp д
д (Pv'y ) . Эр dy ' dz
д
í dv >
|
дУ | дУ
d(pvZvy)
дУ
Постоянство давления по высоте уп-лотнительного зазора позволяет непосредственно проинтегрировать первое и третье уравнения этой системы по координате у. Следуя [5], введем коэффициенты турбулентности Кл и Кг, учитывающие влияние
др =_э_Г э^.
дл ду ^ л ду У В итоге получим основное уравнение для определения поля давления - уравнение Рейнольдса, обобщенное для случая
дополнительной турбулентной вязкости. Коэффициенты турбулентности линейно, но с достаточной точностью аппроксимируют зависимость турбулентных напряжений от вязкого напряжения:
ЭР = 0. ЭР = _Э_
Эу ' ~ ~
(
ITT Л
z dy.
д_ dx
h3p dp
d + —
dz
h3p dp \iKz dz
dz dy
турбулентного течения вязкого сжимаемого рабочего тела:
■ 6—(pUh)- 12р V + 12h dP. dxVF ; F dt
движения ротора в декартовой системе координат:
\Kx dX
где значения скоростей U и V в точках на поверхности ротора находятся по следующим зависимостям при рассмотрении
U =wR + X cos a-Y sin а; V = X sin а + Y cos а. Уравнение баланса энергий для однофазного течения рабочего тела:
íuh h ^
ph
dI dp dT
--- + CP —
dp dt dt
+ p
dp
2 12|Kx dx
dI dp ^ dT
--- + CP —
dp dx dx
ph3 dp
12\Kz dz
dI dp ^ dT
--- + CP —
dp dz dz
dp , Uh dp U2 ^
= h +--- +-\Kx.
dt 2 dx h
1
0,9 0,8 0,7 0,6 0,5 | 0,4 . 0,3 ' 0,2 0,1 0
0 0,2 0,4 0,6 0,8 1
Относительная частота вращения (п/птах)
-Т-20 К —Т-22,5 К ---Т-25 К
а)
1
0,9
ГС
0,8
и
ОС 0,7
ас ? 0,6
и а — 0,5
о с 1 0,4 и.
ОС ГС и^ 0,3
£ 0,2
0,1
и 0
X
О
Ар-10 МПа -5R-0,035 мм ю-5000 рад/с
//
///'
/¿V
/ /'
Й--------
0,2 0,4 0,6
Относительный эксцентриситет
Т-20 К —Т-22,5 К ---Т-25 К
0,8
б)
Рис. 2. Зависимость подъемной силы от радиального зазора (а), эксцентриситета (б) и температуры рабочего тела
1 0,9
(С
| 0,8 и
к 0,7
(С
! 0,6
¡5 _ 0,5 ^ 1)
° Е 0,4
с и.
* ¡Е 0,3
X
£ 0,2 си
£ 0,1 и
5 0
1
0,9
(С
1 0,8 и
к 0,7
(С
1 0,6 и
— 0,5 о = с и.Е 0,4
г ^
? ~ 0,3
Ар-10 МПа 1 / X1
/ У
----
,4
_ _ _ _ ____ ____ ____
_ 5
6^
0 0,2 0,4 0,6 0,8 1 Перекос (0/0тах) а)
0,2 0,1 0
Ар-10 МПа Т-20 К е-0,3 D-100 мм 8R-0,05 мм
/ ✓ / У
е =emax=2SR/L _ У
/
е=о
_ — — "
0,4 0,6
Частота вращения (п/пт
х)
б)
Рис. 3. Влияние перекоса (а) и скорости вращения (б) на подъемную силу: 1 - ш=9000 рад/с; 2 - ш=7000 рад/с; 3 - ш=5000 рад/с; 4 - ш=3000 рад/с; 5 - ш=1000 рад/с; 6 - ш=0 рад/с
Система уравнений - Рейнольдса и баланса энергий - является неопределенной, так как число неизвестных параметров превышает число уравнений. Дополнительными соотношениями, доопределяющими эту систему, являются зависимости теплофизических свойств рабочего тела от давления и температуры, а также уравнение состояния, связывающее плотность, давление и температуру: д = д(р, Т); I = I (р, Т); / (р, р, Т) = 0.
Анализ влияния формы зазора на подъемную силу в уплотнении проводился для уплотнений водородных насосов.
Рис. 2а показывает зависимость подъемной силы от величины радиального зазора 5Л при различных частотах вращения ю и температуре рабочего тела Т. Видно, что увеличение радиального зазора сопровождается снижением подъемной силы. Отрицательное влияние на подъемную силу оказывает и увеличение температуры рабочего тела, сопровождающееся снижением его вязкости.
Увеличение эксцентриситета е приводит к увеличению подъемной силы (рис. 2б). Здесь же видно, что в области малых эксцентриситетов (~ё < 0,5) характеристика уплотнения «подъемная сила - перемещение» имеет почти линейный характер, что позволяет использовать для анализа задач динамики ротора понятие динамических коэффициентов рабочего тела в уплотни-тельном зазоре, линейно аппроксимирующих его реакции. Увеличение температуры рабочего тела сопровождается снижением подъемной силы из-за роста концентрации газовой фазы, падения вязкости и плотности рабочего тела.
На рис. 3а представлено влияние перекоса осей роторного и статорного элементов уплотнения. При малых частотах вращения или неподвижном вале подъемная сила формируется за счет гидростатического эффекта. Перекос приводит к ее снижению вследствие уменьшения гидравлического сопротивления в направлении увеличения зазора и, следовательно, уменьшения градиента давления в нагруженной зоне. С ростом частоты вращения ситуация меняется, при 0> 0,7 можно наблюдать существенный рост подъемной
0
0
0,2
0,8
силы. Преобладание гидродинамического эффекта приводит к большему росту давления в нагруженной зоне вследствие нелинейности характеристики уплотнения. Отметим некоторое уменьшение подъемной силы при малых скоростях по сравнению с неподвижным ротором (рис. 3б).
Влияние отклонений формы роторного и статорного элементов уплотнения на характеристики было исследовано на примере уплотнений с параметрами: О = 0,1 м; I = 0,009 м; 8Я = 0,05 мм; Т = 20 К; рабочее тело - жидкий водород. Результаты расчетов представлены на рис. 4-6.
Уменьшение подъемной силы наблюдается при величине эксцентриситета
0,3 (величина эксцентриситета взята в отношении к радиальному зазору в уплотнении) и выше. Аналогичная картина наблюдается и при изменении скорости вращения ротора, где видимое изменение подъемной силы происходит при относительной скорости 0,3 и выше. Также отмечено, что при изменении давления на входе наблюдается быстрый рост подъемной силы в уплотнении с конусными уплотнитель-ными поверхностями. При относительном давлении больше 0,2 рост подъемной силы замедляется, однако остается в положительной тенденции.
0,2 0,3 0,4
Эксцентриситет
0,6 £
0,4 0,6 0,8
Частота вращения (п/п^,)
0,4 0,6 0,8
Давление на входе (р/рта>)
а) б) в)
Рис. 4. Влияние конусности на подъемную силу в зависимости от эксцентриситета (а), частоты вращения (б), перепада давления (в)
Эллипсность, в отличие от конусности, не оказывает значительного влияния на подъемную силу ни при изменении скорости вращения, ни при увеличении давления на входе в уплотнение. При увеличении скорости вращения ротора подъемная сила уплотнения с эллипсными втулками уменьшилась по сравнению с уплотнением с идеальными втулками. Изменение подъемной силы при увеличении скорости вращения составило не больше 6%.
Аналогичная картина наблюдается при изменении эксцентриситета, где изменение подъемной силы не превышает 5%. Данное обстоятельство объясняется тем, что эл-липсность изменяет геометрию радиального зазора в поперечном направлении, в то время как конусность - в продольном, тем самым уменьшая опорную площадь контакта (т.е. площадь, где создается эффективный гидродинамический клин и достаточный перепад давления).
»_
Др=10 МПа Т=20 К D=100 мм
0/а=1,001 . //
Ш=1,0015
0,2 0,3 0,4 0,5 0,6
Эксцентриситет
0,4 0,6 0,8
Частота вращения (п/птах)
0,4 0,6 0,8 1
Давление на входе (р/р^)
а) б) в)
Рис. 5. Влияние эллипсности на подъемную силу в зависимости от эксцентриситета (а), частоты вращения (б), перепада давления (в)
0,1
0,5
0,2
0,2
0,2
0,1
0,2
При проведении численного эксперимента было выявлено, что влияние, оказываемое на подъемную силу седлообраз-ностью поверхностей уплотнения, идентично влиянию бочкообразности (при сохранении исходной формы одного из эле-
ментов уплотнения). Данное обстоятельство объясняется тем, радиальные зазоры, сформированные этими отклонениями формы, одинаковы. Поэтому в дальнейшем результаты по седлообразности представлены не будут.
1 0,9 § 0,8 к 0,7
Л 0,6
£ _ 0,5
I 3 0,4
л £ 0,3
I °'2
§ о,1 г 0
1 0,9
I 0,
к 0
I 0
Л _ 0 § 10 = ^ 0 ;_ 0 г 0
Т-20 К - Р-100 мм — о/а=1,ооо5 — — -1 к----
1\
-У
/
0,4 0,6 0,8
Частота вращения (п/пта.)
0,4 0,6 0,8
Давление на входе (р/рта.)
а) б) в)
Рис. 6. Влияние бочкообразности на подъемную силу в зависимости от эксцентриситета (а), частоты вращения (б), перепада давления (в)
0,2
0,2
0,1
0,2 0,3
Эксцентриситет
0,4
0,5
0,6
Подъемная сила в уплотнении при бочкообразности его поверхностей имеет существенное изменение лишь при изменении скорости вращения, в то время как при увеличении давления на входе и эксцентриситета это изменение не более 5...6%.
Перекосы плавающих колец относительно уплотнительных поверхностей крыльчаток приводят к снижению устойчивости уплотнений и появлению дополнительных радиальных нагрузок на подшипники. На рис. 7 показаны эпюры давления жидкости в щели уплотнения при отклонениях формы и перекосе осей роторного и статорного элементов уплотнения.
Отклонения формы поверхностей и деформации элементов уплотнения приводят к перераспределению давления в уп-лотнительном зазоре, что, в свою очередь, меняет характер сил, действующих на ротор со стороны уплотнения, и, следовательно, на его динамические характери-
стики. Например, участок щели, показанный на нижней части рисунка (рис. 7е), при перекосе осей становится конфузор-ным, а в верхней части - диффузорным. Как видно из приведенных эпюр, суммарное давление на уплотнительную поверхность конфузорного участка плавающего кольца значительно больше, чем на диф-фузорный участок. Поэтому плавающее кольцо сместится относительно крыльчатки, увеличивая размер щели на конфузор-ном участке и уменьшая на диффузорном. При этом утечка через конфузорный участок будет существенно превышать утечку через диффузорный участок, что приведет к изменению давления в полости крыльчатки (давление в зоне диффузорного участка будет больше давления в зоне конфу-зорного участка), то есть к появлению дополнительной радиальной нагрузки. В то же время перепад давления на поясок крыльчатки действует в противоположную сторону.
\
>4
1 1
Ф) Р,
Р„ Ф) г)
д) е) ж)
Рис. 7. Характерные эпюры давлений вдоль оси уплотнения: а - при отсутствии деформаций; б - при конфузорной форме щели; в - при диффузорной форме щели; г - при бочкообразности роторного элемента уплотнения; д - при седлообразности роторного элемента уплотнения; е - при перекосе осей; ж - при изгибе оси ротора
Можно сделать следующие выводы.
Отклонения формы поверхностей и деформации элементов уплотнения приводят к перераспределению давления в уплотни-тельном зазоре, что, в свою очередь, меняет характер сил, действующих на ротор со стороны уплотнения, и, следовательно, на его динамические характеристики. Результаты, приведенные в статье, позволяют на
стадии проектирования системы «ротор -опоры - уплотнения» определить влияние допусков на изготовление и деформаций элементов уплотнения на динамические характеристики ротора высокооборотной турбомашины, спрогнозировать поведение ротора при различных сочетаниях внешних и внутренних факторов.
СПИСОК ЛИТЕРАТУРЫ
1. Силовые деформации корпуса подшипника при нестационарных тепловых потоках / Ю.А.Равикович [и др.] // Вестник МАИ. - 2009. - Т. 16. - № 3. - С. 65-68.
2. Костюк, А.Г. Динамика и прочность турбома-шин: учеб. для вузов / А.Г .Костюк. - 3-е изд., перераб. и доп. - М.: Изд. дом МЭИ, 2007. - 476 с.
3. Камерон, А. Теория смазки в инженерном деле /
A.Камерон. - М.: Машгиз, 1962. - 296 с.
4. Лойцянский, Л.Г. Механика жидкости и газа / Л.Г.Лойцянский. - М.: Дрофа, 2003. - 840 с.
5. Типей, Н. Подшипники скольжения: расчет, проектирование, смазка / Н.Типей,
B.Н.Константинеску [и др.]. - Бухарест: АН РНР, 1964. - 458 с.
1. Power deformations of bearing body at non-stationary thermal flows / YuA.Ravikovich [et al.] // Bulletin of MAI. - 2009. - Vol. 16. - № 3. - pp. 65-68.
2. Kostyuk, А.G. Turbomachines Dynamics and Strength: textbook for colleges / А.G.Kostyuk. - 3d Edition revised and supplemented. - М.: MEI Publishers, 2007. - pp. 476.
Сведения об авторах:
Иванов Андрей Владимирович, д.т.н., профессор кафедры «Ракетные двигатели» Воронежского го-
3. Cameron, А. Theory of Lubrication in Engineering / А.Cameron. - М.: Machgiz, 1962. - pp. 296.
4. Loitsyansky, L.G. Liquid and Gas Mechanics / L.G.Loitsyansky. - М.: Drofa, 2003. - pp. 840.
5. Tipey, N. Sliding Bearings: Computation, Design, Lubrication / N.Tipey, V.N.Konstantinesku [et al.]. - Bucharest: AS RPR, 1964. - pp. 458.
Статья поступила в редколлегию 18.10.2016. Рецензент: д.т.н., профессор Орловского государственного университета Чернышев В.И.
сударственного технического университета, тел.: (473)-252-27-52, е-mail: iav308@inbox.ru.
Ivanov Andrey Vladimirovich, D. Eng., Assistant Prof., Prof. of the Dep. "Missile Engine", Voronezh
State Technical University, E-mail: iav308@inbox.ru.