Файл: Смирнов, О. Р. Надежность судовых энергетических установок.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 16.10.2024

Просмотров: 137

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

ГЛАВА IV

НАДЕЖНОСТЬ ЭЛЕМЕНТА СЭУ

§14. НАДЕЖНОСТЬ ПРОСТОГО ЭЛЕМЕНТА УСТАНОВКИ

В соответствии с классификацией отказов СЭУ (см. § 2) среди элементов, входящих в состав СЭУ, будем различать простые и слож­ ные. Под простым элементом будем понимать законченную конструк­ цию, предназначенную для выполнения определенной задачи, при­ чем ее отказы являются отказами какой-либо одной из трех ранее указанных групп. Под сложным элементом будем понимать элемент, в котором возможно возникновение отказов различных групп, в том числе и всех трех.

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

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

Пусть величины т и тв распределены показательно с параметрами Я и р соответственно. Тогда вероятность возникновения отказа или окончания восстановления в интервале времени At с точностью до бесконечно малых высшего порядка по сравнению с At равны КAt и рД^ соответственно, и соответствующие вероятности того, что эти события в период At не произойдут, будут: 1 — ХАt и 1 — рД^.

В произвольный момент времени t элемент установки может на­ ходиться в одном из двух состояний: исправном или состоянии отказа. Вероятности этих состояний обозначим соответственно через ёо (О и gi (t)- Рассмотрим момент времени t + At. Событие, заклю­ чающееся в том, что в этот момент времени элемент находится в ис­ правном состоянии, есть сумма двух несовместных событий: 1 ) в мо­ мент t элемент исправен и за At не отказал; 2) в момент t элемент на­ ходится в состоянии отказа, но за At восстановление окончено.

Таким образом, вероятность

g 0 (t +

At) с точностью

до беско­

нечно малых высшего порядка определится соотношением

go (t

+ At) = g 0 (t) (1 -

XAt)

+

gl (t) рДt.

(4.1)

Раскрыв в

(4.1) скобки, перенеся g 0

(t)

в левую часть,

разделив

обе части равенства на At и перейдя к пределу при At —>0, запишем lim go(/+ дУ ~ go {t) = — Xg0 (t) + pgy (0

9*

131


или

go(t) = — Xg0 (t) + Jig! (t).

(4.2)

Рассмотрим теперь событие, заключающееся в том, что в момент t + At элемент находится в состоянии отказа. Оно также представ­ ляет собой сумму двух несовместных событий: 1 ) в момент t элемент исправен, но за At произошел отказ; 2 ) в момент t элемент находится в состоянии отказа и за Аг восстановление не закончено.

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

gx (t + At) = g 0

(t) Ш + gx (t) (1 - jiA*).

(4.3)

Раскрыв скобки в (4.3),

перенеся g x (t) в левую часть,

разделив

на At и перейдя к пределу при А£-—>0 , аналогично предыдущему найдем

 

gi(t) = *go(t) — №(<)■

(4-4)

Уравнения,

аналогичные (4.2) и (4.4), носят в литературе различ­

ные названия:

уравнения Эрланга [85], уравнения процесса

гибели

и размножения

[20], уравнения марковского процесса [70].

Будем

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

Очевидно, что (4.2) и (4.4) не являются независимыми (правые

части отличаются лишь знаком) и что

 

М О + М 0 = 1,

(4-5)

поскольку в момент t элемент обязательно будет либо исправен, либо находится в состоянии отказа.

Присоединив к любому из уравнений (4.2) или (4.4) условие нор­ мировки (4.5), с учетом начальных условий g 0 (0) = 1, g x (0) = 0, имеем:

go ( 0 =

р +

*..

(4.6)

Я + р

 

 

Ы 0

X Хе~

1

(4.7)

Я. + р

 

 

 

 

Определим поведение найденных вероятностей для достаточно большого времени эксплуатации. Для этого перейдем в (4.6) и (4.7)

к

пределу при t —>оо: .j

 

 

 

 

 

lim g0 (/) =

lim

р + Хе~а+Ц) 1

Р

г

 

(4.8)

 

Я -р р

X 4- р

Т + Тв ~

 

 

со

f-»co

£ ° ;

 

 

Пш gl (t) =

lim

Л 1 о—(Л-+М-) t

А

 

(4-9)

 

*

 

= Т Т Г = Sv

не

Таким образом,

финальные значения вероятностей g 0

(t)

и g x (t)

зависят от времени.

 

 

 

 

132


Нетрудно

видеть,

что

g (t)

есть коэффициент готовности эле­

мента g о (t)

— k (t)

и что

его

поведение соответствует характеру

кривой, изображенной на рис. 29.

Представляет интерес исследование вопроса о скорости сходи­ мости величин g 0 (t) и g t (t) к своим предельным значениям. Из рас-

а)

0,9

0,8

fi

0,08

0,08

0,04

0,02

О

fo

1, 000г

0,990\-

Ю

20 40 ВО t, ч

 

 

 

 

20

40 ВО

80 t,4

STo(t)

1

go(t)

 

 

0,98

N . II .

II

Silt)

Ю ?

 

 

 

 

 

 

0,8

 

Рис. 36. Зависимости от времени вели­

0,8

 

чин g0 (t) и gx (t)

при различных зна­

 

 

чениях

К

и р:

а — Х =

10” 2 1/ч,

0,4

 

р = 1 0 - 1

1/ч;

б — Х =

