Управление рисками. Практикум. Учебное пособие
.pdfследующих основных факторов: вид и кратность резервирования, |
|
|||||||
дисциплина обслуживания. |
|
|
|
|
|
|||
|
Повышение |
надежности |
систем |
достигается |
за |
счет |
||
применения двух видов резервирования: с постоянно включенным |
|
|||||||
резервом и по методу замещения. При |
этом |
обслуживание |
|
|||||
системы может осуществляться с двумя видами приоритета- |
|
|||||||
прямым и обратным. При прямом |
приоритете |
система |
||||||
обслуживается в порядке ее поступления в ремонт, при обратном |
|
|||||||
- |
первой |
обслуживается |
система, поступившая |
в |
ремонт |
|
||
последней. |
Наиболее |
эффективным способом |
обеспечения |
и |
||||
повышения надежности и снижения техногенного риска является структурное резервирование с возможностью восстановления отказавших элементов в процессе функционирования системы. Однако применение резервирования повышает стоимость системы
иее эксплуатации. Поэтому кратность резервирования
ограничивается, |
и |
в |
большинстве |
случаев |
применяется |
||
резервирование |
с кратностьюm = |
1 (дублирование). Из двух |
|||||
указанных |
|
видов |
резервирования |
наибольший |
выигрыш |
||
надежности |
достигается |
при |
резервировании замещением. |
||||
Недостатки этого резервирования следующие:
1) физическая реализуемость требует автомат-контроля состояния системы и коммутации при отказе работающей системы;
2) снижается производительность системы, т. к. резервные системы до замещения не работают.
По этим причинам на практике наиболее часто применяется резервирование с постоянно включенным резервом.
Наработка |
на |
отказ |
и |
коэффициент |
готовнос |
|
резервированных |
восстанавливаемых |
систем |
при |
одной |
||
обслуживающей бригаде вычисляются по следующим формулам: |
|
|||||
а) для системы с постоянно включенным резервом: |
|
|
||||
T = T å |
1 |
, |
K = |
T |
; |
|
m |
|
|
|
|
|
|
0 i=0 |
|
|
г |
|
|
|
(i +1)!ri |
|
1 |
|
|||
|
|
|
|
T + |
|
|
|
|
|
|
m |
|
|
61
б) для резервированной системы замещением:
T = T |
m |
1 |
, |
K |
|
= |
T |
|
, |
|||
å |
|
|
|
|
||||||||
|
|
|
1 |
|||||||||
|
|
0 i=0 ri |
|
|
г |
|
T + |
|
||||
|
|
|
|
|
|
|
|
|
m |
|
||
|
|
|
|
|
|
|
|
|
|
|
||
где r = |
l |
|
– |
интенсивность восстановления системы; |
||||||||
m |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|||
срок службы системы до первого отказа.
Показатели надежности Т и Кг зависят обслуживающих бригад. Приведем формулы обслуживающих бригад:
а) для системы с постоянно включенным резервом:
T = T |
æ m -1 |
2i |
|
+ |
2m -1 |
ö |
, Kг |
= |
|
T |
; |
|||
|
å |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
||||||
ç |
|
i |
|
m ÷ |
|
|
1 |
|||||||
0 |
|
(i +1)!r |
|
(m +1)!r |
|
|
|
|
|
|||||
|
è i =0 |
|
|
ø |
|
|
T + |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
2m |
|
||
б) для резервированной системы замещением:
T = T |
æ å |
2 |
i |
+ 2 |
m -1 |
ö , Kг = |
T . |
||||
|
m -1 |
|
|
|
|
|
|
|
|
|
|
|
ç |
|
|
|
|
|
|
÷ |
|
|
|
0 |
|
|
|
|
|
|
|
1 |
|
||
r |
i |
|
r |
m |
|
|
|||||
|
è i =0 |
|
|
|
ø |
T + |
|
||||
|
|
|
|
2m |
|||||||
|
|
|
|
|
|
|
|
|
|||
1
T0 = l –
от числа для двух
Для случая высоконадежных систем, когда r ≤ 0,001, дисциплина обслуживания не оказывает существенного влияния на надежность резервированных восстанавливаемых систем.
Риск системы определяется по формуле:
|
t |
R(t) = rM (t) = rlò pm (t )dt , |
|
|
0 |
где M(t) |
- среднее число отказов системы в течение времениt, |
pm(t) - |
вероятность пребывания системы в пред отказовом |
состоянии в момент времени t. Для расчетов можно использовать |
|
простую приближенную формулу |
|
R(t) = rl pmt , |
|
где pm- |
стационарная вероятность пребывания системы в |
предотказовом состоянии.
62
Порядок выполнения работы
Первое задание. Определить наработку на отказT и коэффициент готовности Кг резервированной восстанавливаемой автоматизированной системы.
Пример выполнения задания. Имеем следующие исходные данные:
-срок службы системы Tс = 5 лет;
-время непрерывной работы t = 1000 ч;
-интенсивность отказа системы λ = 10-5 ч-1;
- интенсивность восстановления системы r = l = 1; 0,1; m
0,05; 0,01;
-кратность резервирования m = 1;
-риск из-за отказа системы r = 150000 усл. ед.:
- |
допустимый |
риск в |
течение времени |
непрерывной |
работы R(1000)=360 усл. ед. |
|
|
||
На |
рис. 1 и 2 |
приведены |
структурные схемы |
и графы |
состояний системы соответственно при общем постоянном резервировании (а) и резервировании замещением (б).
Рис. 1. Структурные схемы резервированных систем
63
Рис. 2. Графы состояний резервированных систем
Расчетные формулы для случая дублированной системы (m=1) имеют вид:
а) дублированная система с постоянно включенным резервом:
одна обслуживающая бригада (n = 1):
æ |
|
1 |
ö |
|
|
|
1 + r |
|
; |
(1) |
T = T0 ç1 |
+ |
|
÷ |
, K г |
= |
|
|
|
||
r |
1 + r + r |
2 |
||||||||
è |
|
ø |
|
|
|
|
|
|||
две обслуживающие бригады (n = 2):
æ |
|
1 ö |
|
|
|
1 + 2r |
|
; |
|
T = T0 ç1 |
+ |
|
÷ |
, K г |
= |
|
|
|
|
|
1 + 2r + 2r |
2 |
|||||||
è |
|
2r ø |
|
|
|
|
|||
б) дублированная система замещением: одна обслуживающая бригада (n = 1):
æ |
|
1 |
ö |
|
|
|
1 + r |
|
; |
T = T0 ç1 |
+ |
|
÷ |
, K г |
= |
|
|
|
|
r |
1 + r + r |
2 |
|||||||
è |
|
ø |
|
|
|
|
|||
две обслуживающие бригады (n=2):
æ |
|
1 |
ö |
|
|
|
1 + r |
|
||
T = T0 ç1 |
+ |
|
÷ |
, Kг |
= |
|
|
|
|
; |
r |
|
|
1 |
|
||||||
è |
|
ø |
|
|
1+ r + |
r2 |
|
|||
|
|
|
|
|
|
|||||
|
|
|
|
|
|
2 |
|
|
||
(2)
(3)
(4)
64
в) нерезервированная система |
|
|||||||
T =T1 =T0 |
= |
1 |
, |
Kг |
= |
1 |
. |
(5) |
l |
|
|||||||
|
|
|
|
1+ r |
|
|||
Из приведенных формул видно, что наработка на отказ и коэффициент готовности резервированной системы являются функциями ρ. Это позволяет автоматизировать расчеты, используя математическую систему Mathcad.
Введем исходные данные:
Tc := 5
t := 1000
l := 10- 5
m := 1 |
i := 0 .. 3 |
r := 15×104
Создадим матрицу формул размерностью9×1 по числу формул (1)-(5) для Т и Кг:
65
é |
|
æ |
1 ö |
|
|
ù |
||||
ê |
|
t×ç1 + |
|
|
|
÷ |
|
|
ú |
|
|
|
|
|
|
||||||
ê |
|
è |
2×r ø |
|
ú |
|||||
ê |
|
1 + 2×r |
|
|
|
ú |
||||
ê |
|
|
|
|
|
|
|
|
|
ú |
|
|
|
2×(r) |
2 |
||||||
ê 1 + 2×r + |
|
|
ú |
|||||||
ê |
|
æ |
1 ö |
|
|
ú |
||||
ê |
|
t×ç1 + |
|
|
|
÷ |
|
|
ú |
|
|
|
|
|
|
||||||
ê |
|
è |
2×r ø |
|
ú |
|||||
|
1 + 2×r |
|
|
|
||||||
ê |
|
|
|
|
ú |
|||||
ê |
1 |
+ 2×r + |
|
2×(r) |
2 |
ú |
||||
ê |
|
|
ú |
|||||||
ê |
|
æ |
1 ö |
|
|
|
ú |
|||
T(r) := ê |
|
t×ç1 + |
|
|
÷ |
|
|
|
ú |
|
|
|
|
|
|
||||||
|
è |
|
r ø |
|
|
|
||||
ê |
|
1 + r |
|
|
|
ú |
||||
ê |
|
|
|
|
ú |
|||||
|
|
|
|
|
|
|
|
|
||
ê |
|
1 + r + |
|
( r) |
2 |
|
|
ú |
||
ê |
|
|
|
|
|
ú |
||||
|
æ |
|
1 |
ö |
|
|
|
|||
ê |
|
|
|
|
|
ú |
||||
ê |
|
t×ç1 + |
÷ |
|
|
|
ú |
|||
|
è |
|
r ø |
|
|
|
||||
ê1 + r ú
êú
ê1 + r + 1 ×(r)2 ú
ê |
|
2 |
ú |
ê |
1 |
|
ú |
êú
ë |
1 + r |
û |
Выполним расчеты для различных ρ
66
Расчеты выполняются для всех значенийr. Сведем полученные результаты в табл. 1.
Таблица 1
|
Результаты расчета наработки на отказ и коэффициента |
||||
|
|
|
готовности системы |
|
|
r |
|
1,0 |
0,1 |
0,05 |
0,01 |
TP1/T0 |
|
1,5 |
6,0 |
11,0 |
51 |
KP1 |
|
0,6 |
0,984 |
0,995 |
1,0 |
TP2/T0 |
|
1,5 |
6,0 |
11,0 |
51 |
KP2 |
|
0,6 |
0,984 |
0,995 |
1,0 |
TZ1/T0 |
|
2,0 |
11,0 |
21,0 |
101 |
KZ1 |
|
0,667 |
0,991 |
0,998 |
1,0 |
TZ2/T0 |
|
2,0 |
11,0 |
21,0 |
101 |
KZ2 |
|
0,8 |
0,995 |
0,999 |
1,0 |
Кг |
|
0,5 |
0,909 |
0,952 |
0,99 |
В табл. 1 приняты следующие обозначения:
- ТР1, ТР2 - наработка на отказ системы с постоянно включенным резервом с одной и двумя обслуживающими бригадами соответственно;
- TZ1, TZ2 - наработка на отказ системы, резервированной
по принципу замещения с одной и двумя обслуживающими бригадами соответственно:
67
- КР1, КР2 - коэффициент готовности системы с постоянно |
|
|
|||||||||||
включенным |
резервом |
с |
одной |
и |
|
двумя |
обслуживающими |
|
|||||
бригадами соответственно; |
|
|
|
|
|
|
|
|
|
||||
- |
KZ1, |
KZ2 |
- |
|
коэффициент |
|
готовности |
системы, |
|
||||
резервированной по принципу замещения с одной и двумя |
|
||||||||||||
обслуживающими бригадами соответственно. |
|
|
|
|
|
||||||||
Анализ |
данных |
табл. 1 позволяет |
|
сделать |
следующие |
|
|||||||
важные выводы: |
|
|
|
|
|
|
|
|
|
|
|
||
1) наработка |
на |
|
отказ |
резервированной |
системы |
с |
|||||||
кратностью m = 1 не зависит от числа ремонтных бригад; |
|
|
|
||||||||||
2) при |
малых |
|
значенияхr |
|
наработка |
|
на |
отказ |
|
||||
дублированной |
системы |
|
замещением |
практически |
вдвое |
|
|||||||
больше, |
чем |
при |
дублировании |
с |
постоянно |
включенным |
|
||||||
резервом; |
|
|
|
|
|
|
|
|
|
|
|
|
|
3) резервирование с восстановлением является средством |
|
||||||||||||
повышения наработки на отказ системы: например, в случае |
|
|
|||||||||||
резервирования |
замещением |
приr = |
0,01 наработка |
на отказ |
|
|
|||||||
TZ1 = TZ2 = Т0×101=1010000 ч, что составляет примерно 115 лет и |
|
|
|||||||||||
в 23 раза превышает срок службы системы (5 лет); |
|
|
|
|
|||||||||
4) число |
ремонтных |
|
бригад |
оказывает |
незначительное |
|
|||||||
влияние на коэффициент готовности дублированной системы, |
|
||||||||||||
если r мало: например, коэффициент готовности дублированной |
|
|
|||||||||||
системы с постоянно включенным резервом и r = 0,05 при одной |
|
|
|||||||||||
и двух |
бригадах |
|
обслуживания |
составляет0,995 |
и |
0,998 |
|
|
|||||
соответственно;
5) при малых r вид резервирования практически не влияет на величину коэффициента готовности: например, коэффициент готовности системы при одной бригаде обслуживания составляет 1 и 1 соответственно для случая резервирования с постоянно включенным резервом и замещением.
Следует иметь в виду, что приведенные расчеты наработки на отказ лишь иллюстрируют эффективность резервирования с восстановлением, но не являются достоверными, т.к. в течение
115 лет работы системы интенсивность отказов не может быть величиной постоянной, как это принято при расчетах.
68
Второе задание. Определить вероятность безотказной работы резервированной восстанавливаемой автоматизированной системы.
Пример выполнения задания. Для определения вероятности безотказной работы необходимо составить и решить систему дифференциальных уравнений функционирования системы:
для схемы рис. 2, а:
ì |
|
dp0 (t) |
||
ï |
|
|
|
|
dt |
|
|||
ï |
|
|
||
ídp (t) |
|
|||
ï |
|
1 |
|
= |
ï |
|
dt |
||
î |
|
|
||
для схемы рис. 2, б:
ì dp0 (t)
ï
ï dt
í
ïdp1 (t) = ïî dt
= -2l p0 (t) + m p1 (t); 2l p0 (t) - (l + m) p1 (t);
= -l p0 (t) + m p1 (t);
l p0 (t) - (l + m) p1 (t).
(6)
(7)
За начальные условия примемp0(0)=1, p1(0)=0. Решим системы дифференциальных уравнений (6) и (7) методом РунгеКутты. Введем правые части уравнений системы (5) в аналитическом виде:
m := l
æ1 ö
p := ç ÷
è0 ø
|
é -2×l×p |
0 |
+ m×p |
1 |
ù |
|
D(t1 , p ) := |
ê |
|
|
ú |
||
ê |
|
(l + m)×p |
ú |
|||
|
ë2×l×p 0 |
- |
1û |
|||
Y := Rkadapt(p ,0, 1000,50,D)
69
Решение дифференциальных уравнений приведено на рис. 3.
|
|
|
0 |
|
1 |
|
2 |
|
|
|
35 |
|
700 |
|
|
0.986 |
0.014 |
|
|
|
36 |
|
720 |
|
|
0.986 |
0.014 |
|
|
|
37 |
|
740 |
|
|
0.985 |
0.015 |
|
|
|
38 |
|
760 |
|
|
0.985 |
0.015 |
|
|
|
39 |
|
780 |
|
|
0.985 |
0.015 |
|
|
|
40 |
|
800 |
|
|
0.984 |
0.016 |
|
|
|
41 |
|
820 |
|
|
0.984 |
0.016 |
|
|
Y = |
42 |
|
840 |
|
|
0.983 |
0.017 |
|
|
|
43 |
|
860 |
|
|
0.983 |
0.017 |
|
|
|
44 |
|
880 |
|
|
0.983 |
0.017 |
|
|
|
45 |
|
900 |
|
|
0.982 |
0.018 |
|
|
|
46 |
|
920 |
|
|
0.982 |
0.018 |
|
|
|
47 |
|
940 |
|
|
0.981 |
0.018 |
|
|
|
48 |
|
960 |
|
|
0.981 |
0.019 |
|
|
|
49 |
|
980 |
|
|
0.981 |
0.019 |
|
|
|
50 |
|
1·103 |
|
|
0.98 |
... |
|
|
Рис. 3. |
Результат |
|
решения |
системы |
дифференциальных |
||||
уравнений методом Рунге-Кутты
Аналогичным образом получается решение системы(7):
(рис. 4): m := l
æ1 ö
p := ç ÷
è0 ø
|
é -l×p |
0 |
+ m×p |
1 |
ù |
|
D(t1, p ) := |
ê |
|
|
ú |
||
ê |
|
|
|
|
ú |
|
|
ël×p 0 |
- |
(l + m)×p |
1û |
||
Y := Rkadapt(p ,0, 1000,50,D)
70
