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

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

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

но стабилизатора с помощью привода 2 на углы, пропорциональные выходным сигналам двойного интегратора 5 во времени. (На фиг. 1.1

обозначено р — — ). dt

Ф иг. 1.1. Упрощенная схема работы одного канала гироориентатора

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

При перемещении системы из положения А в положение В воз­ никают ускорения pV. Двойной интеграл во времени от этих ускоре­ ний равен пройденному пути 5;

р г

р

Ставится задача так подобрать параметры системы, чтобы во время полета положение акселерометра сохранялось горизонталь­ ным. Если передаточное число канала ориентатора, состоящего из последовательно соединенных звеньев — акселерометра, интеграто­ ров и привода, обозначить k, то |условие сохранения горизонтально­ сти акселерометра будет

то есть угол

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

относительно стабилиза­

тора должен

быть равен углу поворота

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

носительно направления вертикали места вылета.

10

Из этого условия и вытекает искомое требование к подбору передаточного числа системы:

к = ±.

R

Если это требование удовлетворено, то во время полета акселе­ рометр сохраняет горизонтальное положение при любых ускорениях летательного аппарата. Следовательно, рассматриваемое устройство приобретает свойства вертикали, инвариантной (невозмущаемой) по отношению к ускорениям движения подвеса.

При отклонении акселерометра от горизонтального положения на угол а (фиг. 1.2,а) в сторону увеличения абсолютного угла поворота происходит следующее.

Ф иг. 1.2. Измерение ускорений при наклонах акселерометра.

Если до наклона акселерометр измерял только горизонтальноеускорение pVr что и обеспечивало поворот акселерометра на угол апер, то в данном случае измеряемое ускорение будет равно:

ах = pV cos a'— gsin а.

Абсолютный угол поворота акселерометра

оа6с = апер + а бу­

дет уменьшаться, поскольку составляющая g sin а

имеет отрицатель^

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

« = -

J R

о

Скорость а достигнет максимума в тот момент, когда угол на­ клона а исчезнет.

В дальнейшем, несмотря на то, что на вход акселерометра не подается сигнал g'sin а , платформа с акселерометром будет продол­ жать движение вверх — в сторону уменьшения угла аа6с. Причи­ ной этого является сигнал первого интегратора, в котором фикси­ руется величина угловой скорости поворота акселерометра. Появ­ ляется наклон акселерометра в другую сторону (фиг. 1.2,6), а составляющая вертикального ускорения g sin а меняет свой знак.

Интегрирование составляющей gsina

приводит к тому, что накоп­

* Ускорение g вызвано подъемной силой

и направлено вверх.

11

ленная в первом интеграторе отрицательная угловая скорость я уменьшается. Увеличение угла наклона акселерометра вверх пре­ кратится, как только в первом интеграторе исчезнет накопленная

 

скорость ж. Поскольку на вход ак­

 

селерометра

поступает положитель­

 

ный сигнал

^sina, то в первом

ин­

 

теграторе

начинается накапливание

 

положительной скорости

а и абсо­

 

лютный угол яабс будет увеличивать­

 

ся, а угол

я — уменьшаться. Далее

 

процесс

повторяется, характеризуя

 

собой

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

колебания

 

платформы

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

На

 

фиг. 1.3 показаны графики переход­

 

ных процессов, возникающих при на­

Ф иг. 1.3. Графики переходных

чальном

отклонении

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

процессов платформы с акселе­

от

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

положения.

 

рометром

ния

Представление о законе движе­

 

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

при

начальных

возмущениях можно получить из рассмотрения процесса его движения.

Угол «абс отработки акселерометра относительно стабилизатора получается в результате интегрирования ускорений и равен:

Лабс = Й*^,

Р2

где

®абс ^пер а.

Полагая угол а наклона акселерометра малым, так что sina —а и cos а = 1, получим следующее значение ускорения:

 

 

ах = pV

g&.

 

 

 

Подставляя значение

ускорения

аг,ша

также учитывая, что

*пер '

, оудем иметь:

 

 

 

 

 

Rp

 

 

 

 

 

р- ж - f