1(Г4 1/ч,

 

 

 

р =

10~2

1/ч;

в Х = 1 0 '3

1/ч,

0,2

 

р =

1 0 '1

1/ч;

г — Х =

10~5

1/ч, '

 

р = 5 -1 0 ~ 2 1/ч; д — Х = р = 5Х

 

 

 

 

Х10‘ 2 1 /ч.

 

 

О

40 80 120 t, ч

 

 

 

 

 

 

смотрения (4.6) и (4.7) можно сделать лишь вывод об экспоненциаль­

ной скорости сходимости.

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

g 0 (t) и g x (t) от

На

рис.

36 изображены

времени при

значениях X от

5-10- 2 до 1 0 ~ 5 1 /ч и значений

р от 1 0 - 1

до 10- 2

1/ч.

Этот

интервал изменения величин

X и р соответствует

среднему времени

безотказной работы элемента

до

1 0 0 0 0 0

ч- и сред­

нему времени ремонта до 100 ч. Если учесть, что 10 000 ч — это при­

133


мерный ресурс многих вспомогательных механизмов, а 1 0 0 0 0 0 ч — предельное значение ресурса главных двигателей и что даже такая трудоемкая операция, как, например, демонтаж поршня малооборот­ ного главного двигателя и дальнейшая его установка на место после ремонта, занимает < = « 2 0 0 чел.-ч (такую операцию выполняют 2 — 5 чел.), то можно считать, что указанный диапазон значений к и р исчерпывает в основном практически возможные случаи.

Рисунок 36 показывает, что через 500—700 ч с точностью до 5% можно считать величины g0 (t) и gy (t) равными своим предельным значениям и пользоваться последними при расчетах, связанных с надежностью СЭУ.

Отметим также, что уравнения (4.2) и (4.4) описывают работу элемента с учетом его восстановления. Если такового не делается, то, полагая в (4.2) и (4.4) р, = 0, получаем уравнения, описывающие работу элемента до первого отказа. Решениями этих уравнений будут являться вероятность безотказной работы и вероятность отказа, которые были определены выше другим способом. Таким образом,, составляя уравнения вида (4.2) и (4.4), можно находить как комплекс­ ные характеристики надежности, так и характеристики безотказ­ ности. Однако здесь следует подчеркнуть, что изложенный метод рас­ чета характеристик надежности элемента может быть использован лишь в случае справедливости экспоненциального закона надеж­ ности.

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

Рассмотрим событие, заключающееся в исправности элемента в момент t. Оно представляет собой сумму двух несовместных собы­

тий:

1 ) элемент не отказал в промежутке (0 ,

t)\ 2 ) элемент отказал

в момент т < t, восстановился в момент т < 0

< t и далее в момент t

был

исправен. Предположив восстановление после отказа

полным

для

определения коэффициента готовности

получим интегральное,

уравнение

 

 

 

t t

 

 

 

k (t) — P (t) -(- J Ja(x) r ( 0 x)k (t 0 )d 0 dx.

(4.10)

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

ности элемента, работающего на переменных режимах.

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

(4.11)

о

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

134


промежутка времени х, для которого выполняется равенство

 

Qi(ti) = QAx),

(4.12)

где Q2 {х ) — вероятность

отказа второго элемента за

время х.

В соответствии с (4.11) и (4.12) этот промежуток

времени есть

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

 

 

h

X

(4.13)

J

ax(t) dt = | a2(t) dt.

оо

Нетрудно видеть, что аналогично может быть найдена величина л: и в случае, когда элементы 1 и 2 одинаковы, но каждый из них рабо­ тает на своем режиме ех и е2, соответственно. Здесь (0, х) есть про­ межуток времени, на котором вероятность отказа и, следовательно, вероятность безотказной работы при эксплуатации на режиме е2 равны этим величинам при работе на режиме ех за время tv

Пусть теперь имеется элемент, проработавший на режиме ех

время

tx и затем на режиме е

2

еще время t2. Найдем вероятность

отказа

элемента за время tx +

 

t2.

В соответствии с изложенным,

промежутку времени tx работы на

режиме ех соответствует, в смысле надежности, некоторый промежу­ ток времени х работы на втором режиме. Таким образом, промежу­

ток (О, t x + t2) работы сперва на режиме ех,

а затем на режиме е2

соответствует времени х

+

t2 работы на втором режиме, т. е.

Q (t 1 +

t2,

s j

, 8 2) =

Q (х +

t2, е2)

И Л И

 

 

X“f“12

 

 

 

 

(4.14)

Q(k +

t2) = j

a2(t)dt,

 

 

 

0

-

 

где величина x определяется из выражения (4.13).

Пример 1. Пусть время безотказной работы на режимах ех и е2 распределено показательно с параметрами Хх и Х2 соответственно:

a1{t) — 'kxe klt; a2 [t)= K 2e %г1 ■

Тогда в соответствии с (4.13)

х

=

к .

 

 

Х2 tl

 

 

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

 

 

 

 

 

 

 

X-\~t 2

 

 

 

 

Q ( ^ 4 - y =

j a2(t)dt =

1 — е - а л + М , ) .

(4.15)

 

 

о

 

 

 

 

В общем виде вероятность Q (tx +

t2) как функция времени может быть запи­

сана в виде

 

 

 

 

 

 

<2(0 =

1 — е- K i t

 

 

о < t < t x

(4.16)

|

g{Kltf^kztt

11)) t1< t < t i

 

135