Сихарулидзе Ю.Г. Баллистика и наведение летательных аппаратов (2-е изд., 2013) (1246775), страница 13
Текст из файла (страница 13)
Уравнения движенияускорения позволяют существенно упростить уравнения движения и задачу ихинтегрирования.Для получения уравнений движения в требуемой системе координат надопроектировать на ее оси векторные уравнения движения центра масс ЛА (1.1.5)и движения относительно центра масс (1.1.12). Уравнение (1.1.5) описывает движение центра масс с учетом его перемещения по отношению к корпусу ЛА, однакопрактический интерес обычно представляет именно движение самого корпуса,рассматриваемого как материальная точка. Заметим, что при полете на активномучастке корпус ЛА перемещается на сотни километров и приобретает скоростьнесколько километров в секунду, в то время как сам центр масс перемещаетсяотносительно корпуса на несколько метров со скоростью, равной несколькимсантиметрам в секунду.
Поэтому во многих задачах баллистики и некоторыхзадачах динамики можно отождествлять движение центра масс ЛА с движениемкорпуса и описывать его уравнениемc = mWF + P.(1.8.1) — главный вектор c — вектор абсолютного ускорения центра масс ЛА, Fгде Wвсех внешних сил, приложенных к ЛА, P — главный вектор реактивных сил (тягадвигателя).1.8.1. Движение центра масс. Получим уравнения движения в начальной стартовой (инерциальной) системе координат. Предположим, что система управленияЛА имеет ГСП с тремя акселерометрами, ориентированными по взаимно перпендикулярным направлениям, и БЦВМ для решения задач навигации, управленияи стабилизации.Если пренебречь уходом осей гироскопов на активном участке, то можно считать, что ГСП определяет направление осей начальной стартовой (инерциальной)системы координат (см.
п. 1.2.2). Проектируя уравнения движения (1.8.1) на осиэтой системы координат, получим(ph)Wx ln = Wx ln + gx ln ,(ph)Wy ln = Wy ln + gy ln ,где(ph)Wz ln = Wz ln + gz ln ,(1.8.2)1(ctr)Px ln + Xln + Xln,m1(ph)(ctr)Py ln + Yln + YlnWy ln =,(1.8.3)m1(ph)(ctr)Pz ln + Zln + ZlnWz ln =m— составляющие кажущегося ускорения; gx ln , gy ln , gz ln – составляющие ускорениясилы притяжения; Px ln , Py ln , Pz ln — составляющие тяги маршевых двигателей;(ctr)(ctr)(ctr)Xln , Yln , Zln — составляющие полной аэродинамической силы; Xln , Yln , Zln —составляющие тяги управляющих двигателей.Акселерометры, ориентированные по осям начальной стартовой системы координат, будут в процессе полета измерять составляющие кажущегося ускорения(ph)Wx ln =.Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»1.8. Уравнения движения в начальной стартовой (инерциальной) системе координат(ph)(ph)53(ph)Wx ln , Wy ln , Wz ln . Если измеряемые величины вводить в БЦВМ как правыечасти уравнений движения (1.8.2) и добавлять к ним составляющие ускорениясилы притяжения gx ln , gx ln , gx ln , то, как уже отмечалось, можно интегрироватьуравнения движения на борту ЛА.
Это позволяет решать навигационную задачупо определению текущего вектора состояния ЛА (включающего его координатыи составляющие вектора скорости), что необходимо для выбора потребного управления на оставшейся части траектории. Начальные условия интегрирования могутбыть нулевыми, например, когда ЛА начинает движение со стартового устройства,или отличаться от нуля, например, при спуске ЛА с орбиты на поверхность Земли.1.8.2. Составляющие гравитационного ускорения.
Наряду с начальной стартовой системой координат 00i xln yln zln , для дальнейших расчетов, связанных с определением текущего радиуса-вектора r и широты ϕ, удобно ввести земную инерциальную систему координат 00 xE yE zE . Ее начало совпадает с исходной, ось00 yE направлена по радиусу-векторуR0, проведенному из центра общего земногоэллипсоида через точку старта 00 в начальный момент времени t = 0, ось 00 xEнаправлена в плоскости местного меридиана к северному полюсу, а ось 00 zEдополняет систему координат до правой (рис. 1.19).
Чтобы перейти от начальнойстартовой системы координат 00i xln yln zln к земной инерциальной 00 xE yE zE , надовыполнить два поворота: на величину угла азимута A вокруг оси 00i yln , чтобысовместить промежуточную ось 00 x с плоскостью местного меридиана, и наугол Δϕ0 вокруг оси 00i zE , чтобы ось 00i yln совпала с направлением радиусавектора R0 , проведенного из центра общего земного эллипсоида через точкустарта в момент времени t = 0 (рис. 1.20).
Здесь Δϕ0 — разность геодезическойи геоцентрической широт в точке старта. Матрица перехода от начальной стартовойсистемы координат к земной инерциальной имеет вид: cos Δϕ0 cos A sin Δϕ0 − cos Δϕ0 sin A N =sin Δϕ0 sin A − sin Δϕ0 cos A cos Δϕ0.sin A0cos AПересчет составляющих произвольного вектора a из начальной стартовойсистемы координат в земную инерциальную осуществляется по формулеaE = Naln .В частности, по известным составляющим текущего радиуса-вектора ЛА в начальной стартовой системе координат xln , yln , zln можно вычислить его составляющие в земной инерциальной системе координат:xE = xln cos Δϕ0 cos A + yln sin Δϕ0 − zln cos Δϕ0 sin A,yE = −xln sin Δϕ0 cos A + yln cos Δϕ0 + zln sin Δϕ0 sin A,zE = xln sin A + zln cos A.Рассмотрим теперь в земной инерциальной системе координат два единичныхвектора, один из которых направлен по вектору угловой скорости вращения Земли E0 = (cos ϕ0 , sin ϕ0 , 0),ω.Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»54Глава 1.
Уравнения движенияРис. 1.19. Земная инерциальная система координата второй направлен по радиусу-вектору r, проведенному из центра общего земногоэллипсоида в центр масс ЛА,r 0 =xE R + yE zE,,rrr,где ϕ0 — широта точки старта,r=x2E + (R + yE )2 + z2E .Скалярное произведение этих векторов позволяет определить текущую геоцентрическую широтуϕ = arcsin(ωE0 ·r0 )илиϕ = arcsinxER + yEcos ϕ0 +sin ϕ0 .rr.Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»1.8. Уравнения движения в начальной стартовой (инерциальной) системе координат55Рис. 1.20. Переход от начальной стартовой системы координат к земнойРасстояние R (вдоль текущего радиуса-вектораr) от центра до соответствующейточки на поверхности общего земного эллипсоида вычисляется по формуле√1 − e2R = a,1 − e2 cos2 ϕоткуда с точностью до членов первого порядка относительно сжатия имеемR ≈ a(1 − α sin2 ϕ).где a — экваториальный радиус, e — эксцентриситет эллипса в меридиональнойплоскости, α — сжатие общего земного эллипсоида.
В частности, при ϕ = ϕ0получим R = R0 .Знание текущих величин геоцентрического радиуса-вектора r и геоцентрической широты ϕ дает возможность вычислить с приемлемой точностью составляющие вектора ускорения силы притяжения, т. е. гравитационного ускорения,направленные соответственно по радиусу (gr ) и перпендикулярно ему (gϕ ) [1.28]:2 2Rav Ravq 21 − 3 sin ϕ ,1+α−gr = −gavrr222Rav q Ravα−sin 2ϕ,gϕ = −gavr2rαμω 2 a3, gav = 2 , q = E .Rav = a 1 −3RavμПереходя к начальной стартовой системе координат, получим:gxln = gϕ cos Δϕ0 cos A − gr sin Δϕ0 cos A,gyln = gϕ sin Δϕ0 + gr cos Δϕ0 ,gzln = −gϕ cos Δϕ0 sin A + gr sin Δϕ0 sin A..Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»56Глава 1. Уравнения движенияПосле подстановки этих величин в уравнения движения (1.8.2) последниеможно интегрировать. Вычисленная величина r позволяет определить текущуювысоту относительно общего земного эллипсоида (рис.
1.19)h = r − R.Если уравнения (1.8.2) интегрируются не на борту ЛА в процессе полета, а припроведении баллистических расчетов с использованием ЭВМ, то составляющиекажущегося ускорения (1.8.3) следует определять как функции текущих параметровдвижения. Выведем необходимые соотношения.1.8.3. Составляющие тяги двигателей и аэродинамической силы.
Составляющие векторов тяги маршевых двигателей P и тяги управляющих двигателей Fctr ,а также составляющие вектора полной аэродинамической силы Ra естественнозадавать в связанной системе координат. Действительно, маршевые двигателинеподвижны относительно корпуса ЛА. Управляющие двигатели или другие «рули»должны поворачиваться относительно корпуса для создания соответствующих сили моментов, а аэродинамические силы для ЛА типа ракеты всегда определяютсяв виде осевой и нормальной составляющих. Для ЛА самолетного типа, имеющихтолько плоскость симметрии, аэродинамические силы обычно задаются в полусвязанной системе координат, но их легко привести к связанным осям, которыеповернуты на угол атаки α относительно полусвязанных осей.Итак, будем полагать, что известны составляющие векторов тяги маршевыхи управляющих двигателей, а также составляющие полной аэродинамической силыв связанной системе координатP = (Px , Py , Pz ) , Ra = (X , Y , Z) , F ctr = (Xctr , Yctr , Zctr ,) ,заданные соотношениями (1.4.2), (1.6.1), (1.6.3), (1.6.4) и условием ориентациимаршевых двигателей относительно корпуса ЛА.
Используя матрицу перехода L−1от связанной к начальной стартовой системе координат, получим составляющиеуказанных векторов в начальной стартовой системе координат:⎡⎤⎡⎤(ctr)Px ln + Xln + XlnPx + X + Xctr⎢(ctr) ⎥⎣ Py ln + Yln + Yln⎦ = L−1 ⎣ Py + Y + Yctr ⎦(ctr)Pz + Z + ZctrPz ln + Zln + Zlnили подробно(ctr)Px ln + Xln + Xln = (Px + X + Xctr ) cos ϑ cos ψ ++ (Py + Y + Yctr ) (sin γ sin ψ − cos γ sin ϑ cos ψ) ++ (Pz + Z + Zctr ) (cos γ sin ψ + sin γ sin ϑ cos ψ) ,(ctr)Py ln + Yln + Yln = (Px + X + Xctr ) sin ϑ ++ (Py + Y + Yctr ) cos γ cos ϑ − (Pz + Z + Zctr ) sin γ cos ϑ,Pz ln + Zln +(ctr)Zln= − (Px + X + Xctr ) cos ϑ sin ψ ++ (Py + Y + Yctr ) (sin γ cos ψ + cos γ sin ϑ sin ψ) ++ (Pz + Z + Zctr ) (cos γ cos ψ − sin γ sin ϑ sin ψ) .(1.8.4).Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»1.8. Уравнения движения в начальной стартовой (инерциальной) системе координат57Уравнения (1.8.4) после деления на массу m будут определять три компонентыкажущегося ускорения.1.8.4.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
