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

книги из ГПНТБ / Золотарь, И. А. Экономико-математические методы в дорожном строительстве

.pdf
Скачиваний:
10
Добавлен:
21.10.2023
Размер:
9.35 Mб
Скачать

Наработка на отказ

2 *

 

4 = — ----- ,

(XII.29)

п

 

где п — число отказов; U — время исправной работы

между

(г— 1)-м и г'-м отказами.

 

В вероятностной форме tcр выражается через со (О-

 

 

( х м » )

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

Для элементов разового (т. е. до первого отказа) использования целесообразно в качестве количественных характеристик надежно­ сти использовать P{t)\ Т\ К(t) и a(t).

Т а б л и ц а 34

Класс элемен­

Наличие ре­

Состояние

тов устройств

зервирования

эле ментов

 

 

устройств

В работе

 

Количественные характеристики

Характер

 

 

надежности

 

 

 

 

 

 

отказов

 

 

 

а СС(0

 

 

P ( t ) Т

*ср

),(/)

 

 

 

 

 

Мгновенные

+

_1_

+ ! -

+

 

 

 

Отсутствует

Постепенные

+

+

+

+

Разового

 

 

В хранении Постепенные

_L

+

 

 

 

 

 

 

 

 

 

пользования

 

Мгновенные

+

+

+

 

 

 

 

 

Имеется

В работе

 

 

 

 

 

Постепенные| +

+

+

 

 

 

В хранении Постепенные

 

 

 

 

 

Мгновенные|

—•

+ *

+

 

+ *

 

 

Отсутствует

В работе

 

 

 

 

 

 

 

 

Постепенные

+

+

Многократ­

 

В хранении Постепенные

+

 

ного пользо­

 

 

 

 

 

 

 

 

вания

 

 

Мгновенные

+

 

 

Имеется

В работе

 

 

 

 

 

 

 

 

Постепенные

 

 

 

+

 

 

 

В хранении Постепенные

| -

|+ * * |

-

| +

 

* Характ еристики Т, а( t)

и Я(£) в да?ihom случае могут

применяться

вместо t cР И

(о(0 лишь в случае распреде ления времени отказов по эксп оненциальному закону,

 

** i

r

данном случае

может приме!чяться, если в

процессе

хранения

устройс тва

гср

в

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

210

Для элементов (устройств) многократного использования тако­ выми являются ^Ср и со (t), причем количественная оценка надежно­ сти в таких случаях более сложна. Элементы (устройства), находя­ щиеся на хранении, в случае их выхода из строя заменяются. Поэ­ тому здесь снова действительны количественные характеристики Др и со(t). Таким образом, выбор количественных характеристик надежности зависит от класса элементов (устройств), наличия или отсутствия резервирования, состояния элементов (в хранении или в работе) и характера отказов. Табл. 34, заимствованная из работы А. М. Половко, дает рекомендации по выбору характеристик для ко­ личественной оценки надежности с учетом упомянутых факторов. Знак « + » в табл. 34 означает пригодность, а знак «—» непригод­ ность данного количественного показателя.

 

 

-I т

" 1 -

 

 

4 а

 

 

 

 

 

<Т)

о

 

 

 

М., ч

д0 ' 4

 

н

 

 

 

•■Ск

*

#

*

 

 

s;

 

 

<1

«=

о.

а

 

 

о

50

 

0,950

0 ,5 0 х

0,514х

1500—1600

о о i

 

 

 

 

 

х ю - з

хЮ -з 1600—1700

100—200

40

 

0,910

0,40

0,430

1700— 1800

200—300

32

 

0,878

0,32

0,358

1800— 1900

300—400

25

 

0,853

0,25

0,289

1900—2000

400—500

20

 

0,833

0,20

0,238

2000—2100

500 -600

17

 

0,816

0,17

0,206

2100—2200

600—700

16

0,800

0,16

0,198

2200—2300

700—800

16

0,784

0,16

0,202

2300—2400

800—800

15

0,769

0,15

0,193

2400—2500

900— 1000

14

0,755

0,14

0,184

2500—2600

000— 1100

15

0,740

0,15

0,200

2600—2700

100— 1200

14

0,726

0,14 0,191

2700-2800

1200— 1300

14

0,712

0,14

0,195

2800—2900

1300— 1400

13

0,699 0,13 0,184 2900-3000

1400— 1500

14

0,685

0,14

0,202

 

 

 

Т а б л и ц а З б

Ап., эле­ ментов

 

"1 =■

- 1 -

а.

«

 

 

*•

 

 

13

0,672

0,13

0,192

13

0,659

0,13

0,195

13

0,646

0,13

0,200

14

0,632

0,14

0,220

12

0,620

0,12

0,192

12

0,608

0,12

0,195

13

0,595

0,13

0,217

12

0,583

0,12

0,204

13

0,570

0,13

0,225

14

0,556

0,14

0,248

16

0,540

0,16

0,290

20

0,520

0,20

0,376

25

0,495

0,25

0,490

30

0,465

0,30

0,624

40

0,425

0,40

0,900

Порядок вычислений количественных характеристик надежности на основе статистических данных испытаний элементов (устройств) можно проследить на следующем примере'. Испытанию подверга­ лись Ао=1000 образцов. Общая длительность испытаний составила 3000 ч, причем выход из строя (отказы) элементов Ащ фиксировал­ ся через каждые Д /,= 100 ч. Данные испытаний приведены в табл. 35. Требуется вычислить характеристики P*(t)', a* (t) и X*(t).

Пользуясь формулой (XII.6), найдем:

Р* (100)=

. 1000 — ?0_= о;95;

' '

1ПЛП

1 Пример заимствован из книги А. М. Половко «Основы теории надежности»

(М , «Наука», 1964, 446 с.).

211

Я*(200) =

10QQ- (50- + .12) = 0 91;

 

 

 

1000

Я*(3000)

=

1000~ 575 = 0,425.

v

'

 

1плл

Величины a* (t) и %* (t) в соответствии с их смыслом будем вычислять, относя величины A/ij к середине соответствующего про­ межутка времени ДU.

Тогда

 

50

\Jy\S

 

 

 

 

 

1000-100

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а*(150)

40

: 0,4

 

 

 

 

 

 

 

1000-100

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а*(2950)

40

0,4

 

 

 

 

 

 

1000-100

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X* (50)

50

.

г

 

 

 

 

 

 

 

1000 + 950

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

--------------- 100

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

X* (150)

 

40

. — с

 

 

 

 

 

 

 

950 +

910

_

0,43 -10~3;

 

 

 

 

 

 

 

--------------100

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

X* (2950)

40

. — Г

 

 

 

 

 

 

465 +

425

=

0 ,9 -10-3.

 

 

 

 

 

 

 

----- z -----• 100

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

На рис. 48 нанесены по полученным

значениям P*(t), а : (/) и

K*(t)

плавные кривые для этих характеристик.

Необходимо отме­

тить,

что характер кривой л* (t)

на рис.

48 является типичным для

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

В начальный период до момента tо значение

 

 

 

 

 

 

X* (/) велико вследствие так

 

 

 

 

 

 

называемых

приработочных,

 

 

 

 

 

 

т.

е. ранних,

отказов,

выз­

 

 

 

 

 

 

ванных зачастую дефектами

 

 

 

 

 

 

производственного

изготов­

 

 

 

 

 

 

ления. Затем до момента t\

 

 

 

 

 

 

X*(/)^const. В дальнейшем

 

 

 

 

 

 

вследствие износа элементов

 

 

 

 

 

 

К* (t) нарастает.

 

 

 

 

 

 

 

 

 

Из данных табл. 35

сле­

 

 

 

 

 

 

дует, что вероятность безот­

 

 

 

 

 

 

казной работы устройств в

 

 

 

 

 

 

течение 3000 ч равна 0,425.

Рис.

48.

Зависимость вероятности

без­

Однако

остается

неясным,

какова достоверность

этого

отказной

работы

P*(t);

интенсивности

показателя и насколько ему

отказов

X*(t) и

частоты

отказов a*(t)

 

 

от времени

 

 

можно доверять.

 

 

212

Т а б л и ц а 36

а

а

ь

С

a

а

ь

С

0,50

0

 

0,290

0,99

2,0206

1,40

1,073

0,75

0,5859

0,58

0,355

0,995

2,2373

1,61

1,250

0,90

1,1131

0,77

0,527

0,999

2,6841

2,09

1,672

0,95

1,4287

0,95

0,681

0,9995

2,8580

2,30

1,857

0,975

1,7023

1,14

0,846

 

 

 

 

Следует учитывать, что показатель надежности Р* (3000) =0,425 получен непосредственно по данным статистической выборки без каких-либо предположений о распределении времени возникнове­ ния отказов. Поэтому доверительный интервал для величины Р*(3000) =0,425 не может быть получен на основе правил, изложен­ ных в гл. IV с использованием распределения Стьюдента. В подоб­ ных случаях для суждения о степени достоверности количественной характеристики надежности можно воспользоваться ^-распределе­ нием Фишера и следующей формулой:

Я (0 = ------------ --------------,

(ХП.31)

п + 1

 

 

1 + --------------------

- V

 

7Va — п

a',2'Vl

 

где P ( t ) — надежность по результатам испытаний с наработкой t

часов

и при коэффициенте доверия

1 —а; п — число полученных

при испытании отказов; N0— количество испытываемых невосста-

навливаемых образцов;

— процентная точка ^-распределе­

ния с соответствующими степенями свободы V2 и vi

(понятие о сте­

пенях

свободы см. в гл. IV),

определяемыми из

соотношений

 

vz= 2ra+2;

(XII.32)

 

vi = 2V0 — 2га.

(ХП.ЗЗ)

В случае объема выборки V0> 30 процентили

/ц,;,,;,, могут

быть с достаточной точностью определены из соотношения:

 

v, ^

а b)~1'2 eg,

{XII.34)

где

h

2nV2

;

(XII.35)

 

 

•ц + v2

 

 

 

g =

v2 — V1

 

(XII.36)

 

-------- — •

 

 

vlv2

 

 

Параметры а, b и с в зависимости от а задаются табл. 36. Если а<0,50, то вычисления ведут следующим образом: для величины 1—а из табл. 36 выбирают а, b и с;

213

находят h и g; при вычислении g из формулы

(XII.36)

величины

у 2—Vi в числителе меняют местами;

 

 

из соотношения (XII.34) находят lg F \-a, Vi;

а затем

F

определяют

с-

1

 

 

 

v, = —------------ •

 

 

* 1—ajvjjv*

На основе изложенной методики оценим достоверность показа­ теля Р* (3000), вычисленного по данным табл. 35:

находим v2= 2 « - |-2 = 2-575--j-2-— 1152;

v1=27V0 —2ti = 2 - 1000 — 2-575 = 850;

зададимся коэффициентом доверия 1-—а=0,9, т. е.

 

 

 

 

а = 1 —0,9 = 0,1;

 

 

из табл. 36 для

1—а = 0,9 найдем а = 1,1131; Ь = 0,77 и с = 0,527;

вычислим

h

2у]У2

2-850-1152

 

 

'Ч + У2

 

977;

 

 

 

850 + 1152

 

 

найдем g =

 

=

850 ~

1152-= -0 ,0 0 0 3 0 8 ;

 

6

у , у 2

 

850-1152

 

 

вычислим lg F\

== 1,1131 (977 — 0,77)-1/2- f 0,527-0,000308 =

= 1,1131-------------1-0,000162=0,0358;

Л _ 01.Ч2=

1,086;

 

1 / 9 7 6 , 2 3

 

 

 

 

 

определим Fa^

t= — — ------- = ^ - ^ " = 0 ’92;

 

 

 

 

 

1—a;vA;v2

*

 

 

по формуле (XII.31) найдем Р*(3000):

 

 

 

Р* (3000)

 

1

0,45,

 

 

1+

575 + 1 0,92

 

 

 

 

 

 

 

 

 

 

 

 

1000-575

 

 

т. е. с коэффициентом доверия

1—а = 0,90

истинная

вероятность

безотказной

работы

элемента

(устройства)

в течение

3000 ч со­

ставляет не менее 0,45.

 

 

 

 

§ 31. Основные статистические модели, используемые в теории надежности

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

Частота отказов a(t) используется в оценке надежности гораздо чаще, чем средняя частота со(/). В большинстве практических слу-

214

чаев n(t)<^N0( см. формулы XII.20 и XII.25), и потому можно пре­ небречь незначительной погрешностью из-за уменьшения количе­ ства N0 испытываемых элементов (устройств) после каждого отка­ за. Поддержание постоянного числа элементов заменой отказавших может обусловить значительно большие погрешности вследствие неодинаковой наработки элементов и отличий характеристик эле­ ментов, вводимых взамен отказавших. К тому же вычисление вели­ чины ©(f) значительно сложнее, чем a(t).

а)

