Dynamic axisymmetric problem of a direct piezoeffect for a round bimorph plate

Abstract


The dynamic axisymmetric problem for a round bimorph structure consisting of the metal substrate and the axially polarized piezoceramic plate is considered. Its bending modes are caused by the action of the mechanical loading (normal tension) on the side surface; the loading is a function of the radial coordinate and time. The rigid and hinge supports of the plate are considered. The initial design relations are formulated for the piezoceramic material with the hexagonal crystal lattice of the 6 mm class. In order to solve the problem of the theory of electrodynamics in a three-dimensional formulation, the finite integral Hankel transformations along the axial coordinate and a generalized transformation (FIT) over the radial variable is used. At each stage, the standardization is carried out which allows implementing an appropriate transformation algorithm. In the first case the boundary conditions are presented in a mixed form; and in the second case they are homogeneous by introducing auxiliary functions. This approach allows gaining precise (within the used models) calculated ratios in a most general form. The built closed solution allows defining the frequency range of the axisymmetric oscillations, the stress-strain state and the nature of the changing induced electric field of the bimorph plate. This makes it possible to establish the conventional solutions of the designed devices, determine a way of fixing the electrical signal, pick up all the geometrical and physical characteristics of the typical elements of the piezoceramic transducers Also, the developed solution allows solving the problems of the elasticity and electroelasticity theory for circular thick and thin plates with an arbitrary number of layers under most general loading conditions without the use of kinematic hypotheses.

Full Text

Введение В различных технических устройствах используются пьезокерамические датчики в виде тонких биморфных пластин [1-8]. Как правило, они состоят из двух жестко соединенных круглых пьезокерамических пластин с противоположной или параллельной направленностью вектора поляризации. Для повышения механической прочности в рассматриваемой конструкции применяется металлическая подложка [9, 10]. Расчет пьезокерамических многослойных конструкций в основном выполняется с помощью прикладных теорий для тонких пластин [11-13], в которых кинематические гипотезы дополняются аналогичными допущениями о характере изменения электрического поля [14-21]. Для более полного учета связанности физических полей в многослойных пьезокерамических пластинах в работах [22, 23] было проведено исследование и получено замкнутое решение нестационарной задачи обратного пьезоэффекта в трехмерной постановке. На основании численных результатов сделан вывод, что при аксиальной поляризации материала напряженность электрического поля по высоте тонкой пластины изменяется по линейному закону. Кроме того, установлено, что при использовании сплошного электродного покрытия касательными напряжениями, возникающими в конструкции, можно пренебречь. В настоящей работе рассматривается нестационарная задача прямого пьезоэффекта для асимметричной (отсутствует симметрия физических свойств материала относительно нейтральной плоскости) сплошной биморфной пластины с жестким, а также шарнирным закреплением ее цилиндрической поверхности. 1. Постановка задачи Пусть круглая двухслойная пластина, занимающая в цилиндрической системе координат область : , состоит из пьезокерамического элемента высотой , выполненного из аксиально поляризованного (электроупругого) материала с гексагональной кристаллической решеткой класса 6mm , и металлической (упругой) заземленной подложки толщиной ( ). Изгибные осесимметричные колебания возбуждаются за счет действия на лицевой поверхности конструкции ( ) механической динамической нагрузки (нормальных напряжений) , являющейся произвольной функцией радиальной координаты и времени (рис. 1). Подключение электродов к измерительному прибору позволяет зафиксировать величину и форму электрического напряжения . Рассматриваются случаи жесткого, а также шарнирного закрепления цилиндрической поверхности конструкции. Рис. 1. Биморфная пластина Fig. 1. Bimorph plate Дифференциальные уравнения движения и электростатики в цилиндрической системе координат и безразмерной форме имеют следующий вид [24, 25]*: , (1.1) , . Граничные условия на цилиндрической поверхности определяются равенствами , , , (1.2) жесткое защемление , (1.3) шарнирное закрепление , (1.4) Первые три неравенства (1.2) являются условиями регулярности решения для сплошной пластины, а последнее учитывает отсутствие электродного покрытия. Граничные условия по аксиальной координате формулируются следующим образом: , (1.5) (1.6) , , (1.7) Равенства (1.6) являются условиями совместности напряжений, деформаций и заземления подложки. Кроме того, последнее равенство (1.7) учитывает подключение электродированной поверхности пьезокерамической пластины к измерительному прибору с большим входным сопротивлением, что соответствует режиму «холостого хода» (отсутствию свободных электрических зарядов). Начальные условия задаются в общем виде и позволяют определить состояние пластины в момент приложения нагрузки: , (1.8) Перевод размерных функций и переменных в безразмерные величины производится по формулам , , , . В равенствах (1.1)-(1.8) приняты следующие обозначения: - компоненты вектора перемещений и тензора напряжений ; - компоненты вектора индукции и потенциал электрического поля; , пьезомодули и коэффициенты диэлектрической проницаемости электроупругого материала ; объемная плотность, модули упругости электроупругого и упругого материалов; известные в начальный момент времени перемещения, скорости перемещений; , . Соотношения (1.1)-(1.8) представляют математическую формулировку рассматриваемой краевой задачи электроупругости. 2. Построение общего решения На первом этапе решения краевой задачи (1.1)-(1.8) используется метод конечных интегральных преобразований Ханкеля по радиальной координате [26]. При этом данное преобразование позволяет удовлетворить только смешанные граничные условия. Для выполнения данного требования необходимо: - для жесткого закрепления первое равенство (1.3) заменить на условие отсутствия на цилиндрических поверхностях пластины ( ) касательных напряжений: (2.1) - для шарнирного закрепления последнее соотношение (1.2) заменить условием наличия потенциала электрического поля на цилиндрической поверхности при . (2.2) Кроме того, вводятся новые функции , , , связанные с , , соотношениями , (2.3) , , . Здесь единичная функция Хэвисайда [27]; - неизвестные функции, определяемые в процессе решения задачи соответственно из условий отсутствия вертикальных перемещений цилиндрической поверхности ( ) пластины при (первое равенство (1.3)) и радиальной компоненты вектора индукции электрического поля при (последнее равенство (1.2)). В результате подстановки (2.3) в (1.1)-(1.8), (2.1), (2.2) получаем новую краевую задачу относительно функций . При этом дифференциальные уравнения (1.1), первое условие (1.5), первые два условия (1.6) и равенства (1.7) становятся неоднородными с правыми частями , , а начальные условия следует заменить на : , , , , а краевые условия (1.2)-(1.4) принимают следующий вид: - жесткое закрепление , , (2.4) - шарнирное закрепление , (2.5) Здесь следует отметить, что в отличие от классической постановки задач теории упругости в настоящем исследовании «жесткое» закрепление (2.4) характеризуется также отсутствием угла поворота. Кроме того, последнее условие (2.5) получается при . Только в этом случае, в дальнейшем, дополнительные внеинтегральные члены в первом трансформированном уравнении (1.1) равны нулю. Применяем к краевой задаче (1.1)-(1.8), (2.1), (2.2), (2.4), (2.5) относительно функций конечные интегральные преобразования Ханкеля, используя следующие трансформанты: (2.6) и формулы обращения (2.7) где положительные нули функций , соответственно при жестком и шарнирном закреплении пластины, расположенные в порядке их возрастания ; при жестком и при шарнирном закреплении. В пространстве изображений получаем новую краевую задачу: (2.8) ; (2.9) (2.10) (2.11) (2.12) где На втором этапе решения используется обобщенный метод конечных интегральных преобразований (КИП) [28] по координате . Предварительно выполняется процедура стандартизации, связанная с приведением граничных условий (2.9)-(2.11) к однородным, с помощью следующих разложений: , , (2.13) . В случае жесткого закрепления, принимая во внимание , , , имеем а при шарнирном закреплении, с учетом получаем , В результате подстановки (2.13) в (2.8)-(2.12), при учете соотношений: при жестком закреплении (2.14) : при шарнирном закреплении (2.15) получаем новую краевую задачу относительно функций , , с однородными граничными условиями по координате . Правые части дифференциальных уравнений (2.8) и начальные условия , заменяются на , , , , , , , , . Функции , входящие в разложения (2.13), определяются с помощью следующих равенств: , (2.16) где постоянные; - параметр, соответствующий количеству граничных условий для функций Подстановка (2.16) в (2.14),(2.15) позволяет определить функции - жесткое закрепление , - шарнирное закрепление , , , . Начально-краевую задачу (2.8)-(2.12) относительно функций решаем, используя обобщенный метод конечных интегральных преобразований (КИП) [28] с использованием неизвестных трансформанты и компонент вектор-функции ядра преобразования - положительные параметры, образующие счетное множество . Общее решение краевой задачи (2.8)-(2.12) для электроупругого и упругого слоев представлено в работах автора [22, 29, 30]. Здесь следует отметить, что в отличие от стандартной процедуры разложения по собственным функциям [26] метод КИП [28] позволяет в процессе решения динамической задачи определить частоты и формы собственных колебаний конструкции. Окончательные выражения функций , получим, применяя последовательно формулы обращения КИП [28] и метод конечных преобразований Ханкеля (2.7). В результате с учетом (2.3), (2.13) имеем , , (2.17) . Заключительным этапом исследования является определение неизвестных функций , . В случае жесткого закрепления пластины, когда функция определяется из условия отсутствия вертикальных перемещений цилиндрической поверхности пластины при и на основании (2.17) принимает вид (2.18) При шарнирном закреплении . Для упрощения расчета функция определяется из условия равенства нулю суммарного значения тока смещения на цилиндрической поверхности пьезокерамической пластины . Для этого представляем ее в виде следующего многочлена: . (2.19) Данное выражение учитывает условия неразрывности деформаций и заземление металлической подложки Постоянная определяется при удовлетворение условия . (2.20) 3. Численные результаты. Выводы В качестве примера рассматривается биморфная пластина, имеющая следующие физические и геометрические характеристики аксиально поляризованных пьезокерамических пластин состава ЦТС-19 [24] и металлической стальной подложки: H/м2, Кл/м2, Ф/м, кг/м3, H/м2, кг/м3, м, м. Рассмотрим случай действия равномерно-распределенной гармонической нагрузки , где - амплитудное значение интенсивности в безразмерной форме; - относительная частота вынужденных колебаний ( - круговая частота вынужденных колебаний). На рис. 2-4 изображены графики, характеризующие изменение напряженно-деформированного состояния и характер распределения электрического поля в рассматриваемой конструкции при ( - положительный параметр, соответствующий первой частоте собственных колебаний). Их анализ позволяет сделать следующие выводы. 1) При жестком закреплении конструкции на лицевой поверхности пьезокерамической пластины одновременно образуются зоны растяжения и сжатия, что приводит к появлению электрических зарядов разных знаков. Поэтому использование сплошного металлического покрытия нецелесообразно в связи с индукцией на электроде электрического потенциала небольшой величины. Для решения данной проблемы, когда частота внешнего воздействия , для регистрации разности потенциалов необходимо использовать два разрезных круговых электрода (радиус раздела электродов определяет нулевое значение функции ). Подключение электродов к измерительному прибору позволяет определить : Рис. 2. Зависимость амплитудных значений от толщины пьезокерамической пластины : а - шарнирное закрепление; б - жесткое защемление Fig. 2. The dependence of the amplitude values of the thickness of the piezoceramic plate: a - securing hinge; b - hard pinching Рис. 3. Изменение амплитудных значений по радиальной координате (1- , 2- ): а - шарнирное закрепление; б - жесткое защемление Fig. 3. The change in the peak values of the radial coordinate (1- , 2- ): a - securing hinge; b - hard pinching Рис.4. Изменение амплитудных значений по высоте пьезокерамической пластины (шарнирное закрепление) Fig. 4. The change in amplitude values, the height of the piezoceramic plate (hinge fixing) При шарнирном закреплении пластины разность потенциалов между электродированными поверхностями пьезокерамических пластин (в данном случае одна плоскость заземлена) определяется равенством 2) Для конструкции с заданной толщиной подложки можно определить оптимальную высоту пьезокерамической пластины, позволяющей наиболее эффективно преобразовать внешнее механическое воздействие в индуцируемый электрический сигнал. На рис. 2 показаны графики изменения амплитудных значений разности потенциалов для различной толщины пьезокерамической пластины. В данном случае для металлической подложки толщиной м необходимо использовать пьезокерамическую пластину высотой м. 3) Амплитудные значения электрического импульса при шарнирном закрепление пластины существенно больше, чем в случае жесткого защемления ее контура (см. рис. 2, а, б). Однако шарнирное закрепление характеризуется также большими нормальными механическими напряжениями в центре пластины (см. рис. 3, а, 1- , 2- ). Данная особенность вводит ограничения на величину интенсивности нагрузки При жестком защемлении контура пластины величина значительно меньше (см. рис. 3, б) ( при ), что приводит к увеличению диапазона изменения интенсивности нагрузки . 4) Амплитудное значение потенциала электрического поля по высоте пьезокерамической пластины при разных случаях закрепления изменяется по параболической зависимости (рис. 4, кривая 1, шарнирное закрепление), соответственно график, описывающий изменение аксиальной компоненты вектора напряженности , представляет прямую линию (см. рис. 4, 2). Аналогичная картина наблюдается при анализе биморфных пластин в задачах обратного пьезоэффекта [22]. В заключение отметим, что разработанный в настоящей работе алгоритм расчета позволяет исследовать биморфные конструкции с произвольным количеством упругих и электроупругих слоев.

