Способ алгоритмической компенсации погрешности гирокомпасирования с применением гироскопического датчика угловой скорости

 

Изобретение относится к области точного приборостроения и может быть использовано при создании гирокомпасов и курсоуказывающих устройств аналитического типа. Предварительно привязывают измерительные оси гироскопа к осям, связанным с объектом. Используют гироскоп в режиме обратной связи по току датчика момента. Определяют сигналы с эталонных сопротивлений датчика угловой скорости. Определяют коэффициенты модели дрейфа гироскопа. Находят погрешности гирокомпасирования в шести пространственных положениях гироскопа, когда курсовой угол равен 0 и 270 ^o, плоскость осей чувствительности горизонтальна, наклонена только по тангажу, наклонена только по крену. В рабочем режиме находят широту местоположения объекта, углы наклона плоскости измерительных осей гироскопа по тангажу и крену. Проводят алгоритмическую компенсацию дрейфа гироскопа. Определяют алгоритмическую компенсацию погрешности в виде аналитического выражения, использующего значения нескорректированного угла истинного курса, погрешностей гирокомпасирования в шести положениях и углов наклона при их предварительном определении, широты местоположения, углов наклона объекта при гирокомпасировании. Обеспечивается повышение точности гирокомпасирования. 2 ил.

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

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

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

За прототип взят способ определения истинного курса с помощью гироскопического датчика угловой скорости, основанный на его аналитическом вычислении (см. статью V.L. Budkin, S.P. Redkin "Gyrocompassing on a Movable Object by Means of Angular Rate Sensor on the Basis of Dynamically Tuned Gyro". 5th Saint Petersbury International Conference on Integrated Navigational Systems, May 25-27,1998).

В этом способе проводят гирокомпасирование с применением датчика угловой скорости на основе динамически настраиваемого гироскопа (ДУС-ДНГ) при бесплатформенной установке его на объекте. Предварительно привязывают измерительные оси гироскопа к осям, связанным с объектом, согласовывают знаки изменения напряжений с эталонных сопротивлений гироскопа с направлением азимутального поворота его корпуса вокруг оси собственного вращения, определяют коэффициенты модели дрейфа гироскопа, широту местоположения объекта, углы наклона плоскости измерительных осей гироскопа по тангажу и крену, измеряют напряжения с эталонных сопротивлений гироскопа, работающего в режиме датчика угловой скорости, а затем азимутальный угол объекта определяют с помощью выражения: где U₁₁, U₁₂ - напряжения с эталонных сопротивлений гироскопа соответственно по первому и второму каналам, К_н11, К_н12 - крутизна гироскопа по напряжению по первому и второму каналам, - угловая скорость вращения Земли, - широта местоположения объекта, , - углы наклона плоскости осей чувствительности гироскопа соответственно по тангажу и крену, ^г1_ду, ^г1_дх - дрейф гироскопа соответственно по первому и второму каналам.

В известном способе производится аналитическая компенсация дрейфа гироскопа. Для этого дрейф гироскопа может быть представлен в виде математической модели. Стационарную модель дрейфа в условиях неподвижного объекта, учитывающую дрейф, не зависящий от ускорения, дрейф от разбаланса вдоль оси собственного вращения и квадратурный дрейф, можно представить в виде: ^г1_дx = ^1x₀+^1x_1xg_x+^1x_1уg_у, ^г1_ду = ^1у₀+^1у_1xg_x+^1у_1уg_у, (2)
где ^1x₀, ^1у₀ - коэфффициенты модели дрейфа, характеризующие дрейф, не зависящий от ускорения,
^1x_1x, ^1у_1у - коэффициенты модели дрейфа, характеризующие дрейф от разбаланса вдоль оси собственного вращения,
^1x_1у, ^1у_1x - коэффициенты модели дрейфа, характеризующие квадратурный дрейф,
g_x, g_у - проекции ускорения силы тяжести g на оси, связанные с гироскопом,
g_x = gsincos, g_у = -gsin.
Для точной аналитической компенсации дрейфа необходимо точное знание коэффициентов модели дрейфа, которые определяются предварительно в результате калибровочных испытаний гироскопа. Однако при калибровочных испытаниях эти коэффициенты определяются с погрешностями, они также изменяются в процессе наработки гироскопа и выполнить их повторные калибровочные испытания без снятия с гироскопического устройства не представляется возможным. Ошибки в задании коэффициентов модели дрейфа обуславливают погрешности гирокомпасирования.

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

