Диагностика и надежность автоматизированных систем. Учебное пособие
.pdfЗадача 10. Для системы, содержащей три основных элемента и один нагруженный резервный, заменяющий любой из основных, составить граф состояний и записать систему дифференциальных уравнений, если вероятность безотказной работы элементов подчиняется экспоненциальному закону, после отказа системы работоспособные блоки отключаются и за время восстановления не отказывают. Задачу решить для случаев:
-интенсивности отказов всех элементов, как и интенсивности восстановлений, одинаковы, восстановление ведут две бригады;
-интенсивности отказов всех основных элементов, как и интенсивности восстановления одинаковы, но отличаются от интенсивности отказов и восстановления резервного элемента, восстановление ведет одна бригада;
-интенсивности отказов всех основных элементов, как и интенсивности восстановления одинаковы, но отличаются от интенсивности отказов и восстановления резервного элемента, восстановление ведут две бригады.
Рекомендации по составлению графа состояний
1.При составлении графа следует соблюдать принцип ординарности, т. е. невозможности возникновения больше одного события (отказа или восстановления) за бесконечно малое время. Это означает, что для двух смежных вершин, соединенных ветвью перехода, состояния системы не могут отличаться состоянием больше чем одного элемента.
2.Граф состояний можно составлять, используя алгоритмы приведенные ниже.
А. Последовательный алгоритм. Начиная с нулевого состояния, наносят все возможные ветви переходов из нулевого состояния и добавляют необходимые вершины, в которых заканчиваются ветви. Затем переходят к первому состоянию, и так последовательно перебирают все добавляемые состояния, пока они не закончатся.
81
Б. Анализируя поведение системы, определяют все возможные её состояния и наносят вершины графа, соответствующие состояниям. После этого прорисовывают ветви переходов между вершинами.
3.С целью уменьшения количества возможных состояний системы для случая, когда после отказа системы работоспособные блоки отключаются и во время восстановления не отказывают, нерезервированные блоки с одинаковой интенсивностью восстановления могут быть объединены в один объект, интенсивность отказа которого равна сумме интенсивностей отказов объединяемых блоков, а интенсивность восстановления объекта равна интенсивности восстановления одного блока.
4.Если две или более вершин графа соединяются только с одной общей для них вершиной и интенсивности переходов в общую вершину (интенсивности восстановления) одинаковы, то эти вершины могут быть объединены в одну. Интенсивность перехода с общей вершины в объединенную будет равна сумме интенсивностей переходов из общей вершины в объединяемые, а интенсивность перехода из объединенной в общую вершину равна интенсивности аналогичного перехода для одной из объединяемых вершин.
5.При определении коэффициента готовности системы с неограниченным количеством восстанавливающих бригад, систему можно разделить на отдельные резервированные и нерезервированные группы, найти коэффициент готовности для каждой группы, а затем, перемножив эти коэффициенты, получить коэффициент готовности системы.
6.Если после отказа резервированной системы ее останавливают и ведут ремонт до полного восстановления, то чтобы не нарушать принцип ординарности вводят дополнительно промежуточные неработоспособные состояния, возникающие при восстановлении, которые могут иметь набор работоспособных элементов, позволяющий системе работать.
82
7. При расчете среднего времени наработки до отказа(до первого отказа), вероятности отказа системы за заданное время, вероятности безотказной работы, гамма-процентной наработки до отказа строят сокращенный граф без выходов из нерабочих состояний.
Пример 19. Для системы, содержащей преобразователь с нагруженным резервным и 5 датчиков температуры, составить граф состояний, если интенсивности отказов преобразователей одинаковы и равны λп, интенсивности отказов первого и второго датчиков равны λ1, остальных – λ2. Интенсивности восстановления (или замены) всех элементов равны μ. Количество восстанавливающих бригад не ограничено, после отказа исправные элементы отключаются.
Объединим все датчики в один блок с интенсивностью отказов равной сумме интенсивностей отказов всех датчиков (все датчики работают одновременно) lд = 2l1 + 3l2 , и интенсивно-
стью восстановления μ. Структурная схема такой системы соответствует схеме, приведенной на рис. 36. Граф состояний показан на рис. 41.
Рис. 41. Граф состояний для примера 19
Если восстанавливает одна бригада, то граф можно упростить, объединив третью и четвертую вершины, т. к. интенсивности восстановления всех элементов одинаковы, а восстанавливать после отказа системы желательно тот элемент, который позволит вернуть систему в рабочее состояние(рис. 42). В третьем состоянии могут быть неработоспособными оба преобразователя либо один преобразователь и один датчик.
83
Рис. 42. Граф состояний при одной восстанавливающей бригаде
При расчете коэффициента готовности или статического коэффициента готовности системы в случае неограниченного количества бригад, систему можно условно разделить на резервированный преобразователь и нерезервированные датчики, для этих двух групп составить отдельно графы, которые показаны на рис. 43, записать системы дифференциальных уравнений, найти коэффициенты готовности для групп, а затем, перемножив их, получить коэффициент готовности системы.
Рис. 43. Графы состояний резервированного преобразователя (а) и нерезервированной группы датчиков (б)
Пример 20. Для расчета коэффициента готовности соста- |
|
||
вить структурную схему и граф состояний системы, состоящей |
|
||
двух основных и одного резервного элемента(рис. 19), если ин- |
|
||
тенсивности отказа и восстановления всех элементов различные. |
|
||
Восстановление ведут две бригады. Резерв нагруженный, после |
|
||
отказа системы исправный блок |
|
||
отключается и в процессе вос- |
|
||
становления |
не |
отказывает. |
|
Включение |
системы |
произво- |
|
дится после восстановления ра- |
|
||
ботоспособного состояния. |
|
|
|
1. Граф состояний по этому |
|
||
условию |
представлен |
на |
|
рис. 44. |
|
|
|
Рис. 44. Граф состояний к примеру 20
84
Построение графа начинают с нулевой вершины, соответствующей всем работоспособным и, в данном примере с нагруженным резервом, работающим элементам. Поскольку отказать может любой работающий элемент и интенсивности отказов у всех элементов различные, а следовательно, и состояния, полученные в результате отказа разных элементов, будут различными, добавляем три новых вершины для состояний 1, 2 , 3 и три ветви перехода с интенсивностями λ1, λ2, λ3.
Впервом состоянии, возникшем в результате отказа первого элемента, работают элементы 2 и 3, которые также могут отказать, а первый элемент может быть восстановлен. В результате восстановления первого элемента система вновь возвратится в нулевое состояние, в результате отказа третьего элемента перейдет в новое, четвертое состояние; в результате отказа второго элемента – в пятое состояние. В четвертом состоянии неработоспособны первый и третий элементы, в пятом – первый и второй. Поскольку число оставшихся работоспособных элементов в этих состояниях меньше количества основных элементов, то эти состояния соответствуют неработоспособной системе и соответствующие им вершины заштрихованы.
Во втором состоянии неработающий второй элемент может быть восстановлен и система возвратится в нулевое состояние; первый элемент может отказать и система перейдет в пятоесо стояние, введенное ранее; может отказать третий элемент и система перейдет в новое, нерабочее шестое состояние.
Аналогичные переходы возможны из третьего состояния.
Вчетвертом состоянии неработоспособны первый и третий элементы и работоспособен второй. Так как по условию задачи при отказе системы исправные элементы отключаются и в- от ключенном состоянии не отказывают, дальнейшего перехода из нерабочего в новое нерабочее состояние не происходит. Если по условию задачи отключение не происходит, граф, как будет показано дальше, продолжает расширяться. Неработоспособные элементы могут быть восстановлены, причем две бригады могут восстанавливать одновременно два элемента, но окончание восстановления все равно происходит в разные моменты времени, и
85
если раньше будет восстановлен первый элемент, система возвратится в третье состояние, а если третий – в первое состояние.
Продолжая анализ всех состояний поочередно, расставляем вершины и ветви графа. Заканчивают построение, когда все введенные вершины будут проанализированы.
2.Если резерв ненагруженный, то при двух восстанавливающих бригадах в графе состояний исчезнет ветвь, выходящая из нулевой вершины в третью, так как в нулевом состоянии третий (резервный) элемент отключен и он не отказывает, но, тем не менее, третье состояние, соответствующее работающим первому и второму элементам и неработоспособному третьему элементу, возникает в результате восстановлений из шестого и четвертого состояний.
3.Если определяется среднее время наработки системы до первого отказа, то все ветви, выходящие из заштрихованных вершин, соответствующих неработоспособным состояниям, следует удалить.
4.Если восстановление производит одна бригада, то из каждых двух ветвей, выходящих из заштрихованных вершин, следует оставить только одну.
5.Если после отказа системы работоспособные блоки не отключаются и восстановление ведут две бригады, граф имеет вид, показанный на рис. 45.
Рис. 45. Граф состояний при возможности отказа всех элементов
86
Седьмое состояние на этом графе соответствует неработоспособным всем трем элементам и переходы в это состояние произошли в результате отказа последнего работоспособного элемента.
6. Если после отказа системы ее отключают и ведут ремонт до полного восстановления, то получим граф, показанный на рис. 46.
Рис. 46. Граф состояний при ремонте до полного восстановления
На данном графе седьмое состояние соответствует только одному неработоспособному элементу, но поскольку ремонт осуществляется до полного восстановления, это состояние, возникшее в результате восстановления одного из двух неработоспособных элементов в нерабочей системе, также является неработоспособным.
Задача 11. Для системы, содержащей контроллер, n основных одинаковых датчиков и m резервных, заменяющих любой из основных, i исполнительных механизмов одинаковых по на-
дежности и с одинаковыми интенсивностями восстановлений μ,
и
составить структурную схему для расчета надежности и граф состояний для расчета коэффициента готовности, если после отказа системы работоспособные блоки отключаются и при восстановлении системы их отказ не происходит. Количество восстанавливающих бригад r. Данные по вариантам приведены в табл. 10.
87
|
|
|
|
|
|
Таблица 10 |
|
№ |
|
|
|
|
Интенсивности вос- |
|
|
n |
m |
i |
r |
становления основных |
Резерв |
||
вар |
|||||||
|
|
|
|
и резервных датчиков |
|
||
|
|
|
|
|
|
||
1 |
3 |
1 |
3 |
1 |
Одинакова |
Нагруженный |
|
2 |
3 |
2 |
2 |
1 |
Различается |
Нагруженный |
|
3 |
3 |
2 |
3 |
2 |
Одинакова |
Ненагруженный |
|
4 |
4 |
2 |
2 |
1 |
Различается |
Ненагруженный |
|
5 |
4 |
2 |
2 |
2 |
Одинакова |
Нагруженный |
|
6 |
4 |
1 |
3 |
2 |
Различается |
Нагруженный |
|
7 |
4 |
2 |
3 |
3 |
Одинакова |
Ненагруженный |
|
8 |
6 |
2 |
3 |
3 |
Различается |
Ненагруженный |
|
9 |
6 |
3 |
2 |
1 |
Одинакова |
Нагруженный |
|
10 |
6 |
2 |
2 |
1 |
Различается |
Нагруженный |
|
11 |
6 |
2 |
2 |
2 |
Одинакова |
Ненагруженный |
|
12 |
3 |
2 |
4 |
2 |
Различается |
Ненагруженный |
|
13 |
4 |
3 |
4 |
3 |
Одинакова |
Нагруженный |
|
14 |
4 |
3 |
4 |
3 |
Различается |
Нагруженный |
|
15 |
4 |
3 |
3 |
1 |
Одинакова |
Ненагруженный |
|
16 |
4 |
3 |
3 |
1 |
Различается |
Ненагруженный |
|
17 |
5 |
2 |
3 |
2 |
Одинакова |
Нагруженный |
|
18 |
5 |
2 |
2 |
2 |
Различается |
Нагруженный |
|
19 |
5 |
2 |
2 |
3 |
Одинакова |
Ненагруженный |
|
20 |
5 |
2 |
2 |
3 |
Различается |
Ненагруженный |
5.4. Определение коэффициента готовности восстанавливаемой системы
Как указывалось ранее, коэффициент готовности, определяющий вероятность работоспособного состояния системы в момент времени t, равен сумме вероятностей работоспособных состояний. Чтобы найти коэффициент готовности необходимо:
-составить структурную схему системы для расчета надеж-
ности;
-разработать граф состояний системы;
-записать систему дифференциальных уравнений;
-решить систему уравнений для выбранных начальных ус-
ловий;
-сложить вероятности работоспособных состояний.
88
Известно несколько методов решения систем обыкновенных дифференциальных уравнений, записанных в форме Коши. В компактной матричной форме такую систему можно записать как
dP |
= A × P(t ), |
(61) |
|
||
dt |
|
|
где А – матрица постоянных коэффициентов дифференциальных уравнений; P(t) - вектор переменных вероятностей.
Если система уравнений простая, можно привести ее к системе независимых уравнений большего порядка и решать уравнения только для работоспособных состояний.
Так, из системы уравнений (60)
dP0 = -2lP0 (t )+ mP1(t ), dt
dP1 = -(m + l)P1 t( +) 2lP0 (t )+ mP2 (t ), dt
dP2 = -mP2 (t )+ lP1(t ), dt
считая, что в нулевой момент времени оба объекта системы, показанной на рис. 33, были работоспособны, и она находилась в нулевом состоянии, необходимо найти функции P0(t) и P1(t) при начальных условиях
P0(0)=1, P1(0)=0, P2(0)=0. |
(62) |
Для этого выразим из первого уравнения системы P1 (t ) и подставим во второе и третье:
P1 (t )= |
1 |
|
[P0' (t )+ 2lP0 (t )], |
|
|
|
(63) |
|||
m |
|
[P0' (t )+ 2lP (t )]+ |
||||||||
|
1 |
[P0'' (t )+ 2lP0' (t )]= -(l + m) |
1 |
|
||||||
|
|
|
(64) |
|||||||
|
m |
|
|
|
m |
0 |
||||
+ 2lP0 (t )+ mP2 (t ), |
|
|
|
|
||||||
P2' (t) = l |
|
1 |
[P0' (t )+ 2lP0 (t )]- mP2 (t ). |
|
(65) |
|||||
|
|
|
||||||||
|
|
|
|
|
m |
|
|
|
|
|
Упрощая уравнение (64) и выразив из него P2 (t )
89
|
|
|
P (t )= |
|
P |
'' (t ) |
|
|
|
|
|
|
|
|
|
|
|
|
P' (t ) |
|
2l2 |
|
|
(t ), |
|
||||||||||
|
|
|
|
0 |
|
|
+ (3l + m) |
0 |
|
+ |
|
|
|
|
P |
(66) |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
2 |
|
|
|
|
m2 |
|
|
|
|
|
|
|
|
|
|
|
|
m2 |
|
m2 |
0 |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
подставим |
|
P2 (t ) в (65) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
[P0' |
(t )+ 2lP0 (t )]- |
|
|||||||||||||
|
P0''' (t) |
+ (3l + m) |
P0'' (t) |
+ |
2l2 |
|
P0' (t )= l |
1 |
|
|
|||||||||||||||||||||||||
|
|
|
|
|
m2 |
|
m |
|
|||||||||||||||||||||||||||
|
m2 |
|
|
|
|
|
|
|
|
|
|
m2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
P '' |
(t ) |
|
|
|
|
|
|
|
|
P |
' |
(t ) |
|
|
2l2 |
|
|
|
|
|
|
|
|
||||||||||
- |
0 |
|
|
|
- (3l + m) |
|
0 |
|
|
|
|
- |
|
|
|
P (t ). |
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
m |
|
|
|
|
|
|
|
|
|
|
|
|
m |
|
|
|
|
|
|
m |
0 |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Объединяя подобные слагаемые, получим уравнение для P0 (t ): |
|||||||||||||||||||||||||||||||||||
|
|
P0''' (t )+ (3l + 2m)× P0'' (t )+ (2l2 |
+ 2lm + m2 )× P0' (t )= 0 . |
(67) |
|||||||||||||||||||||||||||||||
Записываем для уравнения (67) характеристическое уравнение |
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
k 3 + (3l + 2m)× k 2 + (2l2 + 2lm + m 2 )× k = 0 , |
|
||||||||||||||||||||||||||||
находим корни |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
k1 = 0 , |
|
|
|
|
|
|
|
|
|
|
|
|
|
3(l + 2m)2 - 4(2l2 + 2lm + m2 ) |
|
||||||||||||||||||||
k 2,3 = - (3l + 2m)± |
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
||
или k 2,3 = |
- (3l + 2m)± l2 + 4lm |
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Решение уравнения (67) получим в виде |
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
P0 (t )= C1 + C 2 exp (k 2 t )+ C 3 exp (k 3 t ). |
(68) |
||||||||||||||||||||||||
Для определения постоянных Сi |
необходимо найти началь- |
||||||||||||||||||||||||||||||||||
ные значения P ' (t ) |
и P '' |
(t ) при t = 0 , воспользовавшись началь- |
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ными условиями (62) и уравнениями (63), (64) при t = 0 . |
|
||||||||||||||||||||||||||||||||||
Тогда |
' (0)= -2l , |
|
|
|
|
(0)= 2l(3l + m) - 2l2 = 2l(2l + m). |
|
||||||||||||||||||||||||||||
P (0)= 1 , P |
|
P |
'' |
|
|||||||||||||||||||||||||||||||
0 |
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Выражение |
(68) |
|
|
с учетом начальных условий дает три |
|||||||||||||||||||||||||||||||
уравнения для поиска постоянных С: |
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
P0 (0)= C1 + C 2 + C 3 , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
P ' (0)= C |
2 |
k |
2 |
+ C |
3 |
k |
3 |
|
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
90
