- •Глава 1 динамические системы в задачах обработки навигационной информации
- •1.1. Математические модели динамических систем и методы их описания
- •1.1.1. Определение и классификация динамических систем
- •1.1.2. Определение выхода линейных динамических систем с помощью фундаментальной матрицы и весовой функции
- •1.1.3. Передаточная функция стационарных линейных динамических систем
- •1.1.4. Частотная характеристика стационарных линейных динамических систем
- •1.1.5. Взаимосвязь между основными функциями, используемыми при описании линейных динамических систем
- •Основные функции, используемые при описании линейных систем
- •Взаимосвязь различных функций, используемых для описания поведения линейных стационарных систем
- •1.1.6. Определение модели линейной динамической системы в пространстве состояний по заданной передаточной функции
- •1.1.7. Основные свойства линейных динамических систем
- •Задачи к разделу
- •Контрольные вопросы
- •1.2. Линейные стационарные динамические системы
- •1.2.1. Интеграторы
- •1.2.2. Фильтры Баттерворта
- •1.2.3. Модель акселерометра
- •1.2.4. Модель микромеханического гироскопа
- •1.2.5. Простейшая модель ошибок построения вертикали в инерциальной системе
- •Вклад уходов гироскопов в ошибки определения нп в автономном режиме в течение 3 мин
- •Вклад ошибок акселерометров в ошибки определения нп в автономном режиме в течение 3 мин
- •Задачи к разделу
- •Контрольные вопросы
- •1.3. Дискретизация и моделирование линейных
- •1.3.1. Дискретизация непрерывных систем
- •1.3.2. Основные методы описания lti-объектов в Matlab
- •Создание lti-объектов
- •Извлечение информации о моделях
- •1.3.3. Особенности дискретизации стационарных систем
- •Методы дискретизации, используемые в функции c2d
- •Методы вычисления матричной экспоненты
- •Контрольные вопросы
- •1.4. Задание для моделирования
- •Пример выполнения задания в Matlab
- •Заключение к главе 1
1.2.4. Модель микромеханического гироскопа
В настоящее время при построении современных инерциальных навигационных систем широкое применение находят микромеханические гироскопы и акселерометры [19, 53, 54, 58]. Рассмотрим (рис. 1.2.4) модель для так называемого микромеханического гироскопа (ММГ) RR-типа (от англ. rotate-rotate) [59, с. 103].
Рис. 1.2.4. Микромеханический гироскоп RR-типа
1 ротор, 2 неподвижные гребенчатые структуры привода первичных колебаний вокруг оси Z; 3 упругие балки подвеса; 4 датчик угла вокруг оси X; 5 неподвижные гребенчатые структуры датчика угла вокруг оси Z
Принцип действия гироскопа заключается в следующем. При помощи специальной системы управления возбуждаются и поддерживаются на заданном уровне первичные угловые колебания ротора 1 относительно оси Z (главная ось или ось первичных колебаний). После появления подлежащей определению угловой скорости вокруг оси чувствительности Y возникает момент сил Кориолиса, величина которого определяется как , где значение кинетического момента ротора. Момент вызывает так называемые вторичные угловые колебания ротора гироскопа относительно оси X (ось прецессии или ось вторичных колебаний).
Измеряя величину вторичных колебаний ротора, можно определить угловую скорости . Кинетический момент создается за счет возбуждения и поддержания стабильными колебания ротора ММГ при помощи специальной системы управления. Физически момент управления создается датчиком момента или электростатическим приводом (5), имеющим форму гребенчатых структур. Часть гребенчатых структур сформирована в теле ротора, а другая закреплена на корпусе таким образом, чтобы зубцы подвижной и неподвижной частей имели взаимную область перекрытия площадей. В этом случае после подачи напряжения на неподвижную гребенку возникает электростатическая сила, стремящаяся затянуть зубцы подвижной гребенки в пазы неподвижной.
Для измерения угла поворота ротора относительно оси Z используются аналогичные гребенчатые структуры, образующие за счет электрической схемы включения датчик угла (2). При повороте ротора относительно оси Z площадь перекрытия зубцов изменяется, соответственно, изменяется и емкость датчика угла, регистрируемая в дальнейшем специальным преобразователем «емкостьнапряжение».
Измерения поворота ротора относительно оси вторичных колебаний Х также осуществляется при помощи датчика емкостного типа (4), неподвижные электроды которого расположены на корпусе, а в качестве общего электрода используется тело ротора. При повороте ротора относительно оси X изменяется зазор между электродами датчика, а следовательно и емкость, что так же как и в канале первичных колебаний, регистрируется с помощью специального преобразователя «емкость напряжение».
Можно показать, что упрощенная линеаризованная система уравнений движения ротора ММГ, не учитывающая перекрестные связи и вредные моменты, имеет следующий вид [25, 59, с.103]:
; (1.2.27)
, (1.2.28)
где , – углы отклонения ротора ММГ относительно осей первичных Z и вторичных X колебаний; , – показатели затухания колебаний ротора относительно осей Z и X; , – собственные частоты подвеса ротора ММГ относительно осей Z и X; – момент инерции ротора относительно оси X; , – управления, вводимые по осям первичных и вторичных колебаний.
Задачей управления по оси первичных колебаний является возбуждение и стабилизация колебаний ротора с заданной частотой и амплитудой . На практике частота возбуждения первичных колебаний часто совпадает с собственной частотой первичных колебаний ротора , что вызвано необходимостью использования резонансного усиления для достижения заданной амплитуды колебаний. После завершения переходных процессов движение ротора относительно оси Z принимает вид угловых колебаний:
; (1.2.29)
. (1.2.30)
Соответственно выражения для кинетического момента ротора и момента сил Кориолиса принимают вид гармонических и амплитудно-модулированных сигналов:
; (1.2.31)
, (1.2.32)
где – момент инерции ротора относительно оси Z.
С учетом (1.2.32) можно видеть, что уравнение (1.2.28) является нестационарным, поскольку входной сигнал угловой скорости модулируется по величине в соответствии с гармоническим законом с частотой первичных колебаний. Это также означает, что вынужденные колебания ротора по оси вторичных колебаний будут происходить на частоте первичных колебаний, и для выделения полезного сигнала из измеренного сигнала необходимо выполнять процедуру синхронного детектирования.
Также следует иметь в виду, что в общем случае значения моментов инерций , , показателей затухания , и собственных частот , по осям X и Z не равны между собой, что может оказывать существенное влияние на характеристики ММГ.
Вместе с тем нетрудно заметить, что левые части уравнений (1.2.27)(1.2.28), описывающее поведение угла поворота ротора вокруг оси первичных и вторичных колебаний, с точностью до обозначений совпадают с уравнением (1.2.20). Таким образом, собственное движение системы может быть исследовано с использованием полученных в подразделе 1.2.3 функций.
Для анализа движения ротора по оси вторичных колебаний в случае гармонического воздействия, модулированного по амплитуде необходимо привлечение теории анализа систем с амплитудной модуляцией.