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

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

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

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

центральных осей инерции

относительно корпуса

аппарата

[6, 9] дает в уравнениях (2.46) дополнительные

члены,

которые

в большинстве практических

задач являются

пренебрежимо

малыми.

 

 

 

§4. ОСИ КООРДИНАТ, ГЕОМЕТРИЧЕСКИЕ

ИКИНЕМАТИЧЕСКИЕ СООТНОШЕНИЯ

4.1. ОСИ КООРДИНАТ

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

З е м н а я с и с т е м а о с е й к о о р д и н а т

Оси этой системы Охз уз z 3 (рис. 2.3) неподвижно связаны с Землей и участвуют в ее суточном вращении. Сокращенно их называют земными осями.

Рис. 2.3. Земная система осей координат:

М — точка

старта; NABS — начальный

(Гринвичский) меридиан; NMLS — местный

меридиан;

фц— геоцентрическая широта

точки М; к — долгота точки М ; сс — каса­

тельная к местному меридиану в точке М; ф — азимут запуска

Начало земных осей выбирается где-либо на поверхности Земли, например, в точке старта летательного аппарата. Ось Оуз направлена по радиусу-вектору, проведенному из центра общего земного эллипсоида через точку старта; ось Охз ,— го­ ризонтально в каком-либо произвольном направлении; н’апри-

91

мер, ось О хз удобно направить так, чтобы плоскость полета совпадала с плоскостью Охзуз', ось Охз образует с плоскостью местного меридиана угол г|з, называемый азимутом запуска-, ось ÖZ3 направляется так, чтобы система координат была правой.

С т а р т о в а я с и с т е м а к о о р д и н а т

Стартовая система координат Oxcyczc (рис. 2.4) также свя­ зана с Землей и вращается вместе с ней. Начало координат рас­ положено в точке старта; ось Оус направлена вверх по линии отвеса, т. е. противоположно направлению силы тяжести; ось Охс образует с плоскостью местного меридиана угол азиму­ та запуска г|г, ось Огс соответствует правой системе координат.

Рис. 2.4. Стартовая система осей координат:

М — точка

старта; NABS — начальный

(Гринвичский)

меридиан;

NMLS — мест­

ный

меридиан;

ф д — астрономическая

ши­

рота

точки

М;

X — долгота

точки М;

сс

касательная

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

 

ке М; ^ — азимут запуска

 

Рис. 2.5. Ориентация связанных осей на стар­ те летательного аппара­ та

На старте связанные оси баллистической ракеты ориентиру­ ются по осям стартовой системы (рис. 2.5). Продольная ось ра­ кеты Ох 1 совмещается в осью Оус\ поперечная ось Оу\ направ­ ляется в сторону, противоположную оси Охс\ ось Ozi направ­ ляется по оси Ozc.

Н а ч а л ь н а я с т а р т о в а я с и с т е м а о с е й к о о р д и н а т

Оси начальной стартовой системы координат Ox0yoZo в мо­ мент пуска совпадают с осями стартовой системы. В дальней­

92

шем оси начальной стартовой системы не изменяют своего пер­ воначального направления относительно инерциального про­ странства, а оси стартовой системы, жестко связанные с Зем­ лей, поворачиваются за время t на угол о)3 t вокруг оси враще­

ния Земли. Начальная стартовая система координат является инерциальной системой координат.

 

 

 

 

 

 

 

 

 

 

 

/

Г е о ц е н т р и ч е с к а я с и с т е м а к о о р д и н а т

Эта система координат с началом в центре Земли и с осями,

связанными

с Землей, уже использовалась выше при

изучении

поля земного тяготения. Опорны­

 

 

 

 

ми плоскостями

в рассматривае­

 

 

 

 

мой системе

координат

служат

 

 

 

 

плоскости экватора

и начального

 

 

 

 

меридиана.

 

 

 

 

 

 

 

 

 

 

