Пиши Дома Нужные Работы

Обратная связь

ТЕХНОЛОГИЧЕСКИЕ МЕТОДЫ ОБЕСПЕЧЕНИЯ НАДЕЖНОСТИ

ТЕОРИЯ НАДЕЖНОСТИ

 

Методические указания к практической работе №1

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

 

 

Составители:

д.т.н., профессор Богданович В.И.

к.т.н., доцент Плотников А.Н.

 

Самара 2013

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

Введение

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

Например, рассеяние ресурсов по критерию усталости, оцениваемое отношением наибольшего ресурса к наименьшему достигает 40. Значительным является рассеяние ресурсов по износу и коррозии. При выполнении любой технологической операции также возникает рассеяние параметров качества изделий.

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

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



При решении второй задачи надежности, в зависимости от конкретных технических систем и целей, используются следующие методы:

1. расчетный;

2. опытно-статистический;

3. регистрационный;

4. экспертный.

1 Теоретические основы

Пусть технический объект начинает работу в момент времени , а в момент времени происходит его отказ.

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

Следовательно, время – наработка до отказа является случайной величиной.

Множество возможных значений , называемое генеральной совокупностью, совпадает с множеством действительных чисел .

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

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

Рассмотрим методику обработки статистических данных для определения закона распределения наработки до отказа.

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

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

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

, (1)

где – число случаев, благоприятных событию A, общее число случаев.

В нашем примере эмпирическая (или статистическая) оценка вероятности (1) (также называемая частотой или частотностью наступления события) будет для всех событий наступления отказов в моменты времени одинакова и равна

, (2)

так число случаев благоприятных отказу в момент равно одному.

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

Таблица 1 – Ряд распределения

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

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

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

. (3)

В рассматриваемом случае, когда случайная величина является наработкой до отказа , а область возможных реализаций представлена (в общем случае) бесконечным временным интервалом , соотношение (3) для ФРВ наработок до отказа ( ФРВ отказов ) должно быть записано в виде:

. (4)

Из определения ФРВ следует, что для таблицы 1 эмпирическая (или статистическая) ФРВ отказов принимает значения

. (5)

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

Рисунок 1 – Функция распределения вероятностей для таблицы 1

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

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

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

. (6)

Причём при этом естественно будет выполняться условие нормировки статистических вероятностей:

. (7)

Полученная часть таблицы 2, состоящая из трёх рассмотренных строк, носит название группированного статистического ряда. Разделив каждое значение на длину интервала , получим значения плотностей распределения статистической вероятности на каждом интервале:

. (8)

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

Таблица 2 – Построение группированного статистического ряда, статистической плотности распределения, функции распределения вероятностей и интенсивности отказов

Интервал разбиения   …   … k

 

Рисунок 2 – Гистограмма плотности распределения вероятности наработок до отказа

Имея в своём распоряжении данные таблицы 2, строим статистическую функцию распределения вероятности отказов, значения которой на интервалах разбиения в соответствии с (4) рассчитываются по соотношению:

. (9)

Значения приводим в пятой строке таблицы 2. При построения гистограммы учитываем, что при значениях равна нулю, а при значениях единице.

В шестой строке таблицы 2 приводим значения интенсивности отказов для каждого интервала разбиения, значения которых определяем по соотношению:

. (10)

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

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

Задача оптимального выбора выравнивающей кривой является в значительной мере неопределенной и не имеет строгого математического решения. Здесь приходится руководствоваться соображениями физики решаемой задачи, характером эмпирической зависимости, оценкой степени точности результатов наблюдений и т.д. Например, гистограмма, показанная на рисунке 2, явно наводит на мысль о нормальном распределении. Монотонное же убывание свидетельствовало бы о показательном распределении, как наиболее вероятном.

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

Величина представляет собой взвешенную сумму квадратов разностей между гипотетическими и наблюденными частностями попадания случайной величины в -й интервал:

. (11)

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

Алгоритм оценки правдоподобия гипотезы по критерию сводится к следующему.

1.1. В дополнение к группированному ряду по гипотетической плотности распределения вычисляются гипотетические интервальные вероятности:

. (12)

1.2 Вычисляется величина

1.3 По таблице распределения находится вероятность (приложение Б).

1.4 Если найденная вероятность достаточно мала (стандартные уровни значимости ; ), то выдвинутую гипотезу следует отвергнуть, как малоправдоподобную. В противном случае принять, как неопровергаемую результатами наблюдений.

Порядок выполнения работы

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

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

: 943; 7400; 10640; 8850; 12780; 7190; 11600; 8640; 5904; 5996; 16820; 6900; 2900; 10240; 10770; 7100; 9630; 5260; 12350; 10150.

