УДК 681.3 (075)
О.Д. КОЗАЧУК*
ДОСЛІДЖЕННЯ ВПЛИВУ КОЕФІЦІЄНТА НАДМІРНОСТІ НА ПОКАЗНИКИ НАДІЙНОСТІ РЕЗЕРВОВАНИХ СИСТЕМ З ДРОБОВОЮ КРАТНІСТЮ
Національний авіаційний університет, Київ, Україна
Анотація. Розглянуто вплив коефіцієнта надмірності на показники надійності резервованої системи з дробовою кратністю за різними законами розподілу часу безвідмовної роботи елементів системи.
Ключові слова: надійність, резервована система, дробова кратність, закон розподілу, час безвідмовної роботи, коефіцієнт надмірності, раптові відмови, поступові відмови.
Аннотация. Рассмотрено влияние коэффициента избыточности на показатели надежности резервированной системы с дробной кратностью при различных законах распределения времени безотказной работы элементов системы.
Ключевые слова: надежность, резервированная система, дробная кратность, закон распределения, время безотказной работы, коэффициент избыточности, внезапные отказы, постепенные отказы.
Abstract. The influence of redundancy coefficient for reliability redundant system with fractional multiplicity for different distribution laws of no-failure operation time of system elements was regarded. Keywords: reliability, redundant system, fractional multiplicity, distribution law, no-failure operation time, redundancy coefficient, sudden failures, degradation failure.
1. Вступ. Постановка проблеми
Підвищення ефективності та надійності технічних систем за рахунок застосування резервування (додаткової надмірності) у структуру системи - важлива та актуальна проблема. Одним із широко застосованих методів резервування систем є резервування систем із загальним навантаженим резервуванням з дробовою кратністю [1].
Дослідження впливу надмірності на резервовані системи з дробовою кратністю представляє технічний інтерес і має велику практичну значимість.
2. Аналіз досліджень та публікацій
Рішенню цієї проблеми присвячена велика кількість праць [1-3]. Незважаючи на цей вплив коефіцієнта надмірності на показники надійності, дослідження були проведені тільки в резервованій системі елементів з експоненціальним законом розподілу часу безвідмовної роботи.
3. Мета роботи та постановка завдання
Мета роботи - дослідження впливу коефіцієнта надмірності на показники надійності резе-рвованих систем з дробовою кратністю за різними законами розподілу часу безвідмовної роботи елементів системи: експоненціальним, Вейбула, дифузійним немонотонним, узагальненим показовим та суміші розподілів (експоненціальним і Вейбула).
4. Вирішення проблеми
4.1. Під час застосування теореми складання несумісних подій, для випадку однакових елементів системи із загальним навантаженим резервуванням з дробовою кратністю, маємо такий вираз для ймовірності безвідмовної роботи:
© Козачук О.Д., 2015
ISSN 1028-9763. Математичні машини і системи, 2015, № 1
171
Pc(t) = Ci(p(t))N-\l - P(t))' , (1)
де P(t) - імовірність безвідмовної роботи одного елемента;
CN - число поєднання із N елементів за і;
N = п + т - загальне число елементів у резервованій системі; п - число основних елементів; m - число резервних елементів;
Kнадм = m/n ~ коефіцієнтом надмірності резервованої системи (під час загального резервування системи з дробовою кратністю) розуміється відношення числа резервованих (надмірних) елементів системи m до загальної кількості елементів у системі N .
4.2. Дослідження залежності надійності систем із загальним навантаженим резервуванням з дробовою кратністю від коефіцієнта надмірності системи за різними законами розподілу часу безвідмовної роботи елементів, що ураховують випадкові та поступові відмови елементів: експоненціального, Вейбула, дифузійного немонотонного, узагальненого показового та суміші розподілів (експоненціального і Вейбула).
Експоненціальний розподіл часу безвідмовної роботи (однопараметричний - параметр X) ураховує випадкові відмови системи у період нормальної експлуатації системи. Інтенсивність відмов елемента X - стала, не ураховує зношення та старіння у процесі експлуатації.
Розподіл Вейбула часу безвідмовної роботи (двопараметричний - параметри 1 та к ) застосовується під час опису працездатності систем у процесі зношення та старіння.
Дифузійний немонотонний розподіл (ДНР) часу безвідмовної роботи (двопараметричний - параметри а та v ) застосовується для елементів, у яких відомі фізичні процеси деградації, що призводять до відмов у процесі експлуатування, та встановлені визначаючі параметри, що характеризують технічний стан об’єктів [4].
Узагальнений показовий розподіл (УПР) часу безвідмовної роботи - це розподіл тривалості будь-якого складного випадкового процесу, який у початковий момент часу з імовірністю Pi знаходиться на ї-ій стадії, у своєму розвитку може проходити через одну, дві або п стадій, тривалості яких експоненціально розподілені [5].
Суміш розподілів експоненціального і Вейбула часу безвідмовної роботи дозволяє урахувати прояв випадкових відмов та явищ зношення і старіння у процесі експлуатації елементів [5].
4.3. Аналіз впливу коефіцієнта надмірності на показники надійності резервованих систем з дробовою кратністю
4.3.1. У табл. 1 представлені формули для розрахунку показників надійності п’яти різних розподілів часу до відмови елементів. У праці [6] отримані формули та алгоритми для визначення середнього наробітку до відмови систем із загальним навантаженим резервуванням із дробовою кратністю за різними законами розподілу часу безвідмовної роботи елементів.
Для визначення середнього наробітку до відмови систем із загальним навантаженим резервуванням із дробовою кратністю під час дифузійного немонотонного розподілу, узагальненого показового розподілу та суміші розподілів (експоненціального та Вейбула) часу до відмови елементів застосовується рішення трансцендентних рівнянь (2):
1 - K
НАДМ
1
N
Pm (Toc ),
(2)
172
ISSN 1028-9763. Математичні машини і системи, 2015, № 1
де PM (t0) - імовірність безвідмовної роботи елементів:
• для експоненціального розподілу наробіток до відмови елементів:
T 1
■TTL = - ln(1 - KНАДМ - -);
T0M N
для розподілу Вейбула наробіток до відмови елементів:
T
0C
- lnl 1 - K
НАДМ
N
T
0 M
Гі 1 + -
1
b
(3)
(4)
Таблиця 1. Формули для розрахунку показників надійності п’яти різних розподілів часу до відмови елементів
Тип розподілу часу до відмови елементів Формули для ймовірності безвідмовної роботи, щільності розподілу часу до відмови і інтенсивності відмов елементів, L (z ) = fM (z 0) L m{Zq) P (z ) Примітка, z = to ^o t ’ T0 M
Експонен- ціальний розподіл (ЕР) PM (z0 ) = exp(- z0 1 fM (z0 ) = 1M exp(- z0 ), LM (z0 ) = 1 to - час оперативної роботи T =— 1M середній наробіток до відмови
Розподіл Вейбула (РВ) PM (z0 )= exp[-(Kbz0)b 1 fM (z0 ) = bKb (Kbz0)b-1 exp[- (Kbz0)b 1 L m (z0 ) = bKb (Kbz0 )b-1 T0 M = aKb , K Ч1+1) - гамма-функція
Дифузій- ний немонотон- ний розподіл (ДНР) ^1 — z 1 P (z ) Ф 0 exp 2 ІФ ( \ 1 + zo T0M = m 5 V = 1
PM (z0 / _ ф 1— exp Wz0 J 1 fM (z0 )= rz exF mvzj2pz0 2 ф V ) ,[-(1' 2 , Wz0) - Zq )2 1 ^zo j
Узагальне- ний показовий розподіл (УПР) PM (z0) = ^ - Ц j exp(- 2z0 ) +1 exp fM(z>={Mexp(- 2z0 ї-ль+л)' - ^ Zq V h + 1 0 ґ 2h ;xp Л h+1 ) / Y zo )_ T = 2 h + 1 h = 41
Суміш розподілів експоненціального та Вейбула (ЕР+РВ) Pm (Z0) = 0,5exp(-0,625z0) + 0,5exp fM ( Zq ) =l{0,5exp(-0,625 +2bKb (2,5Kb )b-1exp [-(2,5Kb -( 2.5K,Zq )" j zo) + z0)b ]} 5 T0M = 0,5(T0EM + T0WM X T0EM = 1/Л, Tqwm = 0,25/ л, Tqm = 0,625/ Л
ISSN 1028-9763. Математичні машини і системи, 2015, № 1
173
4.3.2. На рис. 1-4 представлені графіки залежності відносного середнього наробітку до відмови від коефіцієнта надмірності резервованої системи за різними законами розподілу часу безвідмовної роботи елементів системи (експоненціальним розподілом (ЕР), розподілом Вейбула (РВ), немонотонним дифузійним розподілом (НДР), узагальненим показовим розподілом (УПР) та сумішшю розподілів експоненціального та Вейбула (ЕР+РВ)).
Відношення середнього наробітку до відмови резервованої системи до середнього наробітку до відмови елементів системи для експоненціального закону розподілу часу безвідмовної роботи елементів,
Т0 C ! То Е
-----N3 = 1024 ел.
-----N2 = 256 ел.
----=N1 = 64 ел.
Коефіцієнт надмірності, Кнадм.
Рис. 1. Графіки залежності відношення середнього наробітку до відмови резервованої системи (під час загального резервування із дробовою кратністю) до середнього наробітку до відмови елементів системи (для експоненціального закону розподілу часу безвідмовної роботи
'0Е
елементів) від коефіцієнта надмірності системи Кна(дм
За графіками на рис. 1 видно, що відношення — поступово зростає:
'0Е
• від 0 до 1,0 під час збільшення Кнадм від 0 до 0,6 для всіх трьох резервованих систем;
• від 1,0 до 6,9 під час збільшення Кнадм від 0,6 до 0,998 для резервованої системи за N = 1024 елементи;
• від 1,0 до 5,8 під час збільшення Кнадм від 0,6 до 0,993 для резервованої системи заN = 256 елементів;
• від 1,0 до 4,16 під час збільшення Кнадмвід 0,6 до 0,969 для резервованої системи заN = 64 элемента.
Графіки відношення — для резервованих систем із дробовою кратністю та чис-
'0Е
лом елементів у системі N = 64, N2 = 256, N3 = 1024 практично співпадають один з одним на інтервалі зміни Кнадм від 0 до 0,834.
З рис. 2 видно, що під час Кнадм ~ 0,80 всі графіки нормованого наробітку до відмови резервованих систем для п’яти різних законів розподілу часу безвідмовної роботи
елементів практично співпадають: — (Кнадм = 0,80) = 1,427 - 1,600.
'0Е
174
ISSN 1028-9763. Математичні машини і системи, 2015, № 1
Відношення середнього наробітку до відмови резервованої системи (N=1024 елементи) до середнього наробітку до відмови елементів системи для різних законів розподілу безвідмовної роботи елементів, Toe ! Т0Е
12
10
8
LD LD Ш
ГПГПГч-,3-Гч-Гч-Гч|і£>СПт
ОіЛО'ГОШО'Г'І^ШСП
гч-гч-гч-сосососпспспсп
6
4
2
0
----ДНР
----УПР
суміш (ЕР+РВ)
Коефіцієнт надмірності, Кнадмір
Рис. 2. Графіки залежності відношення середнього наробітку до відмови резервованої системи (під час загального резервування із дробовою кратністю) до середнього наробітку до відмови елементів системи для різних законів розподілу часу безвідмовної роботи елементів (ЕР, РВ, ДНР, УПР, суміш розподілів експоненціального та Вейбула (ЕР та РВ) від коефіцієнта надмірності
системи (Кнадм, від 0 до 0,997)
На інтервалі зміни Кнадм від 0,80 до 0,997 графіки для відносного середнього наробітку до відмови резервованих систем за різними законами розподілу часу безвідмовної
роботи елементів взагалі не співпадають між собою. Нижче всіх проходить графік для —
ТоЕ
резервованої системи, в якій розподіл часу безвідмовної роботи елементів визначається законом Вейбула ( — {Кнадм = 0,997 ) = 2,97). Вище всіх проходить графік для —
Т0 Е ' Т0е
резервованої системи, в якій розподіл часу безвідмовної роботи елементів визначається сумішшю розподілів експоненціального та Вейбула ( — {Кнадм = 0,997 ) = 9,98 ).
Т0Е
Відношення середнього наробітку до відмови резервованої системи (N = 1024 елементи) до середнього наробітку до відмови елементів системи,
Т0с ! Т0Е
1
0,9
0,8
0,7
0,6
0,5
0,4
0,3
0,2
0,1
0
ІЛ ІЛ ІЛ ГЧІ О Г»-
Q o' Q
U")
о'
Коефіцієнт надмірності, Кнадмір.
----ЕР
----РВ
----ДНР
----УПР
Рис. 3. Графіки залежності відношення середнього наробітку до відмови резервованої системи (під час загального резервування із дробовою кратністю) до середнього наробітку до відмови елементів системи (для різних законів розподілу часу безвідмовної роботи елементів) (ЕР, РВ, ДНР, УПР, суміш ЕР та РВ) від коефіцієнта надмірності системи (Кнадм. від 0 до 0,5)
ISSN 1028-9763. Математичні машини і системи, 2015, № 1
175
Найбільшу цікавість для інженерних додатків мають резервовані системи з загальним резервуванням із дробовою кратністю, які мають коефіцієнт надмірності Кнадм у діапазоні від 0 до 0,25 (рис. 3).
З рис. 3 видно, що криві для ЕР, УПР та РВ розташовані паралельно під одним кутом на осі абсцис ( Кнадм ). Крива для ЕР розташована нижче кривої для УПР та РВ, а крива для УПР розташована нижче кривої для РВ. Криві для ДНР та суміші (ЕР та РВ) є практично паралельними і розташовані під меншим кутом до осі абсцис ( Кнадм) і перетинають криві для ЕР, УПР та РВ.
Відношення середнього наробітку до відмови резервованої системи для неекспоненці-ального закону розподілу часу безвідмовної роботи елементів системи середнього наробітку до відмови резервованої системи для ЕР
S)V O'
----ДНР
----УПР
----Суміш (ЕР+РВ)
Коефіцієнт надмірності, Кнадім
Рис. 4. Графіки кривих залежності відношення середнього наробітку до відмови резервованої системи (під час загального резервування з дробовою кратністю) для різних неекспоненціальних законів розподілу часу безвідмовної роботи елементів системи (РВ, ДНР, УПР, суміш (ЕР+РВ) до середнього наробітку до відмови резервованої системи для експоненціального закону розподілу часу безвідмовної роботи елементів системи від коефіцієнта надмірності
Із графіків на рис. 2 і 3 видно, що крива середнього наробітку до відмови резерво-ваних систем із загальним резервуванням та дробовою кратністю під час неекспоненціаль-них розподілів часу безвідмовної роботи елементів системи значно відрізняється від середнього наробітку до відмови резервованих систем із експоненціальних розподілів часу безвідмовної роботи. На рис. 4 показані графіки кривих відношення середнього наробітку до відмови резервованих систем із дробовою кратністю під час неекспоненціальних розподілів часу безвідмовної роботи елементів системи до середнього наробітку до відмови ре-зервованих систем під час експоненціальних розподілів часу безвідмовної роботи за зміною коефіцієнта надмірності від 0 до 0,25. Графіки кривих мають гіперболічний характер із максимальними значеннями Кнадм = 0. Найбільше значення під час Кнадм = 0 має крива для резервованих систем з часом безвідмовної роботі елементів, які розподілені відповідно дифузійному немонотонному розподілу:
ТТ'ДНР (Кн„„м. = 0) = 81,6.
Т0 ЕР
Графік кривої для резервованих систем з часом безвідмовної роботи елементів, розподілених відповідно закону Вейбула, розташовується нижче (по осі ординат):
Т0РВ
^ СКнадм. = 0) = 36,1.
Т 0 ЕР
176
ISSN 1028-9763. Математичні машини і системи, 2015, № 1
Ще нижче по осі ординат розташовуються графіки кривих для резервованих систем з часом безвідмовної роботи елементів, розподілених відповідно закону для узагальненого показового розподілу:
Т4ПР (К„„„,. = 0) = 3,79
Т0 ЕР
та відповідно закону для суміші (ЕР та РВ):
T°<ff*FB) (К„„а,. = 0) = 4,20.
Т 0 ЕР
5. Висновки
У статті було досліджено вплив коефіцієнта надмірності на показники надійності резервованих систем з дробовою кратністю за різними законами розподілу часу безвідмовної роботи елементів системи: експоненціальним, Вейбула, дифузійним немонотонним, узагальненим показовим та суміші розподілів (експоненціальним і Вейбула). Були побудовані графіки залежностей відношення середнього наробітку до відмови резервованої системи до середнього наробітку до відмови елементів системи для різних законів розподілу часу безвідмовної роботи від коефіцієнта надмірності системи.
Отримані у цій статті результати для резервованих систем з загальним резервуванням із дробовою кратністю можуть знайти широке застосування у технічних додатках під час досліджень надійності фазованих антенних решіток радіолокаційних і гідроакустичних станцій із будь-якою надмірністю для елементів розподільної структури, час безвідмовної роботи яких характеризується неекспоненціальними законами розподілу.
СПИСОК ЛІТЕРАТУРИ
1. Половко А.М. Основы теории надежности: практикум / А.М. Половко, С.В. Гуров. - СПб.: БХВ-Петербург, 2006. - 560 с.
2. Надежность технических систем: справочник / Под ред. И.А. Ушакова. - М.: Радио и связь, 1985. - 608 c.
3. Надежность и эффективность в технике: справочник в 10 т. / Под ред. В.И. Патрушева,
A. Ю. Рембезы. - Москва: Машиностроение, 1990. - Т. 5: Проектный анализ надежности. - 316 с.
4. Азарсков В.Н. Надежность систем управления и автоматики: учебн. пособ. / В.Н. Азарсков,
B. П. Стрельников. - К.: НАУ, 2002. - 164 с.
5. Игнатов В.А. Статистическая оптимизация качества функционирования электронных систем / Игнатов В.А., Маньшин Г.Г., Трайнев В.А.; под ред. Е.Г. Коновалова. - М.: Энергия, 1974. - 264 с.
6. Костановський В.В. Математичні моделі надійності типових апертур фазованих антенних решіток, які враховують раптові та поступові відмови модулів надвисоких частот / В.В. Костановський // Математичні машини і системи. - 2014. - № 2. - С. 142 - 150.
Стаття надійшла до редакції 12.11.2014
ISSN 1028-9763. Математичні машини і системи, 2015, № 1
177