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

книги из ГПНТБ / Динамика полета и конструкция крылатых летательных аппаратов

..pdf
Скачиваний:
52
Добавлен:
30.10.2023
Размер:
21.98 Mб
Скачать

О т к л о н е н и я

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

х а р а к т е р и -

с т и к. При исследовании динамики летательных

аппаратов ис­

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

. фактических аэродинамических характеристик летательного аппа­ рата от аэродинамических характеристик модели объясняется ря­ дом причин.

Во-первых, в процессе аэродинамического эксперимента сами, аэродинамические характеристики модели определяются не вполне точно. При этом различные характеристики определяются с различ­ ной точностью. Если ошибка в определении коэффициентов аэроди­ намических силобычно не превосходит 5%, то ошибка при опре­ делении производных демпфирующих или спиральных моментов ■может достигать 20—40%. Если зависимости аэродинамических коэффициентов от различных факторов носят ярко выраженный не­ линейный характер, то погрешность при определении таких зависи­ мостей резко возрастает. Поэтому для анализа переходных процес­ сов необходимо знать не только аэродинамические характеристики, но и точность их определения.

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

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

рактеристик.

остаточные деформации

элементов

конструкции

В-четвертых,

приводят к тому,

что в процессе

эксплуатации аэродинамические

характеристики

летательного

аппарата

постепенно

изменяются.

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

§ 1. 7. О Т К Л О Н Е Н И Я Х А Р А К Т Е Р И С Т И К А Т М О С Ф Е Р Ы

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

сферы могут довольно сильно отклоняться от их расчетных значе­

ний, причем эти отклонения можно в общем случае считать случай­ ными.

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

Отклонения скорости звука зависят от высоты полета. До высот порядка 80 км эти отклонения обычно не превышают 20% от стан­ дартной величины.

§ 1.8. НАЧАЛЬНЫЕ ОТКЛОНЕНИЯ

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

•называют начальные значения независимых переменных и их про­ изводных.

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

Стартовые возмущения задаются в соответствии с видом стар­ товой установки и условиями старта.

Г Л А В А I I

ПЕРЕХОДНЫЕ ПРОЦЕССЫ И ПЕРЕДАТОЧНЫЕ ФУНКЦИИ ЛЕТАТЕЛЬНОГО АППАРАТА В ПРОДОЛЬНОМ ДВИЖЕНИИ

§ 2. 1. В О З М О Ж Н О С Т Ь Р А З Д Е Л Е Н И Я Л И Н Е А Р И З О В А Н Н О Й

У Р А В Н Е Н И И Д В И Ж Е Н И Я Н А У Р А В Н Е Н И Я П Р О Д О Л Ь Н О И Б О К О В О Г О Д В И Ж Е Н И Я

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

331

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

§ 2. 2 . Л И Н Е А Р И З О В А Н Н Ы Е У Р А В Н Е Н И Я П Р О Д О Л Ь Н О Г В О З М У Щ Е Н Н О Г О Д В И Ж Е Н И Я

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

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

Для исследования устойчивости и характера переходных про­

цессов при продольном движении воспользуемся

общими урав­

нениями движения центра тяжести КЛА в путевой

системе ко­

ординат (см. (1.44), (1.45) и

(1.46)

второго

раздела] и уравне­

нием моментов в связанной

системе

координат

[см. (1.24) вто­

рого

раздела]. Полагая

в этих уравнениях у = 0,

р = 0 и W^ — O

(боковая составляющая

вектора скорости

ветра

отсутствует),

а также принимая (Jy — Jx) ^ x^y ~ 0 (по малости),

имеем

G dV^

Pcos(<pe + а„) + УS i n (а — а„) — Q cos (а — а„) — С? sin 0„,

g

dt

 

 

 

 

 

 

( 2. 1)

 

 

 

 

 

 

 

 

332

( 2.2)

(2.3)

(2.4)

где Mz — сумма продольных моментов, действующих на КЛА

впродольном движении,

