![](/user_photo/_userpic.png)
книги / Управление качеством
..pdf![](/html/65386/197/html_5ycUaVQcyv.6WZX/htmlconvd-jSxjEG261x1.jpg)
Интенсивностьотказаможнотакжеопределитьпоформуле(4.15а):
λ* (3050) = a * (3050) = 0, 0025 = 6, 7 10−31 / ч. P * (3050) 0,375
Средняя наработка до отказа − это математическое ожида-
ние наработки объекта до первого отказа. Средняя наработка до отказа Т1 определяется по формуле:
∞ |
∞ |
|
|
d P (t ) |
|
|
|
|
|||
T1 = ∫t f (t )d t = ∫t |
− |
|
d t. |
||
|
|||||
0 |
0 |
|
|
d t |
Интегрируя это выражение по частям, получим:
|
|
|
∞ |
|
|
|
|
|
|
|
∫ |
|
|
|
|
|
|||
T1 = − t P (t ) |
0∞ |
+ |
0 |
P (t ) d t. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Практически во всех случаях член t P (t ) |
|
0∞ |
обращается в |
||||||
|
|||||||||
|
|
|
|
|
|
|
|
|
|
нуль, так как имеет место соотношение lim tP (t ) = 0, |
следующее из |
||||||||
|
|
|
|
|
t →∞ |
|
|
|
|
того, что срок службы любого объекта ограничен, и потому Р(t) быстрее стремится к нулю, чем t → ∞ . Таким образом, получаем
T1 = ∞∫ P (t )d t, |
(4.18) |
0 |
|
т.е. средняя наработка до отказа численно равна площади под кривой P(t).
Статистическая средняя наработка до отказа группы однотипных объектов определяется как отношение суммарного значения наработки каждого из них до появления отказа к общему числу объектов N, работоспособных в момент t = 0,
|
1 |
n |
|
|
T1 = |
∑ti , |
(4.19) |
||
N |
||||
|
i =1 |
|
где ti − наработка до отказа i-го объекта.
261
![](/html/65386/197/html_5ycUaVQcyv.6WZX/htmlconvd-jSxjEG262x1.jpg)
Имея данные о количестве вышедших из строя элементов ni в каждом i-м интервале времени, среднюю наработку до отказа лучше определять из уравнения
T1 |
≈ |
m |
|
|
|
∑ni tсрi |
/ N. |
(4.19а) |
|||
|
i =1 |
|
|
|
В выражении (4.19а) tсрi m находятся по следующим формулам:
tсрi = (ti–1 + ti) /2, m = tk/∆t,
где ti–1 – время начала i-го интервала; ti – время конца i-го интервала; tk – время, в течение которого вышли из строя все элемен-
ты; ∆t = ti–1 – ti – интервал времени.
Поскольку практически невозможно осуществить испытания всех элементов до отказа, то в первом приближении при большом числе N среднюю наработку доотказа можноопределитьзависимостью
T1 ≈ |
t1 + t2 + ...tN + (N − m)t |
, |
(4.19б) |
|
|||
|
N |
|
где N – число элементов, поставленных на испытания; m – число отказавших элементов; t − время испытания.
Формула (4.19б) справедлива при числе m отказавших элементов, близких к N.
Пример 2. На испытания поставлено N = 10 невосстанавливаемых элементов. Испытания проводились в течение времени t = 100 ч. В процессе проведения испытаний отказало 8 элементов, при этом отказы зафиксированы в следующие моменты времени:
t1 = 20 ч, t2 = 30 ч, t3 = 50 ч, t4 = 30 ч, t5 = 40 ч, t6 = 60 ч, t7 = 70 ч, t8 = 60 ч. Оставшиеся два элемента не отказали. Определить сред-
нюю наработку до отказа.
Решение. Вычислим наработку до отказа для невосстанавливаемого элемента по формуле (4.19б):
T1 = 20 + 30 + 50 + 30 + 40 + 60 + 70 + 60 + (10 − 8)100 = 56 ч. 10
262
Пример 3. Найти интенсивность отказов λ(t) и построить график изменения кривой интенсивности отказов по данным, представленным в табл. 4.1. На испытания поставлено N = 100 элементов, испытания проводились в течение времени t = 100 ч.
Таблица 4 . 1
Исходные данные результатов испытаний
Интервал |
Числоотказов |
ЧислоN (ti) работо- |
|
|
|
времени, |
винтервале, |
способныхэлементов |
∆ti |
n (∆ti) |
Nр (ti) |
∆ t= ti− ti −1 |
n (∆ti) |
кмоменту времениti |
|
|
|
0 … 10 |
10 |
90 |
50 … 60 |
3 |
68 |
10 … 20 |
9 |
81 |
60 … 70 |
2 |
66 |
20 … 30 |
6 |
75 |
70 … 80 |
5 |
61 |
30 … 40 |
2 |
73 |
80 … 90 |
9 |
52 |
40 … 50 |
2 |
71 |
90… 100 |
10 |
42 |
Для построения кривой интенсивности отказов (рис. 4.5) воспользуемся формулой (4.15):
λ(t1 ) = |
n (∆ |
|
t1 ) |
= |
10 |
|
= 1,11 |
10−2 , λ |
(t2 ) |
= |
|
9 |
= 1,11 10−2 , |
||||||
Nр (t1 )∆ |
t1 |
|
|
|
81 10 |
||||||||||||||
|
|
|
|
|
90 10 |
|
|
|
|
|
|
||||||||
λ(t3 ) = |
6 |
|
|
= 0,80 |
10−2 , |
λ(t4 ) = |
2 |
|
|
= 0, 27 |
10−2 , |
||||||||
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
10 |
|
73 10 |
|||||||||||||||
|
75 |
|
|
|
|
|
|
|
|
|
|
||||||||
λ(t5 ) = |
2 |
|
|
= 0, 28 |
10−2 , |
λ(t6 ) = |
3 |
|
|
= 0, 44 |
10−2 , |
||||||||
|
|
|
|
|
|
|
|
|
|||||||||||
71 10 |
|
68 10 |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
λ(t7 ) = |
|
|
2 |
|
|
|
= 0,33 10−2 , |
λ(t8 ) = |
5 |
|
= 0,82 |
10−2 , |
|||||||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
66 10 |
|
|
61 10 |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
λ(t9 ) = |
9 |
|
|
= 1,75 |
10−2 , |
λ(t10 ) = |
10 |
|
= 2,38 |
10−2 . |
|||||||||
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
10 |
|
42 |
10 |
||||||||||||||
|
52 |
|
|
|
|
|
|
|
|
|
|
263
![](/html/65386/197/html_5ycUaVQcyv.6WZX/htmlconvd-jSxjEG264x1.jpg)
Рис. 4.5. Кривая интенсивности отказов
Пример 4. Система состоит из N = 12600 элементов, средняя интенсивность отказов которых λср = 0,32 · 10–6 1/ч. Необходимо определить вероятность безотказной работы в течение 50 ч и вычислить среднюю наработку до первого отказа.
Решение. Интенсивность отказов системы определяют по формуле λс = λcр · N = 0,32·10–6 · 12600 = 4,032 · 10–3 1/ч. Тогда на ос-
новании (4.16) получаем
P(50) = e−4,032 10−3 50 ≈ 0,82
исреднюю наработку до отказа T1 = 1 (4, 032 10−3 ) ≈ 250 ч.
Пример 5. На испытании находилось N = 1000 образцов неремонтируемой аппаратуры. Число отказов n(∆t) фиксировалось каждые 100 ч работы (∆t = 100 ч). Данные об отказах приведены в табл. 4.2. Требуется вычислить количественные характеристики надежности.
Решение. Аппаратура относится к классу невосстанавливаемых изделий. Поэтому критериями надежности будут P*(t), a*(t),
λ*(t), T1*.
Вычислим P*(t). На основании формулы (4.6) имеем:
264
P * (100) |
= 1 − |
50 |
|
|
= |
0,95, |
||
1000 |
|
|
||||||
|
|
|
|
|
|
|||
P * (200) |
= 1 − |
|
90 |
|
|
= |
0,91, |
|
1000 |
|
|||||||
|
|
|
|
|
||||
.................................................. |
||||||||
P * (3000) |
= 1 − |
575 |
= |
0, 425. |
||||
|
||||||||
|
|
1000 |
|
|
Таблица 4 . 2
Данные об отказах
∆ti, ч |
|
n (∆ti) |
|
|
|
|
∆ti, ч |
|
|
|
|
n (∆ti) |
|
∆ti, ч |
n (∆ti) |
|||||
0…100 |
|
50 |
|
|
|
1000…1100 |
|
|
|
|
15 |
|
2000…2100 |
12 |
||||||
100…200 |
|
40 |
|
|
|
1100…1200 |
|
|
|
|
14 |
|
2100…2200 |
13 |
||||||
200…300 |
|
32 |
|
|
|
1200…1300 |
|
|
|
|
14 |
|
2200…2300 |
12 |
||||||
300…400 |
|
25 |
|
|
|
1300…1400 |
|
|
|
|
13 |
|
2300…2400 |
13 |
||||||
400…500 |
|
20 |
|
|
|
1400…1500 |
|
|
|
|
14 |
|
2400…2500 |
14 |
||||||
500…600 |
|
17 |
|
|
|
1500…1600 |
|
|
|
|
13 |
|
2500…2600 |
16 |
||||||
600…700 |
|
16 |
|
|
|
1600…1700 |
|
|
|
|
13 |
|
2600…2700 |
20 |
||||||
700…800 |
|
16 |
|
|
|
1700…1800 |
|
|
|
|
13 |
|
2700…2800 |
25 |
||||||
800…900 |
|
15 |
|
|
|
1800…1900 |
|
|
|
|
14 |
|
2800…2900 |
30 |
||||||
900…1000 |
|
14 |
|
|
|
1900…2000 |
|
|
|
|
12 |
|
2900…3000 |
40 |
||||||
Для расчета характеристик a(t) и λ(t) применим формулы |
||||||||||||||||||||
(4.12) и (4.15), тогда: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
a * (50) |
= |
|
50 |
|
|
= |
0,5 10−3 1 / ч, |
|
|||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
1000 100 |
|
|
|
|
|
||||||||
|
a * (150) |
= |
40 |
|
|
= |
0, 4 10−3 1 / ч, |
|
||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||
|
1000 |
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
100 |
|
|
|
|
|
||||||||
................................................................. |
|
|||||||||||||||||||
|
a * (2950) = |
|
40 |
|
|
|
= |
0, 4 |
10−3 1 / ч; |
|
||||||||||
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
100 |
|
||||||||||||||||
|
|
|
|
|
|
|
1000 |
|
|
|
|
|
||||||||
λ* (50) |
= |
|
|
|
|
50 |
|
|
|
|
|
|
= |
0,514 10−3 1 / ч, |
|
|||||
|
|
|
||||||||||||||||||
[(1000 + 950) / 2] 100 |
|
|||||||||||||||||||
λ* (150) = |
|
|
|
|
40 |
|
|
|
|
|
= |
0, 43 10−3 1 / ч, |
|
|||||||
|
|
|
||||||||||||||||||
|
[(950 + 910) / 2] 100 |
|
265
![](/html/65386/197/html_5ycUaVQcyv.6WZX/htmlconvd-jSxjEG266x1.jpg)
.....................................................................................
λ* (2950) = |
40 |
= 0,9 10−3 1 / ч. |
[(465 + 425) / 2] 100 |
Значения P*(t), a*(t) и λ*(t), вычисленные для всех ∆ti, приведены в табл. 4.3.
|
|
|
|
Таблица 4 . 3 |
|
Вычисленные значения P*(t), a*(t), λ*(t) |
|||
|
|
|
|
|
∆ti, ч |
|
P*(t) |
a*(t), 10–3 ·1/ч |
λ*(t), 10–3 ·1/ч |
0…100 |
|
0,950 |
0,50 |
0,514 |
100…200 |
|
0,910 |
0,40 |
0,430 |
200…300 |
|
0,878 |
0,32 |
0,358 |
300…400 |
|
0,853 |
0,25 |
0,289 |
400…500 |
|
0,833 |
0,20 |
0,238 |
500…600 |
|
0,816 |
0,17 |
0,206 |
600…700 |
|
0,800 |
0,16 |
0,198 |
700…800 |
|
0,784 |
0,16 |
0,202 |
800…900 |
|
0,769 |
0,15 |
0,193 |
900…1000 |
|
0,755 |
0,14 |
0,184 |
1000…1100 |
|
0,740 |
0,15 |
0,200 |
1100…1200 |
|
0,726 |
0,14 |
0,191 |
1200…1300 |
|
0,712 |
0,14 |
0,195 |
1300…1400 |
|
0,699 |
0,13 |
0,184 |
1400…1500 |
|
0,685 |
0,14 |
0,202 |
1500…1600 |
|
0,672 |
0,13 |
0,192 |
1600…1700 |
|
0,659 |
0,13 |
0,195 |
1700…1800 |
|
0,646 |
0,13 |
0,200 |
1800…1900 |
|
0,632 |
0,14 |
0,220 |
1900…2000 |
|
0,620 |
0,12 |
0,192 |
2000…2100 |
|
0,608 |
0,12 |
0,195 |
2100…2200 |
|
0,595 |
0,13 |
0,217 |
2200…2300 |
|
0,583 |
0,12 |
0,204 |
2300…2400 |
|
0,570 |
0,13 |
0,225 |
2400…2500 |
|
0,556 |
0,14 |
0,248 |
2500…2600 |
|
0,540 |
0,16 |
0,290 |
2600…2700 |
|
0,520 |
0,20 |
0,376 |
2700…2800 |
|
0,495 |
0,25 |
0,490 |
2800…2900 |
|
0,465 |
0,30 |
0,624 |
2900…3000 |
|
0,425 |
0,40 |
0,900 |
266
![](/html/65386/197/html_5ycUaVQcyv.6WZX/htmlconvd-jSxjEG267x1.jpg)
Cледует иметь в виду, что в табл. 4.3 данные P*(t) приведены для концов интервалов ∆ti, а данные для a* (t) и λ*(t) – для середины интервалов ∆ti. Поэтому определение λ(t) по формуле (4.15') и данным табл. 4.2 не даст значений λ*(t), указанных в табл. 4.3.
Вычислим среднее время безотказной работы, предположим, что на испытании находились только те образцы, которые отказали.
По формуле (4.19а), учитывая, что в данном случае m = tk/∆t = = 3000/100 = 30 и N = 575, имеем
Т1 = 50 50 + 40 150 + 32 250 + ... + 30 2850 + 40 2950 =1400 ч. 575
Полученное значение средней наработки до отказа является заниженным, так как опыт был прекращен после отказа 575 образцов из 1000, поставленных на испытание.
Средняя наработка на отказ − отношение суммарной нара-
ботки восстанавливаемого объекта к математическому ожиданию числа его отказов в течение этой наработки. В соответствии с определением средняя наработка на отказ Т0 задается формулой
T0 = |
t |
, |
(4.20) |
M { r (t )} |
где t − суммарная наработка; r (t) − число отказов, наступивших в течение времени наработки t; M { r (t )} − математическое ожида-
ние этого числа отказов.
Статистическая оценка средней наработки на отказ вычисля-
ется по формуле |
|
T0 = t / r (t ) . |
(4.21) |
Если на испытании находится N образцов в течение времени t, то наработка на отказ вычисляется по формуле
T0 |
|
N n j |
|
N |
|
= |
∑∑tij |
/ ∑ n j , |
(4.21a) |
||
|
j=1 i=1 |
|
j=1 |
|
267
![](/html/65386/197/html_5ycUaVQcyv.6WZX/htmlconvd-jSxjEG268x1.jpg)
где tij − время исправной работы j-го образца изделия между (i − 1) и i-м отказами; nj – число отказов за время работы tj-го образца.
Если в процессе испытаний или эксплуатации в качестве статистической оценки получена средняя наработка на отказ, то вероятность безотказной работы элемента за время t определяется по формуле
Pi (t ) = e |
− |
tр +tхр 10−3 |
|
|
|
T0 |
(4.22) |
||||
|
|||||
|
, |
где tp − время работы, tхр − время хранения. Например, устройство работает 10 ч, а общее время работы изделия 50 ч. Следовательно, время хранения составит 40 ч.
Пример 6. При испытаниях некоторого устройства, имеющего время работы до отказа, распределенное по нормальному закону, получено 10 реализаций наработки до отказа (в часах): 120, 110, 80, 130, 120, 140, 80, 150, 130, 140. Требуется найти доверительные границы с доверительной вероятностью 0,95 для средней наработки на отказ и дисперсии Дх.
|
|
|
|
1 |
10 |
|
1 |
|
|
|
|
Решение. Имеем |
|
|
= |
∑ti =120 |
ч, Дx = |
|
∑(ti − |
|
)2 = |
||
|
t |
t |
|||||||||
|
|
N −1 |
|||||||||
|
|
|
10 i=1 |
|
|
|
|
||||
= 576 ч2, S = Дx = 24 |
ч. |
|
|
|
|
|
|
|
По формулам (3.22) и (3.23) определим границы средней наработки на отказ:
|
|
− tα( N −1) |
S |
≤ tcp≤ |
|
|
S |
. |
|
|
t |
t+ tα( N −1) |
|||||||
|
|
|
|||||||
|
|
|
N |
|
|
|
N |
||
Согласно табл. П1 |
имеем |
tα( N −1) = t (0,95; 9) =1,83 . Отсюда |
для tср получаем значения нижней и верхней границы средней наработки на отказ соответственно 105,4 и 134,6 ч. Вычисление доверительного интервала для дисперсии ведется по формуле (3.24).
Имеем N = 10, S2 = 576 ч, χкр2 (0,95; 9) = 16,92, χкр2 (0,05; 9) = 3,33 (см. табл. П. 2). Тогда доверительный интервал для дисперсии составит 430 и 1735 ч.
268
Параметр потока отказов − отношение математического ожидания числа отказов восстанавливаемого объекта за достаточно малую его наработку к значению этой наработки. Параметр потока отказов ω(t ) определяется по формуле
ω(t ) = lim |
M { r (t + ∆ t )− r (t )} |
|
|
|
, |
(4.23) |
|
|
|||
∆ t→ 0 |
∆ t |
|
где r (t + ∆ t ), r (t ) − количество отказов к моментам времени t + ∆ t
и t соответственно. |
|
||
Статистическая оценка параметра потока отказов ω (t ) |
опре- |
||
деляется по формуле |
|
||
ω (t ) = |
r (t2 ) − r (t1 ) |
, |
(4.24) |
|
|||
|
t2 − t1 |
|
где r(t2), r(t1) − количество отказов к моментам времени t2 и t1 соответственно.
Для простейшего потока отказов, т.е. для потока, обладающего свойством стационарности, ординарности и отсутствия последствия, параметр потока отказов совпадает с интенсивностью отказов, т.е.
ω(t ) = λ(t ) = λ = const.
Рассмотрим показатели ремонтопригодности.
Вероятность восстановления − это вероятность того, что вре-
мя восстановления τ работоспособного состояния объекта не превысит заданного значения t. По определению
PB (t ) = P{τ ≤ t} , |
(4.25) |
где РВ(t) − вероятность восстановления объекта.
Время восстановления − это то время, которое затрачивается на обнаружение, поиск места отказа и устранение последствий отказа. Вероятность восстановления представляет собой функцию распределения времени восстановления объекта, т.е. интегральный
269
![](/html/65386/197/html_5ycUaVQcyv.6WZX/htmlconvd-jSxjEG270x1.jpg)
закон распределения времени восстановления, а производная по функции РВ(t) по времени дает плотность распределения случайной величины fB(t) (дифференциальный закон распределения) времени восстановления
fB (t ) = d PB (t ) . d t
Интенсивность восстановления − условная плотность веро-
ятности восстановления работоспособного состояния объекта, определенная для рассматриваемого момента времени при условии, что до этого момента восстановление не было завершено. Интенсивность восстановления µ (t ) определяется по формуле:
µ (t ) = |
|
fB (t ) |
= |
|
1 |
|
d PB (t ) |
. |
(4.26) |
|
− PB (t ) |
1− PB (t ) |
|
||||||
1 |
|
|
d t |
|
Интегрируя соотношение (2.26) в пределах (0, t), получаем
|
t |
|
|
|
∫ |
|
(4.27) |
PB (t ) =1 − exp − |
|
µ (t )dt . |
|
|
0 |
|
|
Полученное выражение устанавливает связь между интенсивностью и вероятностью восстановления для любых законов распределения времени восстановления.
Статистическая оценка интенсивности восстановления µ (t ) имеет вид:
µ (t ) = |
n (t + ∆ t )− |
n (t ) |
, |
(4.28) |
NB (t ) ∆ |
|
|||
|
t |
|
где n(t), n (t + ∆ t ) − количество объектов, для которых восстановление закончилось к моментам времени t, (t + ∆ t ) ; NB(t) − количество объектов, ожидающих восстановления к началу промежутка времени (t, t + ∆ t ) .
270