Модель ошибок гироскопа

Автореферат диссертации по теме «Модель погрешностей волоконно-оптического гироскопа»

На правах рукописи

Антонова Мария Вячеславовна

МОДЕЛЬ ПОГРЕШНОСТЕЙ ВОЛОКОННО-ОПТИЧЕСКОГО ГИРОСКОПА

05.11.03 — Приборы навигации

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

005545947

13 ПАР 2014

Москва-2014

005545947

Работа выполнена в Московском государственном техническом университете имени Н.Э. Баумана.

Научный руководитель:

Матвеев Валерий Александрович,

доктор технических наук, профессор

Официальные оппоненты:

Редькин Сергей Петрович,

доктор технических наук,

ЗАО «Инерциальные технологии

«Технокомплекса»

Колбас Юрий Юрьевич,

кандидат технических наук, ОАО «НИИ «Полюс» им. М.Ф. Стельмаха»

Ведущая организация: ОАО «Раменское приборостроительное

конструкторское бюро», Московская область

Защита диссертации состоится 16 апреля 2014 г. в 10 час. 00 мин. на заседании диссертационного совета Д212.141.19 в Московском государственном техническом университете имени Н.Э. Баумана по адресу: 105005, 2-я Бауманская ул., д. 5, стр. 1.

С диссертацией можно ознакомиться в библиотеке МГТУ им. Н.Э. Баумана.

Ваш отзыв в двух экземплярах, заверенных печатью организации, просим направлять по адресу: 105005, г. Москва, 2-я Бауманская ул., 5, стр. 1, МГТУ им. Н.Э. Баумана, учёному секретарю диссертационного совета Д212.141.19.

Автореферат разослан « С1 j » марта 2014 г.

Ученый секретарь диссертационного совета, к.ф.-м.н.

Семеренко Д.А.

Общая характеристика работы

Актуальность темы диссертации. Бесплатформенные командно-измерительные приборы перспективных систем управления для ракетно-космической техники (РКТ) строятся на базе инерциальных чувствительных элементов нового поколения, в первую очередь, на волоконно-оптических гироскопах (ВОГ).

Одной из основных задач на этапе разработки ВОГ является создание математической модели погрешностей (ММП) гироскопа, которая позволяет учесть систематические составляющие ошибок в условиях конкретных применений.

В основе ММП традиционных механических гироскопов, особенно в условиях их применения в ракетно-космической и в авиационной технике, лежит зависимость скорости ухода от перегрузки. Кроме того, ММП дополняется зависимостями её коэффициентов от температуры окружающей среды, напряженности магнитного поля и т.п. Что касается ВОГ, то помимо характерных для этого типа гироскопов источников погрешностей (шумы интенсивности источника излучения, электрооптический эффект Керра, дробовые шумы фотоприёмника и др.), в отечественной и зарубежной литературе имеются сведения о влиянии на выходной сигнал температуры и магнитных полей. В публикациях специалистов фирм «Оптолинк», «Физоптика», «ПНППК», ЦНИИ «Электроприбор», занимающихся разработкой ВОГ, затрагиваются проблемы снижения температурной и магнитной чувствительности ВОГ. При этом представленные в том или ином виде модели погрешностей ВОГ содержат зависимости скорости ухода только от этих факторов.

Информация о влиянии перегрузок на выходной сигнал ВОГ в отечественных публикациях отсутствует, в то время как в зарубежной литературе содержатся сведения о заметной чувствительности ВОГ к перегрузкам. Так, D.H. Titterton и J.L. Weston из Американского института аэронавтики и астронавтики приводят сводку точностных характеристик ВОГ и указывают на скорость ухода, пропорциональную перегрузке, порядка 1 (°/4)/g, и пропорциональную квадрату перегрузки, порядка 0,1 (74)/g2. В стандарте США IEEE Std 952-1997, регламентирующем перечень технических характеристик и методы испытаний ВОГ, чувствительность к линейным перегрузкам среди обязательных паспортизуемых параметров не значится, но указывается в качестве справочной характеристики в подразделе «Чувствительность к другим внешним факторам» для конкретных применений. При этом не описываются механизмы влияния перегрузки, и, тем более, не представляется модель таких погрешностей в аналитическом виде. Для применения ВОГ в системах управления средств выведения космических аппаратов и других объектов, испытывающих большие перегрузки, необходимо изучить влияние линейных ускорений на выходной сигнал ВОГ. Таким образом, создание ММП ВОГ, учитывающей влияние этого возмущающего фактора, является актуальной научно-технической задачей.

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

Для достижения указанной цели были поставлены и решены следующие задачи:

— проведён расчёт напряжений и деформаций типовой конструкции ВОГ при воздействии перегрузок;

— разработана тепловая модель ВОГ, позволяющая получить параметры поля температур элементов конструкции ВОГ;

— разработана ММП, учитывающая влияние на скорость ухода ВОГ линейных перегрузок, внешних магнитных полей, температуры и скорости её изменения;

— разработана методика экспериментального определения параметров ММП;

— проведена проверка адекватности ММП: при действии единичной перегрузки — на основе результатов многопозиционных испытаний ВОГ; при действии многократных перегрузок — на основе результатов испытаний на низкочастотной вибрации с амплитудами виброускорений до 10 g;

— разработана система конструкторско-технологических мер по снижению влияния внешних воздействий на выходной сигнал ВОГ.

Объектом исследования являются ВОГ с интегрально-оптическим фазовым модулятором и эрбиевым волоконным суперлюминисцентным источником излучения.

Предметом исследования является зависимость скорости ухода ВОГ от условий его эксплуатации.

Методы исследований. В диссертации использованы методы прикладной теории гироскопов, теории оптических систем, теории обработки сигналов, методы теории упругости, методы конечно-элементного анализа.

Научная новизна

1. Экспериментально доказана чувствительность выходного сигнала ВОГ к линейной перегрузке на уровне, существенном для применений в PKT.

2. Зависимость скорости ухода ВОГ от проекций перегрузки на оси ВОГ представлена в аналитической форме.

3. Установлен и экспериментально подтвержден факт наличия в ММП ВОГ нехарактерных для ММП традиционных механических гироскопов двух нелинейных составляющих, которые зависят от модуля проекции ускорения на плоскость, перпендикулярную оси чувствительности (04).

4. Дано объяснение одной из причин виброчувствительности ВОГ — возникновение при действии линейной вибрации под углом к ОЧ ВОГ

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

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

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

Реализация и внедрение результатов. Результаты, полученные в диссертации, использованы в ФГУП «ЦЭНКИ» — «НИИ ПМ им. академика В.И. Кузнецова» при создании волоконно-оптических гироскопов КИНД11-240, КИНД11-221, бесплатформенного инерциального блока КИНД34-059, что подтверждено актом внедрения.

Основные результаты и положения, выносимые на защиту:

1. На основе типовой конструкции ВОГ с учётом характеристик материалов её элементов построена расчетная модель ВОГ, позволяющая оценить параметры поля температур и напряжённо-деформированного состояния элементов ВОГ при характерных условиях его эксплуатации как факторы, влияющие на точностные характеристики гироскопа.

2. Разработана ММП ВОГ, учитывающая влияние на скорость ухода ВОГ линейных перегрузок, внешних магнитных полей, температуры и скорости её изменения.

3. Алгоритмическая компенсация выходного сигнала ВОГ с использованием разработанной ММП снижает погрешность определения входной угловой скорости в несколько раз.

4. Наличие в ММП ВОГ нелинейных составляющих, зависящих от модуля проекции ускорения на плоскость, перпендикулярную оси чувствительности, позволяет адекватно описать выходной сигнал ВОГ при действии на него низкочастотной вибрации с амплитудой до 10

Апробация работы. Основные результаты диссертации докладывались на XVIII международной конференции по интегрированным навигационным системам, г. Санкт-Петербург, 2011 г., и 4 всероссийских научно-технических конференциях: II Всероссийской научно-технической конференции «Актуальные проблемы ракетно-космической техники», г. Самара, 2011 г.; Первом Всероссийском космическом инновационном конвенте, г. Углегорск, 2011 г.; II научно-технической конференции «Системы управления беспилотными космическими и атмосферными летательными аппаратами», г. Москва, 2012 г.; Всероссийской конференции молодых специалистов, учёных

и студентов памяти Главного конструктора, академика АН СССР В.И. Кузнецова, г. Москва, 2013 г.

Публикации. Основные результаты диссертации опубликованы в пяти статьях, из которых три — в журналах, входящих в перечень ВАК.

Структура и объем работы. Диссертационная работа состоит из введения, четырёх глав, выводов, списка литературы и двух приложений. Работа изложена на 130 листах машинописного текста, содержит 49 рисунков и 31 таблицу. Список литературы содержит 61 библиографическую ссылку.

Содержание работы

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

В первой главе приводится модель погрешностей ВОГ, построенная исходя из общих соображений о принципе действия, структуре и конструкции прибора, которые сводятся к следующим положениям:

— первоисточником влияния массовых сил на показания ВОГ могут быть только вызываемые ими деформации составных частей ВОГ и, как следствие, напряжения в волокне чувствительной катушки волоконного кольцевого интерферометра и в модуляторе и обусловленные ими оптические эффекты;

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

— диаметральное расположение модулятора на торце катушки и наличие в катушке зоны входа/выхода оптоволокна в какой-то мере нарушают круговую симметрию ВОГ, что не позволяет отрицать полностью возможность зависимости показаний прибора от направления радиальной перегрузки;

— необходимо учесть эффект взаимовлияния деформаций катушки, вызванных осевой и радиальной составляющими «косой» перегрузки, т.е. появление смещения нуля ВОГ, пропорционального произведению этих проекций перегрузки, подобного скорости ухода механического роторного гироскопа из-за неравножесткости подвеса ротора.

Вид модели:

А<о~еот + й>„ ■axlg+at, •о,./г+й)в.-а,/£ + £»„ ■а.^+2еап.-а,-аг1ёг + + 2a>v■a,■ar/g2+2a¡„■a,■a!/gг-6K■mm,

где Дш-скорость ухода;

аг, ау, а. — проекции кажущегося ускорения корпуса ВОГ на его оси;

аг = + а) — радиальная составляющая кажущегося ускорения;

тг — составляющая скорости ухода, не зависящая от ускорения (перегрузки);

тр — составляющая скорости ухода, пропорциональная ускорению по ОЧ (по оси ОХ) при единичной перегрузке;

а>0- составляющая скорости ухода, пропорциональная модулю ускорения в плоскости, перпендикулярной ОЧ при единичной перегрузке;

со^, сав — составляющие скорости ухода, пропорциональные ускорениям по осям О У и ОЪ, соответственно при единичной перегрузке;

£У1Г — составляющая скорости ухода, пропорциональная произведению ускорений вдоль и перпендикулярно ОЧ, при единичной перегрузке, действующей под углом 45° к оси ОЧ;

со,,., <£>_ — составляющие скорости ухода, пропорциональные произведению ускорений по осям ОХ и ОУ, ОХ и ОЪ, соответственно при единичной перегрузке, действующей под углом 45° к оси ОХ в плоскости ОХУ,

охг;

ускорение свободного падения;

Ы— приращение выходного кода ВОГ за такт опроса;

ДК — отклонение масштабного коэффициента от паспортизованного значения;

11ПП.- интервал интегрирования входной угловой скорости, равный такту опроса.

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

Рассматриваются объекты исследования: приборы ВОГ типа КИНД11-221 (далее ВОГ-1), КИНД11-240 (далее ВОГ-2), КИНД11-222 (далее ВОГ-3), разработанные в филиале ФГУП «ЦЭНКИ» — «НИИ ПМ», их основные параметры и конструктивные особенности. Приводятся параметры испытательной аппаратуры и электроизмерительных приборов.

Во второй главе рассмотрены основные механизмы влияния температуры на ВОГ. На основании выражения для скорости ухода (2) рассматриваются способы уменьшения эффекта Шьюппа.

где Дп — малая поправка к основному показателю преломления гг0, постоянному вдоль световода, меняющаяся вдоль световода в зависимости от координаты 2 и момента времени /; Ь — длина волоконного контура; О — диаметр чувствительной катушки.

(2)

Приводится тепловая модель ВОГ, разработанная на основе эскизов ВОГ-1 и примененных материалов с использованием программы Solid Works, на основе следующих методических положений:

— геометрические параметры элементов конструкции тепловой модели соответствуют номинальным размерам элементов;

— геометрическое расположение элементов, порядок сборки и тепловая связь между элементами обеспечивалась командами программы SolidWorks;

— характеристики материалов катушек ВОГ задавались по справочным источникам на примененные материалы;

— тепловые характеристики элементов обмоток световода и компаунда заполнения намотки задавались в виде условного материала (эквивалента волокна);

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

