Научная статья на тему 'ИССЛЕДОВАНИЕ СВОБОДНЫХ КОЛЕБАНИЙ МИКРОМЕХАНИЧЕСКОГО ГИРОСКОПА С УЧЕТОМ НЕОРТОГОНАЛЬНОСТИ ОСЕЙ ТОРСИОНОВ'

ИССЛЕДОВАНИЕ СВОБОДНЫХ КОЛЕБАНИЙ МИКРОМЕХАНИЧЕСКОГО ГИРОСКОПА С УЧЕТОМ НЕОРТОГОНАЛЬНОСТИ ОСЕЙ ТОРСИОНОВ Текст научной статьи по специальности «Физика»

CC BY
69
10
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
Ключевые слова
ГИРОСКОП R-R-ТИПА / ПРЕЦЕССИЯ ГИРОСКОПА / ОЦЕНКА ПОГРЕШНОСТИ ГИРОСКОПА / МИКРОМЕХАНИЧЕСКИЙ ГИРОСКОП / СВОБОДНЫЕ КОЛЕБАНИЯ

Аннотация научной статьи по физике, автор научной работы — Сайпулаев М.Р., Меркурьев И.В., Соловьёв А.В., Тарасов А.Н.

Введение. Исследован процесс свободных колебаний чувствительного элемента микромеханического гироскопа R-R-типа рамочной конструкции разработки Научно-исследовательского института прикладной механики им. академика В. И. Кузнецова с учетом неортогональности осей торсионов. Изучено влияние инструментальной погрешности изготовления на точность гироскопа на подвижном основании в случае свободных колебаний. Целью работы являлось повышение точности прибора посредством разработки математической модели микромеханического гироскопа R-R-типа с учетом неортогональности осей торсионов и исследования влияния указанной погрешности на точность прибора. Актуальность задачи повышения точности микромеханических гироскопов связана с повышением точности инерциальных навигационных систем, основанных на микромеханических датчиках.Материалы и методы. Предложены новая математическая модель, описывающая динамику гироскопа с учетом инструментальной погрешности изготовления прибора, и формула для оценки погрешности гироскопа. Приведены зависимости переменных состояния, полученные по результатам моделирования и на основе эксперимента. В исследовании использованы методы теоретической механики и асимптотические методы, в том числе формализм Лагранжа и методика осреднения Крылова-Боголюбова.Результаты исследования. Разработана новая математическая модель динамики гироскопа с учетом неортогональности осей торсионов. Получены решение уравнений малых колебаний чувствительного элемента гироскопа и оценка угла прецессии для случая подвижного основания. Проведен сравнительный анализ разработанной модели с экспериментальными данными, полученными в случае свободных колебаний чувствительного элемента гироскопа при неподвижном основании. По результатам анализа подтверждена адекватность построенной математической модели. Сформированы аналитические выражения, демонстрирующие тот факт, что неортогональность осей торсионов приводит к перекрестному влиянию амплитуд первичных колебаний на амплитуды вторичных колебаний чувствительного элемента и появлению дополнительной погрешности в показаниях угловой скорости при работе гироскопа в свободном режиме.Обсуждение и заключения. Полученные результаты могут быть использованы для повышения точности прибора с помощью алгоритма аналитической компенсации погрешности гироскопа и методики идентификации параметров математической модели.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Сайпулаев М.Р., Меркурьев И.В., Соловьёв А.В., Тарасов А.Н.

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

STUDY ON FREE OSCILLATIONS OF A MICROMECHANICAL GYROSCOPE TAKING INTO ACCOUNT THE NONORTHOGONALITY OF THE TORSION AXES