В)

6)

0

0,4

0,8

1,2

t_

Рис. 49. Графики зависимостей P(t);

X(t)

и a(t) для

различ­

ных законов распределения:

 

 

а — для нормального закона распределения

Т— Т, + С»а«^<,а

;

б — для

гамма-распределения; в — для распределения Вейбулла

 

В теории надежности часто используются для a(t)

следующие

распределения: экспоненциальное (см.

рис. 48 — на

участке от t0

до 11 ); нормальное (рис. 49,а); гамма-распределение

(рис. 49,6);

бета-распределение (см. рис. 13, 14 и 40);

Вейбулла

(рис. 49, в).

На рис. 48 была представлена характерная кривая для частоты

отказов a(t). Сравнивая рис. 48 и рис. 14,

можно заключить, что

215

для описания a(t) во всем диапазоне времени работы различных устройств (включая период приработки до t0 и период нарастания частоты износовых отказов при t > t \ ) в наибольшей мере пригодно бета-распределение. Однако, как было показано в гл. III, расчеты с использованием бета-распределения достаточно затруднительны. Кроме того, приработка выпускаемых элементов (устройств) часто производится предприятием-изготовителем. К моменту же време­ ни ti (см. рис. 48) устройство часто целесообразно снять с экс­ плуатации и заменить новым или отремонтированным. При таких условиях наибольшую значимость для количественной оценки на­ дежности имеет участок кривой a(t), соответствующий изменению t от t0 до t\. На этом участке в большинстве случаев частота отка­ зов a(t) меняется по экспоненциальной кривой. Вот почему экспо­ ненциальное распределение является в теории надежности основ­

ным. Рассмотрим его более подробно.

Частота отказов описывает­

Экспоненциальное распределение.

ся формулой

a(t) =

'ke

 

(XII.37)

 

 

 

где X— интенсивность отказов.

(XII.21) и (XII.37)

получим:

На основе формул (XII.18),

ч п =

а

 

А

е''>4

Ае- и

:A = const.

 

 

 

 

 

—\t

 

 

1 — | a (t) dt

i

-

j

X 4 te

 

 

о

 

о

 

 

 

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

Учитывая зависимость (XII. 18), получим

 

Э-ХО

(XII.38)

Q W = l

-It

(XII.39)

Математическое ожидание времени безотказной работы Т будет

равно:

 

 

T = ^ P { t ) d t =

e - wrff= — Lie-*<r = - L .

(XII.40)

о

X 1

X

v

 

 

 

 

Очевидно, что tcv = T. В теории вероятностей доказывается, что при экспоненциальном распределении среднеквадратичное откло­ нение случайной величины равно ее математическому ожиданию. В данном случае а(Т) = Т.

216

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

Нормальное распределение. Формула частоты отказов в этом случае имеет вид:

 

(<- ту

 

 

a (J)-

2о2

 

 

(ХН.4П

 

 

 

V 2л

 

 

 

 

где Т — математическое ожидание времени

безотказной

работы;

о2— дисперсия величины Т.

 

 

 

и - т

 

 

 

 

Согласно (XII.18) P{t) = l - \ a { t ) d t = -

. \

— - -е

2ча jj

at.

о

 

о

УТп

 

 

« V

 

Вводя под знаком интеграла новую переменную

 

t — T

■и\

 

 

(XII.42)

 

 

 

а ] / 2

t — T

о VT

получим Р (^ )= Ь — - — Г е~игс/и, что дает (см. гл. III):

V я —т VT

 

 

 

 

 

+

ZL_

(XII.43)

 

 

г

/

2

/

V

0 V 2

Так как обычно Т > о, то Ф ( -----— ) =

1,0. Тогда

 

 

 

а ] /

