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

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

.pdf
Скачиваний:
730
Добавлен:
27.10.2023
Размер:
28.56 Mб
Скачать

Таким образом, квадратичный трехчлен (12.84) в знаменате­ лях передаточных функций (12.80) — (12.82) можно разложить на элементарные сомножители, как показано в табл. 12.3.

 

 

 

 

 

Таблица 12.3

а \ 2 + 011^42

 

а р %+ ft/? + 1

Корни характеристического уравне­

6

 

ния

(12. 85)

 

< 1

П р 2’ + 2 <і Т р + 1

,2 = у - ( — ' ± і Ѵ \ — £2|

> 0

 

р + 1)2

PU2 =

 

1

1

— у

 

 

 

 

 

 

 

 

( Т ' р + 1) {Т " р + 1)

 

1

_

1

 

> 1

р 1

pr ' р%

рп

0

Р (Тр + 1)

Р\ =

0; р2=

- у

 

 

 

 

 

 

 

( Т ' р + \ ) { Т " р + \ )

р1

1 _

1

< 0

Т' < 0

•pt ЛР%

рП

 

 

 

 

 

 

Разложение на множители числителей передаточных функций

Анализ проведем для устойчивых летательных аппаратов (öi2+ önfl42> 0 ), разделяя известные аэродинамические схемы на две группы:

1) летательные аппараты с положительным передаточным ко­

эффициентом по углу атаки:

) > 0 (схемы «утка» и с по-

воротными крыльями);

/уст

 

2) летательные аппараты с отрицательным передаточным ко­

эффициентом по углу атаки: (— ) < 0

(обычная схема и

\ о /у с т

 

«бесхвостка»).

Схемы «утка» и с поворотными крыльями. У статически устой­ чивых летательных аппаратов этих схем передаточный коэффи­ циент К всегда положителен, т. е.

«12«43—«13042>°-

Числитель передаточной функции (12.80), учитывая, что Ö13 <0, можно представить в виде

(7 > + 1 ) (7 > + 1 ),

(12.94)

542

где

T ^ B + V Л + ß 2, Т2 = В — У А + В \

л ____ ^13____ ,2 5 = — «13 —■^13^42 + Öi2a 43

(12.94а)

 

 

а 12^43 — а 13#42

д 12а 43 — 013042

 

 

Числитель передаточной функции (12.81), принимая

во

вни­

мание, что й4з> 0, представим в виде

 

 

 

 

 

Т % р 1-\-2%вТ вр~\~ 1;

(12.95)

Т в =

____ £43____

 

 

 

а 12^43 — «13а 42

 

 

 

 

 

(12.95а)

 

аца43 + й12043 — ^13042

 

 

 

5 0 = --------—

 

 

 

 

 