Introduction. The paper is devoted to the study on free oscillations of the sensing element of a micromechanical R-Rtype gyroscope of frame construction developed by the Kuznetsov Research Institute of Applied Mechanics, taking into account the nonorthogonality of the torsion axes. The influence of the instrumental manufacturing error on the accuracy of a gyroscope on a movable base in the case of free oscillations is studied. The work objective was to improve the device accuracy through developing a mathematical model of an R-R type micromechanical gyroscope, taking into account the nonorthogonality of the torsion axes, and to study the influence of this error on the device accuracy. The urgency of the problem of increasing the accuracy of micromechanical gyroscopes is associated with improving the accuracy of inertial navigation systems based on micromechanical sensors.Materials and Methods. A new mathematical model that describes the gyroscope dynamics, taking into account the instrumental error of manufacturing the device, and a formula for estimating the error of a gyroscope, are proposed. The dependences of the state variables obtained from the results of modeling and on the basis of the experiment are presented. Methods of theoretical mechanics and asymptotic methods, including the Lagrange formalism and the Krylov-Bogolyubov averaging method, were used in the research.Results. A new mathematical model of the gyroscope dynamics, taking into account the nonorthogonality of the torsion axes, is developed. The solution to the equations of small oscillations of the gyroscope sensing element and the estimate of the precession angle for the case of a movable base are obtained. A comparative analysis of the developed model and the experimental data obtained in the case of free oscillations of the gyroscope sensing element with a fixed base is carried out. The analysis has confirmed the adequacy of the constructed mathematical model. Analytical expressions are formed. They demonstrate the fact that the nonorthogonality of the torsion axes causes a cross-influence of the amplitudes of the primary vibrations on the amplitudes of the secondary vibrations of the sensing element, and the appearance of an additional error in the angular velocity readings when the gyroscope is operating in free mode.Discussion and Conclusions. The results obtained can be used to improve the device accuracy using the algorithm for analytical compensation of the gyroscope error and the method for identifying the mathematical model parameters.

Текст научной работы на тему «ИССЛЕДОВАНИЕ СВОБОДНЫХ КОЛЕБАНИЙ МИКРОМЕХАНИЧЕСКОГО ГИРОСКОПА С УЧЕТОМ НЕОРТОГОНАЛЬНОСТИ ОСЕЙ ТОРСИОНОВ»

МЕХАНИКА MECHANICS

H) Check for updates

УДК 531.383 https://doi.org/10.23947/2687-1653-2021-21-3-231-238

Исследование свободных колебаний микромеханического ИН^ЛЕ

гироскопа с учетом неортогональности осей торсионов

М. Р. Сайпулаев И12, И. В. Меркурьев 2, А. В. Соловьёв \ А. Н. Тарасов 1

1 Филиал АО «ЦЭНКИ» — ФГУП «Научно-исследовательский институт прикладной механики им. академика В. И. Кузнецова» (г. Москва, Российская Федерация)

2 ФГБОУ ВО «Национальный исследовательский университет «МЭИ» (г. Москва, Российская Федерация) И saypulaevmr@,mail. т

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

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

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

Ключевые слова: гироскоп R-R-типа, прецессия микромеханический гироскоп, свободные колебания.

гироскопа, оценка погрешности гироскопа,

Для цитирования: Исследование свободных колебаний микромеханического гироскопа с учетом неортогональности осей торсионов / М. Р. Сайпулаев, И. В. Меркурьев, А. В. Соловьёв, А. Н. Тарасов // Advanced Engineering Research. — 2021. — Т. 21, № 3. — С. 231-238. https://doi.org/10.23947/2687-1653-2021-21-3-231-238

© Сайпулаев М. Р., Меркурьев И. В., Соловьёв А. В., Тарасов А. Н., 2021

и

S К

а

*

Study on free oscillations of a micromechanical gyroscope taking into account the nonorthogonality of the torsion axes

M. R. Saypulaev И1'2,1. V. Merkuryev 2, A. V. Solovyev ', A. N. Tarasov 1

1 Kuznetsov Research Institute of Applied Mechanics, TsENKI division (Moscow, Russian Federation)

2 National Research University «Moscow Power Engineering Institute» (Moscow, Russian Federation) И saypulaevmr@mail. ru

Introduction. The paper is devoted to the study on free oscillations of the sensing element of a micromechanical R-R-type gyroscope of frame construction developed by the Kuznetsov Research Institute of Applied Mechanics, taking into account the nonorthogonality of the torsion axes. The influence of the instrumental manufacturing error on the accuracy of a gyroscope on a movable base in the case of free oscillations is studied. The work objective was to improve the device accuracy through developing a mathematical model of an R-R type micromechanical gyroscope, taking into account the nonorthogonality of the torsion axes, and to study the influence of this error on the device accuracy. The urgency of the problem of increasing the accuracy of micromechanical gyroscopes is associated with improving the accuracy of inertial navigation systems based on micromechanical sensors.

