Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Управление качеством

..pdf
Скачиваний:
1
Добавлен:
12.11.2023
Размер:
21.77 Mб
Скачать

Интенсивностьотказаможнотакжеопределитьпоформуле(4.15а):

λ* (3050) = a * (3050) = 0, 0025 = 6, 7 1031 / ч. 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

Имея данные о количестве вышедших из строя элементов 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= titi 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

102 , λ

(t2 )

=

 

9

= 1,11 102 ,

Nр (t1 )

t1

 

 

 

81 10

 

 

 

 

 

90 10

 

 

 

 

 

 

λ(t3 ) =

6

 

 

= 0,80

102 ,

λ(t4 ) =

2

 

 

= 0, 27

102 ,

 

 

 

 

 

 

 

 

 

 

 

10

 

73 10

 

75

 

 

 

 

 

 

 

 

 

 

λ(t5 ) =

2

 

 

= 0, 28

102 ,

λ(t6 ) =

3

 

 

= 0, 44

102 ,

 

 

 

 

 

 

 

 

 

71 10

 

68 10

 

 

 

 

 

 

 

 

 

 

 

 

λ(t7 ) =

 

 

2

 

 

 

= 0,33 102 ,

λ(t8 ) =

5

 

= 0,82

102 ,

 

 

 

 

 

 

 

 

 

 

 

66 10

 

 

61 10

 

 

 

 

 

 

 

 

 

 

 

 

 

λ(t9 ) =

9

 

 

= 1,75

102 ,

λ(t10 ) =

10

 

= 2,38

102 .

 

 

 

 

 

 

 

 

 

 

 

10

 

42

10

 

52

 

 

 

 

 

 

 

 

 

 

263

Рис. 4.5. Кривая интенсивности отказов

Пример 4. Система состоит из N = 12600 элементов, средняя интенсивность отказов которых λср = 0,32 · 10–6 1/ч. Необходимо определить вероятность безотказной работы в течение 50 ч и вычислить среднюю наработку до первого отказа.

Решение. Интенсивность отказов системы определяют по формуле λс = λ· N = 0,32·10–6 · 12600 = 4,032 · 10–3 1/ч. Тогда на ос-

новании (4.16) получаем

P(50) = e4,032 103 50 ≈ 0,82

исреднюю наработку до отказа T1 = 1 (4, 032 103 ) ≈ 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 103 1 / ч,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1000 100

 

 

 

 

 

 

a * (150)

=

40

 

 

=

0, 4 103 1 / ч,

 

 

 

 

 

 

 

 

 

 

1000

 

 

 

 

 

 

 

 

 

 

100

 

 

 

 

 

.................................................................

 

 

a * (2950) =

 

40

 

 

 

=

0, 4

103 1 / ч;

 

 

 

 

 

 

 

 

 

 

 

 

100

 

 

 

 

 

 

 

 

1000

 

 

 

 

 

λ* (50)

=

 

 

 

 

50

 

 

 

 

 

 

=

0,514 103 1 / ч,

 

 

 

 

[(1000 + 950) / 2] 100

 

λ* (150) =

 

 

 

 

40

 

 

 

 

 

=

0, 43 103 1 / ч,

 

 

 

 

 

[(950 + 910) / 2] 100

 

265

.....................................................................................

λ* (2950) =

40

= 0,9 103 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

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

где 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

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

fB (t ) = d PB (t ) . d t

Интенсивность восстановления условная плотность веро-

ятности восстановления работоспособного состояния объекта, определенная для рассматриваемого момента времени при условии, что до этого момента восстановление не было завершено. Интенсивность восстановления µ (t ) определяется по формуле:

µ (t ) =

 

fB (t )

=

 

1

 

d PB (t )

.

(4.26)

 

PB (t )

1PB (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