2016
Математика и механика
№ 6(44)
МЕХАНИКА
УДК 629.76
DOI 10.17223/19988621/44/4
А.В. Азин, С.В. Пономарев, С.В. Рикконен, А.М. Храмцов
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ РЕЖИМОВ РАБОТЫ ПЬЕЗОДВИГАТЕЛЯ1
Разработана одномерная математическая модель пьезодвигателя с учетом механической и акустической передачи энергии. Модель позволяет рассчитать частотные характеристики механоакустической колебательной системы, по которым выбираются материалы элементов и конструкция двигателя, позволяющие согласовать источник энергии с нагрузкой изделия.
Ключевые слова: пьезодвигатель, космический аппарат, колебательная система, динамическая сила, вибросмещение.
В настоящее время миниатюрные двигатели и приводы широко применяются в технике. Практически во всех отраслях промышленности миниатюрные приводы востребованы, а в космической и военно-промышленной отрасли они просто необходимы. Постоянное требование улучшения массогабаритных показателей вынуждает исследователей искать новые принципы преобразования энергии из электрической в механическую, совершенствовать технику и технологию.
При расчете и проектировании пьезодвигателей и пьезоприводов для высокого коэффициента передачи энергии необходимо учитывать условие согласования пьезоактюатора (ПА) (источника механической энергии) с нагрузкой. При неправильном выборе конструкции и ошибочном выборе материалов элементов конструкции возможен вариант полного демпфирования передачи энергии от ПА к нагрузке. В этом случае особо актуальным является математическое моделирование конструкции изделия с правильным учетом конкретной нагрузки для применения в соответствующей отрасли промышленности.
Целью работы является разработка математической модели для расчета режимов работы ПА (согласование излучателя с нагрузкой) с максимальным коэффициентом передачи энергии. Особенность конструирования пьезодвигателей заключается в том, что частотный диапазон возбуждения ПА лежит в области частот, при которых необходимо учитывать совместную передачу энергии механическим и акустическим способами. На рис. 1 представлено схематическое изображение пьезодвигателя, работающего в режиме «короткого замыкания» (КЗ), то есть без механической нагрузки.
ПА с синусоидальной силой F = Fm cos(œt) посредством излучающей плиты массой Мр1 через пружину предварительного пожатия K воздействует на толкатель
1 Работа выполнена при финансовой поддержке Минобрнауки России, уникальный идентификатор ЯРМЕР!57814Х0060.
массой та. Подпружиненная масса толкателя расположена на излучателе, механическая система имеет коэффициент вязкого трения Я. Излучатель является абсолютно жестким, синусоидальная сила равномерно распределена по площади излучения. Это линейная одномерная колебательная система с одной степенью свободы со следующими допущениями:
- колебательная система представлена чисто линейной одномерной системой -все параметры постоянны и независимы от амплитуды колебаний, а также от частоты возбуждения;
- сила возбуждения чисто синусоидальная, излучающая плита создает плоскопараллельные волны в толкателе;
- излучающая плита при любых режимах работы сохраняет контакт с толкателем;
- излучатель является абсолютно жестким, синусоидальная сила равномерно распределена по площади плиты излучения;
- длина толкателя значительно меньше длины волны излучения.
Движение такой чисто механической колебательной системы описывается следующим дифференциальным уравнением:
dX,
. dX,
М-1- + К —1 + RX1 = Fm cos(rat).
dt
dt
(1)
Амплитуда виброскорости (чисто механической системы) массы толкателя определяется по формуле
X =
Fm
R2 + (юМ - К )2 ю
(2)
где масса М включает в себя массу толкателя - та, массу излучающей плиты - Мр\ и 30 % массы ПА - МасЬ К - жесткость пружины предварительного поджатия.
В [1] приведены два случая, описывающие параметры и поведение системы, когда (кг) >> 1 и (кг) << 1, где г - радиус излучателя, м; к = 2п/ X - волновое число; X = с// - длина волны, м.
Рис. 1. Схематическое изображение пьезодви-гателя: 1 - корпус пьезодвигателя; 2 - толкатель; 3 - пружина предварительного поджатия; 4 - излучающая плита; 5 - пьезоактюатор; 6 -винт поджатия
Fig. 1. Schematic representation of a piezo motor: 1, case of the piezo motor; 2, pusher; 3, spring of initial load; 4, emissive plate; 5, PZT stack actuator; 6, preload screw
//////////////////
2
5
1
При низких частотах воздействия длина вибрационной волны много больше размеров излучателя и произведение (кг) << 1. Если диаметр излучателя меньше одной трети длины волны (или (кг) << 1), то дальнее поле излучения можно не учитывать [1]. Частотный диапазон и размеры излучателя в системах пьезоприво-дов имеют произведение близкое к (кг) < 1.5, которое предопределяет учет упругой (акустической) передачи энергии вдоль толкателя. Для этого случая можно применить методику определения параметров колебательной системы, предложенную И. С. Чичининым [2]. Акустическое сопротивление (сопротивление излучения) наиболее просто трактуется, когда все точки излучающей поверхности источника колеблются синхронно, с одинаковой амплитудой (рис. 2). Сопротивление излучения абсолютно жесткой плиты в этом случае определяется отношением вибрационной силы на виброскорость 2п = Е„ / X. Достаточно найти комплексную амплитуду X скорости перемещения поверхности 5". На рис. 2 представлена акустическая подсистема механоакустической колебательной системы, ее механическая и электрическая схемы замещения. На рис. 2 приняты следующие обозначения: 1 - излучающая плита; 2 - присоединенная масса колебательной акустической подсистемы; 3 - масса излучающей платформы в механической схеме замещения - Мр1; 4 - присоединенная масса колебательной акустической подсистемы в механической схеме замещения - пп; 5 - жесткость колебательной акустической подсистемы - Кп; 6 - неупругое сопротивление колебательной акустической подсистемы - Лп; 7 - инерционное сопротивление штампа в электрической схеме замещения - (/шМр1); 8 - неупругое сопротивление акустической подсистемы в электрической схеме замещения - Яп; 9 - упругое сопротивление акустической подсистемы в электрической схеме замещения - (Кп/(/'ю)); 10 - инерционное сопротивление присоединенной массы в электрической схеме замещения - (/юпп).
10
2
Рис. 2. Схема излучения жесткой плиты, расположенной на упругом полупространстве при (kr) < 1.5. Представлена акустическая подсистема механоакустической колебательной системы: а - абсолютно жесткая плита на упругом полупространстве; b - механическая схема замещения колебательной системы; c - электрическая схема замещения излучения волн Fig. 2. Scheme of a rigid emissive plate located on an elastic half-space at (kr) < 1.5. Acoustic subsystem of a mechanical acoustic oscillating system: (a) completely rigid plate on an elastic half-space; (b) mechanical equivalent scheme of the oscillating system; (c) electrical equivalent circuit of the emission of waves
8
3
9
c
Движение такой акустической колебательной системы описывается следующим дифференциальным уравнением (3):
dX2 йХ2
(^р1 + П)-2- + Кп—2 + ад = Рп ес8(Ю/). (3)
ж1 ж
Излучение происходит за счет деформации упругого полупространства под излучателем, процесс излучения зависит от присоединенной массы вещества объ-
екта тп, которая колеблется в фазе с излучателем, акустическая жесткость материала Кп и коэффициент диссипации Яп [3-5]. В процессе излучения участвует масса излучающей плиты - Мр\.
Амплитуда виброскорости жесткой излучающей плиты определяется по формуле
X2 = ■ ¥т 2 . (4)
^ + [а(Мр1 + тп) - К- )
Параметры упругого излучения вычисляются следующим образом [2].
Коэффициент демпфирования (диссипация)
Яп = 7(1 - у2) р У Я2 аг, (5)
Акустическая жесткость материала определяется по формуле
Кп = 8(1 - у2) р У/ Я, , (6)
Присоединенная масса вычисляется следующим образом:
тп = (1 - Т2) Р Я3 от , (7)
где - площадь излучателя, Яi - радиус излучателя, у = V/ Ур, Ув - скорость поперечных звуковых волн, Ур - скорость продольных звуковых волн, р - плотность материала толкателя с массой та, а - поправочные коэффициенты, близкие к единице, но зависящие от частоты возмущения и параметра у.
Колебательную систему пьезодвигателя или пьезопривода можно представить в виде механоакустической системы. Механоакустическая система имеет две колебательные подсистемы, которые имеют свои параметры, частотные характеристики и резонансные частоты. Виброскорости и Х2 как бы «нельзя» складывать, потому что эти колебания имеют разную «природу» и описываются разными математическими моделями. Виброскорость Х2 - это колебание упругого излучения, результатом колебательного процесса является деформация торца массы толкателя та. Виброскорость Х1 - это движение всей массы подвижных частей
конструкции как единого целого объекта. Однако обе подсистемы переносят энергию и участвуют в создании суммарном смещении торца толкателя Х2 = Х1 + Х2.
Аналогии между механическими колебательными системами и электрическими колебательными контурами позволили развить систему электрических аналогий механических систем. Механические аналогии электрических величин приведены в таблице.
Электрическая система Первая система механическая Вторая система механическая Первая система акустическая
Напряжение и Сила ^ Скорость У Звуковое давление Р
Сила тока I Скорость У Сила ^ Объемная скорость (5У)
Индуктивность Ь Масса т Гибкость с Акустическая масса та=(р1/5)
Активное сопротивление Я Сопротивление г Величина, обратная коэффициенту сопротивления Сопротивление акустических потерь га
Емкость С Гибкость с Масса т Акустическая податливость Сс = У/(рС2)
Используя правила преобразования электрических цепей, можно относительно быстро произвести расчеты режимов работы колебательных систем [3-6].
Для расчета механоакустической колебательной системы пьезодвигателя, представленной на рис. 1, должны совместно использоваться дифференциальные уравнения (1) и (3), а совместному решению этих уравнений может соответствовать электрическая схема замещения на рис. 3 [7, 8].
10
1
Т
1
Т
1
Т
Рис. 3. Электрическая схема замещения колебательной системы пьезодвигателя, представленной на рис. 1 Fig. 3. Electrical equivalent circuit of the piezo drive oscillating system presented in Fig.1
На схеме замещения показана масса излучающей плиты Мр1 (1 - инерционное сопротивление излучающей плиты в электрической схеме замещения - ^Мр1), последовательно к ней присоединена механоакустическая колебательная система, с параллельным соединением элементов механической подсистемы 2 - 3 - 4 (2 -инерционное сопротивления толкателя (jrnma); 3 - упругое сопротивление пружины предварительного поджатия толкателя (K/jro)); 4 - неупругое сопротивление механической подсистемы (R)) и акустической подсистемы 5 - 6 - 7 (5 - инерционное сопротивления присоединенной массы акустической подсиситемы (jrom„); 6 - упругое сопротивление акустической подсистемы (Kn/(/ro); 7 - неупругое сопротивление акустической подсистемы (Rn)).Показаны варианты нагрузки колебательной системы - режим холостого хода 8 (Zn = ю), режим короткого замыкания 9 (Zn = 0), инерционная нагрузка без активных потерь 10 (Zn = ^МпаёГ).
Преобразуем электрическую схему замещения колебательной системы пьезо-двигателя, работающей на инерционную нагрузку.
В механоакустической системе чисто механическое сопротивление в символической форме будет выглядеть как
Z = R + j(®ma - -). (8)
ш
Чисто акустическое сопротивление в символической форме
Z2 = Rn + j(шmn - Ш). (9)
ш
Суммарное сопротивление системы будет выглядеть следующим образом:
2Ъ = j&M р1
12
z, + z0
■ + j®M n
(10)
Суммарное вибросмещение подвижной части конструкции определяется по формуле
Р 1
ХЕ= —. (11)
Сила на нагрузке равна
Рп = Р - X]юЫр1 + ■ <12)
Пример расчета
Материал толкателя - оргстекло. Эквивалентный радиус излучающей плиты ^ = 0.4 мм. Жесткость пружины К = 8.57-109 Н/м. Масса излучающей плиты Мр1 = 0.10110-3 кг. Масса толкателя та = 110-3 кг. Коэффициент потерь на трение Я = 525 кг/с. Присоединенная масса толкателя тп = 1.310-7 кг. Жесткость толкателя Кп = 3.63 108 Н/м. Неупругое сопротивление толкателя Яп = 146 кг/с. Масса пригруза МпабГ = 0.18 кг. Сила пьезоактюатора Р = 1082 Н.
При данной конструкции пьезопривода и материале толкателя (оргстекло) система имеет один явный резонанс на частоте 7000 Гц, длина волны равна 0.35 метра (значительно больше длины толкателя), основная часть энергии передается в нагрузку посредством акустической подсистемы - 95 %, 5 % - механической подсистемой (рис. 4).
f кГц
Рис. 4. Частотные характеристики вибросмещений механоакустической колебательной системы, рассчитанные по формуле (11): кр. 1 - вибросмещение механической подсистемы; кр. 2 - вибросмещение акустической подсистемы; кр. 3 - суммарное вибросмещение толкателя
Fig. 4. Frequency-response characteristics of vibration displacements of a mechanical acoustic oscillating system according to formula (11): 1, vibration displacement of the mechanical subsystem; 2, vibration displacement of the acoustic subsystem; 3, total vibration displacement of the pusher
Колебательная система имеет незначительные активные сопротивления, на которых рассеивается энергия и поэтому система может на резонансе существенно раскачаться. Это подтверждается высокими значениями силы на нагрузке (рис. 5).
Рис. 5. Частотная характеристика силы на нагрузке механоакустической системы Fig. 5. Frequency-response characteristic of the force on the load of mechanical acoustic oscillating system
Мгновенная мощность на инерционной нагрузке без потерь изменяется во времени с двойной частотой возбуждения, амплитуда данной мощности равна
Р^п = РптХТ 5Ь(2Ю'), (13)
где Гпт - амплитудное значение силы на инерционной нагрузке, Хт - амплитудное значение виброскорости в нагрузке.
В режиме работы на инерционную нагрузку мгновенная мощность за период колебания дважды меняет знак, временные кривые мощности симметричны относительно оси абсцисс, следовательно, полезной работы система не производит.
Заключение
По результатам проведенного исследования предложен подход к моделированию пьезодвигателя, позволяющий рассчитать по отдельности передачу механической энергии через механическую и акустическую подсистемы с дальнейшим суммированием потоков энергии на нагрузке.
Для решения линейных одномерных задач согласования источника энергии с нагрузкой целесообразно применять математическое моделирование, основанное на аналоговых электрических схемах замещения первого рода.
Механическая система согласования источника энергии с нагрузкой для частот работы пьезодвигателей представляется механоакустической колебательной системой, в которой передача энергии осуществляется через два канала: механической и акустической подсистемами.
В аналоговой электрической схеме замещения колебательной системы механо-акустическая система представляет собой параллельный электрический контур, сопротивление которого зависит от частоты колебаний и может иметь величину от 0 до бесконечности (энергетическая пробка).
Конструкция элементов и материалы, из которых сделаны элементы двигателей, существенно влияют на механическое сопротивление механоакустической системы и, следовательно, на коэффициент преобразования энергии.
Данный математический подход позволяет конструктору в приближенной форме выбрать материал элементов, конструкцию пьезодвигателя в зависимости от того, по какой механоакустической подсистеме намерены передавать энергию (по акустической или механической), с какими виброперемещениями и силами на нагрузке, и окончательно определить частоту работы пьезодвигателя.
ЛИТЕРАТУРА
1. Скучик Е. Основы акустики. М.: Мир, 1976. 542 с.
2. Чичинин И.С. Вибрационное излучение сейсмических волн. М.: Недра, 1984. 221 с.
3. Дьелесан Э., Руайе Д. Упругие волны в твердых телах. М.: Наука, 1982. 424 с.
4. Афонин С.М. Многомерная структурно-параметрическая модель составного пьезодвига-теля наноперемещений // Вестник машиностроения. 2007. № 1. С. 3-13.
5. Богданов Е.П., Рикконен С.В. Экспресс-метод определение параметров нагрузок вибрационных систем // Интернет-журнал «Науковедение». 2013. № 5.
6. Амеличев В.В., Вернер В.Д., Ильков А.В. МЭМС-микрофон. Выбор материалов, конструкции и технологии. Часть I. Электромеханический элемент // Нано-микросистемная техника. 2007. № 2. С. 53-62.
7. Пономарев С.В., Рикконен С.В, Азин А.В. Анализ работы пьезопривода по частотным характеристикам механоакустической колебательной системы // Изв. вузов. Физика. 2014. Т. 57. № 8/2. С. 196-202.
8. Пономарев С.В., Рикконен С.В, Азин А.В. Моделирование колебательных процессов пьезоэлектрического преобразователя // Вестник Томского государственного университета. Математика и механика. 2015. № 2(34). С. 86-95.
Статья поступила 07.11.2016
Azin A.V., Ponomarev S.V., Rikkonen S.V., Khramtsov A.M. MATHEMATICAL MODELING OF PIEZO MOTOR OPERATION MODES. Tomsk State University Journal of Mathematics and Mechanics. 6(44). pp. 45-53
DOI 10.17223/19988621/44/4
At present, micro electromechanical systems appear to be progressively applicable in engineering. The problem of reducing weight-and-dimensional characteristics of the spacecraft (SC) system is of great importance in the space field. One of the solutions is replacement of electromechanical motors in different SC units by piezo motors. It is necessary to take into account the condition of consistency of PZT stack actuators (the sources of mechanical energy) with the load when calculating and designing the piezo motors of high energy transfer coefficient. In this case, the most relevant is mathematical modeling of a product design with correct consideration of the specific drive load in an actual industrial field.
Mathematical modeling by the finite element method in combination with 3D numerical simulation is most accurate and appropriate. However, this method is laborious and time-consuming. In practice, designers and constructors prefer to use simple linear one-dimensional mathematical models, which provide an accurate information about materials and product design. In linear one-dimensional models of oscillatory systems for the range of frequencies applied in piezo motors, it is necessary to consider the acoustic and mechanical principles of energy conversion and transfer from the source to the load. In this paper, a one-dimensional mathematical model of the piezo motor with consideration of mechanical and acoustic energy transfer is developed. This model allows calculating the frequency-response characteristics of mechanical acoustic oscillating system according to which the materials of elements and design of the drive are chosen.
Keywords: piezo motor, spacecraft, oscillating system, dynamic force, vibration displacement.
AZIN Anton Vladimirovich (Candidate of Physics and Mathematics, Tomsk State University, Tomsk, Russian Federation) E-mail: antonazin@niipmm.tsu.ru
PONOMAREV Sergey Vasil'evich (Doctor of Physics and Mathematics, Tomsk State University, Tomsk, Russian Federation) E-mail: psv@niipmm.tsu.ru
RIKKONEN Sergey Vladimirovich (Candidate of Technical Sciences, Tomsk State University,
Tomsk, Russian Federation)
E-mail: rikk2@yandex.ru
KHRAMTSOVAlexeyMikhaylovich (Tomsk State University, Tomsk, Russian Federation)
E-mail: khramtsov.home@gmail.com
REFERENCES
1. Skuchik E. (1976) Osnovy akustiki [Foundations of Acoustics]. Moscow: Mir.
2. Chichinin I.S. (1984) Vibratsionnoe izluchenie seysmicheskikh voln [Vibrational Emission of Seismic Waves]. Moscow: Nedra.
3. Dieulesaint E., Royer D. (2000) Elastic Waves in Solids. Berlin; Heidelberg: Springer-Verlag.
4. Afonin S.M. (2007) Mnogomernaya strukturno-parametricheskaya model' sostavnogo p'ezodvigatelya nanoperemeshcheniy [Multidimensional structural parametric model of a complex piezo motor of nanodisplacements]. Russian Engineering Research. 1. pp. 3-13.
5. Bogdanov E.P., Rikkonen S.V. (2013) Ekspress-metod opredelenie parametrov nagruzok vibratsionnykh sistem [A rapid metod of determining the parameters of loads of vibration systems]. JournalNaukovedenie. 5(18).
6. Amelichev V.V., Verner V.D., Il'kov A.V. (2007) MEMS-mikrofon. Vybor materialov, konstruktsii i tekhnologii. Chast' I. Elektromehanicheskiy element [MEMS-microphone. Choice of materials, designs, and technologies. Part I. Electromechanical sensitive element]. Journal of Nano and Microsystem Technique. 2. pp. 53-62.
7. Ponomarev S.V., Rikkonen S.V, Azin A.V. (2014) Analiz raboty p'ezoprivoda po chastotnym kharakteristikam mekhanoakusticheskoy kolebatel'noy sistemy [Analysis of a piezodrive operation in accordance to the frequency-response characteristics of mechanical acoustic oscillating system]. Russian Physics Journal. 57 (8/2). pp. 196-202.
8. Ponomarev S.V., Rikkonen S.V, Azin A.V. (2015) Modelirovanie kolebatel'nykh protsessov p'ezoelektricheskogo preobrazovatelya [Simulation of oscillatory processes in a piezoelectric transducer]. Tomsk State University Journal of Mathematics and Mechanics. 2(34). pp. 86-95.