About the authors

D A Shlyakhin

Samara State Technical University

References

  1. Piezoelectric bimorph bending sensor for shear-stress measurement in fluid flow / D. Roche, C. Richard, L. Eyraud, C. Audoly // Sensors and Actuators. - 1996. - Vol. 55. - P. 157-162.
  2. Yoo J.H., Hong J.I., Cao W. Piezoelectric ceramic bimorph coupled to thin metal plate as cooling fan for electronic devices // Sensors and Actuators. - 2000. - Vol. 79. - P. 8-12.
  3. Sharapov V. Piezoceramic sensors. - Springer Verlag, 2010. - 498 p.
  4. Jurenas V., Bansevicius R., Navickaite S. Piezoelectric bimorphs for laser shutter systems: optimization of dynamic characteristics // Mechanika. - Kaunas: Technologija, 2010. - No. 5(85). - P. 44-47.
  5. Seung-Bok Choi, Young-Min Han. Piezoelectric actuators: control applications of smart materials. - N.-Y.: CRC Press, 2010. - 280 p.
  6. Шарапов В., Сотула Ж. Пьезокерамические преобразователи. Новые технологи проектирования // Электроника:НТБ. - 2012. - № 5. - С. 96-102.
  7. Ультразвуковые колебательные системы для синтеза полимерных композиционных материалов / Д.А. Негров, Е.Н. Еремин, А.А. Новиков, Л.А. Шестель. - Омск: Изд-во Омск. гос. техн. ун-та, 2012. - 128 с.
  8. Шикина В. Е. Расчет частоты колебаний пьезокерамического первичного преобразователя для массового расходомера жидкостей // Изв. высших учебных заведений. Поволжский регион. Технические науки. -2014. - № 1 (29). - С. 54-63.
  9. Ивина Н.Ф., Тагильцев А.А. Анализ собственных колебаний дисковых изгибных пьезопреобразователей с произвольным соотношением размеров [Электронный ресурс] // Электронный журнал «Техническая акустика». - 2005. - № 2. - С. 1-9. - URL: http:// webcenter.ru/ eeaa/ejta.
  10. Янчевский И.В. Минимизация прогибов электроупругой биморфной пластины при импульсном нагружении // Проблемы вычислительной механики и прочности конструкций. - Харьков, 2011. - Вып. 16. - С. 303-313.
  11. Никофоров С.Н. Теория упругости м пластичности. - М.: Госиздат. по арх. и стр-ву, 1955.-284 с.
  12. Сеницкий Ю.Э. Исследование упругого деформирования элементов конструкций при динамических воздействиях методом конечных интегральных преобразований. -Саратов: Изд-во Сарат. ун-та,1985. - 174 с.
  13. Mohammad Amin Rashidifar, Ali Amin Rashidifar. Vibrations Analysis of Circular Plate with Piezoelectric Actuator Using Thin Plate Theory and Bessel Function // American Journal of Engineering, Technology and Society. - 2015. - No. 2(6). - P. 140-156.
  14. Tsaplev V., Konovalov R., Abbakumov K. Disk bimorph-type piezoelectric energy harvester // J. of Power and Energy Eng. - 2015. - No. 3. - P. 63-68.
  15. Jafar Eskandari Jam, Mahmood Khosravi, Nader Namdaran. An exact solution of mechanical buckling for functionally graded material bimorph circular plates // Metall. Mater. Eng. - 2013. - Vol. 19 (1). - P. 45-63.
  16. Adelman N.T., Stavsky Y. Flexural-extensional behavior piezoelectric cilcular plates // J. Acoust. Soc. Amer. -1980. -Vol. 67. - No. 3. - P. 819-822.
  17. Karlash V.L. Resonance Electro-Mechanic Vibration of Piezo-Ceramic Plates // Int. Appl. Mech. - 2005. - Vol. 41. - No. 7. - P. 535-541.
  18. Ватульян А.О., Рынкова А.А. Об одной модели изгибных колебаний пьезоэлектрических биморфов с разрезными электродами и ее приложениях // Изв. РАН. МТТ. - 2007. - № 4. - С. 114-122.
  19. Wang Y., Xu R.Q., Ding H.J. Analytical solutions of functionally graded piezoelectric circular plates subjected to axisymmetric loads // Acta Mechanica. - 2010. - Vol. 215. - Iss. 1-4. - P. 287-305.
  20. Шляхин Д.А. Вынужденные осесимметричные колебания пьезокерамической тонкой биморфной пластины // Изв. РАН. МТТ. - 2013. - № 2. - С. 77-85.
  21. Shlyakhin D.A., Kazakova O.V. Non-Stationary Flexural Fluctuations of a Round Flat Bimorph Plate with Graded-Varying Thickness // Procedia Engineering. - 2014. - Vol. 91. - P. 69-74. doi: 10.1016/j.proeng.2014.12.014
  22. Шляхин Д.А. Динамическая осесимметричная задача электроупругости для жестко закрепленной биморфной пластины // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2015. - № 2. - С. 164-178. doi: 10.15593/perm.mech/2015.2.11.
  23. Shlyakhin D.A., Kazakova O.V. A dynamic axially symmetric goal and its extended solution for a fixed rigid circular multi-layer plate // Procedia Engineering. - 2016. - Vol. 153. - P. 662-666. DOI information: 10.1016/j.proeng.2016.08.219
  24. Гринченко В.Т., Улитко А.Ф., Шульга Н.А. Механика связанных полей в элементах конструкций. - Киев: Наук. думка, 1989. -279 с.
  25. Партон В.З. Кудрявцев Б.А. Электроупругость пьезоэлектрических и электропроводных тел. - М.: Наука,1988. - 470 с.
  26. Снеддон И. Преобразование Фурье. - M.: Изд-во иностр. лит., 1955. - 669 c.
  27. Владимиров В.С. Обобщенные функции в математической физике. - М.: Наука, 1978. -318 с.
  28. Сеницкий Ю.Э. Многокомпонентное обобщенное конечное интегральное преобразование и его приложение к нестационарным задачам механики // Изв. вузов. Математика. - 1991. - № 4. - С. 57-63.
  29. Шляхин Д.А. Вынужденные осесимметричные изгибные колебания толстой круглой жестко закрепленной пластины // Вестн. Самар. гос. ун-та. Естественно-науч. серия. - 2011. - № 8 (89). - С. 142-152.
  30. Шляхин Д.А. Вынужденные осесимметричные изгибные колебания толстой круглой жестко закрепленной пьезокерамической пластины // Вестн. Самар. гос. ун-та. Естественно-науч. серия. - 2012. - № 6 (97). - С. 124-135.

Statistics

Views

Abstract - 131

PDF (Russian) - 43

Cited-By


PlumX


Copyright (c) 2017 Shlyakhin D.A.

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