Расчет надежности и прогнозирование долговечности деталей газотурбинных двигателей

Расчет надежности и прогнозирование долговечности деталей газотурбинных двигателей

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ УКРАИНЫ

НАЦИОНАЛЬНЫЙ АЭРОКОСМИЧЕСКИЙ УНИВЕРСИТЕТ им. Н.Е. ЖУКОВСКОГО «ХАИ»

Расчетно-пояснительная записка

по дисциплине: «Надежность Авиационных двигателей»

Тема: РАСЧЕТ НАДЕЖНОСТИ И ПРОГНОЗИРОВАНИЕ ДОЛГОВЕЧНОСТИ ДЕТАЛЕЙ ГТД

ВВЕДЕНИЕ

Надежность - это один из основных показателей качества изделия, проявляющийся во времени и отражающий изменения, происходящие в изделии на протяжении всего его периода эксплуатации. Надежность, как свойство изделия, закладывается на этапе проектирования, реализуется при изготовлении и поддерживается в процессе эксплуатации. Следует всегда иметь в виду, что качественно проработанный проект является основой надежности будущего изделия.

Поведение реальных конструкций обусловлено взаимодействием целого ряда факторов, имеющих явно выраженный случайный характер. В связи с этим определение надежности конструкций невозможно без применения методов теории вероятностей и математической статистики.

На основании физики возникновения отказы могут быть разделены на две группы:

внезапные отказы, имеющие характер случайного выброса:

хрупкое разрушение;

превышение предела текучести в какой-либо точке детали, для которой остаточные деформации недопустимы;

возникновение слишком больших упругих деформаций.

постепенные отказы, возникающие в результате необратимого накопления повреждений в детали:

пластические деформации (деформации ползучести);

усталостные повреждения, ведущие к развитию усталостных трещин.

Таким образом, при оценке надежности деталей авиационных ГТД необходимо учитывать внезапные и постепенные отказы. Эти два вида разрушений в первом приближении можно считать независимыми друг от друга.

Изложенные выше допущения позволяют принимать в качестве основного показателя надежности детали вероятность безотказной работы (неразрушения), формула которой имеет вид:

где  и  - вероятность безотказной работы с учетом внезапных и постепенных отказов.

1. ОПИСАНИЕ КОНСТРУКЦИИ КОМПРЕССОРА ГТД

Компрессор авиационного газотурбинного двигателя предназначен для сжатия рабочего тела (воздуха) перед подачей его в камеру сгорания. Прототип проектируемого двигателя - ТВаД ТВ3-117 - оснащен осевым компрессором.

Компрессор - осевой, дозвуковой, двенадцатиступенчатый, состоит из:

корпуса передней опоры;

входного направляющего аппарата (ВНА);

корпусов компрессора с направляющими аппаратами (НА);

спрямляющего аппарата;

ротора компрессора;

первой и второй опор дв-ля.

Для повышения запасов газодинамической устойчивости (ГДУ) дв-ля при работе на нерасчетных режимах, он оборудован следующими конструктивными механизмами:

. Поворотными лопатками ВНА и НА первых 4-х ступеней;

. Двумя клапанами перепуска воздуха (КПВ) из-за 7-ой ступени.

2. РАСЧЕТ НАДЕЖНОСТИ ЛОПАТКИ КОМПРЕССОРА

Рабочие лопатки компрессора подвержены нагрузкам от газовых и цент - рубежных сил, которые вызывают статические и переменные напряжения.

Статические напряжения в лопатках в основном складываются из напряжений растяжения от центробежных сил и изгибных напряжений от газовых и центробежных сил, вызванных наличием выносов центров тяжести для компенсации изгиба от газовых сил.

В качестве исходных данных используем результаты курсового проекта по курсу «Теория и расчет лопаточных машин».

Исходные данные

материал: ВТ3.

температура:- ;

Предел длительной прочности  МПа.

Коэффициент вариаций предела длительной прочности .

