CC BY
35
ЯКІСТЬ, надійність і сертифікація ОБЧИСЛЮВАЛЬНОЇ ТЕХНІКИ І ПРОГРАМНОГО ЗАБЕЗПЕЧЕННЯ
УДК 621.3.019.3 П.В. СТРЕЛЬНИКОВ*
РАЗРАБОТКА ПРОГРАММНОГО ОБЕСПЕЧЕНИЯ ОЦЕНИВАНИЯ ПОКАЗАТЕЛЕЙ НАДЕЖНОСТИ ПРИ МАЛОМ ЧИСЛЕ ОТКАЗОВ
Институт проблем математических машин и систем НАН Украины, Киев, Украина
Анотація. Розроблено програмне забезпечення для обробки результатів спостереження у процесі експлуатації технічних систем з метою визначення показників надійності спостережуваних систем при малій кількості відмов.
Ключові слова: надійність, гамма-відсоткове напрацювання, розподіл наробітку до відмови.
Аннотация. Разработано программное обеспечение для обработки результатов наблюдения в процессе эксплуатации технических систем с целью определения показателей надежности наблюдаемых систем при малом числе отказов.
Ключевые слова: надежность, гамма-процентная наработка, распределение наработки до отказа.
Abstract. The software for processing the results of observations during the operation of technical systems in order to determine the reliability of the observed systems with a small number of failures was developed.
Keywords: reliability, gamma-percent operating time, distribution of operating time.
1. Введение
В ряде работ [1-3] разработаны методики оценки показателей надежности при малом числе отказов при наблюдении в процессе эксплуатации объектов на основе использования в качестве теоретической модели отказов DN -распределения. Типичными планами наблюдений с целью оценки показателей надежности высоконадежных и сложных систем представляются планы с однократным цензурированием (планы [NUr], [NUT]) или планы при многократном цензурировании ([NRr], [NRT]). Следует отметить, что оценки параметров двухпараметрического DN -распределения, в частности, максимально правдоподобные оценки параметра масштаба ц и параметра формы n при усеченных наблюдениях получают из решения достаточно сложных систем уравнений [1]. Ниже представляются программы для вычисления показателей надежности при наблюдении малого числа отказов (но не менее шести).
2. Общие положения
Настоящая программа предназначена для вычисления параметров DN -распределения (m, n) и показателей надежности, когда объем статистических данных об отказах превышает шесть. При числе отказов, равном шести, и при использовании DN -распределения представляется возможным оценить параметры распределения (m, n) с минимальными требованиями по достоверности (q = 0,8; d = 0,4) .
Заметим, что уже много десятков лет уровень доверительной вероятности оценки интенсивности отказов высоконадежных изделий электронной техники (резисторов, интегральных схем и др.) принимается равным q = 0,6.
Программа позволяет получить оценки параметров DN -распределения, вычислить все необходимые показатели надежности при однократном цензурировании и наличии отказов в соответствии с планами испытаний [NUr] и [NUT], а также при многократном цензурировании в соответствии с планами испытаний [NRr] и [NRT]. Все формулы и системы уравнений взяты из работы [ 1].
3. Алгоритм программы
A0. Ввод исходных данных.
A01. План испытаний [NOT]: N=; T=; m = d = ; t1 =; t2 = ;.
A02. План испытаний [NUr]: N=; m = r =; t1 =; t2 = ;.tr
A03. План испытаний [NRr]: N=; r =; t1 =; t2 =;.tr =; n =; t =; t2 =;... tn = (ti - нара-
ботки неотказавших образцов). На A3 .
A„4. План испытаний [NRT]: N=; T=; r = ; t1 =; t2 =;... tr =; n =; t =; t2 =;... tn = . На A3. A2 . Вычисление максимально-правдоподобных оценок параметров DN -распределений при однократном усечении выборки (планы испытаний [NUr] или [NUT]).
Максимально-правдоподобные оценки параметров DN -распределения вычисляют при решении системы уравнений:
m, , (N - m) Г 2 ч m f
v ' n2 dn (tr; mn)+—-mm-
m
. td ; tr td . На A2 .
. На A2 .
t =; t =;... t
n +--mm, r ,
m m [1 - DN (tr; m,n)]
= 0,
2 m, ,
2 + n-----1 — m m_
m
(N - m)
m[1 - DN (tr; m,n)]
(2 + n2) DN (tr; m,n)-
mm-1
m
= 0,
где m = r - число отказов до момента усечения за время испытаний (наблюдений) по плану [NUT] или m = d числу отказов при плане испытаний (наблюдений) [NUr] (при этом
tr = td );
1=—Е t>; m-1 =—Е 7; m-1 = j1 fN (t) dt,
m i=1 m г=! tt 01
•1
lr
mi = j tfN (t) dt; fN (t)
vtJPt
exp
(t -m)2
2v2mt
DN(tr; m,n) = F
n
n
n
На печать //, V . На Д4.
Д3. Вычисление максимально-правдоподобных оценок параметров БЫ -распределения при многократном усечении выборки (планы испытаний [N^.1"] или [N^1]).
Максимально-правдоподобные оценки параметров БЫ -распределения вычисляют при решении системы уравнений:
0
N—r
2 ( mi ............./ 1
m
V2DN(tj ;m,v)
T m,V)+m -mm-ij
1 - DN(t; m,v)
= 0,
2 m' , 1N-
2+v —1—m tn, — у m r1=i
(2 + v2 )DN (tj; m,n)—
mm
—ij
1 — DN tj; m,n)
m
где обозначения m', m—1, ¡u'l ,■, m— ,■ приведены выше;
DN (т; m,n) = F
*1 j ’ ^—1 j
т, —m
+ exp(^ IФ
V
т, +m
. На печать m, v • На A4.
A4. Вычисление доверительных границ параметров DN -распределения.
При достаточной статистике (m > 6) и планах испытаний [NUr] или [NUT] нижние доверительные границы (НДГ) параметров вычисляют по формулам:
(
m=mx
1—q;
v
; v = v<!1 +
(1 + 2~2) Uq — Uq_ 4m 4m
^ [8m + (1 + 2~2)Uq](1 + 2v2) I.
Верхние доверительные границы (ВДГ) параметров БЫ -распределения вычисляют по формулам:
(
m = mx
V
q ;^= I; V = ~ ^1 + (1 + 2 )Uq + [8m + (1 + 2v2)uq ](1 + 2~2) j>.
4m
На печать т, Ц, V , V .
Д5. Вычисление точечных оценок и доверительных границ определяемых показателей надежности на основе БЫ -распределения.
Точечную оценку и доверительные границы средней наработки до отказа (ресурса) вычисляют по формулам:
ср L_ ’ ср
Точечную оценку и доверительные границы гамма-процентной наработки до отказа (ресурса) вычисляют по формулам:
Т7=П ■ х(1 -у; ~); Т7=т-х (1 -у V); Тг = т- х(1 -г; у).
Точечную оценку и доверительные границы вероятности безотказной работы (ВБР) за наработку t вычисляют по формулам:
P(t) = Ф
m—t
V
— exp
2
,v‘
Ф
m+1
V
Jut
; P(t) = ф
( , ) m—t
V
V Vi—
P (t) = Ф
(1~L ) vmt
— exp
( о ) (
Ф
V V J
m+1 yjpt
T4mt
Точечную оценку и доверительные границы вероятности безотказной работы в заданном интервале (¿, ^ + Хзад) вычисляют по формулам:
Ф
m - t - tзад
P (t, t зад ) =
_____V
~V m (t + tзад ) у
- exp( V2
Ф
m+1 +t з
Ф
m-t
V
Jjut
Ф
Ф
P(t, t зад ) = V V~
Ф
f * * m-t - t зад
И Vm (t + t^ad )
V
Í * л
m-1
V
W (2T
Ф
Ф
P (t, t3ad )
m t tзад
n -m+u у
V-
Л f 2 Л
- exp 2 ф
у Чи у
m(t + t зад ) у
m+1 +t зад
m t у m+1
nVm(t + t зад )
N m+1 \
mt y m +1+ / t зад
V —
f +
зад / у
Ф
í - , л
m -1
f 2 )J ”
- exp — Ф
ЧИ у
13).
На печать Tcp, Tcp, Tcp, Tr, Tr, Tg, P(t), P(t), P (t).
Программа была скомпилирована в среде Embarcadero Rad Studio 2009 (C++ Builder
2
2
4. Инструкция, руководство пользователя и пример решения задач с использованием программы n_mater
Программа, реализующая вышеприведенный алгоритм, названа Normat Mater (n_mater). Основная оболочка программы показана на рис. 1.
Инструкция пользователя программы n_mater
1. Запуск программы.
2. Ввод исходных данных в панели «Admission figures».
3. Ввод исходных значений в массив на панели «Arrays data».
4. Запуск программы на выполнение с помощью кнопки «Calculate» в самом низу рабочей формы.
5. Просмотр и анализ результатов в таблице на форме «Normat mater results», вывод на печать с помощью кнопки «Print results».
Примечание. Приняты обозначения: eps - точность вычисления параметров;
N=T count (число образцов, поставленных на испытания); n=Tau count (число цензурирований); T=TNUT (продолжительность испытаний при плане [NUT]); T= NRT (продолжительность испытаний при плане [NRT]); d=dNUT (число отказов при плане [NUT]); r=rNUr (число отказов при плане [NUr]); r=rNRT (число отказов при плане [NRT]); {t1, t2,...,tN} =Array t (count N) = 478, 607, . . . 9250;
{t, t2,..., tn} =Array tau (count n)=.... Если n = 0 (нет наработок до цензурирования), то задают: Tau count =1, а Array tau (count 1) и Value=0. t для P(t, t3ad ) -t ; tзад =tgiv;
q = q (доверительная вероятность оценки параметров);
Пример 1 (план [NUr], DN-распределение)
1. Обнуление старых данных (нажатием кнопки «Clear»).
2. Ввод новых исходных данных: Plan of tests (план испытаний) - NUr.
Ввод исходных данных: T count-32; Array t-478, 607, 770, 860, 990, 1040, 1106, 1180, 1240, 1360, 1441, 1525; Tau count-1; TNUT - 0; TNRT-0; dNUT-0; rNUr-12; rNRr-0; rNRT-0; q-0,9;
3. Если ставится задача оценки таких показателей, как
Tg , P(tзад ), P(t, t + tзад ), то задаются значения: Gamma-0,9; tgiv-200; t-300.
Форма с введенными данными выглядит, как показано на рис. 1.
4. Запуск программы (нажатие S ).
5. Распечатка результатов расчета (нажатие кнопки с пиктограммой «принтер»).
6. Выход из программы (пункт меню Exit или нажатием на крестик в шапке окна).
7. Форма с выведенными данными выглядит, как показано на рис. 2.
g=Gamma.
Программа nmater
Рис. 1. Основная оболочка программы
Рис. 2. Форма с введенными данными
Отметим, что оценка средней наработки до отказа в приведенном примере на основании укороченной выборки (12 наработок), равная Тср = 2373 ч, имеет погрешность менее 5% по сравнению с оценкой средней наработки генеральной совокупности испытанных объектов (магнитофонов) [3] в количестве 96 образцов и равной Т = 2280 ч.
5. Заключение
Разработано программное обеспечение для обработки результатов наблюдения в процессе эксплуатации технических систем с целью определения показателей надежности наблюдаемых систем при малом числе отказов. При этом определяются закон распределения наработки до отказа (на отказ) и все необходимые показатели надежности (средняя наработка до отказа (на отказ), гамма-процентная наработка, вероятность безотказной работы за заданное время или в заданном интервале наработки).
СПИСОК ЛИТЕРАТУРЫ
1. Стрельников В.П. Оценка и прогнозирование надежности электронных элементов и систем / В.П. Стрельников, А.В. Федухин. - К.: Логос, 2002. - 486 с.
2. Стрельников В.П. Оценка надежности по малой выборке с использованием дополнительной априорной информации / В.П. Стрельников, П.В. Стрельников // Материалы Междунар. конф. «Системные проблемы качества, математического моделирования, информационных, электронных и лазерных технологий». - Москва: Радио и связь, 2002. - Ч. 2. - С. 163 - 167.
3. Азарсков В.Н. Надежность систем управления и автоматики: учеб. пособ. / В.Н. Азарсков, В.П. Стрельников. - К.: НАУ, 2004. - 164 с.
Стаття над1йшла до редакцп 08.10.2013