Технический результат достигается тем, что в известном способе гирокомпасирования с применением гироскопического датчика угловой скорости, включающем предварительную привязку измерительных осей гироскопа к осям, связанным с объектом, использование гироскопа в режиме обратной связи по току датчика момента, определение сигналов с эталонных сопротивлений датчика угловой скорости, согласование знаков изменения сигналов с эталонных сопротивлений гироскопа с направлением азимутального поворота его корпуса вокруг оси собственного вращения, определение коэффициентов модели дрейфа гироскопа, а в рабочем режиме нахождение широты местоположения объекта, углов наклона плоскости измерительных осей гироскопа по тангажу и крену, алгоритмическую компенсацию дрейфа гироскопа, аналитическое определение нескорректированного угла К истинного курса объекта, дополнительно предварительно устанавливают гироскоп в положение, при котором угол истинного курса равен 0 ^o, плоскость измерительных осей гироскопа горизонтальна, измерительные оси Х_г1,У_г1 гироскопа направлены соответственно на восток и север, выполняют аналитическое гирокомпасирование и находят погрешность гирокомпасирования ₀₁, потом устанавливают гироскоп в положение, при котором угол истинного курса равен 270 ^o, горизонтально расположенные измерительные оси Х_г1 и У_г1 направлены соответственно на север и запад, проводят гирокомпасирование и находят погрешность гирокомпасирования ₀₂, затем гироскоп устанавливают в положение, в котором угол истинного курса равен 0^o, измерительная ось Х_г1 находится в горизонте и направлена на восток, а измерительная ось У_г1 наклонена по тангажу на угол ₁, проводят гирокомпасирование и находят погрешность гирокомпасирования ₁, после этого гироскоп устанавливают в положение, в котором измерительная ось Х_г1 находится в горизонте и направлена на север, а измерительная ось У_г1 наклонена по тангажу на угол ₁, проводят гирокомпасирование и находят погрешность гирокомпасирования ₂, затем гироскоп устанавливают в положение, в котором угол истинного курса равен 0^o, измерительная ось У_г1 находится в горизонте и направлена на север, а измерительная ось Х_г1 наклонена по крену на угол ₁, выполняют гирокомпасирование и находят погрешность гирокомпасирования ₁, после этого гироскоп устанавливают в положение, в котором угол истинного курса равен 270^o, измерительная ось У_г1 находится в горизонте и направлена на запад, а измерительная ось X_г1 наклонена по крену на угол ₁, проводят гирокомпасирование и находят погрешность гирокомпасирования ₂, а затем в рабочем режиме определяют алгоритмическую компенсацию погрешности гирокомпасирования в виде следующего аналитического выражения:

где
^г1_дх = ^1x₀+^1x_1xgsincos-^1x_1ygsin,
^г1_дy = ^1y₀+^1y_1xgsincos-^1y_1ygsin,
^1x₀ = cossin₀₁,
^1y₀ = -cossin₀₂,
^1x_1y = (^1x₀-cossin₁)(gsin₁)^-1,
^1y_1y = (^1y₀+coscos₁sin₂)(gsin₁)^-1,
^1x_1x = -(^1x₀-coscos₁sin₁)(gsin₁)^-1,
^1y_1x = -(^1y₀+cossin_y2)(gsin₁)^-1,
- угловая скорость вращения Земли,
- широта местоположения объекта,
, - углы наклона плоскости осей чувствительности гироскопа соответственно по тангажу и крену,
g - ускорение силы тяжести.

Найдем аналитическое выражение для погрешности гирокомпасирования с применением датчика угловой скорости на основе динамически настраиваемого гироскопа. Для известного способа аналитическое гирокомпасирование выполняется с использованием выражения (1).

Представим выражение (1) в виде:
tg() = AB^-1, (3)
где

B = (U₁₁K_н11cos+^г1_дуcos-sinsincos).
Положим, что A и B есть отклонения от точных значений А и В. Эти отклонения обусловят погрешность гирокомпасирования и выражение (3) будет иметь вид:
tg(+) = (A+A)(B+B)^-1. (4)
В соответствии с тригонометрической формулой сложения
tg(+) = (tg+tg)(1-tgtg)^-1. (5)
Из выражения (4) с учетом (5) погрешность гирокомпасирования можно представить в следующем виде:

где

U₁₁, U₁₂ - погрешность измерения напряжения с эталонных сопротивлений гироскопа соответственно по первому и второму каналам,
K_Н11, K_Н12 - погрешность задания крутизны гироскопа по напряжению по первому и второму каналам,
- погрешность задания широты местоположения объекта,
, - погрешность задания углов наклона плоскости осей чувствительности гироскопа соответственно по тангажу и крену,
^г1_ду, ^г1_дх - остаточный дрейф гироскопа после его алгоритмической компенсации соответственно по первому и второму каналам,
^г1_дх = ^1x₀+^1x_1xgsincos-^1x_1уgsin,
^г1_ду = ^1у₀+^1у_1xgsincos-^1у_1уgsin.
Для изодромной обратной связи гироскопа измеряемые с эталонных сопротивлений напряжения можно представить в следующем виде:
U₁₁ = a₁₁cos+a₁₂sin+a₁₃,
U₁₂ = a₂₁cos+a₂₂sin+a₂₃, (7)
где a₁₁ = coscosK^-1_н11,
a₁₂ = 0,
a₁₃ = (-^г1_ду+sinsin)K^-1_н11,
a₂₁ = -cossinsinK^-1_н12,
a₂₂ = -coscosK^-1_н12,
a₂₃ = (^г1_дх+sinsincos)K^-1_н12.
Подставляя (7) в (6) и полагая
после преобразований получим:

Полученное выражение (8) характеризует погрешность гирокомпасирования с применением гироскопического датчика угла. Эта погрешность обусловлена остаточным дрейфом ^г1_дх, ^г1_ду после его алгоритмической компенсации, углами наклона по крену и тангажу , широтой местоположения и азимутальным углом . При использовании угла истинного курса К выражение (8) имеет вид:

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

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

Для синуса погрешности гирокомпасирования можно записать следующее выражение:

Используя выражение (10), определим коэффициенты модели остаточного дрейфа через погрешность гирокомпасирования . Для этого введем следующие системы координат: Х_д У_д Z_д, Х_с Z_с У_с, Х_г1 У_г1 Z_г1. Система координат Х_д У_д Z_д является географической. Ось Х_д направлена на север, ось У_д - вертикально вверх, а ось Z_д - на восток. В системе координат Х_д У_д Z_д проекциями угловой скорости Земли являются горизонтальная составляющая cos, направленная на Север, и вертикальная составляющая sin, где - угол широты места. С объектом связана система координат X_c У_с Z_c. Ось X_c направлена вдоль продольной оси объекта, ось У_с перпендикулярна оси У_с и располагается в плоскости симметрии объекта, ось Z_c - перпендикулярна этой плоскости. Система координат X_г1У_г1 Z_г1 связана с гироскопом. Ось X_г1 направлена по оси
чувствительности второго канала гироскопа, ось У_г1 - по оси чувствительности первого канала, а ось Z_г1 - по оси собственного вращения. Вдоль оси Z_г1 направлен вектор кинетического момента Н гироскопа. Пусть привязка осей гироскопа выполнена таким образом, что они совпадают с осями объекта.

