Full Text

Щеточные уплотнения относятся к классу контактных податливых уплотнений и предназначены для уменьшения утечки рабочего тела (жидкости или газа) через зазоры между вращающимися и неподвижными деталями турбомашин [1-3]. Известно об успешном внедрении этих уплотнений в газотурбинные двигатели за рубежом [4; 5], расчетно-экспериментальные исследования ведутся также в России [6-13]. Ключевой характеристикой щеточных уплотнений (рис. 1) является скорость износа щеточного пакета (щетины) при контакте с ротором, которая зависит от рабочих условий, жесткости щеточного пакета, материалов контактной пары и кинетики изменения рабочего натяга/зазора в уплотнении по режимам полетного цикла. Износ щетины ведет к постепенному увеличению утечки газа через уплотнение вплоть до образования радиального зазора на некоторых режимах работы. Именно износ определяет эффективность и ресурс щеточного уплотнения в составе газотурбинного двигателя. Линейная скорость износа щетины может быть оценена на основе закона Архарда: , где - экспериментально определяемый коэффициент пропорциональности; - линейная скорость поверхности ротора; - контактное давление между щеточным пакетом и поверхностью ротора [11]. Контактное давление определяется отношением силы реакции щеточного пакета , зависящей от величины текущего натяга D, к суммарной площади торцов щетинок: . Характерной особенностью зависимости силы реакции щеточного пакета от величины натяга / смещения является наличие гистерезиса в цикле «нагрузка - разгрузка» из-за трения в системе. Зависимость без учета воздействия давления газового потока может быть получена экспериментально при помощи простой лабораторной установки, показанной на рис. 1, в. Использование двух нагружающих колодок, имеющих различную длину по окружности, позволяет исключить краевой эффект и смоделировать ситуацию равномерного натяга между щеточным пакетом и ротором. Наличие осевого перепада давления газа на уплотнении (основной режим работы) приводит к заметному увеличению силы реакции щеточного пакета из-за взаимодействия щетинок с потоком газа. При установке уплотнения с зазором этот эффект приведет к поджатию щеточного пакета к поверхности ротора и частичному (или полному) закрытию зазора. Для испытаний в этом случае разработан стенд с вращающимся ротором и возможностью наддува и нагрева подаваемого воздуха. Cила реакции и жесткость щеточного пакета могут быть определены на основе соотношения для прогиба одной щетинки по модели плоского изгиба консольной прямолинейной балки [14]: , где (1) а б в Рис. 1. Щеточное уплотнение: конструкция, расчет, эксперимент; а - общий вид: 1 - щетинки, 2 - ротор, 3 - опорная пластина, 4 - передняя пластина; б - модель; в - нагружение колодкой для измерения статической жесткости щетины [12] Fig. 1. Brush seal: design, calculation, experiment; a - General view: 1 - bristles, 2 - rotor, 3 - support plate, 4 - front plate; б - model; в - loading with a block for measuring the static stiffness of the bristles [12] Здесь Е - модуль упругости материала щетинки (проволоки); - угол наклона щетинки; - диаметр щетинки; - длина щетинки, которая определяется наружным и внутренним диаметрами щеточного пакета. «Ужесточающее» влияние контактного взаимодействия щетинок в пакете и воздействие газового потока учтено эмпирическими коэффициентами и : , где (2) Здесь nщп - количество щетинок в единице окружной длины щетины. Значения коэффициентов взяты согласно данным работы [14]: ; при несовпадении осей ротора и статора (эксцентричное смещение ротора), при симметричном расширении ротора под действием тепловых и инерционных нагрузок. В настоящей работе выполнено уточнение жесткости щеточного пакета на основе решения задачи пространственного деформирования сегмента. Подобный подход использован в работах [15-21] и обоснован хорошим согласованием расчетных и экспериментальных результатов. Вид расчетных моделей приведен на рис. 2. Расчетное определение силы реакции щеточного пакета основано на решении контактной задачи взаимодействия щетинок между собой, с ротором и с задней (опорной) пластиной корпуса с учетом действия распределенной силы давления газового потока. Области контактного взаимодействия показаны на рис. 3: К1 - контакт между торцами щетинок и поверхностью ротора (коэффициент трения ); К2 - контакт между щеточным пакетом и опорной пластиной (коэффициент трения ); К3 - контакт между щетинками в щеточном пакете (коэффициент трения ). Каждая проволока выполнена в балочном представлении, что позволило создать вычислительно экономичную расчетную модель сегмента щеточного уплотнения. Использован стандартный 2-узловой конечный элемент (КЭ) балки Тимошенко с шестью степенями свободы в узле и кубическими функциями формы. В силу малости окружного размера сегмента кривизна щеточного пакета не учитывалась. Расположение щетинок в сегменте принято регулярным и допускало две предельные конфигурации, рис. 3. Осевое и окружное расстояния между центрами соседних проволок рассчитывались в зависимости от задаваемой характерной величины зазора между щетинками : где для линейного расположения щетинок в пакете; для шахматного расположения щетинок. Влияние деформированного состояния щеточного пакета на перераспределение параметров течения газа в нем на данном этапе работы не рассматривалось. Значения распределенных нагрузок в n-м узле определялись по распределению давления газа в щеточном пакете: ; где , - компоненты вектора нормали к контуру сечения щетинки в узле; , - перепады давления газа на сечении щетинки. Для вычисления вышеприведенных интегралов использовано параметрическое представление контура сечения щетинки и билинейное распределение давления между контрольными точками. Характерный вид нагрузок показан на рис. 5. а б в Рис. 2. Конечно-элементные модели сегмента щеточного пакета: 1 - щетинки, 2 - ротор, 3 - опорная пластина; а - модель на основе объемных КЭ (Adina) [12; 15-17]; б - модель на основе балочных КЭ (Abaqus) [18-20]; в - модель, созданная в настоящей работе (ANSYS), и фото изготовленного пакета Fig. 2. Finite element models of the brush package segment: 1 - bristles, 2 - rotor, 3 - support plate; а - a model based on volumetric FE (Adina) [12; 15-17]; б - a model based on beam FE (Abaqus) [18-20]; в - the model created in this work (ANSYS) and a photo of the manufactured package Рис. 3. Области контактного взаимодействия в щеточном уплотнении Fig. 3. Contact areas in the brush seal а б Рис. 4. Профили безразмерного давления P в радиальном и осевом направлениях: а - осевые сечения; б - радиальные сечения Fig. 4. Profiles of dimensionless pressure P in the radial and axial directions: a - axial sections; б - radial sections Рис. 5. Распределение газовой нагрузки по длине щетинки Fig. 5. Gas load distribution along the bristle length Таблица 1 / Table 1 Моделирование контактного взаимодействия в щеточном уплотнении Modeling contact interaction in a brush seal Взаимодействие (контактная пара) Технология Типы КЭ в ANSYS Щетинки - ротор (К1) Контакт узла и поверхности CONTA175, TARGE170 Щетинки - опорная пластина (К2) Контакт пересекающихся балок CONTA176, TARGE170 Между щетинками (К3) Контакт параллельных балок CONTA176, TARGE170 Для моделирования взаимодействия между элементами уплотнения использованы контактные технологии, реализованные в комплексе ANSYS (табл. 1). Схема назначения контактных пар показана стрелками на рис. 3. Ротор и опорная пластина представлены абсолютно жесткими «целевыми» элементами. Верификация созданной модели выполнена на основе сравнения полученных с ее помощью результатов с экспериментальными и расчетными данными работ [15; 20]. Размеры и параметры тестового щеточного уплотнения (ТЩУ) приведены в табл. 2. Введенные обозначения поясняются рис. 1 и 6. Для моделирования каждой щетинки использовано 20 балочных конечных элементов. Расчеты выполнены в геометрически нелинейной постановке, влияние которой начинает проявляться при смещениях ротора больше 0,6 мм. Для учета контактного взаимодействия с трением использован комбинированный алгоритм (метод множителей Лагранжа для контактных условий по нормали и метод штрафа для условий по касательной); параметр точности - 1 мкм, величина относительной контактной жесткости - 0,01. Зазор между щетинками принят равным , как и в работе [15]. В целом результаты, полученные по созданной расчетной модели, удовлетворительно согласуются с известными экспериментальными и расчетными данными, что позволяет сделать вывод о допустимости использования этой модели для расчетной оценки жесткости щеточного уплотнения. Апробация созданной расчетной модели выполнена на примере расчета жесткости щеточного уплотнения Таблица 2 / Table 2 Основные параметры щеточных уплотнений Main parameters of brush seals Параметр ТЩУ БЩУ ГЩУ Внутренний диаметр щеточного пакета , мм - 435 Наружный диаметр щеточного пакета , мм - 461,4 470,6 Внутренний диаметр опорной пластины , мм - 436,3 Толщина щеточного пакета , мм - 1,85 Длина щетинки , мм 25,4 18,5 24,3 Зазор под опорной пластиной , мм 3,3 0,65 Угол наклона щетинок , 45 45 Диаметр щетинки (проволоки) , мм 0,142 0,100 Модуль упругости материала проволоки , ГПа 207 214 Жесткость щетинки по (4.1.7.1) , Н/м 1,521 0,996 0,440 П р и м е ч а н и е : ТЩУ - тестовое уплотнение [15; 20]; БЩУ и ГЩУ - базовое и гибкое (модифицированное) уплотнения [9]. Рис. 6. Эскиз базового уплотнительного узла Fig. 6. Basic seal assembly sketch большого диаметра (см. табл. 2) с учетом действия давления газового потока. Данное уплотнение было спроектировано в НИИД совместно с ЦИАМ. Коэффициенты трения приняты равными ; ; [15]. Рассмотрен сегмент из пяти щетинок с начальным зазором d между щетинками, равным 5 мкм. Расчеты выполнены в геометрически линейной постановке ввиду небольшой величины смещения ротора согласно внешней кинетике зазора в узле. Рассчитанные зависимости силы реакции щеточного пакета, приходящейся на одну щетинку, приведены на рис. 8. Там же показаны линейные зависимости , рассчитанные по полуэмпирической формуле (2) и используемые в упрощенной полуаналитической постановке задачи определения износа щетины [9]. Как видим, при наличии перепада давления газа на уплотнении имеет место значительное (на порядок) расхождение результатов упрощенной модели и пространственной модели. Этот эффект был отмечен и для тестового щеточного уплотнения (см. рис. 7). а б Рис. 7. Определение статической жесткости щеточного пакета без учета давления газового потока: а - модель из девяти щетинок; б - эксперимент и расчеты [15; 20] Fig. 7. Determination of the static stiffness of the brush pack without taking into account the pressure of the gas flow: а - a model of nine bristles; б - experiment and calculations [15; 20] а б в г Рис. 8. Сила реакции щеточного пакета для двух вариантов ЩУ большого диаметра: а - базовое ЩУ, ; б - «гибкое» ЩУ, ; в - базовое ЩУ, ; г - «гибкое» ЩУ, Fig. 8. Reaction force of the brush pack for two brush options large-diameter seals: а - basic BS ; б - «flexible» BS ; в - basic BS ; г - «flexible» BS На рис. 8 также указаны максимальные значения эквивалентных напряжений в щетинках, которые не превосходят , что подтверждает работу в упругой области. Величина перепада давления характеризует условия работы уплотнения на режиме максимальной скорости износа по полетной программе («малый газ» [9]). Полученные зависимости силы реакции щеточного пакета от величины натяга / смещения ротора (петли гистерезиса) позволяют уточнить жесткость щетины и скорость ее износа при работе в составе двигателя и дать прогноз эффективности щеточного уплотнения с течением наработки. Прогноз эффективности опирается на оценку расходной характеристики щеточного уплотнения, которая стоится по моделям работ [22-27]. Полученные результаты также могут быть использованы для уточнения модели комбинированных щеточных уплотнений [28-30]. Заключение Учет пространственного контактного взаимодействия элементов щеточного уплотнения значительно влияет на величину и характер изменения силы реакции щеточного пакета и, соответственно, на его жесткость. При наличии перепада давления газа на уплотнении результаты расчетов по упрощенной полуаналитической модели и пространственной модели отличаются на порядок. Необходимо учитывать этот эффект при расчете износа щеточного уплотнения и прогнозе его эффективности. Характерной особенностью зависимости силы реакции щеточного пакета от величины смещения ротора является наличие гистерезиса в цикле «нагрузка - разгрузка» из-за трения в системе. Воздействие потока газа на щетинки значительно увеличивает размах петли гистерезиса и меняет ее форму. На размах петли также сильно влияет коэффициент трения пары «щетинки - ротор».

About the authors

V. S Bondar

Moscow polytechnical university

E. V Rodin

Moscow polytechnical university

References

Statistics

Views

Abstract - 381

PDF (Russian) - 362