Определение факторов, оказывающих наибольшее влияние на показатели надёжности диска ГТД с учетом внезапных отказов Текст научной статьи по специальности «Механика и машиностроение»

УДК 621.452.3

ОПРЕДЕЛЕНИЕ ФАКТОРОВ, ОКАЗЫВАЮЩИХ НАИБОЛЬШЕЕ ВЛИЯНИЕ НА ПОКАЗАТЕЛИ НАДЁЖНОСТИ ДИСКА ГТД С УЧЕТОМ ВНЕЗАПНЫХ ОТКАЗОВ

© 2009 А. И. Белоусов, А. В. Грицин Самарский государственный аэрокосмический университет

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

Надежность, интенсивность отказов, параметры надежности, диск турбины, метод конечных элементов

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

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

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

Требуемая надёжность двигателей закладывается на этапе проектирования, доводится до требуемого для начала эксплуатации уровня на этапе опытно-

конструкторской отработки, совершенствуется и обеспечивается при производстве, улучшается и реализуется при эксплуата-ции[1].

Но даже если конструирование, расчёты, изготовление, сборка, испытания отдель-

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

Одной из основных задач проблемы обеспечения безотказности и долговечности является разработка методов определения числовых значений показателей надёжности изделий. Эта задача решается методами обработки данных эксплуатационных наблюдений и ускоренных испытаний, а также расчётным путем. Расчётный метод в настоящее время слабо оформился в инженерный метод из-за сложности задачи, ее новизны и недостаточного количества фактических и опытных данных. Но в настоящее время они наиболее экономически целесообразны и дают возможность судить о надёжности изделий в нормальных условиях эксплуатации по значению соответствующих показателей при форсированных режимах (повышенные нагрузки, скорости, температуры, концентрации реагентов и т.д.).

По характеру потери работоспособности отказы двигателей делятся на две группы. К первой группе относятся внезапные отказы, носящие характер случайного выброса: хрупкое разрушение; превышение предела текучести в какой-либо точке детали, для которой остаточные деформации недопустимы; возникновение слишком больших упругих деформаций и др. Ко второй

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

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

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

К особенностям ДЛА с точки зрения оценки показателей надёжности относятся [2]:

• сравнительно небольшой объём статистической совокупности (в ряде случаев единичные экземпляры);

• трудность получения опытных данных для последующей обработки методом математической статистики;

• сложный характер и взаимосвязанность процессов нагружения;

• повышенное влияние на надёжность конструкций стабильности технологических процессов и сборочных работ. Нагрузки, действующие на конструкцию, как правило, оказывают многократное влияние или развертывают свои вероятностные свойства во времени. Механические свойства материалов, изготовляемых в массовом количестве, могут быть изучены основательно [3]. Таким образом, поведение самой уникальной конструкции определяется случайными факторами массового характера, для каждого из которых допускается статистическое толкование вероятности. Предсказать на основе этого статистического массива данных поведение конструкции - в этом и состоит цель расчета надёжности ДЛА.

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

ствующих нагрузок и фактической прочности. Разрушения из-за нарушения прочности могут возникнуть в результате непредвиденных местных напряжений, вызванных нерасчётными статическими или динамическими нагрузками в сочетании с невыявленными производственными отклонениями. Указанные местные напряжения и нагрузки, как случайные величины, существенно отличаются от средних величин, принимаемых при прочностных расчётах и оценке ресурса изделия.

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

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

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

Введём понятие функции качества (запаса работоспособности) [1]

V = 0пред ~0р , (1)

где о пред - предельные напряжения (предел длительной прочности, предел выносливости и др.); о - рабочее напряжение в расчетном сечении.

Для надёжной работы детали двигателя необходимо выполнение условия

V > 0 . (2)

Предельные напряжения и рабочее напряжение - случайные величины. Они имеют нормальное или логарифмически-нормальное распределение вокруг среднего значения с некоторым стандартным отклонением от него [1,2,4]. Следовательно, функция качества р согласно выражению (1), как композиция распределений о пред и о , также является случайной величиной.

