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

книги из ГПНТБ / Селезнев В.П. Инерциальная навигация летательных аппаратов

.pdf
Скачиваний:
24
Добавлен:
29.10.2023
Размер:
6.91 Mб
Скачать

а также зависимостей (1.13), (1.14), (1.15), (1.16). После алгебраи­ ческих преобразований получим

Р +

<“ о Р == 8 ( " f 'a +

а) +

8 (<“ ув +

я)

 

 

 

Уа

«) +

Л

 

+

 

 

 

2 R

 

*

 

О

 

^

I

/

 

I

*

\

 

£ ха

-

а

ё 2а

(1.49)

 

+ 2“ .гв - д

а -г ® (®.ve

+

а) 0>*а

-

+

 

;

 

5 +

шо 8 = —

fj (<оУа 4-

а) — ^ (u>j

+

 

а) +

шг<г (шУа +

а) +

 

 

 

0 R

 

п R

 

,

,

 

,

\

 

ё*а , ё*а

(1.50)

 

+ 2 Т Г

 

2 —

Шга а +

а (шУа -I-

а) <0 га +

Я

 

R

 

Я

 

 

 

Л

 

 

 

 

 

 

 

 

 

 

 

Уравнения инерциальной вертикали с акселерометрами,

ориентированными

по

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

системе координат, полу­

чаются из уравнений

(1.49)

и (1.50)

путем замены индексов х 0,

Уа>

Za Н3 Ъ Ъ t-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р + шо Р

-

8 ( « f

+

я)

- f 3 (со,, - f a )

(Oj ( u f - f

a)

- f

2 —

mc

- f

 

 

 

+

2

D

 

 

 

 

 

 

PV

 

 

 

 

^ -51)

 

 

 

 

 

+ * («f + «)<ос —

 

 

 

 

 

8 - f <oo p =

. p (со,, - f

a) —

p ( u f - f

a)

+

шс ( " f

+

a) +

2

a>£ —

 

 

 

-

2 R_

 

 

 

 

a) «4

gt

 

 

 

 

 

 

 

 

 

 

R

'4* -Ь “ (ШТ) +

 

+ 'Я a’

 

 

 

 

где u>{, a>4, a>(. и gz определяются по формулам (1.6) и (1.11). Свободная в азимуте ориентировка акселерометров, соответст­

вующая угловой скорости «>,, + a = 0, дает следующие уравнения вер­ тикали:

р + “ о Р — 2 J - - Шга -)- 2я — Юд.д -|- — g x a );

8 + 8 = 2 — co,e — 2a ^ шгй + -I- (g ta + g Xa a).

§ 1.7. МЕТОДИЧЕСКИЕ ПОГРЕШНОСТИ ИНЕРЦИАЛЬНОЙ ВЕРТИКАЛИ

Если направление геоцентрической вертикали принять за истин­ ную вертикаль, то погрешности инерциальной вертикали могут быть найдены как частные решения уравнений (1-49), (1.50), (1.51) или (1.52). Рассмотрим влияние отдельных параметров движения летательного аппарата и особенностей гравитационного поля Зем­ ли на отклонения инерциальной вертикали. Для упрощения анали­ за сделаем допущение, что связь между уравнениями каналов вертикали несущественная и ею можно пренебречь.

44

а)

Погрешность

от

нерадиальности гравитационного поля

Земли.

 

и (1.50) следует, что операционное изобра­

Из уравнений (1.49)

жение погрешности:

 

 

 

й. = - ^ -----8

 

gz- + l_gJ ±_.

(1.53)

Pl

R i p t + vl)

