А. С. Гишваров, А. А. Рыжов, А. Г. Алгушаев, А. Р. Фатыхова
УДК 629.735
ОБОСНОВАНИЕ пРОГРАММ УСКОРЕННЫХ ИСпЫТАНИЙ НА НАДЕжНОСТЬ И РЕСУРС ГАЗОТУРБИННЫХ ДВИГАТЕЛЕЙ И СИЛОВЫХ УСТАНОВОК МНОГОВАРИАНТНОГО пРИМЕНЕНИЯ
Рассматриваются последовательность обоснования программы ускоренных испытаний на надежность и ресурс маршевого авиационного газотурбинного двигателя, устанавливаемого на беспилотном летательном аппарате, предназначенном для эксплуатации по шести профилям полета и обоснование программы испытаний авиационной вспомогательной газотурбинной силовой установки, эксплуатируемой на трех типах самолетов. Газотурбинный двигатель многовариантного применения, ускоренные ресурсные испытания, исследование, оптимальная область, программа испытаний, надежность, ресурс
ОБОСНОВАНИЕ ПРОГРАММЫ УСКОРЕННЫХ ИСПЫТАНИЙ
на надежность И РЕСУРС
АВИАЦИОННОЙ ВСПОМОГАТЕЛЬНОЙ ГАЗОТУРБИННОЙ СИЛОВОЙ УСТАНОВКИ
Авиационная вспомогательная газотурбинная силовая установка устанавливается на трех типах самолетов и эксплуатируется по своим циклам нагружения. Общее число циклов нагружения установки в каждом варианте эксплуатации равно соответственно 2130, 1028 и 900 циклам [1-3].
Требуется определить объем, режимы и длительность ускоренных ресурсных испыта-ний, если по априорным данным известно, что надежность силовой установки определяют:
• малоцикловая прочность дисков ротора Пп;
• длительная прочность рабочей лопатки первой ступени турбины П22;
• контактная прочность радиально-упорного подшипника опоры П33;
• контактная прочность ведущей шестерни редуктора П43.
Обобщенная модель расходования ресурса рассматриваемой силовой установки имеет вид
[4-6]:
Пп [ Р011; г; (т), Тн (т), n (х), NT (т), Nz (т); х ] = 0;
П22 [P0.22; Т; (х), Т (т), n (т), N (О, Nz (х); х ] = 0;
П33 [ P0.33; Т* (О, n (х); х] = 0; п4з[р04з; n (О, NT (х); х] = 0; R(x) = [т;(х), Т(х), n(х), N(т),Nz(х )]т Є Gr;
P0 = idem,
Контактная информация: 8 (347) 273-79-54 Работа выполнена при финансовой поддержке Минобрнауки РФ.
где Тг* - температура газа перед турбиной, К; Тн - температура воздуха на входе в установку, К; п - частота вращения ротора, %; N - загрузка генераторов постоянного и переменного токов, кВт; N - число циклов нагружения; Р0 - область исходного состояния установки, соответствующая техническим требованиям на его изготовление.
Область определения параметров режима нагружения двигателя во всех трех вариантах применения соответствует ограничениям:
650 < Т* < 1080 К; 0 < Ыт < 100 кВт;
Ок% : <¡223 < Г < 320 К; 900 < N < 2556;
60 < п < 105 %.
Уровень адекватности моделей характеризуется остаточной дисперсией, а значения весовых коэффициентов приняты следующими [2]:
b11 = 0,2;
Ь22 = 0,33;
b33 = 0,25; b43 = 0,17.
Значения коэффициентов р характеризующих значимость элементов, выбраны из условия, что отказ элемента ведет к отказу силовой установки в целом [3]:
р*. = 0,25; / = 1,4; ] = 1,3.
Результаты оптимизации параметров ускоренных ресурсных испытаний приведены в табл. 1, откуда видно, что общее число возможных вариантов испытаний силовой установки равно 6. Данные в столбцах 3, 4 и 5 табл. 1 означают номера испытываемых силовых установок. Например, данные второй строки означают, что надежность силовой установки в 1-м и 2-м эксплуатационных вариантах применения проверяется в одних испытаниях (столбец 3), а его надежность в 3-м эксплуатационном варианте применения - в других испытаниях (столбец 4) и т.д.
Для сравнения в строке 6 приведены параметры испытаний, оптимизация которых проведена без учета критерия, характеризующего уровень дифференцированного «зачета» в ускоренных ресурсных испытаниях эксплуатационных вариантов применения силовой установки, начиная от наименее нагруженного и кончая наиболее нагруженным вариантом, и позволяет, даже в случае отказа силовой установки в испытаниях, не браковать всю партию (в зачет которой проводятся испытания), а признавать изделия годными к эксплуатации по тем вариантам, эксплуатационная повреждаемость которых выработана до момента отказа [4, 7]. Видно, что для рассматриваемой силовой установки учет данного критерия позволил существенно повысить уровень дифференцированного «зачета» эксплуатационных вариантов применения (при этом уровень «недифференцируемости» снизился с 47,1 до 0).
Из множества Парето-оптимальных значений N Яик(т), тк за окончательный выбран вариант под номером 5, для которого получены наилучшие значения критериев эффективности Ф Ф2, Ф3 и Ф4.
Согласно данному варианту ускоренные ресурсные испытания проводятся на одной силовой установке с реализацией циклов нагружения, приведенных на рис. 1.
№ N. N. Ти.К. «Зачет» в испытаниях эксплуатационных вариантов при тДІЧц) Общая дли- тель- ность, час Значения критериев эффективности (Ф;= ПІ 111)
к=1 к=2 к=3 к=1 к=2 к=3 1-го 2-ю 3-го ф| ф2 ф, ФJ
1 3 (1) (2) (3) 1674 749 594 1674 (2556) 749 (1235) 594 (1080) 3017 0 0 1 0,74
2 2 (Ь2) (3) - 1674 594 - 1674 (2556) 749 (1235) 594 (1080) 2268 0 0 0,67 0,55
3 2 (М) (2) - 1674 749 - 1674 (2556) 749 (1235) 594 (1080) 2423 0 0 0,67 0,60
4 2 (23) (1) - 749 1674 - 1674 (2556) 749 (1235) 594 (1080) 2423 0 0 0,67 0,60
5 1 (1АЗ) - - 1674 - - 1674 (2556) 749 (1235) 594 (1080) 1674 0 0 0,38 0,41
6 1 (1,23) - - 1674 - - 1674(2556) 1674 0 47,1 0,33 0,41
7 3* (1) (2) (3) 2140 1370 1370 2140 (2556) 1370 (1235) 1370 (1080) 4880 0 0 1 1,20
Примечание: тз - время «зачета» эксплуатационного варианта, час; (Ыц) - число циклов нагружения в испытаниях, соответствующее времени «зачета» тз; (1), (2),..., (1,2,3) - множество эксплуатационных вариантов, проверяемых в каждых испытаниях; * - испытания полностью соответствуют эксплуатационным
Рис. 1. Циклы нагружения авиационной газотурбинной силовой установки в оптимальном варианте ускоренных ресурсных испытаний:
І I - поле допуска на параметр режима;
11’, 12’,... - длительность работы силовой установки на участках нагружения (в мин)
В процессе ускоренных ресурсных испытаний «засчитываются»:
• первый эксплуатационный вариант применения на 1674 часе (2556 циклов нагружения);
• второй эксплуатационный вариант применения на 749 часе (1235 циклов нагружения);
• третий эксплуатационный вариант применения на 594 часе (1080 циклов нагружения).
Таким образом, для рассматриваемой силовой установки, ускоренные ресурсные испытания позволяют обоснованно сократить время на проверку ее надежности в К = 4880 / 1674 « 3 раза.
ОБОСНОВАНИЕ ПРОГРАММЫ УСКОРЕННЫХ РЕСУРСНЫХ ИСПЫТАНИЙ НА НАДЕЖНОСТЬ И РЕСУРС МАРШЕВОГО ГАЗОТУРБИННОГО ДВИГАТЕЛЯ
Маршевый двигатель устанавливается на беспилотном летательном аппарате, эксплуатируемом по двум профилям полета в каждой из трех климатических зон «САУ», «ЗИМА» и «ЛЕТО», т.е. по 6 вариантам применения. Режимы нагружения задаются температурой (Т *) и давлением (Рет*) воздуха на входе в двигатель, частотой вращения ротора (п), а также длительностью нагружения (т). Циклы нагружения двигателя приведены на рис. 2-7.
На предприятии-изготовителе надежность двигателя проверяется по программе серийных ускоренных ресурсных испытаний (рис. 8).
Рассмотрим возможность повышения эффективности программы серийных ускоренных ресурсных испытаний двигателя . т'
и,5 О
■-
/ ч
0,5
0,5
0,5
Рис. 2. Цикл нагружения двигателя в 1-м эксплуатационном варианте применения
т; 1
0.5 /
/
0
р:
0,5
0
0,5
0
■ V
0,5
0,5
0,5
Рис. 3. Цикл нагружения двигателя во 2-м эксплуатационном варианте применения
т;
0,5
О
/
У
0,5
О
0,5
О
0.5
0.5
0,5
Рис. 4. Цикл нагружения двигателя в 3-м эксплуатационном варианте применения
т вх р;
05 0,5
/
0 / 0
0,5
0
0 0,5т
0 0,5:
0 0,5т
Рис. 5. Цикл нагружения двигателя в 4-м эксплуатационном варианте применения
Рис. 6. Цикл нагружения двигателя в 5-м эксплуатационном варианте применения
т;
0,5
0
р.;
0,5
0
п
0,5
0
0 0,5 I
0 0,5 Г
Рис. 7. Цикл нагружения двигателя
в 6-м эксплуатационном варианте применения
Исходными данными для рассматриваемого двигателя являются [3, 4, 8]:
• данные по эксплуатации, приведенные в виде типовых полетных циклов на рис. 2-7;
• типовой цикл нагружения двигателя в стендовых испытаниях (рис. 8);
• результаты термогазодинамических и прочностных расчетов двигателя.
Элементами, определяющими надежность двигателя в эксплуатации, повреждаемость которых обусловлена влиянием режимов и длительности нагружения, являются:
• рабочая (ЛТ) и сопловая (СА) лопатки турбины;
• рабочая (ЛК) и направляющая (НА1) лопатки первой ступени компрессора;
• направляющая лопатка последней ступени компрессора (НА2);
• корпус камеры сгорания (КС);
• радиально-упорный подшипник ротора (Под);
• ведущая шестерня привода агрегатов (Шк).
Определяющими характеристиками расходования ресурса являются длительная статическая, контактная и изгибная прочность.
Область определения параметров режима нагружения соответствует эксплуатационной и задана ограничениями:
240 < Тх* < 623 K; на режиме «О»: 0,7 < P ** < 10 кГ/см2;
’ — вх —
n = 94% при T * < 339 K.
вх
К.
0,5
0,5
Т„
0,5
У V
у
0
0.5
Рис. 8. Цикл нагружения двигателя в серийных испытаниях: P * , T * , п’ - относительные
вХ 7 вх 7
значения соответственно давления и температуры воздуха на входе в двигатель, а также частоты вращения ротора
ФОРМИРОВАНИЕ МАТЕМАТИЧЕСКОЙ МОДЕЛИ ПОВРЕЖДАЕМОСТИ ДВИГАТЕЛЯ
Обобщенная модель повреждаемости двигателя представляет собой систему моделей повреждаемости его элементов: ЛТ, СА, ЛК, НА1, НА2, КС, Под и Шк.
Поскольку Р0 = idem, то формула для расчета повреждаемости примет вид:
п„=1
dt
[Т*,Р*,п]
где і - порядковый номер элемента узла двигателя (і = 1,8);. - порядковый номер характеристики расходования ресурса (/ = 1,3); х*.. - время до разрушения или наступления параметрического отказа элемента узла двигателя.
Длительная прочность элементов узлов (ЛТ, СА, ЛК, НА1, НА2 и КС) описывается параметрической зависимостью Ларсена-Миллера:
с = А - ВР, (1)
где с - предел длительной прочности, кГ/мм2; А, В - коэффициенты; Р - параметр Ларсена-Миллера, определяемый по формуле:
Р = Т (20 + ^ т), (2)
где Т - температура элемента, К; т - длительность работы, час.
При известном напряжении с и температуре Т, с учетом соотношений (1) и (2), время до разрушения определится по формуле:
т* = 10
A -ст.
г___г_
B T
- 20
Коэффициенты А и В. определяются методом наименьших квадратов по кривым длительной статической прочности материалов деталей. Нагружение лопатки характеризуется напряжением слт и температурой Тлт, а повреждаемость лопатки турбины с учетом длительного статического нагружения определится по формуле:
&
209,23 -ст™, 010[ ТД.210^-20
= F [Т*, P*, n ];
= F[Т*, п].
В окончательном виде математическая модель повреждаемости двигателя записывается в виде системы уравнений:
сіі
10
209,23 -а
лт
7,02-10 -Т
-20
лт
ПСА J
о
dt
10
231,50-а
СА
7,709-10-Т,
-20
СА
Плк - | о
dt
10
114,45-а
ЛК
3,29-10 -Т
-20
ЛК
IL
П
НА2
= /■
(
-І
dt
114,45-с
НА1
3 29-10 ~ЪТ 10 А1
-20
dt
104,81 -а
на^-20
10 v 3,39210
+
+
I-
dt
220,41 -о
НА2_-20
10v 1,°98'10 -^^НА2
п
йі
КС
10
164,79-сткс 6,427-10 -3Т
-20 I
КС
П
йі
Под
10
п
60 - п йі
Шк
о 120-106 60 • п • К™ь
ФОРМИРОВАНИЕ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ ЭЛЕМЕНТОВ УЗЛОВ ДВИГАТЕЛЯ
Модели напряженно-деформированного состояния элементов узлов двигателя формируются в виде уравнений:
Я,(г) = [г,(г), Г2(г),...,гДг)]гГ =
—п Т
= Г[Твх,РЙХ,п]
связывающих напряжение и температуру эле-
ментов сЛТ, ТЛТ
CСА, TСА, СЛК
с
Т
сНА2, ТНА2 , сКС , ТКС , й, МКР с параметрами режи-
Г* т-»*
, Р , п :
с і(х) = f [ Т (х) Р*х (х) П (х)£
Т(х) = f [ Т » Г (х), п (х)].
Определение констант моделей напряженно-деформированного состояния проводится по результатам термогазодинамических и прочностных расчетов двигателя [4, 6, 8].
ФОРМИРОВАНИЕ МАТЕМАТИЧЕСКОЙ МОДЕЛИ ЭКСПЛУАТАЦИИ ДВИГАТЕЛЯ
Для рассматриваемого двигателя модель его эксплуатации представляет собой априорно заданные полетные циклы нагружения, изменение параметров в которых приведено на рис. 2-7.
ОЦЕНКА ЭКСПЛУАТАЦИОННОЙ ПОВРЕЖДАЕМОСТИ ДВИГАТЕЛЯ
Результаты оценки эксплуатационной повреждаемости элементов рассматриваемого двигателя, проведенной по «расчетной модели», приведены в табл. 2.
ОПТИМИЗАЦИЯ РЕЖИМОВ И ДЛИТЕЛЬНОСТИ УСКОРЕННЫХ РЕСУРСНЫХ ИСПЫТАНИЙ ДВИГАТЕЛЯ
Предварительными исследованиями было установлено, что рациональным является оптимизация параметров испытаний на шести участках нагружения = 6).
Таблица 2
Относительная повреждаемость элементов двигателя в эксплуатации
Наименование элемента Относительная повреждаемость элементов двигателя в эксплуатационном варианте применения (в %)
1 2 3 4 5 6
Рабочая лопатка турбины 100 44,6 47,4 7,8 40,4 16,4
Сопловая лопатка турбины 100 29,2 32,3 34,9 82,8 36,5
Рабочая лопатка первой ступени компрессора 103 92,1 100 109 10-16 10-10
Направляющая лопатка первой ступени компрессора 100 97,2 105 109 109 10-22
Направляющая лопатка последней ступени компрессора 100 0,01 104 103 0,01 107
Корпус камеры сгорания 0,01 100 0,02 0,01 103 104
Опорный подшипник ротора 78,8 0,06 0,05 18,0 100 0,08
Шестерня привода 100 83,7 85,1 30,6 26,0 48,2
Примечание: значения повреждаемости элементов приведены в % относительно их максимальной повреждаемости в эксплуатации
В процессе оптимизации на участках нагружения варьировались:
• температура воздуха на входе в двигатель Т* , К;
вху 7
• давление воздуха на входе в двигатель Р* х, кГ/см2;
• относительная частота вращения ротора п, %;
• длительность нагружения ХЦ , %.
При оптимизации, дополнительно к ограничениям, учитывались параметрические ограничения, связанные с возможностями испытательных стендов «А» и «Б»:
Стенд "А":
Стенд "Б":
343 < T* < 623 К;
1,3 < P* < 2,5 ^/см2; 80 < n < 102 %.
240 < T* < 623 К;
0,7 < P* < 10 ^/см2. 80 < n < 102 %.
Параметрические ограничения обоих стендов находятся внутри области определения «расчетной модели» двигателя.
В систему ограничений были включены также ограничения, характеризующие скорость изменения параметров в эксплуатационных условиях:
(d T* 75г)„ < max (5 T* 75т)э =1,5 К/с;
(d P* /дт)и < max (5 P* 75т)э =
= 0,025 (кГ/см2)/с;
(5й7дг)„ <max(5n/5x)3 = 0,25 %/с.
Поскольку ранее разработанные на предприятии серийные испытания проводятся на одном двигателе, то оптимизацию параметров опытных ускоренных ресурсных испытаний также проводили для одного двигателя.
По требованию Заказчика в испытаниях должна обеспечиваться проверка всех шести вариантов применения двигателя (означает, что каждый выпускаемый двигатель может быть использован в любом из шести вариантов). Поэтому критерий дифференцированного «зачета» эксплуатационных вариантов применения в ускоренных ресурсных испытаниях (критерий Ф2) в данном случае не рассматривался. По условиям задачи при Nu = 1 Ф3 = 176 = idem, т. е.
целевая функция выбора параметров ускоренных ресурсных испытаний имеет вид:
ф [N*, R*k, X *k, Nk ] = (АФ! + БФ4) ^ min,
Г 1 N« 8 3 \
ф1 = mini-N-t 2 t[5Пj]2• К-pl\;
t N и k=1 i = 1 j=1 J
Пи jjk [P0u , -X J - (max П, jjl [Po3 , ..^3 ])k —______________________l_____________•
j (та^П3 ijl [Po3, Ra ^)д3]) k ’
l
Ф2 e (0..^2 max ); Ф3 = idem
( N
Ф4 = ШШ £x„ k ZT3
^k=1 / v=1
(трПэ ¡лtPos,RfrX*,])k <
I
n„ jjk [ P0„, R ^XxJ <е- aijk; i = 18; j = 1,3; l = 1,6; П = F (P0u, R (x), xj; R (x) e Gr ; P0u = idem; A + D = 1,0; B = С = 0;
bj =pj =idem; о,- =1.
Результаты оптимизации параметров ускоренных ресурсных испытаний показали следующее.
Оптимальным для стенда «А» является проведение испытаний двигателя длительностью 1,94 тс (цикл испытаний показан на рис. 9). При этом обеспечивается гарантированная проверка надежности:
• лопаток соплового аппарата;
• рабочих лопаток первой и направляющих лопаток последней ступеней компрессора;
• корпуса камеры сгорания;
• ведущей шестерни привода агрегатов,
а при этом уровень гарантированной проверки надежности остальных элементов составляет;
• рабочей лопатки турбины - 34,3% (в серийных испытаниях - 16,6%);
• радиально-упорного подшипника ротора -77,7% (в серийных испытаниях - 0,2%).
Оптимальным для стенда «Б» является проведение испытаний двигателя длительностью 2,58 tc. При этом обеспечивается 100%-я проверка надежности всех рассматриваемых элементов двигателя.
Оценка эффективности ускоренных ресурсных испытаний через долю эксплуатационной повреждаемости, вырабатываемой в испытаниях в среднем одним элементом двигателя, показала, что опытные испытания в 4 раза эффективнее серийных, а оценка их эффективности
по количеству эксплуатационных вариантов, проверяемых в испытаниях — в 2,5 раза эффективнее серийных.
Л А
О 0,5
0 0,5 t
п'
0 0,5 г
Рис. 9. Цикл нагружения двигателя при испытаниях на стенде «А»: Р* , Т ,П - относительные
вх 7 вх7
значения соответственно давления и температуры воздуха на входе в двигатель, а также частоты вращения ротора
Таким образом, применение методики выбора оптимальных значений параметров ускоренных ресурсных испытаний для рассматриваемого двигателя многовариантного применения позволило существенно повысить эффективность ранее разработанной серийной программы испытаний.
СПИСОК ЛИТЕРАТУРЫ
1. Гишваров А. С., Идрисов Р. Т. Анализ факторов, влияющих на формирование обобщенного цик-
ла испытаний авиационных ГТД // Испытания авиационных двигателей: межвуз. научн. сб. Уфа: 1986. № 14. С. 24-31.
2. Гишваров А. С., Давыдов М. Н., Агеев Г. К. Основные задачи повышения эффективности ускоренных ресурсных испытаний двигателей и энергетических установок // Вестник УГАТУ Уфа. 2012. Т. 16, № 2 (47). С. 69-75.
3. Гишваров А. С., Агеев Г. К. Исследования авиационных силовых установок с применением планирования эксперимента. Уфа: УГАТУ, 2009. 215 с.
4. Гишваров А. С. Теория ускоренных ресурсных испытаний технических систем. Уфа: Изд-во АН РБ «Гилем», 2000. 350 с.
5. Guishvarov A. S., Kondratieva N. V, Minigaleev S. M. Development technology of the programs of reliability and length of life accelerated tests of aircraft engines // 36th AIAA //ASME/SAE/ASEE Joint Propulsion Conf. Huntsville, Alabama, 2000. P. 36-38.
6. Гишваров А. С., Тимашев С. А. Синтез ускоренных ресурсных испытаний технических систем многовариантного применения. Екатеринбург: УрО РАН, 2012. 188 с.
7. Guishvarov A. S. Development technology of the programs of reliability and length of life accelerated tests of aircraft engines // Proc. of the 1997 Cina-Russia Symp. on Aero-Engines. Nahjing, China, 1997. P. 47-52.
8. Гишваров А. С. Параметрическая оптимизация ресурсных испытаний энергетических установок летательных аппаратов многовариантного применения // Вестник УГАТУ. 2002. Т. 3, № 2. С. 46-66.
ОБ АВТОРАХ
Гишваров Анас Саидович, проф., зав. каф. авиац. двигателей. Дипл. инж.-мех. (УАИ, 1973). Д-р техн. наук по тепловым двигателям ЛА (УГАТУ, 1993). Те-орет. и эксперим. иссл. в обл. надежности, ресурса и испытаний техн. систем.
Рыжов Алексей Андреевич, проф. той же каф. Дипл. инж.-мех. (УАИ, 1959). Д-р техн. наук по двигателям ЛА (УГАТУ, 1997). Иссл.в обл. проектирования и экспериментальной отработки ГТД.
Алгушаев Айнур Гафурович, асп. той же каф. Дипл. инж.-мех. (УГАТУ, 2010). Иссл. в обл. прогнозирования и диагностики технических систем.
Фатыхова Анвария Раисовна, асп. той же каф. Дипл. инж.-мех. (УГАТУ, 2010). Иссл. в обл. прогнозирования и диагностики технических систем.