Materials and Methods. A new mathematical model that describes the gyroscope dynamics, taking into account the instrumental error of manufacturing the device, and a formula for estimating the error of a gyroscope, are proposed. The dependences of the state variables obtained from the results of modeling and on the basis of the experiment are presented. Methods of theoretical mechanics and asymptotic methods, including the Lagrange formalism and the Krylov-Bogolyubov averaging method, were used in the research.

Results. A new mathematical model of the gyroscope dynamics, taking into account the nonorthogonality of the torsion axes, is developed. The solution to the equations of small oscillations of the gyroscope sensing element and the estimate of the precession angle for the case of a movable base are obtained. A comparative analysis of the developed model and the experimental data obtained in the case of free oscillations of the gyroscope sensing element with a fixed base is carried out. The analysis has confirmed the adequacy of the constructed mathematical model. Analytical expressions are formed. They demonstrate the fact that the nonorthogonality of the torsion axes causes a cross-influence of the amplitudes of the primary vibrations on the amplitudes of the secondary vibrations of the sensing element, and the appearance of an additional error in the angular velocity readings when the gyroscope is operating in free mode. Discussion and Conclusions. The results obtained can be used to improve the device accuracy using the algorithm for analytical compensation of the gyroscope error and the method for identifying the mathematical model parameters.

Keywords: gyroscope R-R type, gyro precession, gyro error estimation, micromechanical gyroscope, free oscillations.

For citation: M. R. Saypulaev, I. V. Merkuryev, A. V. Solovyev, A. N. Tarasov. Study on free oscillations of a micromechanical gyroscope taking into account the nonorthogonality of the torsion axes. Advanced Engineering Research, 2021, vol. 21, no. 3, pp. 231-238. https://doi.org/10.23947/2687-1653-2021-21-3-231-238

Введение. Разработка высокоточных микромеханических инерциальных датчиков, в том числе микромеханических гироскопов (ММГ), применяемых для решения задач навигации и управления движением летательных аппаратов и мобильных роботов, является актуальной задачей приборостроения [1]. К достоинствам ММГ относятся малые масса и размеры, а также низкая стоимость по сравнению с гироскопами, основанными на других физических принципах. Однако к главным недостаткам ММГ относятся изменчивость его метрологических характеристик и низкая точность измерений параметров углового движения объекта (угловой скорости и угла поворота). Принцип функционирования вибрационных гироскопов основан на свойстве маятника Фуко сохранять плоскость малых колебаний неподвижной в инерциальном пространстве [2].

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

* ¡^

нелинейность, обусловленная конечными колебаниями чувствительных элементов (ЧЭ) или физической > нелинейностью, связанной с особенностями системы управления колебаниями [2-6].

Исследования, посвященные динамике и конструированию ММГ, были опубликованы и в работах зарубежных авторов [6-8]. Так, например, в публикациях [6, 8] получена формула для оценки уходов гироскопа, основанная на использовании разработанной математической модели движения, которая описывает медленное изменение тороидальных координат колебаний ЧЭ. В работах [0, 7] обсуждаются вопросы изготовления ММГ и анализируются уравнения его малых колебаний. В статье [0] составляются уравнения движения ММГ с угловым (R-R-типа) и линейным (L-L-типа) осцилляторными видами движения ЧЭ. В

й

а £ .й

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

При проектировании ММГ разработчики стремятся использовать явление внутреннего резонанса в системе, обусловленное совмещением собственных частот колебаний ЧЭ [3, 4]. Однако в работах [0, 7] отмечается, что погрешности технологии изготовления, неизвестные и непредсказуемые отклонения элементов конструкции от проектных положений приводят к дополнительным погрешностям в измерениях прибора.

Для повышения точности измерения угловой скорости ММГ поставлена цель: исследовать свободные колебания (при отсутствии управления) ЧЭ ММГ R-R-типа с учетом эффектов, возникающих из-за неортогональности осей торсионов. Этот дефект появляется в виду несовершенства технологии изготовления прибора. Ставятся задачи разработки новой математической модели динамики ММГ с учетом неортогональности осей торсионов, оценки ухода прибора и описания влияния неортогональности осей торсионов на динамику ЧЭ ММГ.

Материалы и методы. Рассматривается модельная конструкция вибрационного ММГ R-R-типа — конструкция с промежуточной рамкой в соответствии с классификацией из источника [3]. Кинематическая схема гироскопа (рис. 1) реализована в виде двухстепенного карданного подвеса ЧЭ.

