![](/user_photo/_userpic.png)
книги из ГПНТБ / Коваленко И.Н. Полумарковские модели в задачах проектирования систем управления летательными аппаратами
.pdfВ процессе работы ЗЭ восстановлений или замен отказавших элементов не производится, а осуществляет ся лишь автоматическое подключение резёрвных эле ментов вместо отказавших основных. Зарезервирован ный элемент считается отказавшим при выходе из строя одного из основных элементов и при отсутствии исправ ных резервных элементов.
Задача |
заключается в |
определении оптимального |
|||||
периода Тк |
профилактического |
обслуживания ЗЭ, при |
|||||
назначении |
которого |
обеспечивается |
максимальная ве |
||||
роятность застать ЗЭ |
в исправном-состоянии в про |
||||||
извольный момент времени t |
(максимальная готовность) |
||||||
и максимальная |
вероятность |
безотказной |
работы ЗЭ |
||||
за время от t до |
t + x |
(максимальная надежность). |
|||||
Обозначим, как и ранее, вероятность того, что ЗЭ ра |
|||||||
ботает к моменту t и продолжает |
работать |
в течение |
времени от t до t+x, через р(х, t), а функцию распре деления времени безотказной работы ЗЭ (случайной ве личины X),через F(t) = Р {X<t} .
Вначале предположим, что время между двумя по следовательными моментами профилактического обслу живания ЗЭ является случайной величиной Y с функци ей распределения G (t) = Р {У </}.
Так как после каждого профилактического обслужи вания ЗЭ возвращается в первоначальное состояние, то времена между двумя моментами окончания профилак тических обслуживании ZKобразуют вновь процесс вос становления.
Величина Z может принимать следующие значения:
|
|
‘У + То, |
если |
k=±0 |
|
|
Y + T» |
если |
k = \ |
|
Z = |
г + т „ |
если |
k— 2 |
|
|
У-\-Тп, если |
k ~ n , |
|
где k = \ , 2,..., |
n — количество отказавших элементов к |
|||
началу профилактического обслуживания. |
||||
В рассматриваемом случае также справедлива фор |
||||
мула (4.6), в которой |
|
|
||
[г] = Ж |
[У] + |
Ро (К) Т0+ рг (У) Гх + . . . + рп (У) Тп, |
140
причем У]/Л-(^) = 1 и P i ( Y ) , і = О, 1, 2,. . п обозначает
; =0 вероятность к началу профилактического обслуживания
иметь і отказавших элементов.
Выберем G ( t ) |
по отношению к Т к так, чтобы выпол |
||
нялось условие |
(4.10). |
Тогда на |
основании выраже |
ния (4.6) |
|
|
|
|
I |
[1 — F (о) -р ж)] dm |
|
Pix, Тк)= |
|
+ Рі ( Л <)Ті + |
. (4.32) |
тк+ Р ъ { Т к) п |
- Р П( Т К) Т П |
Формула (4.32) дает возможность выбирать период Тк профилактического обслуживания ЗЭ, исходя из за данной величины р(х, ТІ(). Однако, наибольший ин терес представляет отыскание оптимальной величины Тѵ, обеспечивающей максимум функции р(х, Тк) и опре деляемой из условий
д р ( х , Т „■) |
(4.33) |
|
д Т к |
||
|
при существовании (для каждого конкретного случая) единственности максимума р{х, Тк).
Для определенности рассмотрим ЗЭ (а для получе ния окончательных результатов перейдем впоследствии к еще более частным случаям), состоящий из одного ос новного элемента (m= 1) и п—1 резервных элементов, которые находятся в состоянии ненагруженного резерва. Тогда время безотказной работы ЗЭ
Х ^ Х 1+ Х 2+ . . . + Х П, |
(4.34) |
:где A4 — время безотказной работы основного элемента; Х2— время безотказной работы первого резервного
элемента; Хп — время безотказной работы (п—1)-го резервно
го элемента. Следовательно,
У-\-Т0, если А'1> К ;
К + 7 \, если А ^ < Г , Х г-\-Х%> У \
Z = \Y-VT2, если Х г+ Х 2< Г , Х ^ Х ъ + Х ^ У і ( 4 .3 5 )
Y -\-Tn, если X < yY .
141
Очевидно, для случая ненагружепного резерва
F{i)=1 |
It- |
(АО2 I |
I (Ар"-1 |
Л—X/ (4.36) |
|
|
21 |
0 - 1 ) 1 |
J |
при условии, что интенсивность отказов каждого элемен та остается постоянной п равной X, а характер потока отказов элемента простейший. В этом случае выраже ние (4.32) принимает вид
Р{х, Тк)=
|
|
\ + \{t + |
x ) + |
[A U + х)]2 |
|
||
|
|
|
- 21— - + • •. + |
|
|||
|
—XT |
—XT |
|
|
(ХТ'к)'Л—1 |
—хг. |
|
Г к + Г 0 е |
т*+Т{кТк е |
* + . . . + Тп^ |
|
|
|||
|
|
|
|
|
|
( И - 1)1 |
|
|
4" |
[ \ ( і + х ) ] п 1 j e- x { t + . r ) dt |
|
||||
|
|
0 - 1 ) ! |
|
|
|
|
(4.37) |
|
|
|
|
|
|
—XT. |
|
+ Тп |
1 - |
1 + XT’к + ... |
+ |
|
|
||
0 - 1 ) 1 |
|
Рассмотрим для примера случай: »z=l, я=2, т. е. когда ів состав системы входят один основной и одни ре зервный элемент и для каждого элемента поток отказов простейший с параметром X. Тогда
’ |
К , |
[А (/ + |
л-)]2 |
) ^ _ Х(І+Л) |
^ |
І1 + |
l ( t + X ) + |
|
1 e“ A(f+'r; d t |
P{X, Тк) = _______ о _ _ ____________________________ ” ______________________________ |
||||
Тк + Т0 |
+ ТХ1ТКе ~ ХТ« + |
Т2 {1 |
- [1 + Х Г К] е _ Х Г к ) ‘ |
|
|
|
|
|
(4. 38) |
Убеждаясь в существовании единственного максиму ма функции р{х, Тк) и используя условие (4.33), полу чаем следующее уравнение:
е - хѴ [27\.(1 + ^ ) + Г2 (1—ХГК)—2 ( 7 \ - 7 ’0) +
+ м ^ 2+ - ^ ) ] - |
(Та- 2 7 \+ Тг) — у = 0 . |
(4. 39) |
Уравнение (4.39) легко решается графически (рис. 4.7). При этом удобно все члены, зависящие от Г1(, обо
И2
значить через функцию Ф(Г,<). В табл. 4.5 приведены ре зультаты вычислении значений ТІ;оѵі для случаев т=1;
/г = 2, |
3, 4 при различных значениях X. |
Вычисления про |
||||||||
изводились с нслользоваиием |
выражения (4.37) |
и ус- |
||||||||
ловия |
(4.33). |
|
|
|
ФШ) |
|
|
|||
Ясно, что при п> 2 |
зна |
|
|
|||||||
|
|
|
|
|||||||
чительно |
усложняются |
вы |
|
|
|
|
||||
кладки |
|
при |
определении |
|
|
|
|
|||
оптимального |
значения |
ве |
|
|
|
|
||||
личины |
Тк, однако, |
принци |
|
|
|
|
||||
пиальных |
затруднений |
не |
|
|
|
|
||||
встречается. |
сложных |
слу |
Рис. |
4.7. Графическое |
опреде |
|||||
В |
более |
ление |
оптимального |
периода |
||||||
чаях |
(нагруженный, |
облег- |
контроля |
САУ с резервом |
||||||
чеиныи резерв, а также при |
|
|
|
|
||||||
^ ^ const) |
вероятности р,- могут быть определены |
с по |
||||||||
мощью |
биномиального |
распределения |
(нагруженный |
резерв) или при применении так называемой «схемы ги- ■бели и размножения» к исследованиям резервных систем,
которое |
было |
предложено |
и |
обосновано |
А. Д. |
|
Со |
|||
ловьевым. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица |
4.5 |
|||
X. •/„ |
т0. ч |
Г,, ч |
Т2 j ч |
То, ч |
Т,, ч |
|
Г к орт’ Ч |
= |
|
|
/1 = 2 |
/1 = 3 |
л |
4 |
|||||||
|
|
|
|
|
|
т - 1 |
ІИ —1 |
т |
|
|
|
|
|
|
|
|
|
|
|
- 1 |
|
0,04 |
0,2 |
0,4 |
0,6 |
0,8 |
1 |
8,7 |
12,7 |
|
17,3 |
|
0,02 |
0,3 |
0,6 |
0,9 |
1,2 |
1,5 |
14,1 |
23,8 |
|
32,9 |
|
0,01 |
0,5 |
1 |
1,5 |
9 |
2,5 |
26 |
45,9 |
|
63,6 |
|
0,005 |
0,7 |
1,5 |
2,4 |
3,2 |
4 |
45,9 |
83,6 |
117,5 |
||
0,0025 |
1 |
9 |
3 |
4 |
5 |
81,8 |
155,9 |
232,2 |
||
0,00125 |
2 |
3 |
4 |
5 |
6 |
163,7 |
316 ,5 |
488,7 |
||
0,00066 |
4 |
О ,0 |
7 |
8,5 |
10 |
316 |
610,2 |
945 |
||
0,00033 |
6 |
7,5 |
9 |
10,5 |
12 |
572 |
1138 |
1823,4 |
||
0,00016 |
8 |
10 |
12 |
14 |
16 |
1014,9 2097,3 |
3469,3 |
|||
0,000083 |
10 |
12 |
14 |
16 |
18 |
1685,8 |
3617,1 |
6174,7 |
4.4. ОПТИМАЛЬНОЕ ОБСЛУЖИВАНИЕ СИСТЕМЫ ПО НАРАБОТКЕ И СОСТОЯНИЮ СОСТАВЛЯЮЩИХ ЭЛЕМЕНТОВ
Опишем сначала кратко модель, подробно изложен ную ів работе [6]. Рассмотрим, пренебрегая простоями системы, ее непрерывную работу.
143
Обозначим через v(t) параметр, характеризующим состояние системы в процессе работы:
,,, _ I 0, если система исправна в момент /;
' |
\ 1, в противном случае. |
к произвольному |
||
Пусть |
0 (t) — наработка системы |
|||
моменту t, тогда состояние системы в момент |
t может |
|||
быть описано вектором {ѵ(/)> |
0(01Считаем, что в об |
|||
ласти А возможных состояний |
систем |
известна |
подоб |
ласть Л0: при {v(t), О(0}<=Л0 система считается отка завшей. Пусть далее в момент ^о = 0 начинается эксплуа
тация новой системы. |
|
2, ... ) |
|
Предполагается, |
что в моменты /і = Ш (і= 1, |
||
наблюдается состояние системы и определяется |
вектор |
||
{''(/;), 0 (fi) >- Если |
в момент U (/ = 1, 2, ...) |
обнаружи |
|
вается, что система |
({ѵ (/,-), 0 (h) }е /4 0) отказала, тоона |
||
подвергается аварийно-профилактическому |
ремонту. |
В противном случае либо принимается решение о прове дении предупредительной профилактики системы, либо о продолжении ее эксплуатации до следующего момен та контроля ti+l = ti+At и т. д. Считаем, что после про филактики или ремонтов система является работоспо собной и ее наработка равна нулю.
Если принять, что аварийно-профилактический ре монт системы занимает постоянное время уа.п с коэффи циентом пропорциональности П О С Т О И М О С Т И Са.тт, то С Т О И М О С Т Ь его проведения
^-а.п Сап\ ап■
Аналогично для стоимости сп.п проведения преду предительной профилактики
^П.Н~ Сі.пУ,,.,,-
Считаем, что са.п>сп.п- Стоимость проверки работо способной системы (контроля ее состояния) полагаем равной нулю. Моменты окончания профилактики и ре монтов являются моментами регенерации (обновления) процесса, описывающего состояние системы. Интервал времени между соседними моментами регенерации назо вем периодом регенерации.
Искомая стратегия обслуживания в период регенера ции полностью определяется разбиением области А воз можных состояний системы {v(4), tK} на три непересекающиеся подобласти Аи А2, А3 {AYUA2UA3=A), при по падании в которые вектора {v(tK), tK} принимается
144
решение о продолжении наблюдений, о проведении пре дупредительного профилактического мероприятия или аівариино-профилактическопо ремонта соответственно
(Л3=Л 0). Наша задача заключается в оптимальном раз биении областей А х и Л2.
Введем критерий оптимальности разбиения области Л. Пусть произведено некоторое разбиение Аи Л2, Л3 и на рассматриваемом периоде регенерации время пребы
вания системы в области Лх равно t. Тогда |
удельные |
||
затраты за период регенерации равны: |
|
|
|
I г |
момент t |
осуществ- |
|
если в |
|||
лен переход |
в область |
Л2; |
|
с (Лх, Л3, Л3)—ct — |
|
|
|
у у , если в момент t осуществ |
|||
лен переход |
в область Л3 = |
Л0. |
Время t пребывания системы в области Лі зависит от разбиения Лі, Л2, Л3. За критерий, характеризующий качество разбиения Ль Л2, Л3, а значит, и качество стра тегии обслуживания в работе [6] были приняты средние удельные затраты, приходящиеся на единицу времени безотказной работы системы за период регенерации:
Мсі = с(Аі, Л2, Л3) .
Задача теперь свелась к определению таких обла стей Лі*, Л2*, Л3*=Л 0(Лі*и Л2* U Л3*=Л ), для которых достигается минимум средних удельных затрат:
с (Лі, |
Л2, Л3) = |
шіп с (Лх, Л2, Ла). |
(4.41) |
|
|
ЛI,Аг, Лз |
|
|
л,ил2ил3=л |
|
|
Строго говоря, |
критерием |
(4.41) можно пользоваться |
при мно |
гократном применении формулируемого ниже оптимального правила
разбиения |
А (, А г, Л3 в |
процессе длительной |
эксплуатации |
системы |
только в |
случае, когда |
возможные изменения |
величины t |
невелики. |
В остальных случаях управление не будет близким к оптимальному, хотя и может быть использовано на практике. Строгая оптимальность правила соблюдается при эксплуатации систем до первого попадания в области Л2 или Л3.
Покажем, когда можно предложенным правилом оптимального разбиения пользоваться на практике для систем с длительным пери
одом |
эксплуатации. Пусть Т — период |
эксплуатации системы, £ — |
случайная длина периода регенерации |
(Л1|=т), ср — математическое |
|
6 |
3353 |
145 |
|
ожидание затрат х\ на периоде регенерации (Мт| = ф), D — множест во стратегии разбиения d. Тогда суммарные затраты за время Г
щт)
;=-і
где .V (Г) — число периодов регенерации до момента Т.
н
мощью теории выпуклых фигур (применяя и рандомизированные привила разбиения) можно показать, что при некотором Q
Если |
Р {д0 < $ < а і \ |
= |
1 и Р |
{*о < т) < |
Ь\\ = 1, гдеЛ 0 >-<р — |
— Дц; |
< ер + Лц; д0 > |
т — Д£; а\ |
< т + Д£ , |
то для параметра Q , |
|
обеспечивающего искомый |
min |
как |
а |
||
rain |
можно |
написать, что |
|
Последнее выражение позволяет обосновать возможность пере-
т. е. использования полу-
ченных выше результатов для систем длительной эксплуатации.
Для определения оптимального разбиения Л Д Л2*, Л3*, обеспечивающего минимум средним удельным за тратам с(Ль Л2, Л3), воспользуемся леммой Дуба (см. п. 1.6). В работе [6] эта лемма применительно к рассмат риваемому случаю сформулирована следующим обра зом.
Пусть {Fs}, s = 1, 2, ... — последовательность а-полей в некотором пространстве W элементарных событий, FsciFs+l, а {г/Д — последовательность случайных вели чин, измеримых относительно Fa, и для которых сущест вует М{уі]. Пусть также задан некоторый класс К пра вил определения моментов k остановки наблюдений та ких, что ЛВД<оо.
146
Если существует правило, определяющее момент ос тановки k*, для которого:
1) Лр*]<оо;
2)
( > Уз-1. 5 > k*
для почти всех элементарных событий a ^ W , где s — моменты остановок наблюдений, определяемые другими правилами из класса К, и если существует константа С<оо такая, что для всех s
3) л Ч[ ^ +1 - ^ ] 1 ^ } < с , то данное правиле является оптимальным, т. е.
М [г/А*]==шіпУИ [ys],
где минимум берется по всем возможным правилам ос тановки из класса К.
В работе [6] было показано, что условия «1)» и «3)» леммы на практике всегда выполняются. Относительно выполнения условия «2)» необходимо сказать следую щее. Для реальной технической системы, не подвергаю щейся обслуживанию, естественным является предло жение о том, что с увеличением наработки надежность системы уменьшается. В этом случае вероятность
|
я | № ) . |
**]€ А,| |
является |
монотонно возрастающей функцией времени |
|
tK = kht, |
что необходимо |
и достаточно для выполнения |
условия «2)», которое в рассматриваемом случае запи сывается так:
< |
с/к_ {, если Р {[ѵ (/й), |
/ft)] € А3 |
||
|
при [ѵ(/,■), |
* = 0 . £ — ! ] } < |
||
|
^ ______^п.п_____ |
• |
||
I |
( с а.п |
^и.п) |
1) |
1 c s j |
М [С( IFt I І
кк I> cti_ v если Р{ѵ(Я>
при [ѵ(/,), і = 0 , k — l ] l >
>
Cll.ll_______
( с а.п Cn.J(Ä-l)
6* |
147 |
|
В данном случае последовательность {Fts} есть по следовательность наблюдаемых состояний системы
=t-t, i'=--0, s),
амомент оптимальной остановки [см. выражение (4.42)] определяется с помощью неравенства
р I[V(#*), Q |
е л , I[V{ t , \t h i = о ; л- i ] } < |
|
< |
__________ gn.n__________ |
(4.43) |
|
( с а.и ^п.ц) |
|
Правая часть выражения (4.43) монотонно |
убывает, |
а левая возрастает при возрастании /г. Поэтому, если в качестве /г* взять такое k, для которого последний раз выполняется неравенство (4.43), то условие 2 будет выполняться.
Таким образом, все условия леммы Дуба примени тельно к данному случаю оказываются выполненными, и оптимальное правило остановки (нахождения оптималь ной границы между областями А{ и Л2) определяется с
помощью выражения |
(4.43), которое в случае оптималь |
|||
ного обслуживания |
системы |
по наработке |
принимает |
|
вид: |
|
|
|
|
|
'к - 1 + Д' |
|
|
|
- |
) |
МОЮ |
K(t)dt |
|
1 — |
|
|
т |
< |
к—1 |
|
|||
|
1' |
Л(ОМ |
|
|
е |
о |
|
|
|
|
< |
^ |
---------- , |
(4. 44) |
|
U a , .. — С п .п)(А — 1) |
’ |
где Л(і) — функция интенсивности отказов системы, оп ределяемая выражением (4.1).
Правило построения области At таково: всякая точна {v(4_j, 4)}, для которой выполняется (4.44), принадле жит области Л).
Обслуживание системы по состоянию составляющих ее элементов
П а р а л л е л ь н о - п о с л е д о в а т е л ь н о е с о е д и н е и и е э л е м е н т о в
На практике во многих случаях соединение элемен тов системы по надежности можно представить как ком
148
бинированное соединение параллельных и последова тельных элементов.
Рассмотрим параллельно-последовательное соедине ние элементов системы по надежности, представленное на рис. 4.8. Пусть число последовательно соединенных разных по надежности элементов системы равно N, а число работающих параллельно и также разных по на дежности элементов — п. Отказ системы (предполага ем, что она работает непрерывно) наступает при отказе
любого из N элементов и при отказе более т (0< /п < я) из п параллельных элементов (резерв нагруженный). Считаем, что подключение резерва автоматическое и происходит достаточно быстро, т. е. временем его под ключения можно пренебречь. Надежность автоматиче ского коммутатора вновь не учитываем. Считаем также, что распределение времени безотказной работы каждого элемента последовательной группы подчиняется некото рому закону Fj(t), у = 1,. . N, а каждого элемента па раллельной группы — закону ФД/), t —L
Для периодического наблюдения за состоянием эле ментов системы назначается некоторый интервал вре мени At и в моменты U = iAt (/ = 0, 1, 2,...) контролиру ется состояние всех элементов системы. Стоимость конт-. роля, как и прежде, учитывать не будем.
Если в момент контроля обнаруживается отказ од ного из N последовательно соединенных элементов или отказ (т+1) - го из параллельных элементов (в любой последовательности), то система подвергается планово му аварийно-профилактическому ремонту. Во время проведения этого ремонта неотказавшие элементы си стемы проверяются в таком объеме (а отказавшие за меняются или капитально ремонтируются и устраняются все последствия имевших место отказов), что после ре-