Подходы к решению задачи стабилизации информационных характеристик микромеханических гироскопов роторного типа. Разработка программы моделирования устройства, позволяющей анализировать его информационные характеристики, инструментальные погрешности.
Аннотация к работе
ММГ и ММА входят в класс так называемых интегральных преобразователей, т.к. электромеханические узлы приборов этих типов формируются из неметаллических материалов (монокристаллический кремний, плавленый кварц, карбид кремния и др.) методами фотолитографии и изотропного или анизотропного травления вместе с элементами электроники возбуждения, датчиками съема и преобразования полезного сигнала, элементами формирования обратных связей и т.д. Единство материала и технологий микроэлектроники позволяет создавать приборы, в которых осуществляется интеграция физико-конструктивных элементов. С этой точки зрения все приборы, в которых осуществляется интеграция воспринимающего измеряемую величину (чувствительного) элемента, всех преобразующих элементов, а также элементов внутренней электроники можно назвать интегральными, т.к. для их изготовления используются технологии микроэлектроники. Микромеханические датчики первичной информации позволяют создавать малогабаритные инерциальные навигационные системы, интегрированные со спутниковыми навигационными системами (GPS, ГЛОНАСС), обеспечивающими необходимую точность определения параметров ориентации и координат подвижных объектов: дистанционно-пилотируемых летательных аппаратах, морских судов, наземных транспортных средств, робототехнических комплексов, медицинской технике, в оборонных областях.
Список литературы
По теме диссертации опубликовано 12 печатных работ, из них 2 в изданиях, рекомендованных в перечне ВАК.
Объем работы. Диссертация изложена на 119 страницах машинописного текста, включая 41 рисунок и 5 таблиц, состоит из введения, пяти глав, заключения, списка литературы, включающего 61 наименование.
На защиту выносятся следующие основные положения: 1. Обобщенная модель динамики микромеханических гироскопов роторного типа;
2. Алгоритм анализа основных информационных характеристик микромеханических гироскопов роторного типа, включающих масштабный коэффициент преобразования и полосу частот пропускания для различных схем обработки вторичных колебаний;
3. Оценка влияния взаимосвязей измерительных каналов на динамику и выходные сигналы двумерного ММГ;
4. Методика стабилизации информационных характеристик ММГ, введением контуров стабилизации амплитуды первичных колебаний и управления частотой расстройки первичных и вторичных колебаний;
5. Оценка квадратурных составляющих погрешностей двумерного ММГ, обусловленных влиянием контура первичных колебаний на измерительные каналы.
Содержание диссертации информационный гироскоп роторный программа
Во введении обосновывается актуальность, формулируются цели и задачи работы, определяются объекты исследования, констатируется научная новизна и практическая ценность полученных результатов.
В первой главе проводится обзор конструктивных схем микромеханических гироскопов с кремниевым чувствительным элементом, и формулируются основные задачи исследования. Конструктивная схема рассматриваемого двумерного ММГ роторного типа приведена на рис. 1. Основу конструкции микрогироскопа с двумя осями чувствительности составляет ротор, подвешенный на симметричном упругом подвесе, содержащем четыре наружных торсиона 2, четыре внутренних торсиона 3 и четыре полукольца 4. Наружные торсионы одним кольцом закреплены в специальных карманах в теле ротора, а другим в сегментах полуколец подвеса. Внутренние торсионы одним кольцом закреплены в сегментах полуколец подвеса, а другим на пятке 5, жестко связанной с корпусом прибора (подложкой). Для размещения электростатического двигателя системы возбуждения колебаний ротора в его теле образованы окна 6, по радиальным образующим которых размещена регулярная структура (зубцы 7 ротора гребенчатого двигателя).
Для упрощения на схеме показаны только единичные элементы гребенчатой структуры правого и левого двигателя в одном окне. Внутри каждого окна размещены два статора 8 электростатических двигателей, которые по радиальным образующим имеют такую же регулярную структуру из зубцов, как и ротор двигателя. Зубцы 9 статора находятся в пазах, образованных зубцами ротора. Емкостные датчики преобразователей первичных колебания ротора по оси возбуждения имеют аналогичную гребенчатую структуру. Емкостные датчики преобразователей вторичных колебаний ротора имеют плоскую структуру образованную свободными поверхностями на теле ротора и четырьмя изолированными обкладками, связанными с корпусом прибора.
Конструктивная особенность рассматриваемой системы упругого подвеса ротора с четырьмя S-образными упругими элементами, связывающими ротор с жестко закрепленной на корпусе центральной пяткой, состоит в практическом исключении деформаций растяжения упругих элементов. Кроме того, выбором отмеченных на схеме (рис. 1б) размеров конструкции гироскопа и использованием двух групп торсионов различной длины обеспечивается практически равенство частот собственных первичных и вторичных колебаний ротора. Для диапазона наружного радиуса ротора R=(1?2,5) мм эти частоты лежат в пределах ?р=(5?0,8) КГЦ.
Вторая глава посвящена составлению уравнений динамики ММГ роторного типа с двумя осями чувствительности, исследованию динамики взаимодействия поступательного и вращательного движения чувствительного элемента и моментов, воздействующих на ротор со стороны электростатической системы возбуждения. С использованием динамических моделей Эйлера, при установившихся первичных колебаний ротора, получена следующая модель динамики вторичных колебаний двумерного ММГ для осей жестко связанных с ротором.
?х, ?у - абсолютные угловые скорости вращения корпуса прибора по экваториальным осям;
МВХР, МВУР - совокупности возмущающих внешних моментов по осям жестко связанных с ротором;
g - относительный угол поворота ротора по оси возбуждения;
A - осевой момент инерции ротора;
Bx, By - экваториальные моменты инерции ротора по осям вторичных колебаний;
ca, cb и ??, ?? - коэффициенты упругих и демпфирующих моментов по осям вторичных колебаний;
МД - момент, создаваемый датчиком момента возбуждения гироскопа.
В ММГ-RR типа выходная информация определяется углами a1, b1 поворота плоскости ротора относительно корпуса прибора. При a<<g и b<<g связь углов поворота ротора соответствует операторам вращения:
(2) где (3)
Таким образом, приведение уравнения (1) к осям корпуса прибора сводится фактически к замене переменной j на переменную j1. Такая замена не является простым кинематическим преобразованием, а принципиально изменяет структуру моментов, поскольку полный набор операторов вращения имеет запись:
Подставляя операторы (4) в уравнение (1), получаем уравнения динамики двумерного ММГ роторного типа в осях жестко связанных с корпусом прибора:
Дальнейшие преобразования сводятся к выделению главной линейной части уравнений ММГ по координатам a1 и b1 и малых нелинейных моментов, порождаемые левыми и правыми частями уравнения (5), обусловленные асимметрией его параметров. В результате получаем:
В уравнениях (6) малые нелинейные моменты, образованные левой и правой частями уравнения (5) изза асимметрии параметров гироскопа:
Дополнительно решена задача динамики взаимодействия поступательного и вращательного движения чувствительного элемента гироскопа при линейных скоростях и ускорениях движения объекта. Результат решения представлен уравнениями:
, , (7)
, , (8)
. где x, y, z - малые линейные относительные перемещения центра масс чувствительного элемента;
m - масса ротора гироскопа;
Сх, Су, Cz и ?x, ?y, ?z - коэффициенты упругих и демпфирующих сил;
Vx, Vy, Vz - абсолютные линейные скорости движения корпуса прибора;
?х, ?у - абсолютные угловые скорости вращения корпуса прибора.
В уравнениях (7), (8) и определяют моменты и силы более высокого порядка малости. Уравнения (7) определяют моменты, которые следует учесть в правых частях уравнений динамики вращательного движения ротора. Основная часть этих моментов, линейным образом связанная с координатами x, y, z и их производными, зависит от линейных скоростей и ускорений движения корпуса прибора. Координаты x, y, z и их производные устанавливаются при интегрировании системы уравнений (8) поступательного движения центра масс чувствительного элемента. Совместно уравнения (7) и (8) определяют основное влияние динамики поступательного движения центра масс чувствительного элемента на динамику ММГ как по оси возбуждения, так и по его выходным осям. Именно эти уравнения позволяют определить возмущающие ММГ моменты, связанные в частности с разножесткостью подвеса чувствительного элемента.
Уравнения (6), (7), (8) образуют замкнутую систему позволяющую исследовать динамику и основные погрешности двумерного ММГ роторного типа.
В третьей главе рассмотрена динамика изолированного измерительного канала гироскопа. Для анализа вторичных колебаний ММГ, образующихся при взаимодействии переносной измеряемой абсолютной скорости вращения объекта и установившихся резонансных первичных колебаний ?(t)=?0sin?t запишем уравнение изолированного движения чувствительного элемента в следующей форме:
(9) где a - угол вторичных колебаний, и - частота и коэффициент затухания вторичных колебаний, С=cg0n - приведенный к размерности угловой скорости постоянный параметр, - измеряемая угловая скорость, ?=A/B
Введя в рассмотрение параметры и , определяющие абсолютную и относительную расстройку частоты вторичных собственных колебаний по отношению к частоте возбуждения, и принимая , <<1, <<1 вынужденная составляющая решения уравнения (9) имеет запись: , (10) где
Коэффициент ?c определяет синфазную составляющую выходных вторичных колебаний, фаза которых в зависимости от знака измеряемой угловой скорости либо совпадает с фазой резонансных первичных колебаний, либо сдвинута по отношению к ним на угол -p. При резонансной настройке вторичных колебаний, когда и , именно эти колебания несут информацию об измеряемой угловой скорости .
Коэффициент ?k определяет квадратурную составляющую вторичных колебаний, фаза которых в зависимости от знака расстройки сдвинут на угол по отношению к фазе первичных колебаний. Очевидно, что коэффициент , определяющий полную амплитуду выходных колебаний
(11)
При l1 преобладает квадратурная составляющая, амплитуда которой при увеличении приближается к полной амплитуде вторичных колебаний; при l=1 амплитуды синфазной и квадратурной составляющих вторичных колебаний совпадают.
Обобщенную структурную схему изолированного канала ММГ можно представить в следующем виде:
Структурная схема изолированного канала ММГ
В этой схеме чувствительный элемент представлен двумя механическими колебательными контурами первичных и вторичных колебаний. В верхней части схемы в упрощенном виде показана автогенераторная система возбуждения гироскопа. В контуре первичных колебаний присутствует положительная обратная связь по знаку относительной скорости . Нижняя часть схемы представляет контур вторичных колебаний. На входе колебательного звена этого контура момент гироскопической реакции Мг образуется при наличии измеряемой угловой скорости W. Момент Мг является амплитудно-модулируемым за счет механической модуляции (ММ) на опорной частоте n колебаний угловой скорости первичных колебаний . Измерительная часть вторичного контура, вырабатывающий выходной сигнал прибора U?, включает емкостной преобразователь угла a, демодулятор ДМ, буферный усилитель (Кк) и два фильтра Ф1 и Ф2. Демодулятор выделяет сигнал UK, пропорциональный квадратурной составляющей колебаний a(t), так как в качестве опорного напряжения использует сигнал, первая гармоника которого соответствует колебаниям напряжения . Фильтр Ф1(s) предназначен для подавления в выходном сигнале высокочастотных колебаний на частотах 2n, 2n±D. Фильтр Ф2(s) предназначен для уменьшения колебательности частотной характеристики прибора. Так же в измерительной части вторичного контура пунктиром показана возможность формирования выходного сигнала через синфазную составляющую колебаний a(t).
Дальнейшее исследование динамика ММГ целесообразно проводить с использованием аналитического аппарата, разработанного для систем управления с амплитудной модуляцией.
При наличии измеряемой угловой скорости, изменяющейся по гармоническому закону
, (12)
уравнение движения в измерительном канале микромеханического гироскопа записывается в виде: , (13) где - модулированный гармонический входной сигнал, в котором W0 и W - амплитуда и частота гармонической измеряемой угловой скорости.
Уравнению соответствует передаточная функция: . (14)
Частотные характеристики для огибающих синфазной и квадратурной составляющих вторичных колебаний будут иметь следующий вид: ; (15)
где и - частотные характеристики правой и левой боковых полос.
На рис. 4 приведены амплитудные характеристики, построенные по выражениям (15) при ?=3000 Гц, ?2=10-3, ?=(10, 50, 100) Гц, ?0=1°, ?=1,6, из которых следует, что резонансные пики характеристик соответствуют частотам ?1=? и ?2=2?, а антирезонансный пик находится на комбинированной частоте
Частотные характеристики по огибающим вторичных колебаний: а) квадратурной составляющей, б) синфазной составляющей
В низкочастотной области, определяющей полосу частот пропускания прибора, амплитуды квадратурной составляющей колебаний для всего диапазона частот расстройки ?=(1?100) Гц на (10?30) Дб превосходит амплитуду синфазной составляющей. Для этого диапазона частот ? характеристики (15) допускают простую для инженерной практики аппроксимацию, соответствующую передаточным функциям: где (16)
Построенные частотные характеристики по передаточным функциям (16), полностью совпадают с предыдущими в области низких частот.
Таким образом, при построении информационного канала ММГ следует отдать предпочтение съему информации по квадратурной составляющей вторичных колебаний. Именно в этом случае частота расстройки ? определяет полосу частот пропускания прибора, а амплитуда квадратурной составляющей вторичных колебаний существенно превосходит амплитуду синфазной составляющей.
Аппроксимированные частотные характеристики по огибающим вторичных колебаний: а) квадратурной составляющей, б) синфазной составляющей
Исследования влияния взаимосвязи симметричных каналов вторичных колебаний ММГ на его характеристики при установившихся первичных колебаниях в соответствии с (6) проведем по уравнениям:
(17) где , , .
Уравнение (17) показывают, что взаимосвязь выходных каналов гироскопа, при описании его динамики в координатах и , определяются гироскопическими моментами , моментами радиальной коррекции , перекрестными моментами со стороны датчика момента возбуждения
Для дальнейшего исследования была составлена схема моделирования в пакете MATLAB 6.5., которая представлена на рис. 6. Исходные числовые параметры модели задаются в процессе эксперимента и формируются в М-файле в виде набора констант, отвечающих условиям эксперимента. Две аналогичные по структуре субсистемы alfa и betta моделируют соответствующие каналы углового движения с учетом их взаимного влияния. Выходами этих субсистем являются углы и угловые скорости колебаний ЧЭ в системе осей 0x и 0y корпуса прибора.
Модельный эксперимент, позволяющий определить уровень полной амплитуды вторичных колебаний в перекрестном канале, по сравнению с колебаниями в основном канале при действии на него измеряемой угловой скорости, проводился при следующих параметрах ?у=1°/с, ?х=0, ?=3000 Гц, ?=50 Гц, ?2=10-3, ?0=1°, ?=1,6, с2. Результаты моделирования приведены в таблице 1, из которых следует, что наибольшее влияние на динамику двумерного ММГ оказывают взаимосвязи каналов по гироскопическим моментам. При этом максимальная перекрестная чувствительность двумерного ММГ составляет 0,47%.
Схема моделирования двумерного ММГ роторного типа
Таблица 1
Моменты и взаимосвязи каналов Амплитуда колебаний в канале ?, угл. сек. Амплитуда колебаний в канале ?, угл. сек. Смещение относительно нуля в канале ?, угл. сек.
Без взаимосвязей 0,16 0 0
0,1587,43?10-4-1,27?10-4
0,1580,95?10-6-0,82?10-6
0,1581,175?10-6-0,99?10-6
Все взаимосвязи 0,159 7,43?10-4 -1,35?10-6
В четвертой главе проведено исследование основных информационных характеристик ММГ, к которым относятся масштабные коэффициента преобразования измеряемой угловой скорости и полоса частот пропускания прибора. При проектировании ММГ необходимо особое внимание уделять обеспечению стабильности этих параметров в процессе эксплуатации. Преобладание той или иной составляющей вторичных колебаний определяется безразмерным параметром ?=?0/?2 .
С точки зрения динамики гироскопа при ?>>1 расстройка ? определяет его полосу частот пропускания, так как ? является сопрягающей частотой колебательного звена соответствующего частотной характеристике прибора. Стабильность механического коэффициента преобразования может нарушаться в процессе эксплуатации за счет изменения параметров g0 , ?, ?, x2. При относительных погрешностях параметров гироскопа получаем следующую оценку погрешности коэффициента преобразования: , где (18)
Индексом «н» отмечены номинальные значения параметров.
При построении прибора со съемом информации по квадратурной составляющей вторичных колебаний, наибольший вес в формировании погрешности коэффициента преобразования измеряемой угловой скорости имеет погрешность амплитуды первичных колебаний ??0 и погрешность частоты расстройки вторичных колебаний DD.
Обеспечение максимальной чувствительности ММГ к измеряемой угловой скорости обеспечивается в области малых расстроек вторичных колебаний при ?<<1, когда амплитуда синфазных вторичных колебаний является преобладающей. В этой ситуации механический коэффициент преобразования определяется выражением:
При относительных погрешностях параметров гироскопа получаем следующую оценку погрешности коэффициента преобразования: , где (19)
При построении гироскопа со съемом информации по синфазной составляющей вторичных колебаний (при lн<<1) стабилизация коэффициента преобразования прибора требует стабилизации амплитуды и частоты первичных колебаний (стабилизации амплитуды кинетического момента гироскопа H0=Ag0n) и стабилизации коэффициента затухания вторичных колебаний x2.
В соответствии с принципом формирования контура автогенераторного возбуждения первичных колебаний чувствительного элемента ММГ, представленном в схеме рис. 3, уравнение его движения можно записать в виде: , где (20) и - собственная частота недемпфированных колебаний ЧЭ и коэффициент затухания; [1/с2В2]; - конструктивный коэффициент датчика момента гребенчатой структуры; b, d0, ri, n - высота, боковой зазор, радиус i-ого элемента, общее число пар элементов гребенчатой структуры; U - напряжение питания.
Для приближенной оценки амплитуды установившихся автоколебаний воспользуемся методом гармонической линеаризации правой части уравнения (20).
В результате получаем следующее условие существования незатухающих автоколебаний ЧЭ:
Таким образом, амплитуды скорости и угла первичных колебаний будут определяться зависимостями: и (21)
В условиях наличия шумов в сигнале момент реверсии двигателя раскачки, определяемый по факту изменения скорости , оказывается неопределенным, поскольку отношение сигнал/шум в управляющем сигнале переключающего устройства стремится к нулю. Улучшение характеристик переключателя и, следовательно, всей схемы автовозбуждения может быть достигнуто включением в схему пороговых элементов. В этом случае уравнение амплитуды установившихся колебаний будут иметь вид: ; , (22) где - относительная величина порога переключения.
Поскольку амплитуда колебаний ЧЭ является параметром, определяющим метрологические характеристики ММГ, важным представляется реализация схемотехнических решений, позволяющих стабилизировать ее на заданном уровне. Решение этой задачи возможно при формировании системы регулирования амплитуды согласно функциональной схеме приведенной на рисунке.
Функциональная схема системы стабилизации амплитуды колебаний ЧЭ ММГ
В диссертации разработана схема моделирования контура первичных колебаний, реализуемая в среде Matlab. Схема предусматривает параллельную работу контуров автогенераторного возбуждения, стабилизации амплитуды первичных колебаний и управления их частотой.
Огибающие первичных колебаний при совместной работе контура автогенераторного возбуждения и стабилизации их амплитуды; б). Огибающие первичных колебаний при параллельной работе всех трех контуров
Рисунок 8а. иллюстрирует форсированное автогенераторное возбуждение первичных колебаний для значение коэффициентов затухания ?=0,5·10-4 (1), ?=10-4 (2), ?=0,25·10-4 (3). Для всех этих вариантах общее время возбуждения не более 0,35 с, а использование в системе стабилизации амплитуды ПИ-регулятора обеспечивает астатизм установившихся колебаний на заданном уровне равном 1 градусу. При заданном напряжении источника питания U=10 В, система возбуждения и стабизации автоматически обеспечивает поддержание установившихся колебаний напряжениями U=4,8 В (1), U=6,7 В (2), U=3,3 В (3).
Особенности технологии объемной обработки кремния приводят к существенному разбросу частоты расстройки первичных и вторичных колебаний в ММГ и, что как следует из выражения (18), к неопределенности масштабного коэффициента прибора. При известных частотах собственных колебаний первичного и вторичного контуров, в диссертации решена задача обеспечения требуемой частоты расстройки ?з, определяющей основные информационные характеритиски ММГ. Решение задачи основано на введение в структуру прибора контуров положительной или отрициательной обратной связи, которая может быть реализована как в канале первичных так и вторичных колебаний. При использовании в ММГ роторного типа в контуре вторичных колебаний датчика момента плоской конструкции не удается обеспечить линейное управление дополнительной электростатической жесткостью подвеса. Однако и в этом случае при создании между обкладками датчика момента постоянного опорного напряжения U0, по оси вторичных колебаний образуется дополнительная отрицательная электростатическая жесткость с коэффициентом, определяемым выражением: [Нм], (23) где S - площадь перекрытия электродов;
r - средний радиус расположения электродов;
d0 - зазор между электродами.
Естесственно что, при различных знаках исходной расстройки, возможности этого принципа управления жесткостью подвеса ограничены.
Для контура первичных колебаний, в котором используется дифференциальный датчик момента гребенчатой структуры, отсутствует зависимость величины создаваемого момента от управляемой координаты g. В этом случае, при линеаризации характеристики датчика момента с помощью подачи на электроды гребенчатой структуры опорного напряжения U0, коэффициент передачи датчика момента по отношению к управляющему напряжению определяется зависимостью: [Нм/В], где b, d0, n - высота, боковой зазор и число элементов гребенчатой структуры;
ri - радиус i-го элемента гребенчатой структуры.
При управлении частотой собственных первичных колебаний дополнительный контур положительной или отрицательной обратной связи должен обладать общим коэффициентом передачи: , где Кпg - коэффициент преобразователя угла g;
n - коэффициент усиления.
Для достижения требуемой частоты расстройки Dз, коэффициент обратной связи выбирается из условия: [Нм], (24) где n - исходная частота первичных колебаний;
Dy - управляемая частота расстройки;
со - коэффициент жесткости подвеса в контуре первичных колебаний.
При этом следует иметь в виду, что незначительно изменяется коэффициент передачи первичного контура и его коэффициент затухания в соответствии с зависимостями:
(25)
В (24) и (25) верхний знак соответствует увеличению частоты первичных колебаний контуром отрицательной обратной связи, а нижний знак - уменьшению частоты первичных колебаний контуром положительной обратной связи. Незначительные изменения коэффициентов передачи и затухания в контуре первичных колебаний не приводит к изменению требуемой амплитуды первичных колебаний, которая стабилизируется параллельно работающим контуром стабилизации этой амплитуды. Работоспособность приведенного принципа управления частотой расстройки подтвержден результатами моделирования при одновременной работе контуров автогенерации, стабилизации амплитуды и регулировки частоты расстройки.
На рис. 8б. для варианта ?=0,5·10-4 показано развитие первичных колебаний при параллельной работе контуров автогенераторного возбуждения, стабилизации амплитуды колебаний и управления их частотой. При исходной частоте ?=3000 Гц ее регулировка в пределах ±50 Гц приводит к незначительному изменению динамики развития колебаний. Однако, эта регулировка оказывает существенное влияние на изменения частоты ? расстройки первичных и вторичных колебаний.
В пятой главе исследованы основные погрешности двумерного ММГ, обусловленные влиянием контура первичных колебаний на его выходные характеристики. Это влияние проявляется в следствии инструментальных погрешностей изготовления чувствительного элемента, таких как его статическая и динамическая неуравновешенность и перекос главной оси гироскопа. При наличии этих погрешностей по выходным осям в двумерной схеме формируются следующие возмущающие воздействия:
(26) где x0, y0, z0 - смещение центра масс чувствительного элемента по выходным осям и по оси возбуждения первичных колебаний; E, D - центробежные моменты инерции чувствительного элемента; ?x, ?y - малые углы перекоса главной оси гироскопа. Из уравнений (26) следует, что только первые слагаемые в правых частях образуют моменты, действующие на частоте ? первичных колебаний.
В этом случае уравнения динамики ММГ по выходным осям принимают вид:
(27)
Из уравнений (26), (27) следует, что погрешности ММГ обусловленные моментами ?МХ, ?МY являются квадратурными по отношению к полезной реакции прибора на измеряемые угловые скорости ?у, ?z. Для основного рассматриваемого варианта получения выходной информации в схеме демодуляции по квадратурной составляющей вторичных колебаний, обеспечивающего компромисс между полосой частот пропускания прибора и величиной масштабного коэффициента преобразования допустимые значения рассматриваемых инструментальных погрешностей соответствующих допустимым погрешностям ?уд и ?хд получены соотношения:
(28) радиус инерции ротора.
Отметим, что оценки (28) применимы и к одномерным ММГ роторного типа с одной выходной осью. Оценки (28) позволяют определить допустимые погрешности изготовления чувствительного элемента. Для примера, для двумерного ММГ при ?хд= ?уд и равным 0,1 град/с или 0,01 град/с, при коэффициенте демпфирования ?2х= ?2у=10-4, частоте расстройки ?х=?у=50 Гц, ?x= ?y=1,6, rи=1 мм и частоте ?=3000 Гц допустимые погрешности изготовления чувствительного элемента характеризуются данными приведенными в таблице 2.
Таблица 2 wхд=wyд=0,1 град/с wхд=wyд=0,01 град/с
Е/А, D/A 1,5?10-5 1,5?10-6
, м20,15?10-121,5?10-14 ?х, ?у, угл. сек 5,09 0,5
Из таблицы 2 следует, что в данном примере допустимые значения смещения центра масс лежат в пределах 0,38?0,12 мкм, допустимые перекосы главной оси гироскопа не должны превышать 5 ? 0,5 угл. сек и допустимые отношения центробежных моментов инерции к главному моменту инерции ротора по оси возбуждения должны быть в пределах 1,5·10-5?1,5·10-6.
Основные результаты и выводы
1. Составлены уравнения динамики двумерного ММГ RR-типа, частным случаем которых являются уравнения ММГ с одной осью чувствительности.
2. Исследована динамика взаимодействия поступательного и вращательного движения чувствительного элемента.
3. Получено описание динамики двумерного ММГ RR-типа как в осях, связанных с ротором, так и в осях корпуса прибора. Это дает возможность расчета упругого подвеса гироскопа и исследования его информационных характеристик.
4. Получены выражения для масштабных коэффициентов преобразования измерительных каналов для вариантов съема информации как по квадратурной так и по синфазной составляющих вторичных колебаний, а так же условия их стабилизации.
5. Получены передаточные функции для синфазной и квадратурной составляющих вторичных колебаний как для систем с амплитудной модуляцией так и их инженерные аппроксимации.
6. Для основного варианта съема информации по квадратурной составляющей вторичных колебаний решена задача стабилизации масштабного коэффициента преобразования.
7. Получены точные выражения для оценки возмущающих моментов, обусловленные влиянием контура первичных колебаний на измерительные каналы. Показано, что такие факторы как статическая и динамическая неуравновешенность ротора и перекос осей двумерного ММГ являются источником квадратурной погрешности. Проведена количественная оценка этих погрешностей.
8. Разработан программный комплекс в среде Matlab двумерного ММГ RR-типа, позволяющий вести исследования его статических и динамических характеристик.
Список публикаций по теме диссертации
1. Овчинникова, Н.А. Обобщенная модель динамики микромеханических гироскопов роторного типа/ Л.А. Северов, Н.А. Овчинникова // Гироскопия и навигация. 2007. №4.
2. Овчинникова, Н.А. Информационные характеристики микромеханических гироскопов на основе кремниевой технологии МЭМС /Н.А. Овчинникова, Л.А. Северов, А.И. Панферов, С.К. Золотарев, В.К. Пономарев // Изв. вузов. Приборостроение. 2011. Т.54, №8. С. 12-22.
3. Овчинникова, Н.А. Динамические модели двумерного микромеханического гироскопа/ Л.А. Северов, Н.А. Овчинникова // Сб. докладов Международного симпозиума аэрокосмических приборных технологий (АПТ «04)/ ГУАП. СПБ. 2003.
4. Овчинникова, Н.А. Метрологические характеристики микромеханического гироскопа роторного типа/ Н.А. Овчинникова // Сб. докл. Восьмой научной сессии аспирантов и соискателей ГУАП, посвященной Всемирному дню космонавтики/ ГУАП. СПБ. 2005.
5. Ovchinnikova, N.A. Methods for increasing the accuracy of micromechanical gyroscopes // L.A. Severov, V.K. Ponomarev, A.I. Panferov, N.A. Ovchinnikova // Proceedings. 2005 International Conference Physics and Control. 2005. Digital Object Identifier 10.1109/PHYCON.2005.1514036, ieeexplore.ieee.org.
6. Овчинникова, Н.А. Влияние перекрестных связей на метрологические характеристики микромеханического гироскопа/ Н.А. Овчинникова // Сб. трудов 8-й конференции молодых ученых «Навигация и управление движением»/ ОАО «Концерн «ЦНИИ «Электроприбор». СПБ. 2006.
7. Овчинникова, Н.А. Формирование сигналов управления в микромеханических гироскопах компенсационного типа/ Н.А. Овчинникова // Сб. докладов XV Международного научно-технического семинара «Современные технологии в задачах управления, автоматики и обработки информации»/ Алушта. 2006.
8. Овчинникова, Н.А. Технология MEMS - основа создания новой линейки сенсоров параметров движения/ Л.А. Северов, Н.А. Овчинникова // Сб. трудов XXV Межведомственной НТК «50 лет космодрому «Плесецк»: история создания, этап становления и перспективы развития»/ Мирный. 2007.
9. Овчинникова, Н.А. Микромеханические электростатические драйверы/ Н.А. Овчинникова // Сб. докладов XVI Международного научно-технического семинара «Современные технологии в задачах управления, автоматики и обработки информации»/ Алушта. 2007.
10. Овчинникова, Н.А. Оценка влияния сигналов управления на метрологические характеристики микромеханических гироскопов компенсационного типа/ Н.А. Овчинникова // Сб. трудов Международного молодежного научно-производственного форума «Будущее высоких технологий и инноваций за молодой Россией», проводимого в рамках XV Международной выставки-конгресса «Высокие технологии. Инновации. Инвестиции»/ СПБ. 2009.
12. Ovchinnikova, N.A. Angular velocity sensors based on ring resonator/ L.A. Severov, V.K. Ponomarev, A.I. Panferov, N.A. Ovchinnikova // Proceedings of 2010 International Symposium on Inertial Technology and Navigation/ Nanjing.