Рис. 1. Конструктивная схема прибора: 1 — основание (корпус); 2 — промежуточная (внешняя) рамка; 3 — чувствительный элемент, состоящий из сбалансированной пластины и инерционной массы; 4 — торсионы

Для описания положения ЧЭ введем системы координат (рис. 2), связанных: с корпусом прибора — 0ХУ2 ; с внешней рамкой упругого подвеса — Ох1у1zx; с сбалансированной пластиной — 0ху2 . Причем 02 является осью чувствительности гироскопа, а система координат Ох2у2 г2 отличается от системы Ох1 у1 zl поворотом на постоянный угол неортогональности всех торсионов. В представленных системах начало координат соответствует точке О и находится в геометрическом центре сбалансированной пластины.

2

3

1

X

Рис. 2. Системы координат

О T3

a £ л

В системе Охху1 z1 зададим осевые моменты инерции промежуточной рамки 31х, 31у, а в системе Охух — осевые моменты инерции ЧЭ 32х, 12у, Отметим, что в данной работе оси систем координат Ох1 у1 и Охух считаются главными центральными осями инерции промежуточной рамки и ЧЭ соответственно.

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

ОХУ2 —— Ох1 у1 ^ —— Ох2 у2 z2 —— Оху2,

X 21 У!

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

Составим уравнения динамики ЧЭ ММГ в форме уравнений Лагранжа 2-го рода [9, 10]:

d I дЬ ) дЬ дФ d

dt I dá) да dá dt

dL ГФ

(1)

Sp) Sp 5(3

где L = T - П — функция Лагранжа; Т и П — кинетическая и потенциальная энергии системы соответственно; Ф — диссипативная функция, характеризующая потери на внутреннее трение. Выражения для этих величин имеют вид:

т = \(J2*al +-/2v®v +J2-®;) + \(JiA2 +(J» sin2 a +Jlz eos2 a)Q2),

1 1 (2)

<s> = \d«á2+\d^2' и = \сУ'+\с^2-

где da, dp — коэффициенты трения; ca, cp — коэффициенты жесткости торсионов.

Выражения для проекций cor, cov, ох угловой скорости ЧЭ на подвижные оси х. v, г имеют вид: сох = ácosPcos¿>-Qcosasinp + Qsinacospsinc»,

cov = p + Qsinacos¿>-c¡:sin¿>, (3)

со. = ásinpcosc) +QcosacosP + QsinasinPsin¿>. Учитывая малость углов а , р и 8 , тригонометрические функции в выражениях (2) и (3) от указанных углов можно заменить их разложениями в ряд Тейлора, ограничившись слагаемыми до первого порядка малости. Тогда получим уравнения малых колебаний из уравнений движения (1), записанных с учетом выражений (2) и (3).

Используя формализм Лагранжа [9] в случае постоянной угловой скорости основания, получаем уравнения малых колебаний ЧЭ, записанные с точностью до слагаемых первого порядка малости в виде:

а + со;а=/РР--+ Р + со^З = -/2£М--^Р + 8а, (4)

(Я, ./2 0/1

где введены следующие обозначения (аналогично тому, как это выполнено в статье [11]):

J2 ж + J 2 y J2 z . J 2 ж + J 2 y J2 z

J1 =-, J2 =-, a>a=„,

J + J J 'a \í j . j

•Лж T 2 ж 2 y 2 ж

. d , dR

O ~ =-a--q -1 = —L_

Ha t j , j \' Hp j

J + J2 ж ) J2 y

Ц2 Здесь ]х, ]2 — безразмерные моменты инерции упругого подвеса; аа, тр и , ^ — собственные частоты

колебаний и добротности по углам а , В соответственно. "¡3 При выводе уравнений колебаний (4) угловая скорость корпуса гироскопа О считалась малой

¡и относительно собственной частоты соа, т.е. << , а также угол 8 полагался малой величиной, т. е. 8 << 1.

/ Отметим, что в уравнениях (4) отброшены слагаемые, обусловленные наличием геометрической нелинейности ММГ. Влияние нелинейности геометрии движения ЧЭ на динамику ММГ Я-Я-типа описано в монографии [5].

Принимая во внимание, что правые части уравнений (4) являются малыми возмущениями, т. е. а + ауаа = 0(е), то с точностью до слагаемых первого порядка малости можно записать: а =-охаа + О (е). 234 Таким образом, из правой части уравнений (4) исключаются вторые производные углов а и р .

