Обоснование долговечности конструкции транспортного самолёта при многоцелевом применении по результатам ресурсных и лётных испытаний Текст научной статьи по специальности «Механика и машиностроение»

УДК 629.7.02.015.4

ОБОСНОВАНИЕ ДОЛГОВЕЧНОСТИ КОНСТРУКЦИИ ТРАНСПОРТНОГО САМОЛЁТА ПРИ МНОГОЦЕЛЕВОМ ПРИМЕНЕНИИ ПО РЕЗУЛЬТАТАМ РЕСУРСНЫХ И ЛЁТНЫХ ИСПЫТАНИЙ

© 2013 П. А. Фомичёв1, В. И. Клепцов2

Национальный аэрокосмический университет им. Н.Е. Жуковского «ХАИ» 2ОАО «Таганрогский авиационный научно-технический комплекс им. Г.М. Бериева»

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

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

Введение

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

В мировой практике существует несколько основных подходов к проектированию авиационных конструкций. Изначально самолёты создавали по методу обеспечения безопасного ресурса (safelife). Этот метод гарантирует конкретный период эксплуатации для всего парка самолётов данного типа при принятых одинаковых условиях без серьёзной вероятности катастрофического или опасного усталостного разрушения. Многолетний опыт создания самолётов показал, что при необходимости увеличения ресурса целесообразно применение метода эксплуатационной живучести, который включает в себя принципы допустимости повреждения (damage tolerance) и безопасного разрушения (fail-safe). Данный метод основан на предположении, что после повреждения либо разрушения одного или нескольких элементов оставшаяся неповрежденной часть конструкции должна выдерживать эксплуатационные нагрузки. В связи

с этим метод эксплуатационной живучести дополняют требованием периодических осмотров в течение всего срока эксплуатации каждого самолёта.

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

В работах [1-5] рассмотрены проблемы проектирования и установления ресурса эксплуатируемых самолётов.

Заключение о ресурсе конструкции получают по результатам лабораторных ресурсных испытаний (РИ). Программу РИ разрабатывают на основе предполагаемых профилей типовых полётов с учётом накопленной ранее статистики для самолётов аналогичной категории. Нагрузки на воздушных участках типовых полётов находят в соответствии с отраслевым стандартом турбулентности атмосферы [6]. В процессе эксплуатации профили типовых полётов уточняют, возможна эксплуатация самолёта как транспортного средства по не предполагавшимся ранее профилям. Для самолёта Бе-200ЧС наиболее повреждающими являются участки

полётов на пожаротушение, для которых характерным является полёт на высотах до 400 м в зоне повышенной турбулентности. Однако отраслевой стандарт не устанавливает характеристики турбулентности на высотах до 300 м [6].

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

Соотношения для расчёта долговечности при многоцелевом применении

Пусть в процессе ресурсных испытаний долговечность до возникновения трещин в регулярных зонах составила 1 исп блоков нагружения (типовых полётов). Усталостное повреждение за один блок программы РИ обозначим Б исп, тогда общее усталостное повреждение в испытаниях составит

Б = Б

ри и

•1 и

(1)

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

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

Усталостное повреждение за 1 полётов в эксплуатации равно

Б э

= Е Б, -11

где 1, - число типовых полётов на 1-м профиле;

1 , = 1-а,,

Б э

= 1-Е Б, -а, ,

(2)

Е а, = 1.

Условие отсутствия разрушения в эксплуатации запишем в виде

Б < Б .

эксп ри

Это соотношение с учётом (1) и (2) примет вид

к Б. 1-Е—— -а, < 1,

; Б ' 1

(3)

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

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

Т = 1 исп - т

исп исп

Т, = 1, - т,

где т исп, т, - время типового полёта, принятого в испытаниях, и время для ,-го профиля. Тогда

1 исп =

Ти

1, =

Т,

После подстановки этих соотношений в (3) имеем

кБ

Е—

,Б

- Т, < Ти

(4)

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

х

х

х

х

повой полёт каждого из рассматриваемых профилей.

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

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

Повреждение от случайных нагрузок

Отраслевой стандарт «Модель турбулентности атмосферы» ОСТ 1 02514-84 устанавливает характеристики турбулентности, влияющие на прочность и выносливость конструкции для всех стадий проектирования, испытаний и эксплуатации [6]. В стандарте регламентирована интегральная повторяемость дискретных вертикальных воздушных порывов, превышающих вертикальную скорость заданной величины на 1 км пролетаемого самолётом пути, в зависимости от высоты полёта