2 У « 4 3 ( а 1 2 а 4 3 —

0 1 3 а 4 2 І

 

 

Разложим на множители числитель

передаточной

функции

(12.82), учитывая, что а4з+ аіз'>0.

Предполагая, что

аіз +

+ йца4з < 0, представим числитель этой передаточной функции в виде

ТеР У 1>

(12.9,6)

где

 

;r6 = _ f f i ±

_ f i £ (12.96а)

а 13 + 0 -1 1 а 4 3

 

Обычная схема и «бесхвостка». Передаточный коэффициент К у статически устойчивых летательных аппаратов этих схем отри­ цателен, т. е.-

аПа4Ъ

013^42

0.

У летательных аппаратов

обычной

схемы коэффициент Сіз'

практически равен нулю, и числитель

передаточной функции

(12.80) запишется в виде

 

 

Т іР+ 1,

(12.97)

где

 

 

— ß 13 + а 12а 43

(12.97а)

«12^43 — «13д 42

 

 

 

В рассматриваемом случае корни квадратичного трехчлена в числителе передаточной функции (12.81) вещественны и его мож­

но записать в виде

 

(7 > + 1 )(7 > 4 -1 ),

(12.98)

543

где

T3 = D + y C + D2,

Т4= Z>— У С D2,

 

с=- — «43

2D:

а 43 ( а 1 1 + ^ 1 2 )

(12.98а)

 

а 12а 43 — а 13а 42

 

а 12а 43 — а 13а 42

 

Так как у аппаратов обычной схемы аіз + 0паіз>0, то в числи­ теле передаточной функции (12_82) получим

где

_____Д-43___

#13 + #11#43

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

7 > 2 + 2£7>+1,

имея в виду, что коэффициент g может быть меньше или больше единицы или равен ей, как показано в табл. 12.3.

Таблица 12.4

Схемы «утка“ и с поворотными

Схемы обычная и „бесхвостка“

крыльями

< (р)

j ,

(Т\Р +

1) (ХчР +

1)

к

Ті Р + 1

 

^

П р 2 + 2\Тр + 1

 

1

 

 

' ч

72/72

+ 2ЧТр +

К ( р )

г,

 

+ 2£ѳТѳр +

1

г

(ТзР +

1) (Т4р +

1)

К

Т2 р 1

+ 2 ЬТр+ 1

 

Т2р2 + 2£7> +

1

 

 

 

 

 

Щ( р )

Kz

 

Т*Р+ 1

 

 

 

Т5р + 1

 

 

 

 

 

 

 

 

T2p2 + 2 Z T p + l

 

T2p2 + 2 t T p +

1

 

 

 

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

Напомним, что полученные передаточные функции справедли­ вы при следующих условиях:

1) скорость полета известна, т. е. рассматривается первый этап возмущенного движения;

2) влияние силы тяжести на возмущенное движение невели­ ко (а44» 0 ) ;

544

В пределах этих ограничений полученные передаточные функ­ ции справедливы для любых летательных аппаратов (обычная схема и «бесхвостка» как с аэродинамическим управлением, так и с управлением поворотом газовой струи, схемы «утка» и с пово­ ротными крыльями). Однако значения параметров передаточных функций (К, Т, %и др.) у аппаратов различных схем могут су­ щественно отличаться, что сказывается на динамических свойст­ вах летательного аппарата.

Маневренные летательные аппараты с малой нормальной силой, создаваемой органами управления

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

мы и схемы «утка», нормальная

сила органов управления Y 8

мала и ее можно не учитывать

при составлении передаточных

функций летательных аппаратов. Такое допущение, очевидно, нельзя принимать для летательных аппаратов схем «бесхвостка» я с поворотными крыльями и для бескрылых аппаратов типа бал­ листических ракет. В этих случаях обязательно приходится учи­

тывать нормальную силу органов управления Е88.

Если нормальная сила органов управления пренебрежимо мала, то в уравнениях возмущенного движения (12.49) или в вы­ ражениях передаточных функций (12.80) — (12.82) можно поло­ жить

(12.99)

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

Принимая во внимание (12.99), получим следующие переда­ точные функции:

(12.100)

( 12. 101)

(12,102)

545

Здесь

К--

012Т"013д 42

 

 

0ЦД42

 

Т! =

1 .

 

042

 

 

 

/ о \ =

— а13 .

(12.103)

«12 + аПаА2

 

7 =

1

 

^ а\2 + апаі2

 

е =

а П + а 12 +

а 42

 

2 1/ а12+ Лца42

Параметры передаточных функций (12.100) — (12.102) зави­ сят от конструктивных и аэродинамических параметров летатель­ ного аппарата, от скорости и высоты полета. Чтобы выяснить влияние этих факторов, упростим выражения (12.103), полагая, что летательный аппарат имеет достаточно большую (по -срав­ нению с демпфированием) степень статической устойчивости. Тогда, пренебрегая величиной ацО-ю по сравнению с ai2, получим:

 

013042 =

mz

P + Y*

'

(12.104)

 

012

 

/и“

гпѴ

 

 

 

 

 

К

У-

ml

Р + у а

 

(12.105)

 

 

G

 

 

g

<

 

 

 

 

 

 

 

 

 

кт, ^

__ 0 1 3 ___

 

m z

 

(12.106)

 

 

а 12

 

к

 

 

 

 

 

 

 

 

