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

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

 

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

Известен способ измерения [RU 2366961 C1, G01P 15/13, 10.09.2009], заключающийся в том, что ускорение преобразуют в отклонение подвижной системы акселерометра с последующим преобразованием отклонения в электрический сигнал, его усиливают и преобразуют в последовательность временных интервалов, сформированных с помощью источника стабилизированного тока и элементов нелинейного звена в виде широтно-модулированного сигнала, который преобразуют в момент импульсной отрицательной обратной связи при помощи датчика момента, а сигналы импульсной отрицательной обратной связи преобразуют в сигнал, модулированный счетными импульсами, поступающими от генератора импульсов стабилизированной частоты f_c, с последующим измерением числа счетных импульсов n_1i и n_2i, которые преобразуют в выходные сигналы акселерометра в виде разности Δn_i=n_1i-n_2i, суммы n_i=n_1i+n_2i=f_cT_i и отношения числа импульсов N_i=Δn_i/n_i, причем величину и знак входного воздействия a_вх=a_i+Δa_i, соответствующего измеряемым ускорению a_i, и погрешностям Δa_i, определяют по окончании каждого i-того периода T_i автоколебаний с помощью выходного сигнала и по формуле N_i=K_i(a_i+Δa_i), (i=1^{--- ∞}), где K_i=ml/K_дмI₀ - коэффициент преобразования входного воздействия, ml - маятниковый момент подвижной системы, K_дм и I₀ - коэффициент передачи датчика момента и поступающий от стабилизированного источника ток, величина которого ограничивает диапазон измеряемого ускорения, при этом коэффициент преобразования K_i, погрешность, например смещение нуля Δa_i=Δa₀, их стабильность, определяют при помощи формулы N_ik=K_aik(a_i+Δa₀), градуировочную характеристику и ее нестабильность (нелинейность) при помощи другой формулы N_ij=К_aij(a_ij+Δa_0j) определяют заранее и каждую раздельно по результатам калибровки в условиях гравитационного поля Земли при вертикальном и горизонтальном положениях измерительной оси акселерометра, причем измерения каждой технической характеристики выполняют на интервале времени сотен периодов автоколебаний, чтобы применить статистическую обработку и фильтрацию результатов измерений, которые используют при измерении ускорения, причем на основе заданного значения линейности градуировочной характеристики, полученного и паспортизованного в результате калибровки акселерометра на заводе-изготовителе в условиях его прецизионного оборудования и заданного диапазона измеряемого ускорения, проводят калибровку акселерометра на том же диапазоне в условиях гравитационного поля Земли при вертикальном положении измерительной оси акселерометра и на интервале времени в несколько сотен периодов автоколебаний путем измерения ускорения известной величины, при этом коэффициент преобразования, погрешность смещения нуля и их стабильность, соответствующие заданному диапазону измеряемого ускорения определяют одновременно при коротком по времени однократном измерении двух известных величин входного ускорения, соответствующих границам заданного диапазона a_вх1=+a_max+Δa₀ и a_вх2=+a_max+Δa₀ с учетом смещения нуля, полученные значения выходных сигналов N_max1=K_a1(+a_max+Δa₀), -N_max2=K_a2(-a_max+Δa₀) суммируют алгебраически и делят на удвоенное значение ускорения границы заданного диапазона (N_max1-(-N_max2))/2a_max=K_a, а затем суммируют арифметически и делят пополам (N_max1+(-N_max2)/2=K_aΔa₀, определяют значения: коэффициента преобразования K_a погрешности смещения нуля на выходе ΔN₀=K_aΔa₀ и на входе Δa₀=ΔN₀/K_a, измеренную погрешность смещения нуля компенсируют алгоритмическим образом до уровня остаточной на выходе ΔΔN₀, соответствующей остаточной погрешности смещения нуля ΔΔa₀=ΔΔN₀/K_a на входе, повторяют измерения ускорения известной величины на границах заданного диапазона вышеуказанным способом с учетом скомпенсированной погрешности смещения нуля и вычисления с использованием полученных значений выходных сигналов $$\left(N_{\max 11}-\left(N_{\max 22}\right)\right)/2a_{\max }=K_a^{*}$$ , $$\left(N_{\max 11}+\left(-N_{\max 22}\right)\right)/2=K_a^{*}\Delta \Delta a_0^{*}$$ , определяют уточненные стабильное значение коэффициента преобразования $$K_a^{*}$$ , значения погрешностей смещения нуля на выходе $$\Delta \Delta N_0^{*}$$ и на входе $$\Delta \Delta a_0^{*}$$ , стабильность полученных погрешностей вычисляют при помощи отношения измеренных и вычисленных алгебраической и арифметической сумм $$\left(N_{\max 11}+\left(-N_{\max 22}\right)\right)/\left(N_{\max 11}-\left(-N_{\max 22}\right)\right)=\Delta \Delta a_0^{*}/a_{\max }$$ , где $$\Delta \Delta a_0^{*}/a_{\max }$$ - относительная величина скомпенсированной погрешности на заданном диапазоне, а представленная в процентах $$\Delta \Delta a_0^{*}\cdot 100\%/a_{\max }={\gamma }_{\Delta \Delta a0}$$ характеризует стабильность измеренной погрешности, стабильность градуировочной характеристики одновременно, причем при стабильном коэффициенте преобразования $$K_a^{*}=const$$ , полученные значения запоминают и используют при формировании полетного задания и для измерения ускорения в заданном диапазоне с помощью выходного сигнала устройства и по формуле

$$N_a=K_a^{*}\left(a_i+\Delta \Delta a_i^{*}\right)$$ .

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

Известен также способ [RU 2140059 C1, G01C 21/18, 20.10.1999], включающий получение данных о местоположении и скорости, определение поправок к выходным параметрам инерциальной гироскопической системы выработку корректирующих сигналов и ввод их в инерциальную гироскопическую систему, при этом непрерывно сравнивают измеряемый и измеренный ранее профили рельсовой колеи в продольном направлении, в вертикальной плоскости и по взаимному положению рельсовых нитей в поперечной плоскости, а поправки к выходным параметрам инерциальной навигационной системы определяют при обработке полученных при сравнении результатов с одновременным выделением деформации рельсовой колеи, при этом корректирующие сигналы вводят в инерциальную гироскопическую систему непрерывно.

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

Наиболее близким по технической сущности к предложенному является способ [RU 2381511 C1, G01P 21/00, 10.02.2010], включающий определение характеристик унифицированных систем позиционирования на основе микромеханических акселерометров и гироскопов в условиях одновременного воздействия на него динамических в виде линейной перегрузки, или синусоидальной вибрации, или одиночных ударных импульсов, а также температурных испытательных параметров, основанный на том, что воздействие всех динамических параметров задают с помощью одного поворотного стенда, при этом воздействие линейной перегрузки обеспечивают созданием нормального ускорения, величина которого определяется соотношением: a_п=ω²·R, где a_п - нормальное ускорение, пропорциональное угловой скорости вращения поворотной платформы стенда, ω - угловая скорость вращения поворотной платформы, R - расстояние от оси вращения поворотной платформы до места установки испытуемого прибора, воздействие синусоидальной вибрации обеспечивают созданием тангенциального ускорения, величина которого определяется соотношением: a_r=∈·R, где a_r - тангенциальное ускорение, пропорциональное угловому ускорению поворотной платформы, ∈_∈(t)=E·sin(2pv·t) - задаваемое угловое ускорение поворотной платформы, R - расстояние от оси вращения поворотной платформы до места установки испытуемого прибора, Е - амплитуда углового ускорения, v - частота углового ускорения, t - время, воздействие кратковременного одиночного ударного импульса обеспечивают созданием углового ускорения поворотной платформы стенда в виде половины периода синусоиды.

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

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

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

Поставленная задача решается, а требуемый технический результат достигается тем, что в способе определения погрешностей инерциальных измерительных приборов (преимущественно гироскопов и акселерометров) при испытаниях на ударные и вибрационные воздействия, основанном на оценке погрешностей измерений по исходным данным, полученным по результатам измерений показаний инерциальных измерительных приборов, согласно предложенному изобретению, измерение показаний инерциальных измерительных приборов производят в три этапа с одинаковой продолжительностью по времени при различной для каждого из этапов пространственной ориентации инерциальных измерительных приборов, причем на первом этапе производят измерение показаний инерциальных измерительных приборов для определения начальной ориентации и интегрирования уравнения навигации, на втором этапе производят измерение показаний инерциальных измерительных приборов для интегрирования уравнений навигации, когда инерциальные измерительные приборы подвергают ударным и/или вибрационным воздействиям, а на третьем этапе производят измерение показаний инерциальных измерительных приборов для начальной ориентации и интегрирования уравнения навигации после ударных и/или вибрационных воздействий, при этом для расчета смещений нулей инерциальных измерительных приборов до ударных воздействий используют величины l_x, l_y, l_z, рассчитанные по уходам координат и углов ориентации относительно начальной матрицы $${\tilde{D}}_0$$ на участке T₁ для расчета смещений нулей инерциальных измерительных приборов после ударных воздействий используют величины l_x, l_y, l_z, рассчитанные по уходам координат и углов ориентации относительно начальной матрицы $${\tilde{D}}_3$$ на участке T₃, а для расчета смещений нулей инерциальных измерительных приборов во время ударных воздействий используют величины l_x, l_y, l_z, рассчитанные по уходам координат и углов ориентации относительно начальной матрицы $${\tilde{D}}_0$$ , где l_x, l_y, l_z, рассчитывают из соотношений l_x=B_0xt³+B_1xt², l_z=B_0zt³+B_1zt², l_y=B_1yt², коэффициенты B_0x,z, B_1x,z, B_1y определяют на каждом шаге интегрирования рекуррентной процедурой метода наименьших квадратов, а смещения нулей инерциальных измерительных приборов определяют по соотношениям:

; , ; Δsα_z=2·B_1z, Δsα_y=2·B_1y'.

На чертеже представлены:

на фиг.1 - циклограмма проведения измерений показаний приборов;

на фиг.2 - зависимость от времени параметра ориентации α;

на фиг.3 - зависимость от времени параметра ориентации β;

на фиг.4 - зависимость от времени параметра ориентации А;

на фиг.5 - временная зависимость ухода линейной координаты 1_x на интервале T₁;

на фиг.6 - временная зависимость ухода линейной координаты 1_y на интервале T₁;

на фиг.7 - временная зависимость ухода линейной координаты 1_z на интервале T₁;

на фиг.8 - временная зависимость ухода линейной координаты 1_x на интервале T₁-Т₃;

на фиг.9 - временная зависимость ухода линейной координаты 1_y на интервале T₁-Т₃;

на фиг.10 - временная зависимость ухода линейной координаты 1_z на интервале T₁-Т₃.

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

Способ позволяет обеспечить компенсацию колебательно-вращательных движений ударных и вибрационных стендов таким образом, чтобы во время механических воздействий оценивались только реальные погрешности инерциальных измерительных блоков (приборов) (ИИБ), в частности смещения нулей лазерных гироскопов (ЛГ) и маятниковых акселерометров (МА), вследствие ударных и вибрационных воздействий, в то время как кажущиеся погрешности из-за паразитных колебаний и вращений практически исключаются.

Измерения проводятся в раздельном масштабе времени: вначале проводится запись показаний ИИБ (далее исходные данные - ИД) по заданной циклограмме, а затем обработка данных.

Функционирующий ИИБ устанавливается и закрепляется на испытательном стенде.

Измерения проводятся по циклограмме, приведенной на фиг.1. Интервалы Т₁, Т₂ и T₃ выбираются равными. Их оптимальное значение 5 мин. Меньшее значение интервала подвержено искажениям из-за случайного характера ошибок измерений, а увеличение интервала приводит к неоправданному увеличению используемых ресурсов. В интервале Т₁=5 мин проводится запись ИД для начальной выставки (определения начальной ориентации) и интегрирования уравнений навигации, в интервале Т₂=5 мин проводится запись ИД для интегрирования уравнений навигации, когда ИИБ подвергается ударным или вибрационным воздействиям, в интервале Т₃=5 мин проводится запись ИД для повторной выставки ориентации и интегрирования уравнений навигации после ударных воздействий. Запись ИД проводится при трех различных положениях ИИБ, при этом оси связанной с корпусом системы координат (СВСК) ориентируются так, чтобы ось СВСК была ориентирована вертикально вверх, а оси X_p и Z_p ориентированы произвольно в плоскости, близкой к горизонтальной.

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

Для проведения расчетов введятся три системы координат.

ИИБ отождествляют с ортогональной декартовой системой координат 0_p, X_p, Y_p, Z_p, жестко связанной с ИИБ (СВСК).

Стенд отождествляют с ортогональной декартовой системой координат 0_ст, X_ст, Y_ст, Z_ст, жестко связанной со стендом - связанная топографическая система координат (СТСК).

Угловые и линейные смещения систем координат СВСК и СТСК определяют относительно геодезической системы координат в данной точке земного сфероида (ГСК), ось X_г которой направлена на Север, ось Y_г - вверх по нормали к земному сфероиду, ось Z_г - на Восток.

Текущая ориентация СВСК, первоначально совмещенной с ГСК, получается путем трех последовательных поворотов на углы, именуемые параметрами ориентации:

- поворот на угол A_z (азимут) вокруг оси Y_p СВСК, совпадающей с осью Y_г ГСК.

- поворот на угол α (крен) вокруг нового положения оси X_p СВСК.

- поворот на угол β (тангаж) вокруг нового положения оси Z_p СВСК.

При таком определении СВСК смена позиций а, б, в означает, что по оси ударных воздействий, ориентированной вдоль оси Y_г, выставляются поочередно все ЛГ, а алгоритмы и программы начальной выставки и интегрирования уравнений навигации, полученные для позиции 6, остаются неизменными.

Для неподвижной системы показания МА складываются из проекций вектора ускорения силы тяжести Земли $$\overline{G}$$ и собственных смещений нулей датчиков.

Проекции вектора $$\overline{G}$$ , измеренные первоначально на осях чувствительности ЛГ и приведенные к ортогональным осям СВСК, описываются соотношением

где α_{xp, yp, zp} - показания акселерометров, считываемые с осей СВСК.

D₀ - истинная матрица ориентации СВСК относительно ГСК.

$$\overline{G}={\left[\mathrm{0,}G\mathrm{,0}\right]}^T$$ - вектор ускорения силы тяжести Земли, заданный проекциями на оси ГСК.

[cα_ij] - матрица ориентации осей чувствительности МА относительно СВСК (далее матрица конструктивной привязки акселерометров).

[Δca_ij] - матрица погрешности измерения матрицы [cα_ij].

$$\overline{\Delta sa}={[\Delta s{a}_x,\Delta s{a}_y,\Delta s{a}_z]}^T$$ - вектор смещения нулей акселерометров.

Интегрирование уравнений навигации, в пределах данной методики, осуществляется на осях ГСК.

Проекции показаний акселерометров на оси ГСК определяются из (1) соотношениями

где D - текущее значение матрицы ориентации СВСК относительно ГСК, вычисленное по показаниям ЛГ.

В алгоритме интегрирования уравнения ориентации, используемом в данной методике, проекции вектора вращения Земли вычитаются из показаний ЛГ, так что матрица D имеет вид

где - кососимметричная матрица, составленная из смещений нулей Δsq_x,y,z

t - текущее время, отсчитываемое от начала интегрирования.

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

Погрешности Δβ и Δα определяются соотношениями

;

Из (2) с учетом (3), (4), (5) следует

Расписав (7) в скалярной форме в линейном приближении, получаем

Из (8) с учетом соотношений (6) получаем

Проинтегрировав уравнения (9) дважды, получаем соотношения для

уходов координат l_x и l_z, обусловленных инструментальными погрешностями ИИБ:

По ИД, полученным в интервале T₁, вычисляются усредненные за время T₁ показания ЛГ и матрица начальной ориентации $${\tilde{D}}_0$$ СВСК относительно ГСК.

Усредненные значения определяются соотношениями:

где $${\tilde{P}}_{x,y,z}$$ - усредненные за время Tl показания ЛГ;

P_x,y,z,i - показания датчиков ЛГ на i-м шаге съема информации;

$${\tilde{a}}_{x,y,z}$$ - усредненные за время Т₁ показания акселерометров;

α_x,y,z,i - показания акселерометров на i-м шаге съема информации;

Δt - временной шаг съема информации.

Параметры ориентации A_z, α, β вычисляются по алгоритмам:

Матрица ориентации $${\tilde{D}}_0$$ через параметры (A_z, α, β) определяется в линейном приближении соотношением

По исходным данным, полученным в интервалах T₁, T₂, T₃, и воспользовавшись начальной ориентацией (14), проводят интегрирование уравнений навигации во всем интервале измерений от t=0 до t=T₁+T₂+T_3.

К результатам интегрирования применяют сглаживающий фильтр:

коэффициенты которого B_0x,z, B_1x,z, B_1y определяются на каждом шаге интегрирования рекуррентной процедурой метода наименьших квадратов (МНК).

Из соотношений (10) и (15) следует:

Соотношения (15) и (16) вычисляются в конце каждого временного интервала.

По ИД, полученным в интервале Т₃, после интервала ударных воздействий, выполняют повторную выставку ориентации по соотношениям (11)-(14) и определяют матрицу $${\tilde{D}}_3$$ .

Для расчета смещений нулей гироскопов и акселерометров до ударных воздействий, во время ударных воздействий и после ударных воздействий используют формулы (15) и (16).

Причем для расчета смещений нулей ЛГ и МА до ударных воздействии используют величины l_x, l_y, l_z, рассчитанные по уходам координат, и углов ориентации относительно начальной матрицы $${\tilde{D}}_0$$ участке T₁.

Для расчета смещений нулей ЛГ и МА после ударных воздействий используются величины l_x, l_y, l_z, рассчитанные по уходам координат, и углов ориентации относительно начальной матрицы $${\tilde{D}}_3$$ на участке Т₃.

Для расчета смещений нулей ЛГ и МА во время ударных воздействий используются величины l_x, l_y, l_z, рассчитанные по уходам координат, и углов ориентации относительно начальной матрицы $${\tilde{D}}_0$$ на полном участке Т₁+Т₂+Т₃ за вычетом уходов на участках T₁ и Т₃.

Все расчеты уходов производятся по модулю.

Экспериментальная проверка способа.

Исходные данные (ИД) при экспериментальной проверке способа записывали с использованием стенда, имеющего разметку для установления ИИБ (СВСК) по известному азимуту А=0,812842 рад.

На фиг.2, 3, 4 приведены графики временной зависимости угловых координат α, β, A, полученных интегрированием уравнения вращения по показаниям датчиков ЛГ. Из графиков видно, что в интервале T₂ имеют место угловые смещения координат ИИБ (т.е. СВСК относительно ГСК) по всем трем координатам. Это обусловлено тем, что ИИБ не закреплен жестко на стенде, и поэтому имеют место смещения изделия в момент ударных воздействий.

Вычисленные смещения, приведенные на графиках фиг.2, соответствуют реальным смещениям оптического элемента (ОЭ), оси которого определяют СВСК, относительно своего первоначального положения.

На фиг.5-10 приведены графики временной зависимости уходов линейных координат l_x, l_y, l_z, полученные интегрированием уравнения навигации по показаниям датчиков ЛГ и МА, с последующим использованием математических фильтров (15). На фиг.5-7 приведены графики, полученные по измерениям в интервале на фиг.8-10 по измерениям во всем интервале T₁+Т₂+Т₃.

Вычисления по соотношениям (16) приводят к результатам на интервале T₁:

на интервале T₃:

на полном интервале T₁+T₂+T₃:

Результаты, полученные по ИД в интервале Т₁+Т₂+Т₃, незначительно отличаются от результатов, полученных в интервале T₁.

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

Способ определения погрешностей инерциальных измерительных приборов при испытаниях на ударные и вибрационные воздействия, основанный на оценке погрешностей измерений по исходным данным, полученным по результатам измерений показаний инерциальных измерительных приборов, отличающийся тем, что измерение показаний инерциальных измерительных приборов производят в три этапа с одинаковой продолжительностью по времени при различной для каждого из этапов пространственной ориентации инерциальных измерительных приборов, причем на первом этапе производят измерение показаний инерциальных измерительных приборов для определения начальной ориентации и интегрирования уравнения навигации, на втором этапе производят измерение показаний инерциальных измерительных приборов для интегрирования уравнений навигации, когда инерциальные измерительные приборы подвергают ударным и/или вибрационным воздействиям, а на третьем этапе производят измерение показаний инерциальных измерительных приборов для начальной ориентации и интегрирования уравнения навигации после ударных и/или вибрационных воздействий, при этом для расчета смещений нулей инерциальных измерительных приборов до ударных воздействий используют величины l_x, l_y, l_z, рассчитанные по уходам координат и углов ориентации относительно начальной матрицы на участке T₁, для расчета смещений нулей инерциальных измерительных приборов после ударных воздействий используют величины l_x, l_y, l_z, рассчитанные по уходам координат и углов ориентации относительно начальной матрицы на участке T₃, а для расчета смещений нулей инерциальных измерительных приборов во время ударных воздействий используют величины l_x, l_y, l_z, рассчитанные по уходам координат и углов ориентации относительно начальной матрицы , где l_x, l_y, l_z, рассчитывают из соотношений l_x=B_0xt³+B_1xt², l_z=B_0zt³+B_1zt², l_y=B_1yt², коэффициенты B_0x,z, B_1x,z, B_1y определяют на каждом шаге интегрирования рекуррентной процедурой метода наименьших квадратов, а смещения нулей инерциальных измерительных приборов определяют по соотношениям
; Δsα_x=2·B_1x, ; Δsα_z=2·B_1z, Δsα_y=2·B_y'.