К 1. 6. Надежность в период постепенных отказов

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

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

Нормальное распределение хорошо описывает распределение вероятностей наработки до отказа, ресурса элементов и других показателей надежности, когда они зависят от большого числа однородных по своему влиянию случайных факторов, влияние каждого из которых по сравнению с совокупностью всех остальных незначительно. Этот закон характерен для постепенных отказов, вызванных износом и старением. Нормальному распределению «подчиняется» наработка до отказа многих восстанавливаемых и невосстанавливаемых объектов, размеры и ошибки измерений деталей и др. Нормальное распределение является наиболее универсальным, удобным и широко применяемым для практических расчетов (рис. 1.5, 1.6).

Плотность распределения

Распределение имеет два независимых параметра: математическое ожидание m_t и среднее квадратическое отклонение S _t. Значения параметров m_t и S_t оценивают по результатам испытаний по формулам (4.5) и (4.7).

Математическое ожидание определяет на графике (см. рис. 1.5) положение петли, а среднее квадратическое отклонение — ширину петли.

Кривая плотности распределения тем острее и выше, чем меньше S_t. Она начинается от t = – ∞ и распространяется до t = + ∞. Это не является существенным недостатком, особенно если m_t≤ 3 S_t, так как площадь, очерченная уходящими в бесконечность ветвями кривой плотности, выражающая соответствующую вероятность отказов, очень мала. Так, вероятность отказа за период времени до m_t –3 S_t составляет всего 0,135% и обычно не учитывается в расче­тах. Вероятность отказа до m_t –2 S_t равна 2,175%. Наибольшая ордината кривой плотности распределения равна 0,399/ S_t.

Интегральная функция распределения .

Вероятность отказа и вероятность безотказной работы соответственно Q(t)=F(t); P(f)= 1– F(t).

Вычисление интегралов заменяют использованием таблиц. Таблицы для нормального распределения в функции (t – m_t) и S_t были бы громоздкими, так как имели бы два независимых параметра. Можно обойтись небольшими таблицами для нормального распределения, у которого m_x =0 и S_x= 1. Для этого распределения функция плотности имеет одну переменную х. Величина х является центрированной, так как т_х = 0, и нормированной, так как S_x= 1. Функция плотности распределения записывается в относительных координатах с началом на оси симметрии петли.

Функция распределения – интеграл от плотности распределения

.

Из этого уравнения следует, F₀ (x) + F₀ (– x) = 1, отсюда F₀ (– x) = 1– F₀ (x)

Для использования таблиц следует применять подстановку x = (t — m_t)/S_t;при этом х называется квантилью нормированного нормального распределения и обычно обозначается u_р.

Плотность распределения и вероятность безотказной работы соответственно f(t)=f_o(x)/S_t Q(t)=F₀(x); P(t) = 1– F₀(x), где f_o(x) и F₀(x) берут по таблицам [5, 50]. Например:

x …………………………
fo(x) ………………………	0,3989	0,2420	0,0540	0,0044	0,0001
F0(x) ……………………..	0,5	0,8413	0,9772	0,9986	0,9999

В табл. 1.1 приведены непосредственно значения P(t) в зависимости от x=u_p=(t – m_t)/S_t в употребительном диапазоне. В литературе по надежности часто вместо интегральной функции распределения F₀(x) пользуются функцией Лапласа:

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

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

; .

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

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

Значения этой наработки (времени) определяют с помощью квантилей нормированного нормального распределения t=m_t + u_p S_t. Значения квантилей даются в таблицах в зависимости от требуемой вероятности, в частности от вероятности безотказной работы.

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

Распределение суммы независимых случайных величин U = X+Y+Z, называемое композицией распределений, при нормальном распределении слагаемых также является нормальным распределением.

Математическое ожидание и дисперсия композиции соответственно равны

; ,

где т_х, т_у, m_z – математические ожидания случайных величин X, Y, Z; S_x, S_y, S_z – дисперсия тех же величин.

