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

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

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

3.2. УПРАВЛЕНИЕ ЗАПРАВКОЙ САМОЛЕТОВ В ВОЗДУХЕ

Вопросам увеличения дальности полета самолета и способности при минимальном запасе топлива доставлять грузы на предельные расстояния по-прежнему уделяется большое вшім.ание. Одним из эффективных средств по­ вышения дальности полета и полетной нагрузки считают системы заправки самолетов в воздухе (см., например,

147]).

В настоящее время системами заправки оснащаются не только бомбардировщики и военно-транспортные са­ молеты, но и истребители. В качестве заправщиков ис­ пользуются как специализированные, так и боевые само­ леты, снабженные подвесным заправочным оборудова­ нием. Специализированные заправщики применяются для дозаправки самолетов ВВС Англии, США, Франции и других стран.

Создание систем, обеспечивающих заправку в полете, идет несколькими путями и диктуется соображениями, связанными с выбором оптимальных параметров для различных целевых назначений.

Основные тенденции в развитии систем заправки та­ ковы:

— расширение диапазона скоростей и высот за­ правки;

—■сокращение времени обнаружения, стыковки и пе­ релива;

снижение веса системы;

обеспечение автономности энергопитания;

снижение нагрузок на топливоприемники заправ­ ляемого самолета;

унификация заправочного оборудования.

Ведутся работы по автоматизации контактирования; при этом используется опыт наведения ракет на цель и стыковки космических аппаратов на орбите. Сколько-ни­ будь детальное рассмотрение применяемых за рубежом систем заправки и организации их эксплуатации выхо­ дит за рамки этой главы. Заинтересованного читателя отсылаем к работе [47].

Если пренебречь инфинитезимальным поведением системы заправки и учитывать лишь временные и отчасти качественные характеристики процесса заправки, то та­ кая система с позиций теории массового обслуживания может быть рассмотрена как СМО смешанного типа.

120

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

Пример 3.3. В районе дозаправки постоянно дежурит п самоле- тов-дозаправщмков. Если дозаправка началась, то опа осуществля­ ется до конца и длится в среднем 10 мин. Если все дозаправивши заняты, то самолет нуждающийся в дозаправке, некоторое время мо­

жет

ожидать

(совершать

полет по кругу в районе

- дозаправки).

Среднее время

ожидания — 20 мин. Если

самолет так

и не дожида­

ется

дозаправки, он садится

на запасной

аэродром.

В среднем в

район дозаправки за час прибывает 24 самолета. Число самолетов, ожидающих дозаправки, ничем не ограничено.

Определить оптимальное число самолетов-дозаправщиков,

при

котором минимизируется платежная функция

 