Плотность материала кг/.

Мощность ступени Вт.

Число лопаток .

Угловая скорость  рад/с.

Высота лопатки м.

Втулочный радиус м.

Средний радиус м.

Периферийный радиус м.

Расход газа кг/с.

Осевые составляющие абсолютной скорости газа на входе:

м/с;

м/с.

Давление на входе и выходе ступени:

Па;

 Па.

Длина хорды лопатки м.

Максимальная стрела прогиба у корневого сечения м.

Максимальная толщина профиля корневого сечениям.

Окружные составляющие скорости газа в относительном движении

м/с;

м/с.

Статическое давление на входе и выходе РК:

Па

Па

Статическая температура на входе и выходе РК:

К

К

Плотность газа на входе и выходе РК:

кг/м3

Расчет интенсивности газовых сил:

в плоскости вращения на среднем радиусе:

в осевой плоскости в корневом сечении:

в осевой плоскости в периферийном сечении:

Исходные данные для расчета на ЭВМ

РАСЧЕТ НА ПРОЧНОСТЬ ПЕРА РАБОЧЕЙ ЛОПАТКИ КОМПРЕССОРА (ТУРБИНЫ)

-------------------------------------------------------------------------------

ВЫПОЛНИЛ(А): koshkina

УЗЕЛ ДВИГАТЕЛЯ: компрессор МАТЕРИАЛ: vt-3

ИСХОДНЫЕ ДАННЫЕ:= 1.000000 CL= 6.200000E-02 RK= 9.300000E-02 RP= 1.550000E-01

VP= 0.000000E+00 UPP= 0.000000E+00 APP= 0.000000E+00= 19718.200000 AA= 0.000000E+00 AU= 0.000000E+00 PU= 1826.270000= 1559.310000 PAP= 2428.830000 RO= 4500.000000= 2.800000E-02 2.800000E-02 2.800000E-02= 7.000000E-02 5.500000E-02 4.000000E-02= 1.960000E-03 1.540000E-03 1.120000E-03= 1.380000 1.180000 7.800000E-01

SPT= 1050.000000 1050.000000 1050.000000 1050.000000

.000000 1050.000000 1050.000000 1050.000000

.000000 1050.000000 1050.000000- указатель типа расчитываемой лопаточной машины (компрессор);, PO9 - плотность материала лопатки и ее отклонение (кг/м³);- предел длительной прочности материала пера лопатки - массив из 11 чисел (МПа); VSPT - коэфициент вариации предела длительной прочности;, OMEGA9 - угловая скорость и ее отклонение (рад/с);, N9 - мощность ступени и ее отклонение (Вт);, MG9 - расход газа и его отклонение (кг/с);A, C2A, CA9 - осевые составляющие абсолютной скорости на входе и на выходе из ступени и их отклонение (м/с);, P2, P9 - давление на входе и на выходе из ступени и его отклонение (Н/м²);- число лопаток;, L9 - длина пера лопатки и ее отклонение (м);, R2, RCP - радиус корневого и периферийного сечений, их отклонение (м);- объем бандажной полки (м³);- вынос центра тяжести полки в окружном направлении (м);- вынос центра тяжести полки в осевом направлении (м);- относительный вынос центра тяжести периферийного сечения в осевом направлении;- относительный вынос центра тяжести периферийного сечения в окружном направлении;, B9 - хорда лопатки и ее отклонение (м);, D9 - максимальная толщина профиля и ее отклонение (м);, H9 - максимальная величина прогибы средней линии профиля (массив из трех чисел) и ее отклонение (м);- угол установки профиля (массив из трех чисел) (рад);- интенсивность окружных газовых сил (Н/м);- интенсивность осевых газовых сил в корневом сечении (Н/м);- интенсивность осевых газовых сил в периферийном сечении (Н/м);

Результаты расчета лопатки на статическую прочность на ПЭВМ

Результаты расчета на прочноcть пера лопаткиX F Jmin Spakt SизгA SизгB SизгCm^2 m^4 МПа МПа МПа МПа

1 .00000 .136E-02 .394E-06 113.265 -.682 .947 -.374

.00620 .130E-02 .324E-06 106.652 -.574 .800 -.343

.01240 .124E-02 .283E-06 99.199 -.462 .647 -.305

.01860 .118E-02 .249E-06 90.857 -.354 .499 -.264

.02480 .113E-02 .219E-06 81.562 -.256 .364 -.221

.03100 .107E-02 .191E-06 71.240 -.172 .247 -.176

.03720 .101E-02 .165E-06 59.798 -.103 .151 -.130

.04340 .951E-03 .141E-06 47.120 -.052 .078 -.085

.04960 .893E-03 .117E-06 33.061 -.020 .031 -.044

.05580 .834E-03 .950E-07 17.435 -.004 .006 -.013

.06200 .776E-03 .735E-07 .000 .000 .000 .000SсумA SсумB SсумC Ka Kb Kc

[МПa] [МПa] [МПa]

112.584 114.212********* 9.326 9.193 9.301

106.078 107.452 106.309 9.898 9.772 9.877

98.738 99.846 98.895 10.634 10.516 10.617

90.502 91.356 90.593 11.602 11.493 11.590

81.306 81.926 81.341 12.914 12.816 12.909

71.068 71.487 71.064 14.775 14.688 14.775

59.695 59.948 59.667 17.590 17.515 17.598

47.068 47.198 47.035 22.308 22.247 22.324

33.041 33.092 33.017 31.778 31.730 31.802

17.431 17.441 17.422 60.236 60.201 60.268

.000 .000 .000************************

Расчет потребного коэффициента запаса прочности производим по методике [1].

Определение среднеквадратических отклонений возмущающих факторов:

Определение номинальных значений коэффициентов влияния на напряжение растяжения:

) Плотности материала лопатки:

)Угловой скорости:

) Высота пера лопатки:

) Среднего радиуса:

Определение дисперсии напряжения растяжения от центробежных сил:

Определение номинальных значений коэффициентов влияния составляющих изгибающих моментов по главным осям инерции на напряжение изгиба от газовых сил:

Где: м

-максимальная толщина профиля во втулочном сечении;=h=0,00127м -максимальная стрела прогиба во втулочном сечении;

Где м;

Где:

Определение номинальных значений коэффициентов влияния на составляющие изгибающих моментов по главным осям инерции:

Коэффициенты влияния на составляющую изгибающего момента:

) Мощности ступени:

) Высоты лопатки:

)Угловой скорости:

)Среднего радиуса:

)Условного радиуса:

) Величины :

)Расхода газа:

)Величины:

)Величины:

Коэффициенты влияния на составляющую изгибающего момента :

) Мощности ступени:

)Высоты лопатки:

)Угловой скорости:

)Среднего радиуса:

)Условного радиуса:

)Величины

)Расход газа:

)Величины:

)Величины:

Определение номинальных значений коэффициентов влияния на главные моменты инерции:

Определение дисперсии составляющих изгибающих моментов по главным осям инерции для газовых сил:

Определение дисперсии моментов инерции сечения:

Определение дисперсии напряжения изгиба от газовых сил:

Определение дисперсии полного напряжения (напряжения растяжения от центробежной силы и напряжения изгиба от газовых сил)

Определение коэффициента вариации напряжения:

Рисунок 1 - Графики функций плотностей распределения напряжений

Силовая установка самолета разбивается на следующие функциональные группы: двигатель, воздухозаборник, топливную систему самолета, систему управления двигателем, реверсивные и шумоглушащие устройства, противопожарную систему.

В свою очередь двигатель состоит из таких основных узлов и систем: компрессора, камеры сгорания, турбины, выхлопной системы с реактивным соплом, топливной системы, системы смазки, системы запуска, трансмиссии, приводов к агрегатам.

Сделав предположение об одинаковой надежности узлов и систем двигателя, находим потребную вероятность безотказной работы двигателя в течении полета.

) Требуемый уровень надежности двухдвигательной силовой установки при продолжительности полета

) Требуемый уровень надежности компрессора:

)Требуемый уровень надежности компрессора высокого давления:

)Требуемый уровень надежности ступени:

) Требуемый уровень надежности РК:

) Требуемый уровень надежности лопаточного венца:

) Требуемый уровень надежности лопатки:

Из таблицы значения функции нормированного нормального распределения определяем Х: Х=5,9

-требуемый коэффициент запаса прочности лопатки;

Таблица 1 - Вероятность неразрушения детали

Рисунок 2 - Зависимость вероятности неразрушения лопатки от коэффициента запаса прочности

Вывод: наиболее напряженным оказался участок в точке В в корневом сечении (согласно таблице 1), где суммарное напряжение равно . Определен требуемый запас прочности  для обеспечения требуемой вероятности неразрушения:

3. РАСЧЕТ НАДЕЖНОСТИ ДЕТАЛЕЙ ПРИ ПОВТОРНО-СТАТИЧЕСКОМ НАГРУЖЕНИИ

Работа узлов двигателя имеет циклический характер, а значит, и напряженно-деформируемое состояние двигателя изменяется циклически. Несмотря на то, что в каждом цикле нагружения носят статический характер, при повторных нагружениях в материале возникают явления, типичные для усталости. Поэтому разрушение деталей при сравнительно небольшом числе циклов (N=10²…10³) называют малоцикловой усталостью, а способность материала сопротивляться такому разрушению - малоцикловой прочностью.

Расчет проводим согласно рекомендациям, представленным в методических пособиях [1,2].

Количество циклов за ресурс работы двигателя

Где -ресурс работы двигателя (ч);

-время полета(ч);

Количество приемистостей за ресурс складывается из:

предполетной проверки

пробы перед взлетом  

взлета

ухода на второй круг

проверки после регламентных работ

количества прерванных взлетов

В эксплуатации за ресурс работы максимальное количество взлетных режимов:

 минимальное

Математическое ожидание взлетных режимов:

Среднеквадратическое отклонение:

Определение параметров кривой усталости по данным [2] для температуры лопатки t=223,35оС:

Определение среднего напряжения  и амплитуды  пульсирующего цикла

Определение величины напряжений

Проанализировав особенности конструкции, технологии изготовления и эксплуатации, определяем коэффициенты кσ, εσ, βσ ,ψσ.

Кσ - эффективный коэффициент концентрации напряжений:

где

 - коэффициент чувствительности материала к концентрации напряжений (для сплавов титана ).

Принимаем

 -теоретический коэффициент концентрации напряжений, равный отношению напряжения при наличии концентратора и напряжению в той же точке в отсутствии концентратора, принимаем  = 1,2.

Тогда

 - коэффициент учитывающий абсолютные размеры детали (масштабный фактор),

где ε∞ =0,5-для деформируемых материалов,

λ=0,02 (1/мм),характерный размер детали, для лопатки - это хорда: d=b=18 (мм);

βσ - коэффициент, учитывающий влияние состояния поверхности и упрочнения,

βσ= βт∙βкор∙βуп, где

βт - отображает влияние шероховатости.

Для шлифования определяем коэффициент βт по графику2 (рисунок 2.1):

βт =0,9;

βкор- коэффициент характеризующий коррозионное повреждение поверхностного слоя.

Определяем βкор для образца с концентраторами напряжений в пресной воде по графику 2 (рисунок 2.2): βкор=0,35;

βуп - коэффициент характеризующий выносливость деталей при упрочняющей технологии или без ее применения ,

Принимаем по данным [2] для специальной термической обработки: βуп=2,5.

Рисунок 3 - Зависимость βт от предела прочности материала

Рисунок 4 -Зависимость βкор от предела прочности материала