Применение нормального закона оправдано тем, что в случае совместного действия достаточно большого числа случайных возмущений, подчиняющихся различным законам распределения (если среди них нет превалирующего) результирующее возмущающее воздействие имеет согласно центральной предельной теореме теории вероятностей распределение, близкое к нормальному. Тогда согласно свойствам композиции законов распределений и функция качества р описывается нормальным законом Гаусса с математическим ожиданием р = о ед - о и среднеквадратическим

отклонением ор =л Юо 2 +оо 2, где о

пред и о р

- математически ожидания, а

оо - среднеквадратические отклонения предельного о пред и рабочего о р напряжений

соответственно.

Вероятность неразрушения детали согласно условию (2)

)>-

Рв{|) = Вер[(о„ред -оР)> °] =

¥ 1 ¥ {р

|/р)р =---------г= |е 2

° р2т °

йр, (3)

¥

Рв {) = Т2Т У

dz =

1 ¥ - 2 1 у - 2

1 е 2 1 е 2 = 0.5 + ¥ {у), (5)

О'Тг » - Опт "

л/2р

л/2Г

где у = = — - коэффициент однородно-

ор

ор

сти (гауссовская мера надёжности); ир =^т

V

- коэффициент вариации функции качества р.

Для интеграла вероятности (4) имеются таблицы в справочниках по математике.

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

и о

р•

пред

-1

У =

о пред ор

2 2 оо + оо

°пред °р

оо

оо

— 2 — ■

о р о пред

— 2 ор

и

к -1

к2 2 + ио2

° пред ° р

(6)

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

чета на прочность) к = —пр

о

пред

?

р

предельных

коэффициент

напряжений

вариации

= ; коэффициент вариации рабо-

о пред

где I - наработка, для которой определяется вероятность неразрушения.

В справочнике [5] приведены значения опред и о для некоторых авиационных материалов.

Введем переменную 2 = р—р . Тогда

ор

йр = орй2 . Выражение (3) для вероятности

неразрушения преобразуется в нормированный нормальный закон распределения

е 2 ; (4)

чего напряжения и = ——.

ор о

р

Вероятность разрушения (внезапного прочностного отказа)

0{1 )= 1-Р. {I). (7)

Показателем для оценки надёжности (безотказности) детали двигателя является интенсивность отказов 1, которая для наиболее опасных нелокализованных разрушений должна быть меньше 1-10-91/ч. Для определения интенсивности отказов необходимо знать время работы двигателя на режимах и вероятность разрушения на каждом режиме.

р

2

Напряжения в детали от центробежных сил на режимах работы двигателя приближённо можно оценить по формуле

(?р І =(?р ),

где

к)

р ^max w W

max

(8)

р,-э и &р.э = ар • ®max

напряжение в

детали и угловая скорость на рассматриваемом режиме эксплуатации; (s р )max и wmax -напряжение и угловая скорость на максимальном (взлетном) режиме; ар - коэффициент режима.

В первом приближении можно считать, что статистические разбросы прочности vs д и рабочего напряжения vs не зависят

от режимов работы двигателя (на всех режимах эксплуатации vs = const и

° пред

vs = const), т.е. функция вероятностей F(x)

зависит только от коэффициента запаса прочности.

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

(ар )

(ар )р.э (ар )max (ар )р

пред

к.

а

р р. э р max р р. э р

где квзл - коэффициент запаса на взлётном режиме.

Определив по зависимости (6) коэффициент однородности g, по таблицам соответствующий интеграл вероятности F (g), по формулам (5) вероятность неразрушения детали и (7) вероятность разрушения на каждом режиме работы двигателя, находим интенсивность отказов детали за 1 ч работы двигателя:

(nv Л

I Q(t)

• t

р. э. р. э

/1 [ч-1],

(9)

где п - количество режимов эксплуатации двигателя; Q(V) рэ - вероятность отказа (разрушения) детали на данном режиме работы V

7 рэ

двигателя; V =— - относительная нара-

t

ботка на рассматриваемом режиме; їрэ - наработка на рассматриваемом режиме; ї -сумммарная наработка изделия в эксплуатации.

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

Р

ли напряжений.

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

можно найти несколькими способами: методом малых возмущений, с помощью модуля “Probabilistic Design System” (PDS) в пакете конечно-элементного анализа ANSYS и д.р. Результатом вероятностного расчета в пакете ANSYS являются интегральные функции распределения возмущающих и результирующих факторов, их номинальные значения и среднеквадратические отклонения, вероятности достижения ими заданных значений, матрицы корреляций и графики вероятностных коэффициентов чувствительности, определяющие рассеивание «выходного» параметра при наличии рассеивание «входных» параметров.

Сравнение двух методов и пример их использования приводятся в [6]. В настоящей работе будем использовать программный пакет ANSYS, а в качестве примера, рассмотри рабочее колесо турбины ТРДД НК-8-4 самолета Ил-62.

На надёжность работы диска ГТД влияют разнообразные факторы. Но среди них можно выделить наиболее очевидные, которые могут оказывать наибольшее воздействие на надёжность работы. При производстве эти параметры невозможно выдержать точными, поэтому на них назначаются допуски с учётом достижений современной технологии, затрат на изготовление и стоимости окончательного изделия, условий сборки, оптимальной работы и т.д. Если бы не было допусков, разброса параметров и условий работы, то изделие работало бы безотказно. Но в действительности допуски на параметры могут сложиться самым неблагоприятным образом, что приводит к потере работоспособности изделия.

Проведя анализ работы диска турбины двигателя НК-8-4, мы выделили следующие параметры, которые наиболее очевидно оказывают влияние на работоспособность диска и являются удобными для замеров:

- геометрические размеры;

- частота вращения;

2

- плотность материала диска.

Рассмотрим влияние допусков на размеры, частоту вращения и плотность материала диска на его показатели надёжности. Для того чтобы иметь возможность проводить вероятностные расчёты, исследовать влияние возмущающих факторов на интересующий нас параметр, сначала необходимо создать параметрическую модель рассчитываемой детали и рассчитать её на прочность. Модель сектора рабочего колеса и графическое отображение результата расчёта на прочность представлены на рис. 1.

Рис. 1. Конечно-элементная модель диска и распределение эквивалентных напряжений

Проведем расчёт показателей надёжности для варианта с установленными допусками. В результате расчёта в программе АЫ-БУБ был получен интегральный нормальный закон распределения рабочего напряжения, среднеквадратическое отклонение которого равно оа = 28.7 МПа.

По формулам, приведенным в [6], найдём интенсивность отказов двигателя

1 = 2.378 -10-10 1/ч. Полученное значение интенсивности отказов меньше величины

1-10"9 1/ч, следовательно, диск отвечает мировым требованиям по надёжности.

Для лучшей обработки результатов и упрощения расчетов были приняты следующие допущения:

1) основные размеры замка и лопатки остаются неизменными;

2) при изменении допуска на один параметр допуски на другие параметры остаются неизменными.

Сначала проанализируем влияние допусков геометрических размеров. Выберем в качестве проверочных следующие размеры диска (рис.2):

• наружный диаметр диска DHap;

• наружный D1 и внутренний D2 диаметры места под болты;

• ширина места под болты H;

• диаметр ступицы Dcm.

I ■ '

