Finite element investigation of the effectiveness of the tubular piezoelectric vibratory gyroscope depending on the type of polarization and boundary conditions

Abstract


In the present paper, a dynamic behavior of piezoelectric vibratory gyroscope in the form of a hollow piezoelectric cylinder with two pairs of electrodes placed crosswise on the outer side surface has been analyzed. In the case of one fixed end, two variants of the piezoceramic material polarization have been considered, namely, full radial polarization and partial radial polarization only under the electrodes on the outer side surface. For a completely polarized material (in addition to the case of a cantilevered end) two other variants of the fixations simulating the conditions of hinged opening were also considered. The behavior of the harmonic oscillation of gyroscope has been studied in the framework of linear theory of piezoelectricity (electroelasticity) which takes into account mechanical damping and rotational effects in the relative coordinate system. All of the investigated variants admit the availability of electrically active bending oscillation modes in two perpendicular planes which can be controlled by given electric potentials on the two pairs of electrodes. In such configurations when the gyroscope operates near to the corresponding resonance frequencies the primary flexural vibrations are generated in one plane and the secondary flexural motions are created due to the axial gyroscope rotation in the perpendicular plane. These secondary oscillations can be a measure to determine the value of the rotation frequency. The finite element method, ANSYS finite element package and specially designed computer programs written in the macro language APDL ANSYS were used for numerical calculations. The results of computational experiments have shown that the variant with one fixed end and with a full radial polarization of the piezoceramic material gives the largest maximum of the output potential in the presence of gyroscope rotation. It has been found that the variant of the gyroscope with the fixations simulating the conditions of hinged opening is also quite promising for application.

Full Text

Введение Пьезоэлектрические вибрационные гироскопы успешно используются в современной технике в качестве миниатюрных датчиков угловой скорости [1-4]. Для эффективной работы вибрационного гироскопа желательно наличие у него по крайней мере двух ярко выраженных мод колебаний в различных, а лучше всего в перпендикулярных плоскостях. Тогда при колебаниях устройства на одной моде гироскопическое вращение за счет эффекта Кориолиса вызовет дополнительные движения на другой моде колебаний. Регистрация этого вторичного движения и позволяет определить частоту вращения. В пьезоэлектрическом вибрационном гироскопе, как правило, и возбуждение первичных колебаний, и регистрация вторичных осуществляются за счет пьезоэлектрического эффекта. Для этого в гироскопическом датчике необходимо иметь электроды как для возбуждения, так и для приема колебательных движений. Пьезоэлектрические вибрационные гироскопы могут быть выполнены полностью из пьезоэлектрического материала, но могут иметь и сложные конструкции, состоящие из упругих и пьезоэлектрических сред. Одними из наиболее известных вибрационных пьезоэлектрических гироскопов являются конструкции камертонного вида с двумя или с четырьмя вилками (H-тип). Несмотря на то, что камертонные гироскопы известны достаточно давно [5], их исследования продолжаются и в настоящее время [6, 7]. Актуальными остаются задачи разработки и анализа других эффективных конструкций вибрационных гироскопов как канонических типов в форме балок [8, 9], пластин [10], цилиндров [11] и других, так и усложненных конфигураций [12-15]. Для расчета пьезоэлектрических вибрационных гироскопов во многих работах с успехом использовался метод конечных элементов, техника которого применительно к пьезоэлектрическим вращающимся телам была описана в [16-18]. При этом если частота вращения много меньше частоты первичных колебаний, то для решения можно также использовать метод конечных элементов, но в совокупности с асимптотическим разложением по малому параметру отношения частот [19]. В [20-26] рассматривался трубчатый вибрационный гироскоп из радиально поляризованной пьезокерамики с двумя парами электродов, расположенными на внешней боковой поверхности, и с одним сплошным электродом на внутренней поверхности. В этих работах трубчатый гироскоп был исследован с использованием приближенных теорий изгиба, метода конечных элементов (МКЭ) и экспериментальных данных. Тем не менее ряд важных вопросов остался невыясненным. Так, были недостаточно изучены модификации устройства при различных граничных условиях на торцах, допускающих наличие электрически активных кратных собственных частот с модами колебаний в двух взаимно перпендикулярных плоскостях. Кроме того, для пьезокерамического материала в полой трубке с электродами, не полностью покрывающими ее боковые поверхности, возможны различные структуры неоднородности, обусловленные технологией процесса поляризации. Последнее, как известно [11], может существенно изменить характеристики гироскопа. Отмеченные вопросы изучаются в настоящей работе на основе [25, 26]. При этом здесь рассматриваются варианты гироскопа с одной и той же геометрией трубки, тогда как в [25, 26] длина трубки подбиралась так, чтобы обеспечить приблизительное равенство основной резонансной частоты изгибных колебаний. Таким образом, в настоящей работе проводятся исследования по двум основным направлениям. Во-первых, изучается эффективность гироскопа с полной радиальной поляризацией пьезокерамики в сравнении с гироскопом с частичной поляризацией только под электродами на внешней поверхности. Во-вторых, анализируется влияние условий закрепления на торцах на возможности регистрации вращения на соответствующих резонансных частотах. Данные вопросы актуальны при использовании гироскопа рассматриваемого типа на практике, поскольку позволяют выбрать наиболее эффективные технологии его поляризации и условия закрепления. Подобные исследования могут представлять интерес и для многих других видов вибрационных пьезоэлектрических гироскопов. 1. Описание задачи Пусть полый цилиндр имеет длину , внешний радиус и толщину . Введем относительную декартовую систему координат , направив ось вдоль центральной оси цилиндра так, чтобы гироскоп занимал положение (рис. 1, а). Будем использовать также цилиндрическую систему координат , связанную с исходной системой соотношениями , , . Примем, что на внешнюю боковую поверхность , , нанесены четыре электрода ( ) с углами раствора , причем , . Расположение электродов в плоскости поперечного сечения показано на рис. 1, б. Наконец, имеется еще электрод с нулевым потенциалом , который покрывает всю внутреннюю боковую поверхность цилиндра: , ; . Рис. 1. Геометрия трубчатого гироскопа (а) и его поперечное сечение (б) Fig. 1. Geometry of the tube gyroscope (a) and its cross section (b) Гироскоп работает в режиме установившихся колебаний на некоторой частоте . При подаче входного потенциала на электроды и соответственно за счет пьезоэффекта и различий в направлениях векторов наводимого электрического поля и поляризации пьезоэлектрический материал под сжимается, а под расширяется или наоборот, и цилиндр совершает изгибные колебания в плоскости (первичное движение). При вращении всего устройства вокруг оси с постоянной угловой скоростью сила Кориолиса, действующая перпендикулярно направлениям первичного движения и вектора угловой скорости , вызывает деформации в направлении (вторичное движение). В режиме вторичного движения материал под электродом сжимается, а под - растягивается или наоборот, что приводит к изгибным колебаниям в плоскости . В итоге на свободных электродах и за счет пьезоэффекта появляются выходные электрические потенциалы , регистрация которых позволяет определить частоту вращения гироскопа . (Отметим, что для более прецизионного анализа нужно, вообще говоря, ввести в рассмотрение внешние электрические цепи [27], позволяющие определять электрические потенциалы и токи на выходных электродах, но эти усложненные постановки задач пьезоэлектричества здесь не рассматриваются.) Для анализа влияния структуры поля поляризации будем считать, что грань гироскопа жестко закреплена, а остальные его грани - свободны от напряжений. Примем также, что полый цилиндр выполнен из пьезокерамики, причем рассмотрим два варианта его поляризации. В варианте 1 примем, что весь пьезокерамический материал цилиндра поляризован в радиальном направлении. Однако в этом случае для радиальной поляризации цилиндра нужен сплошной технологический электрод на всей внешней поверхности . Этот электрод после завершения процесса поляризации должен быть «стравлен» до электродов ( ), что приводит к некоторым дополнительным затратам. Технологически более просто сразу создать четыре электрода ( ) на внешней поверхности и затем их использовать для поляризации пьезокерамического материала цилиндра. Тогда материал будет поляризован в цилиндрической системе координат неоднородно, но можно приближенно принять, что в секторах под электродами ( ) пьезокерамика поляризована в радиальных направлениях, а вне этих секторов - не поляризована. Такой случай кусочно-однородной поляризации назовем вариантом 2 постановки задачи. 2. Математические и конечно-элементные модели Постановка задачи для гироскопа базируется на теории пьезоэлектричества и учете эффектов, вызванных вращением [16-19]. В задаче об установившихся колебаниях пьезоэлектрического тела неизвестными принимаются вектор-функция амплитуды перемещений и функция амплитуды электрического потенциала В неинерциальной (относительной) системе отчета, связанной с вращающимся телом , уравнения теории пьезоэлектричества в отсутствие массовых сил и объемных зарядов можно записать в виде , (1) , (2) (3) , (4) В (1)-(3) использованы следующие обозначения: - плотность; - массив компонент тензора напряжений ; - вектор электрической индукции; - массив компонент тензора деформаций ; - вектор напряженности электрического поля; - матрица упругих жесткостей, измеренных при постоянном электрическом поле, размерностью ; - матрица пьезомодулей, размерностью ; - матрица диэлектрических проницаемостей, измеренных при постоянных деформациях, размерностью ; - параметр демпфирования; - операция транспонирования для матриц и векторов. К системе уравнений (1)-(4) необходимо присовокупить краевые условия на границе , которые подразделяются на механические и электрические. При формулировке механических краевых условий примем, что торец , обозначаемый через , жестко закреплен, а остальная часть границы свободна от напряжений, т.е. , ; , , (5) где - вектор внешней единичной нормали к . Электрические граничные условия включают задание значений потенциала на электродах , и ; нулевые заряды на и ; а также свободную неэлектродированную поверхность ( ). Математически эти условия можно записать следующим образом: , ; , ; , , (6) , ; , , (7) , ; , (8) причем не зависящие от значения изначально неизвестны, а определяются из решения задачи. Для решения задач электроупругости (1)-(8) можно перейти к их слабым постановкам и использовать технику конечно-элементных аппроксимаций. Пусть - область, занимаемая соответствующей конечно-элементной сеткой, , , где - отдельный конечный элемент с номером . Для задач об установившихся колебаниях можно искать приближенное решение { , } в форме , (9) где - матрица базисных функций для перемещений; - вектор-строка базисных функций для электрического потенциала; , - векторы узловых амплитуд перемещений и электрического потенциала. В соответствии с обычной конечно-элементной техникой аппроксимируем континуальную слабую постановку задачи пьезоэлектричества в конечномерных пространствах, связанных с базисными функциями , . Подставляя (9) и аналогичные представления для проекционных функций в слабую постановку задачи пьезоэлектричества для , имеем следующую конечно-элементную систему [18]: , (10) , (11) где , , , (12) , , (13) , , , (14) , (15) Здесь - матрица масс; - матрица демпфирования; - матрица жесткости; - матрица пьезоэлектрической связности; - матрица диэлектрической проницаемости; - кориолисова матрица демпфирования; - матрица жесткости за счет центробежных сил; , - векторы обобщенных узловых сил и зарядов, получающиеся при учете неоднородных главных граничных условий. 3. Варианты постановки задачи с различными граничными условиями Для случая полностью радиально поляризованного пьезокерамического материала трубки рассмотрим еще различные варианты граничных условий на торцах и . Как и раньше, в варианте 1 грань считается жестко закрепленной, а грань - свободной от напряжений. В варианте 3 на каждый торец и жестко прикрепляются упругие стержни, направленные от границ электродов до оси цилиндра (рис. 2, а), а в варианте 4 - стержни, направленные от центральных точек электродов до точек на оси цилиндра (рис. 2, б). Все стержни в вариантах 3 и 4 моделируются как упругие брусы и учитывают задачи растяжения, кручения и пространственного изгиба. Для точек, лежащих на оси симметрии, т.е. в центральных точках пересечения упругих стержней, в вариантах 3 и 4 ставятся условия равенства нулю перемещений и углов поворота по всем трем осям декартовой системы координат. Очевидно, что все описанные выше варианты допускают равноправные изгибные движения системы в плоскостях и . Отметим, что принятые варианты закрепления стержней при нулевых углах поворота в точках на оси цилиндра обеспечивают дополнительную жесткость всей конструкции. Это допущение не дает стержням вращаться в точности, как при их классическом шарнирном опирании, но зато затрудняет случайные вращения гироскопа при изгибных колебаниях, что важно для практики. Тем не менее принятые здесь варианты закрепления 3 и 4 позволяют приближенно реализовать условия шарнирного опирания торцов трубки гироскопа, использовавшиеся ранее в [20, 24] при анализе аналогичного гироскопа в рамках балочной теории изгиба. а б Рис. 2. Варианты крепления стержней на торцах Fig. 2. Variants of the rod end fixation 4. Конечно-элементное моделирование Для решения задачи использовался КЭ пакет ANSYS версии 11.0 [28]. В полом цилиндре создавалась каноническая сетка из пьезоэлектрических конечных элементов (КЭ) SOLID226 (гексаэдры с 20 узлами). Их материальные свойства модифицировались для каждого КЭ SOLID226 так, чтобы векторы поляризации в центрах элементов были направлены по соответствующему радиальному направлению . В варианте 2 данная процедура применялась только для секторов под верхними электродами, а в остальных областях для КЭ SOLID226 принимались материальные свойства неполяризованной пьезокерамики, рассчитываемые по формулам гомогенизации материала кристаллографического класса 6mm до изотропного [29]: ; ; ; , а также пренебрежимо малые значения пьезомодулей. Входные данные, общие при всех расчетах: константы пьезокерамики PZT-5H - кг/м3, Н/м2, Н/м2, Н/м2, Н/м2, Н/м2, Кл/м2, Кл/м2, Кл/м2, , , Ф/м ( - диэлектрическая проницаемость вакуума); геометрические размеры - (мм); (мм); (мм); ; (В). В вариантах 3, 4 упругие стержни считались выполненными из стали плотностью кг/м3, модулем Юнга Н/м2 и с коэффициентом Пуассона Принималось, что поперечные сечения стержней являются квадратными со стороной (мм) для варианта 3 и со стороной (мм) для варианта 4. Стержни моделировались упругими КЭ BEAM4. Для полого цилиндра была построена каноническая КЭ гексаэдральная сетка с разбиением длины на 20 пьезоэлектрических элементов, толщины - на 2 элемента, окружного направления под электродами - на 4 элемента и окружного направления между электродами - на 8 элементов. Стержни по радиальным направлениям от точки на оси до внутренней поверхности трубки разбивались на 10 упругих балочных элементов. Данная густота сетки обеспечивала достаточную стабилизацию и точность КЭ расчетов. Вначале решалась задача модального анализа об определении первых частот электрических резонансов и антирезонансов , причем запоминались первые две частоты изгибных колебаний. Для этого задача (10), (11) относительно частоты , узловых перемещений и потенциалов при , , каждый раз решалась дважды с различными граничными условиями. Для поиска частот электрических резонансов принималось, что на , и , а электроды и считались свободными. Во второй задаче о поиске частот электрических антирезонансов принималось, что на , а все остальные электроды считались свободными. Подчеркнем, что при ненулевой скорости вращения кориолисова матрица также не равна нулю, и поэтому при использовании решателей ANSYS DAMP или QRDAMP [28, 30, 31] в задаче (10), (11) при , , собственные частоты оказывались комплексными. Кроме того, в силу симметрии устройства относительно плоскостей и , а также возможности равноправных изгибных движений в этих плоскостях, первые антирезонансные частоты изгибных колебаний будут кратными при . Если же , то кратные частоты будут трансформироваться в две однократные резонансные частоты, расположенные тем ближе друг к другу, чем меньше частота вращения . На следующем шаге решалась задача (10), (11) о вынужденных колебаниях с, вообще говоря, ненулевыми , , на нескольких интервалах частот, причем стыковочные границы частотных интервалов совпадали с вещественными частями первых частот электрических резонансов и антирезонансов, найденных из задачи модального анализа. Здесь при известных частотах решались задачи (10), (11) при частотно-независимом коэффициенте затухания с добротностью . 5. Обсуждение и результаты Результаты расчетов задачи (10), (11) при (Гц) для двух типов поляризации в виде зависимостей от частоты наведенного потенциала , точнее отношения , приведены на рис. 3, а. Кривые 1 и 2 получены для задач с соответствующими номерами вариантов типа поляризации. Как можно заметить, максимумы достигаются на резонансных частотах изгибных колебаний, причем вариант 1 с полной радиальной поляризацией дает наибольшие максимумы на первых двух частотах изгибных колебаний. Вариант 2 (с частичной поляризацией пьезокерамики под верхними электродами) также дает достаточно большие максимумы, лишь немного уступающие максимумам варианта 1. Для варианта 2 также несколько меньше коэффициенты электромеханических связей, чем для варианта 1. Таким образом, на основании проведенных вычислений можно сделать вывод о предпочтительности гироскопа с полной радиальной поляризацией, хотя и более простой в изготовлении гироскоп с частичной поляризацией под верхними электродами также достаточно эффективен. Рис. 3. Амплитудно-частотные характеристики выходного потенциала Fig. 3. The gain-frequency characteristics of the output electric potential На рис. 3, б приведены аналогичные результаты расчетов задачи (11), (12) при Гц для граничных условий по вариантам 3 и 4, имитирующих условия шарнирного опирания. Кривые 3, 4 получены для задач с соответствующими номерами вариантов граничных условий. Как видно из рис. 3, б, максимумы в вариантах 3 и 4 также достигаются на резонансных частотах изгибных колебаний, причем эти резонансные частоты выше соответствующих резонансных частот для трубки той же длины, но с одним жестко защемленным торцом. Между тем из сравнения рис. 3, а и б, можно заключить, что вариант 1 с жестким закреплением дает существенно большие максимумы сразу на двух соответствующих частотах изгибных колебаний. При этом оба максимума для варианта 1 приблизительно одинаковы, а для вариантов 3 и 4 пары максимумов достаточно существенно различаются. Основным преимуществом вариантов 3 и 4, имитирующих условия шарнирного опирания, по сравнению с вариантом 1 с одним закрепленным торцом является более широкая полоса рабочих частот, для которых значения выходного потенциала достаточны для регистрации частоты вращения. Отметим, что дополнительно были проведены также расчеты гироскопа с граничными условиями по вариантам 3 и 4, но с частичной поляризацией пьезокерамики под электродами на внешней боковой поверхности. Полученные результаты показали, что максимальные значения выходного потенциала на первых частотах изгибных колебаний трубки с частичной поляризацией могут как возрастать, так и быть меньшими по сравнению с приведенными на рис. 3, б. Но даже наибольший максимум отношения , равный 2,93 для варианта 4 с частичной поляризацией, оказался более чем в два раза меньше минимального максимума 5,87 для трубки с одним закрепленным концом. Таким образом, гироскоп с граничными условиями по вариантам 3 и 4, как при полной, так при частичной поляризации пьезокерамического материала трубки, оказался менее эффективным по амплитуде выходного потенциала. Во всех рассмотренных случаях наличие двух близко расположенных резонансных частот обусловлено симметрией конструкции гироскопа. Например, для варианта 1 при отсутствии вращения ( ) при модальном конечно-элементном анализе были найдены две первые частоты электрических резонансов Гц, Гц, и две кратные частоты электрических антирезонансов Гц. На этих частотах трубка изгибается в плоскостях и . В задаче для частот электрических резонансов электроды вдоль оси закорочены, а электроды вдоль оси - свободны. Данное различие в электрических граничных условиях определяет и небольшое отличие собственных частот и . Между тем в задаче для частот антирезонансов все электроды на внешней поверхности считаются свободными. Таким образом, для антирезонансных частот имеется полная симметрия задачи модального анализа относительно плоскостей и , что определяет наличие кратных изгибных мод колебаний. Появление вращения нарушает симметрию задачи, и при (Гц) частоты электрических резонансов и антирезонансов немного сдвигаются и уже не являются кратными: Гц, Гц, Гц, Гц. Именно на частотах и наблюдаются максимумы и минимумы выходного потенциала на свободных электродах, расположенных вдоль оси (кривая 1 на рис. 3, а). Наведенное силами Кориолиса за счет прямого пьезоэффекта электрическое поле вдоль оси в свою очередь, в силу обратного пьезоэффекта вызывает вторичные колебательные движения в плоскости . Эти движения могут даже превосходить по амплитуде первичные колебательные движения в плоскости . Так, на рис. 4 приведены в одном и том же масштабе формы амплитуд смещений для частот (рис. 4, а) и (см. рис. 4, б). Видно, что колебания на обеих резонансных частотах происходят преимущественно в плоскости вторичных движений , и в меньшей степени - в плоскости первичных движений . Этот феномен проявляется в большей степени для второй резонансной частоты (рис. 4, б). Доминирование вторичных движений над первичными позволяет рекомендовать рабочую частоту гироскопа выбрать вблизи второй резонансной частоты , поскольку тогда чувствительность конструкции к частоте вращения будет возрастать. а б Рис. 4. Формы колебаний для варианта 1 вблизи резонансных частот Fig. 4. Vibration modes for case 1 near resonance frequencies Подводя итог сравнения различных вариантов трубчатого пьезоэлектрического вибрационного гироскопа, можно сделать вывод о предпочтительности конструкции с одним жестким торцом, хотя аналогичные устройства со стержнями на торцах, имитирующими условия шарнирного опирания, также вполне работоспособны и имеют некоторые преимущества по ширине полосы частот. Заметим, что аналогичные выводы о преимуществах гироскопа с полной радиальной поляризацией и варианта с жестким закреплением были сделаны в [25, 26] при сравнении трубчатых гироскопов различной длины, но с приблизительно одинаковыми первыми резонансными частотами изгибных колебаний.

