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

книги из ГПНТБ / Коваленко И.Н. Полумарковские модели в задачах проектирования систем управления летательными аппаратами

.pdf
Скачиваний:
9
Добавлен:
23.10.2023
Размер:
7.56 Mб
Скачать

В процессе работы ЗЭ восстановлений или замен отказавших элементов не производится, а осуществляет­ ся лишь автоматическое подключение резёрвных эле­ ментов вместо отказавших основных. Зарезервирован­ ный элемент считается отказавшим при выходе из строя одного из основных элементов и при отсутствии исправ­ ных резервных элементов.

Задача

заключается в

определении оптимального

периода Тк

профилактического

обслуживания ЗЭ, при

назначении

которого

обеспечивается

максимальная ве­

роятность застать ЗЭ

в исправном-состоянии в про­

извольный момент времени 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) - го из параллельных элементов (в любой последовательности), то система подвергается планово­ му аварийно-профилактическому ремонту. Во время проведения этого ремонта неотказавшие элементы си­ стемы проверяются в таком объеме (а отказавшие за­ меняются или капитально ремонтируются и устраняются все последствия имевших место отказов), что после ре-

Соседние файлы в папке книги из ГПНТБ