книги из ГПНТБ / Селезнев В.П. Инерциальная навигация летательных аппаратов
.pdfВ частном случае, когда приращения передаточных чисел инте граторов не меняются во времени, составляющие погрешности счи сления пути
А Х ^ = 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 |
Р |
Х 0— Sx; |
|
|
|
|||
- Ъ ) — |
+Го |
%у |
'о — S y ; |
(1.117) |
Р |
|
Р |
|
|
gz) |
+ ^ 0 |
kЗг |
I z 0— s., |
|
Р |
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 |
1— 3 -—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 (р2— 2D.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