Асимметричные моды колебаний пьезокерамических цилиндров

Аннотация

И.Ф. Бугаян, И.Н Мощенко

В работе исследуются колебания пьезокерамических цилиндров, склеенных из секторов с тангенциальной поляризацией, и электродами, расположенными между секторами. Рассмотрены асимметричные моды колебаний, возбуждаемые синусоидальными пространственными тангенциальными гармониками электрического потенциала. Показано, что ранее разработанная методика решения уравнений пьезоупругости для такого цилиндра при возбуждении косинусоидальными гармониками подходит и для этой ситуации. Определены в аналитическом виде поля перемещений, деформаций и напряжений. Полученные результаты этой и предыдущей работ позволяют получать полные решения при произвольном типе подключения электродов, через ряды Фурье по известным решениям для гармоник. А также определять весь набор резонансов.
Ключевые слова: пьезоупругость, цилиндр, колебания, азимутальные моды, ряд Фурье, синусоидальные моды, бесконечный ряд, аналитическое решение, резонанс, линейное приближение.

01.04.07 - Физика конденсированного состояния

В качестве ультразвуковых излучателей радиальных колебаний иногда используют толстостенные пьезокерамические цилиндры, склеенные из секторов с антифазной (по отношению к ближайшим соседям) тангенциальной поляризацией. На рис. 1 приведено поперечное сечение цилиндра, склеенного из четырех секторов. Стрелками внутри цилиндра показаны направления поляризации в секторах.

Рис. 1. Поперечное сечение цилиндра, склеенного из четырех секторов. Стрелками внутри цилиндра показаны направления поляризации в секторах, полужирными линиями – электроды.

В реальных устройствах число секторов больше (обычно 12 - 24). В статье исследуется общий случай  с (2 m) секторами и  (2 m)  электродами. Для возбуждения механических колебаний цилиндра на электроды подают переменное напряжение со сдвигом по фазе разности потенциалов на соседних секторах. Сдвиг фаз зависит от типа колебаний, которые необходимо возбудить, но чаще всего используют противофазное включение, как показано на рис. 1. При этом потенциалы φ_i на электродах подчиняются следующему условию:
φ_2i - φ_2i-1 = -φ_2i+112i = u₀ cos(ωt), i=1….m;                     (1)
где u₀ – амплитуда напряжения,  ω – частота, t - время.
Для длинных цилиндров характерны моды колебаний, близкие к плоским, когда смещения практически перпендикулярны оси, а их амплитуда и фаза слабо изменяются в осевом направлении. Целью работы является теоретический анализ таких колебаний в двумерном приближении. Рассматривается идеализированный случай полностью плоских колебаний, для которых смещения uz вдоль оси цилиндра отсутствуют, а остальные смещения ur и uq не изменяются в этом направлении. (В работе используется  цилиндрическая система координат (r, θ, z),  с началом координат на оси и координатой z вдоль оси цилиндра).
Задача решается в рамках линейной теории пьезоупругости, в предположении, что по упругим свойствам материал полностью изотропен, а по пьезоэлектрическим – поляризован в тангенциальном направлении и полностью изотропен в перпендикулярном.
В соответствии с этим, упругие свойства описываются двумя коэффициентами Ламе μ и λ, а пьезоэлектрические, в общем случае, пятерыми не равными нулю пьезомодулями. В цилиндрической системе координат это e_r,rθ= e_z,zθ; e_θ,rr= e_θ,zz; e_θ,θθ (в рассматриваемом плоском случае они сведутся к трем пьезомодулям e_r,rθ; e_r,rθ; e_θ,θθ ). 

Отметим, что противофазность как поляризации, так и электрического подключения позволяет редуцировать задачу, рассматривая одинаковую поляризацию всех секторов и их синфазное подключение. При этом разность потенциалов  на соседних электродах:

φ_2i - φ_2i-1 = φ_2i+1 - φ_2i = u₀ cos(ωt), i=1….m.                          (2)

