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

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

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

В частном случае, когда приращения передаточных чисел инте­ граторов не меняются во времени, составляющие погрешности счи­ сления пути

А Х ^ = k \x ( Д ^ 2 л - f e . v 4 " k'ix A k 3x) S x - j - A k 3.ir X Q

 

(Ak2y къу +

Ak3y)Sy 4- Akzy У01\

(1.105)

AZb = k\z (Д^2г k'iz“f- ^2z

“H ^^32^0 ^

 

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

Из рассмотренного следует, что погрешности счисления пути,

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

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

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

где i — порядковый номер погрешности.

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

§1.13. МЕТОДЫ УМЕНЬШЕНИЯ ПОГРЕШНОСТЕЙ ОРИЕНТАТОРОВ БЕЗ АВТОКОМПЕНСАЦИИ УСКОРЕНИЙ ОТ СИЛ ТЯГОТЕНИЯ

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

а) компенсацией погрешностей в начальных условиях (перед вылетом);

б) компенсацией погрешностей с помощью программных уст­ ройств;

в) внешней коррекцией погрешностей с помощью навигацион­ ной информации от других устройств.

Рассмотрим содержание этих методов уменьшения погреш­ ностей.

64

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

Суть метода заключается в том, что перед вылетом осущест­ вляют «обнуление» акселерометров и контроль остаточных по­ грешностей после компенсации. Для этого к выходным сигналам ах , ау , (^акселерометров добавляют компенсационные сигналы ахк, аук, а2К, которые полностью уничтожают выходные сигналы. Усло­ вие компенсации сводится к следующему:

ах -

ахк =

Аах ;

ау аук*~Аау ; аг - агк =

Дяг,

(1.107)

где Да*,

Дау ,

Даг

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

 

 

Поскольку летательный аппарат неподвижен, то

X =

У — Z —

= 0; X =

Хо; Y •= Vo; Z — Zo.

 

 

Из условия

компенсации (1.107) и уравнений (1.88) получаем

следующие значения компенсационных сигналов:

 

 

CLvk ~~

UyK

®zk — -

gx„

? y .

g z 0

+

? 0

а Уо

 

a 0

<!

 

 

+

J 0

^

 

 

( 1 . 1 0 8 )

 

 

 

- ! -

« о

a Xo

-

\ a >\,

A a z ,

 

 

-

 

где

р0, 30 — начальные

углы отклонения акселерометров.

 

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

следует, что компенсационные сигналы

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

Устранение накапливаемой с течением времени погрешности ориентатора на выходах вторых интеграторов позволяет частично скомпенсировать влияние остаточной погрешности «обнуления» акселерометров и инструментальные погрешности интеграторов для начальных условий полета.

Остаточные погрешности компенсации

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

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

 

 

/

t

 

 

 

Д5,к =

f [{Agx + A a x)(dty~,

 

 

 

6

О

 

 

 

 

t

t

 

 

 

-

0f 0J (д£у +

А«у) (^O2;

(1.109)

 

 

t

t

 

 

 

=

f f (Дgz +

даг) (dt)\

 

 

 

o 6

 

 

где

Agx ~ g x gxц;

 

Дg v = gy - g v„; Дgz = g, —- gz.

 

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

 

 

 

65-

— изменения ускорений сил тяготения относительно началь­

ных значений.

 

и среднего

действий

Для

летательных аппаратов ближнего

наиболее

часто применяется такая начальная

ориентация

осей

 

.

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

акселе­

 

 

рометров, при

которой

 

 

ось у

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

 

 

ной

вертикали,

ось

х

 

 

в направлении полета, а

 

 

ось 2

— перпендикулярна

 

 

этому

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

Для

 

 

такого расположения осей

 

 

xyz

в начальный момент

 

 

времени будем иметь

 

 

 

g x0=

 

= 0;

g Va= g t .

 

 

 

При полете в верти­

 

 

кальной

плоскости

ху

 

 

(фиг. 1.32) ускорения сил

 

 

тяготения равны:

 

 

 

 

 

 

gx= -g sin апср;

 

 

 

 

 

g у — geos fltnep,

 

 

 

где

aneD — Угол поворота

 

 

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

 

 

сительно

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

 

 

системы

координат

 

Фи г. 1.32. Схема'поворота

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

 

 

t

 

 

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

вертикали

а лер

Vx dt.

 

 

 

 

 

При дальностях полета, измеряемых десятками и сотнями километ­

ров, можно считать sin апер =

апер и

cosa„ep= 1. В этом случае'

остаточная погрешность компенсации

в счислении

пути будет

 

t

t

 

 

AS,K_ J (VnepW;

 

 

0

6

 

 

 

/

t

 

(1.110)

bSyK=

\ l ( g - g 0){dtf.

 

0.0

 

 

Для высот полета H « Ro, где R0 — радиус Земли, в первом при­

ближении получим

t i

ASVk =

j ^Hi d t f .

(1.111)

0 о и

>ёб

Так, например, при полете с постоянной скоростью Vx — 1000 м/сек на высоте Я •= 20 км в течение 3-х минут накапливаются погреш­ ности = 1520 м и A*SVK = 1020 м.

Компенсация погрешностей с помощью программных устройств

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

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

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

Компенсация погрешностей выполняется в два этапа. На пер­ вом этапе производится «обнуление» акселерометров и компенса­ ция погрешностей в начальных условиях. Второй этап компенса­ ции продолжается в течение всего полета и заключается в следую­ щем. С начала взлета генератор t стабильной частоты (фиг. 1.33) запускает программное счетно-решающее устройство, которое вы­

рабатывает компенсационный

сигнал g x ,

суммирующийся

с вы­

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

показан только

один

канал ориентатора).

сигнал £>мр

уничтожает

дополнительное

Компенсационный

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

gX{i от сил тяготения,

возникающее при изме­

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

Условием

компенсации яв­

ляется

р>

 

= Дач — Дах

,

 

 

 

 

(1.112)

 

олпр

 

л

дПр>

где ах

1

сигнал начальных

ускорений;

 

 

У ( * р - * о ) .

 

 

(1.113)

 

 

 

 

 

 

 

(Япр-До)3

5*

67

X np, RnP — вычисленные

программным

счетно-решающим

 

 

устройством координата X и расстояние R' от ле­

 

 

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

 

 

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

 

 

 

 

ДаЛ-

и ДаЛ.

— погрешности компенсации ускорений в начальных

 

 

условиях и программного компенсирования.

Координаты Х пр и Rnp определяются с погрешностями:

 

 

Х пр =

X + АЛ’пр,

Rnp = R + Д/?Пр;

 

где X и R — истинные значения,

ДЛ'пр и \R np — погрешности.

Погрешность программной

компенсации

ускорений сил тяго­

тения:

 

 

 

 

 

 

 

 

 

 

 

д

Н-' (ХпР-- х 0)

 

(Д' -- А'р)

^ д

/

Д^чр_____ ЗД/^пр \

пр

( / ? п р

р 0)'

~

( / ?

/ ? о ) 8

U А 'о

R - R o l

 

 

 

 

 

 

 

 

 

 

 

(1.114)

Погрешность

счисления

пути,

вызванная

только

погрешностью

программной компенсации, равна:

 

 

 

 

 

 

■t

 

'

Д^пр

 

 

 

 

 

AS.

 

 

З Д

^ п р

 

(1.115)

 

[ ф —

 

 

 

 

 

R -

-Яо

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычисление погрешностей для других каналов ориентатора произ­ водится аналогично.

Программное компенсирование ускорений сил тяготения по­ зволяет существенно уменьшить погрешности счисления пути. Так,

например,

при

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

ДЛ*

10%

погрешность счисления пути уменьшается в 10 раз

------ — в

X — А'о

 

 

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

§ 1.14. ГИРООРИЕНТАТОРЫ С АВТОКОМПЕНСАЦИЕИ УСКОРЕНИЙ ОТ СИЛ ТЯГОТЕНИЯ

Принцип автокомпенсации (или самокбмпенсации) основан на использовании результатов измерения самого ориентатора для фор­ мирования поправок к его показаниям [27].

Полная схема автокомпенсации погрешностей от ускорений сил тяготения представлена на фиг. 1.34. В этой схеме формирова­ ние сигналов коррекции в цепях обратных связей производится

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

снеограниченным изменением высоты полета.

68

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

-

и !1' S*

gy =

, у' Sy

(1.116)

 

 

*4у ^=f-,

И'

где

R = V s l + S 2y + S l

Фиг. 1.34. Структурная схема гироориентатора с аатокомпенсацией ускорений от сил тяготения;

Л |Л-, * iv . z ~ передаточные коэффициенты акселерометров;

fyy* *3у* ^2г* ^Ъг ~~

передаточные коэффициенты интеграторов

 

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

р

J

Р

Х 0Sx;

 

 

 

- Ъ ) —

+Го

S y ;

(1.117)

Р

 

Р

 

 

gz)

+ ^ 0

kЗг

I z 0s.,

 

Р

J

Р

 

 

где A 'o, Yo, Z0, X 0, K 0 , Z 0 — начальные значения координат и ско­ ростей летательного аппарата.

69

Преобразуем уравнения (1.117), для чего выберем передаточ­ ные числа акселерометров и интеграторов из условий:

к

1д- /й2.1- к?,х

к ] v k l y к з у

k \ z

z ^ 3 Z ^ ■

 

Находим

k i x

/ с I

к \ у = к-4 у\

k \ z k i z .

 

 

р 2А Х = Дg x + Да.,-,. ;

 

 

 

 

 

 

 

 

Р 2 A Y - - A g y +

 

д а Уг ;

 

(1.118)

где

 

p*AZ=Agz + Дагр ,

 

 

 

А';

Д К - 5 У

Г; AZ = Sz

(1.119)

Д *= ,S.V

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

 

 

Agx =

g-д- — g x;

Дgy = g.v -

g > ;

Дg z= g2 gz

( 1. 120)

— погрешности компенсации ускорений от сил тяготения. Если осуществлена точная компенсация ускорений от сил тя­

готения

(Agx --- Agy — Дgz = 0)

и отсутствуют уходы стабили­

затора

(ДаЛ-г =

АаУг = AaZr =-0), то

погрешности системы будут

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

только неточностью задания начальных условий.

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

Формируемые вычислителем сигналы компенсации gx , gy , gz могут быть не равны составляющим ускорения g x , g y , g 2- Возни­ кающие погрешности подаются в инерциальную систему и, йройдя замкнутый цикл, опять поступают в вычислитель. Если инерциаль­ ная система устойчива, то погрешности будут затухать. В против­ ном случае они будут, возрастать, что вызовет неточное определе­ ние местонахождения.

Если подстановкой (1.93), (1.94) и (1.116) в (1.120) образовать разности Д^г, Д^у, &gz и разложить их в ряды по приращениям ДА', ДГ, ДZ, то, ограничиваясь линейными членами, получим

Д*, = — * 1- 3 W

AX 3 — ДГ 3 X Z AZ

J

 

R

 

R 2

 

R 2

 

R2

 

Agy =

_g_

 

К2 \

 

3 —

 

YZ

AZ

1

13 -—I д К

ДЛ" - 3 —

( . )

 

R

R2

R 2

 

1 121

 

 

R2)

 

 

 

Д * ,= ~ - £ 1 - 3

Z 2

 

 

 

Z Y

a y

Л

R2

AZ 3 * 1 A X - 3 —

 

*

R

 

 

R2

 

R2

 

 

Подставляя (1.121)

в (1.118), найдем .

 

 

 

 

 

ip2+ r£) ДА' — АДК — BAZ --

А а х

 

 

-

A A X X { p ,2+ yl2y) A Y ~ C ^ Z = A a y v

 

(1 .1 2 2 )

-

ВАХ - СА К+

(р2 Г;2) AZ =

Да,х ,

 

 

70

где

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

Рассмотрим поведение системы (1.122) при условии, что уходы гироскопов отсутствуют. Если система (1.122) устойчивая, то по­ грешности автокомпенсации будут стремиться к нулю. В случае неустойчивости погрешности неограниченно нарастают. Рассмот­ рим частный случай, полагая, что коэффициенты системы (1.122) меняются медленно по сравнению со временем протекания процес­ сов в системе актокомпенсации. Принимая для простоты коэффи­ циенты (1.123) постоянными, можем судить об устойчивости на основе свойств характеристического определителя.

 

Р1-!- \v

А

В

 

D

А

 

С

(1.124)

 

В

С

Р2 + ъ;

 

Если раскрыть этот определитель и приравнять нулю, то полу­ чим характеристическое уравнение системы (1.122):

Рй+ (■»£ + ri), +

п]) Р* + К Tiy а" 42v + ч1 -nl — А2 — В' -

С2)р 2+

:

А*г£-В *г£/ - С 2т?х - 2 А В С = 0.

(1.125)

Поскольку в этом уравнении отсутствуют члены с нечетным» степенями оператора р, то из этого можно сделать заключение о не­ устойчивости системы автокомпенсации.

Преобразуем уравнение (1.125), учитывая (1.116),

 

 

 

 

+

=

'

2

t

ril +

Tj*

r\2z — А 2 — В 2

— С2 — — 3

(1.126):

vl X t -

А* ■?,- В2 Г!2 - С2

- 2

АВС - -

2

 

 

 

ё 2= s i +

g] + g\-

 

 

Подставляя (1.126)

в

(1.125)

и производя

преобразования,,

получим

(р2+

Q2)2 (р22D.2) = 0,

 

(1.127)

 

 

где Q2 —

R

71

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

бания с периодом Т\

, колебания с тем же перио­

дом и нарастающей пропорционально времени амплитудой и, на­ конец, апериодического экспоненциально нарастающего дЕНжения к. постоянной времени

1

(1.128)

7'2 = УТТг

2 V 2‘ т: '

Анализ системы уравнений (1.122) в общем виде затрудняется наличием переменных коэффициентов. Однако рассмотрение част­ ного случая полета по прямой представляется вполне возможным. Допустим, что летательный аппарат движется в направлении оси t/о. В этом случае Н = Y, X — Z =• 0, gx = g%='J, gy -^g м уравнения ( 1..L22) принимают вид:

(р2+ (- ^ Д* = Да,г ;

 

( р * - 2& ) А У = АаУт;

(1.129)

(Р2 --2) \ Z = Д«гг •

 

Отсюда следует, что в направлениях, перпендикулярных на­

правлению поля сил тяготения, погрешности

ДА" и AZ меняются

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

-условий. Так,

например, для начальных условий при t = О, ДА' =

— ДА'0, AZ — ДZ0 и YX —ДА'П, AZ =

AZ0 будем иметь

 

 

ДА' = ДА'0 cos Ш +

sin 9J-

 

 

 

 

 

 

(1.130)

 

 

AZ — AZ0 cos Ш -|-

дг„ ■sin й/.

В направлении поля тяготения погрешность ДY меняется по

закону

 

 

 

Y2~ 9(

 

 

д к = !

АГ> + - £

а -

ДК0-

е V2S', (1.131)

 

2 \

 

/ 2

<>

 

у

2 О

т. е. погрешность нарастает по экспоненциальному закону.

На фиг.

1.35 представлены графики изменения относительно

величины

ДК

в зависимости от времени и высоты полета.

----

ДГа

*Если учесть правые части уравнений (1.122), то характеристическое урав­

нение становится (p2 +Q2) (р2—2Q2) = 0 и нарастающее пропорционально времени движение отсутствует.

72

равен периоду обращения спутника, летящего по кру­

Кривые 1,2 и 3 соответствуют высотам полета 0, 1000 и 2000 км. К увеличением высоты полета погрешность накапливается медлен­

нее, так как влияние сил тя­

 

 

 

 

готения уменьшается.

 

 

 

 

 

При

движении

лета­

 

 

 

 

тельного

аппарата

по

про­

 

 

 

 

извольной

траектории

си­

 

 

 

 

стема автокомпенсации

ус­

 

 

 

 

корений

от

сил

тяготения,

 

 

 

 

как следует

из

уравнений

 

 

 

 

(1.122), неустойчива по всем

 

 

 

 

каналам.

 

неустойчивость

 

 

 

 

Однако

 

 

 

 

не может явиться критерием

 

 

 

 

неприменимости

инерциаль­

 

 

 

 

ных систем. Дело в том, что

 

 

 

 

при нулевых начальных ус­

 

 

 

 

ловиях погрешности

будут

 

 

 

 

отсутствовать, а при началь­

Фиг.

1.35. Графики изменения

началь­

ных погрешностях их нара­

стание происходит

медлен­

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

но. Кроме того, как

это

по­

в зависимости от времени и высоты по­

 

лета

 

 

казано ниже, систему мож­

подачи в систему

дополнительной

но сделать

устойчивой

путем

внешней информации. .

 

 

 

 

 

 

Аналогия

изменения погрешностей системы

 

 

 

‘ автокомпенсации движению спутника

 

 

Отметим аналогию между характером изменения погрешно­

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

системы

навигации и движением

спутника.

Период

изменения

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

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

системы 7\ =

]/" •

говой орбите радиуса R вокруг Земли с первой космической ско­ ростью

 

V, — RQ

2т,R

(1.132)

 

- V g R .

 

 

Л

 

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

тяготения меняется

по экспоненциальному закону с

постоянной

1 .

' ~R

 

 

времени 7 \= ^v = ^|/

— и при ненулевых начальных условиях не­

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

Н2 = У 2 t i R = 2gV R .

1.133)

73