Положение центра масс лета­

 

 

 

 

тельного

аппарата

при этом

мо­

 

 

 

 

жет определяться

или тремя

де­

 

 

 

 

картовыми координатами х, у, z,

 

 

 

 

или, что более удобно, тремя сфе­

 

 

 

 

рическими координатами к, срц, г.

 

 

 

 

Долгота к и геоцентрическая ши­

 

 

 

 

рота фц отсчитываются, как было

 

 

 

 

показано

выше,

на рис. 1.3. Ко­

 

 

 

 

ордината

г

 

представляет

собой

 

 

 

 

расстояние

 

от

центра

Земли

Рис. 2.6.

Местная географиче­

до центра масс летательного ап­

ская система координат

парата.

 

 

 

 

 

 

 

 

 

 

 

М е с т н а я г е о г р а ф и ч е с к а я с и с т е м а

к о о р д и н а т

Начало

этой системы

координат

OxTyTzT (рис.

2.6)

совпада­

ет с центром

масс

летательного аппарата;

ось Охг

проведена

параллельно касательной к меридиану места в северном направ­ лении; ось Оуг направлена по радиусу-вектору г; ось Oz^ парал­ лельна плоскости экватора.

С к о р о с т н а я и п о л у с к о р о с т н а я с и с т е м ы

ко о р д и н а т

Внекоторых задачах динамики уравнения движения центра масс летательного аппарата удобно записывать в проекциях на оси координат, связанные с вектором скорости центра масс V. Начало координат О такой системы располагается в _центре масс аппарата; ось Ох направлена по вектору скорости V, т. е. по касательной к траектории в сторону полета; оси Оу и Oz ле­ жат в плоскости, нормальной к траектории полета. При этом в

93

динамике полета ось Оу выбирают как в плоскости симметрии летательного аппарата Охірі, так и в вертикальной плоскости. Первую систему координат будем называть скоростной Oxyzr вторую — полускоростной Oxy^z*.

Матрицы направляющих косинусов

К о с и н у с ы у г л о в м е ж д у о с я м и с в я з а н н о й и н а ч а л ь н о й с т а р т о в о й с и с т е м к о о р д и н а т

Ориентация летательного аппарата относительно начальной стартовой системы координат определяется тремя углами между

 

 

 

связанной

OxiyiZi

и началь­

 

 

 

ной стартовой

ОхоУо%о систе­

 

 

 

мами

координат

(рис. 2.7):

 

 

 

углом рыскания g — меж­

 

 

 

ду

проекцией

продольной

 

 

 

оси

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

аппарата

 

 

 

Ох1

на

плоскость

Ox0z0

и

 

 

 

осью

Охо;

 

 

 

 

 

 

 

 

 

углом тангажа фт — меж­

 

 

 

ду продольной осью аппара­

 

 

 

та

ОхI

и плоскостью OxoZ0;

 

 

 

углом

крена т) — между

 

 

 

поперечной

осью

 

Оу\

и

 

 

 

плоскостью,

проходящей че­

 

 

 

рез оси

ОхI

и

Оу0.

между

 

 

 

 

Косинусы

 

углов

Рис. 2.7. Ориентация

связанных

осей

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

и

началь­

относительно начальных

стартовых

осей

ной

стартовой

систем коор­

динат приводятся в табл. 2.1. Найдем косинусы углов между осями начальной стартовой и связанной систем координат. Для этого через центр масс лета­ тельного аппарата — начало связанной системы координат, про­ ведем оси OxoyoZo, параллельные осям начальной стартовой си­ стемы координат. Повернем эту систему на угол \ вокруг оси Оу0 так, чтобы плоскость Ох'у' прошла через ось Ох,. Обозна­

чим полученную систему через Ox'y'z'.

Формулы преобразования координат удобно записывать в. матричной форме. Обозначая матрицу-столбец с элементами x', y', z' через [x'], т. е.