Что ведет к следующим граничным условиям для электрического потенциала:

                                      (3)

Очевидно, что эти две задачи (прямая и редуцированная) полностью эквивалентны.  Если же в прямой задаче сдвиг фаз между соседними секторами отличается от 180  градусов, то задачу все равно можно редуцировать к задаче с одинаковой поляризацией, но уже с ненулевым фазовым сдвигом. При этом граничные условия для потенциала изменяться, но только по фазе , амплитуда разности потенциалов на любом секторе останется той же:

                                  (4)

В работе рассматривается задача именно в такой, редуцированной постановке, с граничными условиями типа (4) для потенциала. Распределение потенциала внутри пьезокерамики находится из уравнений электростатики:

                                                                         (5)

здесь  и   - векторы напряженности и индукции электростатического поля, связанные соотношениями теории пьезоэффекта [3]:

                     (6)

где  ε_r и ε_q – соответствующие компоненты диэлектрической проницаемости (для пьезокерамики ε в направлении поляризации и  в перпендикулярном направлении различаются), ε₀ – электрическая постоянная, ε_ij – компоненты тензора деформации. Из первого уравнения системы (5) вытекает потенциальность поля напряженности , где электрический потенциал j должен удовлетворять граничным условиям (4) и уравнениям (5), (6). Приближенное решение для потенциала можно найти, пренебрегаю в (6) обратным влиянием упругой подсистемы на электрическую. Введем вектор :

Тогда (5) можно записать

    (7)

Предположим, что для всех решений  (вернее ). Тогда второе уравнение системы (5) можно приближенно заменить

                                             (8)

Физически мы пренебрегли электрическим полем, создаваемы прямым пьезоэффектом по сравнению с электрическим полем внешних источников в конденсаторе. Подставляя  в (8), получим уравнение для :

,                                               (9)

где  - оператор Лапласа.
Это уравнение имеет частное решение

                                                  (10)

удовлетворяющее нашим граничным условиям.  Отметим, что это решение справедливо только для первого сектора, потенциал для остальных нетрудно получить из него, пользуясь условиями непрерывности и заданными величинами сдвига фаз для разности потенциалов на остальных секторах.
Противофазное включение соседних секторов используется для получения полносимметричных аксиальных колебаний. Теоретически такие колебания пьезокерамических цилиндров рассмотрены нами в двумерном приближении в [1]. Для этого случая найдены аналитические решения в виде сходящихся рядов, и получены уравнения для резонансных частот [1]. Разработанный при этом  метод решения уравнений пьезоупругости в [2] обобщается на случай плоских колебаний  с азимутальной  анизотропией. Такие моды колебаний возбуждаются, во-первых, при включении секторов, отличном от противофазного. Во-вторых, они могут возникать в виде резонансных колебаний и при противофазном включении. Теоретически последнее можно показать, исходя из следующих соображений.
При противофазном включении секторов в электрическую цепь, для вышеупомянутой редуцированной задачи приближенное решение для потенциала (10) будет справедливо для всего цилиндра, а не только на первом секторе. Более того, как показано в [1], для полносимметричной моды колебаний оно выполняется точно. Пространственный спектр распределения типа (10)   содержит азимутально-неоднородные гармоники. Как показано в [2] на примере косинусоидальных (по ) гармоник, для каждой из них существует соответствующее решение уравнений пьезоупругости для нашей задачи. При этом сумма решений по всем гармоникам не совпадает с рядом Фурье (по ) полносимметричного однородного решения [1] для того же потенциала (10). То есть задача даже в приближении линейной теории неоднозначна, Для одного потенциала существует несколько решений. Это не противоречит теореме единственности для линейной пьезоупругости. Теорема единственности выполняется только для односвязных объектов, а наш объект двухсвязный.   
Для многосвязных объектов выбор реализуемого в действительности решения нужно осуществлять на основе анализа динамической устойчивости. Однако это верно лишь для областей вне резонансов. В районе резонанса какого-то одного решения из общих соображений можно предположить, что энергетически выгодным  будет именно резонансное решение. Так как резонансы полносимметричной моды и мод с азимутальной компонентой не совпадают, то следует ожидать, что в резонансном режиме будут наблюдаться все, и однородные, и неоднородные (по ) гармоники.
Как показано в [2], для противофазного включения соседних секторов функцию распределения потенциала по азимутальному углу можно разложить в ряд Фурье либо по синусам (на интервале [0, 2π]), либо по косинусам  (расширяя четным образом область определения с интервала [0, 2π] на интервал [-2π, 2π]). В случае же электрического подключения общего вида, ряд Фурье, для функции распределения электрического потенциала, содержит как синусы, так и косинусы. Возбуждение колебаний пьезокерамического цилиндра косинусоидальными гармониками рассмотрено в [2].  Целью настоящей работы являлось исследование колебаний цилиндра под действием синусоидальных (по азимутальному углу) гармоник электрического потенциала.
Таким образом, в работе рассматривается следующая постановка задачи. Задано распределение электрического потенциала:

                                                                   (11)