В среде ANSYS Workbench проведена имитация распределения температуры по конструкции ЧЭ ВОГ в установившемся тепловом режиме (рабочая температура 55±2 °С; мощность, выделяемая линейным нагревательным элементом, расположенном на верхнем экране ВОГ — 5,5 Вт). Моделирование показало, что при заданной мощности в нагревательных элементах конструкция ВОГ прогревается осесимметрично до температуры 53,9 °С.

Приводится методика определения коэффициентов температурной чувствительности ВОГ. С использованием элементов системы терморегулирования (СТР) — нагревателей и термодатчиков, установленных на блоке катушек ВОГ, и блока усилителя системы терморегулирования задавались температуры термостатирования от 30 до 50 °С с шагом 5 «С. Дополнительно проведён запуск ВОГ при его разогреве от н.к.у. до рабочей температуры 55 °С. Одновременно записывалось сопротивление контрольного термодатчика. Для тепловой изоляции ВОГ помещен в алюминиевый кожух, между ВОГ и кожухом установлены прокладки из пенопласта. Конструкция ВОГ на рабочем месте для тепловых испытаний показана на рис. 1.

Крышка

Катушка внутренняя

Рис. 1. Конструкция типового ВОГ (разрез в плоскости OXZ)

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

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

Аа>т = сол +а)т-АТ+й)т!-АТ1+ солт -¿Т/сИ, (3)

где й>г0 [7ч] — составляющая скорости ухода, измеренная на калибровочных испытаниях при некоторой фиксированной температуре То;

Т— температура ВОГ (по термодатчику); АТ=Т-Т0;

а>т [(7ч)/°С] и а>т, [(7ч)/(°С)2] — коэффициенты влияния на скорость ухода

ВОГ первой и второй степени отклонения температуры ВОГ от ее значения при калибровке (Д7);

со^ [(7ч)/(°С/мин)] — коэффициент влияния на скорость ухода скорости изменения температуры ВОГ.

С использованием тепловой модели проведена оценка перепадов температур между разными зонами волокна внутри катушки при температурах разогрева конструкции, которые задавались в эксперименте. Для этого заданы разные значения мощности тепловыделения нагревательного элемента верхнего экрана ВОГ и определены при этом перепады температур волокна внутри катушки в радиальном направлении.

В третьей главе приводится анализ причин проявления эффекта Фарадея в волоконном кольцевом интерферометре и рассчитывается максимальная ошибка ВОГ-1, обусловленная этим эффектом. Рассмотрены методы уменьшения магнитной чувствительности. Рассчитан коэффициент магнитного экранирования для применённого в ВОГ-1 магнитного экрана из высоконикелевого пермаллоя марки 79 НМ.

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

Приведены результаты испытаний ВОГ-1 (в магнитных экранах) и ВОГ-2 (в магнитных экранах и без них) при воздействии магнитного поля последовательно по трём осям прибора, проекция вектора напряженности (Н) которого варьировалась от Н=-300 мкТл до Н=300 мкТл с шагом 50 мкТл. По (4) определены коэффициенты влияния магнитного поля сот, соНу, сот [°/ч],

определяющие скорость ухода ВОГ при Н=100 мкТл.

й^Да/ /6, (4)

где с/ — индекс оси; До/ = (й^,зо0-ю*_зоо) — разность средних значений угловой скорости, измеренных при воздействии магнитного поля с Н=300 мкТл и Н=-300 мкТл по оси ОХ, ОУ или ОЪ.

Проведённые исследования показали, что у двух испытанных типов ВОГ (ВОГ-1 и ВОГ-2) обнаружена заметная магнитная чувствительность. Применённое экранирование снижает чувствительность ВОГ к магнитному полю, но все равно она остаётся на значительном уровне (до 1,8 °/ч на ЮОмкТл у ВОГ-2).

В связи с этим предложено использование многослойных экранов. В тех же габаритах предлагается использовать трёхслойный экран, в котором в качестве материала промежуточного экрана выбран капролон (ПА6 блочный), который является диэлектриком и облегчает конструкцию. Применение такого экрана снизит магнитное поле, действующее на ВОГ, примерно в 29 раз (по оценке без учёта технологического зазора между нижним и верхним экранами).

По результатам экспериментальных данных разработана ММП, определяющая зависимость скорости ухода от внешних магнитных полей:

Дй>. = а„, •#,/#,+ соНу ■ Н,!НЯ + <аН! ■ Н,1НЕ, (5)

где Нх, Ну, Нг [мкТл] — проекции вектора напряженности магнитного поля на оси ВОГ;

Не — 100 мкТл — нормирующий множитель, близкий по порядку величины к напряженности магнитного поля Земли;