ГX'

ИУ

z'

94

Таблица 2.1

Косинусы углов между осями связанной и начальной стартовой систем координат

Оси

Ох і

Оу і

O z \

О х а

O y a

cos £ cos <рт

sin <pT

— cos £ sin <pTcos 7j +

COS <pr COST]

+ Sin £ Sin Yj

 

cos £ sin y T sin Y] +

— COS <pTSin Yj

+ sin £ cos Yj

 

O z 0

— sin £ cos <pT

 

cos

£ sin 1) +

+

sin

£ sin <pT cos Y)

 

cos

£ cos Yj —

sin

£ sin <pT sin T]

а с элементами х0, уо, z0 через |хо], получим следующую форму­ лу перехода от системы Ox0yoZo к системе Ox'y'z':

[*']= Г е[х0],

(2.47)

где матрица Tg перехода от системы координат Ox0yoZo к систе­ ме Ox'y'z' имеет вид

cos £ 0 — sin £

*

0

1

0

(2.48)

sin t

0

COS £

 

 

 

Повернем систему Ox'y'z' вокруг оси Oz' на угол фт так, чтобы ось Ох" совпала с осью Оху, полученную систему обозна­ чим через Ox"y"z". Соответствующая формула преобразования координат при повороте осей на угол <рт будет иметь вид

[■*"]= ГЛ*'],

 

(2.49)

где

 

 

 

cos <+>T

sin <pT

0

 

sin <PT

COS cpT

0

(2.50)

0

0

1

 

Поворотом вокруг оси Ох" на угол т] совместим оси Оу"' и Oz'" с осями Оу 1 и Oz\. Переход от системы Ox"y"z" к системе OxiyiZi осуществится посредством формулы

[•*і]= Г, [-*"],

(2.51)

95

где

1

0

0

 

0

COS 7j

sin Tj

(2.52)

0

— Sin 7)

COS Tj

 

Подставляя выражения (2.47) и (2.49) в зависимость (2.51), получим матричное уравнение перехода от начальных старто­ вых осей к связанным:

[а:1]= Г Т(Г9Г?[ л:0]= Г [ л:0],

(2.53)

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

cos 5 cos <pr1

sin 9

T

 

— sin 5 cos f T

sin, 5 sin

 

) —

COS yT COS

У)

COS$Sinr;-|-

— cos£ sin<f>TcosT)

 

+ sin£,,sin <fT cos T,

Sin 5

COS

Y

+

— cos<pT

simrj

cos 5 cos г; —

+ cos£ sintf.r sin-fj

— sin$ sin<pTsinij

К о с и н у с ы у г л о в м е ж д у о с я м и с в я з а н н о й и м е с т н о й г е о г р а ф и ч е с к о й с и с т е м к о о р д и н а т

Ориентация

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

относительно местных

географических осей определяется тремя

углами ф, ■0', у, анало­

гичными углам

фт, т).

 

Углом рысканья ф называют угол между местной географи­ ческой осью ОхГ и проекцией продольной оси Oxj на горизон­ тальную плоскость.

Таблица 2.1а

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

Оси

O x i

O y i

O z 1

 

0*г

 

O yT

 

0 z r

cos ф cos ft

 

sin ft

 

sin ф cos ft

— cos ф sin

 

ft cos 7 +

cos ft cos Y

cos Ф sin Y +

+

sin

 

ф sin Y

— cos ft sin Y

+

sin

ф sin ft cos Y

cos ф sin ft sin

Y +

 

cos

4

cos Y —

+

sin

4

cos

Y

 

sin

ф sin ft sin Y

 

 

 

 

 

 

 

 

 

96

Углом тангажа ■й называют угол между продольной осью Охі и горизонтальной плоскостью.

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

Углы ф, & и у показаны на рис. 2.8. Направляющие косину­ сы углов между связанными и местными географическими ося­ ми координат приведены в табл. 2.1а.