kga = pV ( k

 

( 1. 1 )

 

 

 

7?

 

J _

 

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

pV

 

является к-

R '

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

; + х « =- 0.

( 1.2)

R

12

Если начальные возмущения отсутствуют, то угол а равен нулк> в течение всего времени полета. Для начальных условий: при t = 0V

2 = а() и а — а0 движение акселерометра относительно горизонталь­ ного положения будет представлять гармонические колебания вида:

а = у.пcos V Rg_

-°— Sin 1/

Г

 

R

(1.3)

 

X

 

 

Период этих колебаний

R

 

 

R_

 

 

Т

 

(1.4)

 

 

g

для условий Земли (R — 6370 км и g — 9,81 м/сек?) равен 84,4 мин. и носит название периода М. Шулера. Таким же периодом колеба­ ний обладает физический маятник с длиной подвеса, равной ра­ диусу Земли.

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

акселерометра относительно стабилизатора, имея в

виду, что S

R -* пер-

обнаружить

и

Аналогичные динамические свойства можно

у других систем инерциальной навигации, осуществляющих измере­ ния в горизонтальной системе координат.

§1.3. ИЗМЕРЕНИЕ УСКОРЕНИЙ

ВГОРИЗОНТАЛЬНОЙ СИСТЕМЕ КООРДИНАТ

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

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

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

Пусть летательный аппарат перемещается относительно поверх­ ности земли со скоростью W на высоте Я над уровнем моря.

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

Совместим центр горизонтальной

системы

координат £*]С

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

(фиг. 1.4).

Оси координат

направим следующим образом: ось Ok — на север, ось Ог\— по гео­ центрической вертикали вверх и ось_ОС — по параллели на восток (правая система координат). Вектор W находится в плоскости гори­ зонта ЮС и составляет с направлением меридиана Ok курсовой угол

13

Составляющие относительной скорости по осям координат:

 

северная WN—W cos7; восточная WE— W sin -ц вертикальная R.

 

Вследствие вращения Земли летательный аппарат имеет допол­

нительную переносную скорость, равную ю3/?cos®

и направленную

по

оси ОС Составляющие абсолютной

скорости

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

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

vj, С:

 

 

 

 

 

V t = W N, Vn = R, К с = W / f - i - ^ c o s e p ,

( 1.5)

где

u>3 — угловая скорость вращения Земли;

 

R — Rо-f-Я — расстояние

от

центра

Земли

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

 

аппарата;

от

центра

Земли

до

поверхности

 

Ro — расстояние

 

геоида, которое также меняется в зависимости

 

от местонахождения.

 

 

 

скоростей

в

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

 

Ф иг. 1.5. Определение

направления

 

вертикалей

 

системе

координат

 

 

 

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

нат

относительно инерциальной системы координат:

 

 

 

V,

W F .

 

 

 

= —- =

—- + «3cos ф;

 

 

 

 

И

и

 

 

 

 

Vt

w P

(1.6)

 

 

ш7, == — tgcp = 0>з sin cp -f —- tg ? ;

 

 

3

Y i

Ъ . L

 

 

 

 

R

R

 

Зная

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

аппарата

в горизонтальной системе

кординат и угловые скорости этой си­

стемы координат относительно неподвижной (инерциальной) в еоот-

14

ветствии с законами механики, получим следующие значения уско­ рений, измеряемых акселерометрами:

= V%+ «>,, — <az Vч

an ~

+ шс Pj — “eP c+ SV

(1.7)

ac =

Pc + u)£ V,, —

+ go

 

где go gv gc — составляющие ускорения сил тяготения в предполо­ жении, что направление сил тяжести не совпадает с геоцентрической

вертикалью.

 

 

 

е учетом

(1.5) и (1.6):

Рассмотрим состав ускорений (1.7)

1/£ = W N, V n - R,

 

V K=

W N + iPu)3 cos ©-

W N ш3 sin ®;

Дa%■= co-q Vc —(Oj- Vy =

 

 

W l

 

2WE a>3 sin <p-I—

- tg <p-)-

+ ^

32/?sin 2? + - | n V

 

 

 

2

2

 

 

Дач== cot l/€ — l/t =

-

W n ^

- _ 2Н/я о)3со8<р — <o35# c o s2cp;

Дас cot

P ; = P <J>3 COS f R

Ц7д, co3 sin ш

R

W N W E *

— — - tg ф

/?