Рассмотрим случай изотропного упругого подвеса, т. е. равенства собственных частот колебаний и равных добротностей:

Юa=Юp=Юo, <2а= ^р = 0,

где ю0 — характерное значение собственной частоты колебаний; 0 — характерное значение добротности.

Отметим, что случай наличия разнодобротнсти Ф ^) и малой разночастотности (аа Ф «) при

исследовании свободных колебаний ЧЭ ММГ рассмотрен в статье [11]. При введенных обозначениях и принятых допущениях запишем уравнения движения ЧЭ с точностью до слагаемых первого порядка малости в безразмерном виде:

á+co\а = j^fy-O 'со,,« - — ¿>со„р,

Л

(5)

Р + со„Р = -jflá - О 'со0Р - Sala .

Отметим, что система уравнений (5) приводится к стандартной форме записи регулярно возмущенной системы дифференциальных уравнений с одной быстрой угловой переменной [12, 14]. Одним из распространенных способов нахождения решения регулярно возмущенных систем является использование асимптотических методов разделения движений [12-15].

Решение нелинейных уравнений (5) получим с помощью методики осреднения Крылова — Боголюбова [13], а в качестве медленно меняющихся переменных будем использовать переменные Ван-дер-Поля [12] pi, qup2, q<

а = рх sin(co0/) + í71 cos(co0/), á = щр1 cos(a0t)-a0q1 sin(co0/), P = p2 sin (co0/) + í72 cos(co0/), P = (ú0p2 cos(co0/)-co0í72 sin(co0/).

С помощью процедуры осреднения [14, 15] по явно входящему времени, получаем осредненную систему дифференциальных уравнений, разрешенную относительно производных медленных переменных Ван-дер-Поля:

, 1 ^ _ j Ш j8 Л' = --Q Pl +Т— Р2 42,

2' 1

2ю„

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

2 j 2

, 1 _! j, Q j8 ?í = -tQ ъ +T-Pi,

2ю„

, * „ ! j2 Q 8

p2 =--q pipi--4

j i Q

2 j 2 8

(6)

q2 =--Q '42-^r- 4 +T Pv

2 2a0 2 2 2a0 2

Штрихом в уравнениях (6) обозначено дифференцирование по безразмерному времени т = .

Полученная модель в виде линейных дифференциальных уравнений описывает свободные колебания ЧЭ гироскопа на подвижном основании. Решение системы уравнений (6) можно записать в виде:

Г (уРю -№ю) •

P1 (t) = exP| -2Q

41 (т) = exp | - —

P10COS + ^ j Sin

^cos (f^fr)+j Sin (ypTfr)

P2 (t) = exP| - 2Q

42 (t) = exp| - —

P20 COS(sjy2 +Г2*) 420 COS (ф2 +Г2*)

4vl + у2

(VPW + У4ю )

2 + у2

Nú - "УЛо )„

(N+^T)

(7)

где р10 = рх (0), р20 = р2 (0), д10 = ^ (0), д20 = д2 (0) — начальные условия; безразмерная угловая скорость основания прибора V = ^]\]\ о/( 2га0) ; параметр, характеризующий неортогональность осей торсионов

г=4Г^К 2лДГ).

Вторые слагаемые в формулах (7) характеризуют перекрестное влияние первичных колебаний на вторичные колебания и наоборот. Отметим, что в случае ортогональных осей торсионов, когда у = 0,

решение (7) совпадает с результатами работы [5].

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

ей И S К

ей

*

(D

АЛгапсев. Еп&пеегтдResearch 2021. Т. 21, №3. С. 231-238. 2687-1653

Й О

ТЗ ^

и

и

О, С -Й

Результаты исследования. Для валидации разработанной модели сравним результаты моделирования, вычисленные по формулам (7), с экспериментальными данными. Измерительная информация получена с помощью системы наблюдения. В качестве измерительной информации электростатических датчиков имеем переменные Ван-дер-Поля р1, дь р2, д2.

В эксперименте использован образец прибора со следующими параметрами математической модели при неподвижном основании (V = 0): Q = 3856, К = ]'2 = 1, у = 0,2-105, в качестве начальных условий для переменных Ван-дер-Поля выбраны значения равные измерениям в начальный момент времени: р10 = 13,467-10-3, <710 = 20,429-10-, р20 = 0,787-10-3, ц20 = 1,172-10-3.

