НЕСТАЦИОНАРНАЯ ЗАДАЧА ОБРАТНОГО ПЬЕЗОЭФФЕКТА ДЛЯ ДЛИННОГО ПЬЕЗОКЕРАМИЧЕСКОГО ТЕРМОУПРУГОГО ЦИЛИНДРА

Для цитирования М.А. Кальмова. Нестационарная задача обратного пьезоэффекта для длинного пьезокерами-ческого термоупругого цилиндра. Вестник Дагестанского государственного технического университета. Технические науки. 2020;47 (4):57-68. D0I:10.21822/2073-6185-2020-47-4-57-68

For citation: M.A. Kalmova. Dynamic inverse piezo-effect problem for a long piezoceramic thermoelastic cylinder. Herald of Daghestan State Technical University. Technical Sciences. 2020;47(4):57-68.(In Russ.) D0I:10.21822/2073-6185-2020-47-4-57-68

ЭНЕРГЕТИЧЕСКОЕ, МЕТАЛЛУРГИЧЕСКОЕ И ХИМИЧЕСКОЕ МАШИНОСТРОЕНИЕ POWER, METALLURGICAL AND CHEMICAL MECHANICAL ENGINEERING

УДК 539.3

DOI: 10.21822/2073-6185-2020-47-4-57-68

НЕСТАЦИОНАРНАЯ ЗАДАЧА ОБРАТНОГО ПЬЕЗОЭФФЕКТА ДЛЯ ДЛИННОГО ПЬЕЗОКЕРАМИЧЕСКОГО ТЕРМОУПРУГОГО ЦИЛИНДРА

М.А. Кальмова

Самарский государственный технический университет, 443001 г. Самара, ул. Молодогвардейская 244, Россия

Резюме. Цель. Целью работы является решение несвязанной нестационарной задачи термоэлектроупругости для длинного полого пьезокерамического цилиндра при действии на его поверхностях электрической нагрузки в виде разности потенциалов. Метод. Математическая формулировка рассматриваемой задачи термоэлектроупругости включает систему несамосопряженных дифференциальных уравнений. Для ее решения на первом этапе рассматривается связанная задача обратного пьезоэффекта без учета влияния температурного поля, а на следующем - исследуется гиперболическая задача теплопроводности (LS-теория) при заданном (определенном) электроупругом поле. Результат. Построено новое замкнутое решение динамической задачи обратного пьезоэффекта для длинного пьезокерамического термоупругого цилиндра. Рассматривается случай действия на его лицевых поверхностях нестационарной электрической нагрузки в виде разности потенциалов. Заданы температура окружающей среды и закон конвекционного теплообмена (граничное условие 3- рода). Расчетные соотношения, полученные с помощью обобщенного метода конечных интегральных преобразований, дают возможность определить напряженно-деформированное состояние и термоэлектрические поля, индуцируемые в пьезокерамическом элементе при произвольном электрическом внешнем воздействии. Вывод. Построенное решение позволяет определить напряженно-деформированное состояние и электрическое поле в пьезокерамическом цилиндре, а также проанализировать с помощью гиперболической LS-теории теплопроводности влияние индуцируемого температурного поля на электроупругое состояние рассматриваемой системы. Анализ численных результатов позволяет сделать вывод о незначительных потерях энергии, связанных с нагревом электроупругой системы. Разработанный алгоритм расчета находит свое применение при проектировании нерезонансных и резонансных пьезоэлектрических измерительных приборов.

Ключевые слова: длинный пьезокерамический цилиндр, задача термоэлектроупругости, нестационарное воздействие, конечные интегральные преобразования

DYNAMIC INVERSE PIEZO-EFFECT PROBLEM FOR A LONG PIEZOCERAMIC

THERMOELASTIC CYLINDER M.A. Kalmova

Samara State Technical University, 244 Molodogvardeyskaya St., Samara 443001, Russia

Abstract. Objective. The objective of this work is to solve an unrelated dynamic problem of thermoelectroelasticity for a long hollow piezoceramic cylinder under the action of an electric load on its surfaces in the form of a potential difference. Methods. The mathematical formulation of the con-