Мtp — продольный момент от тяги силовой установки.

Силы и моменты, а также параметры движения КЛА можно рассматривать как суммы этих величин в исходном невозмущен­ ном движении и их приращения (возмущения), вызванными дей­ ствиями возмущающих причин, т. е.

.

Р = Р 0 +

ДР,

 

 

Y =

Y 0 - \ - A Y ,

 

 

Q =

Qo +

AQ.

 

Mz=

мг0 +

д м2,

 

M zp = -/Wz p o + A M zp ,

 

 

V n = V M + W „

(2.5)

 

в п —0лО+ Д0я>

 

 

“я = Кл0 +

Аап,

 

« — «л =

[(К)о — ал0] + (Да-

Да«)>

•») = 4о + Д?1

и т. д.

 

Здесь индекс „нулик“ показывает принадлежность соответст­ вующих сил, моментов и параметров движения к исходному ре­ жиму полета. Очевидно что производные параметров с индексом , „нулик" будут равны нулю.

Используя (2.5) и принимая косинусы приращений углов рав­ ными единице, а синусы равными самим приращениям, уравне­ ния (2.1) (2.2), (2.3) и (2.4) представим в следующем виде:

ззз

^ У г = ЛР C0S + Я»)°— Р0 Sin К * + О Ла„—

AQ cos (а — а„)0 + Q0 sin (а — ал)0 (Да — Да„) +

+А к Sin (а — а„)0 + У0 COS (а — а„) (Да — Да„) —

G cos 0ОД0,

- J ^ „ о

— ДР sin (аа, + а„)о +

Р 0 COS (ал + а„)вДал +

+

дQ sin (а — а„)0 +

Q0 COS (а — а„) (Да — Дап) +

 

 

 

 

( 2. 6)

+

Д Y COS (а — а„) — У0 sin (а — ая)0 (Да — Да„) —

G

 

cos0пОД \/e + G ( 1 -----Т7Г - ) sin бпОД0„,

 

 

л-бс. 1

'

кос. 1'

 

 

^Дш,

 

 

Л- ^ г = дл*4+Д;ИгЛ

хг- — sin 0ОД 1/„+ 1/0 cos 0лОДОл.

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

шг0 = 0, 4^0 = 0.

«*0 =

(«)о.

 

 

 

К о =

к .

 

 

 

®л0= ®0-

 

 

 

Поэтому далее принимается

 

 

 

 

Д1Г£ =

1FEj

 

 

 

ДИУ„ =

117,,

 

 

 

Дсог =

Шг.

 

 

 

Вначале найдем зависимость

между Д1/„

и Д1/, Wv

V,, 60

и

между Да„ и Да, Wv U7,, 0О, VQ-

 

 

 

 

Пользуясь формулой (1.40) второго раздела, полагая

U7t =

0,

или непосредственно из фиг. 2.1

получим

,

 

 

V„= V cos р cos (а — ал) + U7, sin 0Л+

cos 0Л.

(2.8)

334

Подставляя значения параметров в виде сумм V„ — 1/л0 + АУ^я»

а = (а)о *Ь А® и т. д. в (2.8) с учетом (2.7), отбрасывая члены

вто­

рого и выше порядка малости, имеем:

 

A V „ ^ A l/+ ^ c o s 0 o+ tt^sin 0О.

(2.9)

При малых 0Оможно считать:

 

W^sinOo^O,

 

Д1/„ x& .V + W v

Далее, как известно [см. (1.42) второго раздела или фиг. 2.1],

ал = а — arc tg Х р - + arc tg

.

Vni

vt

Откуда, переходя к отклонениям, получаем.

Дая = Да — Д^агс tg -уЗ. j + Д ^arctg -^ -j ,

335

ИЛИ

 

 

 

1

Д Vni\ Ун*

Д

VnTfj

 

 

 

Да,, -•= Да —

 

+

 

 

 

n’lO

 

Vko

 

 

 

 

1 +

 

 

 

 

 

 