R ( P 2 + » D

 

Ускорения gxa и g Z(l меняются медленно (в горизонтальном полете) по сравнению с переходными процессами движения вер­ тикали. Установившееся значение погрешностей:

аё*д ~~ ёха

ёга +

*ёха

P l уст

у ст

(1.54)

ёт

ёт

Пренебрегая влиянием азимутального отклонения акселеромет ров от навигационной системы координат (принимаем а =*0), най­ дем установившееся результирующее отклонение вертикали

 

1/ й „ +

<

:Дср2

(1.55)

+1 = К р ? уст

Or

 

1 уст

 

 

 

 

ёт

 

 

 

где Д'-fi, определяется выражением (1.10).

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

шается (ускорение g T от силы тяжести,

что приводит к росту погреш­

ностей о, и fJj.

 

 

ускорений,

вызванных верти­

б)

Погрешности от поворотных

кальной

угловой скоростью

<аУа + а

и горизонтальной скоростью

полета:

 

 

 

 

 

 

 

(°>va + a)(cw>ze — “>ха)

 

+

«) {тга + лшха)-

Р2 =

2

 

 

„о

| 2

 

(1.56)

 

Рг Щ

 

 

р

+ “о

 

 

Результирующее угловое

отклонение

инерциальной

вертикали

 

фг “ V

^2 + &2

 

 

 

(1.57)

В частном случае, когда

летательный

аппарат

неподвижен,

<°уа — шзsin ф, u)r = io3costp установившееся

отклонение

вертика­

ли равно

ла sin фcos <p

 

 

 

-

 

 

-------- 1------—. Максимальное отклонение будет на широте

щ

2

45°, равное —®- = 17-10-3 или 5',9. 2ш?

45

Суммарное установившееся отклонение инерциальной вертика­ ли, вызванное нерадиальностью гравитационного поля Земли и пово­

ротным

ускорением для неподвижного

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

аппарата,

равно (углу отклонения линии отвеса относительно

геоцентрической

вертикали. Следовательно, инерциальная вертикаль с акселеромет­

рами, неподвижными в азимуте (а = 0)

относительно

Земли, при

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

положение, совпадающее с линией отвеса.

 

 

 

Большие скорости полета, как это показано

для одномерной

инерциальной вертикали, могут вызвать погрешности, достигаю­

щие десятков угловых минут.

 

 

 

в)

Погрешности от поворотных ускорений,

вызванных вертикаль­

ной скоростью R и горизонтальными угловыми скоростями вращения

 

К а + «Ч, а)

<*)2

а)

 

Эз

£а

'

(1.58)

Р2 + Шо

Р2+ “о

 

 

 

 

Результирующее отклонение инерциальной вертикали

 

--=Vf3 + о2

 

 

(1.59)

Характерно, что вертикальная скорость R у летательных аппа­ ратов, имеющих ограниченную высоту полета, меняется в незначи­ тельных пределах. Существенные погрешности могут возникнуть при взлете или спуске летательного аппарата,-

Для представления о порядке погрешности рассмотрим пример: летательный аппарат совершает подъем на высоту 30 км в течение 20 мин. с постоянной скоростью.

Полагая u>r =