Сравнивая (12.105) с (8.44) и (12.106) с (5.92), замечаем, что

 

 

57,3ц■у бал 5

 

 

 

кт,

о /б а л

 

 

 

 

 

 

 

Нетрудно видеть,

что

передаточные коэффициенты К и

„ V

К — , характеризующие маневренные свойства аппарата, в силь- g

ной степени зависят от скоростного напора и степени статиче­ ской устойчивости (главным образом от центровки аппарата).

На больших высотах передаточные коэффициенты К я К g

546

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

Передаточный коэффициент по углу атаки КТ\ не зависит от скоростного напора, но весьма существенно зависит от степени статической устойчивости.

Взависимости от скорости и высоты полета коэффициенты К

иК — могут сильно изменяться вдоль траектории. Так, на­

пример, у зенитного управляемого снаряда «Эрликон» коэффици­

ент К изменяется более чем в 20 раз, а К — — более чем в

&

10 раз (рис. 12.7).

к,1/с;к-у/д

Рис. 12.7. Изменение передаточ­ ных коэффициентов снаряда «Эрликон» в течение полета

0 10 20 30 40 50 t,C

Для упрощения требований, предъявляемых к системе управ­ ления в различных условиях полета, желательно, чтобы коэффи­

циент К или К — (в зависимости оттипа системы управле-

g

ния) изменялся бы в возможно меньшем диапазоне.

Допустимые диапазоны изменения передаточных коэффициентов К и

V

К — могут быть установленъ! лишь в процессе разработки системы

g

управления или на основании опыта проектирования однотипных систем. Иногда на начальных этапах проектирования можно использовать некоторые

общие рассуждения. Пусть, например, требуется, чтобы величина

К

менялась возможно меньше. Рассмотрим упрощенное выражение этого пере­ даточного коэффициента, полученное из (12.105) отбрасыванием тяги Р:

К — ж - 57,3с“

qS

57,3ml

qS

G

X j'~ ‘~*xp

^

g

у

547

Очевидно, что величина

Ш,

должна изменяться возможно

 

G

меньше. Поскольку отношение q\G уже определено выбором траектории поле­ та, желательно так подобрать параметры летательного аппарата, чтобы

Щ

~ ----- ~— изменялось примерно пропорционально G/q. Этого можно ДОСТИГ­ ЛАj

нуть:

1)соответствующим перемещением центра масс по мере выгорания топ­ лива (на активном участке полета) или

2)перемещением крыльев вдоль продольной оси летательного аппарата.

Рассмотрим постоянную времени Т\. Пользуясь соотношения­ ми (12.100)—'(12.102), получим передаточные функции, в кото­ рых выходной величиной является угол Ѳ, а входными величина­ ми — углы а и D:

Ѳ(р) =

1

.

(12.107a)

а(р)

Т іР

 

 

 

в (/>) _

1

 

(12.1076)

» (Р)

Т іР +

1

 

Из (12.107a) видно, что величина, обратная постоянной вре­ мени Т1, представляет собой отношение угловой скорости каса­ тельной к траектории к углу атаки

1 _ ѳ

Т і

а

Чем больше Ть тем больший угол атаки требуется для создания нужной угловой скорости Ѳ.

Рис. 12.8. Переходный процесс при ступенчатом отклонении продольной оси летательного аппарата

Из передаточной функции (12.1076) следует, что в результа­ те скачкообразного поворота летательного аппарата на угол AD направление вектора скорости изменяется на такой же угол A0 = AD, но лишь после переходного процесса, описываемого уравнением

t

Д9-

548

и изображенного на рис. 12.8. Чем больше Ть тем больше време­ ни требуется для изменения направления вектора скорости.

В зависимости от скорости и высоты полета постоянная вре­ мени Т1 может значительно изменяться:

j, __

mV

__

1

2m

1

57 ,3c*qS + Р

5 7 ,3с* + ср

?SV '

В качестве примера на рис. 12.6 приведено изменение Тх вдоль траектории снаряда «Эрликон».

Летательный аппарат с поворотными крыльями

