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

2685

.pdf
Скачиваний:
0
Добавлен:
15.11.2022
Размер:
2.06 Mб
Скачать

Рис. 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,

функцией

работоспособности системы (булевой функцией) переменных х12,...,хn

Описанное( , ,… ) =

 

если система функционирут

,

(4.29)

0,

 

.

 

1,

 

еслисистема отказывает

при

выше

 

поведение

систем

 

функционировании или отказе всех ее элементов и свойство монотонности могут быть записаны в виде следующих

соотношений:

 

 

 

 

S(1, ...,1)=1

(

,…,

),

(4.30)

S(0, ...,0)=0

если для( ,…, ) ≤

реализации

 

и

всех

компонент

выполняется условие

,

= 1… .

 

 

 

 

75

 

 

Выражения (4.28)-(4.30) характеризуют модель полностью и однозначно. Однако может оказаться, что в действительности функция работоспособности системы вовсе не зависит от некоторых из n булевых переменных. Подобные элементы для поведения системы являются в таком случае несущественными. Точнее, определим их следующим образом: v-ый элемент системы называется несущественным, если

условие

 

 

+1,…. ) = ( 1, 2,…, −1,1, +1,…. )

, справедливо( 1, 2,…, −1,0,

, ,…,

для

всех

реализаций булевых переменных

,0,

,….

. В противном случае v-ый элемент

 

 

 

называется существенным.

Систему называют нередуцируемой, если она не содержит несущественных элементов. В противном случае система будет допускать редуцирование [52].

Для краткости будем использовать векторную форму

записи. Пусть

 

есть булев вектор-столбец с n

элементами.

Индекс T обозначает операцию транспонирования.

 

= ( ,…, )

 

ТТ

Вчастности, 1=(1, ...,1) и 0=(0,...,0) . Используемые в дальнейшем неравенства между булевыми векторами

определяются следующим образом:

≥ ,если ≥ для = 1,2,…, , > ,если

(4.31)

 

 

 

 

 

 

 

 

 

,

иразличаются по крайней мере одной компонентой

вектора

 

и

. Тогда формулу (4.30)

можно записать в

 

 

 

 

 

 

 

 

 

 

 

 

следующем

 

виде:

 

,если ≤ ,

(4.32)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В качестве вычислительных операций над булевыми переменными ( ,…, ) применяют булево сложение (дизъюнкцию), обозначаемое знаком " ", и булево умножение (конъюнкцию), обозначаемое знаком " ". В таблице 6 приведены результаты сложения и умножения для всех 2n реализаций булева вектора.

Если обе возможные реализации каждой булевой переменной ограничиваются лишь числами 0 и 1, то булево сложение и умножение можно заменить эквивалентными

76

алгебраическими операциями.

Определение операции сложения и умножения в булевой алгебре

Номер

Булевы переменные

 

 

 

 

 

реализации

 

 

 

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

0

0

...

0

0

0

 

0

0

1

 

0

0

0

...

0

0

1

 

0

1

2

 

0

0

0

...

0

1

0

 

0

1

3

 

0

0

0

...

0

1

1

 

0

1

.

 

.

.

.

...

.

.

.

 

.

.

.

 

.

.

.

...

.

.

.

 

.

.

.

 

.

.

.

...

.

.

.

 

.

.

2n-3

 

1

1

1

...

1

0

1

 

0

1

2n-2

 

1

1

1

...

1

1

0

 

0

1

2n-1

 

1

1

1

...

1

1

1

 

1

1

 

С помощью таблицы легко убедиться в справедливости

следующих

 

= min

,

,…

,

 

 

соотношений:

 

= ∏

 

 

= 1−

 

(1−

) = max

,

,…, .

 

 

 

 

 

Часто оказывается целесообразным наряду с булевой переменной х рассматривают и ее отрицаниепеременную ̅, которая оказывается определенной в той же области значений, что и x; при этом в силу определения ̅реализации х и ̅ никогда не совпадают. В нашем случае ̅= 1− , поскольку реализация булевых переменных есть числа 0 и 1 [52].

4.4.2. Булевы функции работоспособности системы и схемы расчета надежности

Для некоторых классов систем структурную функцию работоспособности S(x)=S( , …, ) можно вычислить

77

достаточно просто. Мы будем наряду с функцией работоспособности системы рассматривать также и

структурную функцию отказа системы

(x),определяемую как

x

 

x

x),

(4.33)

(x)=1-S( ) или

( )=1-

̅(

 

Оба определения по

формуле являются эквивалентными.

 

̅

̅

 

Существуют 2 класса систем.

 

 

 

Для первой системы имеем:

 

 

 

-работоспособна только тогда ,когда все n

ее элементов

работоспособны. Функция работоспособности имеет вид в

таком случае:

 

 

 

, (4.34)

 

-̅

=1-

 

(x)=min( , …, )=

 

 

=

 

 

система отказывает, когда хотя бы один из ее n

элементов отказывает. Для функции отказа отсюда следует:

 

 

=1-

 

 

= 1 −∏

, (4.35)

)=max( , …, )=

 

 

 

̅(x

Графически подобную можно представить как систему с

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

х элементов или

 

 

параллельнымсоединением

хэлементов соответственно.

Рис. 4.4. Последовательное и параллельное соединение элементов

Буквы А и Е илиА и Е символически обозначают вход и выход системы [52].

Для второй системы имеем:

-система работоспособна, когда работоспособен хотя бы один из ее элементов. Для функции работоспособности отсюда

следует:

 

 

 

 

̅(x)=max( , …, )=

=

1 −∏

,

(4.36)

 

78

 

 

 

 

-система отказывает только тогда, когда отказывают все n

элементов. Для функции отказа отсюда следует:

 

 

 

̅

)=1-

=

 

,

(4.37)

 

 

(x)=min( , …,

 

 

 

Графически подобную можно представить как систему с

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

 

 

 

соединением

х

 

.

 

 

 

 

 

хэлементов или параллельным соединением

 

элементов соответственно

Рис. 4.5. Последовательное и параллельное соединение элементов

4.4.3. Свойства и представления булевых функций работоспособности систем.

Реализацию х' булева вектора х будем называть путем но работоспособности системы, если S (х') = 1. Реализацию х' булева вектора х будем называть сечением по работоспособности системы, если S(х')=0. Соответственно этому реализация х' булева вектора х называется сечением по

̅'

отказу системы, если (х) = 0. Реализация х' булева вектора х

̅'

называется путем по отказу системы, если (х) =1 [52]. Особое значение имеют так называемые минимальные

пути и минимальные сечения. Реализация х' булева вектора х называется минимальным путем по работоспособности системы, если S(х') =1, но для любой реализации х'', такой, что х'' < х', функция S(х") = 0. Реализация х' булева вектора х называется минимальным сечением по работоспособности системы, если S(х') =0, однако для любой реализации х", такой,

79

что х">х', функция (х') *= 1. Реализация х′ называется

минимальным сечением по отказу системы, если (х′) = 0,

однако для любой реализации х", такой, что х" > х′, функция

(х") = 1. Реализация х′ булева вектора х называется

минимальным путем по отказу системы, если (х′) =1, однако

для любой реализации х" такой, что х" < х" функция S(х") =0. Можно непосредственно проверить, что всякий

минимальный путь (минимальное сечение) по работоспособности оказывается дополнением минимального сечения (минимального пути) по отказу и наоборот. Обозначим общее число минимальных путей и сечений по работоспособности через m и s, а общее число минимальных путей и сечений по отказу — через и ̅. Тогда =s, ̅= m. Знание всех минимальных путей и минимальных сечений позволяет непосредственно получить аналитические выражения для функции работоспособности и отказа системы. Под эквивалентной схемой минимального пути

(эквивалентной схемой дополнения минимального сечения x ) будем понимать последовательное (параллельное) соединение тех компонентов вектора х, которые представлены

вединицами (в x нулями) [52].

Теорема 1а [52]. Эквивалентную схему работоспособности системы можно представить в виде параллельного соединения эквивалентных схем всех

минимальных путей по

работоспособности,

при

этом, если

 

 

минимальные пути по работоспособности и

x (σ = 1,…,s̅= m)

 

 

 

 

 

(

= 1,…,

)− -минимальные сечение по отказу, функция

работоспособности системы имеет вид:

x )

 

 

 

S x =

x = 1 − (1 −

],

(4.38)

 

 

= 1 −

[1 −

(1− x )

80

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]