(ош, аИу, сон, [°/ч] — коэффициенты влияния магнитного поля, действующего по осям ОХ, ОУ, ОЪ, на скорость ухода ВОГ.

В четвертой главе приведен расчёт деформаций и напряжений в корпусе и оптическом волокне ВОГ при действии сил инерции, вызванных ускорением вдоль диаметра корпуса ВОГ, закрепленного по нижней торцевой поверхности. Расчёт проведен с использованием разработанной модели ВОГ-1 с помощью программного продукта АИБУв.

Приведена методика определения параметров модели погрешностей (1) путём проведения многопозиционных испытаний ВОГ в поле действия единичной перегрузки. ВОГ при этом устанавливался на плиту наклонно-поворотного стенда СНП-200, который обеспечивал выставку ВОГ в нужную ориентацию с точностью ±30″. Для определения параметров ММП достаточно выставить ВОГ в 8 позиций. Однако для повышения точности определения

составляющих погрешности прибор выставлялся в 16 позиций по отношению к направлению единичной перегрузки и направлению на Север (см. Таблицу 1).

В таблице 1 введены обозначения: «45° (И; -§)» — 04 ВОГ (ось ОХ) в позициях 13… 16 отклонена на угол 45° от направления на Север вверх; «45° (1чГ; +§)» — ось ОУ в позиции 13 и ось ОЪ в позиции 16 отклонены на угол 45 ° от направления на Север вниз; «45° (Э; -§)» — ось ОУ в позиции 15 и ось ОХ в позиции 14 отклонены на угол 450 от направления на Юг вверх.

Таблица 1.

№ п.п Ориентация осей ВОГ Иллюстрация

ОХ ОУ О Ъ

1 N W -ё N .у _ С

2 N +ё \у Л \У

3 N Е +g

4 N -ё Е

5 +ё N Е

6 +g XV N

7 +ё Б W

8 +ё Е Б

9 -В Е N к е

10 в Е

11 -ё N

12 -ё W 8

13 45° (И; -й 45° (Ы; +й W

14 45° (N5 -8) W 45° (в; -¿)

15 45° (К; -й 45° (Б; -й Е

16 45° (И; -й Е 45° (Ы; +й

Исходя из выражений для выходного сигнала ВОГ в соответствии с ММП в каждой позиции, в результате решения получившейся системы из 16 уравнений (с избыточностью) предложены следующие формулы для расчета параметров ММП:

о, = —

К (Л’6 + ЛГ8+#9+ЛГ12) 4 Т

-6?—

К (Лг’+ЛГ3 + //2+Л’4)

4 Т

к (лг-м1)

К (ЛГ14 + ^16) 0,. +а*+(!>„ +а,

-2Т—-Т2-

2Т 42′

К (Ын-Ы16)

а.

2 Т 42′

где N — приращение выходного кода ВОГ за такт опроса в ¡-ой позиции (/=1…16); Т — интервал интегрирования входной угловой скорости; К — масштабный коэффициент ВОГ, определённый при испытаниях на скоростном стенде.

Углы отклонения ОЧ оцениваются по формулам: К (Ы’-^+Ы9-^2)

п> ___^-1.

АТсо,

К (Л^-ЛГ’+ЛГ»-^10) ог = ——

Приведены результаты многопозиционных испытаний девяти приборов типа ВОГ-1 (выполнялись два — три запуска каждого из приборов). Запуски проводились в одних и тех же условиях без переустановки прибора. Получены параметры математической модели погрешностей в каждом запуске, а также средние значения параметров ММП и их отклонения в каждом запуске от среднего значения. Поскольку параметры й)^, сов, (Оуу, (Оп находятся на уровне нестабильности нулевого сигнала в запуске и их оценки из-за этого меняются случайным образом, то сохранять их в модели погрешностей, призванной учитывать систематические ошибки, нецелесообразно. Получена упрощённая ММП ВОГ:

&а> = сог+а>р-а,/g + a)l:r■arl g + 2colr^aI^ar/gг. (6)

Эта модель применена к имеющимся результатам испытаний. Рассчитанные (средние по нескольким запускам) параметры ММП для всех приборов, а также разброс ошибок (со -со,) выходного сигнала ВОГ в испытательных позициях и разброс ошибок после учёта ММП со -<оиш приведены в таблице 2. Величина ошибки (со — со3 ) рассчитана как разность выходного сигнала ВОГ и проекции угловой скорости вращения Земли на его ось чувствительности в соответствующей позиции; ( со — со иш ) — как разность выходного сигнала ВОГ и величины сигнала, рассчитанной с использованием ММП.

Таблица 2.

Параметры упрощённой ММП

прибора Параметры ММП —^ 001 002 003 004 005 006 007 009 09245

-1,830 -0,053 -0,739 -0,277 0,266 -0,183 0,395 -0,477 -0,061

, °/ч -0,871 -0,363 -1,196 -1,303 -1,290 -1,290 -1,213 -1,335 -1,209

«V , °/ч -0,962 -0,674 -1,184 -1,325 -1,158 -1,339 -1,321 -1,338 -0,846

«V, °/ч 0,229 0,532 0,423 0,513 0,536 0,520] 0,511 0,552 -0,004

9′, рад -0,001 -0,007 0,003 -0,002 0,006 0,000 -0,003 -0,002 0,005

3 > рад 0,004 -0,001 0,004 -0,001 0,004 -0,004 -0,001 0,001 0,002

0) — С03, °/ч от -2,86 -0,81 -2,18 -1,68 -1,13 -1,53 -1,01 -1,98 -1,59

до -0,93 0,34 0,63 1,10 1,66 1,12 1,64 1,04 1Д4

Ы-Ыммп* °/ч от -0,18 -0,05 -0,20 -0,08 0,04 -0,02 -0,05 -0,14 -0,07

до 0,18 0,05 0,23 0,09 0,15 0,02 0,03 0,25 0,01

Как видно из таблицы 2, учёт параметров ММП при расчёте измеренной угловой скорости уменьшает погрешность в определении входной угловой скорости для всех испытанных приборов в несколько раз. Кроме того, для определения четырёх параметров ММП (6) и углов отклонения оси чувствительности от её номинального положения достаточно шести уравнений и, соответственно, шести ориентаций прибора при испытаниях, что значительно сокращает время проведения калибровки по сравнению с использованием ММП (1).

Учитывая, что ММП (6) сформирована по результатам испытаний в поле действия единичной перегрузки, для проверки адекватности ММП в условиях многократных перегрузок, проведены испытания ВОГ на низкочастотной вибрации с виброперегрузками до 10 ед. Приводятся результаты испытаний при действии вибрации: 1) по ОЧ (синусоидальная вибрация на частоте 40 Гц с амплитудами виброперегрузки (Ам), равными 1; 2; 4; 6 и 8 ед.; на частоте 7 Гц с Ам = 0,5 ед. и на частоте 15 Гц с Ам = 1 ед.); 2) под углом 70° к ОЧ (синусоидальная вибрация на частоте 40 Гц с Ам = 1; 2; 4; 6 и 8 ед.); 3) перпендикулярно к ОЧ ВОГ (синусоидальная вибрация на частоте 15 Гц с Ам = 0,5; 1; 1,5; 2; 5; 8 и 10 ед.). Испытания показали, что при действии вибрации вдоль оси чувствительности, где «работает» линейная составляющая ММП Фр-а,!амплитуда переменной составляющей выходного сигнала ВОГ на частоте вибрации пропорциональна амплитуде виброперегрузки, а на постоянную составляющую вибрация практически не влияет. При действии

вибрации под углом к оси чувствительности, наличие в ММП нелинейной

составляющей а>„ -аг1ё (где аг = ^а’ + а]) приводит к появлению постоянной

составляющей скорости ухода ВОГ Д\Уср, линейно зависящей от перегрузки.

Так, при вибрации под углом 70° к

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

Рис. 2. Результаты виброиспытаний ВОГ-1 (1-вибрация по ОЧ; 2 — вибрация под углом 70° к ОЧ)

\ZV2f, ‘/ч

10 Ам, ед-

Рис. 3. Зависимость выходного сигнала

ВОГ-3 на двойной частоте вибровозмущений от виброперегрузки

Д\УС[

возрастает с ростом виброперегрузки от 2 до 8 ед. на 7,3 °/ч, тогда как при действии вибрации по ОЧ не превышает 0,6 °/ч при возрастании перегрузки в тех же пределах (см. рис. 2).

При вибрации по оси ОЪ, перпендикулярной к ОЧ, когда составляющие ММП сор-аг/ё+ +2а>„-ах-аг1 £ обнуляются и «работает» только апеременная составляющая выходного сигнала ВОГ изменяется по гармоническому закону на частоте, вдвое превышающей частоту вибровозмущения (см. рис. 4), причём амплитуда этой составляющей (^¡гЬ растёт пропорционально виброперегрузке (см. рис. 3).

Т.о., проведенные испытания полностью подтверждают разработанную ММП (6). Это свидетельствует о том, что ВОГ реагирует на линейные движения основания по отражённым в этой ММП зависимостям.

Механизмы влияния

перегрузки на ВОГ требуют дополнительных исследований, в

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

Рис. 4. Результаты испытаний ВОГ-3 при действии вибрации по оси ОЪ на частоте 15 Гц с различными амплитудами виброперегрузки: а) 2 ед.; б) 8 ед.

По разработанной ММП (6) сделаны следующие выводы:

— отсутствие в (6) параметров со^, <Оа, > ^ говорит о том, что направление радиальной перегрузки не сказывается значимым образом на величине вызываемой ею деформации катушки, несмотря на нарушающие круговую симметрию особенности конструкции ВОГ (диаметральное расположение модулятора на торце катушки и наличие в катушке зоны входа/выхода оптоволокна);

— ММП содержит линейную составляющую, зависящую от ускорения по оси чувствительности, и нелинейные составляющие, нехарактерные для традиционных механических гироскопов, зависящие от модуля ускорения в плоскости, перпендикулярной оси чувствительности (&>а.), и произведения ускорений вдоль и перпендикулярно оси чувствительности X (¿«V). Наличие этих составляющих в модели было подтверждено как в поле действия единичной перегрузки, так и при многократных вибрационных перегрузках;

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

С учётом влияния на скорость ухода ВОГ всех трёх рассмотренных видов возмущающих факторов, интегральная математическая модель погрешностей ВОГ записана в следующем виде:

Ав> = Фг, + сот ■ ДГ + й),., • ДГ2 + • г?Г / Л + ®„,-Я,/Яя + + а>И!-Ну/Нк +о)шНУНЕ + а>е-аг/£ + й>„ • о, /£ + 2а>, • а, ■ аг / ёг — АК ■ ЛИг„

где АО [°/ч] — скорость ухода;

сол [°/ч] — составляющая скорости ухода, не зависящая от ускорения и от внешнего магнитного поля, измеренная на калибровочных испытаниях при некоторой фиксированной температуре То;

Т— температура ВОГ (по термодатчику); АТ=Т-Т„;

0г[(7ч)ЛС] и агг [(°/ч)/(°С)2] — коэффициенты влияния на скорость ухода ВОГ первой и второй степени отклонения температуры ВОГ от ее значения при калибровке (Д Т)\

соат [(7ч)/(°С/мин)] — коэффициент влияния на скорость ухода скорости

изменения температуры ВОГ;

Нх, Ну, Нг [мкТл] -проекции вектора напряженности магнитного поля на оси ВОГ;

Не — 100 мкТл — нормирующий множитель;

соНх, еаИг, фИх [7ч] — коэффициенты влияния магнитного поля,

действующего по осям ОХ, ОУ, ОХ, на скорость ухода ВОГ;

ах, ау, а, [м/с2] — проекции кажущегося ускорения корпуса ВОГ на его

оси;

аг = (а* + а1.)’2 [м/с2] — радиальная составляющая кажущегося ускорения; g [м/с2] — ускорение свободного падения;

°>р [7ч] — составляющая скорости ухода, пропорциональная ускорению по ОЧ, при единичной перегрузке;

(о^ [7ч] — составляющая скорости ухода, пропорциональная модулю ускорения в плоскости, перпендикулярной ОЧ, при единичной перегрузке;

а>„ [7ч] — составляющая скорости ухода, пропорциональная произведению ускорений вдоль и перпендикулярно ОЧ, при единичной перегрузке, действующей под углом 45° к оси ОЧ;

N [е.м.р.] — приращение выходного кода ВОГ за такт опроса; ДК [«/е.м.р.] — отклонение масштабного коэффициента от паспортизованного значения;

1,„Лс] — интервал интегрирования входной угловой скорости, равный такту опроса.

Основные результаты и выводы диссертационной работы

1. Разработана методика исследования зависимости скорости ухода ВОГ от внешних факторов — линейных перегрузок, магнитных полей и температуры, использование которой позволяет определить параметры ММП.

2. Построена расчетная модель ВОГ, которая позволяет оценить параметры поля температур и напряжённо-деформированного состояния элементов ВОГ при характерных условиях эксплуатации как факторы, влияющие на точностные характеристики гироскопа.

3. Экспериментально подтверждена зависимость выходного сигнала ВОГ от проекций линейной перегрузки на его оси. У девяти испытанных в поле силы тяжести образцов ВОГ составляющие скорости ухода, зависящие от ориентации относительно вертикали, достигают 1,9 °/ч. Такой уровень погрешностей является существенным для применения ВОГ в изделиях PKT.

4. Разработана ММП ВОГ, учитывающая влияние на скорость ухода ВОГ линейных перегрузок, внешних магнитных полей, температуры и скорости её изменения. Алгоритмическая компенсация составляющих ММП с использованием информации от акселерометров, магнитометров и термодатчиков уменьшает погрешность определения входной угловой скорости в несколько раз.

5. В ММП, помимо составляющей, пропорциональной ускорению вдоль оси чувствительности ВОГ, входят нехарактерные для других типов гироскопов нелинейные составляющие, зависящие от модуля ускорения в плоскости, перпендикулярной оси чувствительности. Значимый уровень этих составляющих подтвержден результатами испытаний ВОГ не только в поле силы тяжести, но и при действии вибрационных перегрузок, по крайней мере, до 10 ед.

6. Дано объяснение одной из причин виброчувствительности ВОГ -возникновение при действии линейной вибрации под углом к 04 ВОГ постоянной составляющей скорости ухода, сопровождающейся составляющей на удвоенной частоте вибрации, вызвано наличием в ММП нелинейных составляющих. При этом обе эти составляющие скорости ухода пропорциональны амплитуде виброперегрузки, что подтверждено экспериментально.

7. Разработана система конструкторско-технологических мер по снижению влияния внешних воздействий на выходной сигнал ВОГ.

Публикации по теме диссертации

1. Исследование характеристик бесплатформенного инерциального блока на базе волоконно-оптических гироскопов в процессе наземной отработки / М.В. Антонова [и др.] // Гироскопия и навигация. 2012. №2 (77). С. 98-111.

2. Антонова М.В. Разработка системы термостатирования для волоконно-оптического гироскопа // Вопросы оборонной техники. Серия 9. 2010. Специальный выпуск 3(244) — 4(245). С.163-171.

3. М.В. Антонова, A.B. Корнюхин. Виброустойчивость бесплатформенного инерциального блока на базе волоконно-оптических гироскопов // Вестник МГТУ, серия «Приборостроение», специальный выпуск №5 «Информатика и системы управления». 2012. С. 210-215.

4. М.В. Антонова, A.B. Корнюхин. Повышение виброустойчивости бесплатформенного инерциального блока на базе волоконно-оптических гироскопов // Труды «РКК «Энергия» им. С.П. Королёва» «Расчёт,

проектирование, конструирование и испытания космических систем». 2012. Серия XII выпуск №4. С.45-49.

5. Экспериментальное исследование зависимостей показаний волоконно-оптического гироскопа от его ориентации относительно вертикали: Доклад / М.В. Антонова [и др.] // Научно-технический сборник «Ракетно-космическая техника». Серия XI. Системы управления ракетных комплексов.

2010. Вып.2: в 3-х частях — Часть I. С. 76-84.

6. Работы по уточнению математической модели погрешностей волоконно-оптического гироскопа / М.В. Антонова [и др.] // VIII научно-практическая конференция молодых специалистов и студентов памяти главного конструктора академика В.И. Кузнецова: Сборник докладов. 2010. С. 35-44.

7. Бесплатформенный инерциальный блок на базе волоконно-оптических гироскопов в составе навигационного прибора: наземная отработка и подготовка к лётным испытаниям / М.В. Антонова [и др.] // П Всероссийская научно-техническая конференция «Актуальные проблемы ракетно-космической техники» («II Козловские чтения»): Сборник материалов.

2011. С.377-378.

8. М.В. Антонова, Л.З. Новиков. О некоторых результатах экспериментальной отработки малогабаритного волоконно-оптического гироскопа: Тезисы докл. // II Всероссийская научно-техническая конференция «Системы управления беспилотными космическими и атмосферными летательными аппаратами. 2012. С.81-82.

9. М.В. Антонова, Л.З. Новиков. Результаты испытаний волоконно-оптического гироскопа: Тезисы докл. // Актуальные проблемы российской космонавтики: Труды XXXVII академических чтений по космонавтике. 2013. С.628-629.

10. Некоторые результаты экспериментальной отработки волоконно-оптического гироскопа / М.В. Антонова [и др.] // Всероссийская конференция молодых специалистов, учёных и студентов памяти Главного конструктора, академика АН СССР В.И. Кузнецова: Сборник материалов. 2013. С. 47.

11. Волоконно-оптические гироскопы и приборы на их основе разработки филиала ФГУП ЦЭНКИ» — «НИИ прикладной механики им. академика В.И.Кузнецова» / М.В. Антонова [и др.] // Научно-техническая конференция к 75-летию кафедры «Приборы и системы ориентации, стабилизации и навигации»: Сборник материалов. 2013. С. 18-21.

Подписано к печати 3.03.14. Заказ №123 Объем 1,0 печ.л. Тираж 100 экз. Типография МГТУ им. Н.Э. Баумана 105005, Москва, 2-я Бауманская ул., д.5,стр.1 (499) 263-62-01

@inproceedings{VNovikov2017,
  title={Математическая модель ошибок лазерного гироскопа // Mathematical Error Model of the Laser Gyro},
  author={V. V.Novikov and В. В.Новиков},
  year={2017},
  url={https://api.semanticscholar.org/CorpusID:125866126}
}

In the article presents analysis and systematization of the results of experimental studies of laser gyro (LG) errors is performed. The structure of result LG error in the form of a linear combination of random processes characterizing the natural and technical fluctuations of the difference frequency of counter waves is determined in the presence of a random constant displacement of zero in the sensor readings. The requirements for the structure and form of the analytical representation of the… 

Наука и Образование

МГТУ им. Н.Э. Баумана

Наука и Образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2017. № 07. С. 189-205.

]Э5М 1994-040В

Б01: 10.7463/0717.0001271

Представлена в редакцию: Исправлена:

© МГТУ им. Н.Э. Баумана

УДК 004.942

Математическая модель ошибок лазерного гироскопа

Новиков В.В.1, Енин В.Н

2

Новиков П.В.

1,*

13.06.2017 27.06.2017

етп^ЬтЫхии

:МГТУ им. Н.Э. Баумана, Москва, Россия 2Московский Политехнический Университет,

Москва, Россия

Выполнен анализ и систематизация результатов экспериментальных исследований погрешностей лазерного гироскопа (ЛГ). Определена структура результирующей погрешности ЛГ в форме линейной комбинации случайных процессов, характеризующих естественные и технические флуктуации разностной частоты встречных волн, при наличии случайного постоянного смещения нуля в показаниях датчика. Сформулированы требования, предъявляемые к структуре и форме представления аналитического описания модели ошибок.Получена обобщенная модель флуктуационных процессов ЛГ, на основе которой разработана математическая модель ошибок ЛГ, как инерциального датчика. Модель представлена системой стохастических дифференциальных уравнений и функциональных соотношений, характеризующих результирующую ошибку датчика. Выполнен корреляционный анализ уравнений модели, получены конечные уравнения для среднеквадратичных значений отдельных компонент, позволяющие идентифицировать параметры результирующей ошибки. Параметры модели представлены через значения спектральной плотности мощности отдельных компонент. Рассмотрена дискретная форма модели, показана сходимость непрерывных и разностных уравнений, при выполнении условий предельного перехода. Определены направления дальнейших исследований.

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

Введение

Современный этап развития навигационных технологий открывает широкие возможности для динамичного развития целого ряда направлений данной отрасли. К наиболее перспективным областям отрасли, где достигнут существенный прогресс, следует отнести лазерную гироскопию [1, 8, 9]. Измерительные системы, приборы и отдельные компоненты лазерной гироскопии заняли прочное место, и нашли широкое применение в современных системах управления и навигации.

В конечных технических приложениях лазерных навигационных технологий в качестве инерциального датчика используется лазерный гироскоп (ЛГ). Неоспоримые функциональные и эксплуатационные преимущества ЛГ [1, 8] по сравнению с другими типами гироскопов, наиболее полно проявляются в бесплатформенных инерциальных навигационных системах (БИНС) [2, 4, 5, 7]. Стадия практического применения ЛГ в качестве датчика инерциальной информации БИНС выдвинула в число актуальных двуединую задачу — совершенствование точностных характеристик ЛГ в сочетании с разработкой методов и алгоритмов оптимизации точностных параметров БИНС, построенных на базе ЛГ [10-17].

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

Особую актуальность задача разработки модели ошибок приобретает в связи с появлением на рынке навигационных технологий интегрированных (инерциально- спутниковых) навигационных систем [3,4]. В интегрированных системах инерциальный модуль, выполненный по схеме БИНС, может оснащаться грубыми, малогабаритными ЛГ [1], что требует разработки алгоритмов и методов оптимизации точностных характеристик навигационной системы, требующих определения математической модели ошибок чувствительных элементов.

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

Анализу влияния дестабилизирующих факторов на точностные характеристики ЛГ посвящено значительное число работ, основополагающими из которых являются [9-17]. В большинстве из указанных работ представлены модели ошибок датчика в виде сложных эмпирических зависимостей далеких от задач инженерного анализа точности навигационных систем и мало пригодных для использования современных методов и алгоритмов обработки инерциальных данных [2-6].

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

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

1. Общая характеристика и обобщенная математическая модель

флуктуационных процессов в ЛГ.

Лазерный гироскоп является интегрирующим измерительным устройством с дискретным выходом, измеряющим угол поворота резонатора кольцевого лазера в инерци-альном пространстве за время съема информации [4-7], что позволяет рассматривать ЛГ в качестве датчика угла или (и) угловой скорости. Ошибка (результирующая погрешность, дрейф, уход) определения ЛГ измеряемых параметров порождается флуктуациями разностной частоты, либо разности фаз встречных волн резонатора лазера, неподвижного в абсолютном пространстве.

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

• адекватность предлагаемой математической модели реальному физическому процессу;

• полноту и наглядность математического описания модели;

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

Результаты исследования флуктуационных процессов ЛГ представлены в целом ряде работ, основополагающими из которых являются [6-17]. В основе исследований лежит детальный анализ природы проявления совокупности физически разнородных флуктуаци-онных процессов, порождающих результирующую ошибку в показаниях ЛГ. В соответствии с результатами исследований в спектре флуктуаций разностной частоты встречных волн ЛГ можно выделить две области, характеризующие технические и естественные флуктуации.

Технические флуктуации сосредоточены в низкочастотной части (полоса частот / до

2 1

10- -10- Гц.) спектра, где они превалируют по интенсивности над естественными флук-туациями. Экспериментально показано, что распределение спектральных составляющих технических флуктуаций подчиняется закону близкому к зависимости / 1.

Естественные флуктуации доминируют в высокочастотной области (полоса частот от 10-1 Гц) спектра флуктуаций, имея равномерный характер распределения интенсивности [6, 7, 8-13], при этом уровень интенсивности естественных флуктуаций следует рассматривать как величину адекватную порогу реагирования датчика

Согласно [6, 7, 8-13] основным источником естественных флуктуаций является спонтанное излучение атомов активной среды постоянной интенсивности — процесс с быстро меняющимися величинами, значения которых практически независимы и разделены

весьма малыми промежутками времени. Наиболее реалистичным математическим представлением процесса, обладающего постоянной интенсивностью и малым временем корреляции между отдельными измерениями является случайный процесс типа белого шума (White Noise (WN)).

Подтверждением правомерности математического представления естественных флуктуаций в виде WN служат результаты экспериментальных исследований, представленные в работе [15]. Здесь естественные флуктуации трактуются, как стационарный случайный процесс в эффективной полосе пропускания ЛГ с постоянным уровнем спектральной плотности мощности (Power Spectral Density (PSD)) в пределах 0,016-0,04 2 /.

(град/час) /Гц.

Вопросы количественной оценки интенсивности девиаций разности фаз, обусловленных естественными флуктуациями, рассматривались в работах [6, 7, 14-17]. Было показано, что изменение дисперсии разности фаз за время Т , подчиняется диффузионному закону M{а2ф}= т , где ПАФ коэффициент диффузии, характеризующий уровень

PSD разностной частоты встречных волн. Равномерный закон изменения спектра флук-туаций разностной частоты, в сочетании с линейным законом изменения во времени значений дисперсии разности фаз, является подтверждением на физическом уровне обоснованности математического представления естественных флуктуаций случайным процессом типа WN.

Если естественные флуктуации представляют нестабильность частоты биений ЛГ в виде некоррелированных, быстро изменяющихся величин, то технические флуктуации определяют большие по величине, но сравнительно медленные, хаотические уходы частоты [12, 13].

Источником технических флуктуаций в соответствии с [6, 7] являются нестабильности накачки (тока разряда), длины резонатора, случайные изменения температуры и градиента температур по периметру резонатора. Действия опосредованно связанных мультипликативных помех практически всех известных типов может быть сведено к возмущениям указанных параметров.

Следует отметить, что характер технических флуктуаций зависит от продолжительности интервала измерений, принимая за начало отсчета момент включения прибора. Формат статьи объективно не позволяет свести воедино результаты анализа процесса ухода на различных по продолжительности интервалах наблюдения, поэтому в данной статье рассматривается модель ошибок соответствующая ограниченному интервалу наблюдений (<3 часов).

В соответствии с результатами экспериментальных исследований [12-15], уровень PSD технических флуктуаций внутри временного интервала выхода на стабильный режим работы изменяется обратно пропорционально частоте наблюдения по закону близкому к зависимости (f ). В теории случайных процессов отсутствует строгое обоснование функционального соответствия между спектром f п процесса и его математическим представлением. Приближением к строгому решению указанной задачи может являться

подход, основанный на дробном интегрировании спектра WN [19,20]. В соответствии с данным подходом, «n» кратное интегрирование случайного процесса w(t) типа WN,

имеющего значение PSD равное —w , приводит к появлению случайного процесса v(t) с

PSD величиной Sv (с) = —w . В частности, если n=1, то v(t) следует рассматривать как

ю

случайный процесс на выходе интегрирующего звена с передаточной функцией H (p) = —

Р

, поэтому (= . Такой случайный процесс называется винеровским случайным

v ю

процессом (Wiener Random Progess (WRP)), или случайным блужданием (Random Walk(RW)). Свойства WRP задаются соотношениями [6]:

t

v(t) = J w{r)dr ; R (t, r) = Y min(t, r); Y min(t, r) = •

0

yt t ^r; y

; П v (ю) = -г,

Yt t< t. с

где — w(t) и v(t) соответственно случайные процессы типа WN и WRP, Rv(t,T) — корреляционная функция WRP, Пv (с) — значение PSD процесса v(t), порождаемого процессом

w(t) интенсивностью Y . Значения PSD процесса типа WRP убывают обратно пропорционально частоте наблюдения, а дисперсия процесса возрастает линейно во времени.

В соответствии с результатами экспериментальных исследований [12-15], низкочастотная область спектра флуктуаций разностной частоты аппроксимируется зависимостью —

Sv (ю) =—w , а характер изменения дисперсии флуктуаций определяется наличием прямо-

v с

линейного участка с положительным наклоном на кривой распределения Аллана (Allan Varians) [18], что позволяет обоснованно рассматривать технические флуктуации в качестве случайного процесса типа WRP.

Очевидно, что WRP пригоден для описания флуктуаций на ограниченном интервале наблюдения, так как беспредельное возрастание дисперсии флуктуаций противоречит реальным физическим законам.

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

С целью физически ясного и количественно адекватного представления обобщенной модели флуктуаций необходимо ввести в рассмотрение параметры ухода ЛГ, которые ранее не использовались в теории и практике лазерной гироскопии. К таким параметрам следует отнести — скорость изменения значений технических флуктуаций разностной частоты (rate of frequency changing (RFC)) и ускорение собственного ухода ЛГ. Под RFC следует понимать отношение приращения величины технических флуктуаций за расчетный период к длительности расчетного периода ([Гц/с]). Ускорение собственного ухода ЛГ ([град/с ]) необходимо рассматривать в качестве расчетного эквивалента RFC.

Формализованное описание флуктуационных процессов ЛГ в форме обобщенной математической модели поясняет структурная схема модели, представленная на рис.1.

Рис.1. Структурная схема обобщенной модели флуктуаций разностной частоты

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

• X(t) и v(t) — некоррелированные шумовые процессы типа WN, причем v(t) следует рассматривать, как процесс характеризующий RFC ( [ [Гц/с]) ;

• ç(t) и Ç(t) — шумовые процессы типа WRP ([Гц с] и [Гц]) ;

• y(t ) — случайный процесс изменения разности фаз, порождаемый техническими флуктуациями (Ç(t) ([Гц с]);

• S и G — PSD соответственно процессов A(t) и v(t) ([(Гц)2/Гц] и [Гц/с2)]);

• а — частота наблюдения флуктуаций ([рад/с]).

Полученная обобщенная модель флуктуаций позволяет характеризовать количественно результирующую ошибку ЛГ, используя PSD отдельных составляющих G /а4 и S /а2, что является практически важным результатом.

Наряду с флуктуационными составляющими в результирующем уходе ЛГ присутствует случайная s{t ), но постоянная составляющая, которая меняет свое значение от включения к включению прибора и трактуется как сдвиг нуля в показаниях ЛГ от запуска к запуску прибора [2,3,5-7]. Сдвиг нуля может быть описан случайной величиной с нулевым

средним и дисперсией m2. Сдвиг нуля в показаниях датчика подчиняется очевидному

дифференциальному уравнению s{t ) = 0.

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

Нестабильность выходного сигнала ЛГ следует рассматривать в виде линейной комбинации нескольких случайных процессов — случайного постоянного смещения нуля, естественных и технических флуктуации разностной частоты.

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

Спектр технических флуктуаций сосредоточен в области низких частотах наблюдения, быстро спадая по величине с ростом последних по закону — / 2.

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

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

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

Структурная схема и характеристики математической модели ошибок ЛГ

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

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

Знание математической модели позволяет получать конечные (не содержащие операторов дифференцирования) уравнения для определения параметров модели, что является основной целью исследования погрешностей ЛГ.

Изложение материала раздела содержит две взаимосвязанные части: представление ухода ЛГ непрерывным процессом и в форме дискретного аналога.

Рассмотрим модель ошибок в виде совокупности непрерывных процессов пространства входных (выходных) сигналов и вектора состояния. Структурная схема модели изображена на рис. 2.

™ ( * )

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

и (г) Х2(г)

¥ и ¥

а

4 ‘ 2

а а

Рис. 2. Структурная схема математической модели ошибок.

В структурной схеме модели введены следующие обозначения:

• и (V) и w(t) — случайные, некоррелированные случайные процессы типа WN;

• и ^ — интенсивности процессов и(V) и w(t) соответственно;

• х (V) — угол ухода, х2 (V) — случайное блуждание угловой скорости, обусловленное ускорением дрейфа — и (V);

• х3 (V) — систематическая погрешность смещения нуля;

• С — частота наблюдения.

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

Х(? ) = АХ (V) + BW(t). (1)

Здесь X () = [х ()х2 ()х3 ()] — вектор состояния, Wг (V) =

w(t) u(t) 0

— случайный

вектор шумов , А — постоянная матрица системы А =

0 1 1 0 0 0 0 0 0

, В — единичная матрица.

Компонентами вектора состояния х(/) являются: х (V) — угол собственного ухода ЛГ; х2 (?) — процесс типа WRP; х3 (V) — случайная, но постоянная величина, характеризующая смещение нуля от запуска к запуску в показаниях прибора. Значение смещения нуля меняется по величине и знаку от включения к включению прибора с распределением близком к нормальному [6,7]. Уравнение этой ошибки имеет очевидное представление х3 (V ) = 0, М X (/)}= т2, где т2 — среднеквадратичное значение случайной величины смещения нуля.

Процесс w(t) содержит в качестве составляющих некоррелированные белые гауссовы шумы w(t) и и(г), так что М) и()| = 0.

Вектор математического ожидания случайного векторного процесса w(t) имеет нулевое среднее значение

М ^(г)} = 0 и корреляционную матрицу

М ^ (г)}= Чд^ — г)