2.1 Строим вариационный ряд, располагая значения в порядке возрастания, а результаты заносим в таблицу 3.

Таблица 3

0,943 58,68 8,850 0,06
2,900 32,52 9,630 1,06
5,260 11,16 10,150 2,40
5,904 7,29 10,240 2,69
5,996 6,78 10,640 4,16
6,900 2,89 10,770 4,71
7,100 2,25 11,600 9,00
7,190 1,99 12,350 14,06
7,400 1,44 12,780 17,47
8,640 0,00 16,820 67,57

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

тыс. час. (13)

2.3 Заполняем последний столбец таблицы 3 и вычисляем величину квадрата среднеквадратичного отклонения наработок до отказа (оценка дисперсии) , а извлекая из полученного значения корень квадратный получаем величину статистической оценки стандартного отклонения:

; тыс. час. (14)

2.4 Построение гистограммы.

2.4.1 Определяем количество разрядов по формуле и определяем величину интервала времени для выбранного числа разрядов час. Для удобства дальнейших вычислений в качестве значения выбираем его округленное значение час., а начало первого разряда выбираем не при значении 943 час., а в начальный момент времени. При этом величина интервала времени для (1 – 4) разрядов будет 3000 час., а для последнего 5 разряда этот интервал надо увеличить до 5000 час. чтобы он включал все экспериментальные данные (таблица 4).

2.4.2 Заполняем таблицу 4, которая соответствует ранее рассмотренной таблице 2 с другим расположением строк и столбцов

Таблица 4 – Экспериментальные данные для построения группированного статистического ряда и гистограмм плотности распределения вероятности, функции распределения вероятности и интенсивности отказов

час час час
0÷3 0,100 0,033 0,100 0,037
3÷6 0,150 0,050 0,250 0,067
6÷9 0,300 0,100 0,550 0,220
9÷12 0,300 0,100 0,850 0,667
12÷17 0,150 0,030 1,0 -

При заполнении таблицы 4 , , и вычисляем по формулам (6), (8), (9) и (10): ; ; ; .

2.4.3.Строим три гистограммы, откладывая по оси абсцисс разряды, а по оси ординат значение ,, и .

Рисунок 3 – Гистограммы плотности распределения наработок до отказа, функции распределения вероятности отказов и интенсивности отказов

2.4.4 Построение гипотетического закона распределения.

Исходя из вида полученной гистограммы для можно предположить, что распределение вероятности наработок до отказа подчиняется нормальному закону (усеченное на интервале ) распределения с параметрами :

,

где – нормировочный множитель;

– табличная функция Лапласа.

Вычисляем и по таблице функции Лапласа (Приложение А) находим . Далее вычисляем . Следовательно, без потери точности с погрешностью менее 1 % можно пользоваться стандартным нормальным законом, положив :

; (15)

; (16)

, (17)

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

2.4.5 Для проверки этой гипотезы заполняем таблицу 5, в которую переносим часть данных из таблицы 4.

Гипотетические вероятности вычисляем по формуле:

, (18)

учитывая, что . Значения указаны в таблице Приложения А.

Заполняем столбцы и таблицы 5.

Таблица 5 – Исходные данные для определения величины критерия согласия Пирсона

0÷3 0,100 0,052 0,044
3÷6 0,150 0,175 0,004
6÷9 0,300 0,308 0,000
9÷12 0,300 0,283 0,001
12÷17 0,150 0,169 0,002

2.4.6 Вычисляем Пирсоновскую меру отклонения между гипотетическим и эмпирическим распределениями по формуле:

.

2.4.7 Находим число степеней свободы распределения по формуле , где – число разрядов гистограммы, 3 – число связей, наложенных на величины .

2.4.8 По таблице распределения (Приложение Б) для и находим, что .

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

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

Используем соотношение (17) и получим:

.

Откуда . По таблице Лапласа находим или тыс. час. Поэтому можно назначить ресурс насоса час. с вероятностью безотказной работы не менее 0,99, но не более 0,991.

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

3.1. Почему возникают отказы технических объектов?

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

3.3. Дайте определение понятий средняя наработка до отказа, математическое ожидание, дисперсия, стандартное отклонение и коэффициент вариации.

3.4. Дайте определение понятий функции распределения случайной величины и интенсивность отказов.

3.5. Приведите соотношение для расчёта статистических величин: вероятности, плотности распределения вероятности, функции распределения вероятности и интенсивности отказов случайной величины.

3.6. Как строится функция распределения вероятностей случайной величины?

3.7. Какую характеристику случайной величины отражает гистограмма?