Направления поля тяготения имеют некоторые отклонения от геоцентрической вертикали, вызванные сжатием Земли и неравно­ мерным распределением масс внутри объема Земли (гравитацион­ ные аномалии). Оценим величину отклонения поля тяготения Земли от геоцентрической вертикали, вызванную сжатием Земли. Это от­ клонение происходит в основном только в плоскости географического меридиана. Разность между направлением линии отвеса, совпадаю­ щим с направлении силы тяжести Q, и геоцентрической вертикалью (фиг. 1.5) равна:

Дер = ll',5sin2(p.

(1.9)

Имеются две причины появления угла Д<р:

1)отклонение силы тяжести вследствие действия центробежных сил, вызванных вращением Земли (угол Д<Р]). и

2)отклонение сил тяжести вследствие несовпадения поля тяго­ тения с геоцентрической вертикалью (угол Д©2)

Угол A'pj может быть вычислен следующим образом:

sin Дф] ~ Acpj о>з2 R cos <?' sin ср' = 5',9 sin 2if>',

Si>

где <р и ср' — геоцентрическая и астрономическая широты.

15

Угол отклонения поля тяготения от геоцентрической вертикали

Д<р2 = Д(р — Д<р] = 5',6 sin 2ср".

( 1. 10)

Гравитационные аномалии вызывают отклонение поля тяжести Зем­ ли в пределах единиц угловых секунд. Притяжение Луны незначи­ тельно отклоняет поле тяготения Земли. Для практической навига­ ции, осуществляемой с погрешностями около Г, можно ограничить­ ся учетом угла Д? 2 и считать

gt = g sin Д?2 =s g&?t;

g ч = g cos Д»2 = g:

(1 .П )

gt = °-

Учитывая зависимости (1.8) и (1.11), получим следующие зна­ чения ускорений в географической системе координат:

где gT — укорение от сил тяжести на летательном аппарате.

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

Возьмем систему координат х„уаг,„ ориентированную относи­ тельно географической системы координат следующим образом: ось у а совпадает с вертикальной осью rlt оси ха и z„ лежат в -горизон­ тальной плоскости и повернуты на угол Н относительно осей Е и С

(фиг. 1.6).

Составляющие ускорений по осям х а,у а, z a,вызванных действием негравитационных сил:

(1.13)

16

Из фиг. 1.6 можно установить, что

VXa =

cos 0 — Ис sin в =

W cos + 0 ) — Rw3cos « sin в;

у у* - V

n = R;

 

 

(1.14)

Vza = V( sin 0 +

Кс cos 0 =

W sin (т + 0) + Rm3cos ®cos в ;

 

R

шУа

-

0 ;

 

 

v ,

 

(1.15)

 

 

__л

 

 

 

R

 

 

Составляющие ускорений от сил

тяготения, с учетом

зависимостей

(1.11). равны:

 

 

 

g*a — gi cos 9 — gz sin в =

 

“ =

c°s 0 ;

 

gya= g n ~ g i

0-16)

g*a = ^

sin В -j- gc cos 0 —

 

=

g Atp2 sin B-

 

Ф иг. 1.6. Измерение ускорений

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

Уравнения (1.13) могут быть записаны и в другой форме:

ах —

cos В —

sin В;

а.\ = а^ ё ^

(1.17)

aZa —-

а5 sin 0 -f ас cos И.

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

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

Совместим оси чувствительности

акселерометров с осями х, у,

z системы координат. Эта система координат, связанная с акселеро­ метрами, удерживается в заданном положении с помощью стабили­

затора. Отклонения и уходы стабилизатора приводят к тому,

что

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

(на­

пример, произвольной в азимуте)

системы координат на углы

Р и

8, как это было показано на фиг.

1.7. Используя значения направ­

ляющих косинусов

 

 

COS 0 COS Р

sin 0 sin 8 —

У—cos 0 sin p cos 8

Z . sin 0 cos 8 -j- 4- cos 0 sin P sin 8

sin P cos P cos 8

cos 8 sin 8

в . П. Селезнев

TO v ПУБЛИЧНАЯ

иАУчн-тех! 1ИНЁСКАЯ

 

Б И Б Л И О Т Е К А С Ч С Р

 

t S

 

f i О

— sin 0 cos p