где Ч — неотрицательно определенная матрица интенсивности, £*(•) — дельта — функция Дирака.

Матрица интенсивности Т имеет вид:

Т

0 0 0 ¥ 0 0 0 0

Элементами матрицы Т служат интенсивности и ¥и шумовых процессов )

()

и Щ1) соответственно.

Определим первый и второй моменты для процессов скорости х 1 (/) и угла х1 (г) ухода ЛГ. Выражения для вектора математического ожидания м |х(/)} случайного процесса х(/), заданного стохастическим дифференциальным уравнением (1) в непрерывном варианте имеет вид [4, 5]

Ш

м {х(г)} = ам {х(г)}+вм )}.

(2)

Учитывая, что случайный вектор шума w(t) имеет нулевое среднее м{w(í)}= 0, дифференциальное уравнение (2) распадается на систему скалярных уравнений:

Шм {х1 (г )}=м {х2 (г)}+м {хз (г)}; ш

йм {Х2 (г )}= 0; ш

Шм {хз (г )}= 0.

ш

(3)

Ранее было отмечено, что смещение нуля х3 (?) в исходной реализации (при данном включении датчика) следует рассматривать, как случайную, постоянную величину с законом распределения близком к нормальному [4-6], поэтому осреднение по совокупности реализаций дает м {хз (г)} = 0.