Значение параметра у соответствует углу неортогональности осей торсионов 8 равного одной угловой секунде. Графическое представление зависимостей переменных Ван-дер-Поля (г), ^ (г), р2 (г), 42 (г) медленно меняющихся с течением безразмерного времени, приведено на рисунке 3.

Рг103 - 0

0,6

0,4

0,2

42-10 1,2

1,0

0,8

0,6

0,4

0,2

\ :

чч v л чч

;

5 000 10 000 15 000

20 000

5 000 10 000 15 000 20 000 Рис. 3. Зависимости переменных Ван-дер-Поля: 1 — результаты моделирования; 2 — экспериментальные данные

Графики (рис. 3) демонстрируют значительное совпадение зависимостей для переменных р2, 42 , полученных по результатам моделирования, с экспериментальными данными. Зависимости для переменных р1 , 41 , полученные по результатам моделирования, качественно согласуются с данными эксперимента, а наблюдаемые небольшие количественные отклонения могут быть обусловлены нелинейными эффектами, такими как нелинейность геометрии движения ЧЭ [5], или явлениями разнодобротности, разночастотности и погрешностью смещения инерционной массы [11]. Рассмотрение нелинейных эффектов, влияющих на динамику ММГ, при построении математических моделей колебаний ЧЭ увеличивает точность микромеханических датчиков в составе инерциальных навигационных систем [2].

Уход гироскопа из-за нелинейных эффектов и других инструментальных погрешностей будем оценивать с помощью вспомогательного функционала I [5, 6, 8]:

/= ^ЛК (44 + Р Р2 )

(412 + Р12)- Л (42 + Р22):

К

который связан с углом в через соотношение:

в = — агСап (I).

Причем этот параметр пропорционален интегралу от угловой скорости:

в = -Щ1 Jn(ri ) dTi.

2 0

С помощью формулы (8), учитывая решение (7), можно оценивать уход гироскопа, связанный с неортогональностью осей торсионов, возникшей из-за несовершенства технологии изготовления. На рис. 4. представлены зависимости функционала I от безразмерного времени по результатам эксперимента и расчета по формуле (8).

/

0,120 0,118 0,116 0,114 0,112 0,110 0,i08

0

Рис. 4. Зависимость функционала I (г) : I — результаты моделирования; 2 — экспериментальные данные

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

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

Библиографический список

1. Peshekhonov, V. G. Gyroscopic navigation systems: Current status and prospects / V. G. Peshekhonov // Gyroscopy and Navigation. — 20II. — Vol. 2 (3). — P. III-II8.

2. Журавлев, В. Ф. Управляемый маятник Фуко как модель одного класса свободных гироскопов / В. Ф. Журавлев // Известия РАН. Механика твердого тела. — I997. — Вып. 6. — С. 27-35.

3. Распопов, В. Я. Микромеханические приборы / В. Я. Распопов. — Москва : Машиностроение, 2007.

— 400 с.

4. Неаполитанский, А. С. Микромеханические вибрационные гироскопы / А. С. Неаполитанский, Б. В. Хромов. — Москва : Когито-центр, 2002. — I22 с.

5. Меркурьев, И. В. Динамика микромеханического и волнового твердотельного гироскопов / И. В. Меркурьев, В. В. Подалков. — Москва : Физматлит, 2009. — 228 с.

6. Electronic Gain Error Compensation for Whole-Angle Coriolis Vibrating Gyroscopes with High Q Factor / Yongmeng Zhang, Tongqiao Miao, Kechen Guo [et al.] // In: Proc. IEEE Int. Symposium on Inertial Sensors and Щ Systems (INERTIAL). — 2020. — P. I-4. https://doi.org/I0.II09/INERTIAL48I29.2020.9090062 ¡g Apostolyuk, V. Theory and Design of Micromechanical Vibratory Gyroscopes / V. Apostolyuk // MEMS/NEMS Hand- $ book. — 2006. — Vol. I, chapter 6. — P. I73-I95. https://doi.org/I0.I007/0-387-25786-I 6 S

7. Wei Wang. Design of a Novel MEMS Gyroscope Array / Wei Wang, Xiaoyong Lv, Feng Sun // Sensors.

— 20I3. — Vol. I3 (2). — P. I65I-I663.

