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

книги из ГПНТБ / Рябинин И.А. Надежность судовых электроэнергетических систем и судового электрооборудования учебник

.pdf
Скачиваний:
25
Добавлен:
24.10.2023
Размер:
10.14 Mб
Скачать

i Понятие вероятности отказа в самой своей основе связано с по­ нятием частости отказов х . Частостью отказов Q* (t) в данной серии опытов называется отношение числа опытов, в которых наблю­ дались отказы п (/), к общему числу произведенных опытов т, т. е.

Q* (*) =

( З Л )

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

0 О,

О

Ю .213 ?7\30

40

50

60

70 80

90 100т

Рис. 15.

Верхние Qn

и нижние QH доверительные пределы

 

для

вероятности

отказа

Q.

 

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

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

Как было отмечено выше, при небольшом числе опытов т ча­ стость отказов Q* (I) является случайной величиной. Если к тому же и число наблюденных отказов п (t) невелико (в том числе и равно нулю), то частость отказов оказывается распределенной по биноми­ альному закону. Пользуясь методом доверительных интервалов, можно оценить число опытов т, при котором с доверительной ве­ роятностью б2 можно ожидать, что ошибка от замены вероятности отказа частостью не превзойдет заданного значения.

1 Частость иногда называют еще относительном частотой или просто частотой.

40

На рис. 15 построены

верхние QD и нижние QH

доверительные

пределы для вероятности

отказа Q 1 (с коэффициентом доверия

б 2 = 0,95) при числе отказов п и числе наблюдаемых

элементов т.

Из приведенных кривых видно, что понятие большого числа является относительным. Оно определяется заданной точностью (т. е. шириной доверительного интервала AJ6= Qn — QH), числом наблюдаемых

0,10

0

10

20

30

п

Рис. 16. График

изменения

ширины доверительного

интервала для вероятности

отказа в функции от числа

наблюдаемых отказов и при т = 100 и б 2 =

0,95.

отказов п и коэффициентом доверия б2 . Минимально допустимые значения для т при фиксированных значениях Д/б2 , п и б 2 при­ ведены в табл. 1.

На рис. 16 показан график изменения ширины доверительного

интервала

AJoB5

в функции

числа

наблюдаемых

отказов п при

5, = 0,95

и т =

100.

 

 

 

 

На основании

изложенного можно

заключить, что:

1) «достаточно

большие»

числа

не

являются

астрономически

большими;

 

 

 

 

 

 

2) малое число наблюдаемых отказов п не затрудняет, а облег­

чает оценку

вероятности

отказа;

 

 

 

 

 

 

 

 

Таблица 1

Минимально допустимое число

наблюдаемых

элементов /л,

 

которое можно считать практически большим

 

п

&Jo.n

= 0 ' 2 0

A J 0 . 0 5 = 0 ' 1

0

^0,9 5 = 0 ' 0 5

0

 

14

29

 

58

1

 

21

45

 

90

2

 

27

57

 

113

3

 

31

67

 

150

1 Для упрощения записи здесь опущено указание длительности времени испыта­ ний t.

41

3) для снижения «практически большого» числа опытов т до приемлемой величины не следует стремиться к необоснованно высо­ кой точности и слишком большой доверительной вероятности.

При ширине

доверительного интервала 0,10 и б 2 = 0,95

«доста­

точно большое»

число опытов (при малом п) не превышает

100.

0,31 0,31 0,33 0,34 0,95 0,33 0,37 0,33 В„

