книги из ГПНТБ / Рябинин И.А. Надежность судовых электроэнергетических систем и судового электрооборудования учебник
.pdfМатематическое ожидание длительности импульса потока совпа дений, образованного в результате перекрытия во времени 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 + Р)а _ |
2р |
(7.114) |
|
|
2рХ |
2рХ |
(1 + Р ) |
|
|
2 ' |
что в точности совпадает с формулой, приведенной в табл. 32, для случая нагруженного резерва с неограниченным восстановлением. Все это еще раз подтверждает правильность полученных соотношений.
§ 27. Основные понятия о логико-статистическом методе расчета надежности СЭС с помощью ЭЦВМ
Построение корректной математической модели изучаемых про цессов является основным условием, гарантирующим получение объективных результатов. При этом важно уметь согласовывать способ описания с целью исследования, различать главные и второ степенные факторы.
Появление и развитие электронных цифровых вычислительных машин привело к резкому увеличению возможностей применения математических методов для решения различных теоретических и практических задач. Существенно изменились понятия и об инженер ных методах расчета.
Наиболее перспективным аппаратом, позволяющим стандартизи ровать решение широкого класса задач надежности, является метод статистического моделирования (метод Монте-Карло). Этот метод известен давно и начал свою вторую жизнь благодаря появлению быстродействующих ЭЦВМ.
Сущность статистического моделирования состоит в том, что исследуемый процесс функционирования судовой электроэнергети ческой системы представляется математической вероятностной мо делью, изоморфной во всех существенных чертах рассчитываемому процессу. Данная модель многократно испытывается, в результате чего определяются требуемые статистические характеристики. При
222