С (л) = (цХ^ож 4- с2К „ п 4- сцКзП,

(3. 8)

где С(п) — суммарные потерн в системе дозаправки горючим

в воз­

духе за единицу времени

(1 ч) ее работы;

в воз­

Сі —-ущерб от ожидания одним

самолетом дозаправки

духе за единицу времени; с9 — ущерб от простоя одного дозаправщика за единицу вре­

мени; с:і — затраты на эксплуатацию одного дозаправщика за едини­

цу времени.

 

Принимаем: с, = 50 единиц стоимости;

ся='-25 единиц стоимости;

г3=25 единиц стоимости.

 

Решение. По формулам (1. 13) — (1.19)

с помощью табл. 5 и 6,

приведенных в книге [35]. определяем характеристики гипотетичес­ кой СМО.

1. Вычисляем параметры а и ß:

10

а = 24— = 4 ,

60

_Ю_б0

= 0,5.

6020

2.На основе накопленного опыта в организации систем доза­ правки ориентировочно определяем область допустимых значений

числа каналов обслуживания. Пусть л=1, 2, 3, 4, 5. Для этих значе­

ний я

определяем

характеристики СМО и

результаты

сводим в

табл.

3.1.

значения платежной функции для я —1,...,5 и сво­

3.

Определяем

дим результаты в табл. 3. 2.

следует, что

т іп С (я )

Из приведенных значений функции С(л)

имеет место при я= 4 , т. е. искомое n 0pt= 4.

 

 

На основании изложенного можно отметить, что не­ которые задачи организации и управления воздушным

121

 

 

 

 

 

 

 

 

 

Таблица

3. 1

п X

а

Р

Ро

Л'з

 

Р о5с

Р

 

V 4

*ож’ 4

к »

V

отк

 

1 24

4

0,5

0,0027 0,9973 5,57

0,248 0,71

0,9973

0,00045

0,232

0

2 24

4 0,5

0,0073

1,956

4,05

0,49

0,51

0,977

0,00383

0,168

0

3 24

4

0,5

0,0121

2,8

2,4

0,7

0,3

0,927

0,0113

0,1

0,073

4 24

4

0,5

0,0151

3,3

1,17

0,83

0,17

0,845

0,0305

0,05

0,155

5 24

4 0,5 0,0171

3,7

0,64 0,92

0,08 0,737

0,0593

0,03

0,263

 

 

 

 

Таблица 3.2

П

с' 1'ож

С»Ѵ '

csKsn

C(n)

 

1

278,4

0,075

24,9

303,3

2

201,5

1 .1

49

251,6

3

120

5

70

195

4

60

17,5

82,5

160

5

36

42,5

92,5

171

движением приводят к необходимости изучения явлений ожидания и допускают формулировку на языке теории массового обслуживания.

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

Глава IV

ПЛАНИРОВАНИЕ ОБСЛУЖИВАНИЯ СИСТЕМ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ

4.1. ПЛАНИРОВАНИЕ ЗАМЕНЫ ОТДЕЛЬНЫХ ЭЛЕМЕНТОВ СИСТЕМЫ ПРИ МГНОВЕННОЙ ИНДИКАЦИИ ОТКАЗОВ

В процессе работы САУ возникает необходимость в периодической замене отдельных наиболее важных и наименее надежных элементов.

Рассмотрим элемент с функцией распределения вре­ мени безотказной работы F(t) =Р(Х<1) и монотонно возрастающей интенсивностью отказов (F'(t) = f (t) ).

(4.1)

'1 — F (0

Здесь X — случайная величина времени от момента возвращения элемента ів исходное состояние до отказа. Для элемента назначается профилактическое обслужива­ ние (контроль) через случайное время Y, распределен­ ное по закону

G ( t ) ^ P ( Y < t ) .

I

Выполняется профилактическое обслуживание через min (X, Y) единиц времени. Другими словами, если пос­ ле предыдущего возвращения в исходное начальное сос­ тояние (момент (2 на рис. 4. 1) элемент безотказно рабо­ тал время У (в этом случае X>Y) , то профилактическое обслуживание проводится сразу же по истечении времени Y, причем время проведения профилактики предполага­ ется случайным. Обозначим его символом V; функция распределения величины V Fv{t) произвольна, а матема­ тическое ожидание '.М{У]= Гп<оо.

Если же после предыдущего возвращения в исходное состояние (после моментов времени /= 0 и t=U) элемент отказал до истечения срока, равного У=г/і и У=£/з и т. д.

123

(см. рис. 4. 1), т. е.

при X<Y,

то непосредственно после

отказа происходит

его замена.

Время,

затрачиваемое

на замену, предполагается случайным

и обозначается

величиной U, имеющей произвольную функцию распре­ деления Fc(t) и математическое ожидание M[U] = T3< oo .

Будем считать, что Ти< Т я.

о

Рис. 4. 1. К принципу осуществления замен через слу­ чайное время

Задача заключается в отыскании такого периода профилактического обслуживания, при котором элемент обладает максимальной вероятностью нахождения в ис­ правном состоянии в произвольно выбранный момент t (максимальная готовность) и максимальной вероят­ ностью исправного состояния на отрезке времени (, / + х (максимальная надежность за время работы в воздухе х). Обозначим р(х, і) вероятность того, что элемент окажется исправным в произвольный момент t и безот­ казно проработает в воздухе в течение времени і + х.

Время между очередными моментами возвращения элемента в исходное состояние обозначим через Z. Тогда из рис. 4. 1 ясно, что

(4.2)

Последовательные реализации случайной величины Z обозначим ZA.

Величины Zh образуют последовательность независи­ мых, одинаково распределенных случайных величин, т. е. образуют процесс восстановления с ц= M[Z] и функцией восстановления H(t) =M[Nt\, где Nt — максимальное зна­ чение п, при котором где

Sn —Z1-rZ2+ ... + Zn.

Если случайные величины X и Y обладают плот­ ностью, то dH(() можно рассматривать, как вероятность

124

того, что возвращение элемента в исходное состояние (восстановление) завершится в момент (, так как вероят­ ность такого восстановления на малом отрезке At равна H ( t +A i ) H(t)+0(At). Математическое ожидание чис­ ла восстановлений элемента на малом отрезке времени At равно

 

M[ANt]=pi(At)+o(Al),

(4.3)

где

рі(АІ) — вероятность того, что на отрезке времени

At

завершится одно восстановление.

р(х, t), на

 

Для вывода формулы, определяющей

рис. 4.2 изображена поясняющая диаграмма. Событие,

вероятность

которого

равна р(х, I), может произойти

двумя способами, а имен- -

{

*

 

но: первый

способ — до

ч

>—-

момента t элемент не под-

0

4

-------

вергается профилактике-

 

 

£+*

скому обслуживанию

с

Рис. 4.Ö. К выводу формулы

(4. 4)

вероятностью

1 — G(t)

И

не отказывает до момента

 

второй

способ — в мо­

t + x с вероятностью 1—F(t+x)\

мент t—£ происходит восстановление элемента с вероят­ ностью dH(tI), и далее в течение времени g элемент не подвергается профилактическому обслуживанию с ве­ роятностью 1 —G(%) и не отказывает за время 1+х с ве­ роятностью 1—F(E,+x). Так как £ меняется в интервале от 0 до t, то каждому значению g в данном интервале соответствует определенное значение величин H(t—g),

G(£) и F(£, + x),

поэтому

искомая вероятность

второго

составляющего

 

события

находится интегрированием

по всем %от 0 до

t.

 

 

Таким образом,

 

 

р(х,

0 = [ l - G l O ] [ l — F (*+*)] +

 

t

 

 

l - ^ + JO] d H (*-*).

 

+ J

[1 —0(E)] [

(4.4)

Отметим, что в течение фиксированного отрезка вре­ мени -V (работа ів воздухе) отказ элемента возможен,

а профилактическое обслуживание исключено, что учте­ но в выражении (4. 4). •

125

П р о и з в е д е м

з а м е н у

п е р е м е н н о г о в о в т о р о м с л а г а е ­

м о м п р а в о й ч а с т и ф о р

м у л ы

(4 .4 ): £ = (—о н Т о г д а

р ( х ,

i?) = [-l-G (if)l

[ 1 - / = ■ ( / + jc)] +

 

 

 

+ f [1-0(7-«)]

[1_ ^ (/ + л--со)]г//-/(щ ).

 

(4.5)

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С о г л а с н о

у з л о в о й

т е о р е м е

в о с с т а н о в л е н и я

( с м . г л . I)

с у щ е с т в у е т п р е д е л ь н о е

( п р и

/ — с о )

з н а ч е н и е

в е р о я т н о с ­

т и р ( х ,

t),

о п р

е

д е

л я

е м о е

ф о р м у

л

о й

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

оо

 

 

 

 

 

 

 

 

lim

р

 

(Л',

t) — р

(л-) =

Г[1— 0

(«)] [ 1— F ( ш -)- х)] сіло. (4. 6)

/+«

 

 

 

 

 

 

 

 

(А J

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

О

 

 

 

 

 

 

 

 

О п р е д е л и м

 

в е л и ч и н у

р .

П р и м е н я я

д л я

н а х о ж д е н и я

и с к о м о г о

у с л о в н о г о

м а т е м а т и ч е с к о г о

о ж и д а н и я

р

ф о р ­

м у л у

 

п о л н о й

в е р о я т н о с т и и

у ч и т ы в а я ,

ч т о с р е д н е е

 

в р е м я

р а б о т ы

э л е м е н т а

б е з

о т к а з а

и

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

о б с л у ­

ж и в а н и я

р а в н о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j' [1 — 0

(/)] [ 1— Т 7 (/)] dl,

 

 

 

(4.7)

 

 

 

 

 

 

О

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

оо

 

 

 

 

 

 

 

 

 

 

и м е е м

 

 

 

 

p =

j

U “

*2(01

[1 — Z7 (/)]

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

СО

 

 

 

 

 

 

СО

 

 

 

 

 

 

 

 

 

+ П

j

F

( X ) Ü Q (.X )+ Та j О

И

d F (Л-).

 

(4.8)

 

 

 

 

 

 

о

 

 

 

 

 

 

о

 

 

 

 

 

В

ы

р

а ж

е н и

е

(4 .6)

с

у ч е т о м

 

з а в и с и м о с т и

(4.8)

п р и ­

м е т

в

и

д

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

р(х) =

 

 

 

 

 

 

 

 

 

 

 

 

 

оо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

]■ [1 — о

(и)] [1— С

(w -j- х) ] rfco

 

 

 

,_____________ о________________

 

 

 

 

 

 

о о

 

 

 

 

 

 

 

 

 

о о

 

 

 

 

с о

 

 

*

f [ 1 - е (/)] [ 1 - е ( О] di+ т 3f г

(X) da (X) +

т „ f а

(х) dF (х)

О

 

 

 

 

 

 

 

 

 

0

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4 .9)

12Ö

Ітобы перейти к определению оптимального периода

контроля (профилактического обслуживания), опреде­ лим функцию G(t) следующим образом:

П.

если

(4.

10)

(О,

если

t < 7 'к, I

 

т. е. предположим, что контроль назначается через посто­ янное время Ти после включения элемента в работу.

Обоснование того, что наилучшее распределение

G(t) является вырожденным, дано в работе [6] и п. 4.5. Тогда

 

 

f [I — F (со + х ) ] d a

 

р (ЛГ) = р(х,

Тк)= ^ ---------

5--------------------------

---------•

 

Гh[l -

F (01 d t + T 3F (Т к) + Т ,,

f d F (X )

 

0

 

К

 

 

 

(4.11)

При этом

мы использовали очевидное равенство

 

\ F( x) dG( x) = F ( T K).

(4.12)

 

т

 

 

 

к

 

 

Вычислим производную выражения (4.11) по Тк, введя обозначения:

[К[ 1 - ^ И 1 ^ = А ; F(TK)= F; f ( T K) = f , (4.13)

о

использовав формулу

F (іо - j - х) Ä ; F ( U>) -f-/ (іо) X

(4.14)

и считая, что действительное значение х достаточно мало для возможности такой аппроксимации. Тогда

d p (X , Т к) __ ( 1 - F - f x ) ( А + T 3F + Т П~ T nF )

д Т к

(A + T aF + T „ - T „ F ) 2

 

( Л - Д л - ) ( і - д + r 3/ - r „ / )

 

{ А + T 9F + Т „ - T nF ) 2

127

После проведения элементарных алгебраических преобразовании, приравнивая производную нулю и пре­ небрегая членом Гп/л'-СІ, получим

_ 1и _ = 1 _____________ _____________:

(4.15)

Гз+ *

[1 - Я (Гк)]2

+ /

(Гк)

И] da

 

о [1 - Я

 

или

 

 

1

 

 

Т3-\-X

1 —

 

 

(4.16)

 

 

 

1 — Я (Гк) +

X(Гк) р [1 — Я (<і>)] сіш

 

о

Покажем, что при Тэ> Т п и монотонном возрастании K(t) функция р(х, Т,;) имеет единственный максимум. Для этого достаточно убедиться в том, что выражение

И г г ~ 1 + ------------------------

4 ----------------------

<4 - 17>

1 - f (Г,) + ).<ГК0 f [ l - f <»>]*»

О

при изменении Ти меняет знак плюс на минус, а это возможно лишь тогда, когда функция

Ф(7’к)= --------------------

—1 ------------------

(4.18)

1 — F (Гк) +

\ (Тк) (К[1 — Я (ш)] da

 

 

о

 

монотонно убывает с возрастанием Тк, т. е. когда ее зна­ менатель при этом монотонно возрастает. Следователь­ но, производная от знаменателя выражения (4. 18) дол­ жна быть положительна:

д

'__ 1_ 1

г"

(4.19)

дТк

Ф(Як)

] = ^ , (7’к)\П F (ш)] diо О 0.

 

 

 

В рассматриваемом случае условие (4. 19) всегда вы­

полняется,

так как

X(t)— монотонно возрастающая

функция и производная от нее положительна, а интеграл положителен потому, что интегрируется положительная величина.

При практическом использовании формул (4. 15) или (4. 16), в которых правая часть зависит только от вели­ чины Тк, строят график значений правой части в функции

т

Гк, а по оси ординат откладывают отношение ---- -— и

Тз+ X

128

Рис. 4. 3. К принципу осуществления замен через неслучайное
время

для каждого значения отношения графическим путем

определяют оптимальное время контроля (профилактического обслуживания).

После определения ве­ личины Hjopt организация эксплуатации элемента строится так, как пока­ зано на рис. 4.3, где светлыми кружками обо­ значены моменты выпол­ нения плановых замен, а

темными — неплановых, обусловленных отказами.

Пример 4. 1. Имеется непрерывно работающая система 1 с функ- иней распределенмн времени безотказной работы следующего вида:

F (t) = 1 — е

где 9=0,0005 1/м2.

Время выполнения плановых профилактических работ на систе­ ме постоянно и равно 15 ч (Гп—15 ч); среднее время устранения от­ каза в системе с последующим возвращением ее в исходное состоя­ ние равно 20 ч (Г3= 20' ч); оперативное время работы системы 10 ч (х=10 ч). Определить из условия обеспечения максимальной го­ товности системы к работе в произвольный момент времени t и мак­ симальной надежности системы за время от t до 7+10 ч интервал вы­ полняемых на ней плановых профилактических работ Т„, если эксплу­ атация системы организуется так, как показано на рис. 4.3.

Решение. Убедимся в том, что в нашем случае функция Я(/) мо­

нотонно возрастает во времени. іВ самом

деле,

из

выражения (4. 1)

следует, что K [t)= q t .

соотношением

(4. 16),

обозначив правую

Далее воспользуемся

часть этого соотношения через ф(Гк). Тогда

 

 

 

O.OOOSy

 

Г

0,0005/"-

е

2

+ 0,00057^ f e

 

2 dt

 

 

 

b

 

 

Расчетные данные функции <р(Т,г) для ряда значении ее аргумен­

та сведены

в табл.

4. 1 и иллюстрируются

графиком рис. 4. 4.

 

 

 

 

 

 

 

 

Таблица

4.1

тк, Ч

8

10

20

50

75

100

150

200

?(Тк)

0

0,03

0,09

0,36

0,57

0,65

0,78'

0,83

1 Здесь

система

рассм атривается

как элем ент.

 

 

 

129

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