Анализ безопасности гидравлической системы самолета Ту-154М по результатамвыделения предвестников аварии Текст научной статьи по специальности «Математика»

тающих на щелочных металлах, труды МАИ, Вып. 193. М. :

Машиностроение, 1969. С. 57-102.

4. Краев М. В., Овсянников Б. В., Шапиро А. С. Гидродинамические уплотнения высокооборотных валов. М. :

Машиностроение. 1976.

E. M. Kraeva, I. S. Masich VORTEX STRUCTURE OF TURBULENT FLOWS AND THEIR SIMULATION

5. Хаген Р. Л., Данак А. М. Перенос импульса при турбулентном отрывном обтекании прямоугольной впадины // Тр. Амер. о-ва инженеров-механиков. Сер. Д, Прикладная механика. 1966.№3.С. 189-195.

The work presents analysis of the basic parameters of vortex turbulence in the field of centrifugal forces. The dependencies ofmathematical modeling which allow to calculate turbulent vortex structures, are obtained.

Keywords: model, turbulence parameters, vortex flow.

© Краева Е. М., Macm H. C., 2011

УДК 629.735.017

А. В. Любимцев

АНАЛИЗ БЕЗОПАСНОСТИ ГИДРАВЛИЧЕСКОЙ СИСТЕМЫ САМОЛЕТА ТУ-154М ПО РЕЗУЛЬТАТАМ ВЫДЕЛЕНИЯ ПРЕДВЕСТНИКОВ АВАРИИ

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

Ключевые слова: предвестники аварии, безопасность, летная эксплуатация, рейтинги нарушений, отказы.

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

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

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

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

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

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

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

Рассмотрим особенности расчёта риска и выделения предвестников аварий на примере анализа безопасности гидравлической системы самолетного парка Ту-154М. От гидравлической системы на Ту-154 осуществляется управление гидроусилителями систем управления всеми рулевыми плоскостями, торможением колес, выпуском и уборкой шасси. Отказ этой системы приводит к потере управления лайнером. Согласно статистике Госслужбы Гражданской авиации (ГА) и управления госнадзора за безопасностью полетов к 2009 г. отказы агрегатов в системах воздушных судов (ВС) распределились следующим образом: АиРЭО - 20,6 %, гидросистема - 24 %, системы управления - 17,6 %, двигатель и его агрегаты - 13 %, шасси - 13 %, системы кондиционирования, высотная, противопожарная системы -10%, прочие - 3 %. Наиболее часто отказы агрегатов в гидравлических системах происходили на самолетах Ту-154М. Причем эти отказы продолжают повторяться из года в год, и их число неуклонно растет. В 2008 г. имело место 7 инцидентов, связанных с гидросистемой, в 2009 г. их было уже 14.

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

Анализ развития ИС в аварийную последовательность исходов удобно выполнять с использованием метода «деревьев событий». «Дерево событий» представляет собой граф, включающий исходное событие, промежуточные события и конечные состояния. Граф наглядно показывает возможные пути развития аварийной последовательности. «Дерево событий» строится с целью анализа последствий некоторого исходного события -10, которое изображается в основании дерева, представляющего собой логическую диаграмму, которая определяет множе-

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

