Belokonov
.pdfРис.7. Траекторная система координат.
1.2.5 Пересчет из одной системы координат в другую
Задача заключается в пересчете проекций вектора из одной системы координат в другую.
Рассмотрим это на примере пересчета из нормальной системы координат в связанную систему координат (рис.8)
a axg iд ayg jд azg kд
Рис.8. Вектор a axg ,ayg ,azg
|
|
|
|
ax |
|
|
|
|
|
|
|
|
|
|
axд |
|
|
|
a |
|
|
ay |
|
; |
|
|
a |
|
g |
|
ayд |
|
|||
|
|
св |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
az |
|
|
|
|
|
|
|
|
|
|
azд |
|
|
|
|
|
|
|
XXg |
|
|
|
|
|
|
|
|
|
|
||
a св |
a g |
|
|
|
|||||||||||||
A |
|
|
|
|
|
|
|
ax |
|
11 |
|||
|
a |
|
|
|
|
|
y |
|
21 |
||
|
|
|
|||
a |
|
|
|
31 |
|
|
|
z |
|
|
где
12 |
13 |
|
|
axg |
|
|
|||
22 |
|
|
|
|
|
|
|
, |
|
23 |
ayg |
||||||||
|
32 |
|
33 |
|
|
a |
zg |
|
|
|
|
|
|
|
|
|
11 cos12 cos13 cos
x?; xg , 21 |
cos y?; xg , 31 |
cos z?; xg |
x ?; yg , 22 |
cos y?; yg , 32 cos z?; yg |
|
x?; zg , 23 |
cos y?; zg , 33 |
cos z?; zg |
Матрица перехода от одной системы координат к другой определяется перемножением матриц элементарных поворотов, взятых в последовательности, противоположной последовательности этих поворотов (рис.9).
Рис.9. Переход от одной системы координат к другой
Совместим связанную систему координат с нормальной системой координат. Поворот системы координат на один угол называется элементарным поворотом.
Осуществляется последовательность элементарных поворотов: первый поворот - на угол ψ, второй поворот - на угол υ, третий поворот - на угол γ.
AXXg A A A - матрица перехода
A |
cos |
0 |
sin |
|
|||
|
0 |
1 |
0 |
|
|
- поворот на угол ψ |
|
|
|
|
|||||
|
sin |
0 |
cos |
|
|||
|
|
|
|
|
|
|
|
|
cos |
sin |
0 |
|
|
||
A |
|
|
cos |
|
- поворот на угол υ |
||
sin |
0 |
||||||
|
|
0 |
|
0 |
|
|
|
|
|
|
1 |
|
|
||
A |
1 |
0 |
|
0 |
|
|
|
0 |
cos |
sin |
- поворот на угол γ |
||||
|
|
|
|
|
|
|
|
|
|
sin |
|
|
|
|
|
|
0 |
cos |
|
|
A |
A |
1 |
0 |
0 |
cos |
sin |
0 |
|
|
cos |
|
sin |
0 |
|
|||
0 |
cos |
sin |
sin |
cos |
0 |
sin cos |
cos cos |
sin |
|
||||||||
|
|
|
|
|
sin |
|
|
0 |
|
|
|
sin sin |
cos sin |
|
|
||
|
|
|
|
0 |
cos 0 |
1 |
|
|
cos |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Результирующая матрица перехода принимает вид |
|
|
|
|
|||||||||||||
A |
|
|
|
|
cos |
|
sin |
|
0 |
|
cos |
0 |
sin |
|
|
||
XXg |
|
|
|
|
cos cos |
sin |
|
0 |
1 |
0 |
|
|
|
||||
|
sin cos |
|
|
|
|
||||||||||||
|
|
|
|
sin sin |
cos sin |
|
|
|
|
0 |
cos |
|
|
|
|||
|
|
|
|
cos sin |
|
|
|
|
|
cos cos |
sin |
sin cos |
|
|
|
sin cos cos |
cos cos |
sin cos sin sin cos |
|
|
|
||||
|
|
|
cos sin |
|
|
|
sin sin cos sin cos |
sin sin sin cos cos |
1.3 Уравнения движения летательного аппарата
Предположим, что S – твердая, не деформируемая оболочка (рис.10). Уравнение движения ЛА (системы переменного состава) записывается как уравнение движения твердого тела, в которое входит масса затвердевшего тела, если к силам, действующим на летательный аппарат, добавить вариационные силы и Кориолисовы силы инерции и момент от этой силы, чисто реактивную силу и момент.
Рис.10. Тело переменного состав
dQ Fe Fp Fвар Fкор
dt
dK0 M e MFp Mвар Mкор dt
P Fp Fст.д.
Уравнения движения ЛА в инерциальной системе координат имеют вид
|
|
dVа |
e |
|
|
|||
m |
|
|
|
F |
P |
|||
dt |
||||||||
|
e |
|
|
. |
||||
dK0 |
|
|
|
|||||
|
|
|
M |
Mp |
||||
dt |
||||||||
|
|
|
|
|
При составлении системы уравнений движений для неинерциальной системы координат, добавляются переносная и кориолисова силы инерции, вызванные вращением системы отсчета.
|
dV |
|
e |
|
е |
|
кор |
|
||||
m |
|
|
|
F |
P Fин |
Fин |
|
|
||||
|
|
|
||||||||||
|
dt |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
||||
dK0 |
|
e |
|
е |
|
кор |
||||||
|
|
M |
|
M p M |
ин |
M |
ин |
|
||||
dt |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
Проекции вектора на оси произвольной системы координат записываются в виде:
da |
|
|
d |
|
|
|
|
|
|
dax |
|
day |
|
daz |
|
|
di |
|
dj |
|
dk |
|
da' |
|
|
||||||||
|
|
|
|
|
|
|
axi ay j azk |
|
|
i |
|
j |
|
k |
ax |
|
ay |
|
az |
|
|
|
a |
||||||||||
dt |
|
|
dt |
dt |
dt |
dt |
dt |
dt |
dt |
dt |
|||||||||||||||||||||||
di |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
da |
|
da |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
dj |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
j |
|
|
|
|
|
a |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
dt |
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
dk |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Тогда для вектора скорости и вектора кинетического момента K можно записать векторные уравнения в подвижной системе координат
|
dV |
|
|
|
|
|||
|
|
d V |
|
|||||
|
|
|
|
|
|
|
V |
|
|
|
|
|
dt |
||||
|
dt |
|
|
|
||||
|
|
|
|
|
||||
dK |
d K |
|
||||||
|
|
|
|
|
|
|
K0 |
|
|
|
|
|
dt |
||||
|
dt |
|
|
|
Дифференциальные уравнения движения центра масс в проекциях на оси произвольной подвижной системы отсчета примут вид:
dV |
x |
|
|
|
|
|
||
m |
|
|
yVz |
zVy |
|
Fix |
||
dt |
|
|
||||||
|
|
|
|
|
|
|
|
|
dV |
y |
|
|
|
|
|
||
|
|
|
|
|
|
|||
|
|
|
|
|
|
Fiy |
||
m |
dt |
|
zVx xVz |
|||||
|
|
|
|
|
|
|||
dV |
z |
|
|
|
|
|
|
|
m |
|
|
|
xVy |
yVx |
Fiz |
||
dt |
|
|
|
|||||
|
|
|
|
|
|
|
|
где x , y , z - угловая скорость подвижной системы координат.
Дифференциальные уравнения движения вокруг центра масс в проекциях на главные центральные оси инерции запишутся в виде:
K0 Ix xiгл Iy y jгл Iz zkгл
Ix |
d x |
|
Iz Iy y z |
Mix |
|
dt |
|||||
|
|
|
|
||
Iy |
d y |
|
Ix Iz x z |
Miy |
|
dt |
|
||||
|
|
|
|
Iz d z Iy Ix y x Miz
dt
X, Y, Z – главные центральные оси инерции.
Для осесимметричных или имеющих плоскость симметрии ЛА главные центральные оси совпадают со связанными осями.
Векторные уравнения движение центра масс ЛА относительно неинерциальной гринвичской системы координат имеют вид
|
dV |
|
|
|
е |
кор |
|
m |
|
P Ra |
mg Fин |
Fин |
Fi |
||
|
|||||||
|
dt |
|
|
|
|
|
|
|
|
|
сила тяжести mg |
|
|
Запишем уравнения движения ЛА в проекциях на оси траекторной системы координат (рис.11)
Рис.12. Траекторная система координат
dV |
|
d'V |
|
|
|
|
|
|
V |
||
dt |
dt |
||||
|
|
|
Спроектировав векторное уравнение движения, получим
dV |
|
|
|
|
|
|
|
|
Fixk |
||
m |
xk |
|
ykVzk |
zkVyk |
|
||||||
dt |
|||||||||||
|
|
|
|
|
|
|
|||||
dV |
yk |
|
|
|
|
|
Fiyk |
||||
|
|
zkVxk |
|
|
|
||||||
|
|
|
|
|
|
||||||
m |
dt |
xkVzk |
|||||||||
|
|
|
|
|
|
|
|||||
dV |
|
|
|
|
|
|
|
|
Fizk |
||
m |
zk |
|
xkVyk |
ykVxk |
|
||||||
dt |
|||||||||||
|
|
|
|
|
|
|
|||||
Так как Vxk V , Vyk Vzk 0 , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
m |
dV |
Fixk |
|
|||
|
|
|
|
|
|
|
|||||
|
|
|
|
|
dt |
||||||
|
|
|
|
|
|
|
Fiyk |
|
|||
|
|
|
|
|
mV zk |
||||||
|
то |
|
|||||||||
|
mV yk |
|
|
|
|||||||
|
|
|
|
|
|
|
|||||
|
|
|
|
|
Fizk |
||||||
|
|
|
|
|
|
|
|
|
|
|
Составим кинематические уравнения (рис.12).
д ; кд k ; k д хд илиk k
Рис.11. Нормальная система координат
Для определения проекций угловых скоростей , воспользуемся матричным методом пересчёта
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
xk |
|
xk xд |
|
|
cos |
|
|
|
|
|
A |
|
|
|||||
|
|
|
|
|
|
sin |
; |
|||
|
|
yk |
|
|
|
0 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
zk |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|||
|
xk |
|
|
|
|
|
|
||
|
|
|
|
|
|||||
|
|
Axk xд |
|
|
0 |
|
|||
|
yk |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
||
|
|
zk |
|
|
|
|
|
|
|
где
|
|
|
cos k |
sin cos cos sin |
|
|
|
|
|
cos |
|
|
sin k cos |
|||||
|
|
|
|
|
|
sin k |
|
|
||||||||||
|
yk |
|
|
|
yk |
|
|
zk |
, |
|
|
zk |
|
|
, |
yk |
, |
|
|
|
|
|
|
0, |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
, |
|
|
, |
|
zk |
sin k |
|
|
|
k yk к cos
k zk 0
yk cos k sin cos cos sin sin k sin k cos , zk sin k cos cos k
После преобразований записываем кинематические уравнения движения
|
V |
V cos cos |
k |
|
|
xд |
|
|
|
||
|
r |
|
|||
|
r |
|
|
V V cos sin
zд k
|
|
|
|
|
|
|
|
|
|
|
|
|
rcos |
rcos |
|
||||
Подставляя выражения |
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||
yk |
|
V cos sin k |
cos k sin cos |
V cos sin k |
cos sin |
V cos sin k |
sin k sin k cos |
||||||||||||
|
|
|
|
|
|||||||||||||||
|
|
|
rcos |
|
|
|
|
|
|
|
rcos |
|
|
rcos |
|||||
|
|
|
|
V cos2 sin |
k |
tg |
, |
|
|
|
|
|
|
|
|||||
yk |
k cos |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
r |
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
V cos |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
zk |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
r |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
записываются динамические уравнения движения |
|
|
|
||||||||||||||||
|
|
mV Fixk |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
V2 cos |
|
|
|
|
|
|
|
|
|||
|
|
mV Fiyk |
m |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
r |
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
mV cos k |
Fizk |
|
mV2 cos2 |
sin ktg |
|
|
|
|
|||||||||
|
|
|
|
|
r |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Выполняем проектирование внешних сил, действующих на ЛА на оси траекторной системы координат
Fxk |
|
Axk x |
P |
Axkxa |
xa |
Axkxд |
|
0 |
|
|
Fинкор.xk. |
|||||||||
|
|
F |
|
|
|
0 |
|
|
y |
a |
|
mg |
|
|
Fкор. |
|||||
|
|
yk |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ин.yk |
|||
|
Fzk |
|
|
|
0 |
|
|
|
za |
|
|
|
0 |
|
|
|
кор. |
|||
|
|
|
|
|
|
|
|
|
|
|
|
Fин.zk |
|
|
|
|
|
|
2m |
ik |
jk |
kk |
|
|
|
|
кор |
3xk |
3yk |
3zk |
||||||
|
|
||||||||||
Fин. |
mWk 2m 3 |
V |
2mV 3zk j 2mV 3ykVk |
||||||||
|
|
|
|
|
|
|
V |
0 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Окончательно уравнения движения в траекторной системе координат примут вид - динамические уравнения движения
mV Pcos cos Xa mgsin
|
|
|
|
mV2 cos |
|
|
mV P cos sin sin a sin cos a Ya cos a Za sin a |
mgcos 2mV 3zk |
|
|
3 |
||
r |
||||||
mV cos k |
P cos sin cos a sin sin a Ya sin a |
Za cos a 2mV 3yk |
|
|
||
|
|
m V2 sin k cos2 tan r
- кинематические уравнения движения.
|
V cos cos k |
, |
V cos sin k |
H |
V sin |
|
|
||||||
r |
r |
|||||
|
|
, |
|
Вектором перегрузки n называется отношение суммы внешних сил, исключая силу тяжести (сил тяги и аэродинамической силы), к силе тяжести
n P Ra -перегрузка mg
Проекции вектор перегрузки - на оси траекторной системы координат
nxk |
|
Pcos cos Xa |
|
|
|
mg |
|
|
|||
nyk |
|
|
P cos sin sin a |
sin cos a Ya cos a |
Za sin a |
|
mg |
|
|||
|
|
|
|
|
nzk P cos sin cos a sin sin a Ya sin a Za cos a mg
- на оси скоростной системы координат:
nxa Pcos cos Xa mg
nya Psin Ya mg
nza Pcos sin Za mg
- на е оси связанной системы координат:
nx |
A |
|
|
nxa |
||
|
|
ХХк |
|
|
||
ny |
|
nya |
||||
n |
|
|
|
|
n |
|
|
z |
|
|
|
|
za |
Уравнения движения в перегрузках (в безразмерной форме) могут быть записаны в виде
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
V nxk |
sin |
|
|
|
|
||||
|
|
|
|
|
g |
|
|
|
|
|
|
|
|
|
|
|
V |
|
|
V |
|
|
V2 cos |
||||||||
|
|
|
nyk |
cos 2 |
|
|
zk |
|
|
||||||
|
|
g |
g |
gr |
|
||||||||||
|
V cos |
|
k nzk |
2V |
yk |
V2 sin k cos2 tg |
|||||||||
g |
|
g |
|
|
gr |
При системном анализе космических миссий моделирование движения является главным средством ответа на вопрос об их возможности и реализуемости.
2. Невозмущенное движение космического аппарата
Рассматривается движение летательного аппарата как материальной точки в поле тяготения одного небесного тела (Земли), создающего центральное поле притяжения.
Принимаются следующие допущения: