Переходя к принятой форме записи, главный определитель
Д (р) можно представить в виде
Д ( р ) ^ а , ( 7 > + l)(7Y/7* + 2i:,7> -|- 1), |
(2.65) |
где
( 2. 66)
Уравнения (2.59), (2.60) и (2.61) можно записать через пере даточные функции в виде
ДVn= W l ? ( p ) M p - r & 1 ' ( р ) \ К + К Ч р ) № : + & У Ц р ) \ У ъ, |
(2.67) |
Д0Л= |
^1р [р)Д8Р -г w;°(p)ak+ |
K |
h P) w, + |
w *'\P) w4, |
(2.68) |
Дт)= |
W * P ( р ) Дор + Wj (Р) AS.+ |
w ? |
[p)Wi + |
w p (P ) W ,,. |
(2.69) |
Раскрывая определители в правых частях уравнений (2.59), (2,60) и (2.61), можно получить в развернутой фс>Рме любую из интересующих нас передаточных функций. В качестве примера передаточные функции высоты по углу отлонения руля высоты и по вертикальному ветру приведены ниже:
й * ., „ч ___________К$Р + К в__________ |
(2.70) |
^ W - {TnP + \ ) ( T J p * + . 2 ^ T , p + \ ) |
’ |
где |
|
V ^ + < < 8) , |
|
п г -[л “( _ а као8+ at a v ) + a^ ( a v aoe - al av ) |
(2J1) |
и
(2.72)
Wri4p) (Тмр+1)(Тт* -I- 2 ^ 7 > + 1) ’
где
. (2.73)
§2.4. ПЕРЕХОДНЫЕ ПРОЦЕССЫ В ПРОДОЛЬНОМ ДВИЖЕНИИ
Вцелом переходный процесс в продольном движении, возникаю щий в результате действия возмущения, в общем случае склады вается, мак уже было показано, из колебаний летательного аппара та относительно центра тяжести и колебательного изменения высо ты и скорости .полета, -на которое накладывается медленное апери одическое изменение этих же величии. Частоты этих колебательных
движений .разнесены достаточно далеко. Такой характер переход ного процесса, несмотря на то, что точное заключение об устойчи вости исследуемого движения может быть сделано только на основе анализа полных уравнений движения, дает возможность судить об устойчивости движения и характеристиках качества переходного процесса непосредственно по выведенным в предыдущем параграфе упрощенным передаточным функциям с достаточной для инженер ных целей точностью. Кроме того,' имеется возможность также свя зать устойчивость и качество переходного процесса с параметрами движения и аэродинамическими и .конструктивными параметрами летательного аппарата.
Как было показано в предыдущем параграфе, передаточные функции, соответствующие короткопериодическому движению, имеют в знаменателе либо многочлен вида
F ( p ) = T vp(TSp* + TcaTap + \), |
(2.74) |
либо |
(2.75) |
F (р) = T'v p (TSp* + 2ija Тар + 1). |
Для устойчивости системы необходимо и достаточно, чтобы кор ни знаменателя передаточной функции имели отрицательные веще ственные части.
В рассматриваемом случае знаменатель передаточной функций имеет один нулевой корень. Однако характеристическое уравнение полной системы уравнений продольного движения (2.27) в общем случае .не имеет нулевых корней. . Следовательно, появление нуле вого корня является результатом введенных упрощений.
Практически характер короткопериодического движения опреде ляется корнями квадратного трехчлена, входящего в выражения (2.74) и (2.75). Для того чтобы эго движение не было расходя щимся, необходимо, чтобы эти корни имели отрицательные веще ственные части. Для этого необходимо н достаточно, чтобы коэф- ' фициенты квадратного трехчлена были 'положительными. Поэтому в рассматриваемом случае условия устойчивости принимают вид
|
TJ > |
0, |
(2.76) |
2 |
> |
0. |
|
Если Га9 < 0 и 2£аДа < 0, то имеет |
место |
неустойчивость, а прн |
Та2 = 0 — нейтральность. |
|
|
|
Из формул (2.33) видно, что неравенства (2.76) удовлетво
ряются при выполнении двух условий |
|
» , Ош. |
» 0 |
^ |
• Яш + OqОш- |
Йпао> |
0 |
(2.77)
■ а'“2 > о
Рассмотрим более подробно эти условия. Для этого, подста вляя полученные выше значения коэффициентов в (2.77), по следние условия представим в следующем виде:
|
|
|
c/qSba |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> 0 |
|
|
|
|
|
(2.78) |
|
|
|
VJ, |
|
|
|
|
|
|
|
где |
|
|
|
" ''> и |
|
|
|
|
|
|
|
- |
дсх |
|
|
|
|
gpbg |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
' 1 + |
^ |
|
1 |
|
|
|
|
|
|
|
|
|
|
cv“ |
)2 G S |
|
|
|
|
|
|
Vs Г |
Vs |
\ |
. й |
- |
|
|
|
+ (mfz-r-m |
v2 , |
sin 0 — х т, |
|
|
|
|
|
V2 |
V |
кос. |
1 |
|
|
|
|
|
® — |
1 — |
sin 0 |
|
|
|
|
|
|
|
|
|
V 2 |
, |
|
|
|
|
|
|
|
|
1 + |
|
|
кос. |
1 |
|
cv |
CyQSbJ*- - |
т*г - т\ . |
(2.79) |
|
|
|
|
|
Поскольку |
в |
(2.78) |
сомножители |
cvaqSba |
g s b 2 |
являются |
|
Л |
и |
VJZ |
•существенно |
положительными |
и |
коэффициент |
( всегда |
поло |
жителен, то условием сходимости короткопериодической |
соста |
вляющей продольного движения |
является |
|
|
|
|
|
|
|
|
о», > |
°- |
|
|
|
|
|
(2.80) |
|
|
|
|
|
|
|
|
|
|
Пользуясь выражением (2.79), можно записать
|
= С — Х т |
|
|
|
где |
|
|
|
|
|
— . — |
дсх |
. 1, |
®Cv |
gpbii^a. |
i |
c = x F+ у |
дву |
m " ‘z ( 1+^4 |
2(7/5 |
|
+ {тр + тп*) (1 - |
1/2 |
gba |
sin 6. |
|
Откуда видно, что o„ |
при заданном режиме полета определяется |
величиной продольной центровки х т. |
Очевидно, |
что если Злу = |
— с — хт > 0, имеет место устойчивость, при о„ = с — х т< 0 —
неустойчивость, а при °пу= с — /ст= 0 — нейтральность.
Величину л'т, при которой ап обращается в нуль, называют нейтральной центровкой х ГНевтр ■Ясно, что
Таким образом, величина нейтральной центровки главным образом определяется относительной координатой фокуса KJ1A.
Следовательно, коэффициент о„ представляет некоторую раз
ность между нейтральной центровкой и эксплуатационной цент ровкой КЛА. Поэтому обычно o„v называется запасом устойчи- •
вости КЛА. Так как положение’ центра тяжести меняется по мере выгорания топлива, а положение аэродинамического фо куса зависит от М, то запас устойчивости в процессе полета может изменяться в значительных пределах, что отражается на характере переходных процессов.
Что касается второго из условий (2.78), то единственным от рицательным слагаемым в левой части неравенства может быть
член |
М ,' |
V21/1 |
1 |
' sin 0, который при обычных условиях по |
|
\ |
кос. |
|
лета не может изменить знака суммы. Однако полет по наклон ной траектории по сравнению с горизонтальным полетом при прочих равных условиях по характеристике переходных про цессов несколько будет отличаться.
О качестве переходного процесса можно судить по величине вещественной и мнимой частей корней знаменателя передаточ ной функции, так как вещественная часть определяет затухание
колебаний, а мнимая часть — их частоту. Для принятой формы записи передаточных функций корни квадратного трехчлена, входящего в знаменатель, могут быть представлены в виде
|
Рi,2 |
Са , . У 1 —?а2 |
|
/Г) |
|
|
т |
i J |
-р |
• |
|
(2.81) |
|
|
1 а |
|
* а |
|
|
|
|
Подставляя |
значения |
и |
Та из |
(2.33) |
с |
учетом (2.79), |
по |
лучим |
L _ |
|
рVSba2 |
|
|
|
|
|
|
|
|
(2.82) |
|
Та — |
4Л |
V"y ’ |
|
|
|
|
|
|
|
/ 1 |
— |
|
|
|
|
pSba3vny |
(2.83) |
|
Та |
|
|
|
|
8суаЛо„ |
|
|
|
|
|
|
|
Как видно из формулы (2.82), |
затухание |
короткопериодиче- |
рких колебаний усиливается |
с увеличением |
v„ Поскольку |
эта |
величина пропорциональна |
производной |
и производным т р и |
т \, то рост этих параметров приводит к усилению затухания.
К тому же приводит увеличение геометрических размеров лета тельного аппарата, но увеличение массы и момента инерции дает противоположный эффект. Увеличение высоты полета и связанное с этим уменьшение плотности приводит к ослаблению затухания. Скорость полета на затухании может сказываться различным образом в зависимости от того, влияние каких фак торов преобладает.
Частота короткопериодических колебаний зависит в первую очередь от запаса устойчивости з„ . Увеличение запаса устой
чивости приводит к увеличению частоты. Заданную частоту ко лебаний можно получить, подбирая центровку летательного ап парата. При этом необходимо иметь в виду, что увеличение собственной частоты уменьшает статическую ошибку, но в то же время уменьшает передаточный коэффициент по углу от
клонения руля высоты &ов, ухудшая тем самым управляемость
летательного аппарата.
Увеличение геометрических размеров летательного аппарата обусловливает увеличение частоты короткопериодичеоких колеба ний, но увеличение момента инерции уменьшает ее, поэтому на практике увеличение размеров сопровождается уменьшением ча стоты. Увеличивается частота при увеличении скоростного напора.
Как видно из проведенного анализа, влияние запаса устойчиво сти аналогично влиянию жесткости в простейших упругих колеба тельных системах. Физически это объясняется тем, что увеличение запаса устойчивости приводит к увеличению стабилизирующего мо мента, возникающего при изменении угла атаки.
Из формулы (2.65) видно, что для сходимости длиннопериоди ческого движения необходимо, чтобы выполнялись условия
|
|
Т* > о |
|
|
|
|
|
(2.85) |
|
|
Г,,2 > 0 |
|
• |
|
|
|
|
|
|
|
|
|
|
|
|
26, Г, > 0 |
|
|
|
|
|
Первое из условий (2.85) приводится к виду |
|
|
|
|
c* + c?%PjpVS |
- у ( ] - |
^ r ) si n0> ° - |
(2-86) |
При слабой зависимости тяги |
от |
скорости, |
что характерно |
для ЖРД и РДТТ, и |
малых углах |
наклона |
траектории |
нера |
венство выполняется, так как коэффициент |
сх |
положителен, а |
произведение |
,. М |
хотя и является отрицательным при сверх |
|
звуковых скоростях, |
но по абсолютной величине, как |
правило, |
не превосходит сх . |
|
|
|
заменить |
неравен |
Два последних неравенства (2.85) можно |
ствами вида |
|
|
|
|
|
|
|
|
|
|
^ - > |
0, |
|
|
|
|
|
|
|
агах — а. |
> |
0. |
|
|
|
(2.87) |
Подстановка в левые части неравенств (2.86) значений вхо дящих туда величин показывает, что характер колебательной составляющей длиннопериодического движения зависит от всех факторов, определяющих величину и характер зависимости от условий полета сил, действующих на летательный аппарат.
Система уравнений, описывающих длиннопериодическое дви жение, имеет характеристическое уравнение третьей степени, приближенные значения корней которого можно определить, при равняв многочлен (2.65) нулю. Решая получившееся уравнение, находим
|
Рз ~ |
|
1 |
|
(2.88) |
|
|
т , |
|
|
|
|
1м |
|
|
Р4,5 |
|
^ |
. у "1 - |
1л2 |
(2.89) |
1Ч |
—J |
Т |
|
|
|
1Ч |
|
|
Преобразуя формулу (2.89) в соответствии с формулами (2.66), (2.62) и формулами для коэффициентов уравнений движения, находим
Р3: |
G |
сх + с" М pVS + |
' 1 |
I/ 3 |
sinO. (2.90) |
|
г |
VК2О С . |
11 |
Из выражения (2.90) видно, что затухание апериодической со ставляющей движения определяется характером изменения тяги двигателя л о скорости и силы лобового сопротивления. Чем быстрее уменьшается тяга и увеличивается лобовое сопротивление, тем сильнее затухание. Если тяга двигателя возрастает с увеличением скорости, то некоторые режимы полета могут быть неустойчивыми. Особенно опасен с точки зрения устойчивости полет при малых ско
ростях и больших положительных углах наклона траектории. Такие условия обычно имеют место при взлете летательного аппарата.
При увеличении удельной нагрузки летательного аппарата зату хание уменьшается.
Подобным же образом из формулы (2.89) можно получить вы ражения для вещественной и мнимой частей корней р.\ и р5, опреде ляющих колебательную составляющую длиннопериодического дви жения, но окончательные выражения получаются сложными и не удобными для анализа. Поэтому ограничимся рассмотрением физи ческой картины длиннопериодичеоких колебаний.
.Пусть в результате действия возмущения траектория летатель ного аппарата искривилась и он начал снижаться, сохраняя преж ний угол атаки (фиг. 2.4). В процессе снижения под действием ка сательной составляющей силы веса, направленной вперед, скорость полета начинает возрастать. Кроме того, по мере снижения увели чивается плотность воздуха. В результате подъемная сила увеличи вается, что приводит к новому искривлению траектории, и летатель ный аппарат начинает подниматься. В процессе подъема скорость уменьшается, так как касательная составляющая силы веса направ
лена назад, а плотность воздуха по мере увеличения высоты умень* шается. Происходящее в результате этих двух процессов уменьше ние скоростного .напора приводит к уменьшению подъемной силы, траектория вновь искривляется, летательный аппарат переходит на снижение и т. д.
Из приведенного описания сильно упрощенной картины длин нопериодических колебаний видно, что частота их определяется скоростью изменения скоростного напора, которая, наряду с другими факторами, определяется величиной производной р1) . Частота колебаний тем больше, чем больше по абсолютной ве личине производная рч, т. е. чем меньше высота полета. Увели чение массы летательного аппарата, напротив, должно приво
дить к уменьшению частоты. |
длиннопериоди |
Кроме упомянутых факторов, на характер |
ческих колебаний влияют закон изменения тяги |
в зависимости |
от высоты и скорости полета и др. |
|
ГЛАВА III
ПЕРЕХОДНЫЕ ПРОЦЕССЫ И ПЕРЕДАТОЧНЫЕ ФУНКЦИИ ЛЕТАТЕЛЬНОГО АППАРАТА В БОКОВОМ ДВИЖЕНИИ
§ 3 .1. Л И Н Е А Р И З О В А Н Н Ы Е У Р А В Н Е Н И Я Б О К О В О Г О Д В И
Переходные процессы в боковом движении у таких лета тельных аппаратов, как самолеты и управляемые ракеты, пред
ставляют собой сложные движения, в ходе которых изменяются угол скольжения £), угол крена к, а также угол рыскания и путевой угол Ф/(. На фиг. 3.1 изображена осциллограмма пере ходного процесса малой дозвуковой ракеты, вызванного скачко
образным отклонением руля направления. Как видно из осцил лограммы, процесс складывается нз быстрого изменения угловой
скорости |
крема, |
колебательного изменения угла (5 и угловых ско |
ростей w.v |
и му, |
которое затухает в течение нескольких секунд, |
и медленного изменения угла крена т и угла рыскания ф. |
Такой |
характер переходного процесса позволяет ввести неко |
торые упрощения, позволяющие получить достаточно простые окончательные выражения, отражающие, тем не менее, основные качественные закономерности переходных процессов в боковом движении. Прежде всего, учитывая, что изменение угловой ско рости крена происходит наиболее быстро, будем рассматривать изменение крена независимо от остальных движений.. Система уравнений, описывающих движение крена, записанная в проек циях на связанные оси координат, имеет вид
Л-
|
dj |
|
|
|
|
( 3 .1 ) |
|
= ш д. — |
(By C O S If t g |
& |
|
|
dt |
|
|
|
|
|
|
|
Момент крена |
от двигателя |
считаем |
равным нулю. |
Линеари |
зуя систему (3.1) |
и считая, что в |
расчетном движении |
угловые |
скорости iu.Vi> и (оУа |
равны нулю, т. е., обозначая |
|
|
|
Дц>у = |
(О |
|
|
|
|
|
|
У I |
|
|
систему (3.1) в отклонениях представим в следующем виде:
(3.2)
dky
U).v — cos To tg V y
dt
Если.считать параметры продольного движения неизменными, то ДЖ^ в общем .случае является функцией угла скольжения (5, угловых скоростей <ох и шу , а также углов отклонения элеронов 83 и руля направления оя . Поэтому, раскладывая ДЖ, в ряд и сохраняя только члены первого порядка малости, получаем
ДЖ , = Ж^Др + и , + A Q ев, + Ж‘* Д§3+ Ж Д8„+ ДЖ„ .
(3.3)
Подставляя (3.3) в (3.2), получаем линеаризованные уравне ния движения крена, которые в окончательной форме запишем так:
Р ш* = а "»х |
+д + О< " AА SЗ+« U X |
(3.4)
+“ у
Коэффициенты системы (3.4) определяются формулами (индекс янулик" опускаем):
1.
2. а' |
_ |
1 ДА? |
1 |
(3 р V a с / |
: _ V. |
~ 7 ГП* 2 S ’ |
Л |
(О |
|
1 я жШ |
1 |
Ш Р » |
л >* |
3 - |
Ч = л Л|- ’ = з ; я ' ’ т |
д а ' |
>1 |
о |
|
1 д^о |
1 |
о Р |
Л* |
4- |
a- r T x M! = T , m* ^ S '" - |
р |
О |
|
^ л жо |
^ |
о Р |
/> / |
5. |
|
|
|
|
|
|
6. |
|
= |
cos 7 tg ft, |
|
|
7.и ~ ~ Ш хе.
Производные тхх , т?х , от“v, да°н обычно получаются в
результате аэродинамического эксперимента и задаются как
|
|
|
|
|
|
функции числа М. При |
расчете коэффициентов уравнений (3.4) |
для ракет с крестообразно расположенными |
крыльями |
или бал |
листических1 ракет необходимо в формулах |
для a tv, |
а*а и а (°« |
вместо I |
подставить L, |
а в формулах для |
а “г и а “у |
вместо Р |
подставить |
2ZA |
|
|
|
Для получения системы уравнений движения рыскания вновь |
воспользуемся |
общими уравнениями движения в путевой системе |
координат. |
За |
расчетное движение примем прямолинейный полет |
в спокойной атмосфере с углом крепа, равным нулю. Будем счи |
тать, что угол наклона траектории в расчетном полете не пре- ■ вышает по абсолютной величине 20°, а путевой угол фл0 равен нулю. Изменением модуля скорости под. действием бокового
24 А. Г. Бедункович и др. |
369 |