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

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

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

В соответствии со структурной схемой на фиг. 3.4,6, получаем следующее дифференциальное уравнение:

2 + ki М з Р : k\ k 2 g +

-= pV

ki -

+

кх кг До '

(k4 i

/е, к2 кяр) Д,3.

(3.25)

При соблюдении условия инвариантности k\k2 ---- —*= 0 и вве-

дении обозначений

 

 

 

R

 

 

 

 

к \ кг g к* — т « у-2-

<чо =

 

kyk2k3 = 2dm0х

получим уравнение, совпадающее с (3.23). В зависимости от знака можно увеличить (ki <4 0) или уменьшить (k4 > 0) период астроинерциальной системы по сравнению с периодом М. Шулера, не нарушая условия инвариантности (3.22) с точностью до погреш­

ности навигационного корректора.

Подбором коэффициента k3 может быть получен любой декре­ мент затухания d. Для быстрого затухания колебаний необходимо иметь d = 0,7 : 1,0. Внешняя информация может подаваться в астроинерциальную систему не непрерывно, а только в короткие промежутки времени: например, при начальной выставке (для лик­ видации начальных погрешностей), для ликвидации накопивших­ ся погрешностей в режиме «памяти», когда видимость небесных светил отсутствует. Для того чтобы процесс коррекции был по воз­

можности

коротким, необходимо сокращать

время

переходного

процесса.

С этой целью необходимо подбирать

х2 >

1 за счет из­

менения передаточного коэффициента k4. Но значительное умень­ шение периода инерциальной системы приводит к тому, что она начинает возмущаться помехами от внешней информации.

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

Погрешности возникают от ошибок внешней коррекции Ар и

Ай, поворотных и других ускорений Аа. Рассмотрим два наиболее характерных случая:

1) погрешности

Д(3 и Ар

и ускорение До

не зависят от вре­

мени и являются постоянными величинами;

 

2) погрешности

Д|3 и Д,3

являются случайными функциями

времени.

 

 

 

В первом случае установившаяся погрешность инерциальной

вертикали будет:

 

 

 

Руст-

- ^ 1 д

и - — ДР4-—2.

(3.26)

 

х2

а>0 х

 

Отсюда следует, что при больших значениях коэффициента измене­ ния частоты колебаний х(х > i) величины погрешностей, вызван-, ных ошибками в информации о скорости летательного аппарата

144

(*Р> и поворотными ускорениями (Да), существенно уменье шаются. Преобладающее значение приобретает погрешность, вы­ званная ошибками навигационного корректора в измерении коор­

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

(руст Дз).

О влиянии изменения частоты

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

вертикали

на погрешности коррекции

можно судить из следующего примера. Современные радионавига­ ционные системы определяют местонахождение летательного аппа­ рата с точностью до ДS x = 1 км на дальностях в несколько сот километров, а допплеровские радиолокационные системы измеряют путевую скорость с относительной ошибкой до 0,3 — 0,4%. Если коррекция астроинерциальной системы ведется одновременно: по координатам — от радионавигационной системы, а по скоро­ сти — от допплеровской радиолокационной системы, то

др = Д5,

1

= 0',54;

Др

0,0041/,

R

6,37-10G

 

 

R

Положив скорость летательного

аппарата

V =

1 км/сек, получим

Др =

0,22-10-2 угл. мин/сек.

 

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

для х

1

2d

0,004 V

 

 

 

■3'.7:

 

 

I Р у с т |

R

для х =

10

>uo

 

2d .0,004Ух = 0',54 + 0',37.

: Р у ст I = А ?

 

 

Ю»0

R

 

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

минуты.

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

погрешности Др и Др являются случайными функциями време­ ни. Оценивая эти погрешности их спектральными плотностями Д и Д5 р , можем написать выражение для среднеквадратической по­ грешности вертикали

—2

1

оо

 

 

 

 

 

! (гю) |2

(ш) du

те I'

Ф2 (lot) |2 ДSfi (и>) d и),

(3.27>

е* =

 

оо

 

 

 

 

где

 

Ф, (1ш) ]2 =

 

Ц)р(х2 —

I)2

.

(3.28)

 

 

 

 

 

 

 

 

(<оо X2 —

(О2)2 +

4 d 2 0)0 * 2 “ >2

 

 

 

 

Ф2(/«))р =

 

Ad 2 а>о х2___________

 

(3.29)

 

 

 

 

 

 

(о>0 X2--<02)2 -f- 4d2 Мр X2 0)2

10. В. Г1. Селезнев

145