8. Askari, S. High quality factor MEMS gyroscope with whole angle mode of operation / Sina Askari, Mohammad H. Asadian, Andrei M. Shkel // In: Proc. IEEE Int. Symposium on Inertial Sensors and Systems (INERTIAL). — 2018. — Vol. 12 (6). — P. 141-144. https://doi.org/10.1109/ISISS.2018.8358148

9. Маркеев, А. П. Теоретическая механика / А. П. Маркеев. — Москва, Ижевск : РХД, 2007. — 592 с.

10. Теоретическая механика / С. В. Болотин, А. В. Карапетян, Е. И. Кугушев, Д. В. Трещев. — Москва : Академия, 2010. — 432 с.

11. Сайпулаев, М. Р. Динамика и точность микромеханического гироскопа с учетом смещения инерционной массы / М. Р. Сайпулаев, И. В. Меркурьев // Вестник Бурятского государственного университета. Математика, информатика. — 2020. — Вып. 3. — С. 49-62. https://doi.org/10.18101/2304-5728-2020-3-49-62

12. Asymptotic Methods in Mechanics of Solids / Andrei L. Smirnov, Sergei Filippov, Petr E. Tovstik [et al.] — Birkhauser Basel: International Series of Numerical Mathematics; 2015. — 323 p.

13. Awrejcewicz, J. Introduction to asymptotic methods / J. Awrejcewicz, V. A. Krys'ko. — Boca Raton Fl: Chapman & Hall / CRC; 2006. — 242 p.

14. Burd, V. Method of Averaging for Differential Equations on an Infinite Interval: Theory and Applications / V. Burd. — Chapman & Hall/CRC; 2007. — 343 p.

15. Медведев, Г. Н. Лекции по методу усреднения / Г. Н. Медведев. — Москва : Физический факультет МГУ, 2019. — 109 с.

Поступила в редакцию 28.06.2021 Поступила после рецензирования 19.07.2021 Принята к публикации 19.07.2021

Об авторах:

Сайпулаев Муса Русланович, аспирант кафедры «Робототехника, мехатроника, динамика и прочность машин», ФГБОУ ВО «Национальный исследовательский университет «МЭИ» (111250, РФ, г. Москва, ул. Красноказарменная, д. 14), инженер 1 категории НИИ ПМ им. академика В. И. Кузнецова (111024, РФ, г. Москва, ул. Пруд-Ключики, д. 12А), ORCID: http://orcid.org/0000-0002-5165-654X, saypulaevmr@mail.ru

Меркурьев Игорь Владимирович, заведующий кафедрой «Робототехника, мехатроника, динамика и прочность машин», ФГБОУ ВО «Национальный исследовательский университет «МЭИ» (111250, РФ, г. Москва, ул. Красноказарменная, д. 14), доктор технических наук, доцент, ScopusID: 35422634900, ORCID: http://orcid.org/0000-0001-7682-2228, nir4s@ya.ru

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

Соловьёв Алексей Владимирович, начальник отделения «Инерциальные приборы и датчики», ФГУП «Научно-исследовательский институт прикладной механики им. академика В. И. Кузнецова» (111024, РФ, г. Москва, ул. Пруд-Ключики, д. 12А), кандидат технических наук, ORCID: http://orcid.org/0000-0003-1345-1778, AV. Solovev@russian. space

Тарасов Александр Николаевич, начальник отдела «Инерциальные приборы и датчики», ФГУП «Научно-исследовательский институт прикладной механики им. академика В. И. Кузнецова» (111024, РФ, г. Москва, ул. Пруд-Ключики, д. 12А), кандидат технических наук, ORCID: https://orcid.org/0000-0002-3657-5733, A.Tarasov@russian.space

Заявленный вклад соавторов:

М. Р. Сайпулаев — формирование основной концепции, задачи и цели исследования, построение

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

3 моделирования, подготовка текста, формирование выводов; И. В. Меркурьев — научное руководство, анализ § результатов исследований, доработка текста, корректировка выводов; А. В. Соловьев — обсуждение 2 результатов исследования, доработка текста; А. Н. Тарасов — проверка сравнительного анализа ' ¡3 экспериментальных данных и результатов моделирования, обсуждение результатов исследования.

о, £ л

Все авторы прочитали и одобрили окончательный вариант рукописи.

i Надоели баннеры? Вы всегда можете отключить рекламу.