Тогда коэффициент состояния поверхности

Коэффициент  характеризует чувствительность материала к асимметрии цикла. По данным [2] для сплавов титана: Принимаем:

Так как полученное значение то производим расчет математического ожидания максимального числа циклов до разрушения лопатки по следующей формуле:

Построим кривую усталости в логарифмических координатах для графического определения числа циклов нагружения до разрушения лопатки (Рисунок 5):

Рисунок 5 - Графическое определение

Находим среднеквадратическое отклонение числа циклов до разрушения:

Определяем среднеквадратическое отклонение и математическое ожидание величины накопленных повреждений за ресурс работы двигателя:

Где

Вычисляем вероятность неразрушения детали с учетом величины накопленных повреждений за ресурс работы двигателя:

В расчете лопатки при действии повторно-статических нагружений была получена вероятность безотказной работы  (при ресурсе двигателя 25000 часов). Таким образом, при повторно-статических нагружениях работоспособность лопатки сохраняется гарантированно, однако единичное значение вероятности говорит об излишнем запасе прочности лопатки и чрезмерном расходе материала.

4. РАСЧЕТ НАДЕЖНОСТИ С УЧЕТОМ ВНЕЗАПНЫХ ОТКАЗОВ ДИСКА КОМПРЕССОРА

Диски компрессоров и турбин - это наиболее ответственные элементы конструкции газотурбинных двигателей. От совершенства конструкции дисков зависит надёжность, лёгкость конструкций авиационных двигателей в целом.

Диски находятся под воздействием инерционных центробежных сил, возникающих при вращении от массы рабочих лопаток и собственной массы дисков. Эти силы вызывают в дисках растягивающие напряжения. Кроме напряжений растяжения и сжатия в дисках могут возникать напряжения кручения и изгиба.

Цель расчёта на прочность диска - определение напряжений и запасов прочности в различных сечениях по радиусу диска.

Расчёт производим на ПЭВМ с помощью программы Disk.exe.

При помощи кафедральной программы disk_nad.exe и методических пособий ([1], [4]) проводим расчет диска в 12 сечениях, в результате которого получаем распределение прочностных характеристик и вероятностей неразрушения диска за 1 секунду по сечениям.

Исходные данные для расчета на ЭВМ:- число расчетных сечений- количество скачков на профиле диска- признак постоянства плотности материала диска по его радиусу- признак постоянства температур диска по радиусу- частота вращения диска- радиальное напряжение на внешнем контуре диска, которое эквивалентно распределенной нагрузке от лопаток и замковой части диска, BZ, NZ, QZ - коэффициенты основных формул для нулевого сечения- коэффициент Пуассона(1..K) - радиусы расчетных сечений(1..K) - толщины диска в расчетных сечениях(1..Z) - номера радиусов, на которых профиль диска имеет скачки(1..K) - значения приведенной плотности материала диска при расчете радиальной лопастной машины- плотность материала диска при расчете осевой лопастной машины(1..K) - температура в расчетных сечениях диска(1..K) - модуль упругости(1..K) - коэффициены линейного расширения(1..K) - пределы длительной прочности- предел длительной прочности материала диска

STRENGTH ANALYSIS OF COMPRESSOR AND TURBINE DISCS

*************************************************************: TERNUK IGOR 242 DATA:= 0 DT= 0speed = 14051.2 rpmof design cross-sections = 11of steps in disc profile = 1load = 16.0 MPa=0.0 BZ=0.0 NZ=1.0 QZ=0.0factor = 0.3(1)= 0.1250 R(2)= 0.1340 R(3)= 0.1400 R(4)= 0.1440(5)= 0.1465 R(6)= 0.1480 R(7)= 0.1520 R(8)= 0.1750(9)= 0.2000 R(10)= 0.2000 R(11)= 0.2117 R((1)= 0.0103 B(2)= 0.0103 B(3)= 0.0103 B(4)= 0.0083(5)= 0.0067 B(6)= 0.0054 B(7)= 0.0050 B(8)= 0.0050(9)= 0.0050 B(10)= 0.0200 B(11)= 0.0200 B(factor = 4530.0(Z) 9

Результаты расчетаR(I),M B(I),M SR,MPa ST,MPa SEK,MPa ZAP

0.1250 0.0103 0.0000 519.2178 519.2178 1.8297

0.1340 0.0103 25.9524 478.4047 465.9709 2.0388

0.1400 0.0103 38.1492 455.7270 437.9005 2.1694

0.1440 0.0083 51.9166 444.3472 420.7978 2.2576

0.1465 0.0067 65.1761 439.4681 410.7765 2.3127

0.1480 0.0054 79.4881 438.7797 404.9299 2.3461

0.1520 0.0050 89.2019 429.0701 392.1534 2.4225

0.1750 0.0050 103.7460 366.5777 327.2788 2.9027

0.2000 0.0200 23.8295 293.7911 282.6308 3.3613

0.2117 0.0200 16.0000 270.9117 263.2765 3.6084,MPa DSR,MPa ST,MPa DST,MPa SEK,MPa DSEK,MPa

0.0000 0.0000 519.2178 1214.5548 519.2178 1214.5548

.9524 6.2981 478.4047 1046.4294 465.9709 1052.7275

.1492 14.7856 455.7270 961.5753 437.9005 976.3610

.9166 28.8129 444.3472 927.5266 420.7978 956.3395

.1761 46.7588 439.4681 919.8577 410.7765 966.6165

.4881 70.6508 438.7797 929.4096 404.9299 1000.0604

.2019 92.6420 429.0701 904.6846 392.1534 997.3266

.7460 161.7937 366.5777 732.9059 327.2788 894.6996

.3180 217.7879 315.2376 629.8324 280.0223 847.6203

.8295 13.6113 293.7911 467.9775 282.6308 481.5888

.0000 44.7790 270.9117 425.3217 263.2765 470.1007

Вычисляем коэффициент вариации напряжения:

Вероятность неразрушения детали:

где  - коэффициент запаса прочности;

 - коэффициент вариации предела длительной прочности;

После расчета переменной х, с помощью переводных таблиц [2] найдем значение вероятности неразрушения при принятом коэффициенте запаса прочности.

Определим вероятность неразрушения конструкции:

ЗАКЛЮЧЕНИЕ

компрессор лопатка прочность напряжение

В ходе выполнения данной работы был проведен расчет вероятности безотказной работы лопатки и диска рабочего колеса входной ступени дозвукового осевого компрессора.

Расчет надежности лопатки с учетом внезапных отказов показал, что наиболее напряженным оказался участок в точке В в корневом сечении, где суммарное напряжение равно . Путем расчета был получен требуемый запас прочности , удовлетворяющий нормам прочности, предъявляемым к современным авиационным ГТД.

В расчете лопатки при повторно-статических нагружениях была получена вероятность безотказной работы  (при ресурсе двигателя 25000 часов). Таким образом, при повторно-статических нагружениях за ресурс двигателя должна сохраниться работоспособность лопатки.

Вероятность безотказной работы всей конструкции составила .

ПЕРЕЧЕНЬ ССЫЛОК

1.    Москаленко А.С. «Расчет надежности деталей авиационных газотурбинных двигателей» Харьков «ХАИ» 1985 г.

2.      Москаленко А.С. «Расчет надежности авиационного газотурбинного двигателя» Харьков «ХАИ» 1990 г.

.        Шошин Ю.С., Епифанов С.В., Зеленский Р.Л. «Расчет на прочность рабочих лопаток компрессоров и турбин» Харьков «ХАИ» 2006 г.

.        Шошин Ю.С., Епифанов С.В., Зеленский Р.Л. «Расчет на прочность дисков компрессоров и турбин» Харьков «ХАИ» 2007 г.