CC BY
77
Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2017. Т. 21, № 4. С. 773-785 ISSN: 2310-7081 (online), 1991-8615 (print) d http://doi.org/10.14498/vsgtu1564
УДК 539.3
Вынужденные осесимметричные колебания круглых многослойных биморфных пластин
Д. А. Шляхин, О. В. Ратманова
Самарский государственный технический университет,
Россия, 443100, Самара, ул. Молодогвардейская, 244.
Аннотация
Представлена методика расчета круглых сплошных многослойных биморфных пластин и получены новые замкнутые решения осесиммет-ричных динамических задач прямого и обратного пьезоэффектов. В общем случае при исследовании электроупругого (пьезокерамического) и упругого слоев математическая формулировка рассматриваемых задач включает уравнения движения и Максвелла в пространственной постановке относительно компонент вектора перемещений и потенциала электрического поля, а также соответствующие начально-краевые условия. Рассмотрены случаи шарнирного и жесткого закрепления внешнего контура конструкции. Для исследования связанных линейных задач применяется математический аппарат в виде метода конечных интегральных преобразований Фурье—Бесселя и обобщенного интегрального преобразования (КИП). При этом на каждом этапе исследования использовалась процедура приведения граничных условий к виду, позволяющему выполнять соответствующее преобразование.
Построенные расчетные соотношения позволяют обосновать конструктивные решения многослойных пьезокерамических преобразователей, а именно подобрать геометрические размеры и физические характеристики используемых материалов, определить размеры разрезных круговых электродов, позволяющие наиболее эффективно преобразовать внешнее электрическое воздействия в механические колебания при различной частоте. Кроме этого, появляется возможность проанализировать напряженно-деформированное состояние, характер изменения электрического поля, а также частотный спектр собственных осесим-метричных колебаний рассматриваемых систем.
Ключевые слова: биморфная пластина, электроупругость, нестационарная нагрузка, интегральные преобразования.
Краткое сообщение
3 ©® Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.org/licenses/by/4.0/deed.ru) Образец для цитирования
Шляхин Д. А., Ратманова О. В. Вынужденные осесимметричные колебания круглых многослойных биморфных пластин // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2017. Т. 21, № 4. С. 773-785. doi: 10.14498/vsgtu1564. Сведения об авторах
Дмитрий Аверкиевич Шляхин& http://orcid.org/0000-0003-0926-7388
доктор технических наук, доцент; заведующий кафедрой; каф. сопротивления материалов
и строительной механики; e-mail: [email protected]
Олеся Викторовна Ратманова © http://orcid.org/0000-0001-7591-4921 ассистент; каф. сопротивления материалов и строительной механики; e-mail: [email protected]
Получение: 14 октября 2017 г. / Исправление: 9 декабря 2017 г. / Принятие: 18 декабря 2017 г. / Публикация онлайн: 29 декабря 2017 г.
Введение. В различных электроакустических, оптических и пневматических приборах применяются пьезокерамические преобразователи в виде тонких биморфных пластин [1—3]. Данные конструкции представляют собой симметричные или асимметричные многослойные системы, состоящие из электроупругих и упругих тонких пластин. Их широкое применение объясняется высокими эксплуатационными параметрами, низкой стоимостью и простотой в управлении. Изгибные колебаний в рассматриваемых конструкциях создаются с помощью электрического напряжения, приложенного к электро-дированным поверхностям пьезокерамических пластин (явление обратного пьезоэффекта), или в случае действия механической нагрузки (прямой пье-зоэффект).
Для описания работы и повышения функциональных возможностей би-морфных пластин возникает необходимость в разработке математических моделей и построении общих решений, позволяющих описать сложную картину взаимодействия полей напряжений различной физической природы [4]. В настоящее время для реализации данной проблемы, как правило, применяются прикладные теории для тонких пластин в условиях установившегося режима вынужденных колебаний при использовании различных допущений о характере распределения напряженности электрического поля в пьезокерамиче-ских пластинах [5—11]. Полученные приближенные расчетные соотношения используются в дальнейшем для расчета и конструирования приборов различного назначения [12-14].
Замкнутые решения динамических начально-краевых задач электроупругости в пространственной постановке при аксиальной поляризации пьезоке-рамического материала получены в работах [15-17] только для однородного диска (цилиндра). В [15] рассматривался незакрепленный цилиндр при учете смешанных краевых условий на торцевых поверхностях, а в [16] исследовалось напряженно-деформированное состояние жестко закрепленной толстой пластины в случае действия осесимметричной механической нагрузки. Здесь также можно отметить одну из работ, посвященную определению спектра частот собственных осесимметричных и неосесимметричных колебаний [17] пьезокерамических цилиндров.
Таким образом, в настоящее время отсутствуют замкнутые решения динамических задач электроупругости для многослойных конструкций в пространственной постановке. Для решения этого вопроса в данной работе представлена методика расчета круглых многослойных шарнирно и жестко закрепленных биморфных сплошных пластин в случае действия электромеханической нагрузки. Построенные расчетные соотношения позволяют дать качественную и количественную оценку связанности электрических и механических полей напряжений в многослойных конструкциях. Частные решения данного подхода приведены в публикациях авторов [18,19].
1. Постановка задачи. В общем случае круглая биморфная конструкция, занимающая в цилиндрической системе координат область О = {(г, в, г) : 0 ^ г ^ Ь, 0 ^ в ^ 2п, 0 ^ г ^ Л*}, представляет многослойную систему в ви-
де жестко соединенных между собой электроупругих (пьезокерамических) и упругих тонких пластин. Для определенности принимаем, что она состоит из металлической подложки толщиной Н\ и двух пьезокерамических пластин высотой Н\, выполненных из материала гексагональной системы класса 6шш (Л* = 2Н\ + Л2). Продольно-поперечные осесимметричные колебания индуцируются за счет подведения к торцевым электродированным поверхностям пьезопластин с различным направлением вектора аксиальной поляризации электрического напряжения V* (г*, Ь*) (рис. сверху) или действии на лицевых плоскостях конструкции механических нагрузки в виде нормальных напряжений д*(г*,Ь*) (рис. снизу). В дальнейшем все рассматриваемые величины будут обозначаться с символом *, а безразмерные — без символа. Рассматривается жесткое и шарнирное закрепления внешнего контура электроупругой системы (см. рис.).
Математическая формулировка рассматриваемых осесимметричных начально-краевых задач в безразмерной форме в цилиндрической системе координат включает следующие соотношения [4]:
- дифференциальные уравнения движения и Максвелла относительно безразмерных компонент вектора перемещений и (г, г, Ь), Ш (г, г, Ь), а также
безразмерного потенциала электрического поля ^>(г, г,£):
с(5 д2и с}? + с55 д2Ж
А Уи + С55__
дг с(в) дг2
с11
+
13
я2т
у55
с
М 11
дгдг
+
+
е31 + е15 дУ ф(в)
д 2и
с5?у д^
е33 дгдг д£2
+ с3з) + С# + с55) ^ Уи+
сМ дг сМ дг2 с11 с
с М дг 11 с11
+е*5 У д^ + д^ - Фм
дг дг2 д£2
е33
ез1 + е15 + е^у — + д2^ е33 дг е33 дг дг2
Си£11 у ^ _ С(1)^33 д2^ = 0
33
дг
33
дг2
(1)
где
{и, Ж, г, г} = {и*, Ж*, г*, г*}/Ь, ^ = ^еззДЬС^),
I = ^С(1)/р(1) — время; и*(г*, г*, ¿*), Ж*(г*,г*,£*) —радиальная и аксиальная компоненты вектора перемещений; (г*,г*,£*) —потенциал электрического поля; ет&, вц, в33 — пьезомодули и коэффициенты диэлектрической проницаемости электроупругого материала, т, к = = 1, 2,..., 5; р(а), Стк — объемная плотность и модули упругости электроупругого (8 = 1) и упругого (8 = 2) материалов; Ф(1) = 1, Ф(2) =
С(2)Р(1) я 1
= С11 р = + 1 •
= С(1)р(2) , У = дг + г ;
механические граничные условия:
(0, г, ¿)| < те, |и(0,г,£)| < те;
(2)
жесткое закрепление:
Ж (1, г, ¿) = 0, и (1,г,£) = 0;
(3)
шарнирное закрепление: Ж (1,г,£) = 0,
оу
1г=1
с(1 ди + С2) и + с1з) дЖ + ез1 д^'
(1) дг с11
(1) г с11
С(1) дг с11
езз дг ]
г=1
0;
(4)
и=о
ь=о
с
(1) 13
с с
уи +
С33) дЖ д^
С(1) дг дг с11
2=0
9,
55 + д^ + е15
дг / е33 дг -
(1) дг с11
2=0
0
0
Oz
Or
\z=h
\z=h
13 VU + 33
с
+ df. с(1) dz dz с11
(1) il
с55) ( dW + dU \ + ei5 df
C(1H dr dz ) e33 dr -
z=h
= 0,
z=h
(6)
с
(i) 13
VU +
+ dz-
C33) dW , df
Lcji) C31 -c55) ( dW + dU4 + ei5 df
-с(1) V dr dz ) e33 dr
z=hi
C((3^ C(3) dW
13 VU + 33
C
u) C(1) dz Cii
( dW dU\
z=hi+h2 Lcdr dz J -
11
C
(2) 55
z=hi
z=hi+h2
(7)
Ш (г + 0) = Ш (г - 0), и (г + 0) = и (г - 0); электрические граничные условия: ((0, г, Ь) < ж,
Dr
r=1
Cj1)eii df + ei5 (dW
+ —4
e33 dr ' e33 V dr dz J -
r=1
= 0;
а) задача обратного пьезоэффекта
((г, 0, Ь) = V(г, Ь)/2, ((г, Л, Ь) = -V(г, Ь)/2, ((г, Н\,Ь) = ((г, Н\ + ¡¡2, Ь) = 0;
б) задача прямого пьезоэффекта
(9)
Dz
z={0,hi,hi+h2,h} =
начальные условия:
Cji)g33 df + eai VU + dW e23 dz e33 dz
dU
z={0,h1,h1+h2,h}
U(r,z, 0) = Uo(r,z), — = Uo(r,z),
dt t=o
W(r,z, 0) = Wo(r,z),
dW
dt
(10)
(11)
t=o
= Wo(r, z),
где {ио, ио, Шо, Щ, Л, Нъ ¡2} = [Щ, и*, Ш**, Ш* ,Н* ,Н\,Ц}/Ь; арт — компоненты тензора напряжений, р,т = г, г; Ог, — компоненты вектора индукции электрического поля; и*, и*, Ш0*, Ш* — известные в начальный момент времени перемещения и их скорости; д = о*/С^; V = V *езз/(ЬС(11)).
При исследовании упругой среды (в = 1) система (1) состоит только из уравнений движения, сформулированных относительно компонент вектора перемещений. Выражения (2) и первое неравенство (8) являются условиями ограниченности решения в центре пластины, равенство Ог|г=1 = 0 (8)
0
учитывает отсутствие электродного покрытия на цилиндрической поверхности пластины, а выражения (10) в задаче прямого пьезоэффекта означают подключение электродированной поверхности пьезокерамической пластины к измерительному прибору с большим входным сопротивлением, что соответствует режиму «холостого хода» (отсутствию свободных электрических зарядов). Формулы (7) позволяют учесть условия совместности напряжений и перемещений в плоскости жесткого соединения пластин. Кроме этого, в задаче обратного пьезоэффекта первое равенство (5) приравнивается к нулю.
2. Построение решения. Определение общего интеграла рассматриваемых краевых задач электроупругости выполняется путем последовательного применения конечных интегральных преобразований: Фурье—Бесселя [20] по радиальной переменной и обобщенного преобразования [21] по аксиальной координате. При этом на каждом этапе предварительно выполняется процедура приведения соответствующих граничных условий к виду, позволяющему использовать алгоритм преобразования.
Преобразование Фурье—Бесселя [20] позволяет удовлетворить только смешанные граничные условия. Для выполнения данного требования необходимо произвести следующие замены граничных условий (3), (4):
- при жестком защемлении пластины первое равенство (3) заменяется условием наличия касательных напряжений N1(2, *) на цилиндрической поверхности пластины:
|г=1
С55 /дЖ + + е15 Сп V дг дг / е33 дг
= N1(3,*); (12)
- в случае шарнирного закрепления последнее соотношение (8) заменяется условием наличия потенциала электрического поля ^(г,*) на цилиндрической поверхности при г = 1:
<р(1,г,*) = ^ч(г, *). (13)
Кроме этого, в общем случае вводятся новые функции и(г, г,*), -ш(г, г,*), ф(г, г,*), связанные с и (г, г,*), Ж (г, г,*), ^>(г, г,*) соотношениями
и (г, г, *) = _ ^ г [Я (г _ М + Н (г _ ^ _ ^)] + и(г, г, *),
2езз дг
Ж (г, г,*) = Ж1(*) + Ао Сп£пг2^(г, *) + Цг,г,*), (14)
р(г, г, *) = Леззе15г2^(г, *) + ф(г, г, *), Ао = [2 (С55вц + е25)]-1.
Здесь Н( ■ ) —единичная функция Хэвисайда [22]; Ж1(*), ^(г, *), ^1(г,*) — неизвестные функции, определяемые в процессе решения задачи, с помощью которых удается удовлетворить все условия (3), (4), (8).
В результате подстановки (14) в (1)-(13) получаем новую краевую задачу относительно функций и, -ш, ф. При этом дифференциальные уравнения (1), граничные условия по переменной г (5)-(10) становятся неоднородными с правыми частям Е1, Л2, Е3, В1, В2,..., В1о, начальные условия ио, ио, Жо, И^о следует заменить на ио, йо, -шо, гоо, а краевые условия (2)-(4) принимают следующий вид:
жесткое закрепление:
U (l,z,t) = 0, шарнирное закрепление:
dW
dr
= 0, TT z=i or
= 0;
z=l
W (l,z,t)=0, ф(1^,^ = 0, VU |r=l = 0. (15
Последнее условие (15) получается при = С^-
Применим к краевой задаче относительно функций и, w, ф конечные интегральные преобразования Фурье—Бесселя, используя следующие трансформанты:
ин(jn,z,t)= i u(r,z,t)rJi(jnr)dr, J0
{■WH (jn, z,t), фн (jn ,z,t)} = {w(r,z,t)^(r,z,t)} rJo(jnr)dr
J0
и формулы обращения:
u(r,z,t)=2 ± UHj ^ Ji(jnr), (16)
n=1
jw(r, z, i)Mr, z,t)} = 21 j).
n=0 \J'n>
где jn — положительные нули функций Ji(jn); Jo(jn) соответственно при жестком и шарнирном закреплении пластины, расположенные в порядке их возрастания, n = 0,1, 2,...,jo = 0; S(jn) = Jo(jn) при жестком и S(jn) = Ji(jn) при шарнирном закреплении.
В результате получаем новую задачу относительно трансформант Фурье—Бесселя ин, wh, фн, которую решаем, используя обобщенный метод конечных интегральных преобразований [21] по аксиальной координате. При этом предварительно проводится также процедура стандартизации, связанная с приведением неоднородных граничных условий по переменной z к однородным. Для этой цели вводятся новые функции Uh , Wh , фн, связанные с ин, wh , фн соотношениями
{ин (jnz, t),wH (jnz, t),фн (jn z, t)} =
= {Fl(jnz, t), P{(jnz, t), Рз(jnz, t) } +
+ {Uh(jnz, t), Wh(jnz,t),<pH(jnz,t)}. (17)
В результате подстановки (17) в расчетные соотношения краевой задачи относительно ин, wh , фн получаем новую задачу для функций Uh , Wh , фн. При этом граничные условия с помощью выбора функций Fi, Р2, Р3 становятся однородными.
Начально-краевую задачу относительно функций Uh, Wh, фн решаем, используя обобщенный метод конечных интегральных преобразований (КИП)
[21] с использованием неизвестной трансформанты С(Лгга,п, *) и компонент К1(Лг„,г), К2(Л^га,г), К3(Л^га,г) вектор-функции ядра преобразования (Л^га — положительные параметры, образующие счетное множество, г = 1, 2, 3,...)
С(Лт, п, *) = I (инКИга + Кмга)Жг, (18)
о
и формул обращения
те
{ия, Жн} = ^ Сгга{К1гга,К2г„,Кзг„}укггау-2, (19)
г=1
2^2 2 2 о
Общие решения системы дифференциальных уравнений для электроупругого (8 = 1) и упругого (8 = 2) слоев представлены в работах авторов [16,23].
Подстановка расчетных соотношений К1(Л^га,г), К2(Лг„,г), К3(Л^га,г) в граничные условия по аксиальной координате формирует однородную систему уравнений относительно постоянных интегрирования ..., ^1бт. Разыскивая ее нетривиальное решение, получаем трансцендентное уравнение для вычисления собственных значений Л^га, а также выражения для постоянных , , . . . , ^1бт.
Окончательные выражения функций и (г, г,*), Ж (г, г,*), ^(г, г,*) получим, применяя к трансформанте (18) последовательно формулы обращения КИП (19) и Фурье—Бесселя (16). В результате с учетом (14), (17) имеем
U(r, z, t) = -^ d^i(z,t) r [Я(z - hi) + H(z - hi - h2)] + 2e33 dz
+ 2 V J^rri Pi + V GiraKiirayKiray-2
n=! S(jn) V i=i
W (r,z,t) = Wi(t) + AoCiieiir2Ni(z,i) +
+ 2 V jMfP2 + V С*,*ип||Кп|Г2 ), (20)
n=0 S(jn) V i=i
<p(r,z,t) = Aoe33ei5r2Ni(z,t)+
+ 2 V ( P3 + V ^ггаКзггаУКгпУ-2
n=0 S(jn) V i=i
Заключительным этапом исследования является определение зависимостей Wi(t), Ni(z,t), ^i(z,i). При решении задачи обратного пьезоэффекта первоначально решается задача в случае действия только электрической нагрузки V(r,t) (Ni = 0).
В результате определяется функция Wi(t) из условия отсутствия вертикальных перемещений цилиндрической поверхности пластины при z = h/2, а также характер изменения касательных напряжений arz(1,z,t), который определяет в дальнейшем выражение для Ni(z,t). На втором этапе решения рассматривается действие функции Ni(z,t), которая аппроксимируется параболической зависимостью и определяется из условия равенства нулю перемещений на цилиндрической поверхности пластины в случае действия V(r, t) и Ni (z,t). Суммирование двух расчетов позволяет получить окончательное решение.
Аналогичным образом решается задача прямого пьезоэффекта. Первоначально рассматривается действие механической нагрузки q(r, t), а на следующем этапе функция ф1(г, t) определяется из условия равенства нулю радиальной компоненты вектора индукции электрического поля на цилиндрической поверхности при суммировании двух расчетов.
Полученные расчетные соотношения (20) удовлетворяют дифференциальные уравнения (1) и краевые условия (2)-(11), т. е. представляют замкнутые решения рассматриваемых начально-краевых задач электроупругости для биморфных пьезокерамических пластин.
Заключение. В настоящей работе на основании разработанного в работах авторов [16, 18, 19] алгоритма получены расчетные соотношения для решения частных задач. Представленная методика расчета позволяет получать замкнутые решения связанных осесимметричных нестационарных задач в пространственной постановке для круглых сплошных биморфных пластин на шарнирном и жестком закреплении по контуру.
Допущение об установившемся режиме вынужденных колебаний, используемое при исследовании динамических задач, справедливо только в случае, когда частоты внешнего гармонического воздействия существенно меньше первой частоты собственных колебаний.
Полученные расчетные соотношения позволяют определить напряженно-деформированное состояние, спектр частот собственных колебаний и характер распределения электрического поля в многослойных пьезокерамических преобразователях резонансного и нерезонансного классов в случае действия нестационарной нагрузки.
Конкурирующие интересы. Заявляем, что в отношении авторства и публикации этой статьи конфликта интересов не имеем.
Авторский вклад и ответственность. Все авторы принимали участие в разработке концепции статьи и в написании рукописи. Авторы несут полную ответственность за предоставление окончательной рукописи в печать. Окончательная версия рукописи была одобрена всеми авторами.
Финансирование. Исследование выполнялось без финансирования.
Библиографический список
1. Sharapov V. Bimorph and trimorph piezoelements / Piezoceramic Sensors. Microtech-nology and MEMS. Berlin, Heidelberg: Springer, 2011. pp. 179-229. doi: 10.1007/ 978-3-642-15311-2_6.
2. Sharapov V., Sotula Z., Kunickaya L. Devices to control and diagnose bimorph piezoelements / Piezo-Electric Electro-Acoustic Transducers. Microtechnology and MEMS. Cham: Springer, 2013. pp. 191-212. doi: 10.1007/978-3-319-01198-1_11.
3. Seeley C. E., Delgado E., Kunzman J., Bellamy D. Miniature piezo composite bimorph actuator for elevated temperature operation / International Mechanical Engineering Congress
and Exposition. vol.10 (Seattle, Washington, USA, November 11-15, 2007), Mechanics of Solids and Structures, Parts A and B, 2007. pp. 405-415. doi: 10.1115/imece2007-44088.
4. Гринченко В. Т., Улитко А. Ф., Шульга Н. А. Механика связанных полей в элементах конструкций. Киев: Наук. думка, 1989. 279 с.
5. Smits J. G., Dalke S. I., Cooney T. K. The constituent equations of piezoelectric bi-morphs// Sensors and Actuators A: Physical, 1991. vol. 28, no. 1. pp. 41-61. doi: 10.1016/ 0924-4247(91)80007-C.
6. Ватульян А. О., Рынкова А. А. Об одной модели изгибных колебаний пьезоэлектрических биморфов с разрезными электродами и ее приложениях // Изв. РАН. МТТ, 2007. №4. С. 114-122.
7. Tsaplev V., Konovalov R., Abbakumov K. Disk bimorph-type piezoelectric energy harvester// Journal of Power and Energy Engineering, 2015. vol.3, no. 4. pp. 63-68. doi: 10. 4236/jpee.2015.34010.
8. Цаплев В. М., Коновалов Р. С. Частотные зависимости констант высших порядков пьезоэлектрической керамики// Дефектоскопия, 2017. №7. С. 14-22.
9. Янчевский И. В. Минимизация прогибов электроупругой биморфной пластины при импульсном нагружении // Проблемы вычислительной механики и прочности конструкций, 2011. №16.
10. Шляхин Д. А. Вынужденные осесимметричные колебания пьезокерамической тонкой биморфной пластины// Изв. РАН. МТТ, 2013. №2. С. 77-85.
11. Rashidifar M. A., Rashidifar A. A. Vibrations analysis of circular plate with piezoelectric actuator using thin plate theory and Bessel function // American Journal of Engineering, Technology and Society, 2015. vol. 2, no. 6. pp. 140-156.
12. Jam J. E., Khosravi M., Namdara N. An exact solution of mechanical buckling for functionally graded material bimorph circular plates// Metall. Mater. Eng., 2013. vol. 19, no. 1. pp. 45-63.
13. Jandaghiana A. A., Jafarib A. A., Rahmania O. Vibrational response of functionally graded circular plate integrated with piezoelectric layers: An exact solution // Engineering Solid Mechanics, 2014. vol.2, no. 2. pp. 120-130. doi: 10.5267/j.esm.2014.1.004.
14. Jurenas V., Bansevicius R., Navickaite S. Piezoelectric bimorphs for laser shutter systems: optimization of dynamic characteristics // Mechanika, 2010. no. 5(85). pp. 44-47.
15. Сеницкий Ю. Э., Шляхин Д. А. Нестационарная осесимметричная задача электроупругости для толстой круглой анизотропной пьезокерамической пластины // Изв. РАН. МТТ, 1999. №1. С. 78-87.
16. Шляхин Д. А. Вынужденные осесимметричные колебания толстой круглой жестко закрепленной пьезокерамической пластины// Изв. РАН. МТТ, 2014. №4. С. 90-100.
17. Шульга Н. А., Болкисев А. М. Колебания пьезоэлектрических тел. Киев: Наук. думка, 1990. 228 с.
18. 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. pp. 662-666. doi: 10.1016/j.proeng.2016.08.219.
19. Шляхин Д. А. Динамическая осесимметричная задача прямого пьезоэффекта для круглой биморфной пластины // Вестник Пермского национального исследовательского политехнического университета. Механика, 2017. №1. С. 164-180. doi: 10.15593/ perm.mech/2017.1.10.
20. Sneddon I. N. Fourier Transforms. New York: McGraw-Hill Book Company, Inc., 1950.
21. Сеницкий Ю. Э. Исследование упругого деформирования элементов конструкций при динамических воздействиях методом конечных интегральных преобразований. Саратов: Сарат. ун-т, 1985. 174 с.
22. Владимиров В. С. Обобщенные функции в математической физике. М.: Наука, 1978. 318 с.
23. Шляхин Д. А. Вынужденные осесимметричные изгибные колебания толстой круглой жестко закрепленной пластины// Вестн. СамГУ. Естественнонаучн. сер., 2011. №8(89). С. 142-152.
Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki
[J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2017, vol. 21, no. 4, pp. 773-785
d http://doi.org/10.14498/vsgtu1564
ISSN: 2310-7081 (online), 1991-8615 (print)
MSC: 74F, 74K20
Forced axisymmetric oscillations of circular multilayer bimorph plates
D. A. Shlyakhin, O. V. Ratmanova
Samara State Technical University,
244, Molodogvardeyskaya st., Samara, 443100, Russian Federation.
Abstract
A method for calculating circular multilayer bimorph plates is presented and new analytical solutions of axisymmetric dynamic problems of direct and inverse piezoelectric effects are obtained. The cases of hinged and rigid fixing of the outer contour of the structure are considered. To investigate related linear problems, a mathematical apparatus is used in the form of a method of finite integral transformations. The constructed calculated relationships allow us to substantiate the constructive solutions of piezoceramic transducers.
Built design ratio allows to prove the constructive solutions of multilayer piezoelectric ceramic transducers, namely, to choose the geometrical dimensions and physical characteristics of the materials used, define the dimensions of a split circular electrode, allowing most effectively to convert the external electrical stimulation into mechanical vibrations at various frequency. In addition, it is possible to perform stress-strain state, the nature of the change of the electric field and frequency range of the axisymmetric vibrations of the considered systems.
Keywords: bimorph plate, electroelasticity, nonstationary loading, integral transformations.
Received: 14th October, 2017 / Revised: 9th December, 2017 / Accepted: 18th December, 2017 / First online: 29th December, 2017
Competing interests. We declare that we have no conflicts of interests with the authorship and publication of this article.
Authors' contributions and responsibilities. Each author has participated in the article concept development and in the manuscript writing. The authors are absolutely
Short Communication
Q ©® The content is published under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/) Please cite this article in press as:
Shlyakhin D. A., Ratmanova O. V. Forced axisymmetric oscillations of circular multilayer bimorph plates, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2017, vol. 21, no. 4, pp. 773-785. doi: 10.14498/vsgtu1564 (In Russian). Authors' Details:
Dmitry A. Shlyakhin http://orcid.org/0000-0003-0926-7388
Dr. Techn. Sci.; Head of Dept.; Dept. of Strength of Materials and Structural Mechanics;
e-mail: d- 612-mit2009@yandex. ru
Olesya V. Ratmanova © http://orcid.org/0000-0001-7591-4921 Assistant; Dept. of Strength of Materials and Structural Mechanics; e-mail: [email protected]
responsible for submitting the final manuscript in print. Each author has approved the
final version of manuscript.
Funding. The research has not had any funding.
References
1. Sharapov V. Bimorph and trimorph piezoelements, In: Piezoceramic Sensors. Microtech-nology and MEMS. Berlin, Heidelberg, Springer, 2011, pp. 179-229. doi: 10.1007/ 978-3-642-15311-2_6.
2. Sharapov V., Sotula Z., Kunickaya L. Devices to control and diagnose bimorph piezoelements, In: Piezo-Electric Electro-Acoustic Transducers. Microtechnology and MEMS. Cham, Springer, 2013, pp. 191-212. doi: 10.1007/978-3-319-01198-1_11.
3. Seeley C. E., Delgado E., Kunzman J., Bellamy D. Miniature piezo composite bimorph actuator for elevated temperature operation, In: International Mechanical Engineering Congress and Exposition, vol.10 (Seattle, Washington, USA, November 11-15, 2007), Mechanics of Solids and Structures, Parts A and B, 2007, pp. 405-415. doi: 10.1115/imece2007-44088.
4. Grinchenko V. T., Ulitko A. F., Shul'ga N. A. Mekhanika sviazannykh polei v elementakh konstruktsii [Mechanics of coupled fields in structural elements]. Kiev, Nauk. dumka, 1989, 279 pp. (In Russian)
5. Smits J. G., Dalke S. I., Cooney T. K. The constituent equations of piezoelectric bi-morphs, Sensors and Actuators A: Physical, 1991, vol.28, no. 1, pp. 41-61. doi: 10.1016/ 0924-4247(91)80007-C.
6. Vatul'yan A. O., Rynkova A. A A model of bending vibrations of piezoelectric bimorphs with split electrodes and its applications, Mech. Solids, 2007, vol.42, no. 4, pp. 595-602. doi:10.3103/S0025654407040127.
7. Tsaplev V., Konovalov R., Abbakumov K. Disk bimorph-type piezoelectric energy harvester, Journal of Power and Energy Engineering, 2015, vol. 3, no. 4, pp. 63-68. doi: 10.4236/jpee. 2015.34010.
8. Tsaplev V. M., Konovalov R. S. Frequency dependences of higher-order constants of piezoelectric ceramics, Russ. J. Nondestruct. Test, 2017, vol.53, no. 7, pp. 485-492. doi: 10. 1134/S1061830917070087.
9. Yanchevskiy I. V. Minimizing deflections of round electroelastic bimorph plate under impulsive loading, Problems of computational mechanics and strength of structures, 2011, no. 16 (In Russian).
10. Shlyakhin D. A. Forced nonstationary axisymmetric vibrations of a piezoceramic thin bimorph plate, Mech. Solids, 2013, vol.48, no. 2, pp. 178-185. doi: 10.3103/ S002565441302009X.
11. Rashidifar M. A., Rashidifar A. A. Vibrations analysis of circular plate with piezoelectric actuator using thin plate theory and Bessel function, American Journal of Engineering, Technology and Society, 2015, vol. 2, no. 6, pp. 140-156.
12. Jam J. E., Khosravi M., Namdara N. An exact solution of mechanical buckling for functionally graded material bimorph circular plates, Metall. Mater. Eng., 2013, vol. 19, no. 1, pp. 45-63.
13. Jandaghiana A. A., Jafarib A. A., Rahmania O. Vibrational response of functionally graded circular plate integrated with piezoelectric layers: An exact solution, Engineering Solid Mechanics, 2014, vol.2, no. 2, pp. 120-130. doi: 10.5267/j.esm.2014.1.004.
14. Jurenas V., Bansevicius R., Navickaite S. Piezoelectric bimorphs for laser shutter systems: optimization of dynamic characteristics, Mechanika, 2010, no. 5(85), pp. 44-47.
15. Senitskii Yu. E., Shlyakhin D. A. The nonstationary axisymmetric problem of electroelas-ticity for a thick circular anisotropic piezoceramic plate, Mech. Solids, 1999, vol.34, no. 1, pp. 66-74.
16. Shlyakhin D. A. Forced axisymmetric vibrations of a thick circular rigidly fixed piezoceramic plate, Mech. Solids, 2014, vol.49, no. 4, pp. 435-444. doi: 10.3103/S0025654414040086.
17. Shul'ga N. A., Bolkisev A. M. Kolebaniia p'ezoelektricheskikh tel [Fluctuations of Piezoelectric Bodies]. Kiev, Nauk. dumka, 1990, 228 pp. (In Russian)
18. 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, pp. 662666. doi: 10.1016/j.proeng.2016.08.219.
19. Shlyakhin D. A. Dynamic axisymmetric problem direct piezoeffect for round bimorph plate, PNRPU Mechanics Bulletin, 2017, no. 1, pp. 164-180 (In Russian). doi: 10.15593/perm. mech/2017.1.10.
20. Sneddon I. N. Fourier Transforms. New York, McGraw-Hill Book Company, Inc., 1950.
21. Senitskii Yu. E. Issledovanie uprugogo deformirovaniia elementov konstruktsii pri dinamich-eskikh vozdeistviiakh metodom konechnykh integral'nykh preobrazovanii [Investigation of the Elastic Deformation of Construction Elements under Dynamic Actions by the Method of Finite Integral Transforms]. Saratov, Saratov Univ., 1985, 174 pp. (In Russian)
22. Vladimirov V. S. Obobshchennye funktsii v matematicheskoi fizike [Generalized Functions in Mathematical Physics]. Moscow, Nauka, 1978, 318 pp. (In Russian)
23. Shlyakhin D. A. The compelled axisymmetric bending fluctuations of thick round rigid plate, Vestnik Samarskogo Gosudarstvennogo Universiteta. Estestvenno-Nauchnaya Seriya, 2011, no. 8(89), pp. 142-152 (In Russian).
