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

УчебПособие (Теория надежности)2011

.pdf
Скачиваний:
105
Добавлен:
07.03.2015
Размер:
2 Мб
Скачать
☆

Составим диаграмму пространства состояний для данной системы, которая имеет три состояния: 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