Перейдем к определению среднеквадратических значений скорости м -|х 12 (г)| и угла м {х12 (г)} ухода ЛГ. Уравнение для определения матрицы среднеквадратических значений Vx = м {х(г)хт (г)} векторного случайного процесса х(г) имеет вид [4, 5]

V х (г) = аух (г)+V (г)Ат + вт(г )вт. (4)

Запись уравнения (4) в поэлементной форме приводит к системе дифференциальных уравнений:

dt

dM {x

dt

dM {x

dt

d_ dt

d dt

d dt

dM {x-dt x

dM {x2 dt x

dM (t )}= о

)}= 2[m(xi (t)}x2(t)}+M(xi (t)}хз(t)}]+ ¥wS(t); MX}= Щч3 + m2t2 +

)Ж (t )}=M X2 (t)}+M X (t )Ж (t)};

)Ж (t )}=M {x2 (t)}+M {хз (t )Ж (t)}; )fe (t)}= 0 ;

)}= wus(t);

Ш

M fe(t )X2 (t )}=y- t2; M {x (t)x3 (t)}= m2t ; M {x2 (t )x3 (t )}= 0 ;

MK2(t)}=Ш/ ;

M fe (t)}

=m

Процессы х2 (г) и х3 (г) статистически независимы так, что М {х2 (г )х3 (г)}= 0 .

Это равенство вполне обосновано, поскольку случайное смещение и процесс типа ЖЕР порождаются физически разнородными источниками погрешностей [6, 7].

Относительно начальных значений переменных, входящих в систему уравнений (5) физически обоснованными являются соотношения:

М {х2 (0)} = М {х2 (0)хз (0)} = М {х (0)хз (0)} = М {х (о)х2 (о)} = 0.

Допущение М{х2(0)}= 0 вытекает из свойств процесса типа ЖЕР [4], а допущение М {х2 (0)х (0)} = 0 обусловлено статистической независимостью процессов х2 (г) и х3 (г). Предположения м {х (0)х3 (0)}= 0 и м {х: (0)х2 (0)}= 0 диктуются нулевыми начальными условиями по углу ухода х (0), при наблюдении процесса ухода от момента включения прибора.

Учитывая, что xi (t) = x2 (t) + x3 (t) + w(t)

уравнения для определения среднеквадратич-

ных значений угловой скорости и угла ухода ЛГ соответственно принимают вид:

M

xi

(t)

= ¥ut + ¥wS(t)+m2, M {xi (t )} = Ши t3 + m2t2 + Wwt.

(6)

Важность полученных результатов с практической точки зрения состоит в возможности получения количественной оценки результирующей погрешности ЛГ через значения PSD отдельных компонент математической модели.

2. Математическая модель ошибок ЛГ в дискретной форме

Цифровая форма выходного сигнала ЛГ диктует необходимость представления полученных ранее уравнений в дискретной форме. Непрерывной форме представления динамической системы (1) соответствует векторное разностное уравнение первого порядка

[4]

x[(k + l)T] = ФX(kT)+:rW(kr) .

(7)

2

»

Матрицы Ф и Г , выражаются через матрицу коэффициентов А, единичную матрицу Е размером (з х з) и интервал квантования Т уравнениями вида Ф = Е + ТА и

Г = ТЕ [5, 6]. Следует отметить, что интервал квантования Т , отделяющий (к +1) и (к)

такты проделанных измерений, определяется как величина кратная периоду колебаний частотной подставки ЛГ [6].

Дискретный аналог непрерывного векторного дифференциального уравнения (3) записывается следующим образом [1-3]

V [(к + 1)Т] = Ф^ [кТ]ФТ + ГО[кТ]ГТ . (8)

Здесь В[кТ] — неотрицательно-определенная матрица интенсивности дискретного шумового процесса w[xт ], имеющая вид

^Д., 0 0″

м |>Т]} = ъ[кгдк [(* — } Т ]; в[кТ]=

0 д

0

и 0 0 0

где и Ди — интенсивности гауссовых <<белых>> последовательностей w(kT) и

и(кТ) соответственно. Функция Кронекера бк [(I — j)Т] имеет общепринятое обозначе-

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

ние

Ъ [О — j )Т ] =

11 при I = } [0 при I ф ]

Запись уравнения (8) в поэлементной форме приводит к системе рекуррентных уравнений:

м ^х12 [(к + 1)Т ]}=м к2 [кТ]}+ 2Т

м {х [кТ х [кТ ]}+м {х [кТ]х3 [кТ]}

+

2Т2м х2 [кТ]х3 [кТ]}+ Т2

+ Т2 Д.,

(9)

значений

м {х2 2 [кТ]}+ м {с,2 [кТ]}

