Способ определения истинного курса с помощью двухканального гироскопического датчика угловой скорости

 

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

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

Известен способ определения истинного курса с помощью гироскопического датчика угловой скорости (см., например, книгу Б.И. Назарова и Г.А. Хлебникова "Гиростабилизаторы ракет" М., 1975, стр. 193-196), согласно которому курсовое направление измерительной оси гироскопа определяют аналитически путем вычислений по результатам измерений горизонтальной составляющей угловой скорости вращения Земли. С целью повышения точности определения курса применяют методы, направленные на уменьшение влияния дрейфа гироскопа. Так гироскопический датчик угловой скорости, измерительная ось которого расположена в плоскости горизонта, поворачивают на разные азимутальные углы и снимают на этих углах выходную информацию с гироскопа в виде электрического напряжения с эталонного сопротивления, стоящего в цепи обратной связи. На основе того что ряд моментов гироскопа, не связанных с его корпусом, меняет знак при разворотах его в азимуте, разработаны известные методы, которые позволяют выполнить взаимную компенсацию вредных моментов. Однако эти известные методы автокомпенсации не позволяют полностью исключить погрешность определения курса от вредных моментов гироскопа, а также от других источников, например невыставки осей чувствительности гироскопа в горизонт.

За прототип взят способ определения истинного курса с помощью гироскопического датчика угловой скорости, основанный на его аналитическом вычислении (см. статью S.P.Redkin "Analytical Gyrocompassing by Angular Rate Sensor Based on Dynamically Tuned Gyro", I Saint Petersburg International Conference on Gyroscopic Technology, May, 1994).

В этом способе проводят предварительную настройку по выставке измерительных осей гироскопа в горизонт и привязки их к осям, связанным с объектом, включают гироскоп в режим работы с обратной связью по току датчика момента, определяют сигнал с эталонного сопротивления датчика угловой скорости в первом положении "a" на угле измеряемого курса, поворачивают гироскоп в азимуте во второе положение "c" вокруг собственной оси на 180 градусов, а затем в третье положение "b" на 90 градусов и измеряют сигналы с гироскопа на этих углах, после чего искомый азимутальный угол оси Y первого канала находят аналитически с использованием значений сигналов с эталонного сопротивления, полученных на этих углах с помощью следующего соотношения (1) где U_11a, U_11c, U_11b - напряжения с эталонного сопротивления датчика угловой скорости соответственно в "a", "c", "b".

При наличии практически всегда имеющихся малых углов наклона плоскости расположения осей чувствительности гироскопа в известном способе выражение для определения искомого азимутального угла можно представить в следующем виде: (2) где δ_b,δ_c - погрешность азимутального позиционирования гироскопа соответственно в положении "b" и "c", K_н11 - крутизна гироскопа по напряжению по первому каналу, H₁ - кинетический момент гироскопа,
Ω - угловая скорость вращения Земли,
ϕ - широта местоположения объекта,
ΔM^1x_oc,ΔM^1x_ob - изменение коэффициента момента, независящего от ускорения по оси X гироскопа при его азимутальном повороте из положения "a" соответственно в положение "c" и "b",
M^1x _1yb, M^1a _1ya - коэффициент, характеризующий момент от осевого разбаланса в положении "b" и "a",
ΔM^1x_1yc - изменение моментного коэффициента от осевого разбаланса при повороте гироскопа из положения "a" в положение "c",
M^1x _1xa, M^1x _1xb - квадратурные моментные коэффициенты в положениях "a" и "b",
ΔM^1x_1xc - изменение квадратурного коэффициента при повороте гироскопа из положения "a" в положение "c",
ΔM^1x_1zc,ΔM^1x_1zb - изменение моментного коэффициента, пропорционального ускорению вдоль оси собственного вращения при повороте гироскопа из положения "a" соответственно в положение "c" и "b",
ΔM^1x_1z1yc,ΔM^1x_1z1yb - изменение моментного коэффициента, пропорционального второй степени ускорения вдоль оси собственного вращения z при повороте гироскопа из положения "a" соответственно в положение "c" и "b",
M^1x _1y1ya, M^1x _1z1za, M^1x _1y1yb, M^1x _1z1zb - коэффициенты, обуславливающие момент от неравножесткости соответственно в положении гироскопа "a" и "b",
ΔM^1x_1y1yc,ΔM^1x_1z1zc - изменение коэффициентов, обуславливающих изменение момента от неравножесткости при повороте гироскопа из положения "a" в положение "c",
M^1x _1x1ya, M^1x _1x1yb - коэффициентные моменты от перекрестного влияния ускорений со ответственно в положении "a" и "b",
ΔM^1x_1x1yc - изменение коэффициента момента от перекрестного влияния ускорений при повороте гироскопа из положения "a" в "c".

