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

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

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

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

 

с о

 

Tn.s =

\tUsWx.

(7.85)

 

о

 

Подставив в данное уравнение значение /„,s (т) и

выполнив ин­

тегрирование по частям, получим

 

f „ , s

= 4 ^ - .

(7.86)

Таким образом, имеется принципиальная возможность вычис­

ления математических

ожиданий THiS

для

любых

п и 0 s s=c «.

{j

Y<

Y3

Yj

Yn.t

Yn

Yj

 

Ч^]ЩНЗ

 

 

23.

 

li

Z2

En-f

In.

 

0

 

Гц

 

 

 

 

 

 

 

 

 

 

 

 

 

 

hi

 

 

 

 

 

 

 

 

 

Xi

 

 

 

 

 

 

 

T,

 

 

 

 

 

Хзт

 

 

 

 

 

Рис. 63.

Условия эквивалентирования:

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

 

 

цепи; б—параллельных

ветвей звена

схемы ненадежности.

 

Располагая формулами для

вычисления

T„iS,

напомним содержание

стоящих перед

нами

задач:

Yk и Zk элементов, входящих в

 

1)

зная средние

времена

цепь

звена

схемы ненадежности системы, требуется

определить эти

пара­

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

различных

элементов (рис. 63, а);

 

2)

зная

 

средние времена Y3 и 1Э для т

параллельных ветвей

звена

схемы

ненадежности системы, требуется

определить эти пара-

213

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

На основании анализа модели совпадения импульсов можно утверждать, что решениями этих задач по эквнвалентироваиию

последовательных

цепочек

и параллельных

ветвей схемы

ненадеж­

ности будут

 

 

 

,

 

 

 

У9 = Т„,0;

 

(7.87)

 

 

Z 3 = t f l u k ;

 

(7.88)

 

 

k=i

 

 

 

 

 

Z3 = fm. m ;

 

(7.89)

 

 

Уз = Т м . 0

- 2 3 ,

 

(7.90)

где Г м .0 — с р е д н е е

время

между отказами

звена.

 

Действительно,

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

ожидание времени

исправной

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

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

совпадений Тп^-

Если

в период ремонта одного из элементов

возможны отказы у других,

т. е. возможны совпадения ремонтов по два, по три и т. д., то среднее

время

восстановления цепочки

(обозначим

его в этом случае Z{31))

будет

определяться суммой всех

Tnjr

 

Расчетная

формула для определения Уэ

имеет вид

 

 

V

-J-

S

 

а соответствующая интенсивность отказов элемента хэ будет равна

K = -=-=th.

(7.92)

Расчетное выражение для определения Z3 при произвольном п является весьма громоздким. В связи с этим запишем формулы только для п = 2 и 3:

I f f = То, +

Г 2 2

=

_

Z l Y J +

ZJYx

_ +

- ^ g - ;

(7.93)

 

 

 

Z1 + Z2 + Y1-l-Y2

Zi + Za

 

Z(33

=

Tz,l

+

7*3,2 +

7*3,3

=

 

 

214

 

ZxYtYa

- j - Z2YXY3

+ Z3YXY2

 

 

 

Z i (K s + / , ) + Z 2

+ К,) + Z 3 (Yx +

Уя ) + r 2 r 3 +

K X K 3

+

К 2 Г 2

 

YXZ2Z3

- j - Y2ZXZ3 -\- Y3ZXZ2

 

 

ZXZ2 +

Yi (7-i + Z3) + Y2

(Zx + Z3) + Ya (Zi + Z2 ) + Z 2 Z 3

+ ад,

+

+

(7.94)

 

•^2^3 + ZjZ 3 - f - Z 1 Z 2

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

 

 

ZkTlYt

 

2™ = т„

 

 

(7.95)

П-1

_

 

 

%YtYj...Yi

 

i. j . 1фк

 

 

 

k=i

 

 

 

Соответствующая интенсивность

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

будет равна

 

 

 

k=i

k=i

\t=i

 

(7.96)