V.ntО

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

*V4 ^ - Д

^

1/т,о.

 

 

 

 

 

т]0 1

VI,

 

 

 

 

 

 

1 -+-

 

 

 

 

 

 

£0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если

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

то

 

 

 

 

 

.

.

 

 

 

'

bV,: Va - i V

t V,f

( 2. 10)

Да. = Да-------------- гтт;---------------- [- ----

 

--------

 

 

 

 

•I/и0

 

 

 

Vo2

 

 

Так как

в расчетном

движении

 

 

 

 

 

 

 

 

 

 

V„o= V0,

 

 

 

 

 

 

 

 

 

v„eo= v £0,

 

 

 

 

 

 

 

 

 

Vпт$“ ^т,0 j

 

 

 

 

 

то выражение

(2.10)

принимает вид

 

 

 

 

 

' Да„=Д« +

^ ( - Д У ячУа + ДУя£1/,,о + A^ K

0- A

^ I / t,o) .

 

Откуда,

подставляя

 

Д 1/^= Д К , + WT),

 

 

 

 

 

 

 

 

 

 

 

 

получаем

 

 

+ Wlt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Д а^Д а + ^ ^ ^ ^ о - Г г У а , ) .

 

 

 

Поскольку

 

 

^0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- ^ =

sin0o,

 

 

 

 

 

 

 

 

 

^ о

 

 

 

 

 

 

 

 

 

 

тг =

cos б0,

 

 

 

 

 

 

 

 

 

V о

 

 

 

 

 

 

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

имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

(2. 11)

 

 

Да„= Да + -T7-sin 0 eW,i:--- гу-cosO0\V,

 

 

 

 

 

V0

 

v л

^0 vv V

 

 

Зная

зависимость

Да„ от Да, U7S и

Wq,

легко установить по­

добную связь

между .ДО- и Да,

W W q

и т. д.

 

 

 

336

Так как при продольном движении З' = '0л4- ап, то

 

Д&= ДОя + Дв||.

(2. 12)

Тогда, подставляя (2.11) в (2.12),

находим

 

 

 

 

 

(2.13)

При малых 60, когда

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

sin 0О1#^ можно считать

величиной второго порядка малости, a

v Оп

 

cos60~ l ,

 

 

 

О

 

Теперь, пользуясь полученным соотношением (2.11), из си­

стемы (2.6) исключим Да

и Дкл .

При

этом будем полагать, что

угол (аа„ + а)0 мал. Тогда,

пренебрегая

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

sin(ade+

+ а)0Да„ и APsin (адв + я)0, как

величинами второго

порядка

малости, и считая cos (з.да -|- а) ^ 1 , можно привести систему уравнений (2.6) к следующему виду:

. (2.15)

Входящие в уравнения (2.15)

отклонения сил п

моментов от

их расчетных значений являются

функциями многих

независимых

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

ной силы Д Y зависят от воздушной скорости, числа М, высоты полета, угла атаки и угла отклонения руля высоты. Отклонение

22 А . Г . Б е д у н к о в н ч и д р .

337

продольного момента ДМх зависит от воздушной скорости, числа М, высоты полета, угла атаки, угловой скорости <ог , скорости изме­ нения угла атаки и от угла отклонения органа продольного управления 80. Отклонение момента силы тяги &MzP является

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

дР = Р ''д У + Р5'Лт1+ р !;', Д8р + ДР,

 

A Q ---- Qo Л V +

QoДа +

Q dA-q + Qo"Д °в +

ДQ0

 

 

AY=YoAV+ Уо Да + Уо Дт)+ УЗ"Д8в +

ЛУ„

 

 

A M , = M ~ qA V + Ж “о Да + М ~ о Дт) + М ™ог <», +

 

 

 

 

■dAa

■ ■

(2Л6^

 

+ МЪа- £ + Мя*АЬа + АМ1в

 

 

A M z P = M ^ А V + М 1 ро Да + М ] ро Дт, +

 

 

 

+

ж

