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

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

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

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

Разумеется, часто время восстановления лучше описывается не показательным законом, а, например, законом Эрланга:

где f/e

(/) —

плотность

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

 

v •—• интенсивность восстановления;

 

 

с о

 

Г

(/г) =

| e~4k-ldt

— гамма-функция,

 

 

о

 

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

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

В нашем случае процесс изменения состояний системы будет марковским процессом с непрерывным временем и конечным мно­ жеством состояний: состояние Е0 — система работоспособна, состоя­ ние Ег—система неработоспособна.

Два состояния системы, рассматриваемые как случайные события, свяжем со случайной величиной Е, равной Е0, когда система работо­

способна, и

Ег, когда

система неработоспособна.

Тогда (рис. 56),

если система

в момент

t, находится в состоянии Е0,

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

наступления отказа за последующее время A t равна Q (At) = 1 е~х A t , а если система в момент t} находится в состоянии Ег, то вероятность восстановления за последующее время At равна V(At) — 1 • e~ v A t .

Фиксируя какой-либо произвольный момент времени t, мы можем

застать систему

либо в

состоянии Е0

с

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

перехода

за

время

[t, t + At]

в состояние Ег, равной Q (At),

либо в состоянии

Ег

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

в состояние Е0,

равной V

(At).

Определим

вероятность / > „ ( / + At)

того, что в промежутке времени

[t, t

+

At]

система будет работоспособна, т. е. будет находиться в состоянии

Е0.

Это событие может осуществляться двумя различными

способами:

1 )

в момент t система находится в состоянии Е0

и за

промежуток"

времени At сохраняет это состояние

(событие

А);

 

 

 

 

2)

в момент t система находится в состоянии Ег,

а за

промежуток

времени At восстанавливается, т. е. переходит в состояние Е0

(собы­

тие

В).

 

 

 

 

 

 

 

 

 

13 И. А . Рябимин

193

По теореме сложения вероятностей несовместных событий полу­ чаем

Р0 (t + At) = Р (А) + Р (В).

(7.2)

Свойство отсутствия последействия, которое в данном случае имеет место, поскольку мы условились, что поток отказов является простейшим, позволяет, используя теорему умножения вероятностей, вычислить вероятность каждого из событий А и В. Так, вероятность события А равна

 

 

 

Р (А)

=

Р0 (t) Рю (д 0>

 

 

 

(7-3)

где

Р0 (t) — вероятность застать систему в момент t в состоянии Е0;

Р00

(At)

— вероятность того, что за время

Af система

останется

 

 

в состоянии

Е0.

 

 

 

 

 

 

 

 

 

ш

 

1

 

г

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/

 

 

 

 

 

 

 

 

f

 

 

0

 

т

 

 

:

t

 

 

 

 

 

 

ч

 

 

 

 

 

 

к№,

П)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-

/

/ /

 

 

 

 

 

 

 

 

 

/

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/

 

 

1

 

 

 

 

 

 

О

}

 

 

1

 

 

 

 

 

 

/

 

 

j v a o

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-At

 

 

•At

 

 

t

 

 

Р и с . 56. Г р а ф и ч е с к а я

и л л ю с т р а ц и я

к о п р е д е л е н и ю

п е р е х о д н ы х

 

 

 

в е р о я т н о с т е й системы .

 

 

 

 

Вероятность того, что за время At

не возникает ни одного отказа,

в соответствии с законом Пуассона

равна

 

 

 

 

 

 

Р00

(At) =

 

е~% A t = e-XAt=\—XAt

 

+

o (At)

(7.4)

[здесь через о (At) обозначена величина, бесконечно малая по срав­ нению с At], Следовательно,

P(A)^P0(t)(l~XAt).

(7.5)

Вероятность~события

В

равна

 

Р (В)

=