Z<2>

ft=l

где pk

z^

Yk

Исходя из понятия коэффициента готовности (или простоя) последовательной цепочки из п элементов и допуская совпадение ремонтов, можно получить следующее простое выражение для определения Z 3 (математического ожидания длительности восста­ новления цепочки):

п

 

 

П ( 1 +

Р Й ) - 1

 

2<з) _ k=i

.

(7.97)

Можно показать, что

7(2)

(7.98)

 

215

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

Расчетная формула для определения Z3 имеет вид

2з = - 5 Г

=

- , ! — >

(7-99)

V

J-

S v e /

 

/ = 1

 

 

 

а соответствующая интенсивность

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

элемента х3

будет равна

HI

 

 

 

 

(7.100)

v3

= Sv9 / .

Математическое ожидание времени исправной работы звена Y3 в формуле (7.90) вычисляется с помощью среднего времени между отказами звена Г м .0 , обратно пропорционального средней частоте

следования

импульсов совпадения

r i m

i m :

 

 

 

Л

 

Т~Г Н-Рэ/

Т

U(Y3i +

z3i)

1 J — р — -

P-m, m

 

 

л

м . о

 

 

2 > 9 ,

 

 

 

 

 

/ = 1

где Рэ1 = -^-

 

 

'Э1

звена будет

равна

Интенсивность отказа

 

m

- " m '

+ Pal _ j'

 

7=r

Рэ/

 

 

(7.Ю1)

(7.102)

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

звена

определяется

выражением

9 э ^

= ^-

= - ,

 

<7Л03)

Y3

V3

pj 1 +

рз/ _ 1

 

 

 

/ = 1

Р

э >

 

Если плотности распределения случайных величин Yk и Zk подчиняются экспоненциальному закону (обычно принимаемое до-

216

лущение в теории массового обслуживания), то суммарный поток восстановлений будет простейшим. На этом основании вероятность безотказной работы t'-го звена схемы ненадежности равна

 

Я3 / (*) = е х р ( - Ц * ) ,

(7.104)

а для всей схемы

ненадежности в соответствии

с (7.81) получаем

 

г

 

 