К о с и н у с ы у г л о в м е ж д у о с я м и п о л у с к о р о с т н о й и м е с т н о й г е о г р а ф и ч е с к о й с и с т е м к о о р д и н а т

В полускоростной системе координат ось Оу* находится всегда в вертикальной плоскости, а ось Oz* всегда горизонталь­ на. Поэтому положение полускоростной системы координат от­ носительно местной географической определяется только на­ правлением оси Ох, т. е. следующими двумя углами (рис. 2.9):

Рис.

2.8. Ориентация

связанных

Рис. 2.9. Переход от местной

осей

относительно местных гео­

географической системы

коор­

 

графических

осей

динат к

полускоростной

систе­

 

 

 

 

ме координат

 

1) углом Ч1-2между проекцией

вектора скорости V на гори­

зонтальную плоскость и осью Охти

_

 

2)

углом Ѳ между вектором

скорости

V и горизонтальной

плоскостью.

Угол Ѳ обычно называют углом наклона траектории к гори­ зонту, угол 4я будем называть углом поворота траектории.

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

местной

географи­

ческой системы координат ОхГуГгГ на углы ¥

и Ѳ до совпадения

направления ее осей с полускоростной системой

координат

Oxy*z* (см. рис. 2.9). В результате получаем:

 

 

4-3422

97

[х]= Л [х г];

(2.55)

COS 4’ COS0

sin Ѳ

— sin Ч" cos0

л = — соз 4' sin Ѳ

cos Ѳ

sin Ч- sin 0

sin Ч’

0

cos Ч-

К о с и н у с ы у г л о в м е ж д у о с я м и м е с т н о й г е о г р а ф и ч е с к о й и с к о р о с т н о й с и с т е м к о о р д и н а т

Положение скоростной системы координат относительно местной географической будем определять углами Ч*-, Ѳ и ус, аналогичными углам ф, -д и у. Тогда угол ус будет представлять собой угол между осью Оу и вертикальной плоскостью Охуг, проходящей через вектор скорости полета. Этот угол уе будем называть углом крена. Однако всегда следует помнить, что угол крена ус отличается от угла крена у.

Нетрудно видеть, что положение скоростной системы коор­ динат относительно полускоростной определяется только углом ус. Когда ус= 0, т. е. плоскость симметрии летательного аппара­ та совпадает с вертикальной плоскостью, скоростная и полуско­ ростная системы осей совпадают.

В табл. 2.2 приведены косинусы углов между полускорост­ ными и скоростными осями, а в табл. 2.3 — между полускорост­ ными и местными географическими осями.

 

 

 

Таблица 2.2

 

 

 

Т а б л и ц а 2.3

Косинусы углов между полу-

Косинусы углов между полускорост-

скоростными и

скоростными

ными

и местными

географическими

 

 

осями

 

 

осями

 

Оси

Ох

Оу

Oz

Сси

О х Г

Оут

° z r

O x

1

0

0

О х

cos Ѳ cos ¥

sin Ѳ

— cos Ѳ sin ¥

О у *

0

cos 7с

— sin 7с

Оу*

—sin0 cos'?

cos Ѳ

sin Ѳ sin ¥

О г *

0

sin 7 с

cos 7с

О г *

sin ¥

0

cos ¥

Косинусы углов между скоростными и местными географи­ ческими осями можно определить с помощью табл. 2.1а, -если вместо углов ф, -O', у подставить соответственно углы Ч^, Ѳ, ус-

98

4.2. СВЯЗЬ МЕЖДУ УГЛАМИ Ѳ, ус и ф, ft, у, ß, а

' Положение скоростной системы координат относительно­ местной географической можно определить либо непосредствен­ но углами Чг, Ѳ, ус, либо углами ф, й, у, ß, а. (Углы ф, й, у опре­ деляют положение связанных осей относительно местных геогра­ фических, а углы ß, а — положение скоростных осей относитель­ но связанных).

