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

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

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

точная или упрощенная. При этом в области низких частот рассматриваемые

характеристики могут заметно

отличаться

примере

на

рис. 12.30

при

О)<0,2 1/с).

как влияет

это

различие

на

переходные

про­

Для того.чтобы выяснить,

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

таких характеристик. Они

 

отличаются

только

в

полосе низких

частот

0^co<Qi.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для оценки расхождения переходных функций можно использовать метод,

предложенный в [23]. Составим разность

 

 

 

 

 

 

 

 

где Ро(со)— вещественная

 

р (“ ) =

■ £*(“ ) — Л ) (“ )■

 

 

 

 

строенная без

упрощений,

частотная характеристика замкнутой системы, по­

 

а

Р (м ) — упрощенная

характеристика. Очевидно,

что

р(со) = 0

при со^О ь

 

 

соответствующие

указанным характеристикам,

 

Переходные

функции,

 

 

будут отличаться друг от друга на величину

 

 

 

 

 

 

 

 

 

 

 

 

 

Si

 

sin

t

 

 

 

 

 

 

 

в (t) =

 

2 р

 

da.

 

 

 

 

 

 

-1

р (<•>)------------

 

 

 

 

 

 

 

 

 

It J

 

 

 

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

Учитывая, что время работы системы стабилизации ограничено некоторым

значением і Р, а

наибольшее

 

расхождение переходных

функций

должно

быть

при

больших значениях времени, оценим

e(t)

 

в момент времени

t = tp:

 

 

 

 

 

 

 

 

2,

sin

atp

 

 

 

 

 

 

 

 

 

 

2

г*

tpdu>-

 

 

 

 

 

£( ^ ) = —

\ P(“) —

 

:—

 

 

 

 

 

 

 

 

я

j

 

<1dp

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

Если выполняется условие

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ір

я ,

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

'

а і

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 р

('

 

 

 

 

p F

 

 

 

 

 

 

I £ Ѵр ) \ < — Г - \ Р ( “ ) da = ■

 

 

 

 

 

 

 

 

 

 

я

о

 

 

 

 

я

 

 

 

где

F — площадь, заштрихованная на рис. 12.32.

 

 

 

 

 

Для приближенной оценки можно принять, что

 

 

 

Я (0 )-Р 0 (0)

2і.

поэтому

I * ( * , ) !

я

I Я ( 0 ) - Я о ( 0 ) I •

 

 

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

Qitp 1z3 (0) — Po (0) 1

I Еотн ißp) I

я Я„(0)

582

при

я

Фактическая погрешность, как следует из вывода, будет меньше.