«>*я + (0*o

также постоянной и равной, напри­

мер шэ= 7 ,28 ■10

51/сек,получим следующее значение угла отклоне­

ния:

 

2 /?«ч (cos <»0 t — 1j

 

 

 

(1.60)

 

 

g r

 

 

Поскольку

время

подъема

составляет —- периода

колебаний

вертикали (Т ~

 

 

4

 

84,4 мин.), то можно принять cosm0/=:0. При этих

условиях получаем

= 3,6-10~4 или 1',2.

 

Периодические изменения

высоты полета происходят обычно

с частотной, намного превышающей частоту колебаний инерциаль­

ной вертикали. Возникающие при этом

погрешности

вертикали

обычно пренебрежимо малы.

 

 

 

Из рассмотренных методических погрешностей наиболее суще­

ственное значение имеют погрешности от поворотных ускорений

,32,

82 и от нерадиальности гравитационного

поля Земли

pj и

8j.

У инерциальных систем высокой точности принимаются специальные

меры по компенсации этих погрешностей. В частности, одним из способов компенсации погрешностей от поворотных ускорений ({32, о2) является применение свободной в азимуте платформы о акселе­ рометрами. Погрешности от поворотных ускорений, вызванных вер­ тикальной скоростью, компенсируются у инерциальных систем, при­ меняемых на ракетах с большой высотой подъема.

§1.8. ИЗМЕРЕНИЕ ВЫСОТЫ ПОЛЕТА

ВТРЕХМЕРНОЙ ИНЕРЦИАЛЬНОЙ СИСТЕМЕ НАВИГАЦИИ

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

Фиг. 1.19. Структурная схема канала, намеряющего высоту

полета

фиг. 1.19. В этой схеме имеется обратная связь, позволяющая ис­ ключить влияние ускорения от сил тяготения на измерение рас­ стояния R от центра Земли до летательного аппарата. Из структур­ ной схемы получаем следующее уравнение системы:

( М у Як ~ ■*«„) — + R q

До Д„

(1.61)

Р

 

 

где Яу = R + Я + 4я„ — ускорение, измеряемое вертикальным ак­ селерометром;

ЯК:

к.

— сигнал, компенсирующий ускорение от

 

сил тяготения; ДаК— сигнал, компенсирующий поворотные

ускорения; к2, к3 — передаточные коэффициенты интеграто­

ров;

k\ — передаточный коэффициент акселеро­ метра;

Ro и R0 — начальные значения R и R. Передаточные коэффициенты подбираем из условия:

к\ = 1, к2 = кз = 1, кх р',

где р' — постоянная тяготения.

-17

Обозначим погрешности:

Д/? = /?ИЗм— R> b g = g ~ g K, Да = &ау —Дак.

Уравнение (1.61) преобразуем к виду:

 

В этом уравнении

д £ =

A g

-f Да.

(1.62)

 

 

 

 

. JL

^

* - 2 — Д / ? - °2-^8-AR.

R2

(R + ДИ)2

R3

R

Подставляя значение Д^ в уравнение (1.62), получим

 

Д/? — 2 —

ДЯ = Да.

(1.63)

 

 

R

 

 

Если предположить, что погрешность компенсации поворотных ускорений является постоянной величиной, то изменение погрешно­ сти в измерении высоты при начальных нулевых условиях будет про­ исходить по закону:

ДR = -----—-----(ch 1 f

2 JL t — 1

(1.64)

/? Т

1

 

"

__

 

При t = О гиперболический

косинус ch

/

2 - С - ^ = 1

и ДЯ = 0.

В дальнейшем, с течением времени ДR

 

н

 

возрастает по закону гипер­

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

§ 1.9. МЕТОДИЧЕСКИЕ ПОГРЕШНОСТИ ГОРИЗОНТАЛЬНЫХ ИНЕРЦИАЛЬНЫХ СИСТЕМ НАВИГАЦИИ

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

Полагаем, что курсовая система стабилизирует акселерометры в азимуте с погрешностью, равной углу а . Инерциальная вертикаль имеет погрешности: в вертикальных плоскостях ух и уг соответствен­ но углы р и 8, Измеряемые акселерометрами ускорения опреде­ ляются выражениями (1.42).

Погрешности ориентатора в географической системе координат

Полагаем, что начальные координ^ы местонахождения зада­ ются в инерциальную систему с погрешностями Д<р и ДX, учет ско­ рости вращения Земли происходит с погрешностями Ди>3 и ускоре­

48

ния измеряются с погрешностями: в северном направлении &aN=

= р£т +

Да, и восточном направлении ДаЕ 8gT-f- Да2,

где

 

Да, = а пЕ + g t - а с а;

 

 

Да2 =

ani +

а^ а.

 

Из (уравнения (1.38) можно получить следующие уравнения по­

грешностей:

 

t t

 

 

 

 

 

 

Лш = Д«р0

 

ДаN (dt)2;

(1.65)

 

 

 

Я R

 

 

 

 

 

*Л? sin у W t: dt.

 

 

 

 

I R cos:tp

 

Рассмотрим характерные особенности отдельных погрешностей,

а)

Погрешности при введении начальных координат:

 

 

Д<р, = Дср0,

 

ДХ, = ДХ0.

(1.66)

б)

Погрешности от неточного учета угловой скорости враще­

ния Земли:

 

 

Jt &i»3dt.

 

 

Д^>2 - Oj

ДХ2 =

(1-67)

 

 

 

 

и

 

Для оценки величины этой погрешности допустим, что

Д«>3 =

=const и равняется 0,01 и>3.

Втечение времени t = 3600 сек. накапливается погрешность

ДХ2 = 2,62 • 10~3 =

9'.

 

вертикали:

,

в) Погрешности от инерциальной

 

Л(Рз

Р +

Да,) (dt)2;

 

 

о Sо S

>

 

 

 

ДХ3= J

sec ср j*

( - g r o -f- Да2)

(dt)2.

(1.68)'

 

И

 

 

 

 

В свою очередь угловые отклонения р и 8

вертикали могут быть

записаны в виде операторных изображений, которые получим из уравнений (1.45), или упрощенных структурных схем

— Да,

R

Да,

R

(1.69)

 

 

р 2+ & R

4 В, П. Селезнев

49

Ф иг. 1.20. Схема связи по грешности в измерении дол­ готы от погрешности верти кали

Подставляя значения углов (1.69) в (1.68) и записывая уравне­ ния в операторном виде, получим

 

 

1

baigr

-I-Да,

 

 

Rp'-

Rp2

 

 

 

r (p 2 + g

 

 

 

 

 

R

 

 

 

Да1

= - Р ;

 

(1.70)

 

R \p 2 + ^

 

 

 

 

 

 

 

R

 

 

 

 

AU

(Да2 — g rS)gectP

&

 

 

Rp2

 

COS <P

 

 

 

 

 

 

Из формул (1.70) следует,

что

погрешность в измерении

широты

численно равна углу отклонения инерциальной вертикали от геоцент­ рической вертикали в плоскости геогра­ фического меридиана. Погрешность в из­ мерении долготы численно равна углу по­ ворота инерциальной вертикали вокруг оси мира.

В подтверждение этого рассмотрим отклонение (угол о ) инерциальной вер­ тикали на сфере (фиг. 1.20). Из равно­

бедренного

треугольника PN АВ следует,

sin ДХ/2

1

что --------------------- или для малых

sin 8/2

cos®

%

углов получаем Д

cos <р

Для снижения уровня погрешностей Д<Рз и Да., необходимо уменьшить откло­ нения инерциальной вертикали, вызван­ ные начальными условиями и неполной

компенсацией ускорений Да, и Да2.

г) Погрешности в измерении долготы, вследствие воздействи канала, измеряющего широту:

 

 

 

 

sin ср

 

WEdt.

(1-71)

 

 

д >-4 =

J* А<? R cos2 <р

В полярных районах,

где

cos <р

0 , погрешность А)ч может не­

ограниченно возрастать.

 

 

 

 

 

д) Погрешности от курсовой системы.

 

