УДК 629.7.017
А. Г. Зосимов, О. Г. Бойко, Л. Г. Шаймарданов
МЕТОД РАСЧЕТА НАДЕЖНОСТИ ФУНКЦИОНАЛЬНЫХ СИСТЕМ САМОЛЕТОВ ГРАЖДАНСКОЙ АВИАЦИИ ПО СТАТИСТИЧЕСКИМ МАТЕРИАЛАМ ЭКСПЛУАТАНТОВ
Разработан метод расчета вероятности отказа функциональных систем на 1ч полета по статистическим материалам эксплуатантов.
Статистические материалы эксплуатантов (авиакомпаний, регионов и отрасли гражданской авиации в целом) [1; 2], а также зарубежной авиационной техники [3] убедительно показывают, что средний налет на отказ самолетов в целом, их функциональных систем и отдельных агрегатов для конкретных типов самолетов не зависит от времени (от налета часов самолетов). Постоянство среднего налета на отказ обеспечивается реализацией программ технического обслуживания самолетов, направленных на поддержание уровня летной годности, т. е. надежности агрегатов, систем и самолетов в целом.
Принятый в расчетах надежности функциональных систем метод основывается на теореме умножения вероятностей. В соответствии с ней вероятность безотказной работы последовательно соединенных агрегатов принимается равной произведению вероятностей безотказной работы этих агрегатов. Вероятность отказа параллельно включенных агрегатов определяется в виде произведения вероятностей их отказа.
Сложные функциональные системы содержат большое число агрегатов и имеют многократное общее и индивидуальное резервирование отдельных агрегатов. В связи с этим расчетные зависимости для вероятности отказа либо безотказной работы систем содержат суммы вероятностей отказа агрегатов, суммы их степеней и произведений.
Центральная предельная теорема теории вероятностей [4] утверждает, что распределение вероятности суммы случайных величин неограниченно приближается к нормальному закону распределения при неограниченном увеличении числа слагаемых. В функциональных системах самолетов число агрегатов достаточно велико, но, естественно, ограниченно.
Рассмотрим вид интегральной функции распределения отказа системы и ее плотности вероятности на примере расчета тестовой системы. Пусть система имеет 4-кратное общее резервирование ветви содержащей 4 последовательно соединенных агрегата. Положим, что параметр потока отказа является одинаковым для всех агрегатов и равным 1 • 10-4.
В отраслевой нормали «Надежность изделий авиационной техники» [5] для малых, свойственных агрегатам авиационной техники, значений параметра потока отказов предложено распределение вероятностей отказа агрегатов представить законом равномерной плотности вероятности, т. е.
<2(г) = юг. (1)
Следует отметить, что распределение (1) при малых вероятностях (2(г) хорошо согласуется с экспоненциальным распределением.
При этих предположениях вероятность отказа тестовой системы определяется в виде
ес(0 = [1 - (1 - 1 • 10-4;)4]4. (2)
Рассчитанная по формуле (2) функция вероятности отказа системы приведена на рис. 1. Там же представлена аналогичная зависимость, но построенная в предположении экспоненциального закона распределения вероятности отказа агрегатов.
Вероятность отказа систем авиационной техники нормируется Нормами летной годности самолетов [5] в виде вероятности отказа на 1ч полета, приводящего к последствиям различной степени тяжести. Для рассматриваемой системы эти зависимости приведены на рис. 2. Характер изменения во времени функции распределения вероятности отказа системы (рис. 1) и вероятности ее отказа на 1ч полета близки к интегральной функции распределения вероятности и плотности вероятности нормального распределения.
Здесь следует иметь в виду, что плотность вероятности определяется по формуле
ёд (г)
/ (і )=
а вероятность отказа на 1ч полета в виде
в1с (і ) =
де (і)
м
(3)
(4)
Поскольку Дг = 1 ч, а г изменяется в пределах десятков тысяч часов, допустимо считать Дг величиной 4-го порядка малости. Но тогда зависимости /(г) и Qc1(t) мало чем отличаются друг от друга.
При использовании закона распределения с равномерной плотностью вероятности (1) плотность вероятности равна
/ ч же (і)
/ (і) =--------------= ю = сош!,
у 1 йі
(5)
т. е. в этом случае плотность вероятности и параметр потока отказов тождественны. Тогда, приняв за единицу времени 1 ч, допустимо считать, что параметр потока отказов есть не что иное, как вероятность отказа на 1 ч.
Таким образом, рассматриваемая тестовая система содержит агрегаты, вероятность отказа которых в течение 1 ч не зависит от времени. Естественным было бы ожидать независимость от времени и вероятности отказа системы на 1 ч. Эти представления расходятся с результатом расчета системы, представленном на рис. 2.
В рассмотренной методике задача расчета надежности решена относительно интегральной функции распределения вероятности отказа системы, а вероятность отка-
за на 1 ч определена как производная от нее. При построении интегральной функции распределения в рассмотрение введено время г, которое затем возводится в степени и суммируется многократно в зависимости от сложности системы. При этом остается проигнорированным факт независимости от времени вероятности отказа агрегата на единицу времени.
Ниже предлагается альтернативный методологический подход к построению расчета надежности системы. При этом на первом этапе рассматривается задача расчета вероятности отказа системы на 1ч полета, а функция распределения вероятности отказа находится как результат ее интегрирования.
Рассматривается тестовая система, состоящая из п агрегатов и имеющая т кратное общее резервирование. Вероятности отказов агрегатов приняты одинаковыми и
Рис. 1 Вероятность отказа системы
общее резервирование _ . — . . общее резервирование при
экспоненциальном распределении
Рис. 2 Вероятность отказа системы на 1 ч полета общее резервирование _ л _ я ш общее резервирование при
экспоненциальном распределении
распределенными по закону равномерной плотности вероятности. Вводится понятие суммарного потока отказов системы юс, определяемого суммой потоков отказов агрегатов. Поскольку агрегатов п, а их параметры пока отказов ю одинаковы, то
Юс = п ю. (6)
Реализация отказа с параметром потока юс не отождествляется с отказом системы и не идентифицируется с метом его возникновения.
Необходимым и достаточным условием отказа системы в течении 1 ч является реализация в каждой из т резервированных подсистем с последовательным соединение п/т агрегатов хотя бы одного отказа. Вероятность возникновения в одной из т подсистем отказа одного агрегата в т раз меньше, чем в системе в целом. Поскольку подсистем т, то вероятность отказа одного агрегата в каждой из них в т2раз меньше. С другой стороны, вероятность одновременной реализации т независимых отказов равна произведению вероятности этих отказов.
Поскольку вероятность отказа агрегатов принята распределенной по закону равномерной плотности вероятности, допустимо считать вероятность реализации отказа одного агрегата в системе равной ю . Тогда вероятность отказа системы в течение 1 ч определится в виде
в =
ю"
(п ю)
(7)
т т
Это значение вероятности отказа на 1ч полета не зависит от времени, а интегральная функция распределения вероятности отказов системы Q является линейной функцией времени
(п • ю)т
вс =* ■
(8)
Для рассмотренной выше системы из 16-ти агрегатов вероятность отказа на 1ч будет равна (16 ю)
вс1 =
16
= 163 -ю4.
(9)
В рассматриваемом числовом примере при традиционной методике расчета и ю = 10-4 вероятность отказа на 1ч полета на отрезке времени (0,1) составляет 0,256 • 10-13, среднее значение в интервале от 0до 1 000 ч оставляет 0,14 • 10-4и максимальное значение, достигаемое при г = 3 500 ч, составляет 0,18 • 10-3. При предложенном альтернативном решении задачи вероятность отказа на 1ч полета равна 0,409 6 • 10-12 и не зависит от наработки часов.
Значение интегральной функции распределение вероятности отказа системы в первом случае равно 1 при г = 10 000 ч, во втором случае при этой же наработке часов оно составляет 0, 409 6 • 10-8.
Библиографический список
1. Анализ влияния надежности авиационной техники на безопасность полетов за 2003 г. Минтранс России. М., 2004. 102 с.
2. Анализ инженерно-авиационного обеспечения безопасности полетов в авиапредприятиях Красноярского МТУ ВТ МТ РФ за 2003 год. Красноярск, 2004. 24 с.
3. ДОК 9642-АК/941. Руководство по сохранению летной годности. Международная организация гражданской авиации. 1995 г.
4. Венцель, Е. С. Теория вероятности / Е. С. Венцель. М. : Изд-во Физикоматематической литературы, 1962. 564 с.
5. ОСТ 1 00132-84. Надежность изделий авиационной техники. Методы количественного анализа безотказности функциональных систем при проектировании самолетов и вертолетов.
A. G. Zosimov, O. G. Boyko, L. G. Shaimardanov
THE RELIABILITY CALCULATION METHOD OF THE CIVIL AVIATION AIRPLANES FUNCTIONAL SYSTEMS BY THE STATISTIC MATERIALS OF EXPLOITER
The possibility calculation method of the functional systems breakdown for an hour flight by the statistic materials of exploiter is developed.
УДК 629.7.017
О. Г. Бойко, Л. Г. Шаймарданов
МЕТОДОЛОГИЧЕСКИЕ ОСОБЕННОСТИ РАСЧЕТОВ НАДЕЖНОСТИ ФУНКЦИОНАЛЬНЫГХ СИСТЕМ САМОЛЕТОВ ГРАЖДАНСКОЙ АВИАЦИИ
Рассмотрены! методические противоречия в расчетах надежности сложных систем, обусловленные проявлением центральной предельной теоремы теории вероятностей, приводящим к реализации нормального распределения вероятности отказа в системах со стационарными параметрами потока отказов агрегатов.
Функциональные системы самолетов осуществляют выполнение многих операций, обеспечивающих безопасность полетов. Они представляют собой сложные разветвленные системы, в которых использовано как общее, так и индивидуальное резервирование. Так, гидросистема самолета Ту-154 М содержит 79 агрегатов, имеет трехкратное общее резервирование и 6-кратное резервирование по гидронасосам. В гидросистеме самолетов Ил-86 и Ил-96-300 применено 4-кратное общее резервирование и 8-кратное резервирование насосных агрегатов. Не менее сложными являются системы управления и механизации планера, топливные системы, системы кондиционирования и др.
Показатели надежности самолетов в целом и их функциональных систем задаются государственными и международными требованиями, отраженными в Нормах летной годности самолетов (НЛГС) [1; 2]. В НЛГС требования к надежности определены в виде допустимых на 1ч полета значений вероятностей отказов систем, приводящих к возникновению в полете неблагоприятных ситуаций, классифицированных как сложная, аварийная и катастрофическая ситуация.
В эксплуатационных подразделениях надежность агрегатов и систем контролируется по косвенным показателям, таким как налет часов на 1 отказ по самолету вце-лом, по системам, по агрегатам; налет на отказ, приводящий к определенным неблагоприятным последствиям (инцидент, серьезный инцидент, авиационное происшествие); количество отказов на 1 000 ч налета (К1 000) и др. Надежность агрегатов и систем самолетов иностранного производства в эксплуатации оценивается в косвенных показателях, не имеющих принципиальных отличий от используемых в гражданской авиации России.
Косвенные показатели оценки надежности задаются головными разработчиками типов самолетов и должны отражать требования НЛГС и облегчить эксплуатантам процедуры контроля надежности.
Крупные авиакомпании, эксплуатирующие десятки однотипных самолетов, принципиально имеют право разрабатывать свои программы технического обслуживания самолетов, направленные на уменьшение эксплуатационных расходов и обеспечивающие поддержание надежности на уровне установленных требований. Разработка таких программ сопряжена с высокой степенью ответственности и должна базироваться на анализе надежности агрегатов и систем выполненном в прямых показателях определяемых НЛГС.
Теория надежности является развитым научным направлением, в рамках которого разработаны методики расчета надежности агрегатов и сложных систем в целом. Вместе с этим следует отметить, что предлагаемые методики не доводятся до расчетов надежности на коротком интервале времени, например на 1ч полета, нормируемой в НЛГС.
При анализе надежности агрегатов и систем в пределах малого интервала выявлен ряд некорректностей, присущих принятым методикам, которые никак себя не проявляют при расчете функций распределения вероятности отказов либо безотказной работы. Обсуждению этих некорректностей и посвящена данная работа.
Рассматривается простой случай, когда интенсивность отказов X либо параметр потока отказов ю являются постоянными величинами, не зависящими от времени г наработки агрегата. Это предположение подтверждается опытом эксплуатации агрегатов, систем и самолетов в