ОЦЕНКА ПОКАЗАТЕЛЕЙ НАДЕЖНОСТИ ИСПЫТАНИЙ ПРИ ЭКСПОНЕНЦИАЛЬНОМ ЗАКОНЕ РАСПРЕДЕЛЕНИЯ ОТКАЗОВ СРЕДСТВАМИ VBA

​EVALUATION OF TEST RELIABILITY INDICATORS WITH EXPONENTIAL FAILURE DISTRIBUTION LAW BY vba TOOLS

Vladimir Kovalchuk

professor, Department of Mathematics, Information Systems and Software, Murmansk State Technical University,

Russia, Murmansk

Marina Burzun

graduate student, Department of Mathematics, Information Systems and Software, Murmansk State Technical University,

Russia, Murmansk

АННОТАЦИЯ

В статье проведена оценка показателей надежности безотказной работы системы. На примере показан расчет основных показателей средствами VBA.

.

ABSTRACT

The article evaluates the reliability indicators of the system's uptime. The example shows the calculation of the main indicators by means of VBA.

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

Keywords: uptime, confidence interval, tests, exponential distribution law, number of failures.

Проведение испытаний новых изделий или системы организуется в соответствии с планом, в котором указывается: количество испытуемых изделий (N), будут ли заменяться отказавшие изделия (B) и когда испытания необходимо прекратить (прекращение испытаний после истечения Т часов или прекращение испытаний при возникновении r – го отказа).

Экспоненциальный закон представляет собой постоянную интенсивность отказов, т.к. определяется параметром λ- const при δ=1.

, интенсивность отказов λ(t)≡λ.                            (1)

Если указано время работы каждого изделия от начала работы до его отказа, расчет суммарной наработки всех элементов SБ(r) и SБ(Т) вычисляется с использованием выражения (2).

                                                       (2)

Среднее время безотказной работы – T_ср при экспоненциальном законе распределения равно величине, обратной ИО – 1/λ, т.е.:

                                                                   (3)

-плотность распределения наработки              (4)

=∙(-ln).

Вероятность восстановления P_B(t):

Дисперсия времени безотказной работы:

                                                                      (5)

Условная вероятность безотказной работы устройства на интервале времени t, после того как оно безотказно проработало на интервале τ.

                                                      (6)

Под доверительным интервалом понимается диапазон значений параметра, в пределах которого с некоторой вероятностью γ может находиться его истинное значение. Вероятность γ называется доверительной вероятностью или коэффициентом доверия. Для экспоненциального закона распределения отказов при плане испытаний N, Б, r установлено, что величина U = 2S_Б(r)λ подчиняется χ² – распределению с 2r степенями свободы, где S_Б(r) – суммарная наработка изделий, установленных на испытание (может быть определена из выражения (7), λ – истинное значение интенсивности отказа, r – число отказов (или разрядов, если отказы сгруппированы по разрядам) [1].

                                   (7)

Вероятность того, что величина U находится в пределах χ²₁ и χ²₂, равна

Интеграл табулирован. Поэтому, задавшись значениями λ_min и λ_max и зная из обработки результатов эксперимента суммарную наработку S_Б(r), находим χ²₁ и χ²₂ и по таблице квантилей распределения χ² –квадрат находим коэффициент доверия γ.

Квантилем случайной величины χ называется такое значение случайной величины х_р, для которого с вероятностью 1– р можно утверждать, что полученное значение этой случайной величины попадет в интервал (–∞, х_р).

Однако чаще стоит обратная задача: по коэффициенту доверия γ и суммарной наработке изделий при испытании SБ(r) требуется найти λmin и λmax.

Установлено, что доверительный интервал будет минимальным, если площади под кривой плотности распределения f_2r в интервалах [0, χ²₁] и [χ²₂, ∞] равны.

Тогда значения χ²₁ и χ²₂ ограничивают соответственно площади 0,5(1+γ) и 0,5(1 – γ).

Последовательность определения доверительных интервалов сводится к следующему. Задавшись коэффициентом доверия γ, определяем 0,5(1+γ) и 0,5(1 – γ) и, зная число степеней свободы 2r, по таблице квантилей χ²– распределения находим значения χ²₁ и χ²₂. Доверительные оценки λ_min и λ_max могут быть определены из неравенства

                                                                    (8)

Отсюда найдем

                                                                                            (9)

                                                                                            (10)

Рассмотрим пример расчета показателей, полученных при испытании 100 изделий (из строя вышло 34). Испытания были прекращены после истечения 100 часов.

Для построения статистического ряда время испытаний разобьем на равные интервалы (разряды) продолжительностью 10 часов и для каждого разряда проведем расчет.

Поскольку за время испытания отказало 34 % изделий, оценка интенсивности отказов подсчитывалась с использованием выражений для плана N, Б, Т.

Доверительные границы:

0,2

Частота отказов определялась для каждого разряда из выражения

.

fэ, λэ, Qэ(t) – параметры потока отказов при экспоненциальном законе распределения [1, 2].

Вероятность отказа подсчитывалась по формуле

.

Для определения доверительных интервалов при экспоненциальном законе распределения по таблице квантилей χ² – квадрат распределений. найдем  (20) и (20) [3]:

,

.

Листинг программы расчета показателей при экспоненциальном законе распределения:

'Вычислим 1 строку таблицы(3)=============================  t

Строка = 3

For СтолбецТаблицы = 4 To КоличествоСтолбцовТаблицы + 3

Sheets("ОсновнаяТаблица").Cells(Строка, СтолбецТаблицы).Value = Sheets("ОсновнаяТаблица").Cells(Строка - 1, СтолбецТаблицы).Value * ВременнойИнтервал

Next

'Вычислим 2 строку таблицы(4)=============================   n*

n = 0 'количество вышедших из строя элементов в периоде

СтрокаДанных = 4

СтрокаТаблицы = 4

СтолбецТаблицы = 4

КонтрольноеЗначениеВременногоИнтервала = ВременнойИнтервал

While Sheets("Исходные данные").Cells(СтрокаДанных, 1).Value <> ""

 If КонтрольноеЗначениеВременногоИнтервала > ОбщееКоличествоЧасов Then GoTo конец

    If Sheets("Исходные данные").Cells(СтрокаДанных, 1).Value <= КонтрольноеЗначениеВременногоИнтервала Then

        n = n + 1

    Else

        Sheets("ОсновнаяТаблица").Cells(СтрокаТаблицы, СтолбецТаблицы).Value = n

КонтрольноеЗначениеВременногоИнтервала = КонтрольноеЗначениеВременногоИнтервала + ВременнойИнтервал

 СтрокаДанных = СтрокаДанных - 1

 СтолбецТаблицы = СтолбецТаблицы + 1

 n = 0

    End If

    СтрокаДанных = СтрокаДанных + 1

Wend

If n <> 0 Then Sheets("ОсновнаяТаблица").Cells(СтрокаТаблицы, СтолбецТаблицы).Value = n

конец:

'Вычислим 7 строку таблицы(9)=============================ЛямбдаЭ x 10^3

СтрокаТаблицы = 9

СтолбецТаблицы = 4

СуммаВышедшихЗаВсеВремя = 0

For n = СтолбецТаблицы To (КоличествоСтолбцовТаблицы + СтолбецТаблицы - 1)

СуммаВышедшихЗаВсеВремя = СуммаВышедшихЗаВсеВремя + Sheets("ОсновнаяТаблица").Cells(4, n).Value

Next

СтолбецТаблицы = 4

Sheets("ОсновнаяТаблица").Cells(СтрокаТаблицы, СтолбецТаблицы).Value = Round((1 / 100 * Log(КоличествоЭлементов / (КоличествоЭлементов - СуммаВышедшихЗаВсеВремя))) * 1000, 2)

Sheets("ОсновнаяТаблица").Range(Cells(СтрокаТаблицы, 4), Cells(СтрокаТаблицы, n - 1)).MergeCells = True

Sheets("ОсновнаяТаблица").Range(Cells(СтрокаТаблицы, 4), Cells(СтрокаТаблицы, n - 1)).HorizontalAlignment = xlCenter

'Вычислим 8 строку таблицы(10)=============================fэ х 10^3

СтрокаТаблицы = 10

СтолбецТаблицы = 4

For n = СтолбецТаблицы To (КоличествоСтолбцовТаблицы + СтолбецТаблицы - 1)

Sheets("ОсновнаяТаблица").Cells(СтрокаТаблицы, n).Value = Round(((Sheets("ОсновнаяТаблица").Cells(9, СтолбецТаблицы).Value / 1000) * Exp(-Sheets("ОсновнаяТаблица").Cells(2, n).Value / 1000 * Sheets("ОсновнаяТаблица").Cells(9, СтолбецТаблицы).Value * ВременнойИнтервал) * 1000), 2)

Next

'Вычислим 9 строку таблицы(11)=============================Рэ

СтрокаТаблицы = 11

СтолбецТаблицы = 4

For n = СтолбецТаблицы To (КоличествоСтолбцовТаблицы + СтолбецТаблицы - 1)

Sheets("ОсновнаяТаблица").Cells(СтрокаТаблицы, n).Value = Round(Exp(-Sheets("ОсновнаяТаблица").Cells(9, СтолбецТаблицы).Value * Sheets("ОсновнаяТаблица").Cells(3, n).Value / 1000), 2)

Next

'Вычислим 10 строку таблицы(12)=============================Qэ(t)

СтрокаТаблицы = 12

СтолбецТаблицы = 4

For n = СтолбецТаблицы To (КоличествоСтолбцовТаблицы + СтолбецТаблицы - 1)

Sheets("ОсновнаяТаблица").Cells(СтрокаТаблицы, n).Value = 1 - Sheets("ОсновнаяТаблица").Cells(11, n).Value

Next

Результаты вычислений представлены в таблице Excel (Таблица 1).

Таблица 1.

Результаты расчета основных показателей испытаний

Результаты расчета fэ(t) представлены на рис 1.

Рисунок 1. Гистограмма частоты отказов при экспоненциальном законе распределения

Поскольку в задаче задано время работы каждого изделия до отказа, суммарная наработка всех изделий S_Б(T) подсчитывается по формуле (11)

S_Б(T) =                                                        (11)

S_Б(T) = 8276 ч.

Список литературы:

1. Коваленко, В. Н., Новиков, А. А. Надежность устройств железнодорожной автоматики, телемеханики и связи. учеб. пособие. – Екатеринбург: УрГУПС, 1995. – с. 78.
2. Основы теории надежности автоматических систем управления : учеб. пособие для вузов / Л. П. Глазунов, В. П. Грабовецкий, О. В. Щербаков. – Л. : Энергоатомиздат, Ленинградское отд-ние, 1984. – 208 с.
3. Бронштейн И. Н., Семендяев К. А. Справочник по математике для инженеров и учащихся ВТУЗов. – М. : Наука, 1980. – 976 с.