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

книги / Надежность и диагностика компонентов инфокоммуникационных и информационно-управляющих систем.-1

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

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

На третьем этапе решаются уравнения системы и находятся искомые вероятности пребывания объекта в состояниях его работоспособности. Очевидно, что в рассматриваемом примере коэффициент готовности равен вероятности застать ответ в одном из трех работоспособных состояний:

Kг = P1(t) + P2 (t) + P3 (t).

(2.106)

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

ZP

(Z ) P (0)

= −P

(Z )λ

P

(Z )λ

+ P

(Z )λ

21

+ P

(Z )λ

31

,

 

1

 

1

 

1

 

12

1

 

 

13

2

 

 

3

 

 

 

 

 

ZP2 (Z ) P2 (0) = P1

(Z )λ12 P2

(Z )λ21

P2 (Z )λ24

+ P4 (Z )λ42

,

(2.107)

ZP

 

(Z ) P

(0)

= P

(Z )λ − P

(Z )λ

 

 

P (Z )λ

 

 

+ P (Z )λ

 

 

,

 

 

31

34

43

 

 

3

3

 

1

13

3

 

 

 

3

 

4

 

 

 

 

ZP

(Z ) P

(0)

= P

(Z )λ

24

+ P

 

(Z )λ

34

P

(Z )λ

42

+ P

(Z )λ

 

 

.

 

4

 

4

 

2

 

 

3

 

 

 

4

 

 

4

 

 

43

 

 

Из системы алгебраических уравнений находятся вероятности пребывания системы в состояниях Pi (Z ) . С помощью обратных пре-

образователей Лапласа полученные вероятности Pi (Z ) приводят

кискомому виду Pi (Z ) .

Вслучаях когда вероятности состояний являются постоянными, что характерно для установившегося режима работы, достигаемого в практике сравнительно быстро при существующих соотношениях ω(t) и λ(t) , система уравнений (2.105) становится системой ал-

гебраических уравнений, так как в этом случае dPi (t)dt = 0.

91

0 = −P1(t)λ12 P1(t)λ13 + P2 (t)λ21 + P3 (t)λ31,

 

0 = P

(t)λ − P (t)λ

21

P

(t)λ

24

+ P

(t)λ

42

,

 

 

1

12

2

2

 

4

 

 

 

 

 

0 = P1(t)λ13 P3 (t)λ31 P3 (t)λ34 + P4 (t)λ43,

 

(2.108)

 

0 = P

(t)λ

24

+ P (t)λ

34

P

(t)λ

42

+ P

(t)λ

43

,

 

2

 

3

4

 

 

4

 

 

 

 

 

 

P (t)

+ P (t) + P (t) + P (t) =1.

 

 

 

 

 

 

 

 

1

 

2

 

3

4

 

 

 

 

 

 

 

 

 

 

Добавление

последнего

уравнения

является

 

 

обязательным

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

Решение системы (2.108) позволяет определить установившийся коэффициент готовности. Рассмотрим порядок расчета на более конкретном примере.

Пример 2.6. Пусть необходимо определить надежность изделия, не имеющего резервирования, с заданными интенсивностями переходов – параметров потока отказов

ωω = const и интенсивностью восстанов-

1

µ

2

ления µ. Работоспособность системы

 

 

 

 

 

описывается графом (рис. 2.23): состоя-

 

Рис. 2.23

 

ние 1 – состояние работоспособности,

 

 

 

состояние 2 – состояние отказа.

Описание графа по приведенному правилу дает следующую систему уравнений вида (2.107):

ZP1(Z ) P1(0) = −P1(Z )ω+ P2 (Z )µ,

ZP2 (Z ) P2 (0) = P1(Z )ω− P2 (Z )µ.

Учитывая, что в момент включения t = 0 система должна быть исправна (P1(0) =1, P2 (0) = 0), получаем:

ZP1(Z ) + P1(Z )ω− P2 (Z )µ =1,

ZP2 (Z ) + P2 (Z )µ− P1(Z )ω= 0.

92

Отсюда

P2 (Z ) = PZ1(+t)µω,

(2.109)

= Z

P1(Z ) Z (Z +µ+ω) .

Обратное преобразование вероятности P1(Z ) требует приведения ее к табличному виду. Для этого умножим и разделим P1(Z ) на

(ω + µ):

P (Z ) =

Z

 

 

ω+µ =

µ(µ+ω+ Z ) + Zω

,

 

 

 

 

 

 

 

1

 

Z (µ+ω+ Z ) ω+µ Z (Z +ω+µ)(ω+µ)

 

 

 

P (Z ) =

1

 

µ

+

 

1

 

ω

.

 

 

 

 

 

 

 

1

Z

 

ω+µ

Z +ω+µ ω+µ

 

 

 

 

Отсюда, учитывая, что 1/Z соответствует 1(t), а 1(Z +ω+µ) соответствует e(µ+ω)t , получаем:

P (t) =

µ

+

ω

e(µ+ω)t .

(2.110)

 

 

1

ω+µ

 

ω+µ

 

 

 

 

Анализом полученного выражения устанавливаем, что P1(t) при t не может быть ниже величины µ(ω+µ) . Эта постоянная часть и является стационарным коэффициентом готовности изделия:

 

µ

=

 

1/Tв

=

 

Т0

 

= Kг.

(2.111)

 

ω+µ

1/Тв +1/Т0

Т0

+

 

 

 

 

Тв

 

Постоянная времени экспоненты

Тпэ =1 (ω+µ) .

Переходный

процесс длится 3÷4 Тпэ, после чего наступает установившийся режим.

Пример 2.7. Пусть ω = 10–2 1/ч, а µ = 1 1/ч. Тогда

Тпэ = 1+101 2 1ч.

93

Следовательно, переходный процесс длится 3 4 ч, а далее надежность системы определяется стационарным коэффициентом готовности:

Kг =

1

= 0,99.

1+102

На рис. 2.24 приведен график зависимости коэффициента готовности от времени.

K

1

 

 

ω

Тпэ

 

µ + ω

 

 

 

t

Рис. 2.24. График зависимости коэффициента готовности от времени

2.3.3. Определение среднего времени наработки на отказ системы с восстановлением

Время наработки на отказ (подразд. 2.1.1) определяется как

T = P(t)dt .

0

В то же время преобразование Лапласа определяется следующей формулой [2]:

 

P(Z ) = P(t)eZtdt,

(2.112)

0

 

т.е. при Z = 0 Т = P(Z).

94

Рассмотрим систему из двух восстанавливаемых блоков, один из которых основной, а другой – резервный. Перепишем систему

уравнений (2.107), заменяя

Pi (Z )

на Ti , с учетом того, что состоя-

ние 4 – состояние отказа. В результате P4 (Z ) =T4 = 0 , а также исчеза-

ет строка, соответствующая dP4 (t) / dt ,

 

 

 

1 = −T1λ12 T1λ13 +T2λ21 +T3λ31,

 

 

=T1λ12

T2λ21

T2λ24 ,

(2.113)

0

0

=T λ −T

λ

31

T λ

34

.

 

 

1

13

3

 

3

 

 

Среднее время наработки на отказ всей системы T =T1 +T2 +T3 , так как 1, 2, 3 – состояния работоспособности.

Пример 2.8. Пусть

λ = λ = λ

24

= λ

34

=102

1/ ч,

12

13

 

 

 

λ21 = λ31 = λ42 = λ43 =1 1/ ч.

Тогда

 

 

 

 

 

 

 

1 = −T 2 102 +T +T ,

 

1

 

 

 

2

3

 

0 =T 102

T (1+102 ),

 

1

 

 

2

 

 

 

0 =T 102

T (1+102 ).

 

1

 

 

3

 

 

Решив систему, получаем:

 

 

 

 

T1 5050,

T2 50, T3 50.

Среднее время наработки на отказ

T=T1 +T2 +T3 = 5150 ч.

2.3.4.Расчет надежности систем с восстановлением при основном (последовательном) и параллельном

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

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

95

что отказ любого из них приводит к отказу всей системы (рис. 2.25). Для простоты предположим, что каждый образец имеет одинаковую интенсивность отказов ω и интенсивность ремонтов µ.

Предположим, что у нас имеется один ремонтник. Составим граф переходов системы (рис. 2.26).

 

 

 

 

 

 

2

 

ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b1

 

b2

 

 

0

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

µ

 

2µ

Рис. 2.25. Структура системы

Рис. 2.26. Граф работоспособности

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

 

системы

 

Обозначим:

0 – состояние системы, в котором оба образца исправны; 1 – состояние системы, когда один образец исправен, а другой

ремонтируется; 2 – состояние системы, когда оба образца неисправны, один

ремонтируется.

Из состояния 0 система может перейти в состояние 1 с интенсивностью отказов 2ω. Из состояния 1 система может перейти в состояние 0 с интенсивностью восстановления µ и в состояние 2 с интенсивностью отказов ω. Из состояния 2 система может перейти в состояние 1 с интенсивностью восстановления µ.

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

уравнений:

dP0dP0dP2

/ dt = −2ωP0 (t) P1(t),

 

/ dt = 2ωP0 (t) (ω+µ)P1(t) P2 (t),

(2.114)

/ dt = ωP1(t) −µP2 (t).

 

Будем искать решение только для установившегося значения. Тогда система дифференциальных уравнений перейдет в систему линейных уравнений:

96

0 = −2ωP0 P1,

 

 

 

 

= 2ωP0 (ω+µ)P1

P2

,

 

0

(2.115)

0

= ωP

−µP ,

 

 

 

1

 

2

 

 

 

1

= P + P

+ P .

 

 

 

 

0

1

2

 

 

 

Отсюда коэффициент готовности

Kг = Р0 =

 

 

µ2

.

(2.116)

µ

2

2

 

 

+2ωµ+ 2ω

 

Пример 2.9. Пусть ω = 102 1/ч, µ = 1 1/ч. Определить коэффициент готовности:

Kг = Р0 =

1

= 0,98.

 

1+ 2 102 + 2 104

Вобщем случае, если у нас имеется n образцов оборудования

иодин ремонтник, справедлива формула:

 

µ n

 

 

 

P0 =

 

ω

 

.

(2.117)

n

µ j

 

 

n!j=0 ω

j!

 

В качестве другого крайнего случая рассмотрим систему, когда количество ремонтников равно количеству образцов оборудования. Пусть на оба образца имеется два ремонтника. Составим граф перехода системы (рис. 2.27):

0– состояние системы, в котором оба образца исправны;

1– состояние системы, когда один образец исправен, а другой ремонтируется;

2– оба образца неисправны и ремонтируются.

97

2ω ω

0

1

2

µ 2µ

Рис. 2.27. Граф переходов системы

Запишем по графу переходов систему уравнений для установившегося значения:

0 = −2ωP0 P1,

0 = 2ωP0 (ω+µ)P1 + 2µP2 ,0 = ωP1 2µP2 ,

1 = P0 + P1 + P2.

Решая систему, получим:

K

г

= Р

=

µ2

 

=

 

 

 

µ2

 

.

(2.118)

 

2

 

2

 

 

 

2

 

0

 

(µ+ω)

 

µ

+

2ωµ+ω

 

 

 

 

 

 

 

 

 

 

 

 

Пример 2.10. Пусть ω = 10–2 1/ч, µ = 1 1/ч. Определить коэф-

фициент готовности:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Kг =

 

1

 

 

 

 

 

= 0,98.

 

 

 

 

 

 

 

 

 

 

1+ 2 102 + 2 104

 

 

 

 

Вобщем случае, если у нас имеется n образцов оборудования

иn ремонтников,

K

г

= Р

=

µn

,

(2.119)

(µ+ω)n

 

0

 

 

 

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

98

дать, так как для каждого образца имеется свой ремонтник и Kг каждого образца не зависит от Kг остальных.

Рассмотрим систему из двух образцов оборудования, соединенных параллельно (рис. 2.28). Как уже указывалось, в этом случае отказ системы наступает только при отказе всех элементов системы.

Предположим, что у нас имеется один ремонтник, который сразу начинает ремонтировать отказавший элемент.

Составим граф переходов системы (рис. 2.29):

0– состояние системы, в котором оба образца исправны;

1– состояние системы, когда один образец исправен, а другой ремонтируется;

2– оба образца неисправны, один ремонтируется.

 

 

 

 

 

 

 

 

2ω

ω

 

 

 

1

 

 

 

 

 

 

 

 

 

 

0

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

µ

2µ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 2.28. Структура системы

Рис. 2.29. Граф переходов системы

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

 

 

 

Запишем по графу переходов систему уравнений для установившегося значения:

0 = −2ωP0 P1,

0 = 2ωP0 (ω+µ)P1 P2 ,0 = ωP1 2µP2 ,

1 = P1 + P2 + P3.

Решив систему, получим:

Р0 =

µ2

,

Р1 =

2ωµ2

 

 

.

µ2 + 2ωµ+2ω2

µ2 +2ωµ+2ω2

99

Коэффициент готовности

K

г

= Р + P =

 

 

µ2 +

2ωµ

 

.

(2.120)

 

2

 

 

2

 

0 1

µ

+2ωµ+2ω

 

 

 

 

 

 

 

 

 

 

Пример 2.11. Пусть ω = 10–2 1/ч, µ = 1 1/ч. Определить коэф-

фициент готовности:

 

 

 

 

 

 

 

 

 

 

Kг =

 

1+ 2 102

= 0,9998.

 

 

1+ 2 102 + 2 104

 

Имеется довольно многочисленный класс систем, в которых обслуживание невозможно начать до наступления полного отказа системы. Это может произойти, если контролируется только выход из строя всей системы, а не отдельных образцов оборудования. Допустим, у нас имеется 2 образца оборудования, соединенных параллельно. После того как откажет вся система, два ремонтника начинают ремонтировать каждый свой элемент.

Составим граф переходов системы (рис. 2.30):

0 – состояние системы, в котором оба образца исправны; 1 – состояние системы, когда один образец исправен, а другой

неисправен, но не ремонтируется; 2 – состояние системы, когда оба образца неисправны и ремон-

тируются.

2ω ω

0

1

2

2µ

Рис. 2.30. Граф переходов системы

Запишем по графу переходов систему уравнений для установившегося значения:

100