Поворот

акселерометров

па угол

а относительно географиче­

ской

системы

координат

приводит

к

появлению

погрешностей

Д-fs и

ДХ8.

 

 

 

 

 

 

5 0

Учитывая, что в данном случае

 

Лад,

ас а

и ДаЕ --- а,

получим

t t

 

 

 

 

 

 

 

 

 

 

 

A'f

а с а dt;

Д>..

j

s e c - i / j

 

(1.72)

ояо-

 

 

 

 

 

 

(I

 

/

Так, например, при отклонении курсовой системы на 1° отно­ сительная погрешность в измерении координат местонахождения достигает 1,75%. Следовательно, требования к допустимым по­ грешностям системы должны быть весьма высокие.

Погрешности ориентатора в ортодромической системе координат

Из уравнений (1.39) и (1.40) получаем следующие значения по­ грешности в измерении ортодромических координат:

i

 

 

 

Ь Х = \

ЬУКХ

dt\

(1.73)

о

 

/

 

Лаг dt —ДНК,

dt,

 

где Дах и Да, определяются из формул

(1.42).

 

Если используется второй способ компенсации переносной ско­ рости вращения Земли, то компенсационные сигналы должны быть равны:

VKX= /?<»j — з cos s,

VKZ= /?ш, cos и = /?о)3 sine cos и,

где г — угол наклона круга

ортодромии относительно экватора;

и— угол между вертикалью места и линией пересечения плос­ кости экватора и круга артодромии (аргумент широты).

Учитывая указанные соотношения, получим