Сравнивая выражение (1), используемое в прототипе для определения азимутального угла, с полученным уточненным выражением (2), можно отметить, что в известном способе искомый угол определяется с погрешностью, обусловленной:
- ошибкой азимутального позиционирования δ_b,δ_c гироскопа в положениях "b" и "c",
- нескомпенсированными остаточными моментами ΔM^1x_a,ΔM^1x_1a,
- проекциями на измерительную ось вертикальной составляющей угловой скорости вращения Земли при наклоне плоскости измерительных осей гироскопа на углы γ и ν.

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

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

Целью настоящего изобретения является обеспечение возможности определения истинного курса с помощью гироскопического датчика угловой скорости в условиях неподвижного наклонного основания.

Поставленная цель достигается тем, что в известном способе определения истинного курса с помощью датчика угловой скорости, включающем предварительную привязку измерительных осей гироскопа к осям, связанным с объектом, применение гироскопа в режиме обратной связи по току датчика момента, определение сигналов с эталонных сопротивлений датчика угловой скорости при операции предварительной настройки проводят фазировку системы измерения напряжений с эталонных сопротивлений гироскопа, обеспечивающую в горизонтальном положении измерительных осей отрицательные приращения сигналов с первого и второго каналов гироскопа при азимутальном повороте его корпуса против часовой стрелки на 90 градусов из положения, при котором измерительная ось у первого канала направлена на север, а измерительная ось x второго канала на восток, определяют коэффициенты модели дрейфа гироскопа, а при гирокомпасировании в искомом курсовом положении объекта определяют углы наклона плоскости измерительных осей гироскопа по тангажу и крену, угол широты местоположения объекта, напряжения с эталонных сопротивлений гироскопа по первому и второму каналам, а затем определяют значение истинного курса измерительной оси первого канала гироскопа по следующей формуле:
K = 2π - arctg a, если b₁ > 0, b₂ > 0;
K = π - arctg a, если b₁ > 0, b₂ < 0;
K = π - arctg a, если b₁ < 0, b₂ < 0;
K = - arctg a, если b₁ < 0, b₂ > 0;
где a = b₁b₂ ^-1;

b₂= U₁₁K_н11Cosγ+ω²_gyCogγ-ΩSinϕSinνCosγ,
ω²¹_gx,ω²¹_gy - модель дрейфа гироскопа по его осям X и Y,


ω^1x_o,ω^1y_o - коэффициенты модели дрейфа, независящие от ускорения,
ω^1x_1x,ω^1x_1y,ω^1y_1x,ω^1y_1y,ω^1y_1z,ω^1x_1z - коэффициенты модели дрейфа, пропорциональные ускорению в первой степени,
- коэффициенты модели дрейфа, пропорциональные ускорению во второй степени,
U₁₁, U₁₂ - напряжения с эталонных сопротивлений датчика угловой скорости соответственно по первому и второму каналам,
K_н11, K_н12 - крутизна гироскопа по напряжению для первого и второго канала,
ν,γ - углы наклона плоскости осей чувствительности гироскопа соответственно по тангажу и крену,
Ω - угловая скорость вращения Земли,
ϕ - широта местоположения объекта.

Схема двухканального гироскопа в режиме обратной связи по моменту показана на фиг. 1. Первый канал гироскопа содержит датчик угла 1, измеряющий поворот корпуса гироскопа 9 относительно измерительной оси Y_г1, усилитель 2, эталонное сопротивление 3, датчик момента 4, создающий момент вокруг оси X_г1. Второй канал содержит датчик угла 5, измеряющий поворот корпуса гироскопа относительно измерительной оси X_г1, усилитель 6, эталонное сопротивление 7, датчик момента 8, создающий момент вокруг оси Y_г1.