sidered problem of thermoelectroelasticity includes a system of non-selfadjoint differential equations. At the first stage, the authors consider the associated inverse piezoelectric effect problem without taking into account the influence of the temperature field, and at the next stage, study the hyperbolic heat conduction problem (Lord-Shulman theory) for a given (defined) electroelastic field. Result. A new closed solution to the dynamic inverse piezoelectric effect problem for a long piezoceramic thermoe-lastic cylinder is constructed. The case of the action of a dynamic electric load in the form of a potential difference on its front surfaces is considered. The ambient temperature and the law of convection heat transfer (3-kind boundary condition) are set. The calculated relations obtained using the generalized method of finite integral transformations allow determining the stress-strain state and thermoelectric fields induced in a piezoceramic element under an arbitrary electrical external influence. Conclusion. The constructed solution allows determining the stress-strain state and electric field in a pie-zoceramic cylinder, as well as analyzing the effect of the induced temperature field on the electroelas-tic state of the system under consideration using the hyperbolic Lord-Shulman theory of thermal conductivity. Analysis of the numerical results allows concluding that there are insignificant energy losses associated with heating the electroelastic system. The developed calculation algorithm is used in the design of non-resonant and resonant piezoelectric measuring devices.

Keywords: long piezoceramic cylinder, thermoelasticity problem, dynamic action, finite integral transformations

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

Математическая формулировка начально-краевых задач термоэлектроупругости включает систему несамосопряженных дифференциальных уравнений, исследование которых, как правило, проводится при использовании численных методов [7-9]. Однако достаточно слабые эффекты взаимодействия полей различной физической природы удается проанализировать только с помощью замкнутых аналитических решений. Проблема интегрирования исходных расчетных соотношений и построение общего решения приводит к проведению расчетов в упрощенной постановке, а именно, исследуются несвязанные задачи [10,11] или анализируются бесконечно длинные тела [12-14]. При этом в большинстве работ, в качестве внешнего воздействия рассматривается нестационарная температурная нагрузка. В свою очередь, при решении задач прямого и обратного пьезоэффектов, электроупругая конструкция рассматривается, как консервативная система без потери энергии на ее нагрев [15-18].

Постановка задачи. Целью настоящей работы является решение несвязанной нестационарной задачи термоэлектроупругости для длинного полого пьезокерамического цилиндра при действии на его поверхностях электрической нагрузки в виде разности потенциалов. Построенное решение позволяет определить напряженно-деформированное состояние и электрическое поле в пьезокерамическом цилиндре, а также проанализировать с помощью гиперболической LS-теории теплопроводности влияние индуцируемого температурного поля на электроупругое состояние рассматриваемой системы.

Методы исследования. Пусть полый длинный незакрепленный пьезокерамический цилиндр занимает в цилиндрической системе координат (г,, 6, z) область Q:

{a < < b, 0 <6 < 2ж, -да < z < да} . Рассматривается случай действия на лицевых электродиро-

ванных поверхностях нестационарной нагрузки в виде разности потенциалов 2V* (t,) , явля-

ющейся произвольной функцией времени t*. При этом известны температура окружающей

среды 3* и закон конвекционного теплообмена (граничное условие 3- рода).

В общем случае дифференциальные уравнения термоэлектроупругости и краевые условия рассматриваемой задачи в безразмерной форме записываются следующим образом [19, 20]:

V

ди

и

дф 1 дф

- пх— + V--п2

дг r дг r дг

-V0+п--

0 д2и

dtг

= 0

-vd п V дг

vd-

дг

ди

дг

д ад2' — + ß—т

дt дt2

+ п

1 ди

г дг

+ n V0 = 0

П0 + пVU - п9 —

дг j

= 0

г = R,1

ди

дг

+ п10

дф

г дг

0 = 0 , ф = +у

д0 дг

+ х0 = пх 3 ;

д

