книги из ГПНТБ / Динамика полета и конструкция крылатых летательных аппаратов
..pdfгде ftizб.г.о—коэффициент продольного момента без горизонталь ного оперения. тгг.о— — Лкг.0аг.0(а + <?ст —е+ пвЪв)—коэффициент продольного момента от горизонтального оперения.
Под продольной |
балансировкой в |
установившемся горизон |
||
тальном полете |
в дальнейшем будем |
понимать равновесие про |
||
дольных моментов в соответствии с уравнением (2.15). |
||||
Отклонение |
руля |
высоты (аэродинамического или газового), |
||
при котором |
обеспечивается равновесие моментов, будем назы |
|||
вать балансировочным отклонением. |
|
|||
Равновесие |
моментов будет обеспечено, если |
|||
. |
• |
|
т : .б.г.о + П 1 ;.го = 0 . |
Решая последнее уравнение для определения балансировочного отклонения руля высоты, имеем:
( тг.о.г.о |
а — ?ст + £ |
(2.16) |
|
Акгм аг.0 |
|||
|
|||
|
|
Пользуясь этой формулой, можно построить балансировоч ную кривую для отклонения руля высоты, т. е. зависимость балансировочного отклонения руля высоты от скорости или числа М полета. Для этого при различных заданных скоростях полета определяются т-в.г.о, - угол атаки а, скос потока s и по формуле (2.16) вычисляется балансировочное отклонение руля высоты ов. На фиг. 2.8 приведена типичная балансировочная кри вая для дозвуковых самолетов.
На вед ‘балансировочной кривой сильное влияние оказывают по ложение центра тяжести, фокуса КЛА, угла установки стабилиза тора и изменение аэродинамических характеристик, вызванное влия нием сжимаемости воздуха.
На фиг. 2.9 и 2.10 приведена качественная зависимость баланси ровочной кривой от угла установки стабилизатора и центровки, а на фиг. 2.11 показан типичный вид балансировочной кривой совре-
^ менного самолета.
230
Балансировочные кривые могут быть построены для различных режимов работы силовой установки и высот полета КЛА.
Балансировочное отклонение руля высоты, определяемое из (2.16), может быть использовано при программировании установив
шегося горизонтального полета.
В качестве программы установившегося движения КЛА по дуге
.большого круга вокруг Земли (особенно при значительных дально стях полета) может быть рекомендован следующий вариант:
а) Н = const, т. е. постоянная высота полета. Эта программ может быть осуществлена обычным высотомером в комплексе с ап
паратурой управления рулем высоты;
б) при движении по дуге большого круга курсовой угол ф (угол между связанной осью х х и касательной к меридиану)
Фиг. 2.10
будет меняться по определенному закону 6= ili(£), который может быть заранее рассчитан по известной скорости установившегося полета. Рассчитанный закон изменения курсового угла может быть выдержан, например, гирокомпасом, связанным с рулем направления;
в) угол крена у=0, т. е. полет без крена. Эта программа движения может быть реализована, если использовать гировер тикаль, показания которой могут быть связаны с отклонениями элеронов;
г) угол скольжения £1—0. При выполнении данной программы может быть использован обычный флюгер.
Весьма важно, и это представляет значительную трудность, что бы была согласована работа отдельных Чувствительных элементов, выдерживающих ту или иную программу, в особенности в тех слу чаях, когда будут отклонения от заданной программы из-за возму щающих причин (поток воздуха, ударная волна и т. д.).'
Отклонение руля высоты, найденное из (2.16) при данной скоро сти полета, должно быть меньше максимального конструктивного значения с тем, чтобы иметь запас хода руля примерно 4—6° на уп равление.
2Л
Приведенная программа движения не является единственной • и могут быть предложены другие программы, учитывающие специ фику задачи и особенности аппаратуры управления.
§ 2.7. РАЗГОН И ТОРМОЖЕНИЕ ПО ПРЯМОЛИНЕЙНОЙ ТРАЕКТОРИИ. ПРОГРАММА ДВИЖЕНИЯ КРЫЛАТОГО ЛЕТАТЕЛЬНОГО АППАРАТА
В динамике боевого маневрирования разгон и торможение КЛА по прямолинейной траектории занимает важное место. Основными характеристиками разгона и торможения являются время и прой денный путь при выполнении данного маневра.
Для расчета характеристик разгона и торможения по наклон ной прямолинейной траектории пз (1.41), имея в виду, что
6=const и dO — 0, получим следующие основные уравнения:
- j |
= Р соя (а + ?ae) - Q - G sin 0 , |
(2.17) |
||
Р sin (а + tpas) + Y = G cos 0 f 1 — |
— ] . |
|
||
|
|
\ |
V Z k o c ) |
|
Как известно, |
|
|
|
|
|
dH |
V sin 0 |
|
(2.18) |
|
d t ~ |
|
||
и в общем случае |
|
|
|
|
|
G = G0- |
Qcext' |
|
(2.19) |
где qeeK— постоянный секундный расход топлива.
Таким образом, для расчета прямолинейного маневра следует
воспользоваться уравнениями (2.17), |
(2.18) и (2.19). |
|||
Как |
это |
следует из (2.17), |
при |
разгоне Pcos (a-t-cpa) — Q— |
—О sin |
9 > 0, |
т. е. суммарная |
результирующая сила направлена |
|
по скорости. |
В противном случае движение совершается с тор |
можением.
При подъеме КЛА по прямолинейной траектории составляю щая силы веса GsinO направлена против движения, а при сни жении—наоборот. Следовательно, при подъеме характеристики разгона ухудшаются, а характеристики торможения улучшаются.
Очевидно, что4 при |
снижении |
будет иметь |
место обратное яв |
|
ление. |
|
|
и (2.19) является нелинейной |
|
Система уравнений (2.17), (2.18) |
||||
и в общем виде не |
интегрируется. |
Поэтому |
при расчете харак |
|
теристик разгона и |
торможения |
приходится |
пользоваться мето |
дом численного интегрирования или вычислительными машинами.
232
Для выполнения численного интегрирования уравнений двш жения следует воспользоваться известными соотношениями
d V _ .
(2.20)
dt J
и
(2.21)
Здесь у—касательное ускорение, которое определяется из пер-5 вого уравнения системы (2.17), т. е.
. |
d V |
р |
Q |
(2. 22) |
] ~ |
dt ~ g |
~q ~cos (а+'?а)— -q — sin* |
Тогда последовательность расчета характеристик разгона и
торможения сводится к следующему: |
|
|
||||
|
1) задаются интервалом времени Дt\ |
(G)0, (а)0, V0i |
||||
Я 0 |
2) |
по начальным |
значениям параметров (Я)0, |
|||
и при заданных 0 |
и <рйв определяют ускорение |
(у ),; |
||||
|
3) |
принимая (у’)0 = |
const в заданном диапазоне |
изменения вре |
||
мени At, вычисляют приращение скорости A V —(j)t A t; |
||||||
т. |
4) |
определяют величину |
скорости |
в момент времени t ^ A t , |
||
е. |
|
V i = V 0 + AV\- |
|
|||
|
|
|
|
|||
|
5) |
вычисляют приращение |
пути |
|
|
|
где |
|
AL — VtpA t, |
|
|||
|
v |
Ve+ V t. |
|
|||
|
|
|
|
|||
|
|
|
v ср |
2 |
3 |
|
изменение высоты |
ДЯ = |
Vep sin |
, |
|
||
|
|
|
|
изменение веса
=ЯсенЫ
исоответствующие значения высоты
Я 1= |
Я 0 + ДЯ |
|
и веса |
- |
■ |
Gi = G0 — ДО; |
|
233
6) |
|
зная |
величины |
параметров |
VX) H lt |
Gx в момент врем |
||||||||
ни tx аналогичный расчет повторяется. |
Однако предварительно |
|||||||||||||
необходимо |
определить |
величину коэффициента сW, соответст |
||||||||||||
вующую |
параметрам |
Vx и Н х. Для этого из второго уравнения |
||||||||||||
системы (2.17), |
принимая |
sin (а + |
|
^ |
(a+ cpd) = |
+ |
+ |
на |
||||||
ходится |
коэффициент |
подъемной силы |
|
|
|
|
||||||||
|
|
|
|
<М 1 |
|
|
cos О— Я, (а0 + |
<рД |
|
|
||||
|
|
|
|
|
|
Р г |
|
| |
Р . У ,2 |
с |
|
|
|
|
|
|
|
|
|
|
с* |
"г- |
о |
^ |
|
|
|
||
|
|
|
|
|
|
су, |
|
|
z |
|
|
|
|
|
и из поляры |
объекта |
в зависимости |
от |
сб> и М, — — ( а ,— ско- |
||||||||||
рость |
звука |
на |
|
|
|
|
|
|
|
у |
а, |
4 |
|
|
высоте И х) определяется с<!>. |
|
|
|
|||||||||||
В |
случае |
горизонтального |
|
разгона |
или торможения |
при |
||||||||
G= £onst |
касательное |
ускорение |
является функцией только |
ско |
рости, так как коэффициенты су , сх , с* и а0 зависят только от числа М = . Имея в виду, что
VdV
dL
j
получим
(2.23)
(2.24)
I
Следует иметь в виду, что при торможении двигатель дрос селируется и могут быть использованы тормозные устройства.
В соответствии с изменением. V, Н и а будет меняться угол отклонения руля высоты. При разгоне и торможении отклоне ние руля высоты нельзя определять из (2.16), ибо продольные моменты не будут уравновешены в силу наличия вращения КЛА, вызванное изменением угла атаки. Таким образом, для обеспе чения полученного закона изменения угла атаки а руль высоты должен отклоняться в соответствии с уравнением моментов си стемы (2.1), т. е.
234
qSba ( т 2б.г.о + n i Z!.0).
Решая последнее дифференциальное уравнение относитель но получим необходимый закон отклонения руля высоты при разгоне и торможении по прямолинейной траектории:
(2.25)
Таким образом, изложенный метод позволяет определить все характеристики разгона и торможения, в том числе параметры, необходимые для программирования движения. Например, если движение программируется по углу тангажа, то для обеспече ния прямолинейности движения угол тангажа &должен меняться по закону
%{t) = 0 + «(*), |
(2.26) |
где 0—заданный постоянный угол траектории с горизонтом (он равен нулю при горизонтальном полете),
я (?)—закон изменения угла атаки от времени, определяемый
у
Программа движения г>(t) может быть осуществлена, если в качестве чувствительного элемента, реагирующего на откло нения от заданной программы, будет использована гировер тикаль.
В качестве программы движения в данном случае может быть принят также закон изменения высоты от времени, который получится в результате решения уравнений движения, т.' е.
H = H (t) . |
(2.27) |
В качестве чувствительного элемента, обеспечивающего про грамму движения (2.27), можно использовать чувствительный высотомер, работающий без запаздывания, или соответствующую инерционную систему.
Кроме того, на КЛА должна устанавливаться аппаратура, позволяющая выдерживать полет без крена и скольжения, т. е.
7=0 и |0=О. |
(2.28) |
ГЛАВА III
ДВИЖЕНИЕ КРЫЛАТОГО ЛЕТАТЕЛЬНОГО АППАРАТА В ВЕРТИКАЛЬНОЙ ПЛОСКОСТИ
§3.1. УРАВНЕНИЯ ДВИЖЕНИЯ КРЫЛАТОГО ЛЕТАТЕЛЬНОГО АППАРАТА В ВЕРТИКАЛЬНОЙ ПЛОСКОСТИ
Движение крылатых летательных аппаратов (.КЛА) в вертикаль ной плоскости является одним из мастных случаев общего случая движения КЛА по пространственным траекториям, когда вектор ско рости полета лежит в плоскости симметрии КЛА, совпадающей с
вертикальной плоскостью, т. е. движение совершается без крена и скольжения. Примерами таких движений могут служить движение
крылатой ракеты на активном участке траектории |
полета, подъем' |
||
и снижение оамолета, взлет и посадка, |
движение |
по |
траектории, |
представляющей замкнутую петлю в |
вертикальной |
плоскости, |
и т. д. |
Типичные траектории движения КЛА в вертикальной пло |
|||
скости |
изображены на фиг. |
1.1 раздела II. |
Силы, действующие |
|
на |
КЛА при движении в |
вертикальной |
плоскости, показаны |
|
на |
фиг. 3.1. |
|
|
Движение КЛА в вертикальной плоскости определяется двумя уравнениями сил и одним уравнением моментов. Так как в общем случае движение КЛА в вертикальной плоскости является криволи нейным и неустановившимся, то в проекциях на скоростные оси Ох и Оу уравнения движения с учетом кривизны земли были записаны в следующем виде (1.41):
236
V |
= Р cos (а + <р5в) — Q — G sin 0 |
|
|
о l / ^ - = P sin(a-b(Pde)+}/- ^1 — |
j G cos 0 |
(3.1) |
, d 2(а + 0 — <р)
Мг
d t2
Уравнения (3.1) в общем случае являются нелинейными и могут быть решены только при заданной программе движения либо с помощью электронных вычислительных машин, либо приближенными методами интегрирования. Ниже рассмотрим некоторые виды движений КЛА. в вертикальной плоскости и общие основы их расчета.
§ 3.2. ОСНОВЫ РАСЧЕТА ДВИЖЕНИЯ КРЫЛАТОГО ЛЕТАТЕЛЬНОГО АППАРАТА В ВЕРТИКАЛЬНОЙ ПЛОСКОСТИ
Расчет движения центра тяжести
Для расчета движения центра тяжести КЛА воспользуемся первыми двумя уравнениями системы (3.1), которые запишем в следующем виде:
|
dV_ |
COS (a + |
<PJ |
Q |
|
|
|
|
dt |
■g | -Q |
G |
Sln |
|
||
|
|
|
|
(3.2) |
|||
dO |
JL |
P . |
|
Y |
|
V2' |
|
|
1 - |
cos 6 |
|||||
d t |
V |
- q s in (a + |
?<3e) + |
~q — |
V2 |
||
|
|
|
|
|
|
V1 KOt |
|
В правые части уравнений (3.2) входят сила тяги Р, аэроди намические силы Q и Y и переменная в общем случае сила веса О.
Тяга Р в общем случае определяется скоростью V, высотой полета Н и режимом работы двигателя. При заданном режиме работы двигателя, как известно,
P = P ( V , Н ).
Аэродинамические силы Q и Y определяются тремя пара метрами: a, V и Н, т. е.-
Q = Q («■ V, Я) = с , К М, Re)p(//2) V/2S,
Y = Y (* t V, Н ) —су (а, М, Ке )р(Я9)У 25 .
237
При этом |
аэродинамические характеристики сх=сх (а, М, Re) и |
|
Су = су (а, |
М, Re) должны быть |
известны. |
Для определения закона изменения веса КЛА необходимо |
||
знать секундный расход топлива |
qеек — qceK ( t ) . Тогда |
|
|
|
(3.3) |
где G0—начальный вес КЛА. |
КЛА в вертикальной плоскости |
|
Кроме того, при движении |
имеет место кинематическое соотношение
(3.4)
В уравнениях (3.2) и (3.4) содержатся четыре независимых параметра: а, V, Н и 0. Следовательно, без дополнительного условия эти уравнения не решаются. Поэтому для расчета дви жения центра тяжести КЛА необходимо задать программу дви жения. В рассматриваемом случае программой движения КЛА может служить, например, одна из зависимостей
|
|
|
а = a (t) , |
|
|
|
|
е = 0 (о , |
|
|
|
» = |
а -f 0 = 0 (t) . |
|
Кроме того, необходимо задать начальные условия, |
т. е. V = V 0, |
|||
Н — Н0, G = G0 при t= t0. |
уравнений (3.2) вместе |
с уравнения |
||
Таким |
образом, |
система |
||
ми (3.3), |
(3.4) при |
заданной программе движения |
и начальных |
условиях позволяет определить все параметры, характеризующие движение центра тяжести КЛА.
Рассмотрим метод численного интегрирования системы урав
нений (3.2), |
которую вместе с уравнениями (3:3) и (3.4) запишем |
в конечных |
разностях: |
АН — Vsin ОДt
AG— qeeKAt
238
В качестве программы движения примем зависимость
»=& (*),
тогда начальные условия будут определяться значениями еле*
дующих параметров при t = t 0 : 0 = 0О, |
V — V0, # |
= # 0, G = G0< |
|||||
Расчет движения центра тяжести .КЛА выполняется в следую |
|||||||
щем порядке: |
|
времени |
. |
|
|
|
|
1. |
Задается приращение |
|
qceK, |
а = & — 0, |
|||
2. |
Определяются значения величин |
Р, Q, Y, |
|||||
& при |
t = t 0. |
|
|
|
|
|
|
3. По формулам (3.5) вычисляются AVy, AOj, АН и AGx .„ |
|||||||
4. |
Находятся значения: |
Vx— V0-f AVy, |
01 = Oo + |
A0i |
H 1= H 0 + |
||
+ AHlt Gy = G0-f- AG, соответствующие |
моменту |
времени tx — |
|||||
—t0 + |
A^i . |
расчет |
выполняется |
|
аналогичным |
образом, |
|
5. |
Дальнейший |
|
|||||
т. е. задается At2, |
и по значениям указанных в пункте 2 вели |
||||||
чин при tx = t0¥ Atx вычисляются V2, 02, |
Н 2, G, |
и т. д. |
Расчет углов отклонения рулей
Для выполнения криволинейного движения в вертикальной пло-' скости л-1еобходимо отклонять руль высоты или иной орган продоль ного управления. При этом создается продольный момент, под дей ствием которого КЛА приобретает угловое ускорение, изменяет угол' тангажа и угол атаки, а следовательно, и нормальную перегрузку. В результате изменения перегрузки траектория движения центра тяжести КЛА искривляется, причем кривизна траектории изменяется в зависимости от отклонения органа продольного управления.
Возможность осуществления того или иного движения в верти кальной плоскости определяется не только величиной располагаемой перегрузки или коэффициента су, но и величиной потребных для вы полнения этого движения -углов отклонения рулей. При недостаточ ной эффективности рулей может оказаться,'что потребные для вы полнения заданного движения углы отклонения рулей превосходят их предельные значения и, очевидно, в этом случае выполнить за данное движение будет невозможно. Поэтому необходимо произве сти расчет потребных углов отклонения рулей. Знание этих углов представляет значительный интерес и для пилотируемых летчиком КЛА, так как они характеризуют его управляемость и маневренные
В О З М О Ж Н О С Т И . |
а |
В качестве одного ив основных органов продольного управления при движении КЛА в вертикальной плоскости служит руль высоты или управляемый стабилизатор. Возможно также управление с по мощью газовых рулей и изменения направления тяги двигателя. По этому в зависимости от способа управления необходимо прежде все го определить либо потребные углы отклонения руля высоты или стабилизатора, либо потребные углы отклонения газовых рулей,
239