В установившемся режиме для изодромной обратной связи напряжения с эталонных сопротивлений датчика угловой скорости можно представить в следующем виде
U₁₁= -ω²¹_gyK_н12+ω²¹_yK_н12;
U₁₂= ω²¹_gxK_н11-ω²¹_xK_н11, (3)
где ω²¹_x,ω²¹_y - проекции абсолютной угловой скорости корпуса гироскопа соответственно на оси X_г1 и Y_г1.

Свяжем с гироскопом систему координат O₂₁X₂₁Y₂₁Z₂₁. Оси O₂₁X₂₁ и O₂₁Y₂₁ направлены по измерительным осям гироскопа, а ось O₂₁Z₂₁ вдоль оси собственного вращения гироскопа. С объектом свяжем систему координат O_cX_cY_cZ_c. Ось O_cX_c направлена вдоль продольной оси объекта, ось O_cY_c перпендикулярна оси O_cX_c и располагается в плоскости симметрии объекта, ось O_cZ_c - перпендикулярна этой плоскости.

Система координат O_дX_дY_дZ_д является системой координат Дарбу. Взаимное расположение систем координат показано на фиг. 2.

В соответствии с взаимным расположением выбранных систем координат проекции абсолютной угловой скорости на оси гироскопа в случае неподвижного относительно Земли объекта можно представить в следующем виде:
ω²¹_x = CosψSinγSinνΩCosϕ+SinψCosγΩCosϕ;
ω²¹_y = CosψCosνΩCosϕ. (4)
Подставив выражение (4) в выражение (3), можно получить следующую систему уравнений:
C₁₁Cosψ+ = d₁;
C₂₁Cosψ+C₂₂Sinψ = d₂, (5)
где
C₁₁= -ΩCosϕCosν;
C₂₁= -ΩCosϕSinγSinν;
C₂₂= -ΩCosϕCosγ;
d₁= -U₁₁K_н11-ω²¹_gy+ΩSinϕSinν;
d₂= U₁₂K_н12-ω²¹_gx-ΩSinϕSinγCosν.
Решая уравнение (5), определяем выражение для cosψ и sinψ. В практике определения курсового угла с помощью датчика угловой скорости его вычисления проводят через функцию арктангенс. Тогда аналитическое выражение для определения азимутального угла ψ можно представить в следующем виде:
ψ = arctg a, (6)
где a = b₁b₂ ^-1;

b₂= U₁₁K_н11Cosγ+ω²¹_gyCosγ-ΩSinϕSinνCosγ.
При использовании функции тангенса квадрант, в котором расположен азимутальный угол ψ, может быть вычислен по соотношению знаков b₁ и b₂. Зная квадрант, формулу для вычисления угла истинного курса можно представить в виде
K = 2π - arctg a, если b₁ > 0, b₂ > 0;
K = π - arctg a, если b₁ > 0, b₂ < 0;
K = π - arctg a, если b₁ < 0, b₂ < 0;
K = - arctga, если b₁ < 0, b₂ > 0. (7)
При изодромной обратной связи в установившемся режиме стационарную модель дрейфа гироскопа можно представить в виде


где f_x, f_y, f_z - составляющие удельной силы вдоль положительных направлений осей, связанных с гироскопом.

Для неподвижного объекта уравнение (8) может быть выражено через составляющие ускорения силы тяжести путем подстановки
f_x = g_x, f_y = g_y, f_z = g_z. (9)
С учетом взаимного расположения осей координат Дарбу и системы координат, связанной с гироскопом, показанных на фиг. 2, проекции ускорения силы тяжести на оси, связанные с гироскопом, можно представить в виде
g_x= gSinγCosν, g_y= -gSinν, g_z= -gCosγCosν; (10)
где g - ускорение силы тяжести.

С учетом (9) и (10) модель дрейфа гироскопа будет иметь вид


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

На фиг. 1 показана схема двухканального гироскопа в режиме обратной связи по моменту.

На фиг. 2 показано взаимное расположение систем координат.

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

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