0,33 0,3320,391 0,333 0,33! („

Рис. 17. Графики минимального числа необходимых наблю­

дений /л^ в зависимости

от нижней

границы доверительной

вероятности безотказной

работы изделия

i?H для двух зна­

чений коэффициента доверия

6j.

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

для того

чтобы

подтвердить

с

определенной

достоверностью

бх ,

что надежность данных изделий

не ниже RH =

1 Q„, необходимо

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

 

На

рис. 17

представлены

графики /щ, = /

(R„) для 6Х =

0,95

и 0,99.

Из

рис. 17 хорошо видно,

как быстро возрастает число необ­

ходимых наблюдений /Пб, с увеличением надежности контролируе­ мого изделия.

42

Итак, чем надежнее изделие, тем труднее становится статистиче­ ская проверка гарантированного уровня надежности RH. Тем не менее ее приходится проводить, когда дело касается весьма ответ­ ственных систем. Так, например, на космическом корабле «Апол- лон-8» для проверки надежности пуска двигателей последние приш­

лось

в

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

3000 раз, что, как видно из

рис.

17,

гарантирует

с достоверностью

в 95% безотказную их ра­

боту

с вероятностью

не ниже 0,999.

 

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

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

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

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

43

 

Математические методы оценки показателей надежности изделий

 

по результатам их испытаний в отечественной литературе довольно

 

подробно изложены в книгах

[8, 14,18, 35]. Поэтому не будем оста­

 

навливаться

на данном

вопросе.

 

При оценке показателей надежности судового электрооборудова­

 

ния по результатам его

работы в реальных условиях необходимо

I

учитывать

следующие

особенности:

 

1) поток

отказов судового электрооборудования за все время

 

эксплуатации, как правило,

является нестационарным;

2)в связи с ограниченным числом объектов наблюдения п сравни­ тельно небольшой длительностью эксплуатации приходится иметь дело со всей совокупностью судов, а следовательно, необходимо учитывать, что число объектов наблюдения за время эксплуатации может изменяться;

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

4)почти все изделия судового электрооборудования при их

отказе либо ремонтируются, либо заменяются новыми из ЗИПа;

5)при отказе некоторых изделий вместе с ними по техническим причинам иногда снимаются и исправные изделия, проработавшие значительный срок;

6)у целого ряда элементов судового электрооборудования отказы могут возникнуть как при работе, так и в нерабочем состоянии (например, из-за снижения сопротивления изоляции).

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

Сформулируем

некоторые общие принципы

подхода к работе

со статистическим

материалом, собранным на

судах:

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

ственную информацию можно из него извлечь.

 

2.

Если данных мало и нужно объединить информацию из раз­

личных источников, то абсолютно необходимо проверить

однород­

ность

объединяемых выборок.

 

3.

При одной и той же точности и достоверности для

оценки

какой-либо фиксированной вероятности [например, Q (^),

R (^)]

требуется значительно меньший объем информации, чем для оценки любой из функций, характеризующих надежность в интервале времени [0, tx ].

44

4. Если данных достаточно, то в первую очередь целесообразно определить характер изменения интенсивности потока отказов со (t)

или

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

отказов

% (/),• в

том числе

проверить гипотезу

о стационарности

потока

отказов.

 

 

5.

Стремление

в первую

очередь

оценить

среднюю наработку

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

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

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

§ 9. Оценка вероятности отказа по частоте

Как мы видели выше, для получения многих характеристик надежности нужно знать вероятность отказа за какой-то промежуток

времени

[t, t-\- At]1.

 

Обозначим

эту

вероятность Q [t,

t +

At]

или, несколько короче, Q (t,

At).

Простейшая задача, с которой стал­

киваются

на

практике,

состоит

в

оценке

неизвестной

вероятности

отказа

Q (t,

At)

по наблюденной частоте

его

появления:

 

 

 

 

 

 

 

 

Q

4 t i A

t )

=

J l ( ^ l )

 

 

 

 

 

( 3 2 )

где

п

(t,

At)

— число

наблюденных

отказов

в

интервале

времени

 

 

 

 

It, t +

 

АЛ;

 

 

 

 

 

 

 

 

 

 

 

 

 

т

(t) — общее

число

изделий,

за

которыми

установлено

 

 

 

 

наблюдение,

в

 

момент

времени

t.

 

Q (t,

At)

 

В

качестве точечной

оценки

неизвестной

вероятности

во

всех

случаях

разумно

принимать

частоту

Q* (t,

At).

Однако