где n – любое целое.
Необходимо определить в двухмерном приближении распределение смещений, деформаций и напряжений в пьезокрамическом цилиндре, обусловленных потенциалом (11).
Для решения поставленной задачи необходимо найти решение уравнений движения:

                  (12)

здесь ρ – плотность; - i-я компонента вектора перемещения;  - компоненты тензора напряжения, связанные с ε_ij – компоненты тензора деформации Гука [3]:

                                   (13)

здесь - относительное изменение объема; α и μ – коэффициенты Ламе.
Напомним, что задача рассматривается в полярной системе координат. В соответствии со сделанными допущениями двухмерности следует двухмерность уравнений (12) (u_z=0, u_r= u_r(r, θ) и u_q= u_q(r, θ)), равенство нулю z компонент векторов напряженности и индукции электростатического поля, а также iz компонент тензора деформации (i - любое) и rz и qz компонент тензора напряжений.
Определяя по (11) напряженность электрического поля, подставляя ее в (13), а (13) в (12) и выражая деформации εij через перемещения , получим обычные уравнения движения в перемещениях с эффективной объемной силой:

      (14)

здесь E и σ – модуль Юнга и коэффициент Пуассона (), а массовая сила имеет следующие компоненты:

                               (15)

Для решения уравнения (14) воспользуемся методом Ламе. Представим где  т.е.  (). Отметим, что в двухмерном случае векторный потенциал  сводится к скалярному  и компоненты вектора смещения следующим образом выражаются через оба скалярные потенциалы:

                                                          
                                                       (16)

Выразим так же массовую силу из (15) через соответствующий скалярный потенциал:

    .        (17)

Подставим и в уравнение движения (14) и получим два скалярных уравнения для потенциалов  и :

                                                   (18)
                                                                  (19)

где  и  - продольная и поперечная скорости звука соответственно. 
Уравнение для поперечной составляющей (19) является чисто волновым, его решения известны:

,  (20)

где  - функции Бесселя первого и второго рода, k1 – волновое число,  , α₂, β₂, D₂ и G₂ – константы, определяемые из граничных условий на боковых поверхностях цилиндра и на торцах первого сектора. Так как эффективная сила, связанная с пьезоактивностью, не входит ни в уравнения (19), ни в граничные условия, то поперечная составляющая за счет линейных эффектов не возбуждается. В принципе, она может возбудиться в импульсном режиме или на резонансе этой моды  за счет нелинейного взаимодействия с продольной составляющей, но рассмотрение таких эффектов выходит за рамки статьи.
Уравнение для продольной  составляющей (18) представляет собой волновое уравнение с эффективной массовой силой . Ранее [2] нами было получено аналогичное уравнение при исследовании колебаний такого же пьезокерамического  цилиндра, возбуждаемых косинусоидальными гармониками. Единственное отличие состоит в том, что для синусоидального возбуждения в эффективную силу  входит , а не , как в [2]. Методика решения уравнений такого типа, предложенная в  [2], легко обобщается и на рассматриваемый случай.
Будем искать решение уравнения (18) в виде тогда оно переходит в

