CC BY
130
УДК 629.73.083(075.8)
ТЕМПЕРАТУРНАЯ ДИАГНОСТИКА РОТОРНЫХ ПОДШИПНИКОВ ГАЗОТУРБИННЫХ ДВИГАТЕЛЕЙ
Р.В. БИРЮКОВ, Ю.В. КИСЕЛЕВ
Статья представлена доктором технических наук, профессором Кузьмичевым В.С.
В статье рассмотрены вопросы диагностирования подшипников газотурбинных двигателей (ГТД) с использованием штатно измеряемых тепловых параметров. Представлена математическая модель теплового состояния ра-диально-упорных подшипников газотурбинного двигателя наземного применения, и показаны примеры ее применения для диагностирования.
Ключевые слова: диагностирование, подшипники, критерии подобия, тепловыделение.
По опыту эксплуатации известно, что одним из компонентов, ограничивающих ресурс ГТД, являются подшипники роторов, отказы которых приводят к серьезным последствиям. Поэтому актуальной является задача обеспечения эффективного контроля и диагностирования подшипниковых опор ГТД в эксплуатации с учетом ограниченного перечня измеряемой на двигателе информации.
Развитие неисправностей подшипников сопровождается изменением вибросостояния, повышенным тепловыделением и изменением в масле концентрации металлов [1].
В данной статье рассматриваются вопросы создания методики контроля и диагностирования подшипников по штатно измеряемым функциональным параметрам, связанным с температурным состоянием подшипников, запись которых осуществляется в систему регистрации все время работы двигателя.
Перечень измеряемых параметров, отражающих тепловое состояние подшипника, различен для разных ГТД и, как правило, включает в себя температуры масла на входе и выходе из двигателя (иногда - на выходе из опор (отдельно по опорам)) и в некоторых случаях - непосредственно измерение температур подшипников (наружных обойм). Такой замер, в частности, реализован для наиболее нагруженных шариковых подшипников на двигателях НК-36СТ и ряде других типов двигателей [2]. По опыту эксплуатации непосредственно измеренная температура подшипника является весьма чувствительной к изменению технического состояния.
В системах контроля и диагностирования применяется допусковый контроль температур, однако он позволяет выявить изменение состояния лишь на поздней стадии развития неисправности. В ряде систем диагностирования [2] применяются алгоритмы трендового анализа температурных параметров. Однако применение алгоритмов трендового анализа затруднено ввиду довольно существенных колебаний температуры подшипника при сменах режима и условий работы (т.е. не вызванных изменением технического состояния).
Если учесть каким-либо образом влияние на измеряемую температуру режима и условий работы, то изменение температуры может быть объяснено изменением технического состояния. А по величине отклонения измеренной температуры от некоторого эталона исправного состояния, с учетом накопленной статистики и опыта эксплуатации конкретного двигателя можно судить о фактическом состоянии подшипника.
Рассмотрим влияние режима и условий работы на температуру подшипников.
Наиболее полно вопрос температурного состояния подшипников ГТД рассмотрен в работе [3]. Температура подшипников зависит от нескольких основных факторов:
- потребляемой мощности, иначе говоря, тепловыделения в подшипнике;
- внешнего подвода тепла (через элементы конструкции опоры или из-за обдува подшипника горячим воздухом);
- прокачки масла и температуры масла на входе.
Рассмотрим вопрос определения потребляемой шариковым подшипником мощности. В работе [3] принято разделение этой мощности Q^ на две составляющие: потери на трение Q тр и гидродинамические потери (на перемешивание смазки) QraaF,. Тогда потребляемая (рассеиваемая) мощность
QX = Qтр + Qгидр •
Гидродинамические потери в работе [4] рассчитываются с использованием критериального уравнения
Q гидр=16,6405^е-1'25^г-1^ Vu3 ,
n(d -dm ■ cosy) u ■ da
где u - окружная скорость сепаратора u = —--- ■ n ; Re =--— критерий Рейнольдса;
120 v
v
Pr =--критерий Прандтля; - число тел качения; ёш - диаметр тела качения; р - плотность
a
смазки.
Выражение для гидродинамических потерь преобразуем следующим образом Q^ = 16,6-105-Re-1'25-Pr-1-z ш•dш2•р•u3 = abRe-1,25^-1 •u3,
5 2
где а1 - постоянный коэффициент, равный а1 = 16,6^10 zш•dш •р (не учитывая зависимость р = Д1;м вх)). Выражение для окружной скорости сепаратора перепишем следующим образом
п( d4rn - du ■ c0sy)
u = —--- ■ n = 51 ■ n,
120
где s1 - постоянный коэффициент, описывающий геометрию подшипника.
С учетом этого перепишем выражение для критерия Рейнольдса, дополнительно «внеся» в значение s1 умножение на диаметр шарика
„ u ■ du 51 ■ n
Re =-ш =-.
v v
Вязкость масла можно учесть через температуру масла на входе, тогда с учетом критерия Прандтля (который в нашем случае перепишем через температуру масла на входе) и видоизмененного выражения для окружной скорости сепаратора можно записать выражение для гидродинамических потерь
Q гидр = a1 •Re-1, 25Фг-1 •u3 = a1^((s1^n)/t мвхо1)в1- t мвхо2 (sbn)3,
где о1, о2, в1 - некоторые постоянные коэффициенты.
Тогда, раскрыв скобки и внеся значения всех постоянных величин в коэффициент а1, показатель степени при оборотах в1 и при температуре масла с1, получим выражение для гидродинамических потерь
Qгидр = a1 • t мвх • n •
Потери на трение определяются следующим образом:
Qтр = С1 • z ш •d ш2•р•u3;
С1 =9(Re, Eu, Pr)= М^Ш-5^214- Eu0,287 •Pr0,44,
Pc ш а
где Eu =-ср-Ш—^ - критерий Эйлера; Р =--средняя нагрузка на одно тело качения;
p\u ■ dm) - Zrn
А - осевая сила, действующая на подшипник (если подшипник нагружен прежде всего осевой силой).
Для радиально-упорных подшипников нагрузку на тело качения определяет главным образом осевая сила, зависящая от конкретной геометрии исполнения системы осевой разгрузки и режима работы двигателя. Нетрудно отметить, что это выражение может быть преобразовано подобно выражению для гидродинамических потерь, и с учетом того, что критерий Эйлера
также может быть выражен через обороты (так как от них зависит и окружная скорость, и осевая сила), можно записать аналогично
Отр =а2 • 1 мвх 11 •
В результате рассеиваемая подшипником мощность
ГЛ лс с 1Г>5 т> -1,25 т» -1 I 1 Л т 1П-5 т> 0,214 ^ 0,287 „ 0,44ч А 2 3
О^ = (16,6-10 -ке •Рг +14,7-10 -ке • Би -Рг ' )• ъш • аш -р-и .
С учетом проведенных преобразований это выражение запишется в виде
О х = а1 • 1мвх • п + а2 • 1мвх • п .
С учетом прокачки и входной температуры масла (но без учета внешнего подвода тепла) температура подшипника может быть определена формулой
3600• От
1 = 1 +-.
п м вх /Т
Ср • дм
Здесь я м [кг/ч] - прокачка масла.
Для учета внешнего подвода тепла к подшипнику по методике [4] определяется добавка Д1п
Д1 п = Х-О вн-1 ,
где х - экспериментально определенный коэффициент; Овн - тепловой поток внешнего подвода тепла; I - характерный размер.
Проделаем для итоговой формулы 1п преобразования, аналогичные по своей идее проведенным ранее для потерь на трение и гидродинамических потерь. Тогда можно записать выражение для температуры подшипника в виде
1 п = 1 м вх+(а3^м> (а1 • 1мвхс1- пв1+ а2 • 1 мвхс2- пв2)+а4- пв4.
Учтя значение а3 в а1 и а2
1 п = 1м вх+ (а1 • 1 мвхс1- пв1+ а2 • 1мвхс2- пв2)/ ям +а4- пв4. (1)
Здесь с помощью слагаемого а4-пв4 учитывается внешний подвод тепла; а3, а4, в4 -постоянные коэффициенты.
Таким образом, фактическая температура подшипника будет зависеть от входной температуры масла, прокачки масла через опору (или связанной с ней величины) и частоты вращения ротора двигателя, причем эта зависимость степенная. Однако напрямую использовать эту формулу для определения некоторой «эталонной» температуры, с которой следует сравнивать измеренную 1 п не представляется возможным. Это объясняется тем, что не известны точно прежде всего действующая осевая сила (как указано в начале статьи) и внешний подвод тепла.
В связи с этим предложена тепловая модель подшипника, учитывающая общие законы влияния различных факторов на 1п и индивидуальные особенности конкретного типа и экземпляра двигателя.
Прорабатывались различные виды формул, связывающих температуру подшипника с влияющими факторами. Наилучший результат получен для формулы, соответствующей преобразованному виду формулы (1)
г = п а • г ь • в а • с, (2)
п отн _ р мвх отн _ р м _ отн '
где 1 п - расчетная температура соответствующего подшипника; п отн_р - физические обороты соответствующего ротора, отнесенные к некоторому выбранному характерному режиму работы; 1 мвхотн_р - температура масла на входе, отнесенная к некоторому выбранному характерному режиму работы; Ом_отн - относительная прокачка масла; А, Ь, с, ё - эмпирически получаемые коэффициенты.
Эта модель позволяет вводить поправки к измеренной температуре для учета влияния режима и условий работы (т.е. приводить температуру к некоторому режиму или, иначе говоря, корректировать ее). Кроме того, разработана методика автоматизированного получения числовых значений коэффициентов в этой модели, и таким образом она может быть адаптирована к любому двигателю (авиационному или индустриальному), где осуществляется измерение необходимого набора параметров.
В данной работе расчет коэффициентов проводился по данным, зафиксированным при эксплуатации трехроторного ГТД наземного применения. Ротор НД, СД газогенератора - тре-хопорные, ВД - двухопорный. На каждом из радиально-упорных подшипников (НД, СД, ВД) осуществляется измерение температуры наружной обоймы.
В расчетах использовались штатно измеряемые параметры: частоты вращения роторов, входная температура масла и два параметра, определяющие фактическую прокачку масла через опоры - давление масла на входе в двигатель Рмвх и давление в объединенной системе суфлирования Рсуфл. Величина Рмвх-Рсуфл напрямую связана с перепадом давления на форсунках подачи масла, т.е. при заданных гидравлических характеристиках системы форсунок определяет подачу масла в опору
О = ц- Г • у]2р-АР или Ом = к1 • у]Рмвх - Рсуфл .
Относительные значения параметров п
отн_р Л м вхотн_р
вычислялись следующим образом
Ппр_р Пзад/Пвыб, (3)
где Пзад - задаваемое значение параметра; Пвыб - базовое значение параметров выбранного для приведения наиболее характерного режима работы.
Относительная прокачка масла определялась из следующего соотношения
С*™ = ^^ =0=412 - ^(Рмвх-Рсуфл)
Значения показателей степеней а, Ь, с, с! рассчитаны методом наименьших квадратов по данным параметров работы конкретного экземпляра двигателя.
При выполнении расчетов отмечено, что полученные формулы дают наилучший результат на стационарных режимах и примерно дают значения температуры на всех режимах работы, в том числе на запуске и на переходных режимах.
Коэффициент С в формулах отражает некоторую «базовую» температуру, которую подшипник будет иметь при указанных выше в формулах отнесения выбранных базовых значениях параметров.
При увеличении режима работы двигателя или температуры масла на входе от выбранного базового режима отнесенные значения параметров буду становиться больше 1, и температура подшипника по модели будет увеличиваться. По отношению к прокачке масла будет наблюдаться обратная зависимость.
С использованием выражения (2) проведено сравнение получаемых по модели значений параметров с измеренными.
На рис. 1 приведено сравнение температур подшипника НД для конкретного экземпляра двигателя за период в 14 дней.
Ш нд °С
105 100 95 90 85 80
Измеренное значение
ЛИ.11^1111^ -
Расчетное по модели значение
75
70 ■
№ замера
ННННННННН
Рис. 1. Измеренное и модельное значения температуры упорного подшипника ротора НД
На рис. 1 видно хорошее совпадение модельного и измеренного значений.
Также проведено сравнение измеренных температур подшипников с расчетными температурами для подшипников НД, СД и ВД нескольких экземпляров двигателей. Проведенные расчеты и сравнения показали хорошее совпадение расчетных и измеренных температур, модель адекватно отражает влияние на температуру режима и условий работы. Следовательно, модель можно применить для целей диагностики.
Для решения задач диагностики можно проводить анализ тенденции изменения разницы между расчетным (по модели) и измеренным значением, однако более удобно использовать полученную модель для «корректировки» или, другими словами, приведения измеренных значений в соответствие с фактическими условиями работы. В этом случае скорректированное (приведенное) значение в идеале при неизменном техническом состоянии должно оставаться постоянным.
Корректировка значений температур подшипников 1;пНд, 1;пСд, 1пВд осуществляется по формуле, полученной преобразованием формулы (2)
Г = Г • п ~а ■ г ~ъ ■ О ~а ■ с (4)
п кор п измер отн _ р мвх отн _ р м _ отн ' V /
При использовании данной формулы устраняется недостаток модели в части не учета индивидуальности измеренного температурного состояния подшипников для разных двигателей, а используется часть формулы, отражающая влияние на температуру режима и условий работы. Очевидно, что при возникновении неисправностей и повышенном тепловыделении в подшипнике скорректированное значение температуры не будет маскировать неисправность.
На рис. 2, 3 приведен пример применения формулы для коррекции значений температур подшипников двигателя по данным, полученным за период 14 дней.
За указанный период времени температура воздуха колебалась в пределах +3,5...+18,5°С, температура масла - в пределах +55,6...+74,2°С, частота вращения ротора НД - в пределах 4670.. .4860 об/мин.
Ш нд °С
85
80
75
/и
Рис. 2. Измеренное и скорректированное значения температуры упорного подшипника ротора НД
Ш сд °С
- Измеренное
« \
Скорректированное значение
85
№ замера
Рис. 3. Измеренное и скорректированное значения температуры упорного подшипника ротора СД
Как видно из рис. 2, 3 скорректированное по внешним условиям и режиму работы значение температур в исправном состоянии имеет существенно меньший разброс по сравнению с измеренной температурой. Это дает основание для использования скорректированной температуры для трендового анализа.
На рис. 4 приведем пример изменения температуры подшипника одного из двигателей в ходе развития неисправности (разрушение подшипника ротора ВД). Можно отметить два участка на графике - участок плавного роста температуры перед аварией (с 120°С до ~ 160°С, начало плавного роста за 14 минут до разрушения) и участок скачкообразного увеличения температуры при разрушении подшипника (с ~160°С до 310°С). Из рис. 4 хорошо видно существенное изменение скорректированного значения температуры при развитии неисправности. Также изменялось и измеренное значение температуры, однако по этому изменению нельзя было однозначно судить о развитии аварии на участке плавного роста температур, так как температура оставалась в пределах допуска. Введение контроля скорректированной температуры позволяет надежно определить начальный момент разрушения подшипника путем использования тренд анализа или введением более жесткого допуска.
Рис. 4. Изменение температуры подшипника ротора ВД в ходе развития аварии
Подобное рассмотрение можно провести не только для температуры подшипника, но и для температуры масла на выходе из опоры, для теплоперепада (разности температур на входе и выходе) на опоре, для теплоотдачи двигателя в масло.
Таким образом, полученная эмпирическая модель может служить основой для построения алгоритмов диагностирования с использованием измеряемых температурных параметров опор роторов ГТД.
Работа выполнена при финансовой поддержке Правительства Российской Федерации (Минобрнауки) на основании постановления Правительства РФ № 218 от 09.04.2010 г.
ЛИТЕРАТУРА
1. Люлько В.И. Эксплуатация авиационных двигателей по техническому состоянию (теория и практика). - М.: МГУ, 2002.
2. Зарицкий С.П. Автоматизированная система диагностирования двигателя НК-36СТ / С.П. Зарицкий, В.Н. Исланов, В.Б. Коротков, В.Н. Михнович, Ю.Н. Тарасенко // Энергодиагностика / Диагностика и надежность: сб. тр. междунар. конф. - М.: Изд-во «Энергия», 1995. - Т. 2.
3. Демидович В.М. Исследование теплового режима подшипников ГТД. - М.: Машиностроение, 1978.
TEMPERATURE DIAGNOSTICS OF GAS TURBINE ENGINE BEARINGS
Birukov R.V., Kiselev Y.V.
The objective of the current research was to develop methodology for diagnosing industrial gas turbine engine bearings using the standard performance parameters. This paper presents mathematical thermal model of combined thrust and radial bearing and provides the model application examples for diagnostics.
Key words: diagnostics, bearings, similarity criterion, thermal emission.
Сведения об авторах
Бирюков Роман Владимирович, 1990 г.р., окончил СГАУ (2013), аспирант СГАУ, область научных интересов - диагностика и регулирование авиационных двигателей.
Киселев Юрий Витальевич, 1948 г.р., окончил КуАИ (1972), доцент, кандидат технических наук, доцент кафедры эксплуатации авиационной техники СГАУ, автор более 60 научных работ, область научных интересов - диагностика авиационной техники.