About the authors

A V Nasedkin

Southern Federal University

E I Shprayzer

Don State Technical University

References

  1. Apostolyuk V. Coriolis vibratory gyroscopes: Theory and design. - Cham, Switzerland: Springer International Publishing, 2016. - 117 p.
  2. A review of vibratory gyroscopes / J.S. Burdess, A.J. Harris, J. Cruickshank, D. Wood, G. Cooper // Eng. Sci. Educ. J. - 1994. - Vol. 3. - No. 6. - P. 249-254. doi: 10.1049/esej:19940603
  3. Singh A.K. Piezoelectric gyro sensor technology // Def. Sci. J. - 2007. - Vol. 57. - No. 1. - P. 95-103. doi: 10.14429/dsj.57.1735
  4. Yang J. A review of analyses related to vibrations of rotating piezoelectric bodies and gyroscopes // IEEE Trans. Ultrason., Ferroelect. and Freq. Control. - 2005. - Vol. 52. - No. 5. - P. 698-706. doi: 10.1109/TUFFC.2005.1503958
  5. Hunt G.H., Hobbs A.E.W. Development of an accurate tuning fork gyroscope // Proc. Inst. Mech. Eng., Conf. Proc. - 1964. - Vol. 179. - No. 5. - P. 129-139. doi: 10.1243/PIME_CONF_1964_179_116_02
  6. Fan M., Zhang L. Research progress of quartz tuning fork micromachined gyroscope // Proc. 2015 Int. Conf. Artificial Intelligence and Industrial Engineering, Atlantis Press, 2015. - P. 361-364. doi: 10.2991/aiie-15.2015.100
  7. Quartz rate gyro sensor for automotive control / Y. Nonomura, M. Fujiyoshi, Y. Omura [et al.] // Sensors and Actuators A. - 2004. - Vol. 110. - P. 136-141. doi: 10.1016/j.sna.2003.10.048
  8. Bhadbhade V., Jalili N., Mahmoodi S.N. A novel piezoelectrically actuated flexural/torsional vibrating beam gyroscope // J. Sound and Vibration. - 2008. - Vol. 311. - P. 1305-1324. doi: 10.1016/j.jsv.2007.10.017
  9. Seok J., Scarton H.A. Dynamic characteristics of a beam angular-rate sensor // Int. J. Mech. Sci. - 2006. - Vol. 48. - P. 11-20. doi: 10.1016/j.ijmecsci.2005.09.015
  10. Seok J., Tiersten H.F., Scarton H.A. An analysis of a vibratory angular-rate gyroscope using polarized piezoceramic bimorph plates. Part 2: Solution procedure for the gyroscope with superposed angular velocity // J. Sound and Vibration. - 2005. - Vol. 280. - P. 289-310. doi: 10.1016/j.jsv.2003.12.031
  11. Effect of polarized electric field on piezoelectric cylinder vibratory gyroscope / Y. Ting, J.-S. Huang, J.-L. Huang, C.-M. Yang // Sensors and Actuators. A. - 2006. - Vol. 128. - No. 2. - P. 248-256. doi: 10.1016/j.sna.2006.02.034
  12. Optimization and analysis of novel piezoelectric solid micro-gyroscope with high resistance to shock / Y. Lu, X. Wu, W. Zhang, W. Chen, F. Cui, W. Liu // Microsyst. Technol. - 2010. - Vol. 16. - No. 4. - P. 571-584. doi: 10.1007/s00542-009-0963-6
  13. Vibration analysis of a piezoelectric micromachined modal gyroscope (PMMG) / X.S. Wu, W.Y. Chen, Y.P. Lu, Q.J. Xiao, G.Y. Ma, W.P. Zhang, F. Cui // J. Micromech. Microeng. - 2009. - Vol. 19. - No. 12. - P. 125008. doi: 10.1088/0960-1317/19/12/125008
  14. Wu X., Chen W., Reza A. Optimal design of piezoelectric micromachined modal gyroscope (PMMG) with modes matched: modal analysis // Microsyst. Technol. - 2014. - Vol. 21. - No. 1. - P. 101-108. doi: 10.1007/s00542-014-2156-1
  15. Yang J.S, Fang H.Y., Jiang Q. A vibrating piezoelectric ceramic shell as a rotation sensor // Smart Mater. Struct. - 2000. - Vol. 9. - P. 445-451. doi: 10.1088/0964-1726/9/4/307
  16. Kagawa Y., Tsuchiya T., Kawashima T. Finite element simulation of vibrator gyroscopes // IEEE Trans. Ultrason., Ferroelectr., Freq. Control. - 1996. - Vol. 43. - P. 509-518. doi: 10.1109/58.503709
  17. Kagawa Y., Tsuchiya T., Sakai T. Three-dimensional finite element simulation of a piezoelectric vibrator under gyration // IEEE Trans. Ultrason., Ferroelectr., Freq. Control. - 2001. - Vol. 48. - P. 180-188. doi: 10.1109/58.895931
  18. Наседкин А.В. Моделирование пьезоэлектрических преобразователей в ANSYS: Учеб. пособие. - Ростов н/Д: Изд-во Юж. федер. ун-та, 2015. - 176 с.
  19. Nasedkin A.V. Some finite element methods and algorithms for solving acousto-piezoelectric problems // Piezoceramic materials and devices, ed. I.A. Parinov. - N.-Y.: Nova Science Publishers. - 2010. - P. 177-218.
  20. Yang J.S., Fang H.Y. A new ceramic tube piezoelectric gyroscope // Sensors and Actuators. A. - 2003. - Vol. 107. - No. 1. - P. 42-49. doi: 10.1016/S0924-4247(03)00209-7
  21. Левицкий А.А., Маринушкин П.С. Численное моделирование пьезоэлектрического вибрационного гироскопа // Датчики и системы. - 2009. - № 9. - С. 11-14.
  22. Левицкий А.А., Маринушкин П.С. Исследование режимов работы пьезоэлектрического вибрационного гироскопа осцилляторного типа // Датчики и системы. - 2011. - № 3. - С. 55-59.
  23. Левицкий А.А., Маринушкин П.С. Анализ погрешностей миниатюрного пьезоэлектрического гироскопа // Успехи современной радиоэлектроники. - 2012. - № 9. - С. 34-38.
  24. Marinushkin P.S. Dynamic analysis of tubular piezoelectric gyroscopes // XI International Conference and Seminar on Micro/Nanotechnologies and Electron Devices EDM'2010: Conference Proceedings. - Novosibirsk, 2010. - P. 168-171. doi: 10.1109/EDM.2010.5568838
  25. Наседкин А.В., Шпрайзер Е.И. Конечно-элементный анализ пьезокерамического вибрационного гироскопа в форме полого цилиндра при различных условиях закрепления // Инновационные технологии в науке и образовании. ИТНО-2015: сб. науч. тр. междунар. науч.-метод. конф. Ростов-на-Дону - пос. Дивноморское, 7-10 сентября, 2015. - Ростов-на-Дону - Зерноград, 2015. - С. 203-208.
  26. Наседкин А.В., Шпрайзер Е.И. Анализ вибрационного гироскопа в форме полого цилиндра из неоднородно поляризованной пьезокерамики // Фундаментальные проблемы радиоэлектронного приборостроения. - 2015. - Т. 15, № 4. - С. 57-60.
  27. Юрлов М.А., Юрлова Н.А. Стратегии демпфирования колебаний конструкций с пьезоэлементами и внешними электрическими цепями и их экспериментальная иллюстрация // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2014. - № 4. - С. 240-270. doi: 10.15593/perm.mech/2014.4.10
  28. ANSYS. Theory reference for ANSYS and ANSYS Workbench. Rel.11. Ed. P. Kothnke. - Canonsburg, ANSYS Inc., 2007. - 1072 p.
  29. Шермергор Т.Д. Теория упругости микронеоднородных сред. - М.: Наука, 1977. - 400 с.
  30. Numerical recipes: The art of scientific computing / W.H. Press, S.A. Teukolsky, W.T. Vetterling, B.P. Flannery. - New York, Cambridge University Press, 2007. - 1256 p.
  31. Peters G., Wilkinson J.H. The calculation of specified eigenvectors by inverse iteration // Handbook for Automatic Computation. Vol. 2. Linear Algebra, еds. J.H. Wilkinson, C. Reinsch. - New York, Springer-Verlag, 1971. - P. 418-439. doi: 10.1007/978-3-642-86940-2_29

Statistics

Views

Abstract - 89

PDF (Russian) - 84

Cited-By


PlumX


Copyright (c) 2016 Nasedkin A.V., Shprayzer E.I.

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies