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

V. A. Fel’k, Yu. S. Vorontsov, A. N. Fomin, A. I. Stupnikov MEASUREMENT OF DISTANCE WITH THE HELP OF THE MODIFIED DOPPLER EFFECT

The method of high-precision measurement of distance is suggested with the use of the modified Doppler effect. Keywords: Doppler effect, phase self-tuning of frequency.

© Фельк В. А., Воронцов Ю. С., Фомин А. Н., Ступников А. И., 2011

УДК 629.735.064

Л. Г. Шаймарданов, О. Г. Бойко

МЕТОД РЕШЕНИЯ ЗАДАЧИ РАСЧЕТА НАДЕЖНОСТИ СЛОЖНЫХ СИСТЕМ ПРИ ПЕРЕМЕННЫХ ПАРАМЕТРАХ ПОТОКОВ ОТКАЗОВ АГРЕГАТОВ

Предложен метод эквивалентной линеаризации переменного параметра потока отказов агрегатов.

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

В работах [1-3] предложен метод расчета надежности сложных функциональных систем самолетов гражданской авиации при стационарности потока отказов ю. В соответствии с ним вероятность отказа функциональной системы рассчитывается на дискретных отрезках времени [0, /].

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

q(t) = ю • t. (1)

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

Приняв t = 1, из выражения (1) можно определить вероятность отказа за 1 ч, как этого требуют нормы летной годности самолетов:

q(1) = ю • 1 = Юь (2)

Для принятия дискретного отрезка времени [0, {] протяженностью т как единичного (например, 3 600 с или продолжительность типового полета самолета) в указанных работах предложено приводить к т параметр потока отказов, при этом

ют = ю • т. (3)

Тогда выражение (2) запишется в виде

q(т) = юх • 1. (4)

При определении вероятности отказа агрегата на отрезке т, принятом за единичный, интегральная

функция определена уже как дискретное событие. Это

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

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

0С = {1 - [1 - q (т)]Т. (5)

В отмеченных работах по выражению (5) определяется вероятность отказа системы только на дискретных отрезках времени [0, /]. Ему не приписываются свойства интегральной функции вероятности отказа системы.

В монографии [4] отмечается, что даже у обслуживаемых (восстанавливаемых) агрегатов потоки отказов могут быть медленно меняющимися во времени, т. е. ю = var. При работе необслуживаемых систем, имеющих различные виды резервирования (горячий резерв, холодный резерв, скользящее резервирование), потоки отказов агрегатов, естественно, возрастают по мере увеличения наработки. В связи с этим задача построения решения для расчета надежности системы при переменных по времени параметрах потоков отказов агрегатов является актуальной.

Для решения поставленной задачи необходимо заменить на единичном отрезке времени [0, /] переменный параметр потока отказов эквивалентным постоянным. Для этого можно использовать метод, разработанный в [1-3]. В теории колебаний для решения задачи с нелинейным трением широко применяется метод эквивалентной линеаризации. В ряде случаев нелинейные дифференциальные уравнения с нелинейным диссипативным звеном не имеют решения. Для приведения таких уравнений к линейным дифференциальным уравнениям и используется этот метод.

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

В рассматриваемой работе предлагается метод замены на отрезке времени [0, г] переменного параметра потока отказов ю(г) постоянным юэ при условии равенства вероятности отказа агрегата q(t) при постоянном и переменном параметрах потоков отказов (рис. 1).

Рис. 1. Замена переменного параметра потока отказов ю(г) постоянным юэ

Интегральная функция вероятности отказа агрегата вида (1) предполагает параметр потока отказов ю численно равным плотности вероятности распределения с равномерной плотностью. Тогда условие эквивалентности ю(г) и юэ на отрезке [0, г] определится в виде

откуда

= | ю(т^т .

0

г

| ю( т)d т

_0_________

г

(6)

(7)

ю(г)

0,018

0,016

0,014

0,012

0,01

0,008

0,006

0,004

0,002

0

10 20 30 40 50 60 70 80 г

Рис. 2. Пример зависимостей параметра потока отказов от времени:

— а = 0,01, к = 0;-----а = 0, к = 0,0002;-----а = 0,004,

к = 0,00012;......а = 0,008, к = 0,00004

ю(г)

0,9

0,8

0,7

0,6

0,5

0,4

0,3

0,2

0,1

/Г

^у,

■У'

/ г

**'/ ✓

^ 4 /

' , ►

* О

♦ ♦ ✓

0

10 20 30 40 50 60 70 80 г

В левой части равенства (6) юэ принимается постоянным на отрезке времени [0, г]. Вместе с этим при изменении протяженности отрезка [0, г] изменяется и величина юэ. Но в пределах рассматриваемого отрезка [0, г] запись ю э • г означает не что иное, как определение вероятности отказа на отрезке при распределении с равномерной плотностью вероятности. Это обеспечивает возможность построения для агрегата интегральной функции распределения вероятности отказа, если задаться различными значениями протяженности отрезка [0, г].

В качестве иллюстрации предлагаемого метода (рис. 2) приведены различные реализации переменного значения параметра потока отказов ю(г) в виде

ю(г) = а + кг (9)

для ряда значений а и к.

Для рассматриваемого примера приведены интегральные функции распределения вероятности отказа агрегата, построенные в соответствии с предложенным методом эквивалентной линеаризации (рис. 3).

Рис. 3. Интегральные функции распределения вероятности

отказа агрегата, построенные методом эквивалентной линеаризации:

---а = 0,01, к = 0;-а = 0, к = 0,0002;-а = 0,004,

к = 0,00012;..а = 0,008, к = 0,00004

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

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

1. Бойко О. Г., Шаймарданов Л. Г. Математическое моделирование схемной надежности сложных систем // Проблемы безопасности и чрезвычайных ситуаций. 2010. Вып. 3. С. 82-88.

2. Бойко О. Г. Надежность функциональных систем самолетов гражданской авиации : монография // Избранные труды Российской школы по проблемам науки и технологий / РАН. М., 2009.

3. Бойко О. Г., Шаймарданов Л. Г. Моделирование надежности агрегатов функциональных систем самолетов // Проблемы машиностроения и надежности машин. М., 2010. Вып. 5. С. 35-39.

4. Новожилов Г. В., Неймарк С., Цесарский Л. Г. Безопасность полетов, концепция и технологии. М. : Машиностроение, 2003.

L. G. Shaymardanov, O. G. Boyko

METHOD FOR SOLUTION OF A PROBLEM OF CALCULATION OF RELIABILITY OF COMPLICATED SYSTEMS WITH VARIABLE PARAMETERS OF BREAKDOWN OF AGGREGATES

The method of equivalent linearization variable parameters of stream of breakdown of aggregates is presented in the article

Keywords: functional systems reliability, breakdown probability integral function, breakdown stream parameter.

© Шаймарданов Л. Г., Бойко О. Г., 2011