эта оценка в ряде случаев оказывается недостаточной для харак­ теристики надежности судового электрооборудования. Действи­ тельно, довольно часто за промежуток времени At отказов вообще можно не наблюдать, т. е. п (t, At) = 0, а это приводит к тому, что вероятность отказа за время At в соответствии с формулой (3.2) также следует считать равной нулю. При этом никак не учитывается

1

Промежуток времени

принято обозначать указанием либо

его границ [tit

li+i],

либо длины Ati = ti+1

— ti, либо одной границы н длины (t(,

At[). В данной

книге в основном используется последнее обозначение. При цифровой записи аргу­ менты /(• и Ati разделяются точкой с запятой.

45

общее число изделий т {t), за которыми установлено наблюдение. В тех же случаях, когда число п {t, At) =4= 0, но мало (что весьма характерно для судового электрооборудования), величина оценки на разных судах резко меняется и не может служить устойчивой характеристикой надежности элементов. Это приводит нас к выводу

о целесообразности использования для

оценки вероятности

отказа

за время At метода доверительных интервалов.

 

Пусть частота отказов за промежуток

времени It, t + At]

равна

Q* (t, At). Мы хотим оценить возможную ошибку и указать какие-то границы для неизвестной вероятности отказов Q (t, At).

Абсолютно достоверными границами для Q (t, At) являются числа 0 и 1. Указание всяких других границ сопряжено с риском совершить ошибку, вероятность которой будем называть уровнем значимости. Вероятность противоположного события, а именно вероятность не совершить ошибку, обычно называют доверительной вероятностью или коэффициентом доверия.

Выбор величины доверительной вероятности в значительной степени зависит от той цели, которую мы перед собой ставим. Жела­ ние лучше застраховать себя от возможной ошибки приводит обычно к выбору весьма больших доверительных вероятностей (порядка 0,99 и более). Однако следует иметь в виду, что всякая перестраховка в статистических исследованиях имеет своп отрицательные стороны, так как чем больше доверительная вероятность, тем шире границы для неизвестной вероятности. Опыт показывает, что выбор довери­

тельных

вероятностей

0,95

и даже 0,90

вполне достаточен для

прак­

тических

целей.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Двусторонним

 

доверительным

интервалом

для

вероятности

Q

(t,

At)

с коэффициентом

доверия,

не меньшим, чем

б, называется

слу­

чайный интервал

J6 =

IQI{

(t,

At),

Qa

(t,

At)],

концы которого (дове­

рительные границы) QH (t, At) <

QB (t, At) зависят только от исхо­

дов

наблюдения

[т. е. п

(t,

At)

и т (t)]

и для

любого

Q (t,

At)

 

 

 

 

 

Р {QH

(t,

At)

<

Q (t,

At)

<

QB

(t,

 

At)\ ^

6.

 

 

(3.3)

Верхним

[0, QB

{t,

At)]

и нижним

[QH (t, At),

1 ]

односторонними

интервалами

называются

такие случайные

интервалы

J6 ,

для

кото­

рых

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P

{0 <

Q (t,

At)

<

QB

(t,

At)] ^

6;

 

 

 

 

(3-4)

 

 

 

P\QH

(t, At)

<Q(t,

At)

<

1} ^

6.

 

 

 

 

(3.5)

Дадим еще одно истолкование

понятия

доверительного

интер­

вала для вероятности отказа Q (t,

At): его можно рассматривать как

интервал значений вероятности Q (t,

At),

совместимых

с

опытными

данными и не противоречащих

им. Если высказывается

 

утверждение,

что

неизвестная

вероятность

Q (t,

At)

заключена между

вероятно­

стями QH (t, At) и QB

(t,

At),

это не означает,

 

что данное

утвержде­

ние справедливо с вероятностью б

в каждом отдельном

случае

при

заданных

п

(t, At)

и т (f).

 

В любом случае

 

вероятность

Q (t,

At)

имеет определенное неслучайное, хотя и неизвестное нам

значение,

и

46

утверждение

QH

(t, At) <

Q (t,

At) < QB (t,

At)

либо справедливо,

либо ложно. Однако если имеется несколько

последовательностей

наблюдений,

то,

вычислив

для

каждой из

них

значения Q„ (t, At)

О» •'Qi

Рис. 18. Графическая иллюстрация двустороннего доверительного интервала Jc2 Для вероятности от­ каза Q.

и Qa (t, At),

можно утверждать,

что в 1006 случаях из ста

приведен­

ное выше

высказывание будет

справедливым. Поэтому

величину б

нужно

 

толковать

не как вероятность попадания точки

Q (t,

At)

в

интервал

J6,

 

а

как

ве­

 

 

 

 

 

 

 

роятность того, что случай­

 

 

 

 

 

 

 

ный

интервал

J6

 

накроет

 

 

 

 

 

 

 

точку

Q (t,

At)

 

(рис.

18).

 

 

 

 

 

 

 

 

На

 

рис.

 

19

 

построено

 

 

 

 

 

 

 

100x3 = 300 двусторонних

 

 

 

 

 

 

 

доверительных

интервалов

 

 

 

 

 

 

 

Л.95 =

[Qn. Qn

1 Для вероят­

 

 

 

 

 

 

 

ностей

отказа

Q,

равных

 

 

 

 

 

 

 

Qi

=

 

0,01;

Q2

 

=

0,14

и

 

 

 

 

 

 

 

Q3

= 0,40.

Для

 

их

по­

 

 

 

 

 

 

 

строения

были

 

использо­

 

 

 

 

 

 

 

ваны

результаты

статисти­

 

 

 

 

 

 

 

ческого

моделирования на

 

 

 

 

 

 

 

БЭСМ-2

 

соответствующих

 

 

 

 

 

 

 

вероятностей,

 

каждая

из

 

 

 

 

 

 

 

которых

оценивалась

по

 

 

 

 

 

 

 

т =

100

опытам.

 

Таким

 

 

 

 

 

 

 

образом,

всего

 

было

про­

 

 

 

 

 

 

 

ведено три серии по 10 000

 

 

 

 

 

 

 

опытов

 

в

каждой.

 

Как

 

 

 

 

 

 

 

видно из рис. 19, интер­

 

 

 

 

 

 

 

валы

 

/„,95=

 

[Q„>

Q b 1

Дей­

 

 

 

 

 

 

 

ствительно

являются

слу­

 

 

 

 

 

 

 

чайными (по длине и место­

 

 

 

 

 

 

 

положению),

но

более 95

 

 

 

 

 

 

 

из

 

100

 

каждой

се­

 

 

 

 

 

 

 

рии) накрывают

оценивае­

 

 

 

 

 

 

 

мую

неслучайную

вероят­

 

 

 

 

 

0,5

8 0,6

ность

Q.

 

 

 

 

 

 

 

 

Рис. 19. Триста вычисленных доверительных

 

Интересно

 

 

отметить

 

 

 

интервалов

для трех вероятностей отказа: Qx =

здесь

 

и

другие

особен­

 

=

0,01; Q2 =

0,14; Q3

=

0,40.

 

ности

построенных

интер

 

 

 

отказа

Q,

тем

более

валов.

Чем

 

меньше

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

асимметрично относительно нее расположены доверительные интер­ валы. Действительно, при Q3 = 0,4 примерно половина интер-

47

валов оказалась смещенной в сторону нижней доверительной гра­ ницы, а половина — в сторону верхней. При Q: = 0,01 все интер­ валы, накрывающие искомую вероятность, оказались смещенными вправо. Иначе говоря, истинное значение малой вероятности всегда ближе к нижней доверительной границе, чем к верхней. Это полезно знать при оценке малых вероятностей отказа, тем более что длина самого доверительного интервала в данном случае может изменяться в достаточно широких пределах (в несколько раз).

Описанное толкование доверительного интервала одинаково спра­ ведливо как для двустороннего, так и для одностороннего (верхнего или нижнего) доверительного интервала. Поэтому вероятность того, что случайный интервал Уб накроет точку Q (/, At), в формулах (3.3), (3.4) и (3.5) обозначена одной и той же буквой б. Однако при поль­

зовании статистическими таблицами часто

приходится

сталкиваться

с

различным обозначением доверительных

интервалов,

что

связано

с

использованием одних и тех же таблиц

для разных по

величине

коэффициентов доверия. Чтобы исключить возможные ошибки, будем

обозначать

доверительную вероятность одностороннего

интервала

 

 

6i = Р,

 

(3.6)

а уровень

его

значимости

 

 

 

 

Yi = 1 б]. =

1 — Р.

(3.7)

Вероятность

Yi характеризует

вероятность выхода

Q (t, At)

из доверительного интервала только за одну границу. При двусто­ роннем интервале всегда существует вероятность ошибки в обе

стороны:

доверительный

интервал / б может

оказаться

целиком

либо левее, либо правее (см. рис. 19) искомой

вероятности

Q (t,

At).

Приняв

для вероятности

односторонней ошибки обозначение

Yi

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

 

 

 

Y 2 = - 2 Y I -

 

 

 

(3.8)

Тогда доверительная вероятность двустороннего интервала

б 2 ока­

зывается связанной с величинами 8Х

и

Yi следующими соотноше­

ниями:

 

 

 

 

 

 

1 = 2P 1;

 

 

 

fig =

1 _ V a

= 1 _

 

= 26i —

(3.9)

 

 

 

61

=

l ± h . ;

 

 

(ЗЛО)

 

 

 

7 l

=

- L = A . .

 

 

(3.11)

 

Перейдем к вопросу о нахождении

доверительных

границ QH

(t,

At) и QB

(t,

At).

элементов

т (г),

установленных

на

судах

и

При

малом числе

находящихся под наблюдением в момент времени t, а также при малом числе отказов /г (t, At), которые удается наблюдать за про-

48

межуток времени

[t,

t + At], частота

отказов

Q* (/,

At) имеет

биномиальное

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

 

 

 

 

 

F(Q*

I Q) = F(ii

I m,Q) = P\N^n

I m,Q\ =

 

=

%P{N

= k\

m, Q\ =

£

C£Q*(1 -

Q)'"-*,

(3.12)

где

 

 

С,* — число сочетаний

из т элементов по k;

F (Q* | Q) =

F (п | /«, Q) — условная функция распределения слу­

 

 

 

 

чайной величины Q* (или N) при усло­

 

 

 

 

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

 

 

 

 

равна

Q,

а общее

число

наблюдае-

,мых элементов равно т.

Здесь и ниже для сокращения записи опущены временные аргу­ менты.

Исходя из биномиального распределения, английские стати­ стики Клоппер и Э. Пирсон указали правило, которое дает гарантию того, что вероятность непокрытия доверительным интервалом неиз­ вестной вероятности Q не превосходит 2yv Перефразируем это общее правило следующим образом. Пусть при работе т независимых изделий, имеющих постоянную вероятность отказа Q, за время At наблюдалось п отказов. Тогда в качестве верхней границы довери­ тельного интервала следует взять единственное решение уравне­ ния

t cm<# 0 - Q B r - * = vi,

(3.13)

а в качестве нижней доверительной

границы — единственное реше­

ние уравнения

 

 

I C * Q * ( l - Q , . ) m - * = Yi-

(3.14)

Точное решение уравнений (3.13)

и (3.14) описано в

[27]. Для

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

О

. у[100 У 1 о/о . 2(/г+1)]

'

 

, „ , , - .

V d ~

2/и —/i+0.5x[I00Yi %,2 (71+1)]

 

\ 0 - 1 0 )

п

х [ЮОбд, о/о, 2я]

 

- fi

 

27п — 7 г + 1 +0,5x11006! %,2п]

'

к

'

в которых

/ (Q, г)

 

(3.17)

 

 

— так называемая Q-процентная точка ^-распределения с г сте­ пенями свободы ( Q = 1007!% и л н Q— 1006;!%). Этот параметр при заданном числе степеней свободы г и уровне вероятности Q легко определяется по табл. 2.2а, приведенной в работе [5]1 .

1

Эта таблица приведена также в работе

[27].

4

И. А. Рябшшн *

49

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