Рг (t) Р 1 0 (АО,

(7-6)

194

где

Рх

(t)

— вероятность застать систему в момент t в состоянии £ \ ;

 

Р10

(At)

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

того,

что за время А? система

перейдет

 

 

 

 

из

состояния Ег

в

состояние

Е0:

 

 

 

 

 

р10

(At) = V (At) =l—e-vAt

 

=

vAt

+

o

(At).

(7.7)

Следовательно,

 

Р

(В)

Рг

 

(t) v

At.

 

 

 

 

(7.8)

Отсюда

 

 

 

 

 

 

 

 

Р0

(t +

ДО =

Р0

(О Р 0

 

(At)

+

Рг

(t)

P1Q

(At)«

 

 

0

 

 

 

 

 

Р0 (t)

(1 — К At)

+ Ру

(t) v

At.

 

 

(7.9)

 

Перенося Р0

(t) в левую часть, деля

на

и переходя

к пределу

при

At—>0, получим дифференциальное

уравнение для Р0

(t):

*Ejdu. = -XPo{t)

+ vP1(t).

(7Л0)

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

Pi (*):

-?EjP- = XP0(t)-vP1(t). (7.П)

1-Ш

ш

Рис. 57. Схема возможных измене­

ний состояний системы.

Схема возможных переходов исследуемой системы из одного состояния в другое дана на рис. 57.

Таким образом, для определения вероятностей состояний системы Ро (0 и Рг (t) в произвольный момент времени t мы имеем систему дифференциальных уравнений

Pu(t) =

--kP0(t)

+ vPx

(0;

P[(t)

=

 

(7.12)

KP0(t)-vP1(t).

Эта система при начальных условиях выполнении требования PQ (t) + Рг (0 = ние:

V

p0(t) = X + v

X

Pl(t)-- + v

X

Р0 (0) =

1, Рх (0)

= 0 и

1 имеет

следующее

реше­

 

 

(7.13)

Величина P0 (t) представляет собой не что иное, как вероятность застать систему исправной в любой произвольный момент времени t, т. е. функцию готовности КГ (t), которую и требовалось найти:

 

* г ( 9 = Ро(*) =

Т ^ 7 7 [1 + Pe~iX+V)

(7-1 4 )

где

 

1 +

р

 

 

 

 

 

 

Р

А

А>

(7.15)

 

V

 

— коэффициент

неисправности.

 

 

13*

195

Зависимость функции готовности Кг (/) от значения коэффициента неисправности р показана на рис. 58. При увеличении р надежность

системы

снижается.

 

 

Рассмотрим частные случаи формулы

(7.14):

 

1) р =

О, X = const. Так как здесь Ть

= 0, то этот случаи

может

быть интерпретирован как автоматическое резервирование,

когда

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

новый

элемент (пли мгновенно устраняется состояние отказа). Из

(7.14)

получаем

Кг (/) = 1.

 

 

2) р = оо, ?i = const. Так как здесь Тп = оо, то этот случай соответствует отсутствию восстановления. Из (7.14) получаем Кг (t) =

 

 

О

tc

 

t

 

 

Рис. 58. Зависимость функции готовности Кг

(0 от значения

 

 

 

коэффициента неисправности р.

 

=

е~,

т. е. формулу для оценки

надежности

невосстанавливаемой

системы

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

работы);

 

 

3) р = const и X = const. В этом случае существует такой момент

времени

tc,

после которого устанавливается стационарное значение

Кг

(tc) = const, не зависящее от времени и равное

Величина Кг представляет собой постоянную стационарную вероятность того, что система будет находиться в исправном состо­ янии в любой момент времени после tc, т. е. коэффициент готовности системы (2.38). Из выражения (7.16) видно, что можно достичь высокого уровня готовности систем даже при низкой исходной их надежности, для чего необходимо добиться условия

о

Вероятность непрерывной безотказной работы системы в течение времени т, начиная с момента t, в условиях рассматриваемой задачи

196

легко найти по формуле экспоненциального

закона надежности:

 

R(x)

= e-Kx.

(7.18)

Таким

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

того, что система будет исправна

в момент

/ и будет дальше безотказно работать

в течение времени т,

равна

 

 

 

 

Я &

^) = T T 7 [ 1 + P e - ( M - V ) ' ] e " u .

(7.19)

Если t достаточно

велико

(t >

/ с ), то

 

 

R(oo,

х) =

^ — е - х х .

(7.20)

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

R(t, х) = е~м е~Кх =

а + т ) .

(7.21)

Эффективность восстановления можно количественно оценить коэффициентом Э (/), равным отношению вероятностей, определя­ емых по формулам (7.19) и (7.21):

3(t) = T ^ [ e X

'

+ pe-vt\

(7.22)

или при t >> tc

 

 

 

5 ( c o ) = T T 7 e W

=

X T V e 4

(7.23)

Как и следовало ожидать, эта эффективность возрастает с увеличе­ нием интенсивности восстановления v, интенсивности отказов X и времени t.

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

Для составления уравнений, описывающих поведение такой системы, воспользуемся допущениями о справедливости марковского случайного процесса, о равенстве параметров Хх — Х2 — X; \>г =

=v 2 = v и о невозможности отказов у неработающих генераторов. Рассмотрим четыре возможных варианта этой задачи:

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

происходит без ограничений); 2) оба генератора находятся в действии, а в случае их отказа

ремонтируются по одному (ограниченное восстановление);

371

197

3)один генератор находится в действии, второй—в ненагруженном резерве; восстановление происходит без ограничений (как в первом варианте);

4)то же, что п в третьем варианте, но восстановление ограни­ ченное.

Дифференциальные уравнения, описывающие возможные состоя­ ния исследуемой системы ( £ 0 , Ех и Е2) для названных четырех вариантов задачи, имеют следующий вид:

для первого варианта

 

 

 

 

Р0

(t) = -2kP0(f)

+

vP,(/);

 

 

 

 

 

 

Р\ (0 =

2кРй (I) -

-'- v) Pi (0 -|- 2vP2 (0;

 

(7.24)

 

 

 

 

 

Pi(t)

=

XPl

(l)-2vPz(l);

 

 

 

 

для

 

второго

варианта

 

 

 

 

 

 

 

 

 

 

 

 

P'0(l) =

— 2\PQ(t)

-\- vPi (t)\

 

 

 

 

 

 

P[ (t) =

2kPD

(t) -

(k -\

v) P, (0 + vP, (/);

 

(7.25)

 

 

 

 

 

 

Po(t)

=

 

kP1(t)-vP,(t);

 

 

 

 

для

 

третьего

варианта

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ед

= _ Я Р 0 ( / ) + г Л ( 0 ;

 

 

(7.26)

 

 

 

Pi (t) =

кРй (0 -

+ v) Pi (/) +

vP2

(/);

 

 

 

 

 

 

 

 

 

 

 

P2 (0

=

^Pi (0 2vP2 (0;

 

 

 

 

для

 

четвертого

варианта

 

 

 

 

 

 

 

 

 

 

 

P'o(t) =

bP0(t)+vP1{i);

 

 

 

 

 

 

Pi

(t) = kP0

(t) -

+

v) P,(0 + vP2 (/);

 

(7.27)

 

 

 

 

 

 

 

 

 

 

 

P2(/) =

^ P i ( ^ ) - v P 2 ( 0 ,

 

 

 

где

P 0

(I) — вероятность

состояния E0,

при котором

отсутствуют

 

 

 

отказавшие генераторы;

 

 

 

 

 

 

Pi(t)

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

состояния Ех,

при котором

имеется один

 

 

 

отказавший

генератор;

 

 

 

 

 

 

Р 2

(t) — вероятность состояния Е2,

при котором оба генератора

 

 

 

отказали.

 

 

 

 

 

 

 

 

 

 

Интегрирование уравнений (7.24)—(7.27) при начальных усло­

виях

Р„ (0) = 1; Рг

(0) =

Р 2 (0) = О

и

выполнении

требования

Р0

(t)

+

Р х (I) +

Р 2

(0 =

1

дает

зависимость

Рк [i)

для

любого

к (к = 0 , 1 , 2). Указанные уравнения

можно

решать

при

помощи

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

Лапласа.

С

учетом принятых

условий имеем

 

 

 

 

KrV)

= P*{t) + Pi(t)=

 

1-РЛ*)-

 

(7.28)

198

В табл. 32 представлены функции готовности К г 1 (I) K r i (t) для четырех вариантов задачи, а на рис. 59 они представлены в гра­ фическом виде для следующих значений параметров:

% = 0,002 1/ч (Т0 = 500 ч), v = 0,01 1/ч (Тв = 100 ч), р = 0,2.

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

1,00КгW

0,99

 

 

/

 

Kfilt)

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

0,98

 

 

 

 

 

 

0,97

 

 

 

 

 

 

 

0,9В

 

 

 

 

Kn\t)

 

 

 

 

 

 

 

 

 

0,95]

 

 

 

 

 

 

 

0,9^

 

 

 

 

 

 

 

0,93

 

 

 

 

 

 

 

ОД

 

 

 

 

 

 

 

0

50

100

150

200

250

300

t.4

Рис.59.

Кривые изменений функции

готовности Кг

(0, '"построен­

 

 

ные по данным табл. ЗвТ"

 

 

а.

Стационарные значения коэффициентов готовности исследуемых систем равны

К л = 0,9722, Д'г 2 = 0,9459, /<г 3 = 0,9836, K T i = 0,9677.

При уменьшении коэффициента неисправности р коэффициент готовности системы Л"г возрастает. Уже при р — 0,1 соответствующие К г будут равны 0,9917; 0,9836; 0,9954; 0,9909.

При составлении дифференциальных уравнений для определения функции надежности системы R (t) следует «запретить» переходы системы из состояния Е% в состояние Ех, т. е. рассмотреть марков­ ский процесс, у которого состояние Е2 является поглощающим. Дифференциальные уравнения в этом случае примут вид:

для первого и второго вариантов

• P0(t) =

-2XP0(i)-]-vPl(t);

 

Р[ (t) = 2XP0(t)-(X

+

v)Pl(l);

(7.29)

P'2(i) =

XPl(t)-

 

 

199

 

 

Таблица функций

готовности

 

В а р и а н т

Ф у н к ц и и г о т о в н о с т и

д у б л и р о в а н н о й

К о э ф ф и ц и е н т ы

готов ­

з а д а ч и

системы K r i

(t)

ности систем

К ^

Я3

1 —

2

(X + v)* + X2

Я2

(A, + v ) 2 - f v2

X2

(Я + v)2 —Яv

1

1 + V~№ + 4Яv( z ^ ' - z ^ )

/ V 2 +

4 ^ V

+ 1 П 7 й Г ( г " " ' - * " г " )

1 + 2 р ( 1 + Р ) 2

1 + 2 p

(1 + p ) » + p*

2(1 + P ) (1 - Ь Р ) 2 + 1

1 + p

( 1 - Ь Р ) 2 - Р

Таблица 32

Ко р н и х а р а к т е р и с т и ч е с к о г о

ур а в н е н и я

z1

=

-(X

+

v)

z2

=

— 2

(Я. + v)

ЗЯ-1 - 2v— V%2 + 4Ъ>

ЗА, + 2v +

J / ^ 2 + 4A,v

2Я + 3v — Vv2 + 4A.V

~2

2Я + 3v + Vv2 + 4Xv

2i = — (Я-f-v — j/"A.v)

г2 = - (Я + v +

для третьего

и четвертого

вариантов

 

 

 

P'0(t)

=

—JJ>o(t)+vPl(t);

 

 

P'i(t) = M3o(t)-(b

+

v)Pl(t);

(7.30)

 

 

p'2(t)

= XPl(t).

 

 

При начальных

условиях

Р0

(0) =

1 и Рх

(0) = Р2

(0) = 0 вероят­

ность безотказной работы дублированной системы при нагруженном

резервировании (первый

и второй

варианты)

будет определяться ^

выражением

 

# 1 . 2 ( 9 = -

2 еа >

 

•«1

 

 

 

(7.31)

где

 

•0,1

_ -..(1

+ Зр)н=Кр' +

6р-

А 1 , 2 —

 

 

(7.32)

а при ненагруженном резервировании (третий и четвертый вари­ анты) — выражением

1

Р 2 - Р

(7.33)

где

(1 + 2 р ) т К 4 р + 1

Pi.3 =

-

(7.34)

0,5\

0 i f

«3

о/

в/

л' да т боо

Рис. 60. Графическая иллюстрация функции надежности электростанции.

Для принятых выше в качестве примера параметров, получим следующие численные выражения этих функций:

 

R1<2

(t) = Q 1

(0,015485e-Q 0051 s< -

0,000515е-°Л 1 5 4 8 5 ( );

 

 

R3A

(0 =

- о 1 1 Г з Г ( 0 , 0 1 3 7 е - о ^ — О.ОООЗе-0 .0137 ').

Графически они показаны на рис. 60. Здесь же приведены для

сравнения

еще две функции надежности:

 

 

1)

R^

(t) — для дублированной

системы

без

восстановления;

2)

R(1)

(t)—для

случая одного

нерезервированного генератора.

Из рис. 60 наглядно видна роль резервирования,

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

и допущения

о невозможности отказов

у

неработающих генера­

торов.

 

 

 

 

 

 

 

 

201

В случае различных параметров Хг =/= Х2 и vx ф v 2 описание дублированной системы намного усложняется, хотя методика в прин­

 

 

 

 

 

 

ципе остается той же самой.

 

 

 

 

 

 

 

 

 

 

Рассмотрим для примера

нагружен­

 

 

 

 

 

 

ное

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

восстановлением,

 

 

 

 

 

 

которое

начинается

немедленно

 

после

 

 

 

 

 

 

отказа основного ( я ^ и л и резервного (х2 )

 

 

 

 

 

 

элемента

(генератора).

 

 

 

 

 

 

 

 

 

 

 

Диаграмма

возможных

изменений

 

 

 

 

 

 

состояний

такой

системы

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

 

 

 

 

 

 

на

рис

61.

 

Как

видно

из

рисунка,

 

 

 

 

 

 

в этом случае число состояний дубли­

 

 

 

 

 

 

рованной

 

системы

равно

четырем

0,

Рис. 61. Схема возможных изме­

Ег,

 

Е2

п Е3),

 

причем отказовым (погло­

нений состояния системы.

 

 

щающим)

 

состоянием

является

только

 

 

 

 

 

 

состояние

Е3.

 

 

 

 

 

 

 

 

 

Запишем систему дифференциальных уравнений для определения

функции надежности

R

(t):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ро = -

+

К) Ро (t)

+

V j P i

(0 +

V2P2 (t);

 

 

 

 

 

 

-Pi (0 = № ( 0 - ( ^ + v i ) / M 0 ;

 

 

 

 

(7.35)

 

 

P2(t) = bP0{t)-(h

 

+

 

 

Vi)P2(t);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Pz(t)

=

hPi

(0 +

hPi

(t).

 

 

 

 

 

 

 

Применяя преобразование Лапласа к системе уравнений (7.35),

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

zPQ(z)

-

Р0 (0) =

-х

+

л,) Р0 (2) -(- vА

(2) + v2 P2 ( г

) ;

|

 

 

 

(z) -

Рх

(0) =

ХгР0

(2) -

 

2 + Vl)

Рх (г);

 

 

}(7.36)

 

zP2(z)-P2(0)

 

= Х2Р0

(2) -

 

х

+

V

2 ) Ра

(2)-

 

 

 

 

 

 

 

 

 

 

 

zP3(z)

-

Р3

(0) = А А

(2) ХР2

 

(2),

 

 

 

 

где Р А (z) =

J е _ 2 ( PA

(t) dt — изображение

по

Лапласу

вероят­

 

 

 

 

 

 

 

ности

 

Pk(t);

 

 

 

 

 

 

 

 

P'k (0~г*- zPk (z) — P/j (0) — преобразование

 

Лапласа

для

про-

 

 

 

 

 

 

изводной от функции Р

(t).

 

 

 

Перепишем (7.36) в следующем виде:

 

 

 

 

 

 

 

 

 

 

(2 +

Хх + Х2) Р0

(2) -

v A i (2) -

(2) =

Р0 (0);

 

 

 

 

-

КРо (г) +

(2 +

Х2

+

V l ) Р х

(2) -

 

Р х

(0);

 

 

(7.37)

 

-

* А

(г) +

(г -Ь Ях

+

v2 ) Р 2

(г) -

 

Р2

(0);

 

 

 

 

 

 

 

 

 

 

- > A i ( z ) - Х Х Р 2 ( 2 )

+

A ( z ) =

Р3

(0).

 

 

 

 

202

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