Следовательно, углы Т, Ѳ, ус можно выразить через углы ф, й, у, а, ß. Зная из табл. 2.1а и 1.1 косинусы углов между свя­ занными и местными географическими осями, между связанны­ ми и скоростными осями, можно получить следующие три соот­ ношения:

sin Ѳ= sin & cos а cos ß— cos ft cos у sin а cos ß —

 

—cos ft sin у sin ß;

 

sin T" cos Ѳ= sin фсов &cos а cos ß -f- cos ф sin у sin а cos ß-J-

 

-f- sin ф sin ftcos у sin а cos ß — cos ф cos у sin ß -f-

J. (2 57)

-j- sin ф sin &sin Y sin ß;

 

sin Yccos Ѳ —sin 8 cos а sin ß — cos ft cos у sin а sin ß 4- -f- cos 8 sin у cos ß.

Из последнего равенства системы (2.57) следует, что при ма­ лых углах Ѳ, •O', а, ß, у и уо углы крена у и ус приблизительно равны: Y~Yc-

Рассмотрим некоторые частные случаи. Пусть полет совер­ шается в вертикальной плоскости без крена (у = 0). Так как при этом ß=0, то из первого равенства системы (2.57) получим

Ѳ —ft — а .

Пусть теперь полет происходит в горизонтальной плоскости

без крена. Положим c o s a ^ l, sin a ~ 0 , COS8 Ä ;1. Тогда из вто­ рого равенства системы (2.57) следует

Ф’ ^ ф - ß .

4.3. КИНЕМАТИЧЕСКИЕ УРАВНЕНИЯ, ОПИСЫВАЮЩИЕ ВРАЩЕНИЕ ЛЕТАТЕЛЬНОГО АППАРАТА ВОКРУГ ЦЕНТРА МАСС

Общий случай

Для исследования полета летательного аппарата необходи­ моиметь кинематические уравнения, описывающие изменение угловых координат летательного аппарата фт, £ и г] в зависимо-

4*

99

сти от проекций вектора угловой скорости летательного аппара­ та на связанные оси юх, tov, coz.

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

скорости I направлен по оси Оуо, вектор <рт — по оси Oz', а век­

тор т] — поосиОяі.

Вектор угловой скорости аппарата со можно представить как сумму

w= £4 CPT+ Yi-

(2.58)

Косинусы углов между векторами £, <рт и г) и связанными осями даны в табл. 2.4. Пользуясь этой таблицей, получим, что проекции угловой скорости аппарата на связанные оси равны:

 

 

d т)

d j

sin f T;

 

 

 

 

 

 

dt

dt

 

 

 

 

/

 

 

dj

COS tpT COS Tj-

 

 

 

 

 

 

Sin YJ;

 

(2.59)

 

 

dt

dt

 

 

CO

dv?

COS T j------------- COS cp

S in T j.

 

 

 

Z

—1—

 

 

 

 

dt

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 2.4

Оси

 

 

k

 

 

'P-r

 

TI

О х 1

sin <рт

 

 

0

 

1

 

 

 

 

sin 1}

1

0

 

 

COS cpTCOS

 

 

О у і

 

 

 

I

O z

!

— COS <pTsin

f]

 

cos тд

 

0

 

 

 

 

Составим теперь кинематические уравнения, устанавливаю­ щие связь между производными углов ф, Ф и у по времени и проекциями вектора угловой скорости летательного аппарата на связанные оси сож, соу, coz-

Вектор угловой скорости летательного аппарата ш можно представить как сумму угловых скоростей поворота связанных осей относительно географических, географических осей относи­ тельно некоторых земных осей и земных осей относительно инер­ циальных, например, начальных стартовых осей. В качестве осей, связанных с Землей, удобно взять географические оси OoXroi/roZro в некоторый начальный момент полета.

100