Голуб М.В., Шпак А.Н., Buethe I., Fritzen C.P. «Моделирование гармонических колебаний и определение резонансных частот полосового пьезоэлектрического актуатора методом конечных элементов высокого порядка точности» Вычислительная механика сплошных сред, 8, № 4, с. 397-407 (2015)

Моделируется динамическое поведение прямоугольного полосового пьезоэлектрического актуатора с помощью метода конечных элементов высокого порядка точности. Строится гармоническое решение в частотной области, что позволяет применить преобразование Лапласа для построения решения во временной области. В качестве аппроксимационных полиномов и тестовых функций используются полиномы Гаусса–Лежандра–Лобатто. Рассматривается два варианта граничных условий. В первом случае все границы пьезоактуатора свободны от напряжений, электрический потенциал задан на нижней границе, в то время как на верхней границе он равен нулю, также нулевыми являются электрические перемещения на боковых границах. Во втором случае на нижней границе приложены некоторые нормальная и касательная нагрузки, при этом на левой боковой границе пьезоактуатор жестко защемлен. Составляется система линейных алгебраических уравнений относительно функций перемещений и электрического потенциала в узловых точках. Вектор-столбец правой части системы формируется в соответствии с граничными условиями из электрических потенциалов на границе, а в случае смешанных граничных условий, еще и из нормальных и касательных напряжений на нижней границе актуатора. Результаты моделирования сравниваются с результатами, полученными с помощью пакета Comsol Multiphysics: проводится сопоставление функций перемещений, напряжений, электрического потенциала и электрических перемещений, а также максимальных и минимальных значений этих величин. Анализируется поведение пьезоактуатора в зависимости от граничных условий и частоты гармонических колебаний. Рассчитываются резонансные частоты колебаний актуатора для различных граничных условий, изучается соответствующие резонансным частотам собственные формы колебаний.