Спектральные плотности ошибок Д53(ш) и \Sa («4 определяются н основном свойствами измерительных систем, дающих величиныS x и V и колебаниями объектов, на которых установлены системы Влияние колебаний объектов на величины погрешностей особенно значительно в тех случаях, когда результаты измерения зависят от положения объекта по отношению к вертикали места. В таких слу­ чаях колебания, являющиеся случайными функциями времени, определяют в значительной мере спектральные характеристики по­ грешностей Д^е и Д . Допплеровский измеритель путевой ско­ рости и является такой системой, погрешности которой определя­

ются колебаниями объекта. При из­ мерении координат местонахожде­ ния объекта при помощи радиона­ вигационных систем влияние коле­ баний менее значительно. Учитывая

сказанное, можно взять

для спект­

ральных

плотностей

Д5р(ш)

и

Д5р(<»)

следующие выражения

[8]:

Фиг. 3.5. Графики спектральных плотностей

4х, а»

(3.30)

4/., 32 (г2 +

- 2)

Д.^э (и>)

(3.31)

( г 2 ю2)2

« 2

где

,2 ^ >4 + Й2,

)-i — параметр затухания; й — резонансная частота.

На фиг. 3.5 даны графики функции ДS$ и ДS? .

Подставляя выражения (3.28) — (3.31) в (3.27), получим после интегрирования и преобразований

 

 

 

~2

“ 2 —2

5

(3.32)

 

 

 

 

*р 'Г

где еЛ.р и

среднеквадратические значения

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

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

системы, проистекающие от неточной

информации

о координатах и скорости объекта. При этом имеем:

 

72

‘"о

(*2 — 1 )2 (2й?% -л +

),,)

(3.33)

Едгр —

 

a»o*2 + 2</Q0xX1)

 

 