м {х [(к + 1)т х [(к + 1)т ]}=м {х [кТ [кТ ]}+ ТЫ {х2 2 [кТ ]}+ ТЫ {х3 [кТ^2 [кТ]}; м х [(к + 1)Т Зх, [(к + 1)Т ]} = м ^ [кТ ]х, [кТ ]}+ ТИ {х3 2 [кТ]}+ ТИ {х3 [кТ ]^2 [кТ]}; м{х22 [(к + 1)Т]}= м{х22 [кТ]}+ Т2Д ;

м х2 [(к + 1)Т ]} = м {^з2 [кТ ]}. Рассмотрим последовательно каждое из уравнений системы. В силу неизменности » последовательности х3 [кТ], выражение для м{х-2 [кТ]} будет иметь тот же вид,

у м {хз2 [кТ ]}= т2

что и в непрерывном варианте записи, поэтом

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

м {х22 [(к + 1)Т ]}, м Ы(к + 1)Т ]х2 [(к + 1)Т ]}, м {х! [(к + 1)Т [(к + 1)Т ]}, м{х2 [(к + 1)Т х [(к + 1)Т]},

а именно:

МК2 [(к + 1)т]}= (к + 1)т2Би; М{х[(к + 1)т]к [(к + 1)т]}= ^ т3Д,;

М{к [(к + 1)тК [(к + 1)т]}= (к + 1)Тт2 ; М{х2 [(к + 1)тК [(к + 1)т]}= 0.

Подставляя полученные соотношения в первое уравнение системы (9) окончательное выражение для определения угла ухода ЛГ в дискретной форме примет вид

М{к2 [(к + 1)т]} = (к +1) Т2Б* + (1 + 22 + ••• к2 ) Т4Би + 2 к(1 + 2 — к) Т2т2 =

= (к +1) Т2Б* + к (к + 1)(2к +1 Т4Ои + к (к +1) Т2 т2. (10)

6

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

Нетрудно убедиться в корректности выполненных преобразований путем сопоставления непрерывной (6) и дискретной (10) форм записи уравнения ошибки.

Учитывая, что связь между дискретными и непрерывными процессами устанавливается предельным переходом, когда выборки берутся достаточно часто к ^ , Т ^ 0, кТ ^ t можно записать:

Ь1М к (к + 1)(2к +1) т 4В = = ¥/ . к-^мо 6 «33′

Ы .М. к (к + 1)т2 т2 = t2 т2′ Ы .М. (к + 1)т2 О = = Ж. (11)

Объединяя уравнения (11) получаем

Ы М. М К [(к + 1)Т ]}=^ + т2t2 + ¥ t.

При выводе соотношений учитывалась связь между интенсивностями непрерывных и дискретных шумовых процессов типа ШЫ через период дискретизации, а именно

¥ = ТБ ¥ = ТБ

* ^ и и ■

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

Заключение

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

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

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

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

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

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

Список литературы

1. Лукьянов Д., Филатов Ю., Голяев Ю., Курятов В., Виноградов В., Шрайбер К.-У., Перлмуттер М. 50 лет лазерному гироскопу // Фотоника. 2014. № 1(43). С. 42-61.

2. Матвеев В.В., Распопов В.Я. Основы построения бесплатформенных инерциальных навигационных систем: учеб. пособие / Под общ. ред. В.Я. Распопова. СПб.: ОАО «Концерн «ЦНИИ Электроприбор», 2009. 278 с.

3. Бендат Дж., Пирсол А. Прикладной анализ случайных данных: пер. с англ. / Под. ред. И.Н.Коваленко. М.: Мир, 1989. 540 с. [Bendat J.S., Piersol A.G. Random data: analysis and measurement procedures. 2nd ed. N.Y.: Wiley, 1986. 566 p.].

4. Salychev O.S. Applied inertial navigation: Problems and solutions. Moscow: BMSTU Press, 2004. 302 p.

5. Salychev O.S. MEMS-based inertial navigation: Expectations and reality. Moscow: BMSTU Press, 2012. 208 p.

6. Jazwinski A.H. Stochastic processes and filtering theory. N.Y.: Academic Press, 1970. 376 p.

7. Басараб М.А., Кравченко В.Ф., Матвеев В.А. Методы моделирования и цифровая обработка сигналов в гироскопии. М.: Физматлит, 2008. 248 с.

8. Применения лазеров: пер. с англ. / Под ред. В.П. Тычинского. М.: Мир, 1974. 445 с. [Laser applications / M. Ross, F. Aronowitz, B.J. Thompson, J.C. Owens a.o. N.Y.; L.: Academic Press, 1971].

9. Aranda J., De la Cruz J.M., Ruiperez H., Dormido S. Calibration and stochastic modeling of a laser-gyro for laboratory testing // Mathematical and Computer Modelling. 1990. Vol. 14. Pp. 231-236. DOI: 10.1016/0895-7177(90)90181-L

10. Mark J., Brown A., Matthews T. Quantization reduction for evaluating laser gyro performance using a moving average filter // Intern. conf. on acoustics, speech and signal processing: ICASSP 1984 (San Diego, CA, USA, March 19-21, 1984): Proc. 1984.

Pp. 614-617. DOI: 10.1109/ICASSP.1984.1172629

11. Coccoli J.D., Helfant S. RLG evaluation: Complementary modeling and testing // Proc. of the NAECON. 1979. Vol. 1. Pp. 14-21.

12. Chow W.W., Gea-Banacloche J., Pedrotti L.M., Sanders V.E., Schleich W., Scully M.O. The ring laser gyro // Reviews of Modern Physics. 1985. Vol. 57. Iss.1. Pp. 61-104. DOI: 10.1103/RevModPhys.57.61

13. Giardina C.R. Effect on navigation performance by noise present in RLG // Guidance and control conf. (San-Diego, CA, USA, August 16-18, 1976): Proc. N.Y.: AIAA, 1976.

Pp. 169-179.

14. Baum R.A., Slabinski R.J., Sturner B.A. A low-noise, high-bandwidth precision gyro for space pointing // Guidance and control 1980: Annual Rocky Mountain conf. (Keystone, Colo., USA, February 17-21, 1980): Proc. San Diego, 1980. Pp. 501-550.

15. Hummel S.G., Tesch E.L. Laser gyros in precision spacecraft attitude determination systems // Guidance and control conf. (Boulder, Colo., USA, August 6-8, 1979): Collection of technical papers. N.Y.: AIAA, 1979. Pp. 143-154. DOI: 10.2514/6.1979-1716

16. Oravetz A.S., Sandberg H.J. Stationary and nonstationary characteristics of gyro drift rate. AIAA J. 1970. Vol. 8. No. 10. Pp. 1766-1772. DOI: 10.2514/3.5988

17. Sudhakar M. Paniit, Wwibang Zhang. Modeling random gyro drift rate by data dependent systems // IEEE Trans. on Aerospace and Electronics Systems. 1986. Vol. 22. No. 4.

Pp. 455-460. DOI: 10.1109/TAES.1986.310781

18. Bilger H.R. Low frequency noise in ring laser gyros // Proc. of the Society of Photo-Optical Instrumentation Engineers (SPIE). 1984. Vol. 487. DOI: 10.1117/12.943248

19. Lawrence C.Ng, Darryll J. Pines. Characterization of ring laser gyro performance using the Allan variance method // J. of Guidance, Control, and Dynamics. 1997. Vol. 20. No. 1. Pp. 211-214. DOI: 10.2514/2.4026

20. Алалуев Р.В., Иванов Ю.В., Матвеев В.В., Орлов В.А., Распопов В.Я. Измерительный модуль микросистемной бесплатформенной инерциальной навигационной системы // Нано- и микросистемная техника. 2007. № 9. С. 61-64.

21. Букингем М. Шумы в электронных приборах и системах: пер. с англ. М.: Мир, 1986. 398 с. [Buckingham M.J. Noise in electronic devices and systems. Chichester; N.Y.: Halsted press, 1983. 372 p.].

Science ¿Education

of the Baurnan MSTU

Science and Education of the Bauman MSTU, 2017, no. 07, pp. 189-205.

DOI: 10.7463/0717.0001271

Received: 13.06.2017

Revised: 27.06.2017

© Bauman Moscow State Technical Unversity

Mathematical Model of the Laser Gyro Errors

V.V. Novikov1, V.N. Enin1’*, «enm^hmstuju

P.V. Novikov2

1Bauman Moscow State Technical University, Moscow, Russia 2Moscow Polytechnic University, Moscow, Russia

Keywords: laser gyro, random drift, systematic drift, natural and technical fluctuations, mathematical

error model, white noise, random walk

The paper presents the analysed and systemised results of the experimental study of laser gyro (LG) errors. Determines a structure of the resulting LG error, as a linear combination of the random processes, characterizing natural and technical fluctuations of difference frequency of the counter-propagating waves, with a random constant zero shift available in the sensor readings. Formulates the requirements for the structure and form of the analytic description of the error model. Shows a generalized model of the LG fluctuation processes, on the basis of which a mathematical model of LG errors was developed as an inertial sensor.

The model is represented by a system of the stochastic differential equations and functional relationships to characterize a resulting error of the sensor. The paper provides a correlation analysis of the model equations and final equations obtained for the mean-square values of the particular components, which allow us to identify the resulting error parameters. The model parameters are presented through the values of the power spectral density of the particular components. The discrete form of the model is considered, the convergence of continuous and difference equations is shown in fulfilling conditions of the limiting transition. Further research activities are defined.

References

1. Luk’ianov D.P., Filatov Yu.V., Goliaev Yu.D., Kuriatov V.N., Vinogradov V.I., Schreiber K.-U., Perlmutter M. Laser gyroscope is 50 years old. Fotonika [Photonics], 2014, no. 1(43), pp. 42-61 (in Russian).

2. Matveev V.V., Raspopov V.Ia. Osnovy postroeniia besplatformennykh inertsial’nykh navigatsionnykh system [Basis of strapdown inertial navigation systems buildings]: a text-

book / Ed. by V.Ia. Raspopov. S.-Petersburg: OAO «Konzern «TSNII Elektropribor» Publ., 2009. 278 p. (in Russian).

3. Bendat J.S., Piersol A.G. Random data: analysis and measurement procedures. 2nd ed. N.Y.: Wiley, 1986. (Russ. ed.: Bendat J.S., Piersol A.G. Prikladnoj analiz sluchajnykh dannykh. Moscow: Mir Publ., 1989. 540 p.).

4. Salychev O.S. Applied inertial navigation: Problems and solutions. Moscow: BMSTU Press, 2004. 302 p.

5. Salychev O.S. MEMS-based inertial navigation: Expectations and reality. Moscow: BMSTU Press, 2012. 208 p.

6. Jazwinski A.H. Stochastic processes and filtering theory. N.Y.: Academic Press, 1970. 376 p.

7. Basarab M.A., Kravchenko V.F., Matveev V.A. Metody modelirovaniia i tsifrovaia obrabotka signalov v giroskopii [Methods for modeling and digital signal processing in gyroscopy]. Moscow: Fizmatlit Publ., 2008. 248 p. (in Russian).

8. Laser applications / M. Ross, F. Aronowitz, B.J. Thompson, J.C. Owens a.o. N.Y.; L.: Academic Press, 1971. (Russ. ed.: Primeneniia lazerov. Moscow: Mir Publ., 1974. 445 p.).

9. Aranda J., De la Cruz J.M., Ruiperez H., Dormido S. Calibration and stochastic modeling of a laser-gyro for laboratory testing. Mathematical and Computer Modelling, 1990, vol. 14, pp. 231-236. DOI: 10.1016/0895-7177(90)90181-L

10. Mark J., Brown A., Matthews T. Quantization reduction for evaluating laser gyro performance using a moving average filter. Intern. conf. on acoustics, speech and signal processing: ICASSP 1984 (San Diego, CA, USA, March 19-21, 1984): Proc. 1984. Pp. 614-617. DOI: 10.1109/ICASSP.1984.1172629

11. Coccoli J.D., Helfant S. RLG evaluation: Complementary modeling and testing. Proc. of the NAECON, 1979, vol. 1, pp. 14-21.

12. Chow W.W., Gea-Banacloche J., Pedrotti L.M., Sanders V.E., Schleich W., Scully M.O. The ring laser gyro. Reviews of Modern Physics, 1985, vol. 57, iss.1, pp. 61-104.

DOI: 10.1103/RevModPhys.57.61

13. Giardina C.R. Effect on navigation performance by noise present in RLG. Guidance and control conf (San-Diego, CA, USA, August 16-18, 1976): Proc. N.Y.: AIAA, 1976.

Pp. 169-179.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

14. Baum R.A., Slabinski R.J., Sturner B.A. A low-noise, high-bandwidth precision gyro for space pointing. Guidance and control 1980: Annual Rocky Mountain conf. (Keystone, Colo., USA, February 17-21, 1980): Proc. San Diego, 1980. Pp. 501-550.

15. Hummel S.G., Tesch E.L. Laser gyros in precision spacecraft attitude determination systems. Guidance and control conf. (Boulder, Colo., USA, August 6-8, 1979): Collection of technical papers. N.Y.: AIAA, 1979. Pp. 143-154. DOI: 10.2514/6.1979-1716

16. Oravetz A.S., Sandberg H.J. Stationary and nonstationary characteristics of gyro drift rate. AIAA J.., 1970, vol. 8, no. 10, pp. 1766-1772. DOI: 10.2514/3.5988

17. Sudhakar M. Paniit, Wwibang Zhang. Modeling random gyro drift rate by data dependent systems. IEEE Trans. on Aerospace and Electronics Systems, 1986, vol. 22, no. 4,

pp. 455-460. DOI: 10.1109/TAES.1986.310781

18. Bilger H.R. Low frequency noise in ring laser gyros. Proc. of the Society of Photo-Optical Instrumentation Engineers (SPIE), 1984, vol. 487. DOI: 10.1117/12.943248

19. Lawrence C. Ng, Darryll J. Pines. Characterization of ring laser gyro performance using the Allan variance method. J. of Guidance, Control, and Dynamics, 1997, vol. 20, no. 1,

pp. 211-214. DOI: 10.2514/2.4026

20. Alaluev R.V., Ivanov Yu.V., Matveev V.V., Orlov V.A., Raspopov V.Ia. Measuring module of micro system platformless inertial navigation system. Nano — i mikrosistemnaia tekhnika [Nano- and Microsystems Technology], 2007, no. 9, pp. 61-64 (in Russian).

21. Buckingham M.J. Noise in electronic devices and systems. Chichester; N.Y.: Halsted press, 1983. (Russ. ed.: Buckingham M.J. Shumy v electronnykh priborakh i systemakh. Moscow: Mir Publ., 1986. 398 p.).

Каталог статей

    • Датчик IMU
      • Кинематика вращения
      • Модель измерения и модель движения
        • МЭМС акселерометр
        • Гироскоп
      • Модель ошибки IMU
        • Детерминированная ошибка
        • Метод калибровки детерминированной ошибки — шестисторонний метод
        • Температурная корреляция
        • Случайная ошибка
        • Метод калибровки случайных ошибок — метод дисперсии Аллана
      • математическая модель
        • Акселерометр
        • Гироскоп
        • Модель в VIO
      • Модель движения
      • Дискретный метод интегрирования
        • Метод Эйлера
        • Метод средней точки (средняя точка)
        • Рунге Кута четвертого порядка (RK4)
      • Общий процесс калибровки

Датчик IMU

Кинематика вращения

[ , , (img-dUA3Tb8z-1584641636300)(img/Untitled/1560661643256.png)]

Линейная скорость вращения составляет:

r

˙

=

(

a

θ

˙

sin

θ

,

a

θ

˙

cos

θ

,

0

)

=

[

0

θ

˙

0

θ

˙

0

0

0

0

0

]

[

a

cos

θ

a

sin

θ

h

]

=

ω

×

r

(1)

begin{aligned} dot{mathbf{r}} &=(-a dot{theta} sin theta, a dot{theta} cos theta, 0)^{top} \ &=left[begin{array}{ccc}{0} & {-dot{theta}} & {0} \ {dot{theta}} & {0} & {0} \ {0} & {0} & {0}end{array}right]left[begin{array}{c}{a cos theta} \ {a sin theta} \ {h}end{array}right] \ &=boldsymbol{omega} times mathbf{r} end{aligned} tag{1}

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

r

˙

=

ω

r

sin

ϕ

=

a

θ

˙

(2)

|dot{mathbf{r}}|=|boldsymbol{omega}||mathbf{r}| sin phi=a|dot{theta}| tag{2}

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

r

˙

I

=

R

I

B

r

˙

B

+

R

˙

I

B

r

B

=

R

I

B

r

˙

B

+

[

R

I

B

ω

b

]

×

r

=

R

I

B

v

b

+

ω

×

r

(3)

begin{aligned} dot{mathbf{r}}_{I} &=mathbf{R}_{I B} dot{mathbf{r}}_{B}+dot{mathbf{R}}_{I B} mathbf{r}_{B} \ &=mathbf{R}_{I B} dot{mathbf{r}}_{B}+left[mathbf{R}_{I B} boldsymbol{omega}_{b}right]_{ times} mathbf{r} \ &=mathbf{R}_{I B} mathbf{v}_{b}+boldsymbol{omega} times mathbf{r} end{aligned} tag{3}

v

I

R

I

B

v

b

+

ω

×

r

R

I

B

v

b

v

I

ω

×

r

(4)

mathbf{v}_{I} equiv mathbf{R}_{I B} mathbf{v}_{b}+omega times mathbf{r} Leftrightarrow mathbf{R}_{I B} mathbf{v}_{b} equiv mathbf{v}_{I}-boldsymbol{omega} times mathbf{r} tag{4}

Производная скорости в формуле (3) для получения ускорения:

r

¨

I

=

R

I

B

v

˙

B

+

R

˙

I

B

v

B

+

ω

×

r

˙

+

[

R

˙

I

B

ω

b

+

R

I

B

ω

˙

b

]

×

r

=

R

I

B

a

B

+

2

ω

×

v

+

ω

×

(

ω

×

r

)

+

ω

˙

×

r

(5)

begin{aligned} ddot{mathbf{r}}_{I} &=mathbf{R}_{I B} dot{mathbf{v}}_{B}+dot{mathbf{R}}_{I B} mathbf{v}_{B}+boldsymbol{omega} times dot{mathbf{r}}+left[dot{mathbf{R}}_{I B} boldsymbol{omega}_{b}+mathbf{R}_{I B} dot{boldsymbol{omega}}_{b}right]_{ times} mathbf{r} \ &=mathbf{R}_{I B} mathbf{a}_{B}+2 boldsymbol{omega} times mathbf{v}+boldsymbol{omega} times(boldsymbol{omega} times mathbf{r})+dot{boldsymbol{omega}} times mathbf{r} end{aligned} tag{5}

Преобразуйте формулу (5), чтобы получить:

a

=

a

I

2

ω

×

v

Филиал

Ши

сила

ω

˙

×

r

Европа

Тянуть

сила

ω

×

(

ω

×

r

)

из

сердце

сила

(6)

mathbf {a} = mathbf {a} _ {I} — underbrace {2 boldsymbol { omega} times mathbf {v}} _ { text Кох Force} — underbrace { dot { boldsymbol { omega}} times mathbf {r}} _ { text Euler Force} — underbrace { boldsymbol { omega} times ( boldsymbol { omega} times mathbf {r})} _ { text центробежная сила} tag {6}

где:

v

=

R

I

B

v

b

mathbf{v}=mathbf{R}_{I B} mathbf{v}_{b}

a

=

R

I

B

a

b

mathbf{a}=mathbf{R}_{I B} mathbf{a}_{b}

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

Модель измерения и модель движения

МЭМС акселерометр

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

  • Масса акселерометра не движется с высокой скоростью и не подвержена влиянию силы Кориолиса.

Гироскоп

  • Вибрирующий гироскоп использует силу Кориолиса для измерения вращения, ось главного движения + ось чувствительности, а сила Кориолиса применяется в направлении, перпендикулярном направлению движения (направлению оси чувствительности).

[ , , (img-sVL2ow4K-1584641636301)(img/VIO-DB/1560669359061.png)]

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

Модель ошибки IMU

Подразделяется на: детерминированная ошибка, случайная ошибка. Первое можно откалибровать заранее, а второе моделируется как распределение Гаусса.

Детерминированная ошибка

  • Смещение
  • Масштаб
  • Отклонение оси (неортогональная ошибка)
  • Коэффициент смещения / масштабирования между запусками означает, что значение меняется каждый раз при включении.
  • В рабочем режиме (стабильность) коэффициент смещения / масштабирования также изменяется во время работы.
  • Температурно-зависимый коэффициент смещения / масштабирования — это величина изменения, на которую влияет температура

Коэффициент масштабирования и ошибка асимметрии:

[

l

a

x

l

a

y

l

a

z

]

=

[

s

x

x

m

x

y

m

x

z

m

y

x

s

y

y

m

y

z

m

z

x

m

z

y

s

z

z

]

[

a

x

a

y

a

z

]

(7)

left[begin{array}{l}{l_{a x}} \ {l_{a y}} \ {l_{a z}}end{array}right]=left[begin{array}{lll}{s_{x x}} & {m_{x y}} & {m_{x z}} \ {m_{y x}} & {s_{y y}} & {m_{y z}} \ {m_{z x}} & {m_{z y}} & {s_{z z}}end{array}right]left[begin{array}{c}{a_{x}} \ {a_{y}} \ {a_{z}}end{array}right] tag{7}

Метод калибровки детерминированной ошибки — шестисторонний метод

Для акселерометра;

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

b

=

l

f

u

p

+

l

f

d

o

w

n

2

S

=

l

f

u

p

l

f

d

o

w

n

2

g

(8)

begin{aligned} b &=frac{l_{f}^{u p}+l_{f}^{d o w n}}{2} \ S &=frac{l_{f}^{u p}-l_{f}^{d o w n}}{2 cdot g} end{aligned} tag{8}

При рассмотрении межосевой ошибки соотношение между фактическим значением и измеренным значением

[

l

a

x

l

a

y

l

a

z

]

=

[

s

x

x

m

x

y

m

x

z

m

y

x

s

y

y

m

y

z

m

z

x

m

z

y

s

z

z

]

[

a

x

a

y

a

z

]

+

[

b

a

x

b

a

y

b

a

z

]

(9)

left[begin{array}{l}{l_{a x}} \ {l_{a y}} \ {l_{a z}}end{array}right]=left[begin{array}{lll}{s_{x x}} & {m_{x y}} & {m_{x z}} \ {m_{y x}} & {s_{y y}} & {m_{y z}} \ {m_{z x}} & {m_{z y}} & {s_{z z}}end{array}right]left[begin{array}{c}{a_{x}} \ {a_{y}} \ {a_{z}}end{array}right]+left[begin{array}{c}{b_{a x}} \ {b_{a y}} \ {b_{a z}}end{array}right] tag{9}

становится однородными координатами, чтобы получить:

[

l

a

x

l

a

y

l

a

z

]

=

[

s

x

x

m

x

y

m

x

z

b

a

x

m

y

x

s

y

y

m

y

z

b

a

y

m

z

x

m

z

y

s

z

z

b

a

z

]

[

a

x

a

y

a

z

1

]

(10)

left[begin{array}{l}{l_{a x}} \ {l_{a y}} \ {l_{a z}}end{array}right]=left[begin{array}{lll}{s_{x x}} & {m_{x y}} & {m_{x z}} &{b_{a x}} \ {m_{y x}} & {s_{y y}} & {m_{y z}} &{b_{a y}} \ {m_{z x}} & {m_{z y}} & {s_{z z}} &{b_{a z}}end{array}right]left[begin{array}{c}{a_{x}} \ {a_{y}} \ {a_{z}} \ 1end{array}right] tag{10}

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

Для гироскопа;

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

Температурная корреляция

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

  • Метод замачивания: контролируйте различные значения постоянной температуры и измеряйте данные датчика.
  • метод линейного изменения: записывайте данные датчика во время линейного нагрева и охлаждения в течение определенного периода времени.

Случайная ошибка

Гауссовский белый шумМоделируется как среднее значение 0, дисперсия

σ

sigma

Независимый гауссовский процесс в каждый момент, как правило, является внешним шумом, вызванным аналого-цифровым преобразователем.

E

[

n

(

t

)

]

0

E

[

n

(

t

1

)

n

(

t

2

)

]

=

σ

2

δ

(

t

1

t

2

)

(11)

begin{array}{l}{E[n(t)] equiv 0} \ {Eleft[nleft(t_{1}right) nleft(t_{2}right)right]=sigma^{2} deltaleft(t_{1}-t_{2}right)}end{array} tag{11}

где

δ

(

)

delta()

Представляет функцию Дирака.

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

n

d

[

k

]

n

(

t

0

+

Δ

t

)

1

Δ

t

t

0

t

0

+

Δ

t

n

(

τ

)

d

t

E

(

n

d

[

k

]

2

)

=

E

(

1

Δ

t

2

t

0

t

0

+

Δ

t

t

0

t

0

+

Δ

t

n

(

τ

)

n

(

t

)

d

τ

d

t

)

=

E

(

σ

2

Δ

t

2

t

0

t

0

+

Δ

t

t

0

t

0

+

Δ

t

δ

(

t

τ

)

d

τ

d

t

)

=

E

(

σ

2

Δ

t

)

(12)

n_{d}[k] triangleq nleft(t_{0}+Delta tright) simeq frac{1}{Delta t} int_{t_{0}}^{t_{0}+Delta t} n(tau) d t \begin{aligned} Eleft(n_{d}[k]^{2}right) &=Eleft(frac{1}{Delta t^{2}} int_{t_{0}}^{t_{0}+Delta t} int_{t_{0}}^{t_{0}+Delta t} n(tau) n(t) d tau d tright) \ &=Eleft(frac{sigma^{2}}{Delta t^{2}} int_{t_{0}}^{t_{0}+Delta t} int_{t_{0}}^{t_{0}+Delta t} delta(t-tau) d tau d tright) \ &=Eleft(frac{sigma^{2}}{Delta t}right) end{aligned} tag{12}

Итак, дискретный гауссовский белый шум:

n

d

[

k

]

=

σ

d

w

[

k

]

(13)

n_{d}[k]=sigma_{d} w[k] tag{13}

где:

w

[

k

]

N

(

0

,

1

)

w[k] sim mathcal{N}(0,1)

σ

d

=

σ

1

Δ

t

sigma_{d}=sigma frac{1}{sqrt{Delta t}}

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

b

˙

(

t

)

=

n

(

t

)

=

σ

b

w

(

t

)

(14)

dot{b}(t)=n(t)=sigma_{b} w(t) tag{14}

То же преобразование между дискретным и непрерывным:

b

d

[

k

]

b

(

t

0

)

+

t

0

t

0

+

Δ

t

n

(

t

)

d

t

E

(

(

b

d

[

k

]

b

d

[

k

1

]

2

)

=

E

(

t

0

t

0

+

Δ

t

t

0

t

0

+

Δ

t

n

(

t

)

n

(

τ

)

d

τ

d

t

)

=

E

(

σ

b

2

t

0

t

0

+

Δ

t

t

0

t

0

+

Δ

t

δ

(

t

τ

)

d

τ

d

t

)

=

E

(

σ

b

2

Δ

t

)

(15)

begin{aligned} b_{d}[k] triangleq & bleft(t_{0}right)+int_{t_{0}}^{t_{0}+Delta t} n(t) d t \ Eleft(left(b_{d}[k]-b_{d}[k-1]^{2}right)right.&=Eleft(int_{t_{0}}^{t_{0}+Delta t} int_{t_{0}}^{t_{0}+Delta t} n(t) n(tau) d tau d tright) \ &=Eleft(sigma_{b}^{2} int_{t_{0}}^{t_{0}+Delta t} int_{t_{0}}^{t_{0}+Delta t} delta(t-tau) d tau d tright) \ &=Eleft(sigma_{b}^{2} Delta tright) end{aligned} tag{15}

означает:

b

d

[

k

]

=

b

d

[

k

1

]

+

σ

b

d

w

[

k

]

(16)

b_{d}[k]=b_{d}[k-1]+sigma_{b d} w[k] tag{16}

где:

w

[

k

]

N

(

0

,

1

)

w[k] sim mathcal{N}(0,1)

σ

b

d

=

σ

b

Δ

t

sigma_{b d}=sigma_{b} sqrt{Delta t}

Метод калибровки случайных ошибок — метод дисперсии Аллана

Конкретный процесс выглядит следующим образом:

  1. Положите гироскоп на время

    T

    T

    , Единый период выборки

    τ

    0

    tau_0

    Из

    N

    N

    Значение выборки группы.

  2. Рассчитайте выходной угол одного образца

    θ

    theta

    И средний коэффициент

    m

    m

    m

    m

    Постарайтесь добиться ровности.

θ

(

n

)

=

τ

0

Ω

(

n

)

d

t

m

=

Rand

(

1

,

N

1

2

)

(17)

begin{aligned} theta(n) &=int^{tau_{0}} Omega(n) d t \ m &=operatorname{Rand}left(1, frac{N-1}{2}right) end{aligned} tag{17}

  1. Рассчитайте дисперсию Аллана, когда

    m

    m

    При выборе разных значений будут разные значения дисперсии Аллана.

θ

2

(

τ

)

=

1

2

τ

2

(

N

2

m

)

k

=

1

N

2

m

(

θ

K

+

2

m

2

θ

K

+

m

+

θ

K

)

2

(18)

theta^{2}(tau)=frac{1}{2tau^{2}(N-2m)} sum^{N-2m}_{k=1}(theta_{K+2m}-2theta_{K+m}+theta_{K})^2 tag{18}

среди них:

τ

=

m

τ

0

tau=m tau_{0}

  1. Обычно квадратный корень из дисперсии Аллана используется при построении кривой дисперсии Аллана, и берется корень из формулы (18).

[ , , (img-O62HubuW-1584641636303)(img/VIO-DB/1560696133790.png)]

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

[ , , (img-a0wDmW7b-1584641636304)(img/VIO-DB/1560696833421.png)]

[ , , (img-unMFIj7e-1584641636307)(img/VIO-DB/1560696853535.png)]

[ , , (img-pfNnnVkz-1584641636309)(img/VIO-DB/1560697012276.png)]

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

τ

tau

Сила.

Учитывая, что эти источники шума не зависят друг от друга, дисперсия Аллана представляет собой сумму квадратов различных типов ошибок:

σ

 total 

2

(

τ

)

=

σ

Q

2

(

τ

)

+

σ

N

2

(

τ

)

+

σ

B

2

(

τ

)

+

σ

K

2

(

τ

)

+

σ

R

2

(

τ

)

=

3

Q

2

τ

2

+

N

2

τ

+

2

B

2

π

ln

2

+

K

2

τ

3

+

R

2

τ

2

2

(19)

begin{aligned} sigma_{text { total }}^{2}(tau) &=sigma_{Q}^{2}(tau)+sigma_{N}^{2}(tau)+sigma_{B}^{2}(tau)+sigma_{K}^{2}(tau)+sigma_{R}^{2}(tau) \ &=frac{3 Q^{2}}{tau^{2}}+frac{N^{2}}{tau}+frac{2 B^{2}}{pi} ln 2+frac{K^{2} tau}{3}+frac{R^{2} tau^{2}}{2} end{aligned} tag{19}

Организуйте формулу (19):

σ

t

o

t

a

l

2

(

τ

)

=

i

=

2

2

C

i

τ

i

(20)

sigma_{mathrm{total}}^{2}(tau)=sum_{i=-2}^{2} C_{i} tau^{i} tag{20}

где:

i

=

(

2

,

1

,

0

,

1

,

2

)

i=(-2,-1,0,1,2)

Используйте метод наименьших квадратов, чтобы подобрать значения, рассчитанные по формуле (20) и формуле (18), чтобы получить

C

i

C_i

Value, а затем получить соответствующий коэффициент ошибки:

{

Q

=

C

2

3600

3

(

)

N

=

C

1

60

(

)

/

h

B

=

C

0

0.664

(

)

/

h

K

=

60

3

C

1

[

(

)

h

1

]

/

h

R

=

3600

2

C

2

[

(

)

h

1

]

/

h

(21)

left{begin{array}{l}{Q=frac{sqrt{C_{-2}}}{3600 sqrt{3}}left(^{circ}right)} \ {N=frac{sqrt{C_{-1}}}{60}left(^{circ}right) / sqrt{h}} \ {B=frac{sqrt{C_{0}}}{0.664}left(^{circ}right) / h} \ {K=60 sqrt{3 C_{1}}left[left(^{circ}right) cdot h^{-1}right] / sqrt{h}} \ {R=3600 sqrt{2 C_{2}}left[left(^{circ}right) cdot h^{-1}right] / h}end{array}right. tag{21}

Данные гироскопа преобразуются в град / ч,

τ

tau

Блок Still sec, акселерометр по-прежнему м / с2

imu_utilsСреди них — использование цереры для решения наименьших квадратов для получения соответствующих коэффициентов, процесс перехода от формулы (18) к формуле (21), обратите внимание на преобразование единиц.

Я не понимаю, что эти коэффициенты ошибок применяются к каждому сегменту отдельно? Тест, найденный в коде, должен сделать

τ

tau

Это полученное значение. Нарисуется общая кривая дисперсии Аллана.Почему ее можно подогнать для представления только части шума? Как формула (19) получает только шум?

Найдите объяснение: ошибка акселерометра в диапазоне среднего времени порядка 1 с — это в основном случайное блуждание по скорости, а диапазон выше 100 с — в основном нестабильность смещения. Для более длительного среднего времени он ограничивается сбором данных Длина и количество независимых подмножеств, используемых в дисперсионном анализе Аллана, ограничены, поэтому его достоверность очень низкая.

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

математическая модель

Акселерометр

a

m

B

=

S

a

T

a

R

B

G

(

a

G

g

G

)

+

n

a

+

b

a

(22)

mathbf{a}_{m}^{B}=mathbf{S}_{a}mathbf{T}_{a} mathbf{R}_{B G}left(mathbf{a}^{G}-mathbf{g}^{G}right)+mathbf{n}_{a}+mathbf{b}_{a} tag{22}

Гироскоп

ω

m

B

=

S

g

T

g

ω

B

+

s

g

a

a

B

+

n

g

+

b

g

(23)

boldsymbol{omega}_{m}^{B}=mathbf{S}_{g} mathbf{T}_{g}boldsymbol{omega}^{B}+mathbf{s}_{g a} mathbf{a}^{B}+mathbf{n}_{g}+mathbf{b}_{g} tag{23}

S

mathbf S

Масштабный коэффициент,

T

mathbf T

Асимметричная ошибка,

n

mathbf n

Гауссовский белый шум,

b

mathbf b

Это случайное блуждание с нулевым смещением,

s

g

a

mathbf{s}_{g a}

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

Модель в VIO

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

ω

~

b

=

ω

b

+

b

g

+

n

g

a

~

b

=

q

b

w

(

a

w

+

g

w

)

+

b

a

+

n

a

(24)

begin{aligned} tilde{omega}^{b} &=omega^{b}+mathbf{b}^{g}+mathbf{n}^{g} \ tilde{mathbf{a}}^{b} &=mathbf{q}_{b w}left(mathbf{a}^{w}+mathbf{g}^{w}right)+mathbf{b}^{a}+mathbf{n}^{a} end{aligned} tag{24}

здесь **

g

g

Знак плюс минус ** обратите на это внимание, а

g

g

Определенные положительные и отрицательные значения связаны с положительным и отрицательным ускорением, выдаваемым акселерометром.

Модель движения

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

Дискретный метод интегрирования

p

w

b

k

+

1

=

p

w

b

k

+

v

k

w

Δ

t

+

1

2

a

Δ

t

2

v

k

+

1

w

=

v

k

w

+

a

Δ

t

q

w

b

k

+

1

=

q

w

b

k

[

1

1

2

ω

δ

t

]

(25)

begin{aligned}mathbf{p}_{w b_{k+1}}&=mathbf{p}_{w b_{k}}+mathbf{v}_{k}^{w} Delta t+frac{1}{2} mathbf{a} Delta t^{2} \ mathbf{v}_{k+1}^{w}&=mathbf{v}_{k}^{w}+mathbf{a} Delta t \ mathbf{q}_{w b_{k+1}}&=mathbf{q}_{w b_{k}} otimesleft[begin{array}{c}{1} \ {frac{1}{2} omega delta t}end{array}right]end{aligned} tag{25}

Для нелинейных функций существуют следующие методы.

Метод Эйлера

X

n

+

1

=

X

n

+

K

Δ

t

K

=

f

(

t

n

,

X

n

)

(26)

mathbf X_{n+1} =mathbf X_{n} + mathrm KDelta t\ K=f'(t_{n}, mathbf X_n) tag{26}

Где в формуле (24)

a

=

q

w

b

k

(

a

b

k

b

k

a

)

g

w

ω

=

ω

b

k

b

k

g

begin{aligned} mathbf{a} &=mathbf{q}_{w b_{k}}left(mathbf{a}^{b_{k}}-mathbf{b}_{k}^{a}right)-mathbf{g}^{w} \ boldsymbol{omega} &=boldsymbol{omega}^{b_{k}}-mathbf{b}_{k}^{g} end{aligned}

Метод средней точки (средняя точка)

X

n

+

1

=

X

n

+

K

2

Δ

t

K

2

=

f

(

t

n

+

1

2

Δ

t

,

X

n

+

K

1

Δ

t

2

)

K

1

=

f

(

t

n

,

X

n

)

(27)

begin{aligned} mathbf X_{n+1} &= mathbf X_{n}+mathrm K_2Delta t \ mathrm K_2 &= f'(t_n+frac{1}{2}Delta t, mathbf X_n+mathrm K_1frac{Delta t}{2} ) \ mathrm K_1 &=f'(t_n, mathbf X_n) end{aligned} tag{27}

Где в формуле (24)

a

=

1

2

[

q

w

b

k

(

a

b

k

b

k

a

)

g

w

+

q

w

b

k

+

1

(

a

b

k

+

1

b

k

a

)

g

w

]

ω

=

1

2

[

(

ω

b

k

b

k

g

)

+

(

ω

b

k

+

1

b

k

g

)

]

begin{aligned} mathbf{a} &=frac{1}{2}left[mathbf{q}_{w b_{k}}left(mathbf{a}^{b_{k}}-mathbf{b}_{k}^{a}right)-mathbf{g}^{w}+mathbf{q}_{w b_{k+1}}left(mathbf{a}^{b_{k+1}}-mathbf{b}_{k}^{a}right)-mathbf{g}^{w}right] \ boldsymbol{omega} &=frac{1}{2}left[left(boldsymbol{omega}^{b_{k}}-mathbf{b}_{k}^{g}right)+left(boldsymbol{omega}^{b_{k+1}}-mathbf{b}_{k}^{g}right)right] end{aligned}

Рунге Кута четвертого порядка (RK4)

В дифференциальном уравнении метода Рунге-Кутта должна быть исходная величина, а именно

y

˙

=

f

(

t

,

y

)

dot y=f(t,y)

Форма применяется.

X

n

+

1

=

X

n

+

Δ

t

6

(

K

1

+

2

K

2

+

2

K

3

+

K

4

)

K

1

=

f

(

t

n

,

X

n

)

K

2

=

f

(

t

n

+

Δ

t

2

,

X

n

+

Δ

t

2

K

1

)

K

3

=

f

(

t

n

+

Δ

t

2

,

X

n

+

Δ

t

2

K

2

)

K

4

=

f

(

t

n

+

Δ

t

,

X

n

+

Δ

t

K

3

)

(28)

begin{aligned} mathbf X_{n+1} &= mathbf X_{n}+frac{Delta t}{6}(mathrm{K_1+2K_2+2K_3+K_4}) \ mathrm K_1 &=f'(t_n, mathbf X_n) \ mathrm K_2 &=f'(t_n+frac{Delta t}{2}, mathbf X_n +frac{Delta t}{2} mathrm K_1) \ mathrm K_3 &= f'(t_n+frac{Delta t}{2}, mathbf X_n +frac{Delta t}{2} mathrm K_2) \ mathrm K_4 &= f'(t_n+Delta t, mathbf X_n +{Delta t} mathrm K_3) \ end{aligned} tag{28}

где в формуле (24)

a

=

1

2

[

q

w

b

k

(

a

b

k

b

k

a

)

g

w

+

q

w

b

k

+

1

(

a

b

k

+

1

b

k

a

)

g

w

]

begin{aligned} mathbf{a} &=frac{1}{2}left[mathbf{q}_{w b_{k}}left(mathbf{a}^{b_{k}}-mathbf{b}_{k}^{a}right)-mathbf{g}^{w}+mathbf{q}_{w b_{k+1}}left(mathbf{a}^{b_{k+1}}-mathbf{b}_{k}^{a}right)-mathbf{g}^{w}right] end{aligned}

template <typename _T>
 inline void imu_tk::quatIntegrationStepRK4(
                 const Eigen::Matrix< _T, 4, 1> &quat,
                 const Eigen::Matrix< _T, 3, 1> &omega0,
                 const Eigen::Matrix< _T, 3, 1> &omega1,
                 const _T &dt, Eigen::Matrix< _T, 4, 1> &quat_res )
{
  Eigen::Matrix< _T, 3, 1> omega01 = _T(0.5)*( omega0 + omega1 );
  Eigen::Matrix< _T, 4, 1> k1, k2, k3, k4, tmp_q;
  Eigen::Matrix< _T, 4, 4> omega_skew;
  
  // First Runge-Kutta coefficient
  computeOmegaSkew( omega0, omega_skew );
  k1 = _T(0.5)*omega_skew*quat;
  // Second Runge-Kutta coefficient
  tmp_q = quat + _T(0.5)*dt*k1;
  computeOmegaSkew( omega01, omega_skew );
  k2 = _T(0.5)*omega_skew*tmp_q;
  // Third Runge-Kutta coefficient (same omega skew as second coeff.)
  tmp_q = quat + _T(0.5)*dt*k2;
  k3 = _T(0.5)*omega_skew*tmp_q;
  // Forth Runge-Kutta coefficient
  tmp_q = quat + dt*k3;
  computeOmegaSkew( omega1, omega_skew );
  k4 = _T(0.5)*omega_skew*tmp_q;
  _T mult1 = _T(1.0)/_T(6.0), mult2 = _T(1.0)/_T(3.0);
  quat_res = quat + dt*(mult1*k1 + mult2*k2 + mult2*k3 + mult1*k4);
  normalizeQuaternion(quat_res);
}

Общий процесс калибровки

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

[ , , (img-OB3gwTq0-1584641636310)(img/IMU /1560749248488.png)]

[ , , (img-437mRqJR-1584641636311)(img/IMU /1560749264561.png)]

[ , , (img-W0rtTOQ9-1584641636313)(img/IMU /1560749280077.png)]

[ , , (img-xY3VeGG8-1584641636316)(img/IMU /1560749299011.png)]

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

1.
Геометрические (карданные) погрешности.

Определение
положения летательного аппарата
относительно оси ротора, производится
посредством намерения углов поворота

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

2.
Скоростные кинематические пог­решности.

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

3.
Кинематические погрешности.

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

4.
Инструментальные погрешности.

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

15.3. Радиокомпас, принцип действия, погрешности

15.3.1. Основные определения

Курсовым
углом радиостанции (КУР)

называется угол, заключенный между
продольной осью самолета и действительным
направлением на радиостанцию. КУР
отсчитывается от продольной оси самолета
по ходу часовой стрелки до направления
на радиостанцию от 0 до 360° (рис.15.5).

Отсчетом
радиокомпаса (ОРК)

называется угол, заключенный между
продольной осью самолета и измеренным
направлением на радиостанцию. ОРК
отсчитывается от продольной оси самолета
до измеренного направления на радиостанцию
от 0 до 360°.

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

Пеленгом
радиостанции (ПР)

называется угол между меридианом начала
отсчета курса и направлением от самолета
на радиостанцию. ПР отсчитывается от
северного направления меридиана по
ходу часовой стрелки до направления на
радиостанцию от 0 до 360°.

В
зависимости от начала отсчета курса
самолета пеленги радиостанции могут
быть истинными ИПР и магнитными МПР.

Пеленги
радиостанции рассчитываются по формулам
:

МПР = МК + КУР;

ИПР=ИК+КУР

(15.7)

ИПР
= КК + ∆ К
+ ∆ М
+ КУР.

Рис.15.5.
Основные радионавигационные элементы

Пеленгом
самолета ПС

называется угол между меридианом
радиостанции и направлением от
радиостанции на самолет. ПС отсчитывается
от северного направления меридиана по
ходу часовой стрелки до направления на
самолет от 0 до 360°. В зависимости от
начала отсчета курса самолета пеленги
самолета могут быть истинными ИПС и
магнитными МПС. Пеленги самолета
рассчитываются по формулам:

МПС = МК + КУР ±
180°;

ИПС
= ИК- КУР ± 180°;
(15.8)

ИПС
= КК +∆ К
+ ∆ М
+ КУР ± 180°.

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

ИПС
= ИК + КУР ± 180″ + (± 8).

(15.9)

Для
карт видоизмененной поликонической
проекции поправка на угол схождения
меридианов будет равна:

(15.10)

где
-долгота
радиостанции,
долгота самолета,


средняя широта листа карты.

Соседние файлы в папке книга посл.редакция

  • #
  • #
  • #
  • #
  • #
  • #
  • #
  • #
  • #
  • #
  • #

Понравилась статья? Поделить с друзьями:
  • Модель ошибки наемных работников
  • Моделирование ошибок бинс
  • Могут ли быть ошибки при узи молочной железы
  • Модель коррекции ошибок эконометрика
  • Множество фильмов рождались на глазах актера ошибка