t = 0 и = ф = 0 = 0 , — {U, ф, 0} = {U0, фо, 00}; (3) где { U, U0, r, R} = { U*, Uo, r*, a}/b, {0, 0o, 0} = —, {0*, 0o, (0 - To}

^33

{ф, фо, V} = , {ф\ Фо, V* }, t = | Cl,

^33° b V p

C

n

11

C

n

33

'31

"33

n3 =

Y11

Yi3

e

n4 =

33

C s

33 33

П5 =

e e

C s

33 33

n

g3e:

3^33

S33Y33

n7 = ¿ Cl ,

7 M p

n8 =

bYi23T 0

A^cii

p

n9 =

bY33 g,T0 C:

33

C

e33A

p

n

13

b

10

C

n

11

33

a— ß = ßß-k b \

C2i , V = ^ +1, p дг г

и (г;, 4) - радиальная составляющая вектора перемещений; ф (г;, — потенциал электрического поля;

©;= Т — Т0 ; 0;, Т, Т0 — приращение и текущая температуры, а также температура первоначального состояния тела, при котором отсутствуют механические напряжения; Ст5, , ¿Г33 — модули упругости, пьезомодули и коэффициент диэлектрической проницаемости электроупругого материала (т, 5 = 1,3^; А, к — коэффициенты теплопроводности и объемной теплоемкости материала; /11, /33 — компоненты тензора температурных напряжений; ^3 — компонента тензора пирокоэффициентов; — время релаксации;

и — коэффициент теплоотдачи; и°, фд, 0 известные в начальный момент времени скорости перемещения, потенциала и температуры.

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

Построение общего решения задачи электроупругости. Математическая формулировка рассматриваемой задачи на основании расчетных соотношений (1)-(3) имеет вид:

2

^dU U ^дф 1 дф d2U л

V--n — + V — - n2—---- = 0

dr r dr r dr dt

1 dU

r dr

r = R,1 t = 0

^ дф ^dU

+ n4V--h n

dr dr

dU U дф

— + nio — + ^ = 0 • ф = ±у ;

dr

= 0

г дг ди

и = 0 , — = Ис;

дЬ

На первом этапе решения выполняется процедура приведения расчетных соотношений (4) - (6) к виду, позволяющему в дальнейшем использовать метод конечных интегральных преобразований [21]. Для этого вводятся новые функции и(г,t), %(г,t) , связанные с и (г, t) ,

ф( т, t) следующими соотношениями:

и

где {И, (г, t), И2 (г, t)} = {Л (г),Л (г)}У (t).

r = R,1

Подстановка (7) в (4) - (6), при удовлетворении условий

df

+ П10 ^ + = 0 , f2 (r) = ±1;

dr r dr 2

позволяет получить краевую задачу относительно функций и (г,t) , %(г,t) с однородными граничными условиями:

^ du u ^dr 1 dr d2u

V--nx — + V^-n2—^--- = R

dr r dr r dr dt

1•

„^ dr ^ d 1 du ^ V— + n.V— + n--= R •

/-s 4 /-s 5 ^ 2 1

dr dr r dr

r = R,1

t=0

du d r fn10 u + dr r dr = 0 , X = 0 ;

dU Ht = Uo - дНг dt ;

dH —1 + n dr H _ 1 r2 Y7dH 2 -V—2 + dr 1 dH? n2--2 r dr d2 H + d • R =VdH 2 dr оr 1 dH - n5--1 r dr

Начально-краевую задачу (9)-(11) решаем, используя структурный алгоритм вырожденного обобщенного конечного интегрального преобразования (КИП) [21]. Для этого вводим на сегменте [К, 1] КИП с неизвестными компонентами вектор-функции ядра преобразований

К (Л, г), К2 (Л, г) :

1

О

{и (г, г), %{г, 0} = £ О (Л, 0{К (Л, г), К 2 (Л, г)}КЦ2, 11К1Г = п4 / К1Л г)2г^г ,

1=1 Я

где Л ~ собственные значения.

В результате использования алгоритма КИП получаем задачи для трансформанты

О (Л, О :

R

t = 0

Gi\t=o = G0

= RH

dt2 г г н dG dGn

dt \t=o dt

и компонент ядра преобразований Kl ( Д, r), K2 (Д, r ) :

+ Д2 к = o

^dK K ^dK 1 dK „2.

V—1 -щ-f + V—2-щ--2 '

dr r dr r dr

r = R,1

_ dK ^ dK 1 dK r,

V-2 + П4 V-1 + Щ--1 = 0 ;

dr dr r dr

dKi K1 >dK2_o K = o •

dr

+ Що — +

r dr

где Rh = -j(R\K\ + RK)rdr, Go = -jH^K.rdr, ^ = £ (щ - | t=o Krdr

Общий интеграл уравнения (13), с учетом начальных условий (14), имеет вид:

G(Л,t) = G0 cos A,t + —0-+ Л-1 JRH (t) sin Л (t-t)dr .

dt Л r.

Система (15) приводится к следующему разрешающему уравнению относительно функции К (Л, г) :

d 2K, 1 dK

г

dr2

+—

r dr

+

.2 Л

A2 - —

Ai 2 r2

V ' У

K = Щ12D

r2

решение которого, имеет вид

К

Здесь Qv,—1 (...), Jv (...),(•••)— неэлементарные функции Ломмеля, обыкновенные функции Бесселя 1-го и 2-го родов порядка V [22] ; —1г ...—Зг — постоянные интегрирования;

2 Д

A =

2 щ + п2п5

v =

П12 =

п

1 + щ ' 1 + щ 1 + щ

Интегрирование второго уравнения

Подстановка К (Л, г), К2 (Л, г) в граничные условия (16) формирует систему алгебраических уравнений, решение которой позволяет определить постоянные интегрирования —i и собственные значения Л .

Окончательные выражения функций и (г, t), ф( г, ^) получим, применяя к трансформанте (17) формулы обращения (12). В результате, с учетом (7), имеем:

да

{и

г=1

Функции (г), /2 (г) определяются при решении следующих дифференциальных уравнений

v ш _ n,m+v ш _ „ I m=о

dr

V

df2

_ nV

dr

r dr

dfi

— n

0

^ <3г г с1г

и удовлетворение граничных условий (8), что позволяет существенно упростить правые части расчетных соотношений (9).

Построение общего решения задачи теплопроводности. Математическая формулировка данной задачи записывается следующим образом:

V--пп

Г d „d2^

dr

^r + ß-

ydt dt j

0 = R

-з ;

r = R,1 t = 0 где R3 =

d0 гл Ü

— + Пц0 = nii_& ;

dr

dQ _ 0 = 0 , Tt = ;

d+ßd

2 Л

dt dt1

rs 1

n8VU _ П9 dT

dr

\

j

Процедура стандартизации граничных условий (24) выполняется с помощью замены ©(г, t) новой функцией Ь (г, t) с помощью такого разложения:

0

Подстановка (26) в (23) - (25) позволяет получить краевую задачу относительно функций Ь (г, t) с однородными граничными условиями:

XI dL

V--n

dr

dL

Г d „d2^

-r + ß-

у dt dt j

L = R

-3 ;

r = R,1 — + nL = 0

dr

ii-*

dL

х = 0 Ь = -3 , — = ;

Использование обобщенного конечного интегрального преобразования (КИП) [21] путем введения неизвестной функции ядра преобразований N (д, г) :

Y (д. t ) = J L (r, t) N (д. r )rdr,

r

1

L(r,t) = jj Y(д,t)N(д,r)\\N,l2, ||N| = J N(A,r)2rdr

i=1 R

позволяет получить задачи для трансформанты Y (д, t) :

dt+ßdt2

2 Л

Y +—Y. = F

n

h .

t = 0 Y) = _^f Nrdr , \leQrNdr;

J dt JR 0

R

и функции ядра преобразований

d2N 1 dN 2лт n

—T + - — + A2N = 0,

dr r dr

dN

r = R,1 -+ nnN = 0;

dr

1 f

где FH =--I RrNdr, Д - собственные значения.

п ^

7 R

Общий интеграл уравнения (31), с учетом начальных условий (32), имеет вид:

Я

I ¿JY t I

Y

0 J

а решение равенства (33) записывается следующим образом:

N

-P 1 ±4P~2 - 4w2 P

ГДе Pl,2 =-"--- •

Подстановка (36) в граничные условия (34) формирует выражения для постоянных интегрирования Е1г , Е2г:

Ei = niiYo (wr) - WYi (Wr), E2, = WJi (wr)- niiJo (WR),

и трансцендентное уравнение для определения W :

ЕЬ [niiJ0 (w ) - WJi (w )] + ЕЪ [niiY0 (W ) - M (W )] = 0 •

Окончательные выражения функции ©

формулу обращения (30).

В результате, с учетом (26), имеем:

да

©

i=i

Численный анализ результатов. В качестве примера рассматривается радиально поляризованный пьезокерамический цилиндр (R = 0.8, b = 0.02 м) состава PZT-4, имеющего следующие физические характеристики: {e3i,e33} = {-5.2, i5.i} Кл/м2 ,

£33 = 5.62xi0-9ф/м, {Cii,С33,Ci3}= {i3.9, ii.5, 7.43,}xi0i0 H/м2 , p = 7500 кг/м3 ,

У,/зз} = {4.6, 3.9}xi05 Н/(м2 К), g3 = 2xi0-4 Кл/(м2 К), k = 3xi06 Дж/(м3 К ), Л = i.6

2

Вт/(м К), а= 5.6 Вт/(м 2 К) [16].

Рассматривается случай действия гармонической электрической нагрузки на поверхностях цилиндра:

V (t ) = V sin 0t,

где V, 6 - амплитудное значение потенциала и частота вынужденных колебаний в безразмерной форме.

Обсуждение результатов. На рис. 1,2 представлены графики изменения перемещений U(r, t) по времени t, а также потенциала ф(r,t) (кривая 1) и напряженности Er (r, t) (2) электрического поля по радиальной координате r (6 = 0.9\). Пунктирной линией обозначена осциллограмма внешней электрической нагрузки.

Анализ представленных графиков позволяет сделать следующие выводы:

1. При высокочастотном внешнем воздействии (0 = 0.9\ ) вследствие наложения волн

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

2. Использование гипотез о линейной зависимости потенциала ф( r, t ) и постоянном зна-

дф

чении напряженности Er ( r, t ) =--—- электрического поля при исследовании задач

дг

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

и (1, t ) / V

200

- 200

0

t

0 10 20 30 40 50

Рис.1. График изменения перемещений U (l, tj по времени t Fig.1. Time displacement graph

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

{ф( г, г), вг (г, г)/5} / V

1

Рис.2. Графики изменения амплитудных значений потенциала r, tj (кривая 1) и напряженности

Er (r, t j (2) электрического поля по радиальной координате r

Fig.2. Graphs of changes in the amplitude values of the potential (curve 1) and intensity (2) of the electric

field along the radial coordinate

2

0

r

0

0.9

На рис. 3 представлены графики изменения температуры ©*(1, t ) во времени t. Очевид-

д©*

но, при t > 10 условие --—- = 0 выполняется при 0О = 19/VQ (рис. 3, график 1).

dt

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

©•(l, t) ,0 С

О.О8

О.Об

0.04

0.02

/ 1

—^ т 2

l0

l5

20

Рис.3. Графики изменения температуры 0* (1, t j по времени t

Fig.3. Temperature graphs over time

Принимая T0 = 295 К (22 0 С), 00 = 1.9/V0, V* = 10 B , получаем следующую картину изменения температурного поля по радиальной координате в цилиндре (рис. 4). Цифрами 1,2 обозначены результаты для времени t = 3 0,7 .

0* (r, t jС

О.О7б

0.075

-----

1

2

Рис.4. Графики изменения температуры 0* (i, t) по радиальной координате

r

Был выполнен эксперимент, где измерения температуры проводились с помощью прибора Center 301. Результаты показали, что температура на лицевых поверхностях исследуемого

элемента при внешнем воздействии повысилась на 0.1 0С , что практически повторяет результаты теоретических расчетов (рис.3,4).

Для оценки влияния температуры на электроупругое поле правые части дифференциальных уравнений R, R (9), с учетом равенств (1), записываются следующим образом:

t

О

О

5

r

О.8

0.9

1

0 QH R =V0-n--V—1 + n

H

д r

R = -n6 V0 + V

1 2 r

1 -V

dHn

dH.

2 n4V

д r

dH,

+n

n

1 dH2 d2 h

--2 +-2

r dr dt

1 dHx

дг дт г дт

На рис.5 представлены зависимости "и(1,X) — X" с учетом (сплошная линия) и без учета (пунктирная линия) влияния температуры на радиальные перемещения. Численный анализ результатов показывает, что температурное поле практически не влияет на величину и х) , т.е.

в задачах обратного пьезоэффекта можно пренебречь потерями энергии на нагрев электроупругой системы.

и (1, х) / V

200

- 200

10

20

30

40

50

Рис.5. Графики изменения температуры U (1, t) по времени t (сплошная и пунктирная линии-с учетом и без учета влияния температурного поля) Fig.5. Temperature graphs U (1, t) over time t

Вывод. Следует отметить, что при решении динамических задач прямого и обратного пьезоэффекта пьезокерамические элементы можно рассматривать, как консервативные системы.

Однако при решении нестационарных задач термоэлектроупругости в случае действия температурной нагрузки необходимо исследовать расчетные соотношения с учетом связанности упругих, электрических и тепловых полей [13].

Библиографический список:

Бобцов А.А. Исполнительные устройства и системы для микроперемещений. СПб ГУ ИТМО, 2011.131 с. Джагуров Р.Г. Пьезоэлектронные устройства вычислительной техники, систем контроля и управления. СПб.: Политехника, 1994. 608 с.

Евсейчик Ю.Б., Медведев К.В. Чувствительность гидроакустического датчика давления // Гидравлика и гидротехника. Науч.- техн. сб. Киев: НТУ. 2008. Вып. 62. С. 10-16.

Янчевский И.В. Управление колебаниями изгиба круглого асимметричного биморфного пьезопреобразователя с разрезными электродами // Пробл. машиностроения. 2012. Т.15, № 2. С. 37-43.

Abedi M., Jafari-Talookolaei R., Valvo P. A new solution method for free vibration analysis of rectangular laminated composite plates with general stacking sequences and edge restraints // Computers & Structures. 2016. Vol. 175. рр. 144-156. Berndt E.A., Sevostianov I. Action of a smooth flat charged punch on the piezoelectric half-space possessing symmetry of class // International Journal of Engineering Science. 2016. Vol. 103. рр. 77-96.

Abbas I.A., Youssef H.M. Finite element analysis of two-termoperature generalized magneto-thermoelasticity/ Arch Appl Mech. 2009.79. 917-925.

He T. et al. A generalized electromagneto-thermoelastic problem for an infinitely long solid cylinder/ European Journal of Mechanics A-Solids. 24 (2005). рр.349-359.

Youssef H.M. Theory of two-temperature generalized thermoelasticity/ IMA J.Appl.Math. 71(3) (2006)/383-390. 10. Фирсанов В.В., Нгуен, Ле Хунг. Напряженно-деформированное состояние произвольных оболочек с учетом термоэлектрического воздействия на основе уточненной теории// Тепловые процессы в технике. 2010. №3. С.110-117.

9.

r

0

t

0

11. Abbas I.A., Zenkour A.M. LS model on electro-magneto-thermoelastic response of an infinite functionally graded cylinder/ Composite Structures. 96. (2013) 89-96.

12. Ватульян А.О., Кирютенко А.Ю., Наседкин А.В. Плоские волны и фундаментальные решения в линейной термоэлек-троупругости// ПМТФ. 1996. Т.37. №5. С.135-142.

13. Ватульян А.О., Нестеров С.А. Динамическая задача термоэлектроупругости для функционально-градиентного слоя// Вычислительная механика сплошных сред. 2017. Т.10. №2. С.117-126.

14. Белянкова Т.И., Калинчук В.В. К моделированию преднапряженного термоэлектроупругого полупространства с покрытием // Изв. РАН. МТТ. 2017. № 1. С. 117-135.

15. Ватульян А.О., Рынкова А.А. Об одной модели изгибных колебаний пьезоэлектрических биморфов с разрезными электродами и ее приложениях// Изв. РАН. МТТ. 2007. №4. С. 114-122.

16. Бардзокас Д.И. Математическое моделирование в задачах механики связанных полей. Т.II: Статические и динамические задачи электроупругости для составных многосвязных тел. М.: Комкнига, 2005. 376 с.

17. Шляхин Д.А. Вынужденные осесимметричные колебания пьезокерамической тонкой биморфной пластины // Изв. РАН. МТТ. 2013. №2. С. 77-85.

18. Шляхин Д.А. Вынужденные осесимметричные колебания толстой круглой жестко закрепленной пьезокерамической пластины // Изв. РАН. МТТ. 2014. №4. С. 90-100.

19. Коваленко А.Д. Введение в термоупругость. Киев: Наук. думка, 1965. 204 с.

20. Гринченко В.Т., Улитко А.Ф., Шульга Н.А. Механика связанных полей в элементах конструкций. Киев: Наук. дум-ка,1989. 279 с.

21. Сеницкий Ю.Э. Метод конечных интегральных преобразований - обобщение классической процедуры разложения по собственным вектор-функциям // Изв. Саратовского ун-та. Новая серия. Матем., механ., информатика, 2011. № 3(1). С. 61-89.

22. Камке Э. Справочник по обыкновенным дифференциальным уравнениям. М.: Наука,1965. 703 с.

References:

1. Bobtsov A.A. Ispolnitelnyie ustroystva i sistemyi dlya mikroperemescheniy. SPb GU ITMO, 2011. 131 s. [Bobtsov A.A. Actuators and systems for micromovements. - SPb GU ITMO, 2011. 131 p. (In Russ)]

2. Dzhagurov R.G. Pezoelektronnyie ustroystva vyichislitelnoy tehniki, sistem kon-trolya i upravleniya. SPb.: Politehnika, 1994. 608 s. [Dzhagurov R.G. Piezoelectronic devices of computer technology, control and management systems. SPb .: Polytechnic, 1994. 608 p.(In Russ)]

3. Yevseychik Y.B., Medvedev K.V. Sensitivity of the hydroacoustic pressure sensor // Hydraulics and hydraulic engineering. Scientific and technical Sat. Kiev: NTU. 2008. Issue. 62 . рр. 10-16

4. Yanchevskiy I.V. Upravlenie kolebaniyami izgiba kruglogo asimmetrichnogo bi-morfnogo pezopreobrazovatelya s razreznyi-mi elektrodami // Probl. mashinostroeniya. 2012. T.15, № 2. S. 37-43 [Yanchevsky I.V. Control of bending vibrations of a circular asymmetric bi-morphic piezoelectric transducer with split electrodes // Probl. mechanical engineering. 2012. T.15, No. 2. рр. 37-43 (In Russ)]

5. Abedi M., Jafari-Talookolaei R., Valvo P. A new solution method for free vibration analysis of rectangular laminated composite plates with general stacking sequences and edge restraints // Computers & Structures. 2016. Vol. 175. рр. 144-156.

6. Berndt E.A., Sevostianov I. Action of a smooth flat charged punch on the piezoelectric half-space possessing symmetry of class // International Journal of Engineering Science. 2016. Vol. 103. рр. 77-96.

7. Abbas I.A., Youssef H.M. Finite element analysis of two-termoperature generalized magneto-thermoelasticity/ Arch Appl Mech. 2009.79. 917-925.

8. He T. et al. A generalized electromagneto-thermoelastic problem for an infinitely long solid cylinder / European Journal of Mechanics A-Solids. 24 (2005).рр.349-359.

9. Youssef H.M. Theory of two-temperature generalized thermoelasticity/ IMA J.Appl.Math. 71(3) (2006)/383-390.

10. Firsanov V.V., Nguen, Le Hung. Napryazhenno-deformirovannoe sostoyanie proiz-volnyih obolochek s uchetom termoel-ektricheskogo vozdeystviya na osnove utochnennoy teorii // Teplovyie protsessyi v tehnike. 2010. - №3. рр.110-117.

11. Abbas I.A., Zenkour A.M. LS model on electro-magneto-thermoelastic response of an infinite functionally graded cylinder/ Composite Structures. 96. (2013) 89-96.

12. Vatulyan A.O., Kiryutenko A.Yu., Nasedkin A.V. Ploskie volnyi i fundamentalnyie resheniya v lineynoy termoelektroupru-gosti // PMTF. -1996. -T.37. -№5. - S.135-142 [Vatulyan A.O., Kiryutenko A.Yu., Nasedkin A.V. Plane waves and fundamental solutions in linear thermoelectroelasticity // PMTF. 1996. T.37. №5. рр.135-142 (In Russ)]

13. Vatulyan A.O., Nesterov S.A. Dinamicheskaya zadacha termoelektrouprugosti dlya funktsionalno-gradientnogo sloya// Vyi-chislitelnaya mehanika sploshnyih sred. 2017. T.10. №2. S.117-126 [Vatulyan A.O., Nesterov S.A. Dynamic problem of thermoelectroelasticity for a functional-gradient layer // Computational mechanics of continuous media. 2017. Т.10. №2. pp.117-126 (In Russ)]

14. Belyankova T.I., Kalinchuk V.V. K modelirovaniyu prednapryazhennogo termoelekt-rouprugogo poluprostranstva s pokryit-iem // Izv. RAN. MTT. 2017. № 1. S. 117-135 [Belyankova T.I., Kalinchuk V.V. On modeling a prestressed thermoelectroe-lastic half-space with a coating. Izv. RAS. MTT. 2017. No. 1. рр. 117-135 (In Russ)]

15. Vatulyan A.O., Ryinkova A.A. Ob odnoy modeli izgibnyih kolebaniy pezoelektri-cheskih bimorfov s razreznyimi elektrodami i ee prilozheniyah// Izv. RAN. MTT. 2007. №4. S. 114-122 [Vatulyan A.O., Rynkova A.A. On one model of bending vibrations of piezoelectric bimorphs with split electrodes and its applications // Izv. RAN. MTT. 2007. No. 4. рр. 114-122 (In Russ)]

16. Bardzokas D.I. Matematicheskoe modelirovanie v zadachah mehaniki svyazannyih po-ley. T.II: Staticheskie i dinamicheskie zadachi elektrouprugosti dlya sostavnyih mnogo-svyaznyih tel. - M.: Komkniga, 2005. 376 s. [Bardzokas D.I. Mathematical

modeling in problems of mechanics of coupled fields. Vol. II: Static and dynamic problems of electroelasticity for compound multi-connected bodies. M .: Komkniga, 2005. 376 p. (In Russ)]

17. Shlyahin D.A. Vyinuzhdennyie osesimmetrichnyie kolebaniya pezokeramicheskoy ton-koy bimorfnoy plastinyi // Izv. RAN. MTT. 2013. №2. S. 77-85 [Shlyakhin D.A. Forced axisymmetric vibrations of a piezoceramic thin bimorph plate // Izv. RAS. MTT. 2013. №2. рр. 77-85. (In Russ)]

18. Shlyahin D.A. Vyinuzhdennyie osesimmetrichnyie kolebaniya tolstoy krugloy zhestko zakreplennoy pezokeramicheskoy plas-tinyi // Izv. RAN. MTT. 2014. №4. S. 90-100 [Shlyakhin D.A. Forced axisymmetric vibrations of a thick circular rigidly fixed piezoceramic plate // Izv. RAS. MTT. 2014. №4. рр. 90-100. (In Russ)]

19. Kovalenko A.D. Vvedenie v termouprugost. -Kiev: Nauk. dumka, 1965. -204 s. [Kovalenko A.D. Introduction to thermoelas-ticity. Kiev: Nauk. dumka, 1965. 204 p. (In Russ)]

20. Grinchenko V.T., Ulitko A.F., Shulga N.A. Mechanics of linked fields in structural elements. Kiev: Nauk. dumka, 1989. 279 р.

21. Senitskiy Yu.E. Metod konechnyih integralnyih preobrazovaniy - obobschenie klassicheskoy protseduryi razlozheniya po sobstvennyim vektor-funktsiyam // Izv. Sara-tovskogo un-ta. Novaya seriya. Matem., mehan., informatika, 2011. № 3(1). S. 61-89. [Senitsky Yu.E. The method of finite integral transformations is a generalization of the classical procedure for expansion in eigenvector functions // Izv. Saratovsky University. New series. Mat., Mechan., Informatics, 2011. No. 3 (1). рр. 61-89. (In Russ)]

22. Kamke E. Spravochnik po obyiknovennyim differentsialnyim uravneniyam. M.: Nauka,1965. 703 s. [Kamke E. Handbook of ordinary differential equations. M .: Nauka, 1965. 703 р. (In Russ)]

Сведения об авторе:

Кальмова Мария Александровна, старший преподаватель кафедры строительная механика, инженерная геология, основания и фундаменты, e-mail: kalmova@inbox.ru Information about the authors:

Kalmova Maria Aleksandrovna, senior lecturer of the department of structural mechanics, engineering geology, foundations and foundations, e-mail: kalmova@inbox.ru

Конфликт интересов: Conflict of interest.

Автор заявляет об отсутствии конфликта интересов. The author declare no conflict of interest.

Поступила в редакцию 12.10.2020. Received 12.10.2020.

Принята в печать 14.11.2020. Accepted for publication 14.11.2020.