/?е .„(0 = П в - Ч ' .

(7.Ю5)

Последнюю вероятность удобнее вычислять из соотношения

 

1 п # с . н ( 9 = - * £ Ч -

( 7 - 1 0 6 )

 

1=1

 

Коэффициент

готовности можно определить

по формуле

=

1

~

7 ,

(7.107)

П(1 + р3 .)

I + SP3 .

 

(=i

 

1

i=i '

 

г

которая при малой величине 2 р 3 . может быть еще упрощена:

tfr^l-£p3..

(7.108)

i=l

Пример 26. Оценим надежность рассмотренной выше системы (пример 25) с учетом восстановления, приняв следующие параметры для элементов

П = 10,000 ч, Z f t = 100 4 , рА = 0,01 (fe= 1, 2, ... . 8).

Для того чтобы в дальнейшем можно было сравнить надежность вос­

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

системы с надежностью системы без учета восста­

новления (для элементов которой ранее

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

Rk = 0,9), расчет

произведем для времени

t — 1050 ч, так как

 

Rk

(1050) = е х р (

-

^ ) =

0,9003.

В соответствии с рис. 64 и формулой

(7.92)

находим

= К*. = 2А* = 0,0002 1/ч;

Yx,Xt

= YXlX, = 5000 ч;

Км.

= SXk = 0,0003 1/ч,

YXtXtXt

= 3333,34 ч,

217

По формуле

(7.93)

определяем

 

 

 

Z[\\

=

Z[\

-

3 7 * _

+ J - =

JkY,l

+ 0.5Z,

-

1 0 0 ' 1

0 0

0

• +

0,5-100 =

i™<™ +

50 = 90,9 + 50 =

140,9 ч.

100+

1000

'

'

 

1100

 

 

 

Для цепочки

из трех элементов по формуле (7.94)

имеем

 

~ ( 1 )

_

 

 

3-100-ю ооо-ю ооо

._

 

6х.х<х,—

з.юо(10000+ 10000)+ 3-10000-10000 1

 

 

 

 

_.

3-10000-100-100

._

 

 

 

 

 

'

3-10000(100+ 100)+ 3-100-100 '

 

+

1 з°100°100° = 9 8

- 1 + 4 9 > 7 + 3 3 - 3 =

181,1 ч.

 

Зная параметры всех эквивалентных цепей, перейдем к расчетам

звеньев.

По

формуле

(7.99)

находим

 

 

 

 

 

 

 

Z3,_ =Z3i

=

 

 

1

 

 

1

 

 

 

 

 

1

_ j _

l

 

 

1

, 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ZSl

'

ZXtXt

 

100

1

140,9

 

 

 

 

 

 

1

 

 

1

 

58,4 ч;

 

 

 

 

0,01 + 0,0071

 

0,0171

 

 

 

 

 

 

 

 

 

 

 

^з,

=

;

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~ZXi

 

 

Z W

t

 

100 1

181,1

 

 

 

 

0,01 +

1

0,0055

= - nnW = 6

4 > 5 4 ;

 

 

 

 

 

 

0,0155

 

 

 

 

2 3 o = Z3, = _ ! _ +

 

1

 

 

 

'

1

 

 

 

_ J _ + _ L

_ L , ^ , _ 1 _

 

 

 

ZXl

^

 

Z,„ ^

Zx.

100+

100 ^

100

 

 

 

 

 

 

 

1

 

33,3 ч.

 

 

 

 

 

 

 

 

 

 

 

0,03

 

 

 

 

 

 

По формуле (7.101)

определяем

 

 

 

 

 

 

71

 

~f

 

 

 

(Yxt

+

Zx,)

(Vx.Xj +

Zx.x4)

 

 

1

М.0З2

1 м-оЗ»

 

 

 

 

- ~n

 

 

 

_

(10 ОООЧ- 100) (5000+ 140,9)

ZA.-, + j£x2xt

 

_

9 ,R n nn

_ 10 100-5140,9

~

 

100+ 140,9

 

 

_

240T9

 

/

1 D

U U U 4 ;

 

 

f<

 

(Yx3 +

ZXa) (Yxax,x,

+ Zx,xtx,)

 

 

 

 

' ( l , 0 3 l

 

 

 

"v

r "v

 

 

 

 

 

 

 

 

 

 

 

 

 

Zx,

+

£хгх,х,

 

 

 

 

218

_

(10 000+

ЮР) (3333,3 + 181,1)

_

10 100-3514,4 _

.9 R

„„„

~

100+181,1

 

_

 

281,1

 

 

f

= f

 

-

( ^ • + ^ , ) ( K , v a + Z Q ( F ^ + Z v a ) __

 

м. оЗ«

м.оЗо

 

A - =

.

= -

- = - =

 

~

 

 

 

 

 

^*o^*e -J- /Lx^i. xb -f- бх^х,

 

 

 

10 100»

34 300 000 ч.

 

 

 

 

 

3-1002

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким

образом,

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

отказов

звеньев

исследуемой

схемы ненадежности

будут равны

 

 

 

 

 

 

 

 

 

2 1 6 000 1 - 58, 4

= 0.0463.10-* 1/ч;

 

 

 

,26300-64,5

~

0.079810~4 1/4;

 

 

^ =

**• - 34300000-33,3

=

0,00029-Ю-4

1,ч.

По формуле (7.106)

получаем

 

 

 

 

 

 

 

In RC.H (1050) =

—1050 (1 +

2-0,0463 + 0,0798 +

2-0,00029) х

 

 

 

X Ю-1

= —0,1336,

 

 

 

откуда

 

 

Rc.u

(1050) -

0,875.

 

 

 

 

 

 

 

 

 

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

довательных цепей и всех звеньев (Y3., Z3.y

Эти характеристики

дают наглядное представление о тех звеньях

схемы ненадежности,

которые в первую очередь и определяют вероятность ее отказа. Та­ кими звеньями являются, во-первых, З х , во-вторых, 33 , и, в-третьих, 3 2 и 34 . Звенья 3 5 и 3 6 практически не снижают надежности системы, и их можно было бы совсем не учитывать. Надежность данной си­ стемы в основном определяет элемент х7 (звено 3\), который не имеет резервирования.

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

219

ю

ю

о

Рис. 64. Схема ненадежности системы, изображенной на рис. 43.

Вероятность безотказной работы резервированной части системы достаточно высока:

Д с ( 1 0 5 0 ) = fx*\\oS)

 

= 0 ' 9 7 2 -

 

Без учета восстановления она равна 0,939.

 

Ввиду того что на практике Z3

Y3,

формулу

(7.102) можно

еще несколько упростить, приняв

Y3 ^

Тм, 0 , и,

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

m

 

 

l3~—i=±

.

 

(7.109)

П (1

+ Р э / )

 

 

Для вычисления коэффициента готовности системы определим сначала по формуле (7.103) коэффициенты неисправности всех звеньев:

Р 3 ' =

100

n

 

m

 

 

 

1оооб =

0 ' 0 1 :

 

 

 

Рз, = Рз< = 2 -foT2 =

0,000271;

p^ = r2W5 =

0

- 0

0 0

5

1 6

;

Рз. = Рэ. =

й | ж

=

0

' 9

6

- 1

0 " в -

По формуле (7.107)

определим

Кг:

 

 

 

 

 

1

 

 

1

г — 1 +(10000 + 2-271 +516

+ 2-0,96)-10-°

1+0,011058

 

1

 

1 —0,011058 =

0,988942.

,011058

 

 

 

 

Итак, заметим, что

Rc, н (1050) = 0,875, а

Кг

= 0,989.

Завершая данный

 

параграф,

рассмотрим

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

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

 

системы с разнотипными

элементами. Найдем

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

V - Т

_ 7

-

J Z k + Z i ) i E a ± ! s )

'

_

'

1 м.о

'-•

 

Z j +

Z a

4 - + 4-

 

 

 

 

 

 

YJt

+ YiZi+YiZt

+

ZJ, _ _ZtZL

=

Y,Y2 + К А +

YaZt

 

Z\ + Z 2

 

Z j + Z 2

Zx + Z 2

 

 

 

 

 

' + P l +

p 3 4 .

 

(7.110)

 

 

 

 

PlP2 ( V i + v2)

 

 

221

По внешнему виду эта формула существенно отличается от вы­ ражения (7.60), которое запишем здесь с иной группировкой членов

1 + Pi +

Ра

+ ^ ( 1

+ P i ) + Ь - (!4 - р я ) +PiPa

 

7\ = -

 

 

Vet

 

 

(7.111)

 

 

PlPa ( V i +

V 2 )

4- Pip, (kx + Л.)

 

 

 

 

Определим

разницу между

7\

и Y для случая идентичных эле­

ментов. Формула (7.110) при этом

примет вид

 

 

 

 

Y

1 + 2 р

(7.112)

 

 

 

 

2\р

 

 

 

 

 

 

Учитывая

(7.62)

п

(7.112),

получим

 

T1

— Y — 1 +

1 + 2р

 

2рк

2рЛ.

1

Т

(7.113)

Ж ~ Т

0

 

Коэффициент готовности дублированной системы из идентичных элементов равен

1 4- 2р , (1 + Р)а _

(7.114)

 

2рХ

2рХ

(1 + Р )

 

2 '

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

§ 27. Основные понятия о логико-статистическом методе расчета надежности СЭС с помощью ЭЦВМ

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

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

Наиболее перспективным аппаратом, позволяющим стандартизи­ ровать решение широкого класса задач надежности, является метод статистического моделирования (метод Монте-Карло). Этот метод известен давно и начал свою вторую жизнь благодаря появлению быстродействующих ЭЦВМ.

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

222

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