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

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

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

Ростовский строительный университет, Северо-Кавказский научный центр высшей

школы ЮФУ

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

2

О

Цсс^иЛ)

о

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

В реальных устройствах число секторов больше (обычно 12 - 24). В статье исследуется общий случай с (2 т) секторами и (2 т) электродами. Для возбуждения механических колебаний цилиндра на электроды подают переменное напряжение со сдвигом по фазе разности потенциалов на соседних секторах. Сдвиг фаз зависит от типа колебаний, которые необходимо возбудить, но чаще всего используют противофазное включение, как показано на рис. 1. При этом потенциалы ф; на электродах подчиняются следующему условию:

ф21 - ф21-1 = -фя+1 + ф21 = ио соб(ю^, ¡=1.. ..т; (1)

где и0 - амплитуда напряжения, ю - частота, 1 - время.

Для длинных цилиндров характерны моды колебаний, близкие к плоским, когда смещения практически перпендикулярны оси, а их амплитуда и фаза слабо изменяются в осевом направлении. Целью работы является теоретический анализ таких колебаний в двумерном приближении. Рассматривается идеализированный случай полностью плоских колебаний, для которых смещения и2 вдоль оси цилиндра отсутствуют, а остальные смещения иг и ие не изменяются в этом направлении. (В работе используется цилиндрическая система координат (г, 9, ъ), с началом координат на оси и координатой ъ вдоль оси цилиндра).

Задача решается в рамках линейной теории пьезоупругости, в предположении, что по упругим свойствам материал полностью изотропен, а по пьезоэлектрическим -поляризован в тангенциальном направлении и полностью изотропен в перпендикулярном.

В соответствии с этим, упругие свойства описываются двумя коэффициентами Ламе ц и X, а пьезоэлектрические, в общем случае, пятерыми не равными нулю пьезомодулями. В цилиндрической системе координат это ег,г9= еъ,ъ9; е9,гг= е9,ъъ; е9,99 (в рассматриваемом плоском случае они сведутся к трем пьезомодулям ег,г9; е9,гг; е9,99).

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

ф21 - Ф21-1 = Ф21+1 - ф21 = uo cos(rot), i=1.. ..m. (2)

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

в = +0; ф0 = -шщ cos ot,

в = -0; ср2т = +шщ cos ot,

— .

в = — Г;ф = + u0 cos ot,

ш

i = 1...2ш . (3)

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

вг =—, Ф = Ф + u0 c0s(ot + X); (АЛ

ш (4)

i = 1...Ш.

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

\rotE = 0,

\ - (5)

[divD = 0,

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

D = ££E + 4— e £ й,

r r 0 r r,гв гв ’

D„ = ££,E + 4— e„ E + 4— e„ „,

в в 0 в e,rr rr в,вв вв? (6)

где 8r и 8е - соответствующие компоненты диэлектрической проницаемости (для пьезокерамики 8 в направлении поляризации и в перпендикулярном направлении различаются), 80 - электрическая постоянная, 8ij - компоненты тензора деформации. Из первого уравнения системы (5) вытекает потенциальность поля напряженности

E = - sradrn „

* ^, где электрический потенциал ф должен удовлетворять граничным

условиям (4) и уравнениям (5), (6). Приближенное решение для потенциала можно найти,

пренебрегаю в (6) обратным влиянием упругой подсистемы на электрическую. Введем

вектор А :

Ar =

Ав = 4—eeJrr£rr + 4—ee,z£z + 4—eе,ееEее,

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

- - дЕ -

divD = £в£divE + (e - £ ) —- + divA = 0.

r dr (7)

Предположим, что для всех решений (ЛгуА ~ 0 (вернее

- - дБ

ЛгуА << (£в£аЛгуБ + (£г - ев)—L)

дг ). Тогда второе уравнение системы (5) можно

приближенно заменить

ЛгуБ + (£' -£<в) дБ- = 0.

БвБо дг (8)

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

Е = - ятайф ,оч т

ь ^ в (8), получим уравнение для ^ :

Ат-(е- -£‘) дт о

£в£о д г (9)

где А - оператор Лапласа.

Это уравнение имеет частное решение

т = | Г ^ \0,ш), (Ш)

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

Противофазное включение соседних секторов используется для получения полносимметричных аксиальных колебаний. Теоретически такие колебания пьезокерамических цилиндров рассмотрены нами в двумерном приближении в [1]. Для этого случая найдены аналитические решения в виде сходящихся рядов, и получены уравнения для резонансных частот [1]. Разработанный при этом метод решения уравнений пьезоупругости в [2] обобщается на случай плоских колебаний с азимутальной анизотропией. Такие моды колебаний возбуждаются, во-первых, при включении секторов, отличном от противофазного. Во-вторых, они могут возникать в виде резонансных колебаний и при противофазном включении. Теоретически последнее можно показать, исходя из следующих соображений.

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

Для многосвязных объектов выбор реализуемого в действительности решения нужно осуществлять на основе анализа динамической устойчивости. Однако это верно лишь для областей вне резонансов. В районе резонанса какого-то одного решения из общих соображений можно предположить, что энергетически выгодным будет именно резонансное решение. Так как резонансы полносимметричной моды и мод с азимутальной

компонентой не совпадают, то следует ожидать, что в резонансном режиме будут наблюдаться все, и однородные, и неоднородные (по в ) гармоники.

Как показано в [2], для противофазного включения соседних секторов функцию распределения потенциала по азимутальному углу можно разложить в ряд Фурье либо по синусам (на интервале [0, 2п]), либо по косинусам (расширяя четным образом область определения с интервала [0, 2п] на интервал [-2п, 2п]). В случае же электрического подключения общего вида, ряд Фурье, для функции распределения электрического потенциала, содержит как синусы, так и косинусы. Возбуждение колебаний пьезокерамического цилиндра косинусоидальными гармониками рассмотрено в [2]. Целью настоящей работы являлось исследование колебаний цилиндра под действием синусоидальных (по азимутальному углу) гармоник электрического потенциала.

Таким образом, в работе рассматривается следующая постановка задачи. Задано распределение электрического потенциала:

(р = А §Ыпв)со$,(Ш), (11)

где п - любое целое.

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

Для решения поставленной задачи необходимо найти решение уравнений движения:

дагг , 1 дагв , а гг — авв_д 2иг

дг г дв г ді2 ’

дугв , 1 2 д ч

дг г дв г гв ді2 ’

и • а..

здесь р - плотность; 1 - і-я компонента вектора перемещения; ,к - компоненты тензора

напряжения, связанные с 8у - компоненты тензора деформации Гука [3]:

агг = 2/^гг +^0 - Єв,ггЕв , авв = ^^вв + ^0 — ев,ввЕв, агв = 2М£гв — Єг,гвЕг ,

0

(13)

е — £ + £

здесь 0 гг вв- относительное изменение объема; Л и ^ — коэффициенты Ламе.

Напомним, что задача рассматривается в полярной системе координат. В

соответствии со сделанными допущениями двухмерности следует двухмерность

уравнений (12) (иъ=0, иг= иг(г, 9) и ив= ив(г, 9)), равенство нулю ъ компонент векторов напряженности и индукции электростатического поля, а также \ъ компонент тензора деформации (1 - любое) и гъ и въ компонент тензора напряжений.

Определяя по (11) напряженность электрического поля, подставляя ее в (13), а (13)

в (12) и выражая деформации 81] через перемещения и, получим обычные уравнения движения в перемещениях с эффективной объемной силой:

Л 2и Б Б

р—и —------------Аи +---------------grad (Лгуи) + Р, (14)

Ли2 2(1 + а) 2(1 + а)(1 - 2а)

здесь Е и о - модуль Юнга и коэффициент Пуассона

Б Ба

/л =-----------, А —----------------

( 2(1 + а) (1 2а)(1 + а) ), а массовая сила имеет следующие компоненты:

Рг = —Ап ■ псо%(пв)со%(аг);

/ (15)

= —0_0^ а _ п2 оо$,(п6)оо$,(Ш). г

Для решения уравнения (14) воспользуемся методом Ламе. Представим и = V + й ,^ сИуЙ = 0, гогУ = 0, и = gгadФ + го№

где т.е.

У = gradФ, й = го№ ч ~

( ). Отметим, что в двухмерном случае векторный потенциал

¥ = (0,0, ¥)

сводится к скалярному 4 ’ ’ 7 и компоненты вектора смещения следующим

образом выражаются через оба скалярные потенциалы:

дФ 1

и =-----------1----------

г дг г дв

1 дФ д¥

иг

в г дв дг (16)

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

р = рgгadv, у = ее.ввА ■пСО^со^,). (17)

р г

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

1 д2 Ф 1

ДФ--------------------=------------у'

2 Я+2 2 У’

С, д С,

'I

1 д2 ¥ Д¥ = 1 д т

(18)

(19)

С с,

где 1 и і - продольная и поперечная скорости звука соответственно.

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

Ч* = (Ц,./„ (кхг) + С2Мп (кхг)){а2 софів) + у32 5т(пв))со$(саі), (20)

J N кі = ~

где ” ’ ” - функции Бесселя первого И второго рода, кі - волновое ЧИСЛО, с*, 0(2, (Зг,

о2 и а2 - константы, определяемые из граничных условий на боковых поверхностях цилиндра и на торцах первого сектора. Так как эффективная сила, связанная с пьезоактивностью, не входит ни в уравнения (19), ни в граничные условия, то поперечная составляющая за счет линейных эффектов не возбуждается. В принципе, она может возбудиться в импульсном режиме или на резонансе этой моды за счет нелинейного взаимодействия с продольной составляющей, но рассмотрение таких эффектов выходит за рамки статьи.

Уравнение для продольной составляющей (18) представляет собой волновое уравнение с эффективной массовой силой у. Ранее [2] нами было получено аналогичное уравнение при исследовании колебаний такого же пьезокерамического цилиндра, возбуждаемых косинусоидальными гармониками. Единственное отличие состоит в том, что для синусоидального возбуждения в эффективную силу у входит СОБ(пв) , а не БІп(пв), как в [2]. Методика решения уравнений такого типа, предложенная в [2], легко обобщается и на рассматриваемый случай.

Будем искать решение уравнения (18) в виде Ф = Фх соб^, тогда оно переходит в

Дф + е 0.а0пАп соб(п^) = _«! Ф

c2 г c2

с е ' се

Или

Дф + А„, = _k 2ф . (21)

г

Ю е^а^мА

1 1ЛУ а У.УУ П /~ч

где к = —, А5п =--------- ---. Сделаем вложенную цепочку замен переменных:

Се С1

Дф , А5п с°<п0) _ к2ф , , А5соъ(нв) ■ г ^ , А5саъ{п@) ■ г

1 г 1 ’ V

^ к2А, со^Спв) ■ г

дф---------5п---— = -к 2фп;

Ф1 + ' 5 4 2 7 I = -к2Ф1; Ф2 = Ф1 + -?------------------—

1 1 - п2 J 1 2 1 (1 - п2)

Д

(1 - п2) 2

к 1 А5п С°Ъ(пв) ■ г3 ] = _к2ф ,ф =ф к2 А5С°$(пв) ■ гГ