Для каждого нарушения, неисправности или отказа агрегата гидросистемы, выделим N возможных аварийных последовательностей (АП), пронумерованных как Е1,..., EN, и для каждой из них рассчитаем соответствующую условную вероятность Qi (I =1, .■■,N перехода нарушения в аварию. Другими словами каждому нарушению, которое является исходным событием, поставим в соответствие N вероятностей перехода этого нарушения в аварию. Авария может наступить при реализации любой из N несовместных аварийных последовательностей. Таким образом, авария - это событие (в вероятностном смысле), которое является суммой несовместных событий Е1,..., Е1. Следовательно, вероятность аварии Qio описывается формулой

Qio= ^ Е//)^=1,2,..., N (1)

/=1

где Qi(Ei / I) - вероятность реализации ^й АП для данного ИС. °

Тогда полная (безусловная) вероятность аварии при наступлении исходного события 1о может быть рассчитана по выражению:

Я(Г)=Р(Г)Е//),I =1,2, ...,^ (2)

где Р(1) - вероятность наступления исходного события I.

Графически «дерево событий» изображается в виде таблицы состояний и дихотомического правостороннего дерева.

В качестве примера построим «дерево событий» и рассчитаем вероятность исходного события на примере отказа, связанного с многоочаговым усталостным разрушением корпуса гидравлического насоса НП-89 первого двигателя. Тяжесть данного нарушения надо оценивать в виде вероятности возможного перехода его в аварию. Проанализируем потенциальные аварийные последовательности, которые могли бы реализоваться в процессе развития такого нарушения и рассчитаем их соответствующие риски.

Гидронасос НП-89 служит для создания и поддержания рабочего давления в гидросистеме. Гидронасос плунжерной конструкции приводится в действие от коробки приводов двигателя [3]. Более 50 % всех инцидентов, связанных с отказами гидросистем самолетов Ту-154, были обусловлены низкой надежностью гидронасосов НП-89. За последние годы произошло более 100 случаев отказов насосов. Отказы отмечаются как на новых, так и на отремонтированных агрегатах в пределах существующего межгарантийного и назначенного ресурсов.

В основание дерева поместим наше исходное событие «многоочаговое усталостное разрушение корпуса гидравлического насоса НП-89». Аварийное состояние, связанное с разрушением гидронасоса, обусловлено следующими промежуточными событиями:

- ^: значительной потерей гидрожидкости из бака первой и второй гидросистем с последующим падением давления в них до 150 кг/см2 и одновременным загоранием светосигнализатора «падения давления в ГС»;

- !2: работа второго НП-89 первой гидросистемы на втором двигателе в условиях масляного голодания, вызывающего кавитационный износ его узлов, создание условий для отказа второго НП-89 из-за заклинивания, и срез шлицевой муфты привода гидронасоса из-за масляного голодания;

- !3: утечки масла из коробки приводов, в результате разрушения корпуса гидронасоса, что приводит к отказу двигателя;

-14: возникновение автоколебаний давления первой гидросистемы линии нагнетания с размахом, значительно превышающим допустимые из-за неравномерной работы НП-89 и разрушение его корпуса.

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

Таким образом, для рассматриваемого исходного события !0 выделено девять несовместных аварийных последовательностей Е1, ., Е9. Вероятности реализаций каждого пути развития аварии определяем по формуле вероятности произведения независимых событий:

Р(Е1/7) = 0,99994 • 0,99999 • 0,999999 • 0,999999 * 0,99999.

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

Р(Е2/I) * 9 • 10-7; Р(Е3/I) * 3- 10-15;

Р(Е4/7) * 6 • 10-6; Р(Е5 //) * 1,2 • 10-11;

° P(E6/Io) * 1,2 • 1011;

Р(Е7 //)* 1,3 • 10-17;Р(Е8/Io)*5 • 10-17;

0 P(E9/Io) * 6 • 10-17.

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

- для первой аварийной последовательности

Q1(Io) = 9•10-7+6•10-6+1,2•10-11 +

+1,2 • 10-11 + 3^ 10-15= 0,000 006 09.

Тогда полная (безусловная) вероятность аварии при наступлении исходного события ^ примет вид

R1(Io)=P1(Io)•Q1(Io)=10-s•6,09•10-6=6,09•10-11.

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

- повышенные по сравнению с установленными техническими условиями силы трения между золотником и гильзой тормозного клапана УГ-92 - R2(I) = 1,6 • 10-11;

- нарушение целостности цепи питания электромагнитного крана ЭК-47 от электрошины -Ю(/)) = 2 • 10-12;

- кратковременное заклинивание золотника золотникового пульта РГ-16А в результате попадания посторонних частиц в рабочую полость - _К4(! ) = 2^ 10-11;

- разрушение возвратной пружины золотника крана подачи давления на выпуск-уборку внутренних интерцеп-торов-_К5(/)) = 3 • 10-11;

- обрыв штока-компенсатора усилия рулевого привода РП-59 управления элеронов-интерцепторов -^(/^•Ю12;

- трещина в корпусе гидромотора управления зак-рылками-Л7(/) = 3 • 10-12;

- разрушение трубопровода подвода гидрожидкости к фильтру-R8(Io) = 4,7• 10-12;

- разрушение дросселя постоянного расхода -^(/^¿•Мг11;

- усталостное разрушение гидрофильтра - Ш0(!) =

= 8 • 10-12. 0

Предполагается, что рейтинги нарушений рассмотренных выше отказов Я1 = 5(1). ...К 10 = 5(10) образуют выборку значений рейтинга из бесконечной генеральной совокупности. Полагая, что рейтинг 5 подчиняется нормальному распределению, в [1] предложено по выборочным характеристикам 5 и Ds (5-оценка среднего значения рейтинга, Ds-оценка дисперсии рейтинга) найти такой предел, что с доверительной вероятностью Y > 0,5 можно гарантировать попадание большей доли совокупности значений рейтинга в интервал [О, 5 ], где 5пр - в математической статистике определен как верхний толерантный предел. Рассчитывается он по формуле для нормального распределения:

5 =5 + k^Ds, (3)

пр ---------- ’ 4 '

где к - толерантный множитель.

Значения толерантного множителя выбирают из таблицы в зависимости от доверительной вероятности - Yи объема наблюдений г. Выбор значений Yи г обусловлен задачами анализа безопасности. Для обеспечения безопасности технических систем значения Yцелесообразно выбирать на уровне 0,9 и выше. Количество наблюдений примем равным г =10, тогда толерантный множитель к =1,842.

Значения 5 и Ds вычисляют по формулам [ 1 ]:

Ds = Е[5(0-5]2/г, (4),

5 = £5(0 / г. (5)

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

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

5 = (¿0/10 = 0,77 • 10-11, Ds = (^г)-5)2/ 10= 1,2 • 10-11.

Тогда значение верхнего толерантного предела 5пр = = 5+к • * 2,98 • 10-11. Таким образом, с высокой вероятностью можно отнести нарушения ^(1), ^(5), Я(9) к предвестникам аварии, так как значение их рейтингов наруше-

Исходное событие Промежуточные состояния Конечное состояние вероятно ст :о стояния

Многоочагово е устало стно е разрушение корпуса гидравлического насоса НП-89 первого двигателя Значительная потеря гидр о жидкости из бака первой и второй гидросистемы до 4 литров с последующим падением давления в гидросистемах до 150 кг/с м2 с одновременным загоранием красного свето сигнализатора “падение давления в первой ГС” Работа второго НП-89 первой гидросистемы На втор ом двигателе в условии масляного голодания, вызывающего кавитационный износ его узлов (поверхностей). Созданы условия для отказа второго НП-89 из-за несовершенства конструктивной схемь со единения насосов, заклинивание и срез шлицевой иуфты привода второго НП-89 из-за масляного голодания При разрушении корпуса гидр о н ас о са н ачин аются утечки масла из коробки приводов, что приводит к отказу двигателя. давления первой гидросистемы линии нагнетания с размахом значительно превышающим д опусти мые из -з а н ер авн о мер н ой работы НП-89 и разрушения его корпуса

Нет

Нет

1о

Исходное событие

Нет 0.999999

0.99994 Нет 0.99999 0.999999 Да

Да ю-6 Нет 0.999999

Да

Нет Нет 0.999999

Да б* ю'3 Нет Ла

0.999998 Да 2*10'“ РТ^т 0 999999

1

Нет Дп Ё1ет 0.999999

да 2*10“ Па

Да 0.5*10° п чччччч

Да

0.99999

9*10'''

з*ю^

6*10° 1.2*10 11

1.2*10'“

1.3*10-*''

5*1 О"1' б* 10'1'

Рис. 1. «Дерево событий», реализующееся при разрушении корпуса НП-89

нии резко выделяется из других значении, которые лежат ниже £ .

пр

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

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

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

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

1. Александровская Л. Н., Аронов И. З. Статистические методы анализа безопасности сложных технических систем. М.: Логос, 2001.

2. Аронов И. З. Современные проблемы безопасности технических систем и анализа риска // Стандарты и качество. 1998. №3.

3. Безопасность полетов : учебник для вузов / Р. В. Сакач, В. Н. Голего, М. Ф. Давиденко и др. М.: Транспорт, 1989.

Рис. 2. Контрольная карта безопасности гидравлической системы самолета ТУ-154 по результатам выделения предвестников аварии

Таблица 1

Значения толерантного множителя к

У Р г= 10 г= 15 г =20 г= 30

0,80 0,80 1,096 1,035 1,002 0,965

0,90 1,665 1,573 1,523 1,467

0,95 2,140 2,022 1,957 1,886

0,99 3,024 2,845 2,766 2,668

0,90 0,80 1,210 1,124 1,076 1,024

0,90 1,842 1,710 1,637 1,558

0,95 2,365 2,196 2,102 2,000

0,99 3,345 3,105 2,973 2,828

0,95 0,80 1,316 1,205 1,143 1,075

0,90 2,003 1,833 1,739 1,636

0,95 2,572 2,353 2,234 2,100

0,99 3,638 3,328 3,159 2,970

0,99 0,80 1,539 1,373 1,282 1,181

0,90 2,342 2,089 1,950 1,797

0,95 3,007 2,683 2,504 2,308

0,99 4,252 3,794 3,540 3,263

A. V. Lyuibimtsev

SECURITY ANALYSIS OFAIRPLAIN TU-154M HYDRAULIC SYSTEM AS A RESULT OF DETECTION OF ACCIDENT MESSENGERS

In the article the authors analyze statistical performance data on failures and incidents of airplane Tu-154M hydraulic system, taken as basis for identification of violations with the highest probability values of transition into accident, and blowouts as accident messengers.

Keywords: breakdown hurbingers, safety, flight exploitation, violations ratings, breakdowns.

©Любимцев А. В., 2011