3.8. Плотности распределения вероятности наработок до отказа имеет максимум при часов (рисунок 3).Означает ли это, что при наработке часов изделия отказывают реже?

3.9. В чем заключается геометрический смысл меры отклонения между гипотетическим и эмпирическим законами распределения, определяемой по критерию Пирсона?

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

3.11. Приведите выражение для плотности распределения вероятности, функции распределения вероятности наработок до отказа и функции распределения вероятности безотказной работы для нормального закона распределения вероятности.

4 Контрольные задания

Вариант №
943.263 1.075∙104 1,186∙104 2,972∙103 9,233∙103 2.972∙103 5.619∙103 5,079∙103 5.619∙103
7.403∙103 1.188∙104 8,248∙103 6,201∙103 1.006∙104 6.201∙103 4.416∙103 6,007∙103 4.416∙103
1.064∙104 1.055∙104 1,142∙104 7.822∙103 7.648∙103 7.822∙103 5.474∙103 5,406∙103 5.474∙103
8.847∙103 6.911∙103 1,177∙104 6.923∙103 7.357∙103 6.923∙103 5.589∙103 6,003∙103 5.589∙103
1.278∙104 9.434∙103 6,52∙103 8.889∙103 9.136∙103 8.889∙103 3.84∙103 3,773∙103 3.84∙103
7.186∙103 1.013∙104 1.292∙104 6.093∙103 8.657∙103 6.093∙103 5.973∙103 4,9516∙103 5.973∙103
1.166∙104 1.204∙104 1,013∙104 8.332∙103 1.159∙104 8.332∙103 5.043∙103 4.435∙103 5.043∙103
8.461∙103 7.126∙103 9,442∙103 6,731∙103 7.925∙103 6.731∙103 4.814∙103 3,921∙103 4.814∙103
6.004∙103 9.939∙103 1.492∙104 5.502∙103 6.564∙103 5.502∙103 6.638∙103 5,973∙103 6.638∙103
6.856∙103 1.157∙104 1,307∙104 5.928∙103 8.992∙103 5.928∙103 5.125∙103 5,254∙103 5.125∙103
1.682∙104 6.859∙103 9,787∙103 1.091∙103 7.025∙103 1,091∙104 4.929∙103 4.334∙103 4.929∙103
6.461∙103 6.781∙103 1.307∙104 5.731∙103 8.19∙103 5.731∙103 6.023∙103 2.056∙103 6.023∙103
2.894∙103 1.151∙104 9,669∙103 3.947∙103 3.94∙103 3.947∙103 4.89∙103 5,864∙103 4.89∙103
1.024∙104 9.444∙103 1.636∙104 7.619∙103 6.607∙103 7.619∙103 7.119∙103 4,196∙103 7.119∙103
1.077∙104 1.55∙104 1,192∙104 7.887∙103 7.833∙103 7.887∙103 5.641∙103 5,134∙103 5.641∙103
7.093∙103 6.932∙103 7,432∙103 6.046∙103 8.977∙103 6.046∙103 4.144∙103 3,793∙103 4.144∙103
9.629∙103 1.277∙104 1.298∙104 7.314∙103 7.443∙103 7.314∙103 5.992∙103 5,681∙103 5.992∙103
5.26∙103 7,381∙103 1,001∙104 5.13∙103 8.482∙103 5.13∙103 5.002∙103 5,109∙103 5.002∙103
1.235∙104 1.271∙104 4,243∙103 8.675∙103 7.342∙103 8.675∙103 3.081∙103 4.841∙103 3.081∙103
1.015∙104 6.961∙103 1,055∙104 7.575∙103 8.486∙103 7.375∙103 5.183∙103 5,514∙103 5.183∙103

 


 


Приложение А

Значения функции Лапласа


 

Приложение Б

Значение в зависимости от и






ТОП 5 статей:
Экономическая сущность инвестиций - Экономическая сущность инвестиций – долгосрочные вложения экономических ресурсов сроком более 1 года для получения прибыли путем...
Тема: Федеральный закон от 26.07.2006 N 135-ФЗ - На основании изучения ФЗ № 135, дайте максимально короткое определение следующих понятий с указанием статей и пунктов закона...
Сущность, функции и виды управления в телекоммуникациях - Цели достигаются с помощью различных принципов, функций и методов социально-экономического менеджмента...
Схема построения базисных индексов - Индекс (лат. INDEX – указатель, показатель) - относительная величина, показывающая, во сколько раз уровень изучаемого явления...
Тема 11. Международное космическое право - Правовой режим космического пространства и небесных тел. Принципы деятельности государств по исследованию...



©2015- 2017 pdnr.ru Все права принадлежат авторам размещенных материалов.