2

 

 

 

 

 

 

1 — Ф

t — T

 

(XII.44)

 

 

 

 

а V2

 

 

Внося (XII.41) и (XI 1.44) в (XI 1.21)

получим

 

 

 

 

 

 

t — T

\2

 

_

 

e x p

I

-

 

 

V 2

 

 

 

>1/2

 

(XII.45)

4 0 =

 

1 —Ф

t — T

 

а 1 /

я

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

0 V 2

 

 

Зависимость (XI 1.45) можно записать в виде:

 

х (0 =

l/o

 

 

 

 

 

— -?(0-

 

 

а | / Я

8—1092

217

На графике рис. 49, а представлена зависимость K(t). Как видно из графика, интенсивность отказов K(t) с течением времени увели­

чивается в отличие от

экспоненциального

распределения,

когда

л (/) = const. Рост K(t)

по мере увеличения

длительности

работы

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

Гамма-распределение. Частота отказов при этом распределении

выражается формулой

 

а (()

(XII.46)

Г)!

где Ко и к соответственно параметры масштаба и формы кривой плотности гамма-распределения, причем А о > 0 и к>0.

Используя зависимости (XII.18), (XII.21) и (XI 1.46), можно по­ лучить следующие формулы:

 

 

к— 1

^

 

 

 

Р(^) = е - Х°< V

-

(XII.47)

 

 

л5)

 

 

 

или

Р(/) = е~х°*

V с- \.t

(У )2

с-хч |

 

 

1

1-2

 

 

 

4 0 =

Ко(У )к 1

 

(XII.48)

 

к — 1

( W

 

 

 

 

 

 

 

п = О

п\

 

 

 

 

Т = — .

 

 

(XII.49)

 

 

Ко

 

 

 

При к 1 из (XII.46)

получаем

 

 

 

а(П = ).0е_Хо<

и Т = ——

, т. е. формулы для

экспоненциального

распределения [см. (XII.37) и (XII.40)].

Таким образом, гамма-распределение есть обобщение экспонен­ циального распределения для случая, когда отказ системы проис­ ходит при появлении ровно к независимых событий (отказов эле­ ментов). Это позволяет использовать гамма-распределение для ко­ личественной оценки надежности резервированных систем с одним действующим и к 1 резервными элементами (общее число элемен­ тов 1 + к—1 = к ). Каждый новый элемент включается в работу после отказа предыдущего, причем время безотказной работы каждого элемента распределено по экспоненциальномуУзакону. Отказ систе­ мы происходит после последовательных отказов всех к элементов. Гамма-распределение, как ясно из рис. 49, в, позволяет оценивать количественно надежность устройств и в период приработки их.

При к = 1 гамма-распределение переходит в экспоненциальное. Кроме того, частным случаем гамма-распределения (при к — целом

218

положительном числе) является распределение Эрланга, широко применяемое в теории массового обслуживания.

Бета-распределение. Это распределение было подробно рассмот­ рено в гл. III. Плотность его, а следовательно, и частота отказов выражается следующей формулой:

a гЛгГГ^ 1(1 (ХП'50)

г (Т) Г (1))

г д е 0 < / < 1 ; у > 0 ; Л > 0 .

Здесь весь интервал времени работы устройства принят за еди­ ницу. В соответствии с (XI 1.18) будем иметь:

t

P { t ) = 1 — j"a{t)dt

о

 

 

 

 

 

t

 

 

 

 

или

P { t ) =

1 -

— <T+.^-

ГfT-i(1 _

 

 

(XII.51)

 

 

 

Г (т)г Н)

J

 

 

 

 

где символ Г означает гамма-функцию

(см. гл. III. § 7)

от соответ­

ствующего аргумента, а интеграл в правой

части — неполную бе­

та-функцию [Bt (у, г))].

 

 

 

 

 

 

 

Интенсивнсть отказов h(t)

может

быть

вычислена

на основе

соотношения (XII.21).

 

 

 

 

 

 

 

Математическое ожидание времени безотказной работы в соот­

ветствии с (XII.28) определится:

 

 

 

 

 

T = [ p ( t ) d t = [ \

1 --------L 1 I ±

j !!2 _ С / т-

1 ( j _

ty -'d t

dt ,

J

 

J

г ^ г ^ )

,)

1

;

 

о

 

о

 

1

t

 

 

 

 

что дает

т== I

Г (тг+.Tj)

 

 

 

 

■§dt j

 

 

 

(XII.52)

 

 

г (7)г (1)

о

о

 

 

 

 

т. е. отыскание Т связано с интегрированием неполной бета-функ­ ции, что весьма затруднительно.

Частными случаями бета-распределения являются равномерное,

треугольное и параболическое распределения.

 

 

Так, при у = г) = 1

из (XII.50) получаем a(t) = 1,

т. е. плотность

равномерного распределения (см. рис. 14, в).

 

 

При у = 1 и 17 =

2

будем иметь:

a ( t ) = ---- ---------- (1— /).

. Так как Г(х) = (дс— 1 )-Г(дс— I) *,

 

 

Г (2) Г (1)

V

то Г (3)=2Г (2)

и

Г(1) = 1,0 и

далее a ( t ) —2(1—t),

т. е. плотность

треугольного

распределения

(см. рис. 14, в).

 

 

 

 

 

 

 

* Б р о н ш т е й н

И.

Н. и С е м е н д я е в

К-

А. Справочник по математи­

ке для инженеров и учащихся втузов. 3-е изд. М.,

Гостехтеориздат, 1953, с. 608.

8*

219

 

Соседние файлы в папке книги из ГПНТБ