2 (1 -п2)(32 -п2) ^ - 3 2 (! -п2)(32 -п2)

ДФ +■

(-1)т-1 к2(т-1)А5„ со8(п0) ■ г2т (1 - п2)(32 - п2)...((2т - 3)2 - п2)

Д

Ф_ +

(-1т-1 к2(т-1)А5и со8(п0) ■ г2

(1 - п2)(32 - п2)...((2т - 1)2 - п2)

= -к2ф ;

Ф , = Ф --

т+1 т

(-1)тк2(т-1) А5п соБ(п^) ■ г2т-1

(1 -п2)(32 -п2)...((2т-1)2 -п2)

22

(22)

И так далее, вплоть до т ^ <х>.

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

Ф , = Ф _ В

т+1 т т ’

Вт =

Вт Вт-1

(-1)тк2(т-1) Аи СОБ(И^) ■ г2т-1

(1 -п2)(32 - п2)...((2т -1)2 - п2): к2 г2

(23)

((2т_ 1)2 _п2)'

При этом для любых к и г можно найти N. такое, что при т > N. выполняется:

к2 Г 2

уг;-----"Т2-----Л < 1.т.е. |вт < |Вт_^ и в пределе К1 ^ 0. при т (24)

^(2т _ 1] _ п ]

При этом мы перейдем к уравнению