;>

г +

ж ^

д в р + д мгРв

 

 

где

символы вида Р о ,

M z

и

т. д. означают частную

производ­

ную

соответственно

 

 

> ( “5 ^ ) и т - Д-»

>AQ0>ДУв, bMze,

ДМ,Рв — слагаемые

первого

порядка малости, учитывающие влия­

ние возмущений,

не

вошедших

в остальные члены

разложения

(отклонения параметров атмосферы, отклонения характеристик летательного аппарата и др.).

При исследовании переходных процессов в продольном дви­ жении за независимые переменные удобно выбрать путевую

скорость V„,

угол наклона траектории 0„, угол

тангажа & и вы­

соту полета

т|. Составляющие скорости ветра

и Wn являются

возмущениями. Если рассматривается летательный аппарат без системы управления, то управляющие воздействия 8р и 80 также

следует считать возмущениями.

Воспользовавшись соотношениями (2.9) и (2.13), из (2.16)

исключим A V и Да. Тогда имеем

 

АР=Р%А Va-+РЗ Дт)+ Рор Д8я + Ро cos 0ОЩ - Р о sin 001 +

ДРв

AQ= QoA V „ - Ql Дб„+ Qg Д» + Q? Дц+ Q0S"ДВ,■+ (■-Q<Tsin 60

—Qo y - Sin e0j + ( Qo sin 0O+ Qo - y cos 0O\ Wn +

AQe (2.17)

AY—Yo AVn—УоД0„+ Уо*Д9 + У^Д7) + Уз®Д8в+ (— У^cos0o— .

- У0“ ^ sin0o^H75 + / —Ko^sin 0О+ У0* cos 0оj W, + ДУ.

338

Д/И2= Л ^Д Vn- M o ~

- Mlo^n + (m %+ M*o)

+ Mzo A9 + '

+

Mo&y H- Л У Д8,— Mo у

sin 0o —jjj- +

(—vWj£cos0o—

-

Mlo -^-sin0o'j U7f + Mlо-1

cos 60^ +

(■- M

l sin0O+

 

+ ^ o -^ -c o sO ^ lF , + AM2e

 

(2.17)

т г р = < poА ц . - м ; * , Д0л + У И > 2+ ж ;ро д 0 + ж ^ 0 Ar, +

 

+ М ^ М р+ ^ _

Mvzpacos е0 - м :ро1 - sin 0oj

W,+

+ ^— М],р0 sin 0o -f- M гРд -у- cos 60^ W-^ -f- АМгРй

Так как уравнения проекций сил и моментов (2.15) позволяют определить только четыре независимые переменные: А1/л , Д0П, Дт), шг , то необходимо добавить к системе (2.15) уравнение, свя­ зывающее Д0 с остальными неизвестными, входящими в систему. Таким уравнением является кинематическое соотношение

 

dt

'Z

 

(2-18)

 

 

 

 

Подставляя в уравнения (2.15)

выражения (2.17) и обозначая

d

 

 

-ij£ через оператор

р, линеаризованную систему уравнении про­

дольного движения с учетом (2.18) запишем в общем виде:

рА V/, о.'уA Vn+ а^Д0 + а^у АН + а'J, Atj+ а у Аор-f-

+ а у AS,+ a y ЕIFS+

a v 11117, + f v

р Д0 = a]f A V„ + a{jД0 4- аЦД&'+ао‘Дт| + а^Д80 +

 

+ a ^ W , + a ^ W , + / e

ри>г — a v A Vn-f- аа_ рА6 + аш^АО -j- amz шг -)- ат А& + . (2.19)

 

 

О

,

А С*

+ а»

Atj + d t f А8р +

 

*

аюв До, + а-™*1р И7£+

 

 

>

 

 

 

 

7

 

 

+

1 ^ + а<> pWv +

 

W4 + / щ

рД0— (1>г

 

 

 

 

рАч= а^ А^ л+ а®Д0

 

 

 

22*

 

 

 

339

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