Пример. Оценить вероятность Р (t) безотказной работы в течение t =1,5-10⁴ ч изнашиваемого подвижного сопряжения, если ресурс по износу подчиняется нормальному распределению с параметрами m_t = 4-10⁴ч, S_t =10⁴ ч.

Решение. Находим квантиль ; по табл. 1.1

определяем, что Р (t) =0,9938.

Таблица 1.1

Нормальное распределение
Квантиль up	Вероятность безотказной работы P (t)	Квантиль up	Вероятность безотказной работы P (t)	Квантиль up	Вероятность безотказной работы P (t)
0,000 –0,1 –0,126 –0,2 –0,253 –0,3 –0,385 –0,4 –0,5 –0,524 –0,6 –0,674 –0,7 –0,8 –0,842 –0,9 –1,0 –1,036	0,5 0,5398 0,55 0,5793 0,6 0,6179 0,65 0,6554 0,6915 0,7 0,7257 0,75 0,7580 0,7881 0,8 0,8159 0,8413 0,85	–1,1 –1,2 –1,282 –1,3 –1,4 –1,5 –1,6 –1,645 –1,7 –1,751 –1,8 –1,881 –2,0 –2,054 –2,1 –2,170 –2,2 –2,3	0,8643 0,8849 0,9 0,9032 0,9192 0,9332 0,9452 0,95 0,9554 0,96 0,9641 0,97 0,9772 0,98 0,9821 0,985 0,9861 0,9893	–2,326 –2,4 –2,409 –2,5 –2,576 –2,6 –2,652 –2,7 –2,748 –2,8 –2,878 –2,9 –3,0 –3,09 –3,291 –3,5 –3,719 –4,0	0,99 0,9918 0,992 0,9938 0,995 0,9953 0,996 0,9965 0,997 0,9974 0,998 0,9981 0,9986 0,999 0,9995 0,9998 0,9999 0,99997

Поскольку длительность безотказной работы объекта не может быть отрицательной, нормальное распределение в общем виде к задачам надежности может применяться только при S << m_t. Если условие не выполняется, то применяется усеченный нормальный закон. Сущ­ность усечения заключается в том, что из совокупности значений случайной величины исключаются все значения t <0 и f(t)= 0при t< 0.

Функция плотности усеченного нормального закона распределения

,

где А — нормирующий постоянный множитель. Значения А вычисляют по формуле

,

где Ф(m_t / S_t) — функция Лапласа.

Усеченное нормальное распределение с хорошей точностью описывает результаты анализа надежности сложных систем с учетом ухода параметров элементов в процессе эксплуатации (точности, прочности, вибраций, температуры и др.) за допустимые пределы. Тем не менее при m_t / S_t >2, что имеет место в большинстве случаев при расчетах надежности рассматриваемых машин с нормальным распределением времени безотказной работы, множитель А мало отличается от единицы, и усеченное нормальное распределение достаточно точно описывается обычным нормальным законом.

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

Плотность вероятности наработки до отказа для распределения Вейбулла

,

где a — параметр масштаба (задает масштаб кривой распределения по оси абсцисс);

b — параметр формы (определяет остроту и асимметрию кривой плотности распределения).

Вероятность безотказной работы

.

Вероятность возникновения отказа

Интенсивность отказов .

Кривые изменения плотности вероятности и вероятности безотказ­ной работы представлены на рис. 4.4, а, б.

Величины a и b всегда положительны. При b = 1распределение Вейбулла переходит в экспоненциальное, которое, таким образом, является частным случаем распределения Вейбулла. При b < 1 интенсивность отказов f(t) становится убывающей функцией времени, поэтому закон Вейбулла с параметром b < 1 можно использовать для оценки надежности объектов в период их приработки. При b > 1 распределение Вейбулла характеризуется возрастающей функцией интенсивности отказов, и его удобно использовать для оценки надежности износовых отказов. При b = 3,5...4,0 распределение Вейбулла близко к нормальному, и функция плотности вероятности приобретает колоколообразную форму.

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

Для расчета параметров a и b по опытным данным удобно пользо­ваться методом наименьших квадратов.

Дата добавления: 2015-07-15