Для упрощения анализа рассмотрим случай, когда можно

пренебречь влиянием запаздывания скоса потока (ап ~ а із~ 0 ). Тогда на основании формул (12.94) — (12.96а) получим, что

— «із

«12«43 — «13«42

7’2= 0 ;

*43

Ts

Y «12«43— «13«42

«11«43

$Ѳ:

2 У «43 («12«43 — «13«42)

-Г'

«43

7 6--------;

«Ц «43

«13 +

 

Реакция летательного аппарата на отклонение органов уп­ равления будет описываться передаточными функциями [см. табл. 12.4], справедливыми также и для схемы «утка»:

&(/>)■_._ is

Т\Р + 1

(12.108а)

Ч р )

р(Т2р2 + 2 і Т р +

1) ’

ѳ (Р) _ ts

Т\р2 + 2$07 > +

1

р (Т2р2 4 . 2£7> +

(12.1086)

Н р )

1)

а (Р) _

тъР + <

(12.108B)

в (р)

а T2 p2 + 2 t T p + \

Нетрудно видеть, что передаточные функции (12.108) и (12.100) — (12.102) аппаратов с поворотными и закрепленными крыльями отличаются друг от друга коэффициентом оиз-

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

549

мулам гл. Ill и V) имеет место следующее соотношение:

« із

т*

(12.109а)

т,

«12

— —— «0,1 ч-0,2.

 

 

Откуда получаем, что в среднем у аппаратов с поворотными крыльями

„5

343

: 0,3 4-0,4;

«42

043

013

я 42

«12

Из (12.109а) следует, что

«12343 — ß i3ß42 = “ З12З42.

(12.1096)

При приближенном анализе естественно также пренебречь произведения­ ми аца42 и апа4з по сравнению с а12 и аіз (соответственно). В результате, используя (12.1096), получим упрощенные выражения параметров передаточ­ ных функций аппаратов с поворотными крыльями:

^f l 3 _

 

012

 

 

2

013 .

Т1 1,

---Ä У;

 

 

а 42 012

( 12. 110)

■ л / 2 - ^ - Г ;

У«42

 

З ц

1 /

343

 

 

К

«42 ’

 

й4з

 

 

 

013

 

 

Хотя передаточная

функция (12.108а) аппарата с поворотными крыльями сов­

падает по форме с

аналогичной передаточной функцией (12.100) аппарата

с закрепленными крыльями, значения параметров К и Т, этих передаточных функций существенно различаются. Для примера представим, что аппараты с поворотными и закрепленными крыльями имеют одинаковые динамические коэффициенты ап, а 12 и «42- Тогда параметры Т и | у этих аппаратов будут равны, но

Кпов.кр'

 

Кз.кр

 

 

 

2 «13 \

 

 

 

 

 

 

 

 

«12 /з-кр

 

Іпов.кр'

2

а 13

I

13.кр-

 

1

1

 

 

0 1 2

/пов.кр

 

.550

4.3. ПЕРЕХОДНЫЕ ПРОЦЕССЫ ПРИ СТУПЕНЧАТОМ ОТКЛОНЕНИИ ОРГАНОВ УПРАВЛЕНИЯ

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

Представим установившееся криволинейное движение лета­

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

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

аппарата параметры #, Ѳ, пу и а изменяются в течение некоторо­ го промежутка времени. Этот процесс называется переходным. По окончании переходного процесса устанавливаются новые зна­

чения параметров “О, Ѳ, пу и а, соответствующие новому положе­ нию органов управления.

Изменение угла атаки, перегрузки и угловой скорости касательной к траектории

(летательный аппарат с закрепленными крыльями)

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

Выражения (12.101) — (12.102) показывают, что по отношению к выходным величинам Ѳ, а, пу летательный аппарат является звеном второго порядка, передаточную функцию которого можно записать в виде

Z L E L = --------- а----------,

(12.111)

»(/»)

7 Ѵ + 2£7>+ 1

 

где X— любая из величин: Ѳ, а, пу;

К — соответствующий передаточный коэффициент:

К, K — , К Т Ѵ g

551