Дп = СУ Р кр - К. у "

где V. - индикаторная скорость полёта на высоте Н., р0 - плотность воздуха на уровне моря, О. - вес самолёта на .-м режиме полёта, К. - коэффициент ослабления порыва [6].

Каждый этап типового полёта, включающего набор высоты, крейсерский полёт, снижение и посадку, разобьём по времени на несколько режимов с постоянными для каждого режима скоростью, высотой, весами самолёта и топлива в крыле. Общее число режимов за весь типовой полёт обозначим через г. Определим высоту Н., скорость V., веса самолёта О. и топлива О, пролетаемый путь Ь. для

каждого .-го режима полёта.

С учётом пути Ь., пролетаемом на м режиме, число порывов, превышающих скорость составит

= ..е

Введём обозначение

Су р0 ^ ^р т--К..

2О. .

(5)

Тогда

= р - е

"е-

ДПу

Р.(ДПу) = ..е

где Б0, С№ - параметры интегральной повторяемости порывов, зависящие от высоты полёта; р - общее число порывов на 1 км пути; Б(-) - число порывов на 1 км пути с эффективной вертикальной скоростью, превышающей

Приращение вертикальной перегрузки Апу в центре масс самолёта зависит от вертикальной скорости порыва:

Вероятность превышения перегрузки Апу на .-м режиме

Р.(ДПу)

Р (ДП у )

Дпу

Плотность вероятности распределения приращений перегрузки

С

С

С

е

фДАпу;> =

¿РДЛпу) _ 1

АПу

алп.

сп

Вероятность попадания в интервал ёДпу составит

Ф ДЛПу)ёЛПу.

Поскольку общее число перегрузок на режиме равно LJ • F0J, то приращение

числа циклов нагрузок dn _ LjF0jj (Лпу)dЛny

или

L:F.

Лпу

ап _,,е с^аЛпу

сп: у'

Накопленное на ,-м режиме усталостное повреждение в соответствии с линейной гипотезой суммирования усталостных повреждений равно

Лпу

Б

ап _ ,, Г 1е -

N С - N

е п аЛпу .

этом моменте, определённое по результатам расчёта общего напряжённого состояния крыла методом конечных элементов; М, - изгибающий момент в сечении крыла при единичной перегрузке на ,-м режиме полёта; О ■ , От, - веса самолёта и

Ч Ц

топлива в крыле на ,-м режиме полёта; Осо , Ото - веса самолёта и топлива в крыле, при которых выполнен расчёт изгибающих моментов; а^), Ь^), 0(2) - параметры зависимости изгибающего момента от весов самолёта, топлива в консоли и веса крыла при единичной перегрузке для сечения крыла с координатой 2.

Степенное уравнение кривой усталости примем в виде

N •ат _ С,

где N - число циклов до разрушения при регулярном нагружении; а - максимальное напряжение отнулевого цикла нагрузок; с, т - параметры уравнения.

С учётом степенного уравнения кривой усталости получим

С

(2ЛЛу (1 + ЛЛу ))

а,

В соответствии с формулой Одинга эквивалентное напряжение отнулевого цикла нагрузок на ,-м режиме полёта

аэкв _ а, ^ 2Му(1 + Лпу)

где а , - напряжение при единичной перегрузке на ,-м режиме полёта,

Зависимость для расчёта накопленного повреждения на ,-м режиме полёта для дискретной схемы атмосферной турбулентности примет вид

Б _ ^

с • сп

_ т -Лпу

•ат | (2Лпу(1 + Лпу))2 е ^ аЛпу

(7)

_ о ГТ

о з _--Мз,

Mэ 3

— О о

Мз _ а{г) • -3 - Ъ{2) • О3

ОО

(6)

- Ф).

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

Усталостное повреждение от случайных нагрузок на всех режимах 1-го типового полёта

Б с

(8)

с

е

т

2

Расчёт повреждения от цикла ЗВЗ

В соответствии с рекомендациями ЦАГИ для определения максимального приращения перегрузки цикла ЗВЗ необходимо принимать

Рс(Дп » ах) = 0,694.

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

Р(ДПутах) = е-Рс( ДПу).

В качестве вероятности превышения пиковых значений предложено принимать Р(Дпутах) = 0,5 . Условие для нахождения