.
Или

                                                (21)

где . Сделаем вложенную цепочку замен переменных:

 ; ;  ;                
; ;;
.
.
.
;                                            
;                                          
;                                       (22)

                
И так далее, вплоть до

Отметим, что в этой цепочке

                      (23)

При этом для любых  и  можно найти  такое, что при   выполняется: т.е.  и в пределе при      (24)
При этом мы перейдем к уравнению , т.е. к уравнению Гельмгольца, его частное решение

,             (25)

где  - функция Бесселя первого рода n-го порядка,  - функция Бесселя второго рода n-го порядка,   - константы.
Обращая цепочку замен (22), определим Ф₁

      
где   и члены ряда Bi определены в (23). Отметим, что в соответствии с (23), ряд S мажорируетсягеометрической прогрессией, а значит, сходится абсолютно. Более того, его можно почленно дифференцировать

,                                                  (26)

Для этих рядов также нетрудно показать, что и они сходится абсолютно и их можно почленно дифференцировать

   (27)

и ряды также сходится абсолютно.
Таким образом, потенциал Ф₁ выражается через известный потенциал  (25) и бесконечный абсолютно сходящийся ряд
.                                                                                                 (28)
Это дает возможность определить потенциал Ф поля перемещений, а по нему определить сами перемещения,  деформации и напряжения. Последние при этом будут выражаться через производные от   и абсолютно сходящиеся ряды (26) и (27). Этим решается задача об описании колебаний, возбуждаемых одной гармоникой (11) электрического потенциала.  Ввиду линейности используемой теории, колебания цилиндра, возбуждаемые всем электрическим потенциалом (10), описываются суперпозицией найденных решений. При этом каждое отдельное решение содержит восемь неизвестных констант.
Найденные решения уравнений (12 - 13) должны удовлетворять граничным условиям для упругих переменных:

                                                        (29)

Граничных условия у нас всего два и им должно удовлетворять только суммарное решение. Мы можем разложить каждое из (29) в ряд Фурье и потребовать удовлетворения своему граничному условию для каждой гармоники. Однако для нулевых граничных условий все гармоники будут также нулевые, и получается, что условию (29) должно удовлетворять каждое решение по отдельности. В частности, это будет выполняться, если для каждой найденной моды колебаний, взять  , а D₁ и G₁ определять из (29).  Таким образом, всегда существует, по крайней мере, одно решение для каждой синусоидальной гармоники (11), удовлетворяющее граничным условиям. Ранее [1] мы нашли аналогичные решения для косинусоидальных гармоник. При любом типе подключения электродов электрический потенциал может быть разложен в ряд Фурье по азимутальному углу, содержащий в общем случае как синусоидальные, так и косинусоидальные гармоники. Таким образом, полученные результаты позволяют получать полные решения при произвольном типе подключения электродов, как сумму найденных решений по всем гармоникам.

Литература

1. Мощенко И.Н. и др., Расчет осесимметричных колебаний пьезокерамического цилиндра с тангенциальными поляризацией и возбуждением //Инженерный вестник Дона, 2009. №1. http://www.ivdon.ru/magazine/archive/n1y2009/117/ (доступ свободный) — Загл. с экрана. — Яз. рус.
2. Мощенко И.Н., Товаровская Н.М. Радиально - круговые моды  колебаний пьезокерамического цилиндра с тангенциальными поляризацией и возбуждением //Инженерный вестник Дона, 2009. №1. http://www.ivdon.ru/magazine/archive/n2y2009/125/ (доступ свободный) — Загл. с экрана. — Яз. рус.
3. Сиротин Ю.И., Шаскольская М.П. Основы кристаллофизики – М.: Наука, 1979. – 639 с.