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

ЯКІСТЬ, надійність і сертифікація обчислювальної

ТЕХНІКИ І ПРОГРАМНОГО ЗАБЕЗПЕЧЕННЯ

УДК 621.3.019:658.562 П.В. СТРЕЛЬНИКОВ

ЭКСПЕРИМЕНТАЛЬНАЯ ОЦЕНКА НАДЕЖНОСТИ ИЗДЕЛИЙ В УСЛОВИЯХ МАЛОГО ЧИСЛА ОТКАЗОВ

Анотація. Представлено результати аналізу існуючих методів і методик експериментальної оцінки показників надійності. Розглядаються методики оцінки показників надійності в умовах малої статистики відмов на основі використання двопараметричних імовірнісно-фізичних моделей надійності (дифузійних розподілів).

Ключові слова: надійність, відмова, експериментальна оцінка показників надійності, статистика відмов.

Аннотация. Представлены результаты анализа существующих методов и методик экспериментальной оценки показателей надежности. Рассматриваются методики оценки показателей надежности в условиях малой статистики отказов на основе использования двухпараметрических вероятностно-физических моделей надежности (диффузионныхраспределений).

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

Abstract. The results of the analysis of existing methods and techniques of experimental reliability evaluation are presented. Methodologies for the assessment of reliability in low failure statistics based on the use of two-parameter probabilistic-physical models of reliability (diffusion distributions) are considered. Key words: reliability, failure, experimental evaluation of reliability, failure statistics.

1. Введение

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

2. Оценивание показателей надежности изделий в условиях малого числа отказов

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

оценивания показателей надежности на основе использования диффузионного немонотонного распределения ( DN -распределения).

Использование двухпараметрических вероятностно-физических моделей отказов позволяет, определив априори параметр формы диффузионных распределений, преобразовать двухпараметрическую функцию распределения в однопараметрическую функцию,

которая при этом имеет типичную форму закона распределения наработки до отказа - одномодальную кривую с левосторонним смещением моды. На рис. 1 показаны кривые плотностей однопараметрического экспоненциального распределения и однопараметрического DN -распределения. Как видно из приведенных на рисунке плотностей распределения наработки до отказа (принято условно математическое ожидание наработки до отказа равным единице), расхождение между экспоненциальным законом распределения и DN -распределением весьма существенно. При этом очевидно, что закон распределения наработки до отказа для наиболее распространенного коэффициента вариации наработки (кривая f3) представляется наиболее реальным. Многочисленные исследования многих авторов подтверждают этот вывод. Установлено, что использование однопараметрического экспоненциального распределения приводит к большим погрешностям в оценке искомых показателей надежности, а также к существенному увеличению объема испытаний в 1,5 и более раз по сравнению с использованием более адекватных двухпараметрических функций распределения наработки до отказа.

Благодаря многолетнему опыту эксплуатации и испытаний, значения коэффициентов вариации наработки до отказа (ресурса) для типовых объектов электронного, электротехнического и механического оборудования априори определяют коэффициент вариации наработки на основании анализа типовых процессов деградации, приводящих к отказам. Значения коэффициентов вариации наработок и видов процессов разрушений для типовых объектов приведены в ряде стандартов [2-4]. Учитывая объем и число рассмотренных данных (выборок), на основании которых определены диапазоны коэффициентов вариации (десятки выборок), можно утверждать, что доверительная вероятность указанных интервалов существенно выше 0,9.

Принятие численного значения коэффициента вариации из указанных диапазонов выше приведенных в стандартах в каждом конкретном случае диктуется соображениями общего характера: увеличение отношения нагрузки к пределу выносливости (прочности) относительно среднего статистического приводит к уменьшению коэффициента вариации и наоборот, то есть, чем меньше коэффициент нагружения, тем больше коэффициент вариации.

Оценка показателей надежности при планах испытаний (наблюдений) [NUr] и [NUT], приводящих к однократному цензурированию статистических данных об отказах При малом числе отказов (r < 6, d < 6) значения параметра формы (i~,v,v) DN -распределения определяют согласно вышеприведенным рекомендациям.

