УДК 534.11 Вестник СПбГУ. Математика. Механика. Астрономия. 2024. Т. 11 (69). Вып. 2 МБС 74Н60
Об эволюции кососимметричных изгибных колебаний круглой пластинки при ее вращении вокруг оси, расположенной в плоскости пластинки*
Н. Ф. Морозов1'2, А. В. Лукин3, И. А. Попов3
1 Санкт-Петербургский государственный университет,
Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9
2 Институт проблем машиноведения Российской академии наук, Российская Федерация, 199178, Санкт-Петербург, Большой пр. В. О., 61
3 Санкт-Петербургский политехнический университет Петра Великого, Российская Федерация, 195251, Санкт-Петербург, ул. Политехническая, 29
Для цитирования: Морозов Н. Ф., Лукин А. В., Попов И. А. Об эволюции кососимметричных изгибных колебаний круглой пластинки при ее вращении вокруг оси, расположенной в плоскости пластинки // Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия. 2024. Т. 11 (69). Вып. 2. С. 354-370. https://doi.org/10.21638/spbu01.2024.209
В настоящей работе выполняется построение и исследование модели связанных плоско-поперечных колебаний круглой тонкой пластинки с концентричным отверстием при действии кориолисовых и центробежных сил инерции, вызванных вращением системы вокруг оси, расположенной в плоскости пластинки. Уравнения колебаний в частных производных получены с применением вариационного принципа Гамильтона — Остроградского. В предположении малости угловой скорости вращения по отношению к частоте рабочей кососимметричной изгибной формы колебаний пластинки найдено приближенное аналитическое решение как для радиальной и окружной, так и для поперечной компонент поля перемещений в режиме свободных колебаний. С помощью проекционного метода Галеркина задача была сведена к системе двух линейных дифференциальных уравнений второго порядка для модальных координат взаимно ортогональных базисных кососимметричных форм колебаний пластинки. Обнаружено, что режим начально возбужденных гармонических колебаний при наличии вращения преобразуется в режим амплитудно-модулированных биений. Найдены аналитические выражения как для частоты медленной огибающей биений, так и для относительной глубины их амплитудной модуляции. Показана принципиальная возможность определения модуля проекции вектора угловой скорости на плоскость пластинки по измеряемой величине частоты огибающей. Исследована задача о выборе оптимальной геометрической формы резонатора с точки зрения максимизации чувствительности системы к изменениям величины угловой скорости вращения. Рассмотрен вопрос об определении направления проекции вектора угловой скорости на плоскость пластинки по измеряемой глубине амплитудной модуляции режима биений.
Ключевые слова: микромеханический гироскоп, твердотельный волновой гироскоп, многоосевой МЭМС гироскоп, инерциальная навигация, параметрические колебания.
* Исследование выполнено за счет гранта Российского научного фонда №21-71-10009, ЫЛрБ://rscf.ru/project/21-71-10009/.
© Санкт-Петербургский государственный университет, 2024
1. Введение. Одной из важнейших задач измерительного приборостроения на современном этапе является создание малогабаритных инерциальных навигационных систем, обеспечивающих определение параметров ориентации и координат летательных, надводных и подводных аппаратов, наземных транспортных средств, роботов и др. Ключевым направлением разработок здесь является создание микроэлектромеханических датчиков инерциальной информации — акселерометров и гироскопов тактического и навигационного классов точности [1]. Наряду с высокими требованиями к точностным характеристикам датчиков, определяющее значение имеют параметры температурной стабильности и устойчивости систем к внешним воздействиям (вибрационным, ударным, акустическим) [2].
В последнее десятилетие в мировой индустрии микросистем интенсивное развитие получило направление разработки многоосевых датчиков — вибрационных микромеханических систем с архитектурой чувствительного элемента (ЧЭ), позволяющей одновременно определять более одной компоненты измеряемой характеристики движения (векторов переносного ускорения или угловой скорости). В отношении микромеханических гироскопов здесь следует указать работы передовых инженерных и исследовательских центров [3-8]. Одним из перспективных решений при разработке конструкции ЧЭ подобных сенсоров является использование особых свойств режимов упругих колебаний в телах с осевой симметрией — кольцах, пластинах, дисках [4, 8, 9]. Данный подкласс датчиков в известном смысле примыкает к обширной группе разработок в области микромасштабных твердотельных волновых гироскопов с осесимметричным резонатором [10-12], хотя и опирается на существенно иные принципы съема выходного сигнала, чем в системах со свободной гироскопической прецессией рабочей формы колебаний.
Вопросам колебаний осесимметричных упругих тел (колец, пластинок, оболочек) в неинерциальных системах отсчета посвящено значительное число работ. Кроме уже упомянутых выше задач вибрационной гироскопии, подобные постановки широко известны в роторной динамике [13-18]. В подавляющем большинстве исследований этого класса рассматривается случай совпадения направления вектора угловой скорости с осью симметрии упругого тела (диска), что приводит к возникновению в нем предварительного напряженного состояния за счет действия центробежных сил инерции и изменению характеристик свободных и вынужденных изгибных колебаний.
Отдельную группу исследований составляют работы по динамике тонких круглых пластинок с различными типами граничных условий на внешней и внутренней (при наличии у пластинки отверстия) кромках. Рассматриваются как линейные, так и нелинейные геометрические постановки [19, 20] для осесимметричных и кососим-метричных колебаний [21, 22], а также бегущих упругих волн [23].
Насколько известно авторам предлагаемой работы, в научной литературе ранее не исследовались задачи динамики круглых пластинок (в частности, их кососиммет-ричных колебаний) при совместном действии гироскопических и центробежных сил инерции, вызванных вращением системы вокруг некоторой оси, не совпадающей с осью симметрии пластинки. Наиболее близкой к названной постановке является работа [24], в которой изучаются аналогичные вопросы для тонкого упругого кольца.
В настоящей работе исследуются особые свойства гироскопически связанных плоско-поперечных колебаний тонкой круглой пластинки и рассматривается возможность их использования для определения величины и направления проекции на плоскость пластинки вектора угловой скорости подвижного объекта, на котором
установлен резонатор. Иными словами, изучаются перспективы применения чувствительного элемента указанной формы в новых архитектурах двухосевых микроэлектромеханических гироскопов.
2. Математическая модель. Рассматривается задача об эволюции кососим-метричных изгибных колебаний круглой пластинки с концентричным отверстием, защемленной по внутреннему краю и свободной по внешнему краю, в неинерциаль-ной системе отсчета, движение которой представляет собой вращение вокруг некоторой оси, расположенной в плоскости пластинки. Графическая схема исследуемой задачи приведена на рис. 1.
Рис. 1. Графическая схема рассматриваемой задачи.
Для материала пластинки принимается модель однородного изотропного линейно-упругого тела с объемной плотностью р, модулем Юнга Е и коэффициентом Пуассона V. Толщину, внутренний и внешний радиусы пластинки обозначим Ь,Ъ и а соответственно. Примем, что вращение системы совершается с постоянной угловой скоростью О, направленной вдоль оси Ох декартовой системы координат Охух, связанной с пластинкой; в силу осевой симметрии конструкции ориентация этой оси в плоскости пластинки может быть задана произвольно.
В качестве механической модели пластинки будем рассматривать геометрически линейную бессдвиговую модель.
Для вывода уравнений движения системы применим интегральный вариационный принцип Гамильтона —Остроградского:
0= [ (ёТ - ёП) (И, (1)
Jo
где Т и П — кинетическая и потенциальная (упругая) энергии пластинки; ( и ё — операторы дифференцирования и варьирования соответственно; £ — время.
Вариация кинетической энергии может быть записана следующим образом [25]:
ёТ = - ! ! рп • ёпалах, (2)
где и — вектор перемещений материальных точек пластинки; интегрирование выполняется по объему тела (¿Л — элемент площади поверхности пластинки). Вариация упругой энергии в общем виде может быть представлена как
т= / («711 fell + «225£22 + «33^33 + «23^23 + «13^13 + «12^12) dAdz, (3) JzJa
где 7ij, £ij — компоненты тензоров напряжений и упругих деформаций в выбранной системе координат. Конкретный вид потенциальной энергии с учетом допущений принятой модели будет приведен далее.
Уравнения движения будем записывать в связанной с пластинкой цилиндрической системе координат r, в, z.
Векторное поле перемещений имеет вид
z
и = (и - zwr)ji + (v--we).72 + wj3, (4)
r
где u,v,w — радиальная, окружная и поперечная компоненты вектора перемещений; J1,32,33 — базисные векторы цилиндрической системы координат; нижние индексы обозначают дифференцирование по соответствующей пространственной координате. Вектор сообщаемой пластинке угловой скорости имеет следующий вид:
О = Í1 (cosej1 — sinej2) . (5)
Вектор абсолютной скорости материальных точек пластинки, записанный в базисе используемой неинерциальной цилиндрической системы координат, определяется выражением
du du ^
где оператор d/dt обозначает покомпонентное дифференцирование вектора без учета фактора подвижности рассматриваемого базиса.
С учетом (6) в базисе {31,2,3} вектор скорости покомпонентно имеет вид
du z
— = \ й — zwr — ílw sin 6,v--we — Oto cos в, w+
dt r
(v--we j cos в + (и — zwr ) sin в j . (7)
Поступая аналогичным образом, найдем следующее представление для вектора абсолютного ускорения в подвижном цилиндрическом базисе:
d2u dt2
= {ü — zwr — 2Otó sin в — Q2 sin в ^v--we j eos в + (и — zwr) sin в
z
z2 v--we — 2Otó cos в — ir cos в
WW + 2Ü
(v--we j cos в + (и — zwr) sin в , (8)
--we j cos в + (й — zwr) sin в — Q2w}T.
r
r
Подстановка (8) в (2) и интегрирование полученного выражения по толщине пластинки позволяют записать вариацию кинетической энергии в виде
г2п ra
ST = - / {I0 [ü - 20w sin в - О2 sin в (v cos в + и sin в)] Su+
J 0 Jb
+Io [v — 20w sin в — О2 cos в (v cos в + и sin 0)] Sv+ +I0 [w + 2О (ú sin в + V cos в) - О2w] Sw+
+I2
wr — Í2 sin в ( —we cos в + wr sin в r
+ -
1 2 /1
— Í} COS в -We cos в + wr sin в rr
Swr+ Swg }гс1гс!в,
(9)
где /0 = А 12 =
Вариация упругой энергии может быть записана как [25]
/•2п ta
sn = - / {[(rNi)r + Nee - N2] Su + [М2в + (rNe)r + N] Sv+ 0b
+ [(rQi)r + Q2e] Sw - [(rMi)r + Mee - rQi - M2] Swr -
— [М2в + rM6r ~ rQ2 + 2M6] &/;e}drd0+
r
r2n
+
0
+ N6Sv + ( Qi + -M6e ) ¿w - Miówv
rdB,
r=b
(10)
где силы и моменты в поперечных сечениях элементов пластинки определяются выражениями
rn/2
{Ni,N2,Ne} = Wii,*22,012} dz,
J-n/2
Г n/2
{Mi,M2,Me} = z {an,a22,02} dz.
n/2
(11)
Таким образом, применение вариационного принципа Гамильтона — Остроградского (1) приводит к следующей системе уравнений движения пластинки:
1 Н---—^ Н----- = /о \й — 2Qw sin в — Q2 sin в (v cos в + и sin 0)] ,
ОГГО0Г
6 Н---—- Н--- = /о íi) — 2Qw cos в — Q2 cos в (v cos в + и sin 0)] ,
ОГГО0Г
+ + — = 1о [w + 2Q(üsme + icose)-Q2w] , (í2)
or r дв r v '
dMi 1 dMe Mi - M2 --1---яа" ---Qí — h
dr r дв r
dM6 1 дМ2 2Me
dr
дв
+ Q2 = I2
21
—wr + Q sin в —We cos в + wr sin в r
1 2 (1
—we — Í2 cos 0 -we cos в + wr sin в rr
2
r
r=a
r
r
с рядом возможных комбинаций граничных условий вдоль кромок г = Ь,а:
6т = 0
5и = 0 или N = 0,
5у =0 или N6 = 0,
или + 1 ЭМ6 г дв = 0, (13)
6»г = 0 или М1 = 0.
Для изотропного материала внутренние силы и моменты связаны с полем перемещений следующим образом:
Ек 1-1/2
1 V 0
V 1 0
0 0 (1 - V) /2
иг
(уд + и) /г (ид + гуг — у) /г
(14)
М1
М2) = — Мб
ЕН3
12(1 — V2)
1 V 0
V 1 0
0 0 (1 — V)
(тдд + г»г) /г2 (гшгд — »д) /г2
(15)
С учетом (14), (15) система уравнений (12) может быть сведена к трем уравнениям свободных связанных плоско-поперечных колебаний пластинки:
иГг--2 («е + и) Н— (Угд + иг)
1 — V . ,1 — V
- 2г2 (идд + ГУгд - 3Уд - 2и) Н---—иг
= и — 20>и} эт в — О2 эт в (у сов в + и эт в).
1 — V , \ 1 , \ V 1
2 2 + гмг - г;) + - (игд + ГУГГ) + -игд + — + ид)
= у — 2О» соэ в — О2 соэ в (у соэ в + и эт в),
(16) (17)
ДУ4» + /от — /2У2» + 2/0О (их эт в + у соэ в) — /0О2т+
+ /гО^эт0 ^—сов$ Н—соэ0 + тгг втА ) —
--вт в ( —иид соэ в + уиг вт в I +
гг
Н--СОЭ в (^—Юдд СОЭ в--Шд ЭШ в + иигд ЭШ в + тг СОЭ в ) +
Н—втА ( —иид сов в + уиг втА )] = 0, (18)
Ек3
где введены обозначения с2 = , О — щг-и2) •
Предметом дальнейшего исследования является изучение качественных особенностей решений системы уравнений (16)—(18) для частного случая эволюции начально сообщенной пластинке некоторой кососимметричной формы изгибных колебаний
гг
2
С
2
С
при действии гироскопических сил (сил Кориолиса), вызванных малой по отношению к частоте возбужденных изгибных колебаний угловой скоростью вращения системы.
3. Анализ гироскопически связанных плоско-поперечных колебаний пластинки. С учетом вышеуказанных особенностей рассматриваемой постановки выполним ряд упрощений в общей системе (16)—(18):
1) в предположении малой толщины пластинки по отношению к ее характерным размерам в своей плоскости пренебрежем наличием сил инерции на поворот поперечных сечений (отбросим слагаемое ¿у|-У2и>);
2) имея в виду малость величины угловой скорости по сравнению с рабочей частотой возбужденной формы изгибных колебаний и гироскопическую (пропорциональную П) связь компонент вектора перемещений с поперечным перемещением й, пренебрежем центробежными (пропорциональными П2) слагаемыми в уравнениях (16)-(17);
3) с учетом пп. 1, 2 будем считать величиной старшего порядка малости и отбросим из рассмотрения последнюю группу слагаемых в уравнении (18), пропорциональную множителю ^-И2;
4) рассматривая низкочастотные (по сравнению с формами плоских движений пластинки) формы изгибных колебаний, пренебрежем силами инерции в радиальном и окружном направлениях (слагаемые и и V в уравнениях (16) и (17) соответственно).
Кроме того, перейдем к следующим безразмерным величинам, отмеченным сим-волом~:
и = Ни, V = НЪ, й = Нй,
~ (19)
г = аЪ, г = Т~, П = ТП,
где Т2 =
Отбрасывая далее символ ~ в обозначениях безразмерных величин, получим следующую упрощенную систему уравнений колебаний пластинки:
v , ч ч1 — V
2 г2
игг--ñ (ve + и) + - (vrg + ur) + „ (uge + rvrg - 2>ve - 2u) +
+ -—-ur = —2a2ílw sin в, (20) r
2r2 (ue + rur -v) + - (urg + rvrr) + -ure + (vee + ив) = -2а2Пъисовв, (21) yAw + w + 2Q (u sin в + V cos в) — ü2w = 0, (22)
2 a2
где a2 = ф^.
Исследование эволюции некоторой начально возбужденной формы изгибных колебаний выполним с помощью приближенного (проекционного) метода Бубнова — Галеркина. Будем искать решение для поперечной компоненты поля перемещений пластинки в виде
w (r, в, t) = Rnk (r) [n(t) cos пв + Z(t) sin пв], (23)
где п,к — индексы осевой симметрии и радиальной изменяемости рассматриваемой формы колебаний; п(^),С(^) — подлежащие определению модальные координаты; Кпи (г) — функция радиальной координаты, определяемая решением задачи на собственные значения для уравнения свободных изгибных колебаний пластинки на неподвижном основании:
Е(г) = С1 Ъ(кг) + С2 /п(кг) + С3 Уп(кг) + С\Кп(кг). (24)
Здесь 1п,/п и Уп,Кп — обычные и модифицированные функции Бесселя I и II рода соответственно. Параметр к связан с размерной круговой частотой свободных колебаний ш соотношением
л рНа4ш2
«4 = (25)
Значения коэффициентов С^ и параметра к определяются из условий удовлетворения решением (24) граничных условий, которые в рассматриваемом случае пластинки, защемленной по внутреннему краю и свободной по внешнему (рис. 1), имеют вид
г = в : Я = 0, Я' = 0, 2
г = 1 : Д" + -Д'-^Д = 0
— (26)
Я'" + \я" - 1 + 2п*-т2я + (3-;)га2д = о,
где введено обозначение в = Ь/а.
Вид первых двух форм изгибных колебаний пластинки из семейств п = 0 и п =1 для значения параметра в = 0.2 показан на рис. 2 (коэффициент Пуассона V принят равным 0.22).
Примечательной особенностью принятой комбинации граничных условий является пересечение частотных ветвей (п,к) = (0,1) (осесимметричные формы с наименьшей радиальной изменяемостью) и (п,к) = (1,1) (исследуемые кососим-метричные формы) при достижении геометрическим параметром в определенного значения: при в < 0.312, т. е. для полых круглых пластинок с достаточно малым внутренним радиусом, наинизшей собственной частотой обладает кососимметрич-ная форма колебаний. Этот эффект отражен на рис. 3.
Везде далее будем рассматривать кососимметричную форму колебаний с наименьшей окружной и радиальной изменяемостью: (п,к) = (1,1).
Подстановка (23) в (20)—(21) приводит к системе неоднородных уравнений в частных производных относительно и и у, точное решение которой может быть представлено в форме
u (r,9,t) = a2Q -CUi(r)+ ((cos29 - r¡sin2^ U2(r) , v (r,e,t) = a2Q -f¡V1 (r) - (r¡ cos 29 + С sin 29^ V2 (r) ,
Вестник СПбГУ. Математика. Механика. Астрономия. 2024■ Т. 11 (69). Вып. 2
(27) 361
Рис. 2. Формы изгибных колебаний круглой пластинки с отверстием.
Рис. 3. Ветви безразмерного частотного параметра к.
где функции и\{г),и2(г),У1 (т),У2(г) являются решениями следующих граничных задач:
Щ - -Vi + -и[ = Ди(г),
Г.У.: г = ß : r = 1 :
Ui =0,
(28)
U[ + -U! = 0,
ir
Vi" +
1 -V
2r
V -
1 - v
2 r2
Г.У.: r = ß : r = 1 :
Vi = Rn(r),
Vi =0,
(29)
rVi' - Vi = 0,
1 :
Г.У.:
v 2v
Щ + -U2--v2 = o,
1 + vTr, 3 - vTr
+ -T-V2
r2 r2
1-v 9 - v
Vo 2 г 2 2 г2 2
=ß: U2 =0,
Rii(r), Rii(r),
V2 =0, rV2' - V2 + 2U2 = 0.
(30)
Решение приведенных граничных задач может быть получено численно с помощью соответствующих методов вычислительной математики (в частности, методом коллокаций1).
Подставим найденные выражения для радиальной и окружной компонент поля перемещений (27) в уравнение изгибных колебаний (22). Учитывая искомый вид поперечной компоненты поля перемещений (23), запишем проекционные условия метода Галеркина: поэтапно умножим (22) на Дц (г)еов в, Дц(г)вт в и проинтегрируем полученные уравнения по площади пластинки. Таким образом, придем к следующей системе обыкновенных дифференциальных уравнений для модальных координат ф),((1):
(1 + 2a2tt2Ai) п + («и - О2) П = 0, (1 + 2а2ü2Bi) С + (к4и - О2) С = 0,
(31)
где введены обозначения:
Ai =
¿гДп (V1 + lU2 + ly2)dr L1 rR2ndr
Bi = -
jp гДц (Uj + \U2 + \V2) dr
Si rRiidr
(32)
Задача (31) позволяет провести качественный анализ эволюции начально возбужденной (с некоторыми начальными условиями по, Со) кососимметричной формы колебаний при действии гироскопических и центробежных сил.
Примечательно, что гироскопическая по своей физической природе связь движений пластинки в своей плоскости и вне плоскости приводит для рассматриваемого
1 MathWorks MATLAB version: 9.13.0 (R2022b). Natick, Massachusetts: The MathWorks Inc. URL: https://www.mathworks.com (дата обращения: 05.01.2023).
r
r
r
r
случая ориентации вектора угловой скорости к появлению слагаемых, пропорциональных квадрату скорости вращения (см. множители при вторых производных в системе (31)).
Уравнения (31) могут быть записаны в виде
v + rfv = 0, С + = 0, (33)
где безразмерные частоты модальных осцилляторов с окружной изменяемостью по cos в и sin в (см. представление решения для прогиба (23)) определяются выражениями 4 2 4 2
2 К41 — О 2 Кц — О / \
^ = 1 + 2a2íl2Ai' ^ = l + 2a2íí2Bl- ^
Считая угловую скорость вращения пластинки малой по отношению к рассматриваемой круговой частоте изгибных колебаний системы на неподвижном основании, запишем в первом приближении по О2:
ШГ1 = К2п - + J
) К11 N (35)
Для дальнейшего исследования будет удобно перейти от частот к средней
частоте шт и частотному расщеплению 2ДШ по формулам
1 , 2 1 2 К + = «и " 2
4- + а2к2п (25-Ма + < К11
Ü2
(36)
где с учетом структуры коэффициентов Ai,Bi (32) введены обозначения:
¡i R (U2 + V2) dr ¡i rRiiVidr ¡i rRiiUidr S=--¿-j-, 5a = -^-i-, 5b = --^-l-. (37)
2 ¡13 rRiidr J ¡3rRhdr lв rR21dr
Предположим, что в начальный момент времени пластинке сообщено начальное поперечное отклонение по рассматриваемой кососимметричной форме изгибных колебаний, причем угловая ориентация этой формы по отношению к связанной с пластинкой осью Ox (см. рис. 1) определяется некоторым углом 7:
w0 (r, в) = CRii(r) cos (в + 7) = CRii(r) [cos 7 cos в — sin 7 sin в]. (38)
Параметр C характеризует амплитуду начального прогиба; в силу линейности задачи везде далее примем C = 1.
Подчеркнем, что угол 7 определяет, таким образом, ориентацию возбужденной формы колебаний по отношению к направлению вектора угловой скорости Q. Здесь необходимо отметить, что одним из возможных технических применений рассматриваемой задачи является предложение нового принципа работы и конструктивного решения для датчиков угловой скорости (прежде всего — микромеханических вибрационных гироскопов с электростатическими возбудителями колебаний и системами
m
емкостного съема выходного сигнала). В контексте этих приложений направление и величина вектора Q (точнее, его проекции на плоскость пластинки) подлежат определению, а ориентация начально возбужденной формы колебаний чувствительного элемента является заданной. Из сказанного следует, что с учетом осевой симметрии геометрии пластинки мы можем интерпретировать угол 7 как параметр, определяющий искомое направление вектора угловой скорости по отношению к известной начальной ориентации возбужденной кососимметричной формы колебаний.
Сопоставление (38) с (23) показывает, что начальные значения модальных координат в рассматриваемом случае равны по = cos 7, Со = — sin 7, а значит,
Таким образом, свободные изгибные колебания пластинки при действии гироскопических и центробежных сил определяются выражением
которое, с учетом введенных обозначений (36), может быть преобразовано к форме
w (r, в, t) = Rn(r) [cos (в + 7) cos Дшt cos umt + cos (в — 7) sin Дшt sin umt] . (41)
Необходимо отметить, что отдельного рассмотрения требует техническая и алгоритмическая реализация задачи сообщения резонатору заданных начальных условий и поддержания свободных колебаний пластинки с течением времени в условиях физически неизбежной диссипации механической энергии, с одной стороны, и недопустимости внесения возбудителем колебаний окружной асимметрии в систему — с другой. Принципиальное решение этой группы вопросов может быть найдено путем расположения пластинки между двумя эквидистантными плоскими электродными структурами с осевой симметрией, состоящими из электростатических приводов с переменным межэлектродным зазором и емкостных датчиков поперечных перемещений. При работе подобной электродной системы в режиме возбуждения параметрических колебаний резонатора указанные выше условия будут выполнены [11].
Имея в виду разработку модели датчика угловой скорости, рассмотрим зависимость от времени движения точки на внешней (свободной) кромке пластинки в принимаемом известном направлении пучности начально возбужденной формы колебаний в = —7:
W(t) = w(1, —7,t) = Rii(1) [cos Дшt cosumt + cos 27sin Дшt sin umt] . (42)
Как видно из полученного выражения, колебания в точке съема выходного сигнала представляют собой амплитудно-модулированные биения с частотой огибающей Дш и относительной глубиной модуляции 1 — | cos 271. Характер этих колебаний показан на рис. 4.
Из сказанного следует, что возникает принципиальная возможность определения амплитуды и направления вектора угловой скорости Q (или, как отмечено выше, в общем случае — его проекции на плоскость пластинки) по измеряемым характеристикам колебаний в некоторых точках пластинки: частота биений Дш, как видно из (36), пропорциональна квадрату модуля вектора угловой скорости О2; параметр глубины амплитудной модуляции cos 27 связан с направлением этого вектора.
n(t) = cos 7 cos unt, Z(t) = — sin 7 cos uzt.
(39)
w (r, в, t) = Дц (r) [cos 7 cos в cos unt — sin 7 sin в cos uzt],
(40)
О 100 200 300 400 500 600 700 800 900
г
Рис. 4. Вид зависимости от времени колебаний в точке съема выходного сигнала.
Исследуем подробнее выражение (36) для частоты биений Дш. Видим, что коэффициент при О2 есть произведение а2/2 и к21 (5а — 5ъ), где первый множитель определяется конкретными физико-механическими и геометрическими параметрами пластинки, а второй представляет собой безразмерную величину, зависящую лишь от формы пластинки — параметра в отношения ее внутреннего радиуса к внешнему. Таким образом, возникает задача о нахождении оптимальной формы пластинки (оптимальной величины параметра в), соответствующей наибольшему возможному значению множителя чувствительности к22 (6а — 6ь). Результаты расчета зависимостей названных величин от параметра в приведены на рис. 5.
о
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
0 0
а б
Рис. 5. Зависимость параметров влияния вращения от /3: а — множители к^ и ¿а — 3ь; б — полный множитель чувствительности.
Как видно из рисунка, эффект гироскопического расщепления кратных собственных частот рассматриваемых взаимно ортогональных кососимметричных
форм колебаний резонатора проявляется тем сильнее, чем меньше внутренний радиус пластинки по отношению к внешнему. Этот результат можно считать ожидаемым и связанным с пропорциональностью сил Кориолиса эффективной массе, вовлеченной в колебания, которая, в свою очередь, растет с уменьшением радиуса отверстия.
Кратко рассмотрим и вопрос об определении направления проекции вектора Q на плоскость пластинки по измеряемой глубине амплитудной модуляции, связанной с величиной cos27. Данная задача составляет известные трудности, поскольку наблюдение медленной огибающей сигнала в виде биений (см. рис. 4) позволяет определить лишь | cos27|. Кроме того, функция cos27 сама по себе является двузначной на интервале (—■§, ■§) (одному ее значению соответствуют два возможных угловых направления, симметричных относительно 7 = 0). Дополнительная информация, необходимая для однозначного определения 7, вероятно, может быть получена путем синхронных измерений сигналов колебаний в нескольких точках на кромке пластинки и их фазовой демодуляции по высокочастотной гармонике. Конкретное исследование этого вопроса составляет предмет дальнейших исследований.
4. Заключение. В настоящей работе была впервые построена математическая модель связанных плоско-поперечных колебаний круглой тонкой пластинки с концентричным отверстием при действии кориолисовых и центробежных сил инерции, вызванных вращением системы вокруг оси, расположенной в плоскости пластинки. Уравнения колебаний в частных производных получены с применением вариационного принципа Гамильтона — Остроградского. В предположении малости угловой скорости вращения по отношению к частоте рабочей кососимметричной изгибной формы колебаний пластинки найдено приближенное аналитическое решение как для радиальной и окружной, так и для поперечной компонент поля перемещений в режиме свободных колебаний. С помощью проекционного метода Галеркина задача была сведена к системе двух линейных дифференциальных уравнений второго порядка для модальных координат взаимно ортогональных базисных кососимметрич-ных форм колебаний пластинки. Обнаружено, что режим начально возбужденных гармонических колебаний при наличии вращения преобразуется в режим амплитуд-но-модулированных биений. Найдены аналитические выражения как для частоты медленной огибающей биений, так и для относительной глубины их амплитудной модуляции. Показана принципиальная возможность определения модуля проекции вектора угловой скорости на плоскость пластинки по измеряемой величине частоты огибающей. Исследована задача о выборе оптимальной геометрической формы резонатора с точки зрения максимизации чувствительности системы к изменениям величины угловой скорости вращения. Рассмотрен вопрос об определении направления проекции вектора угловой скорости на плоскость пластинки по измеряемой глубине амплитудной модуляции режима биений.
Литература/References
1. Peshekhonov V. G. The Outlook for Gyroscopy. Gyroscopy Navig. 11, 193-197 (2020). https://doi.org/10.1134/S2075108720030062
2. Indeitsev D. A., Belyaev Y. V., Lukin A. V., Popov I. A., Igumnova V. S., Mozhgova N. V. Analysis of imperfections sensitivity and vibration immunity of MEMS vibrating wheel gyroscope. Nonlinear Dynamics 105 (2), 1273-1296 (2021). https://doi.org/10.1007/s11071-021-06664-0
3. Tsai D.H., Fang W. Design and simulation of a dual-axis sensing decoupled vibratory wheel gyroscope. Sensors and Actuators, A: Physical 126 (1), 33-40 (2006). https://doi.org/10.1016 /j.sna.2005.09.004
4. Johari H., Shah J., Ayazi F. High frequency xyz-axis single-disk silicon gyroscope. Proceedings of the IEEE International Conference on Micro Electro Mechanical Systems, 856—859 (2008). https://doi.org/10.1109/MEMSYS.2008.4443791
5. Vigna B. Tri-axial MEMS gyroscopes and six degree-of-freedom motion sensors. Technical Digest — International Electron Devices Meeting IEDM) (2011). https://doi.org/10.1109 /IEDM.2011.6131635
6. Cao H., Zhao R., Cai Q., Shi Y., Liu L. Structural Design and Simulation Analysis of Silicon Micro Triaxial Wheel-ring Gyroscope. International Conference on Sensing, Measurement and Data Analytics in the Era of Artificial Intelligence ICSMD 2020 — Proceedings, 130-132 (2020). https://doi.org/10.1109/ICSMD50554.2020.9261699
7. Takahashi H., Abe K., Takahata T., Shimoyama I. MEMS triaxial gyroscope using surface and sidewall doping piezoresistors. Journal of Micromechanics and Microengineering 30 (10) (2020). https://doi.org/10.1088/1361-6439/ab9e4d
8. Sedebo G.T., Shatalov M.Y., Joubert S.V., Shafi A. A. The Dynamics of a Three-Dimensional Tuning Functionally Graded Plate Gyroscope. Mechanics of Solids 57 (6), 1577-1589 (2022). https://doi.org/10.3103/S0025654422060267
9. Smirnov A. L. Free Vibrations of the Rotating Shells of Revolution. ASME. J. Appl. Mech. 56 (2), 423-429 (June 1989). https://doi.org/10.1115/L3176100
10. Hu Z., Gallacher B.J. Effects of Nonlinearity on the Angular Drift Error of an Electrostatic MEMS Rate Integrating Gyroscope. IEEE Sensors Journal 19 (22), 10271-10280 (2019). https://doi.org/10.1109/JSEN.2019.2929352
11. Indeitsev D.A., Zavorotneva E.V., Lukin A. V., Popov I. A., Igumnova V. S. Nonlinear Dynamics of a Microscale Rate Integrating Gyroscope with a Disk Resonator under Parametric Excitation. Rus. J. Nonlin. Dyn. 19 (1), 59-89 (2023).
12. Maslov A. A., Merkuryev I.V., Maslov D.A., Podalkov V. V. Scale Factor of the Wave Solid-State Gyroscope Operating in the Angular Velocity Sensor Mode. 29th Saint Petersburg International Conference on Integrated Navigation Systems, ICINS (2022). https://doi.org/10.23919/ICINS51784.2022.9815350
13. Luo A. C.J., Mote C.D., Jr. Nonlinear Vibration of Rotating Thin Disks. ASME. J. Vib. Acoust. October 2000 122 (4), 376-383 (2000). https://doi.org/10.1115/L1310363
14. Nayfeh A. H., Jilani A., Manzione P. Transverse Vibrations of a Centrally Clamped Rotating Circular Disk. Nonlinear Dynamics 26, 163-178 (2001). https://doi.org/10.1023/A:1012957024898
15. Hamidzadeh H. R. Non-linear free transverse vibration of thin rotating discs. Proceedings of the Institution of Mechanical Engineers, Part K: Journal of Multi-Body Dynamics 221 (3), 467-473 (2007). https://doi.org/10.1243/14644193JMBD46
16. Angoshtari A., Jalali M. A. On the existence of chaotic circumferential waves in spinning disks. Chaos 17 (2) (2007). https://doi.org/10.1063/L2735813
17. Hu Y., Wang T. Nonlinear resonance of the rotating circular plate under static loads in magnetic field. Chinese Journal of Mechanical Engineering (English Edition) 28 (6), 1277-1284 (2015). https://doi.org/10.3901/CJME.2015.0720.097
18. Kang J. H. Axisymmetric Vibration of Rotating Annular Plate with Variable Thickness Subjected to Tensile Centrifugal Body Force. International Journal of Structural Stability and Dynamics 17 (9), 1750101 (2017). https://doi.org/10.1142/S0219455417501012
19. Vasiliev G. P., Smirnov A. L. Free vibration frequencies of a circular thin plate with nonlinearly perturbed parameters. Izvestiya of Saratov University. New Series. Series: Mathematics. Mechanics. Informatics 21 (2), 227-237 (2021). https://doi.org/10.18500/1816-9791-2021-21-2-227-237
20. Asher A., Gilat R., Krylov S. Natural Frequencies and Modes of Electrostatically Actuated Curved Bell-Shaped Microplates. Applied Sciences 12 (5), 2704 (2022). http://dx.doi.org/10.3390 /app12052704
21. Touze C., Thomas O., Chaigne A. Asymmetric non-linear forced vibrations of free-edge circular plates. Part 1: Theory. Journal of Sound and Vibration 258 (4), 649-676 (2002). https://doi.org/10.1006/jsvi.2002.5143
22. Thomas O., Touze C., Chaigne A. Asymmetric non-linear forced vibrations of free-edge circular plates. Part II: Experiments. Journal of Sound and Vibration 265 (5), 1075-1101 (2003). https://doi.org/10.1016/S0022-460X (02)01564-X
23. Nayfeh T. A., Vakakis A. F. Subharmonic travelling waves in a geometrically non-linear circular plate. Pergamon Int. J. Non-Linear Mechanics 29 (2), 233-245 (1994). https://doi.org/10.1016/0020-7462(94)90042-6
24. Eley R., Fox C. H. J., Mcwilliam S. The dynamics of a vibrating-ring multi-axis rate gyroscope. Proceedings of The Institution of Mechanical Engineers Part C-journal of Mechanical Engineering Science 214, 1503-1513 (2000). https://doi.org/10.1243/0954406001523443
25. Nayfeh A.H., Pai P. F. Linear and Nonlinear Structural Mechanics. Wiley Series in Nonlinear Science. Wiley-Interscience (2004).
Статья поступила в редакцию 1 февраля 2023 г.;
доработана 7 октября 2023 г.; рекомендована к печати 9 ноября 2023 г.
Контактная информация:
Морозов Никита Федорович — д-р физ.-мат. наук, проф.; n.morozov@spbu.ru Лукин Алексей Вячеславович — канд. физ.-мат. наук, доц.; lukin_av@spbstu.ru Попов Иван Алексеевич — popov_ia@spbstu.ru
Gyroscopically coupled in-plane and transverse vibrations of an annular free-clamped microplate*
N. F. Morozov1'2, A. V. Lukin3, I. A. Popov3
1 St. Petersburg State University, 7—9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation
2 Institute for Problems in Mechanical Engineering of the Russian Academy of Sciences, 61, Bolshoy pr. V. O., St. Petersburg, 199178, Russian Federation
3 Peter the Great St. Petersburg Polytechnic University,
29, ul. Politekhnicheskaya, St. Petersburg, 195251, Russian Federation
For citation: Morozov N.F., Lukin A. V., Popov I. A. Gyroscopically coupled in-plane and transverse vibrations of an annular free-clamped microplate. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy, 2024, vol. 11(69), issue 2, pp. 354-370. https://doi.org/10.21638/spbu01.2024.209 (In Russian)
In this work, we construct and study a model of coupled plane-transverse oscillations of a circular thin plate with a concentric hole under the action of Coriolis and centrifugal inertia forces caused by the rotation of the system along an axis located in the plane of the plate. Equations of vibrations in partial derivatives are obtained using the variational principle of Hamilton — Ostrogradsky. Assuming the smallness of the angular velocity of rotation with respect to the frequency of the working skew-symmetric flexural form of the plate oscillations, an approximate analytical solution is found for both the radial and circumferential, and transverse components of the displacement field in the free oscillation mode. Using the Galerkin projection method, the problem was reduced to a system of two second-order linear differential equations for modal coordinates of mutually orthogonal basic skew-symmetric modes of the plate vibrations. It is found that the regime of initially excited harmonic oscillations in the presence of rotation is transformed into the regime of amplitude-modulated beats. Analytical expressions are found both for the frequency of the slow beat envelope and for the relative depth of their amplitude modulation. The fundamental possibility of determining the modulus of the projection of the angular velocity vector onto the plane of the plate from the measured value of the envelope frequency is shown. The problem of choosing the optimal geometric shape of the resonator from the point of view of maximizing the sensitivity of the system to changes in the value of the angular velocity of rotation is studied. The question of determining the direction of the
*The research was funded by the Russian Science Foundation grant no. 21-71-10009, https://rscf.ru/project/21-71-10009/.
projection of the angular velocity vector onto the plane of the plate from the measured depth of the amplitude modulation of the beat regime is considered.
Keywords: micromechanical gyroscope, rate-integrating gyroscope, triaxial MEMS gyroscope, inertial navigation, parametric ocsillations.
Received: February 1, 2023 Revised: October 7, 2023 Accepted: November 9, 2023
Authors' information:
Nikita F. Morozov — n.morozov@spbu.ru Alexey V. Lukin — lukin_av@spbstu.ru Ivan A. Popov — popov_ia@spbstu.ru
ХРОНИКА
18 октября 2023 г. на заседании секции теоретической механики им. проф. Н. Н. Поляхова в Доме ученых им. М. Горького (Санкт-Петербург) выступили кандидат физ.-мат. наук, доцент А. С. Кулешов и аспирант Н. М. Видов (механико-математический факультет МГУ им. М. В. Ломоносова, Москва) с докладом на тему «Нелинейные эффекты вблизи многообразия равновесий неголономных систем».
Краткое содержание доклада:
В 1986 г. Я. В. Татаринов заложил основы теории слабо неголономных систем — механических систем с неголономными связями, содержащими малый параметр е и таких, что: 1) при е = 0 связи интегрируемы, т.е. получается семейство голоном-ных систем, зависящее от постоянных интегрирования; 2) эти голономные системы — вполне интегрируемые гамильтоновы системы. При е = 0 поведение таких систем можно рассматривать при помощи асимптотических методов, представляя его как трансгрессию: сочетание движения слегка модифицированной голономной системы с медленным изменением былых констант. В докладе представлено исследование эффекта трансгрессии в двух задачах неголономной динамики: 1) движение тяжелого тонкого твердого стержня по поверхности прямого кругового цилиндра; 2) качение тяжелого однородного шара по неподвижной поверхности в окрестности наинизшей точки (эллиптического типа) данной поверхности.