1д

г дг

дФ

= -к2 Ф

т.е. к уравнению

К дг у

Гельмгольца, его частное решение

фда = {.^П (кг) + °1 КП (кгЖ°а с°8(п0) + Рх 81П(п^)) , (25)

где ,/и - функция Бесселя первого рода п-го порядка, N п - функция Бесселя второго рода

п-го порядка, 01гС1г а± и 0± - константы.

Обращая цепочку замен (22), определим Ф1

г

Ф , =Ф + в ,

т-1 т т ~

ад

ф=ф_ + в =ф+ в + в =Ф+ в + В + в =ф + £ в =Ф + s,

1 22 323 4234 ад / л , ад 5

г=2

ад

где 5 = £ в и члены ряда Б1 определены в (23). Отметим, что в соответствии с (23), ряд

г=2

s мажорируется геометрической прогрессией, а значит, сходится абсолютно. Более того, его можно почленно дифференцировать

дя ” дБ дя ” дБ

— = £—Ч — = £—L, (26)

дг 1=1 дг до 1=1 дО

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

почленно дифференцировать

д2 я _ “ д2 в д2 5 _ ^ д2 в д2 5 _ “ д2 в

дг2~£~дТГ’ ~дгдО~£дгдоо до~£Ю’ (27)

и ряды также сходится абсолютно.

Таким образом, потенциал Ф^ выражается через известный потенциал Фш (25) и бесконечный абсолютно сходящийся ряд

Ф1 =Фад+ 5 . (28)

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

Найденные решения уравнений (12 - 13) должны удовлетворять граничным условиям для упругих переменных:

г = г{, г2° ^гг = 0- (29)

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

Литература

1. Мощенко И.Н. и др., Расчет осесимметричных колебаний пьезокерамического цилиндра с тангенциальными поляризацией и возбуждением //Инженерный вестник Дона, 2009. №1. http://www.ivdon.ru/magazine/archive/n1v2009/117/ (доступ свободный) — Загл. с экрана. — Яз. рус.

2. Мощенко И.Н., Товаровская Н.М. Радиально - круговые моды колебаний пьезокерамического цилиндра с тангенциальными поляризацией и возбуждением //Инженерный вестник Дона, 2009. №1. http://www.ivdon.ru/magazine/archive/n2v2009/125/ (доступ свободный) — Загл. с экрана. — Яз. рус.

3. Сиротин Ю.И., Шаскольская М.П. Основы кристаллофизики - М.: Наука, 1979. - 639 с.