Экспериментальные исследования проводились с помощью датчика угловой скорости на базе динамически настраиваемого гироскопа ГВК-6. Предварительно определялись коэффициенты модели дрейфа гироскопа. Затем динамически настраиваемый гироскоп ГВК-6 устанавливался на платформу стенда таким образом, что плоскость измерительных осей гироскопа имела наклон по тангажу ν = 10 град., по крену γ = 10 град. Стенд также обеспечивал задание различных курсовых углов наклонного гироскопа с точностью 10 угл. секунд.

Динамически настраиваемый гироскоп работал в режиме датчика угловой скорости путем включения его в систему обратных связей по моментам. С эталонных сопротивлений, стоящих в цепях обратных связей, измерялись напряжения. Гироскоп из исходного положения, при котором курсовой угол равен 0, поворачивался против часовой стрелки на курсовые углы 90, 135, 180, 225, 270, 315, 360 градусов. На каждом из этих углов производилось измерение напряжений по двум каналам гироскопа, на базе которых с помощью программ, реализующих известный и предлагаемые способы, определялись значения курсовых углов на устанавливаемых азимутальных положениях гироскопа.

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

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

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

Из зависимости, представленной на фиг. 3, видно, что применение предлагаемого способа позволяет выполнять гирокомпасирование с применением гироскопического датчика угловой скорости в условиях наклонного объекта. В то время как применение известного способа в этих условиях не позволяет выполнить гирокомпасирование из-за неприемлемо больших погрешностей. В предлагаемом способе погрешность определения курса в диапазоне курсовых углов от 0 до 360 градусов не превысила 1,284 град.

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

Формула изобретения

Способ определения истинного курса с помощью двухканального гироскопического датчика угловой скорости, включающий предварительную привязку измерительных осей гироскопа к осям, связанным с объектом, применение гироскопа в режиме обратной связи по току датчика момента, определение сигналов с эталонных сопротивлений датчика угловой скорости, отличающийся тем, что предварительно проводят фазировку системы измерения напряжений с эталонных сопротивлений гироскопа, обеспечивающую в горизонтальном положении измерительных осей отрицательные приращения сигналов с первого и второго каналов гироскопа при азимутальном повороте его корпуса против часовой стрелки на 90^o из положения, при котором измерительная ось Y первого канала направлена на север, а измерительная ось X второго канала на восток, определяют коэффициенты модели дрейфа гироскопа, а при гирокомпасировании в искомом курсовом положении объекта определяют углы наклона плоскости измерительных осей гироскопа по тангажу и крену, угол широты местоположения объекта, напряжения с эталонных сопротивлений гироскопа по первому и второму каналам, а затем определяют значение истинного курса измерительной оси первого канала гироскопа по следующей формуле:
k = 2π - arctg a, если b₁ > 0, b₂ > 0;
k = π - arctg a, если b₁ > 0, b₂ < 0;
k = π - arctg a, если b₁ < 0, b₂ < 0;
k = - arctg a, если b₁ < 0, b₂ > 0;
где
a = b₁b₂ ^-1;

b₂= U₁₁K_H11Cosγ+ω²¹_дyCosγ-ΩSinϕSinυCosy;
ω²¹_дx, ω²¹_дy - модель дрейфа гироскопа по его осям Х, Y;


ω^1x_o, ω^1y_o - коэффициенты модели дрейфа независящие от ускорения;
ω^1x_1x, ω^1x_1y, ω^1x_1z, ω^1y_1x, ω^1y_1y, ω^1y_1z, ω^1x_1z - коэффициенты модели дрейфа пропорциональные ускорению в первой степени;
- коэффициенты модели дрейфа, пропорциональные ускорению во второй степени;
U₁₁, U₁₂ - напряжения с эталонных сопротивлений датчика угловой скорости соответственно по первому и второму каналам;
K_H11, K_H12 - крутизна гироскопа по напряжению для первого и второго каналов;
υ, γ - углы наклона плоскости осей чувствительности гироскопа соответственно по тангажу и крену;
Ω - угловая скорость вращения Земли;
ϕ - широта местоположения объекта.

РИСУНКИ

Рисунок 1, Рисунок 2, Рисунок 3