О т к л о н е н и я |
а э р о д и н а м и ч е с к и х |
х а р а к т е р и - |
с т и к. При исследовании динамики летательных |
аппаратов ис |
пользуются аэродинамические характеристики, получаемые, как правило, в результате аэродинамического эксперимента. Отличие
. фактических аэродинамических характеристик летательного аппа рата от аэродинамических характеристик модели объясняется ря дом причин.
Во-первых, в процессе аэродинамического эксперимента сами, аэродинамические характеристики модели определяются не вполне точно. При этом различные характеристики определяются с различ ной точностью. Если ошибка в определении коэффициентов аэроди намических силобычно не превосходит 5%, то ошибка при опре делении производных демпфирующих или спиральных моментов ■может достигать 20—40%. Если зависимости аэродинамических коэффициентов от различных факторов носят ярко выраженный не линейный характер, то погрешность при определении таких зависи мостей резко возрастает. Поэтому для анализа переходных процес сов необходимо знать не только аэродинамические характеристики, но и точность их определения.
Во-вторых, за счет технологических допусков аэродинамические характеристики летательных аппаратов одного типа отличаются от их расчетных значений. Такие отклонения аэродинамических харак теристик также носят случайный характер (по ие меняются в про цессе полета). Зная допуски на изготовление летательного аппа рата, можно определить вероятные отклонения аэродинамических характеристик, вызванные этой причиной.
В-третьих, деформации летательного аппарата в процессе по лета приводят к изменениям аэродинамических характеристик. Эти изменения могут носить колебательный характер. Такие изменения аэродинамических характеристик особенно сильно влияют на устой чивость системы, особенно при наличии в ней управления нелиней ных элементов. Для точного учета влияния деформаций необходимо совместное решение уравнений движения летательного аппарата и уравнений, определяющих деформации элементов конструкции и вызванные этими деформациями изменения аэродинамических ха
рактеристик. |
остаточные деформации |
элементов |
конструкции |
В-четвертых, |
приводят к тому, |
что в процессе |
эксплуатации аэродинамические |
характеристики |
летательного |
аппарата |
постепенно |
изменяются. |
Этот фактор необходимо учитывать при анализе динамики летатель ных аппаратов многократного применения.
§ 1. 7. О Т К Л О Н Е Н И Я Х А Р А К Т Е Р И С Т И К А Т М О С Ф Е Р Ы
При анализе динамики летательных аппаратов пользуются так называемой стандартной атмосферой, отражающей усредненную за кономерность изменения основных параметров атмосферы с высо той. В реальных условиях фактические значения параметров атмо
сферы могут довольно сильно отклоняться от их расчетных значе
ний, причем эти отклонения можно в общем случае считать случай ными.
В уравнения переходных процессов из многочисленных парамет ров атмосферы входят только плотность воздуха р и скорость зву ка а, необходимая для определения числа М полета. Как показы вают метеорологические исследования, отклонения плотности воз духа зависят от времени года, а также изменяются в течение суток. Кроме того, изменение плотности может быть вызвано рядом дру гих причин. Градиент плотности по высоте тоже может меняться в широких пределах.
Отклонения скорости звука зависят от высоты полета. До высот порядка 80 км эти отклонения обычно не превышают 20% от стан дартной величины.
§ 1.8. НАЧАЛЬНЫЕ ОТКЛОНЕНИЯ
За начальные отклонения можно принять любые возмущения в начале переходною процесса. Обычно начальными отклонениями
•называют начальные значения независимых переменных и их про изводных.
Особый вид начальных возмущений составляют так называемые стартовые возмущения, цли отклонения условий старта от расчет ных. Типичными стартовыми возмущениями являются отклонения' координат места старта, отклонения углового положения летатель ного аппарата в момент старта, а также угловая скорость и угло вое ускорение стартового устройства в момент старта. Наиболее тя желые для переходных процессов комбинации стартовых возмуще ний могут возникнуть при старте летательного аппарата с подвиж ного стартового основания, расположенного, например, на корабле или на самолете.
Стартовые возмущения задаются в соответствии с видом стар товой установки и условиями старта.
Г Л А В А I I
ПЕРЕХОДНЫЕ ПРОЦЕССЫ И ПЕРЕДАТОЧНЫЕ ФУНКЦИИ ЛЕТАТЕЛЬНОГО АППАРАТА В ПРОДОЛЬНОМ ДВИЖЕНИИ
§ 2. 1. В О З М О Ж Н О С Т Ь Р А З Д Е Л Е Н И Я Л И Н Е А Р И З О В А Н Н О Й
У Р А В Н Е Н И И Д В И Ж Е Н И Я Н А У Р А В Н Е Н И Я П Р О Д О Л Ь Н О И Б О К О В О Г О Д В И Ж Е Н И Я
Одним из дополнительных упрощений, используемых при ана лизе динамики летательных аппаратов, снабженных системами уп равления, является разделение уравнений движения на уравнения продольного движения и уравнения бокового движения. При этом
под уравнениями продольного движения понимают уравнения, опи сывающие движение летательного аппарата, при котором его центр масс все время находится в вертикальной, точнее в ортодромической плоскости. Возможность такого разделения обеспечивается тем, что большинство летательных аппаратов (самолеты, ракеты и др.) можно считать аэродинамически симметричными относитель но ветрикальной плоскости. Это означает, что при отсутствии кре на, т. е. тогда, когда плоскость аэродинамической симметрии совпа дает с плоскостью движения, возмущения, входящие только в урав нения продольного движения, не могут вызвать изменений пара метров бокового движения. Аналогичное допущение для бокового движения было бы неверным, так как возмущения, входящие в уравнения бокового движения, в реальных условиях вызывают из менение параметров продольного движения, например . скорости. Тем не менее при качественном исследовании переходных процессов этим обстоятельствам часто пренебрегают для того, чтобы наиболее простым путем получить основные закономерности переходных про цессов в боковом движении.
§ 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) |
|
|
|
|
|
|
|
|
|
( 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) |
Подставляя значения параметров в виде сумм 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 ,
Так как при продольном движении З' = '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, + ДУ.
|
|
|
|
|
|
Д/И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 |