Установим гироскоп в положение 1, как показано на фиг.1. В этом положении угол истинного курса равен 0^o, плоскость измерительных осей гироскопа горизонтальна, измерительные оси X_г1, У_г1 гироскопа направлены соответственно на восток и север. Выполним аналитическое гирокомпасирование известным методом и найдем погрешность гирокомпасирования ₀₁. Подставив ₀₁ у (10), в этом положении будем иметь
^1x₀ = cossin₀₁. (11)
Установим гироскоп в положение 2, как показано на фиг.1. В этом положении угол истинного курса равен 270^o, горизонтально расположенные измерительные оси X_г1,У_г1 направлены соответственно на север и запад. Проведем гирокомпасирование известным методом и определим погрешность гирокомпасирования ₀₂. Подставив ₀₂ в (10), в этом положении получим
^1у₀ = -cossin₀₂. (12)
Установим гироскоп в положение 3, как показано на фиг.1. В этом положении угол истинного курса равен 0^o, измерительная ось X_г1 находится в горизонте и направлена на восток, а измерительная ось У_г1 наклонена по тангажу на угол ₁. Выполним гирокомпасирование известным методом и определим погрешность гирокомпасирования ₁. Подставив ₁ в (10), в этом положении получим:
^1x_1у = (^1x₀-cossin₁)(gsin₁)^-1. (13)
Установим гироскоп в положение 4, как показано на рис. 1. В этом положении угол истинного курса равен 270^o, измерительная ось X_г1 находится в горизонте и направлена на север, а измерительная ось У_г1 наклонена по тангажу на угол ₁. Проведем гирокомпасирование известным методом и найдем погрешность гирокомпасирования ₂. Подставив ₂ в (10), в этом положении будем иметь:
^1у_1у = (^1у₀-coscos₁sin₂)(gsin₁)^-1. (14)
Установим гироскоп в положение 5, как показано на фиг.1. В этом положении угол истинного курса равен 0^o, измерительная ось У_г1 находится в горизонте и направлена на север, а измерительная ось X_г1 наклонена по крену на угол ₁. Выполним гирокомпасирование известным методом и найдем погрешность гирокомпасирования ₁. Подставив ₁ в (10), в этом положении будем иметь:
^1x_1x = -(^1x₀-coscos₁sin₁)(gsin₁)^-1. (15)
Установим гироскоп в положение 6, как показано на фиг. 1. В этом положении угол истинного курса равен 270^o, измерительная ось Y_г1 находится в горизонте и направлена на запад, а измерительная ось X_г1 наклонена по крену на угол ₁. Проведем гирокомпасирование известным методом и найдем погрешность гирокомпасирования ₂. Подставив ₂ в (10), в этом положении будем иметь:
^1у_1х = -(^1у₀+cossin₂)(gsin₁)^-1. (16)
Таким образом, предлагаемый способ алгоритмической компенсации погрешности гирокомпасирования с применением гироскопического датчика угловой скорости имеет следующие отличия от известного способа:
- в предварительной операции находят коэффициенты модели остаточного дрейфа, обусловленного несовершенством модели дрейфа, используемого в известном способе;
- в операции гирокомпасирования вводится новое действие связанное с алгоритмической компенсацией непосредственно погрешности гирокомпасирования с использованием новой аналитической зависимости;
- коэффициенты модели остаточного дрейфа находятся новым способом на основе полученных связей между коэффициентами остаточного дрейфа и соответствующих им погрешностям гирокомпасирования.

На фиг. 1 показаны системы координат: географическая, связанная с объектом гироскопом.

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

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

Экспериментальные исследования проводились с помощью датчика угловой скорости на основе динамически настраиваемого гироскопа ГВК-6 392116. Предварительно в калибровочных испытаниях определялись коэффициенты модели дрейфа гироскопа. Для испытываемого гироскопа коэффициенты модели дрейфа имели следующие значения:
^1x₀ = 0,532 град/час,
^1у₀ = -0,659 град/час,
^1x_1у = 0,827 град/часg,
^1у_1у = 0,094 град/часg,
^1x_1x = 0,275 град/часg,
^1у_1x = 0,833 град/часg.
После этих калибровочных испытаний и наработки 500 ч гироскоп был установлен на платформу наклонно-поворотного стенда. Для проведения испытаний по точности гирокомпасирования с использованием ранее полученных коэффициентов и модели дрейфа плоскость измерительных осей гироскопа наклонялась на угол тангажа =10^o и угол крена =10^o. Известным способом выполнялось аналитическое гирокомпасирование на разных углах истинного курса и находилась погрешность гирокомпасирования. Погрешность гирокомпасирования, полученная при использовании известного способа, представлена зависимостью 1 на фиг.2.

После этого с помощью предлагаемого способа определялись коэффициенты остаточного дрейфа. При испытаниях делались наклоны по тангажу и крену на углы ₁= 10^o, ₁=10^o. В результате были получены следующие коэффициенты остаточного дрейфа:
^1x_0/ = -0,055 град/час,
^1у₀ = -0,167 град/час,
^1x_1у = -0,276 град/часg,
^1у_1у = -0,205 град/часg,
^1x_1x = 0,135 град/часg,
^1у_1x = -1,215 град/часg.
На разных углах истинного курса при наклоне гироскопа по крену и тангажу на углы = =10^o выполнялось гирокомпасирование с применением предлагаемой алгоритмической компенсации погрешности гирокомпасирования и скорректированное значение истинного курса определялось с помощью выражения
K_к = K-, (17)
где К - значение истинного курса, определенное известным методом;
- предлагаемая алгоритмическая компенсация погрешности гирокомпасирования, определяемая выражением (9).

