Алгоритм обработки результатов и расчета показателей надежности

Постановка задачи

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

По результатам испытаний N невосстанавливаемых одинаковых объектов получена статистическая выборка – массив наработки (в любых единицах измерения) до отказа каждого из N испытывавшихся объектов. Выборка характеризует случайную величину наработки до отказа объекта T = {t}.

Необходимо выбрать закон распределения случайной величины T и проверить правильность выбора по соответствующему критерию.

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

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

2.2.1. Формирование статистического ряда

При большом числе испытываемых объектов полученный массив наработок {…, t_i, …} является громоздкой и мало наглядной формой записи случайной величины T.

Поэтому для компактности и наглядности выборка представляется в графическом изображении статистического ряда – гистограмме наработки до отказа.

Для этого необходимо:

- установить интервал наработки [t_min, t_max] и его длину

,

где

, ;

- разбить интервал наработки [t_min, t_max] на k интервалов равной ширины Δt – шаг гистограммы

, .

- подсчитать частоты появления отказов во всех k интервалах

,

где Δn(t_i , t_i + Δt) – число объектов, отказавших в интервале [t_i , t_i + Δt].

Очевидно, что

,

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

Рис. 5.2. Гистограмма.

2.2.2. Расчет эмпирических функций

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

- функция распределения отказов (оценка ВО)

;

 

;

 

;

 

;

 

- функция надежности (оценка ВБР)

;

 

;

 

Рис. 5.3. Графики ВО и ВБР.

 

- плотность распределения отказов (оценка ПРО)

;

 

Рис. 5.4. График ПРО.

 

- интенсивность отказов (оценка ИО)

.

 

Рис. 5.5. Графики ИО.

На рис. 5.3, 5.4, 5.5 приведены соответственно графики статистических оценок , , .

Правила построения графиков ясны из приведенных выше расчетных формул. Каждый из графиков имеет свой масштаб.

2.2.3. Расчет статистических оценок числовых характеристик

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

Оценки характеристик определяются:

- оценка средней наработки до отказа (статистическое среднее наработки):

;

 

- оценка дисперсии наработки до отказа (эмпирическая дисперсия наработки):

;

 

где – середина i-го интервала наработки, т. е. среднее значение наработки в интервале.

2.2.4. Оценка СКО .

Целесообразно рассчитать оценки и некоторых вспомогательных характеристик рассеивания случайной величины T:

- выборочный коэффициент асимметрии наработки до отказа

;

 

- выборочный эксцесс наработки до отказа

;

 

Эти характеристики используются для выбора аппроксимирующей функции.

Так коэффициент асимметрии является характеристикой «скошенности» распределения,

например, если распределение симметрично относительно МО, то A = 0.

На рис. 5.6, а распределение f₂(t) имеет положительную асимметрию A > 0, а f₃(t) – отрицательную A < 0.

Эксцесс характеризует «крутость» (остро- или плосковершинность) распределения.

Для нормального распределения E = 0.

Кривые f(t), более островершинные по сравнению с нормальной, имеют
E > 0, а наоборот – более плосковершинные, E < 0 (рис. 5.6, б).

а)	б)
Рис. 5.4. График ПРО.

2.2.5. Выбор закона распределения.

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

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

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

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

,

где a, b, c, … – неизвестные параметры распределения.

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

На графике вместе с строится теоретическая ПРО f(t), что позволяет визуально оценить результаты аппроксимации (расхождения между и f(t). Поскольку эти расхождения неизбежны, то возникает вопрос: объясняются ли они случайными обстоятельствами, связанными с тем, что теоретическое распределение выбрано ошибочным? Ответ на этот вопрос дает расчет критерия согласия.

2.2.6. Расчет критерия согласия.

Критерий согласия – это критерий проверки гипотезы о том, что случайная величина T, представленная своей выборкой, имеет распределение предполагаемого типа.

Проверка состоит в следующем. Рассчитывается критерий, как некоторая мера расхождения теоретического и эмпирического распределений, причем эта мера является случайной величиной.

Чем больше мера расхождения, тем хуже согласованность эмпирического распределения с теоретическим, т. е. меньше мала, то гипотезу о выборе закона распределения следует отвергнуть, как мало правдоподобную.

В противном случае – экспериментальные данные не противоречат принятому распределению.

Из известных критериев наиболее применяемый критерий согласия χ² (хи-квадрат) Пирсона.

Проверка согласованности распределений по критерию χ² производится следующим образом:

1) рассчитывается критерий χ² (мера расхождения)

,

где – теоретическая частота (вероятность) попадания случайной величины в интервал [t_i, t_i + t];

2) определяется число степеней свободы R = k – L ,

где L – число независимых условий, наложенных на частоты _i ,

например:

а) условие ;

б) условие совпадения ;

в) условие совпадения и т. д.

Чаще всего L = 3. Чем больше число степеней свободы, тем больше случайная величина ² подчиняется распределению Пирсона;

3) по рассчитанным χ² и Rопределяется вероятность P того, что величина, имеющая распределение Пирсона с R степенями свободы, превзойдет рассчитанное значение χ².

Ответ на вопрос: насколько мала должна быть вероятность P, чтобы отбросить гипотезу о выборе того или иного закона распределения – во многом неопределенный.

На практике, если P < 0,1, то рекомендуется подыскать другой закон распределения.

В целом, с помощью критерия согласия, можно опровергнуть выбранную гипотезу, если же Pдостаточно велика, то это не может служить доказательством правильности гипотезы, а указывает лишь на то, что гипотеза не противоречит данным эксперимента.

 

 

Контрольные вопросы:

1. Что представляет математическая модель, и для каких целей она используется в задачах надежности?

2. Из каких условий выбирается закон распределения наработки до отказа объекта?

3. В чем заключается постановка задачи при испытаниях объектов на надежность?

4. Что представляет собой процедура формирования статистического ряда по результатам испытаний?

5. Какие эмпирические функции рассчитываются при обработке результатов испытаний?

6. В чем заключается выбор закона распределения наработки до отказа по результатам испытаний?

7. Что представляет собой критерий согласия?