Хз(х? +

 

“ 2.

16rf2o2 [«ш2'/.2 Дд + Kxz 2(Дс 4- «>2 Х 2 Д0)

 

 

 

 

пг2 ш0 у.Д

 

где

 

 

 

 

 

 

у-4 ) 2 +

(2А — v2 ) (2/гг2 — v 2 io 2 -a 2) ;

д =

( 24

Да = 22[—\z*— u^x4) J- (2Л — v*)(z2-— v2)];

146

Лс - -

{z* -

*«) -J- (2 h - v-) -f (2hz2 -

®« xM; '

 

=

Z i

<■>< 7* - I - z2 (2h

V2 );

 

 

/г == V 2 d2 1— 2 d Y d 2— 1 + \' 2 d — 1 I - 2d Y d 2 - T ;

2A --

2<o2x2 (2a!2 - 1),

 

 

 

v2- 2 z 2 -

4ш*.

 

 

 

 

 

Проведем

анализ

погрешностей гх? и гх%.

Из выражения

(3.33) следует,

что при

/•=

1, т. е. в случае системы с периодом

Шулера, позиционная погрешность

еХд обращается в нуль.

С воз­

растанием х погрешность

гх?

возрастает и при большом

сокра­

щении периода инерциальной системы стремится к своему макси­

мальному значению s ^ max -- 2з|з.

 

 

Скоростная погрешность также зависит от -/..

Если предполо­

жить, что d

1 и X, < (о0, то при изменении / от

х = 1 до •/. -> оо

погрешность

 

 

 

 

 

 

« Ь - 4 .

 

 

 

°>о 2:2

 

при х = 1 начинает возрастать, достигая максимума

 

- 2

при

х ——

 

 

 

 

шо

 

При дальнейшем возрастании коэффициента погрешность на­

чинает уменьшаться, достигая

нуля при х-^со.

Такой характер

зависимости погрешности ех^

от х

объясняется тем, что с ростом

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

ную частоту 12 .

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

стемы появятся

неучитываемые погрешности от угловых колеба­

ний самолета.

 

 

 

 

Выбор

коэффициента х

можно

обусловить

требованием

«>0 у. < 12.

12 •=

0,02 11сек, то

х < 16,3. По-видимому, десятикрат­

Если

ное уменьшение периода будет приемлемым.

 

Если параметр затухания

> 12,

то функция зд.,

не имеет

максимума по х и монотонно убывает при возрастании х.

Ю*

14/

§3.5. КОМПЛЕКСНЫЕ АСТРОИНЕРЦИАЛЬНЫЕ СИСТЕМЫ НАВИГАЦИИ

ВТРЕХМЕРНОМ ПРОСТРАНСТВЕ

Для навигации летательного аппарата в трехмерном простран­ стве необходимо знать три координаты местонахождения (X, Y, Z) и составляющие скорости (X, Y, Z) относительно избранной системы координат.

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

Структурные особенности комплексной астроинерциальной системы

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

Однократное интегрирование сигналов акселерометров ах, ау, аг дает скорости полета, а двухкратное интегрирование — мгно­ венные значения координат. Полученные значения координат по­

даются в

вычислитель Bi

для

формирования

компенсационных

сигналов

g XK, g yK, g z к

и в

вычислитель В2,

дающий угловые

координаты для управления телескопами Ti, 'Гг, Т3._Сигналы рас­ согласования AS, Avj, АС между задаваемыми S, т], С и фактиче­

скими S, т), С угловыми координатами пеленгуемых небесных тел или звезд поступают в вычислитель В3 поправок АА\ АУ, AZ к ко­ ординатам местоположения. Эти поправки передаются через уси­ лительные звенья цепи обратной связи в схему измерения и инте­ грирования ускорений.

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

1) Б о д н е р В. А., С е л е з н е в В. П. К теории невозмущаемых систем с тремя каналами автокомпенсации ускорений от сил тяготения, корректируемых внешней информацией. Известия АН СССР, ОТН, «Энергетика и автоматика», 1960, № 3.

148

 

 

 

 

 

Кгх

Кгх

IX

0

В г

 

 

 

 

 

 

 

 

 

 

 

9 ^

Р

Р

 

SX

 

 

d Y

 

 

 

 

 

 

 

 

Нзу

 

 

%

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

л

 

 

5,

 

MZ

 

*4;

»Ук.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х?я

I Kzi

 

 

 

 

 

 

 

3z ZL

tE

 

 

 

 

 

 

 

 

 

 

 

Ф ii г.

3.6. Структурная схема комплексной аетроинерциальной систе­

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

Т ,, Т 2, Т3 — телескопы:

Bj, В2,

Вя вычислители; k\x , k<iK f

ф t

k^z — пере­

 

 

 

 

даточные коэффициенты

 

 

 

Вычислитель Bt должен формировать сигналы компенсации

ускорений от сил тяготения по формулам:

 

 

 

 

П

 

 

 

П

 

 

Я х к

Ц

/

' т

( S ,

• - * / ) .

V

/ - ^ < - S

' y - К.),

 

.

,

Ri

 

 

.

,

Ri

 

 

 

 

 

 

 

(■=. 1

 

(3.34)

 

 

 

 

■я, к

У /

^ A S z - Z i ) ,

 

 

 

 

 

 

 

 

 

 

 

 

i-1

R i

 

 

 

 

 

 

 

 

 

 

 

 

149

где т„ Х (, Yt, Z —

соответственно

массы- и координаты

центров

/?, —

небесных тел;

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

аппарата до

не­

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

S x, S y , S z —

бесных тел;

местоположения,

измеренные

на­

координаты

 

вигационной

системой;

 

 

 

f - - постоянная тяготения.

 

S 4, S z

и

Вычислитель В2 на основе входных координат 5,.,

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

? -

F x [X ,

У ,

Z,

Л',,

У/,

Zf);

 

rt = F 2 (X,

У,

Z,

X i% УР Z ()\

 

C=

f'8(AC

У,

Z,

X it

Kf,

Z,).

(3.35)

Примеры этих уравнений можно найти в § 2.11. Поскольку вместо

координат X, Y,

Z вычислитель В2 получает измеренные координа­

ты S x, S y, S z,

то вычисленные

им

координаты $, т(, " отлича­

ются от измеренных телескопами координат £, rt, С

На основе рассогласования

Д$ =

f — Дт) = ^’—v(, Л^==^_Л-*

в вычислителе В3 производится вычисление поправок ДА', ДК, \Z .

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

_ d F x\ Х +

dFJ Д У + dFx ДА;

 

дХ

дУ

l z

 

 

_ д Г 2 ДА'4-

ДК +

dF, ДА;

 

~ дХ

дУ

dZ

 

 

dFs

дУ ДГ-(-

dF,

ДА.

 

дХ

dZ

 

Отсюда

Dv

 

 

Г)

 

д7

D

У,

-JL

D

^

D

где

 

 

 

 

dF\

dX l

dP\

 

 

dX

 

 

дУ

dZ

 

 

dF2 dF,

dF2

 

 

dX

дУ

Tz

 

 

dFt

dF,

dF,

 

 

aX

dY

dZ

 

 

(3.36)

(3.37)

150

а определители Dx, Dy, D z получаются из D заменой соответствую­ щих столбцов правыми частями уравнений (3.36). Таким образом, уравнения (3.37) определяют программу работы вычислителя Вз.

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

Предположим, что система координат xyz, по осям которой на­ правлены оси чувствительности акселерометров, повернута отно­ сительно инерциальной системы координат x0y0Zo на углы а, {3, у (фиг; 3.7). Измеряемые акселерометрами ускорения будут:

ах =

+ Аа*г :

 

«у= (hu +

А%г :

’ (3.38)

а* -

+

дя*г •

 

где

 

 

 

ax,= X + g x\

 

=

Г + g,r

 

a2o =

’Z - ' - g z.

(3.39)

Ф и г. 3.7. Отклонение системы координат (xyz) акселе­ рометров относительно инерциальной (.vof/nzo) системы * координат

Для малых углов я, |3, у получаем

к а л-г — 8а —

ааг„;

Аа>у =

Т«-’„ — Ря»-.,;

Да*г =

вл.,-0-

( 3. 40)

151

Из структурной схемы (фиг. 3.6) получаем выражения для из­

меряемых скоростей:

 

 

 

 

 

= (ах ёх К Кх АХ) — - Кх Кх АХ + Х 0->

 

 

 

 

Р

 

 

 

Sy = (а, -

g y к -

Л3у ДК) ^

— Лг3у kiy ДК + К0;

(3.41)

 

 

 

Р

 

 

 

S z = [ a z -

g 2 к -

А3г AZ)

-

k3z kiz AZ + А

 

и координат:

 

 

Р

 

 

 

 

 

 

 

 

 

 

5 , =

^ 5 , +

^

;

 

 

 

 

D

 

 

 

 

5 у =

^ - 5 у +

Г0;

(3.42)

 

 

 

Р

 

 

 

 

S z = - ^ S z + Z Q,

 

 

 

 

р

 

 

 

где Хо, Уо, Z0 и j 0,

Y0, Z n— начальные значения координат и ско­

ростей.

 

 

 

 

 

 

Выражения (3.42) являются уравнениями астроинерциальной

системы. Запишем эти уравнения в приращениях. Полагая

 

АХ = Sx X;

ДГ = 5 У- К ;

A Z = S z - Z ;

(3.43)

Аёх== ё х ~

ёхк,

Аёу = ё у ~ ё У«> Аёг = ё г ~ ё г

(3.44)

и Кх Кх — КуКу —

получаем вместо (3.42)

 

2 +

2dx о хр +

22) ДА"=

Д g x +

АаХг;

 

2 +

2dy ilyр +

22) ДК =

Agy +

ДаУг ;

(3.45)

(p2 + 2 dz Qzp + U2z) b Z = b g z + baZT,.

 

где

 

 

 

 

 

i l l ^ k s x ,

il2y

 

ill

- К г

 

2 dx 2 ^ = Кх kix,

2 dy 2 у = /?3у kiy,

2 d ZQ г = k%z kiz.

 

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

'от сил тяготения

Погрешности Agz, Д^у, Д^г, формируемые на основе сигналов S x, S v, S z, отличны от нуля и равны:

A g x ^ — ^ A X + AAY+BAZ-,

Аёу — ААХ — Т)2 ДГ + CAZ;

(3.46)

Д g z = ЯАХ + C A Y - y]2 AZ,

152

где

■-Е/5 1 - 3

'X - X t 2

 

/=■1

R l

Ri

 

 

lTli 1 - 3

Y - Y, 's

 

/«1

R l

Ri

 

П

 

z —zt 2

 

Е Л 1 - 3

 

 

я /

Д,

П« , ( X - X t) ( Y - Y t) .

л= з £ /

ы\

Ri

R l

 

 

n

m, (АГ-Л-,) (£ -■ ? ,).

B =

Л

 

- » E

R l

r-1

Ri

 

 

П

 

{Y - Уi)(Z ~ Z , )

c-«S/

я?

i-i

 

 

 

Подставляя (3.46) в (3.45), получим

{p2 + 2d, Й,/> + <4 + -^) ДА' -

ЛАК - /^ Z = Aa,r ;

(3.47)

(3.48)

— ЛДА" Ч- (/>•’ -!- 2dу 9.ур

 

| (->2т)р ДГ — CAZ =

Д%г ;

(3.49)

- 5ДЛ" - СДГ + (р 2 f 2

^

‘2 , р 4- <Z2 -г,2) AZ -

Aa*r .

 

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

Погрешности в измерении ускорений Ай*,., АаУг , ДaZ{ , вызван­

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

Рассмотрим устойчивость системы автокомпенсации ускоре­ ний от сил тяготения в предположении, что отклонения телескопов отсутствуют (До*г = ДаУг = а2г = 0), а коэффициенты системы

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

F A P )

-

В

 

А (Р) =■ — А

-

С - F x (p)Fy(P)Fz (p)~

~ В

F Z(P)

 

~ C ' F X{P)

&Fy(p) - A2 Fz { р ) ~ 2 А В С = 0 ,

(3.50)

153