cos 0 sin 8 + (1.18) + sin 0 sin pcos8

cos 0 cos 8 — SlJj 0

•'йУ 'пе - ТОЧНИЧ» кся

CCP

 

9;шсмпляр

СТАЛЬНОГО ЗАЛА

а также допуская малость углоз

J3 и

8,

получим следующие вели­

чины составляющих ускорений по осям

системы координат xyz:

CL% —

tlx g " h

/ П у г

=

\j/_\^ “i-

 

"H

 

ay — a v — oaz — $ax = g

+ A a + V v;

(1.19)

az =

a2n ЬаУа — Vz —8gT-f Дйг,

 

где

Aa.v -

a„

+ gx ;

 

 

 

 

 

 

 

 

Actv =

a

+

H

pa* ;

(l .20)

 

a„

laz -

 

 

•'a

 

 

a

 

Aaz = a ttza + g 2" .

Ф иг. 1.7. Измерение

ускорений в системе координат,

связанной

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

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

Vx= V xa + pv.vj

VyS=Vya' ^ W *a ~ $ V *a'

(1.21)

Абсолютные угловые скорости системы координат xyz относи­ тельно инерциальной системы координат, с учетом относительных угловых скоростей р, 5 и (-):

ш = (и* -j- 8 + Р<о v ;

 

ШУ=

*% +

(0)zo + ?)8 —

(1-22)

= “Ч, + Р— "Ч,8-

 

18

 

 

 

f

Ъ *

s

 

§ 1.4. ОСНОВНЫЕ СВОЙСТВА ОДНОМЕРНЫХ ИНЕРЦИАЛЬНЫХ СИСТЕМ

Инерциальные системы навигации по степени сложности и объ­ ему решаемых задач могут быть условно разделены на одномер­ ные, двухмерные и трехмерные.

Под этим подразумевается, что инерциальные системы могут из­ мерять ускорения, скорости и расстояния по направленным одной, двум и трем осям горизонтальной системы координат. Чем больше объем измерений, тем сложнее конструкция системы навигации. Из методических соображений рассмотрение систем навигации будем

производить с простейших одномерных систем.

 

 

 

 

 

В табл. 1.1

помещены схемы, на которых изображены возмож­

ные простейшие инерциальные системы навигации.

 

 

 

 

На схемах используются общие обозначения:

 

координат;

 

 

 

XoijoZo — инерциальная система

 

Лд-> Лу„,

Аг, и

A v

Ег(! — горизонтальная система координат;

Л- — акселерометры, ориентированные соот­

 

 

 

 

ветственно по осям инерциальной и го­

 

 

Т\,

Т2

ризонтальной систем координат;

 

 

 

телескопы,

следящие

за

небесными

S х„, S Vll,

 

 

 

светилами С! и С2;

 

 

ориентато-

Sz„ и S(_, S T:, Sc — измеренные

инерциальным

 

 

 

 

ром составляющие

пройденного

рас­

 

 

 

 

стояния по осям инерциальной и гори­

 

 

 

 

зонтальной

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

 

 

 

с.р.у. — счетно-решающее устройство;

 

 

 

 

 

П — привод, осуществляющий поворот аксе­

 

 

 

 

лерометров

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

стабилиза­

 

 

 

 

тора.

 

1.1, фиг.

1.2) счет­

В схемах с аналитической вертикалью (табл.

но-решающее устройство вычисляет углы поворота

(япер)

местной

вертикали относительно стабилизатора и пройденные расстояния

5'е,

S v Sc относительно

горизонтальной системы координат.

 

 

В табл. 1.1

(фиг.

3—8)

акселерометры поворачиваются относи­

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

(угол а„ер)• Поворот осуществляется приводом П (в табл. 1.1 показаны схе­

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

Гироинерциальный ориентатор (табл. 1.1, фиг. 5) содержит горизонтальные акселерометры и гиростабилизатор, жестко свя­ занные между собой. Счетно-решающее устройство вычисляет коор­ динаты местонахождения относительно горизонтальной системы координат и вырабатывает управляющие сигналы для поворота ги­ ростабилизатора и (удержания его в горизонтальной плоокости. Ана­

логичные управляющие сигналы используются в схемах в табл.

1.1

фиг. 6—8.

19

2*