приращения перегрузки цикла ЗВЗ имеет вид

Эквивалентное напряжение цикла ЗВЗ найдём по формуле Одинга

_ /-.^.звз / звз .-г-звз \

экв Л/ тах V тах тт '

или

а звз = _ экв ср

- /путах (путах + 0.

Число циклов до разрушения при регулярном нагружении с этим напряжением равно

N =

ЗВЗ

С

/_звз \ П

(СТэкв )

Усталостное повреждение от цикла ЗВЗ за один типовой полёт

Б звз = 1Ш звз .

г ДП звз

Е ..ехр( --С^) = 0,694.

.=1 Сп.

(9)

Это уравнение легко решить численно относительно Дп утах .

Максимальная перегрузка цикла ЗВЗ

пзвз = 1 + Дпзвз

у тах у тах •

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

_ звз звз

а = а - п

тах ср у тах ■

где аср - среднее значение напряжения

при единичной перегрузке на режимах полёта.

В качестве минимального напряжения цикла ЗВЗ принимают значение, соответствующее наземным режимам эксплуатации или пропорциональное а ср :

а.

- * -а ср .

Суммарное усталостное повреждение за ,-й типовой полёт составит

Б, = Б звз + Б сл . (10)

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

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

^□ах 1 т1п 1

моментов и число циклов их повторения п1 в одном блоке программы РИ. Эквивалентный отнулевой момент равен

М01 =

Мтах1 (Мтах1 - МШп1 ) .

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

м экв = (£ п; • мтт

Напряжение в элементе конструкции, соответствующее этому моменту, равно

0 экв

о

М

- • Мэкв.

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

Б исп =

1 С

о

(11)

Соотношения (10), (11) позволяют решить задачу об обосновании долговечности регулярных зон крыла самолёта согласно (3) или (4) при многоцелевом применении по результатам ресурсных испытаний. Пределы применимости обусловлены принятой гипотезой о линейном суммировании усталостных повреждений, что исключает необходимость в использовании дополнительных (зачастую противоречивых) эмпирических коэффициентов. Преимущество состоит в исключении параметра С, который зависит и от технологии изготовления натурной конструкции. В случае, когда условия эксплуатации самолёта, например Бе-200ЧС при полётах на пожаротушение, не соответствуют ограничениям, указанным в [6], целесообразно применение следующего подхода.

Усталостное повреждение (7) на ]-м режиме типового полёта зависит от параметров Б0 и Сп интегральной повторяемости приращения перегрузок. Эти параметры могут быть найдены по результатам обработки лётных измерений прира-

щений перегрузок. При этом важно синхронизировать значения приращений перегрузок с весами самолёта и топлива в крыле. Напряжения при единичной перегрузке должны быть найдены согласно (6) для рассматриваемого режима типового полёта.

На рис. 1 представлены результаты расчёта приращений перегрузок в центре масс самолёта Бе-200ЧС при выполнении транспортных полётов. Для сопоставления представлены интегральные повторяемости приращений перегрузок за один полёт самолётов КШ-95ЬК и Ту-134. Данные по этим самолётам приведены в статье [7]. В ней отмечено, что повторяемость Ту-134 соответствует лётным измерениям, а ЯК1-95ЬЯ - расчётам по ОСТ 1 02514-84 [6].

На этом же рисунке пунктирной линией представлена интегральная повторяемость приращений перегрузок самолёта Бе-200ЧС при полётах на пожаротушение по типовому профилю, рассчитанная согласно ОСТ [6]. Сплошная линия соответствует данным лётных измерений эквивалентных моментов в четырёх сечениях крыла, по которым восстановлена интегральная повторяемость приращений перегрузок в центре масс самолёта при таких полётах.

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

на пожаротушение превышает соответствующую повторяемость согласно [6].

Величина усталостного повреждения (7) регулярной зоны конструкции за один типовой полёт с учётом Б0 и Сп, найденных по результатам измерений перегрузок при полётах на пожаротушение, оказывается выше в 2,75 раза.

у

'к як X \

к х — V \

00 2 0 4\ 0 ] 1 1 2 1

¡Лч

"ч?

—■

от о

н ш с о с

га

п

л н о о 2 ш к о. о н ш о

-2

-5

Приращение перегрузки пу

-♦-Бу-200ЧС транспортная функция -■-Ту-134