Найденная погрешность гирокомпасирования с применением предлагаемой алгоритмической компенсации представлена зависимостью 2 на фиг.2.

На фиг.2 показаны экспериментальные зависимости погрешностей гирокомпасирования на разных азимутальных углах при использовании известного и предлагаемого способов. Из зависимостей, представленных на фиг.2, видно, что применение предлагаемого способа позволяет повысить точность гирокомпасирования. Так, погрешность гирокомпасирования с применением известного способа в диапазоне курсовых углов от 0 до 360^o не превысила 2,6^o, а при использовании предлагаемого способа погрешность в этих условиях не превысила 0,3 ^o. Применение предлагаемого способа позволило уменьшить погрешность гирокомпасирования в 8,7 раза.

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

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

Способ алгоритмической компенсации погрешности гирокомпасирования с применением гироскопического датчика угловой скорости, включающий предварительную привязку измерительных осей гироскопа к осям, связанным с объектом, использование гироскопа в режиме обратной связи по току датчика момента, определение сигналов с эталонных сопротивлений датчика угловой скорости, согласование знаков изменения сигналов с эталонных сопротивлений гироскопа с направлением азимутального поворота его корпуса вокруг оси собственного вращения, определение коэффициентов модели дрейфа гироскопа, а в рабочем режиме нахождение широты местоположения объекта, углов наклона плоскости измерительных осей гироскопа по тангажу и крену, алгоритмическую компенсацию дрейфа гироскопа, аналитическое определение нескорректированного угла К истинного курса объекта, отличающийся тем, что предварительно устанавливают гироскоп в положение, при котором угол истинного курса равен 0^o, плоскость измерительных осей гироскопа горизонтальна, измерительные оси Х_г1, У_г1 гироскопа направлены соответственно на восток и север, выполняют аналитическое гирокомпасирование и находят погрешность гирокомпасирования ₀₁, потом устанавливают гироскоп в положение, при котором угол истинного курса равен 270^o, горизонтально расположенные измерительные оси Х_г1 и У_г1 направлены соответственно на север и запад, проводят гирокомпасирование и находят погрешность гирокомпасирования ₀₂, затем гироскоп устанавливают в положение, в котором угол истинного курса равен 0^o, измерительная ось Х_г1 находится в горизонте и направлена на восток, а измерительная ось У_г1 наклонена по тангажу на угол ₁, проводят гирокомпасирование и находят погрешность гирокомпасирования ₁, после этого гироскоп устанавливают в положение, в котором угол истинного курса равен 270^o, измерительная ось Х_г1 находится в горизонте и направлена на север, а измерительная ось У_г1 наклонена по тангажу на угол ₁, проводят гирокомпасирование и находят погрешность гирокомпасирования ₂, затем гироскоп устанавливают в положение, в котором угол истинного курса равен 0^o, измерительная ось У_г1 находится в горизонте и направлена на север, а измерительная ось Х_г1 наклонена по крену на угол ₁, выполняют гирокомпасирование и находят погрешность гирокомпасирования ₁, после этого гироскоп устанавливают в положение, в котором угол истинного курса равен 270^o, измерительная ось У_г1 находится в горизонте и направлена на запад, а измерительная ось Х_г1 второго канала наклонена по крену на угол ₁, проводят гирокомпасирование и находят погрешность гирокомпасирования ₂, а затем в рабочем режиме определяют алгоритмическую компенсацию погрешности гирокомпасирования в виде следующего аналитического выражения

где ^г1_дx = ^1x₀+^1x_1xgsincos-^1x_1ygsin;
^г1_дy = ^1y₀+^1y_1xgsincos-^1y_1ygsin;
^1x₀ = cossin₀₁;
^1y₀ = -cossin₀₂;
^1x_1y = (^1x₀-cossin₁)(gsin₁)^-1;
^1y_1y = (^1y₀+coscos₁sin₂)(gsin₁)^-1;
^1x_1x = -(^1x₀-coscos₁sin₁)(gsin₁)^-1;
^1y_1x = -(^1y₀+cossin₂)(gsin₁)^-1;
- угловая скорость вращения Земли;
- широта местоположения объекта;
, - углы наклона плоскости осей чувствительности гироскопа соответственно по тангажу и крену;
g - ускорение силы тяжести.

РИСУНКИ

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