УчебПособие (Теория надежности)2011
.pdf
Составим диаграмму пространства состояний для данной системы, которая имеет три состояния: S5, S6, S7 (рис 5.11).
1 (q6 7 |
q6 5) t |
t |
|
|
q6 7 |
S7 |
|
|
|
|
|
|
S6 |
q7 6 t |
|
1 q7 6 t
q6 5 t
q5 6 t
1 q5 6 t
S5
Рис. 5.11. Диаграмма пространства состояний системы
Определим интенсивности переходов из состояния в состояние.
|
|
|
|
|
|
|
|
|
Pj ( ) qJ i |
|
|
|||||
qi J |
= qi j ; qJ i = |
|
j J |
|
|
|
; |
|
||||||||
|
|
Pj ( ) |
||||||||||||||
|
|
j J |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
j J |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
Pi ( ) qi j |
|
|
Pj ( ) q j i |
||||||||||
qI J = |
i I |
|
j J |
|
; qJ I |
= |
j J |
i I |
. |
|||||||
|
|
Pi ( ) |
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
Pj ( ) |
|||||
|
|
|
|
|
i I |
|
|
|
|
|
|
|
|
j J |
||
|
|
|
|
P ( ) q |
P ( ) q |
|
|
|
|
|
||||||
q |
|
1 |
|
1 5 |
2 |
|
|
2 5 |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|||||||
6 5 |
|
|
|
|
P ( ) P ( ) |
|
|
|
|
|
||||||
|
|
|
|
|
|
1 |
2 |
|
|
|
|
|
|
|
|
|
|
0,514 3,891 0,194 10,309 |
5, 65 год-1; |
||||||||||||||
|
||||||||||||||||
|
|
|
|
0,514 0,194 |
|
|
|
|
|
|
|
|
||||
q |
q |
|
23, 256 год-1; |
|
|
|
|
|
||||||||
5 6 |
|
|
5 2 |
|
|
|
|
|
|
|
|
|
|
|
||
81
|
|
|
P ( ) q |
P ( ) q |
|||||
q |
|
|
1 |
1 3 |
2 |
2 3 |
|
||
|
|
|
|
|
|
||||
6 7 |
|
|
P ( ) |
P ( ) |
|
|
|
|
|
|
|
|
|
1 |
2 |
|
|
|
|
|
0, 514 5,832 0,194 30, 928 |
12, 709 год-1; |
|||||||
|
|
0,514 0,194 |
|||||||
|
|
|
|
|
|
|
|||
q |
|
|
P3 ( ) q3 2 P4 ( ) q4 2 |
|
|||||
|
|
|
|||||||
7 6 |
|
|
P3 ( ) P4 ( ) |
|
|
|
|
||
|
|
|
|
|
|
|
|
||
|
0, 051 117, 647 0, 068 44,118 |
75, 63 год-1. |
|||||||
|
|||||||||
|
|
|
|
0, 051 0, 068 |
|
|
|
|
|
Определим вероятности нахождения системы в каждом из трех состояний в стационарном режиме работы.
Система уравнений Чемпена – Колмогорова:
P5 |
( ) |
=0=( q5 6 ) P5 ( ) q6 5 P6 ( )+0 P7 ( ); |
|||
|
t |
||||
|
|
|
|||
P6 |
( ) |
|
|
||
|
=0=q5 6 P5 ( ) (q6 5 q6 7 ) P6 ( )+q7 6 P7 ( ); |
||||
|
t |
||||
|
|
|
|||
P7 |
( ) |
|
|
||
|
=0=0 P5 ( ) q5 6 P6 ( )-q7 6 P7 ( ). |
||||
|
t |
||||
|
|
|
|||
Условие нормировки: |
|||||
P ( ) P ( ) P ( )=1. |
|||||
5 |
|
6 |
7 |
||
Тогда |
|
|
|||
dP5 |
( ) |
0 |
( q5 6 ) P5 ( ) q6 5 P6 ( ) 0 P7 ( ); |
||
|
dt |
||||
|
|
|
|||
|
( ) |
|
|
||
dP7 |
0 |
0 P5 ( ) q6 7 P6 ( ) q7 6 P7 ( ); |
|||
|
|
|
|||
dt |
|||||
|
|
|
|||
P ( ) P ( ) P ( ) 1. |
|||||
|
5 |
|
6 |
7 |
|
|
|
|
|
|
|
Подставив в полученную систему уравнений значения интенсивностей переходов из состояния в состояние, определим вероятности нахождения системы в состояниях S5, S6 и S7.
Вероятность нахождения системы в работоспособном состоянии для стационарном режиме работы равна P6( ) = 0,709, в состоянии аварийного ремонта – P7( ) = 0,119 и в состоянии плановопредупредительного ремонта – P5( ) = 0,172.
Ответ:
1. 1 = 9 год-1; 2 = 17 год-1; 3 = 9 год-1; 4 = 3 год-1; 5 = 4 год-1;
82
P1( ) = 0,514; P2( ) = 0,194; P3( ) = 0,051; P4( ) = 0,068; P5( ) = = 0,172;
q1 2 = 7,782 год-1; q2 1 = 46,392 год-1; q1 3 = 5,832 год-1; q2 3 =
= 30,928 год-1; q1 5 = 3,891 год-1; q2 5 = 10,309 год-1; q3 2 = 117,647 год-1; q4 2 = 44,118 год-1; q3 4 = 58,824 год-1; q5 2 = 23,256 год-1.
2. P6( )=0,709; P7( )=0,119; P5( )=0,172.
Пример 4. Диаграмма пространства состояний системы представлена на рис. 5.12.
1 (q1 2 q1 3
q1 4 q1 5 ) t
|
S1 |
|
t |
t |
|
|
|
|
1 |
2 |
|
|
2 |
1 |
q |
q |
|
|
S2 |
1 (q2 1 q2 3
q2 4 q2 5 ) t
|
|
|
2 |
|
|
S3
|
|
|
q |
4 |
|
t |
|
|||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
q |
|
|
|
t |
|||
|
4 |
1 |
||||||||
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5 |
|
|
|
|
|
|
|
|
|
q |
q |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5 |
|
|
5 |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
1 |
|
|
t |
||
2 |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
||||
5 |
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
2 5 |
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
S5 |
|
|
|
|
|
|
|
|
|
|
|
|
1 (q3 1 q3 2
q3 4 q3 5 ) t
1 (q4 1 q4 2
q4 3 q4 5 ) t
S4
|
t |
|
|
|
|
|
|
|
|
|
|
|
5 |
|
|
|
|
|
|
q |
4 |
|
t |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
4 |
|
|
|
|
|
q |
5 |
|
|
|
||
1 (q5 1 q5 2
q5 3 q5 4 ) t
Рис. 5.12. Диаграмма пространства состояний системы
Известно следующее:
-продолжительность пребывания системы в каждом из состояний подчинена экспоненциальному закону;
-все интенсивности переходов из состояния в состояние имеют одинаковые значения:
83
q1 2 = q1 3 = q1 4 = q1 5 = q2 1 = q2 3 = q2 4 = q2 5 = q3 1 = q3 2
=q3 4 = q3 5 = q4 1 = q4 2 = q4 3 = q4 5 = q5 1 = q5 2 = q5 3 = q5 4 =
=2 год-1;
-S1 S2 S5=S6 – работоспособное состояние; S3 S4=S7 – неработоспособное состояние;
-в начальный момент времени система находится в состоянии S1.
Задание:
Преобразовав исходную диаграмму пространства состояний
вдиаграмму с двумя состояниями, определить:
1)вероятность нахождения системы в работоспособном состоянии в момент времени t = 0,1 года;
2)вероятность безотказной работы в момент времени t = 0,5
года.
Решение:
1. Определение вероятности нахождения системы в работоспособном состоянии в момент времени t = 0,1года.
Преобразуем исходную диаграмму пространства состояний в диаграмму с двумя состояниями (рис. 5.13).
1 q7 6 t
1 q6 7 t
S6= |
|
q6 7 t |
||
|
|
|
||
=S1 S2 |
S5 |
|
S7=S3 S4 |
|
q7 6 |
t |
|||
|
|
|||
Рис. 5.13. Диаграмма пространства состояний для определения вероятности нахождения системы в работоспособном состоянии в момент времени t = 0,1 года
На рис. 5.13 S6 – работоспособное, а S7 – неработоспособное состояния системы.
Определим интенсивности переходов из состояния в состояние.
84
|
|
P (t) (q |
q |
) P (t) (q |
|
q |
) P (t) (q |
q |
) |
|
|||||
q |
|
1 |
1 3 |
1 4 |
2 |
2 3 |
2 4 |
5 |
5 3 |
|
5 4 |
|
; |
||
|
|
|
|
|
|
|
|
|
|
|
|
||||
6 7 |
|
|
|
|
P (t) P (t) P (t) |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
1 |
2 |
5 |
|
|
|
|
|
|
|
q |
|
P3 (t) (q3 1 q3 2 q3 5 ) P4 (t) (q4 1 q4 2 q4 5 ) |
. |
|
|
||||||||||
|
|
|
|
|
|||||||||||
7 6 |
|
|
|
|
P3 (t) P4 (t) |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
Так как условие объединяемости выполняется, то |
|
|
|
|
|
|||||||||
|
q6 7 q1 3 q1 4 |
2 2 4 ( |
1 |
); |
|
|
|
|
|
|
|
||||
|
год |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
q7 6 q3 1 q3 2 |
q3 5 2 2 2 6 ( |
1 |
). |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
год |
|
|
|
|
|
|
Составим систему дифференциальных уравнений Чемпена – Колмогорова с учетом условия нормировки и определим P6РС(t = 0,1).
Уравнения Чемпена – Колмогорова:
|
Р6РС |
(t) |
q6 7 |
P6РС (t) q7 6 |
P7НС (t); |
||
|
t |
|
|||||
|
|
|
|
|
|
|
|
|
Р |
(t) |
|
|
|
|
|
|
7РС |
|
q6 7 P6РС (t) q7 6 |
P7НС (t). |
|||
|
t |
|
|||||
|
|
|
|
|
|
|
|
Условие нормировки: P |
(t) P |
|
(t) 1. |
||||
|
|
|
|
6РС |
7НС |
|
|
Тогда |
|
|
|
|
|
|
|
|
Р6РС (t) |
q6 7 |
P6РС (t) q7 6 P7НС (t); |
||
|
|
|
|||
t |
|
||||
|
|
|
|
|
|
|
Р |
(t) P |
(t) 1. |
||
|
6РС |
|
7РС |
|
|
|
Р6РС (t) |
q6 7 |
P6РС (t) q7 6 (1 Р6РС (t)); |
||
|
|
|
|||
t |
|
||||
|
|
|
|
|
|
P |
(t) 1 Р |
|
(t). |
||
|
7НС |
|
6РС |
|
|
Из первого уравнения найдем P6(t) методом разделения переменных:
85
|
dP6РС (t) |
q6 7 P6РС (t) q7 6 (1 |
P6РС (t)); |
|
||||||||||||||||
|
|
dt |
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dP6РС (t) |
|
|
|
dt; |
|
|||||||||||
|
P6РС (t) (q6 7 q7 6 ) q7 6 |
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
1 |
|
|
|
ln |
|
P6РС (t) (q6 7 q7 6 ) q7 6 |
|
t C; |
||||||||||||
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|||||||||||||||
|
(q6 7 q7 6 ) |
|
|
|
|
|
|
(q6 7 q7 6 ) (t C); |
||||||||||||
ln |
P6РС (t) (q6 7 q7 6 ) q7 6 |
|||||||||||||||||||
|
P6РС (t) (q6 7 q7 6 ) q7 6 |
|
e |
(q6 7 q7 6) (t C) |
. |
|||||||||||||||
|
|
|||||||||||||||||||
|
|
В этом случае уравнение имеет два решения: |
|
|||||||||||||||||
|
|
- первое общее решение: |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
P |
|
(t) |
q7 6 e |
(q6 7 q7 6) (t C) |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
; |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
6РС |
|
|
|
|
q6 7 q7 6 |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
- второе общее решение: |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
P |
|
(t) |
q7 6 e |
(q6 7 q7 6) (t C) |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
. |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
6РС |
|
|
|
|
q6 7 q7 6 |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
Найдем частные решения, |
|
используя начальное условие |
||||||||||||||||
P1(t = 0) = P6РС(t = 0) = 1. Первого частного решения не существует. Второе частное решение находится следующим образом:
|
P6РС (0) 1 |
q7 6 e |
(q6 7 q7 6 ) (0 C) |
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
; |
|
|
|
|
|
|
|||||||||
|
|
|
|
q6 7 |
q7 6 |
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
e |
(q6 7 q7 6 ) C |
(q6 7 q7 6 ) q7 6 ; |
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
e |
(q6 7 q7 6 ) C |
q6 7 ; |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
(q6 7 q7 6 ) C ln(q6 7 ); |
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
C |
|
ln(q6 7 ) |
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
(q6 7 q7 6 ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
Тогда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
(q6 7 q7 6) (t |
ln(q6 7 ) |
) |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
(q6 7 q7 6) |
|
|
|
|
|
|||||||||||
P |
|
|
(t) |
q7 6 e |
|
|
|
|
|
|
|
|
; |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
6РС |
|
|
|
|
|
|
q6 7 q7 6 |
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
q7 6 |
|
|
|
|
|
1 |
|
(q |
q |
) t |
|
ln(q |
) |
|
|||
P |
|
|
(t) |
|
|
|
|
(1 |
|
|
e |
6 7 |
|
7 6 |
|
|
e |
6 7 |
|
) ; |
||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
6РС |
|
|
q6 7 q7 6 |
|
|
q7 6 |
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
86 |
|
|
|
|
|
|
|
|
|
|
|
q7 6 |
|
q6 7 |
|
(q |
q |
) t |
|
|
Р6РС (t) |
|
|
(1 |
|
e |
6 7 |
7 6 |
|
) – частное |
(q6 7 |
q7 6 ) |
q7 6 |
|
||||||
|
|
|
|
|
|
|
|||
решение.
Вероятность нахождения системы в работоспособном состоянии в момент времени t = 0,1 года равна
P |
(0,1) |
|
6 |
(1 |
4 |
e (4 6) 0,1) 0,747 год-1 . |
|
|
|
||||
6РС |
4 |
6 |
6 |
|
||
|
|
|||||
2. Определение вероятности безотказной работы в момент времени t = 0,5 года.
Составим диаграмму пространства состояний для определения вероятности безотказной работы в момент времени t = 0,5 года.
При определении вероятности безотказной работы системы принимают, что восстановление работоспособности системы из неработоспособного состояния невозможно. Диаграмма пространства состояний в этом случае выглядит следующим образом (рис. 5.14):
1 q6 7 t
S 6 |
= |
q6 7 t |
|
S7= S3 S4 |
|||
= S1 S2 S5 |
|||
|
|||
Рис. 5.14. Диаграмма пространства состояний для определения вероятности безотказной работы системы в момент времени t = 0,5 года
На рис. 5.14 S6 – работоспособное, а S7 – неработоспособное состояния системы, q6 7 = 4год-1.
Составим систему дифференциальных уравнений Чемпена – Колмогорова с учетом условия нормировки и определим P7БР(t = 0,5 лет).
87
|
P6БР |
(t) |
q6 7 |
P6БР (t); |
|
|
t |
|
|
||
|
|
|
|
|
|
|
P |
(t) |
|
|
|
|
7ВО |
|
|
q6 7 P6БР (t). |
|
|
t |
|
|
||
|
|
|
|
|
|
Так как при одной неизвестной имеем одно уравнение, то условие нормировки можно не использовать.
P6БР (t) q6 7 P6БР (t);
t
P6БР (t) q6 7 t;
P6БР (t)
P6БР (t) q6 7 t;
P6БР (t)
ln P6БР (t) q6 7 t C;
P6БР (t) e q6 7 t C .
В этом случае уравнение имеет два решения: - первое общее решение:
P (t) е q6 7 t C ;
6БР
- второе общее решение:
P |
(t) е q6 7 t C (уравнение решений не имеет, так как |
6БР |
|
вероятность должна быть положительной).
Найдем частное решение, используя начальные условия
P1(t = 0) = P6БР(t = 0) = 1:
P |
(0) 1 е q6 7 0 C |
С 0. |
6БР |
|
|
Тогда P6БР (t) е q6 7 t |
– частное решение. |
|
Вероятность безотказной работы в момент времени t = 0,5 лет:
P (0,5) е 4 0,5 0,135.
6БР
Ответ:
1.Р6РС (2) 0,747 ;
2.P (0,5) 0,135.
6БР
Пример 5. Два одинаковых объекта (А, В) образуют дублированную систему с резервированием замещением и
88
ограниченным аварийным ремонтом. Планово-предупредительные ремонты на объектах не проводятся.
Известно следующее:
-объекты имеют независящие от времени интенсивности
переходов: интенсивности отказов: λА = λВ = λ = 2 год-1; интенсивности восстановления из аварийного ремонта: μАР,А = μАР,В = µАР = 10 год-1;
-в начальный момент времени система находится в состоянии S1.
Задание:
1)Определить для стационарного режима работы: вероятность работоспособного состояния (РС), вероятность неработоспособного состояния (НС) и вероятность промежуточного состояния (ПС) системы, в котором один из элементов находится в работоспособном, а другой – в неработоспособном состоянии.
2)Составить диаграмму пространства состояний и систему уравнений Чемпена – Колмогорова с условием нормировки для определения вероятности безотказной работы системы.
Решение:
Составим диаграмму пространства состояний для определения вероятностей нахождения системы в работоспособном PРС( ), неработоспособном PНС( ) и промежуточном PПС( ) состояниях (рис.
5.15).
1 t 1 ( АР ) t 1 АР t
S1 a b |
t |
|
t |
S2 a b |
S5 a b |
|
|
|
|
|
|
|
|
Р |
|
|
|
|
|
|
|
|
|
|
|
O |
АР |
|
|
O’ |
|
|
|
|
|
|
|
|
|
|
|
|
|
Р |
|
S3 a b |
|
|
S4 a b |
S b a |
t |
|
|
||
|
|
t |
6 |
|
|
|
|
|
|
1 t |
|
|
1 ( АР ) t |
1 АР t |
|
|
|
Рис. 5.15. Диаграмма пространства состояний для определения вероятности нахождения системы в работоспособном состоянии (дублированная система с резервированием замещением и ограниченным аварийным ремонтом)
89
Здесь приняты следующие обозначения: а, b – рабочие
состояния объектов А или В; a , b – неработоспособные состояния А или В (в этом состоянии на объекте проводится аварийный ремонт либо объект ожидает его начала); a, b − режим ожидания объекта А или В.
Система может находиться в шести состояниях: S1 = a b – состояние, в котором А находится в рабочем состоянии, а B – в режиме ожидания; S2 = a b – состояние, в котором B находится в рабочем
|
|
|
|
|
|
|
|
|
состоянии, а A – в режиме ожидания; |
S3 a b – |
состояние, |
в |
|||||
котором А |
находится в неработоспособном состоянии (состоянии |
|||||||
|
|
|
|
|
|
|
||
аварийного |
ремонта), а B – в рабочем |
состоянии; |
S4 a b |
– |
||||
состояние, в котором B находится в неработоспособном состоянии (состоянии аварийного ремонта), а A – в рабочем состоянии;
S5 a b – состояние, в котором система находится в
неработоспособном состоянии. При этом на объекте A производится аварийный ремонт, так как он раньше отказал. Объект B подвергнется
аварийному ремонту после окончания работ на объекте А; S6 b a –
состояние, в котором система находится в неработоспособном состоянии. При этом на объекте B производится аварийный ремонт, так как он раньше отказал. Объект A подвергнется аварийному ремонту после окончания работ на объекте B.
При этом
P ( ) 1 P |
( ) P ( ) P ( ) P ( ) P ( ); |
||||
РС |
НС |
1 |
2 |
3 |
4 |
PПС ( ) P2 ( ) P4 ( );
PНС ( ) P5 ( ) P6 ( ).
Составим систему дифференциальных уравнений Чемпена – Колмогорова с учетом условия нормировки для стационарного режима и определим вероятности нахождения системы в работоспособном РРС(∞), неработоспособном РНС(∞) и промежуточном РПС(∞) состояниях.
90