-X-RRJ-95LR —*— Бе-200ЧС, пожарная функция, ОСТ

-Ж-Бе-200ЧС, пожарная функция ЛИ

4

3

0

4

Рис.1. Интегральная повторяемость приращения перегрузки в центре масс самолётов на воздушном режиме за один полёт

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

1. Зависимость ресурсных характеристик среднемагистральных пассажирских самолётов от веса конструкции планера [Текст] / В.С. Дубинский, М.Д. Зайцев, Е.И. Крючков [и др.] // Тр. ЦАГИ. - Вып. 2658. -М., 2002. - С. 207 - 214.

2. Машиностроение. Энциклопедия. Т. 4-21. Самолёты и вертолёты. Кн. 1. Аэродинамика, динамика полёта и прочность [Текст] / ред. В.Г. Дмитриев. - М.: Машиностроение, 2002. - 799 с.

3. Нестеренко, Г.И. Ресурс и живучесть самолётных конструкций [Текст] / Г.И. Нестеренко // Проблемы машиностроения и надёжности машин. - 2005. -№1. - С. 106 - 118.

4. Левин, А. С. Уточнение безопасного ресурса критических элементов конструкции крыла «стареющего» пассажирского самолёта [Текст] / А.С. Левин, В.Е.

Стрижиус // Науч. вестн. Моск. гос. техн. ун-та гражд. авиации. Сер. Аэромеханика, прочность, поддержание лётной годности ВС.- 2004. - Вып. 84 (2). - С. 35 - 43.

5. Стрижиус, В.Е. К расчету эквивалентов программ усталостных испытаний крыльев транспортных самолётов [Текст] / В.Е. Стрижиус // Ученые записки ЦАГИ. -2000. - № 3-4. - С. 178 - 191.

6. ОСТ 1 02514-84 Модель турбулентности атмосферы. - Введ. 01.01.1986. - 13 с.

7. Стрижиус, В.Е. Типизированная программа нагружения элементов конструкции крыла регионального транспортного самолёта [Текст] / В.Е. Стрижиус // Науч. вестн. МГТУ ГА. Сер. Аэромеханика и прочность, поддержание лётной годности ВС.- 2006. - № 103. - С. 145154.

ASSURING THE DURABILITY OF TRANSPORT AIRCRAFT STRUCTURE IN MULTIPURPOSE USE BASED ON THE RESULTS OF ENDURANCE AND FLIGHT TESTS

© 2013 P. A. Fomichev1, V. I. Kleptsov2

1National Aerospace University named after N.E. Zhukovsky «KhAI» Beriev Aircraft Company plc, Taganrog

An approach is proposed to determine the durability of transport aircraft wing structure in multi-purpose use. Relationships are obtained for durability calculations based on fatigue test results depending on typical flight profiles and their relative flying time in service. Dependences are obtained for the calculation of damages within the framework of the known discrete model of atmospheric turbulence. The possibility of refining the fatigue damage calculation using the experimental cumulative probability of flight overload increment is shown. Applying to Be-200ChS aircraft, the fatigue damage in firefighting flights obtained from flight test data was found to be 2.75 times greater than the value calculated by the known method.

Durability of structure, fatigue damage, flight profile, cumulative probability, overload, tests.

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

Фомичёв Петр Александрович, доктор технических наук, профессор, заведующий кафедрой прочности летательных аппаратов, Национальный аэрокосмический университет им. Н.Е. Жуковского «ХАИ». E-mail: nil prochnost@khai.edu. Область научных интересов: методы расчёта долговечности элементов конструкций при нерегулярном нагружении.

Клепцов Виктор Иванович, заместитель начальника отдела прочности, ОАО «Таганрогский авиационный научно-технический комплекс им. Г.М. Бериева». E-mail: kleptsov@beriev.com. Область научных интересов: наземные и лётные прочностные испытания, ресурс воздушных судов.

Fomichev Pyotr Aleksandrovich, doctor of technical science, professor, head of the department of aircraft strength, National Aerospace University named after N.E. Zhukovsky «KhAI». E-mail: nil prochnost@khai.edu. Area of research: methods of calculating durability under irregular loading.

Kleptsov Victor Ivanovich, deputy chief of the strength department, Beriev Aircraft Company. E-mail: kleptsov@beriev.com. Area of research: ground and flight strength tests, durability of aircraft.