Если относительная погрешность, вычисленная по этой формуле, оказы­ вается в допустимых пределах, то при анализе свойств рассматриваемой си­ стемы можно не учитывать большие постоянные времени в передаточной функции летательного аппарата, которым соответствуют сопрягающие частоты (Ü< Q I, так как учет этих постоянных времени окажет несущественное влияние на переходную функцию. С помощью аналогичных приемов можно достаточ­ но просто обосновать возможность тех или иных упрощений рассматриваемых передаточных функций в любом конкретном случае. Пример оценки влияния вариаций скорости на процессы стабилизации ракеты можно найти в [14].

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

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

ГЛАВА XIII

БОКОВОЕ ВОЗМУЩЕННОЕ ДВИЖЕНИЕ

§ 1. ХАРАКТЕР СВОБОДНОГО ВОЗМУЩЕННОГО ДВИЖЕНИЯ

Выясним характер свободного возмущенного движения. Для этого, положив в уравнениях (11.48) бн= б э= 0 и Z'&— MXB = MV&= = 0, рассмотрим полученную систему однородных линейных диф­ ференциальных уравнений.

При изучении свободного бокового движения систему (11.48) можно упростить, понизив порядок системы до четырех.

Сначала исключим неизвестную гТ с помощью соотношения

cos 0-'F = cos Ѳ-ф — ß-j-ay.

Продифференцируем это равенство

и отбрдсим

произведения

малых величин

dQ

1ТР

de

,

da

 

 

 

----- dt

Фу

------dt

%

■— у. Тогда, используя четвертое

 

 

 

 

dt

 

 

 

 

уравнение системы

(11.48), получим

 

 

 

 

о äW

 

о

di/

йф

I

d~{

cos Ѳ

<u„ —

d$

, d~(

cos Ѳ -----

= cos Ѳ —i-------

—4- a —- =

---------

dt

-4- a ——.

dt

 

 

dt

dt

 

dt

cos ft y

dt

Следовательно, третье уравнение системы (11.48) можно перепи­ сать так:

Ьцшу4" ~ ~ +(^42 — ^42) ß a

-j- bi6y —0,

где

cos Ѳ

cos &

При изучении свободного движения уравнение

d'b

1

(13.1)

—- =

------о)„

dt

cos ft у

 

584

можно отбросить, так как в остальные уравнения приращение ф не входит *. После решения основных уравнений легко найти приращение ф:

Ф—------ ( “Ѵ W

(13.2)

cos а ,і

а

 

о

 

 

Таким образом, будем исследовать следующую систему урав­ нений:

+ с\і®*+ сп°Ѵ + сі2?= 0 ;

Ьіошх-\-~7~

bnw bl2^-\- bn —^- = 0;

 

 

dt

dt

 

 

 

(13.3)

^41ЮІГ~Ь ~ 7' + (^42~ ^42) ß— а —у-+^4бУ = 0;

 

я

dt

dt

wx — tg &• (O---- = 0.

x

s

y

dt

Пусть невозмущенное движение представляет собой прямо­ линейный установившийся полет. Тогда коэффициенты уравнений (13.3) будут постоянными.

Как и в случае продольного возмущенного движения, частные решения однородной системы (13.3) ищем в виде

wx= A e lt, <Oy=Belt, ß = C ew, y — Dlt,

(13.4)

Тде А, В, С, D — постоянные.

Подставляя эти частные решения в однородную систему (13.3) и производя сокращение на множитель еХг, получаем систему че­ тырех линейных однородных уравнений относительно неизвест­ ных А, В, С и D. Эта система имеет ненулевое решение, если

Х-фсц

 

С12

'

0

 

*10

X-f Йц

b\ik-\-bn

 

0

(13.5)

0

 

"X"I- *42 — Ьіч

— aX-)-è4S

*41

 

1

— tg»

 

0

 

— X

 

Раскрыв этот определитель, получим характеристическое

уравнение системы (13.3)

в

виде алгебраического

уравнения

4-й степени:

 

 

 

 

 

 

 

*4+ Л хз+ М

2+ / ^ + Л

= 0-

(13.6)

* Это уравнение используется при исследовании системы стабилизации угла рыскания ф.

20— 8422

585

Коэффициенты ри Р2, Рз, Рі этого уравнения зависят от коэффи­

циентов уравнений

(13.3), т. е. от производных

коэффициентов

аэродинамических

/3

3

3

ei)

to

си.,

сил и моментов (с$,

тх, ту, тхх, тух, тху5

туу, т1)>

от конструктивных параметров летательного

аппа­

рата (G, S,

I, Іх, Іу)

и от режима установившегося полета

(F*, Я*,

Ѳ*, Ф», а*).

Корни характеристического уравнения определяют характер свободного возмущенного движения летательного аппарата. Если корень К будет вещественным, то соответствующее ему частное решение будет определять апериодическое движение, затухаю­ щее или возрастающее в зависимости от знака К. В случае комп­ лексного корня Хі = х + іѵ характеристическое уравнение будет иметь и сопряженный ему корень Кг== >с—іѵ. Паре комплексных сопряженных корней соответствует колебательное движение с пе­ риодом Г= 2я/ѵ и коэффициентом затухания х. В зависимости от знака х амплитуды колебаний с течением времени возрастают или убывают.

При исследовании бокового возмущенного движения лета­ тельного аппарата наиболее часто встречается случай, когда ха­ рактеристическое уравнение (13.6) имеет одну пару комплексных сопряженных корней х±іѵ и два вещественных корня. При этом общее решение уравнений (13.3) имеет вид

шх= Л1еХі<4~ ^ 2еХа< +

^ /е’І<sin K + 'h);

Ш у = В хеХі< -f-Я2е;'2* +

B ' e xt sin (^ + ф2);

ß =

C1eV + C2ex'^ -fC /ex< sin (ѵ/+фз);

у =

D 1eXi< -f-

-j- D'ext sin (vtf-j- ф4)

и свободное боковое возмущенное движение складывается из

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

Обычно из двух вещественных корней один по абсолютной

величине получается очень

большим, а второй — очень малым,

комплексные же корни имеют промежуточные значения.

1 гл. XII.

Для примера рассмотрим тот же самолет, что и в §

При полете этого самолета

на

высоте 12 км

со скоростью Ѵ=

— 222 м/с

(М = 0,75) коэффициенты уравнений

(13.3) равны:

Яю = 0,0198

1/ с;

Ьц = 0,19

1/с;

Ьи' = 0; Ъі2 = 2,28

1/с2;

ö4i = —1;

*42 = 0,059 1/с; *46 = 0,0442 1/с; *42/ = 0; с„ = 1,66

1/с;

с„' = 0,56 1/с;

£12 = 6,2 1/с2. Здесь для простоты принималось ft = 0 и сс = 0.

Характеристическое уравнение в рассматриваемом

примере

имеет следующий вид:

 

 

 

 

 

 

Я +

1,909X3_|_2,69Х2-f 3,95Х- 0,00437 =

0.

 

586

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

\ = — 1,695; Х2 = 0,001105; Я3і4 = -0 ,1 0 7 ± 1,525/.

Как видно, свободное боковое движение самолета характеризует­ ся двумя вещественными корнями (Яі — большой корень харак­ теристического уравнения, Яг — малый корень) и парой комплекс­ ных корней Я3 и Я4.

Решение уравнений возмущенного движения для начальных

условий (при / —0) cö;c = ft)2/ = ß = Y = 0 и бн=5,73° имеет следую­ щий вид:

со^ = — 0,0085с-1•695< -f 0,0881 е°-001105< -f

Ц-0,1026е-°'1075< cos (87,5/ + 140,8°);

ші/= 3,502-0,00015е-1’695і-3 ,5 1 0 е0'00110«4-

+ 0,0364e-°-107« cos (87,5/ ■+ 92,5°);

ß = —0,3261+ 0,2923е0’00110^ + 0,0353е-°-1075<cos (8 7 ,5 /- 3,1°);

y = — 79,95 + 0,0050e-1'695< + 79,85e°-0011051 +

+ 0,0668e-°’10751 cos (87,5/+ 47,5°).

На рис. 13.1 приведены графики результатов расчета. Расче­ ты показывают, что свободное боковое движение летательного аппарата самолетной схемы (с плоскими крыльями и обычной

о

1

1 3

4

5

6

7

8

9

10

11

12 13

14

15 16

77 18

і,с

Рис. 13.1. Пример бокового возмущенного движения Само­ лета

20*

587

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

Частное движение, соответствующее большому вещественно­ му корню, называют движением крена. Это движение апериоди­ ческое, оно продолжается обычно не больше секунды (в рассмат­ риваемом примере время его затухания вдвое составляет 0,41 с). Обычно в выражениях (13.7) произвольные постоянные В і и Сі во много раз меньше, чем другие постоянные. Это говорит о том, что рассматриваемое апериодическое движение в основном за­ ключается в изменении угла крена и угловой скорости со*.

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

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

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

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

Характер свободного бокового движения на втором и третьем этапах определяется соотношением между коэффициентами Ь\г

и Си-

Если коэффициент Ьи достаточно велик по сравнению с си, колебательное движение быстро затухает (по истечении несколь­ ких секунд) и на третьем этапе остается только спиральное дви­ жение с параметрами со*, ѵ, ß, у либо медленно затухающими, либо медленно возрастающими (в случае спиральной неустойчи­ вости) .

Если же коэффициент ЬІ2 имеет сравнительно малые значе­ ния или отрицательные, то на втором этапе затухает спиральное движение, а колебательное движение продолжается на третьем этапе.

588

§ 2. ПРИБЛИЖЕННЫЕ УРАВНЕНИЯ И ПЕРЕДАТОЧНЫЕ ФУНКЦИИ ДЛЯ БОКОВОГО ВОЗМУЩЕННОГО ДВИЖЕНИЯ ДИНАМИЧЕСКИ ОСЕСИММЕТРИЧНОГО

ЛЕТАТЕЛЬНОГО АППАРАТА

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

Для упрощения кинематических и геометрических соотноше­ ний в уравнениях возмущенного движения летательного аппара­ та удобно использовать наклонную земную систему осей коорди­ нат OxByaza, которая получается из обычной земной системы по­ воротом последней относительно оси Oz3 на некоторый угол Ѳн. Если при этом направление земной оси Охн будет близко к на­ правлению полета, то можно считать малыми не только углы Ч'1 и ф,но и углы Ѳ и -O'*.

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

Если измерять углы ф, #, у, Ч*, Ѳ, ус в наклонной системе ко­ ординат, когда траектория полета близка к оси Охп этой систе­ мы, то существенно упрощаются кинематические и геометриче­ ские соотношения в системе уравнений бокового движения за счет замены косинусов углов § и Ѳ единицами, синусов этих же углов — самими углами и отбрасывания членов с произведения­ ми малых величин, например ау и tg^-coy.

Тогда последние три уравнения системы (11.48) примут прос­ той вид:

dfy

 

dt

 

dl

(13.8)

dt

 

W = ф — ß.

 

Заменив в первых двух уравнениях системы (11.48)

и щ

на у и ф и отбросив в третьем уравнении произведение sin 0 -ß, получим

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

589

J * L . , . ' u J S L + &

“*

■f" c\$ —

^із^э Сіз8н-}- с15М Хьф,

dt2

dt

 

dt

 

 

+ »I. —

+ «И - ^

+ M* + »Й —

dt*

1 11 dt

1 10

dt

1 12 1

dt

 

=

— ^13&H— ^13 [dt

(13.9)

 

 

c№

dt

* 4 2 P -M = M H+ *45^b;

 

ф — ЦГ — ß = 0 .

Дальнейшие упрощения обусловливаются конкретными усло­ виями рассматриваемой задачи — аэродинамической компонов­ кой аппарата и типом системы управления. Уравнения (13.9) ис­ пользуются обычно при исследовании стабилизации летательных аппаратов самолетной схемы. У таких аппаратов наблюдается сильное взаимное влияние движений рыскания (вращения отно­ сительно оси Оу\) и крена (вращения относительно Охі). Эти

движения взаимно связаны благодаря влиянию членов Сц'ір,

Ci2ß, сіз'бв, öioY, ^4бУ и, строго говоря, одно движение нельзя от­ делить от другого.

2.1. ДВИЖЕНИЕ РЫСКАНИЯ

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

Угловая скорость крена приводит к появлению момента рыс­ кания М “*(ох, которым в уравнениях бокового возмущенного

движения осесимметричных аппаратов можно пренебречь *. Кро­ ме того, влияние крена на рыскание проявляется в том, что при накренении летательного аппарата возникает боковая сила, вы­ зывающая скольжение (рис. 13.2). Если система стабилизации достаточно быстро ликвидирует крен, то влияние указанной бо­

ковой силы на движение рыскания мало и в первом из уравнений р а |_ у

(11.33) можно опустить член -— усЕсли принять указанные

допущения, то из системы уравнений бокового движения можно

* Влияние отклонения элеронов на момент рыскания в большинстве слу­

чаев мало даже для аппаратов самолетной схемы, т. е. Л1^Э8Э« 0 .

590

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

аГ/2 dt dt dt

+ bi s ^ . ;

(13.10)

dl”

b4$ b43^H4bv/-v'

d t

ф „ 1 г _ р = о

Часто опускают члены b[$ и бі38н, учитывающие запаздыва­ ние скоса потока (см. разд. 3.2 гл. XII). В результате получают более простую систему уравнений:

і - ьп ~ ~ -\-ьі $ ~ ~ ьпК-\-ьхъМув',

 

 

 

 

~~Т,-----Ь4$ = 6438н “Ь Ь45^в\

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

ф —

 

— р= 0.

 

 

 

 

Легко

проверить, что

уравнения

(13.10), характеризующие

движение рыскания, ничем не

 

 

отличаются по форме от урав­

 

 

нений

 

продольного

движения

 

 

(12.50).

 

 

 

 

 

 

 

 

 

 

 

 

Углу

■тангажа

O'

соответ­

 

 

ствует угол

рыскания

ф,

углу

 

 

наклона

траектории

Ѳ — угол

 

 

поворота

траектории

'Ф’,

углу

 

 

атаки

 

а — угол

 

скольжения ß,

 

 

углу

отклонения

органов уп­

 

 

равления

 

тангажом

б — угол

 

 

отклонения

 

органов

управле­

 

 

ния

рысканием

бн,

динамиче­

 

 

ским

 

коэффициентам

ощ, aJ2,

 

 

öi3, о42, CLkz— коэффициенты bn,

 

 

b b13,

 

 

6 4 3

 

 

 

 

можно

 

G

12,

Первые

 

уравнения

 

 

 

642,

 

 

.

 

 

 

 

 

 

 

получить

из

последних,

заме-

Рис.

13.2. Образование скольже-

нив коэффициенты йш коэффи-

ния

при крене летательного аппа-

циентами Ьш. При этом совпа-

 

рата

дает не только

форма

уравне­

 

 

ний, но и числовые значения коэффициентов, так как для дина­ мически осесимметричных летательных аппаратов оц = 6ц; аі2=

— Ь і 2', ci\2 — b i 2 \ йіз = ^із; a l 3' = b i 3' ; 015 = 615; 042= 642; 043= 643;

0 4 5 = 6 4 5 .

591