Рис. 1. Графики плотностей распределения наработки до отказа (А -экспоненциальное распределение; £2, £3 - DN -соответственно для коэффициента вариации, равного 1 и 0,7)

Определение параметра масштаба /Л БЫ -распределения осуществляется следующим образом.

По результатам наблюдений формируют вариационный числовой ряд по неубыванию суммарных наработок до отказа (предельного состояния): г1, г2,..., гг (г1, г2,..., г й) . Вычисляют значения эмпирической функции распределения наработки в каждый момент г вариационного ряда (] = 1, 2,..., т; т = г или т = ё) по формуле

F, = j / N .

(1)

Вычисляют точечную оценку параметра т и доверительные границы параметра масштаба (л, Л) по формулам

1

m=-Е t—

т j=i

xl —; v

VN’ .

f

m = mx

1 - q;

V

m = mx

q;

V

4m

(2)

где m = r (r - число отказов при плане испытаний [NUr]) или т = d (d - число отказов при плане испытаний [NUT]).

Значение величины x (F; v), представляющей собой относительную наработку при DN -распределении для вероятности отказа F при коэффициенте вариации наработки V, определяют по соответствующим таблицам DN -распределения или решая следующее уравнение относительно x:

Ф

+ exPl -т-1'Ф

x j

у

x +1

■F.

где ф(-) - функция нормированного нормального распределения.

Определив оценки параметров (л, Л, Л, П, ~, Л), вычисляют точечные оценки и доверительные границы соответствующих показателей надежности по формулам табл. 1. Таблица 1. Формулы для вычисления оценок показателей надежности при БЫ -

Оценки Средняя наработка (ресурс) Г амма-процентная наработка до отказа(ресурс) Вероятность безотказной работы за наработку t

Точечная оценка m m x (1 -y;v) Ф ^ v л m-1 V ^4^ у -exp(il ]Ф m+1 ^ V y4n~tj

НДГ уровня q m mx (1 -g; v) Ф С , 1 m-1 V № j •eXP^V2 )Ф m+1 ^ ( v № j

ВДГ уровня q m mx (1 -g; v) Ф ' m-t ) (21Ф 1— 1 -exP 2 Ф Vv4mtJ Vv J r m+11 V vjpt J

Оценивание показателей надежности при планах испытаний [ЫИг] и [ЫЯТ], приводящих к статистическим данным с многократным цензурированием

Параметр формы (коэффициент вариации наработки) определяют аналогично предыдущему случаю.

Определение параметра масштаба т БЫ -распределения.

По результатам наблюдений формируют вариационный числовой ряд по неубыванию суммарных наработок до отказа и до цензурирования: (ї1, їїг, т1, т2,..., тп) . Определяют значение эмпирической функции распределения в каждый момент ї. (т.) вариационного ряда по формуле

р. = Р- +(1 - ^м) г.

7-1

N + г _ X (Г + П)

¿=1

-1

(3)

где г

7(г)’ П7(¿)

- соответственно число отказов (полных наработок) и число неполных наработок в ] (/') интервале (j = 1,2,..., w), г - число замененных образцов. Число интервалов и их граничные значения принимают из соображений удобства расчета.

Вычисляют точечную оценку параметра масштаба р в общем случае по формуле

т

X к.

7=1

X кЛ /N; ~)]1 ’

7=1

где kj - число совпадающих наработок на j -ом интервале.

Доверительные границы параметра масштаба (р, р) вычисляют по формулам (2).

Точечные оценки и доверительные границы соответствующих показателей надежности вычисляют согласно табл. 1.

Оценка показателей надежности при отсутствии отказов (малая выборка)

В теории биномиальной схемы испытаний с остановкой в случае безотказных испытаний установлено (в частности, Клоппером и Пирсоном [5]), что значение нижней доверительной границы отсутствия отказа при испытании N образцов может быть оценено выражением

Р = (1 -¡)1/ N, (4)

где Ь - односторонняя доверительная вероятность оценки вероятности безотказной работы системы (совокупности N образцов).

Связь между односторонней доверительной вероятностью Ь и двусторонней q доверительными вероятностями: ¡5 = (1 + q) / 2.

В случае, когда на момент контроля (цензурирования) группа идентичных образцов (N > 4) не имела отказов, параметр формы V ( V , у ) определяют априори согласно известным рекомендациям.

Нижнюю границу вероятности отсутствия отказа за интервал испытаний (эксплуатации) ?и , если на испытании (эксплуатации) находилось N образцов и отказ не зафиксирован, вычисляют по формуле

\1/N

1 - q 1

Р(К) =

2

Вычисляют нижнюю доверительную границу параметра масштаба р DN -распределения, решая уравнение

Р(!и) = Ф

т_ їи ( 2 т + їи

1 П^т?и і _ ехр1 ^ ] ф і

1

Используя таблицы DN -распределения, по которым установленным способом (по значениям Е = 1 -Р(1и) и V = V) определяют х[1 -Р(1и),у], выражение для оценки р можно записать в следующем виде:

р = —,----------л = 1и ■ К(Р,V),

- х[1 - Р(1и); у] " '•

где К* (Р, V) - поправочный коэффициент, учитывающий эмпирическую вероятность отсутствия отказа; х[] - относительная наработка по DN -распределению, определение которой указано выше.

Используя оценку р, соответствующую доверительной вероятности q, получают выборочную среднюю оценку параметра р, :

р = (,р ~) = !„ ■ К*(Р, V)■ К2~),

X ( 1 - ^ V )

где К2 ^, ~) - поправочный коэффициент, учитывающий вид распределения и доверительную вероятность оценки параметра.

Оценку верхней доверительной границы параметра масштаба р вычисляют по формуле

р=р х (q; ~ ) = р К 2* (^ ~).

Определив оценки параметров (р, р, р, V , ~, р), вычисляют точечные оценки и доверительные границы соответствующих показателей надежности согласно табл. 1.

3. Заключение

Применение в качестве теоретической модели диффузионных распределений за счет использования дополнительной априорной информации (об ожидаемом значении коэффициента вариации распределения наработки) позволяет существенно сократить объем испытаний.

Предложенные методики на основе вероятностно-физических моделей отказов (диффузионных распределений), в условиях малой статистики отказов, на основе использования дополнительной априорной информации о коэффициенте вариации распределения наработки до отказа позволяют оценить параметры закона распределения наработки до отказа, что дает возможность определить все необходимые показатели надежности исследуемых изделий (среднюю наработку, гамма-процентную наработку, вероятность безотказной работы за заданную наработку, остаточный ресурс и др.). То есть позволяют решить задачи оценки надежности в условиях, когда двухпараметрические, строго вероятностные, модели отказов (Вейбулла, логарифмически нормальное и др.) не работают ввиду малой статистики отказов, необходимой для оценки двух параметров.

СПИСОК ЛИТЕРАТУРЫ

1. Стрельников В.П. Оценка и прогнозирование надежности электронных элементов и систем / В.П. Стрельников, А.В. Федухин. - К.: Логос, 2002. - 486 с.

2. ГОСТ 27.005-97. Надежность в технике. Модели отказов. Основные положения. - Введ. 01.01.99.

- 43 с.

3. ГОСТ 27.506-2000. Надежность в технике. Планы испытаний для контроля средней наработки до отказа (на отказ). Ч.: 2. Диффузионное распределение. - Введ. 01.07.2001. - 36 с.

4. ДСТУ 2862-94. Надежность техники. Методы расчета показателей надежности. Общие требования. - Введ.01.01.96. - 39 с.

5. Надежность и эффективность в технике: справочник в 10 т. / Под ред. Р.С. Судакова, О.И. Тес-кина. - М.: Машиностроение, 1989. - Т. 6: Экспериментальная отработка и испытания. - 376 с.

Стаття надійшла до редакції 26.01.2011