■( I

1.. /

I

^ I ^ I

Рис. 2. Геометрические размеры, выбранные для анализа

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

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

Теперь рассмотрим влияние частоты вращения ротора высокого давления на показатели надёжности диска. Частота вращения ротора ВД двигателя НК-8-4 на крейсерском режиме составляет 6850-8° об/мин.

Результаты расчёта представлены графически на рис.4.

Рис. 3. Изменение интенсивности отказов при изменении допусков 5г на размеры

Рис. 4. Изменение интенсивности отказов при изменении допуска д„ частоты вращения

Аналогичным образом рассмотрим влияние плотности материала, на который также существует допуск. Результаты расчёта представлены на рис. 5.

Ы010,чл

Рис. 5. Изменение интенсивности отказов при изменении допуска др на плотность материала

Как видно из результатов расчёта, влияние допусков всех параметров на пока-

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

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

4, а 100 или 200, то все вместе они оказали бы колоссальное воздействие на работоспособность изделия.

Библиографический список

1. Белоусов, А.И., Биргер И. А. Прочностная надежность деталей турбомашин: Учеб. Пособие / А.И. Белоусов, И. А. Биргер. - Куйбышев: Ку-АИ, 1983. - 75с.

2. Косточкин, В. В. Надежность авиационных двигателей и силовых установок: Учебник для студентов / В.В. Косточкин. - М: Машиностроение, 1976. - 248с.

3. Акимов, В.М. Основы надежности газотурбинных двигателей: Учебник для студентов машиностроительных специальностей вузов / В.М. Акимов. - М.: Машиностроение, 1981.207с.

4. Скубачевский, Г.С. Авиационные газотурбинные двигатели. Конструкция и расчет деталей: Учебник для студентов авиационных вузов / Г.С. Скубачевский.- М.: Машиностроение, 1974.

- 520 с.

5. Вероятностные характеристики прочности авиационных материалов и размеров сортамента: Справочник / Под ред. С.О. Охапкина. - М.: Машиностроение, 1970. -568 с.

6. Белоусов, А.И. Оценка параметров прочностной надежности деталей ДЛА на этапе проектирования: Учебное пособие. / А.И. Белоусов, А.В. Грицин. - Самара: СГАУ, 2006 - 115с.

References

1. Belousov A.I., Birger I.A. Strength reliability of turbomachine parts: Tutorial. - Kuybyshev: KuAI, 1983. - 75 pages.

2. Kostochkin V.V. Reliability of aircraft engines and power plants: Tutorial for students. -Moscow: Mashinostroyeniye (Machine building), 1976.-248 pages.

3. Akimov E.M. Fundamentals of reliability of gas turbine engines: Tutorial for student of machine building specialty. - Moscow: Mashi-

ESTIMATION THE FACTORS, MAKE THE GREATEST INFLUENCE ON RELIABILITY PARAMETRES OF A GTD DISK WITH A GLANCE SUDDEN FAILURES

© 2009 A. I. Belousov, A. V. Gritsin Samara State Aerospace University

In this work the methodology of an estimation of parameters of reliability GTD with a glance sudden failures is presented. By this methodology the analysis of influence of the limit of the geometrical sizes, rotating frequency and disk material density on failure rate is carried out.

Reliability, failure rate, reliability parameter, turbine disk, finite element method

Информация об авторах

Белоусов Анатолий Иванович, доктор технических наук, профессор кафедры Конструкции и проектирования двигателей летательных аппаратов Самарского государственного аэрокосмического университета. E-mail: aibelousov@mail.ru. Область научных интересов: динамика и прочность, надежность авиационных и ракетных двигателей, конструкция ракетных двигателей.

Грицин Алексей Валерьевич, инженер-конструктор, аспирант Самарского государственного аэрокосмического университета. E-mail: grialeksey@gmail.com. Область научных интересов: надежность двигателей летательных аппаратов.

Belousov Anatoliy Ivanovich, doctor of technical sciences, professor of department of the Construction and designing of the aircraft engines of Samara State Aerospace University. E-mail: aibelousov@mail.ru. Area of research: dynamics and durability, reliability of the aircraft and rocket engines, construction of the rocket engines.

Gritsin Aleksey Valerevich, engineer-designe, postgraduate student of Samara State Aerospace University. E-mail: grialeksey@gmail.com. Area of research: reliability of the aircraft and rocket engines.

nostroyeniye (Machine building), 1981.-207 pages.

4. Skubachevskiy G.S. Aircraft gas turbine engines. Construction and calculation parts: Tutorial for students of aircraft institute of higher education

- Moscow: Mashinostroyeniye (Machine building), 1974. - 520 pages.

5. Probabilistic characteristics of aircraft mate-

rials strength and range sizes: Reference book. Edited by S.O. Ohapkina - Moscow: Mashino-

stroyeniye (Machine building), 1970. - 568 pages.

6. Belousov A.I., Gritsin A.V. Estimation strength reliability parameters of aircraft engine units on the design stage: Tutorial. - Sama-ra:SSAU, 2006 - 115 pages.