Рис. 4.3. Общая живучесть для d=6 в двух вариантах
Условием, что рассматриваемые числа реализуются в виде графа, является отношение (4.11). Следствие из него – формула (4.12).
∑ d(i) = 2m, |
(4.11) |
||
< d >= |
|
, |
(4.12) |
|
|||
Используя дисперсионные соотношения для случайной величины d(i):
S |
|
∑ (d(i)−< >), |
(4.13) |
|
Получают верхнюю границу вероятности связного графа:
P(p = 1) ≤ 1 −nq [1− |
|
q ( |
|
)], |
(4.14) |
|
|
||||
|
Нет необходимости искать возможные значения Р, т.к. на практике можно ограничиться вполне приемлемым значением Р(р =1) ≥3/4. Для таких значений Р множитель в квадратных можно заменить значением 0,5. Тогда критерий живучести стохастического однородного графа с заданными n, d, q, d<<n:
71
P(p = 1) ≤ 1 − |
|
q , |
(4.15) |
|
Если задано Р( р= 1), то:
d ≥ |
|
ln(1 − ) − |
|
, |
(4.16) |
|
|
т.е. необходимое условие для степени вершин, чтобы граф был связным, т.е. из каждой вершины в среднем должно исходить d ребер при заданных qи n.
Исследуем изменения структуры СИС при внезапном воздействии ЧС.
Введем условия:
1.Каждая вершина имеет в среднем d исходящих
ребер;
2.Каждое ребро, инцидентное вершине Vi также инцидентноVjсвероятностью:
P = |
|
,ij=1…n, |
(4.17) |
|
3.Всевершиныиребраидентичны;
4.Воздействие ЧС направлено в некоторую случайно выбраннуюобластьw,тогдавероятностьпораженияравна—w/D,где D-площадь,занимаемаяСИС;
5.Частота возникновения последствий воздействия ЧС составляетm(t)единицнаединицуплощади;
6.ВоздействиеЧСнаСИСодномоментное;
7.Вершина разрушается в том случае, если подвергается воздействиюЧС,мощностьюнеменееKb S ;
8.Ребро разрушается в том случае, если подвергается воздействиюЧС,мощностьюнеменееKPS ;
9.ВосстановительныеработынаСИСнепроводятся. Если m(t)велико, то ожидаемое число разрушенных вершин
(или ребер) будет соответствовать вероятности того, что любая даннаявершинаможетбытьразрушена.
72
Пусть gkb (m)и gkP(m)обозначают ожидаемые части вершин и ребер, которые подверглись воздействию ЧС с частотой
возникновенияпоследствийm. |
|
||
Вероятность того, |
что |
некоторая вершина ks |
не будет |
разрушенаравна: |
g |
(m), |
|
∑ |
(4.18) |
||
Вероятность того, что некоторое ребро ki не будет разрушено, равна:
∑ |
g (m), |
(4.19) |
Пусть πbи πP- вероятности того, что воздействиеЧСразрушает данные вершину и ребро, тогда, согласно распределению Пуассона длявероятностей:
g (m) = e π |
(π ! |
) |
, |
(4.20) |
g (m) = e π |
(π ! |
) |
, |
(4.21) |
Предположим, что после удара осталось в среднем а неповрежденных ребер, направленных из Vi, причем Vj считается не разрушенной. Для графов, удовлетворяющих условию (4.4) существуетсоотношение:
b = 1 −e |
,a.= |
( ) |
, |
(4.22) |
|
|
(4.23) |
||
P{Q } = q |
|
|
|
|
Следовательно,учитываявыражения(4.18)–(4.23): |
||||
b = 1− exp{−d ∑ |
g (m) [∑ |
g (m)]b}, (4.24) |
||
|
73 |
|
|
|
Таким образом, выражение (4.24) определяет среднее количество вершин, не разрушенных воздействием ЧС,
аb∑ g (m)–среднее число вершин первоначального графа, которые могут установить связь с другими после воздействия ЧС и выражение (4.24) можно использовать для определения коэффициентаживучести.
Компонентаa ,учитывая(4.18),(4.19),(4.20)будетиметьвид:
d ∑ |
g |
(m) |
∑ |
|
g |
(m) |
= |
( |
) |
, |
|
(4.25) |
||||
Задавb>b0 |
каккоэффициентживучести: |
|
|
|
|
|
|
|||||||||
d ∑ |
g |
(m) |
∑ |
|
g |
(m) |
≥ |
( |
) |
, |
(4.26) |
|||||
Учитываявыражения(4.20)–(4.22)получим: |
|
|
|
|
|
|||||||||||
d ∑ |
(π |
! |
) |
∑ |
(π |
! |
) |
≥ (π π ) |
( |
( |
|
) |
), (4.27) |
|||
Принимаявовнимание,чтоks иkL—целые числа,ограничение некоторымпредельнымзначением,находимминимальноезначениеd
[20].
4.4.Булевы модели оценки живучести
4.4.1.Сущность булевых моделей
Предполагается, что исследуемую систему можно достаточно хорошо описать следующим образом. Система состоит из конечного числаэлементов(подсистем,элементарныхстандартныхблоков).Для каждого элемента допускаются лишь два возможных состояния — полной работоспособности и полного отказа. Точно так же предполагается, что и система может находиться лишь в двух состояниях — полной работоспособности и полного отказа. Всякая возможность частичного функционирования всей системы или ее
74
элементовисключается[52].
Работоспособность или отказ системы определяются состояниями ее элементов. В предельном случае одновременного отказа или одновременного функционирования всех элементов система должна соответственно отказывать или работать. Предполагается
рует),тосистемадолжнафункционировать
выполнение следующего свойства монотонности: если система функционирует, когда отказало некоторое подмножество М1, ее элементов (а дополнительное множество М элементов функционитакжеивтомслучае,если
1
отказалолишьподмножествоМ2 М1 элементов.
Названная модель функционирования системы явно приспособлена к применению аппарата булевой алгебры. Состояния системы и каждого элемента можно описать с помощью булевых переменных, которые принимают значение 1 (в случае работоспособности)и0(вслучаеотказа).
Пусть система состоит из n элементов, при этом v-муэлементу
соответствует булева |
переменная |
возможными |
реализациями |
|||
|
|
1или0,т.е. |
|
|
|
|
которойявляютсячислаесли система функционирут, |
|
(4.28) |
||||
здесь v=1, = |
1, |
если система отказывает |
, |
|||
2, ..., n. |
определяется структурной. |
|
||||
Состояние системы0, |
функцией |
|||||
работоспособности системы (булевой функцией) переменных х1,х2,...,хn
Описанное( , ,… ) = |
|
если система функционирут |
, |
(4.29) |
||
0, |
|
. |
||||
|
1, |
|
еслисистема отказывает |
при |
||
выше |
|
поведение |
систем |
|
||
функционировании или отказе всех ее элементов и свойство монотонности могут быть записаны в виде следующих
соотношений: |
|
|
|
|
S(1, ...,1)=1 |
( |
,…, |
), |
(4.30) |
S(0, ...,0)=0 |
||||
если для( ,…, ) ≤ |
||||
реализации |
|
и |
всех |
компонент |
выполняется условие |
, |
= 1… . |
|
|
≤ |
|
|||
|
|
75 |
|
|