ДКK.V — Д^s R cos г —Де/?(0я Sin е;

д ^ к,

 

 

(1.74)

д<°3 R sin е cos и 4- Д£ш8 R cos £ cos и Д« “ з R sin г sin и,

где

Дг — ошибка в задании угла наклона

круга ортодро­

 

 

мии относительно экватора;

 

Дм=Д«0-|-—-

— ошибка в измерении аргумента широты;

 

Дм0

— ошибка в установке начального

аргумента ши­

 

 

роты.

 

а) Погрешности от инерциальной вертикали:

 

 

 

ДЛ ^-РЯ ; дZ, = S#.

(1.75)

б)

Погрешности от курсовой системы:

 

t t

длг2 = \.Г V ia*(dtfr

о о

■t

(1.76)

AZ2= \ j v ra*(dtf.

и о

 

А

51

При

а == const погрешности равны:

 

 

 

\ Х 2 =

— Zo.\ AZ2=

Л’а,

(1.77)

t

t

i

t

 

 

где A = j

JVxa{dt)2\ Z =

\

\Vza {dt)2 — составляющие истинного

o o

o o

 

 

пройденного пути.

 

 

 

 

в) Погрешности от неточного введения угловой скорости вра­

щения Земли:

 

 

 

 

 

t

 

 

t

 

ДА', — — j Дс1>3 R cos е dt\

AZS =

J Аш3/? sine cos udt.

(1.78)

 

о

 

6

 

г) Погрешности от неточного задания угла наклона круга орто­ дромии относительно экватора:

t

>

 

A X i = \As/?«o3 sin е d t \

ДZ i — ■— j" AeRio3 COS e COS ll d t .

(1-79)

0

U

 

д) Погрешность в измерении бокового смещения от заданной ортодромии, возникающая вследствие взаимной связи измеритель­ ных каналов:

AZr,=

А Х со3/? sine sin udt.

(1.80)

 

~R

 

Погрешность Azb (как и погрешность Д>.ч) , является следст­ вием перекрестной связи между каналами, измеряющими X и Z, необходимой для компенсации вращения Земли.

§ 1.10. ДЕМПФИРОВАНИЕ КОЛЕБАНИЙ ИНЕРЦИАЛЬНОЙ ВЕРТИКАЛИ

Отсутствие демпфирования колебаний инерциальной вертикали является весьма существенным ее недостатком. Вследствие этого возникают жесткие требования к точности начальной установки платформы с акселерометрами и задания начальной абсолютной скорости в первые интеграторы. Так, например, при возникновении в первых интеграторах погрешности в измерении скорости, равной 2,3 м/сек, возникают незатухающие колебания вертикали с ампли­

тудой в 1'.

Демпфирование колебаний можно осуществить с помощью пря­ мых и обратных связей и фильтров в схеме инерциального ориентатора.

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

52

рактеристику, уменьшающую воздействие на

вертикаль

внешних

периодических возмущений.

 

 

Рассмотрим простейшие связи и фильтры, применяемые в инер­

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

 

 

Демпфирование с помощью прямой связи,

 

охватывающей один из интеграторов

 

(фиг. 1.21)

 

 

Полагая измеренное ускорение в виде а =

pV — $gr,

где 13-

угол наклона оси чувствительности акселерометра (все другие составляющие ускорения будем считать скомпенсированными), по­ лучим следующее уравнение движения системы в одной плоскости:-

2 + ьх h kt gTp + (*! k2 ks +

kx k3 kt p) gT] p =

 

= p v { k xk2k3 -

kikt kt p' I/.

(181)

Фиг. 1.21. Структурная схема демпфирования прямой связью, охватывающей интегратор

Условие невозмущаемости от ускорений р V возможно при k\k2k% =

— — . Обозначая я>о

^4 Р) ё"т»

^4

R

 

 

нолучим

 

 

(p*+2rlo>0p + <»l)V =

^ P 2V.

(1.82)

 

g Т

 

Собственное движение вертикали

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

однородным

уравнением:

 

 

Р ~г 2?)ш0 ^ -J- сор р —0,

(1.83)

где и>0 в общем виде может быть переменной величиной. В горизон­ тальном полете pV — 0 и решение уравнения (1.82) при начальных

условиях (при t — 0, р = ро и Р = Ро) »меет вид:

р =

(A sin о>0 / 1 — r(2 t + В cos о)0 V 1 - т]2 t),

(1.84)

53