Шляхин Д.А., Кальмова М.А. Нестационарная задача термоэлектроупругости для длинного пьезокерамического цилиндра // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2021. - № 2. С. 181-190. DOI: 10.15593/perm.mech/2021.2.16
Shlyakhin D.A., Kalmova М.А. The nonstationary thermoelectric elasticity problem for a long piezoceramic cylinder. PNRPU Mechanics Bulletin, 2021, no. 2, pp. 181-190. DOI: 10.15593/perm.mech/2021.2.16
ВЕСТНИК ПНИПУ. МЕХАНИКА № 2, 2021 PNRPU MECHANICS BULLETIN
https://ered.pstu.ru/index.php/mechanics/index
Б01: 10.15593/регш.шесЬ/2021.2.16 УДК 539.3
НЕСТАЦИОНАРНАЯ ЗАДАЧА ТЕРМОЭЛЕКТРОУПРУГОСТИ ДЛЯ ДЛИННОГО ПЬЕЗОКЕРАМИЧЕСКОГО ЦИЛИНДРА
Д.А. Шляхин, М.А. Кальмова
Самарский государственный технический университет, Самара, Россия
О СТАТЬЕ АННОТАЦИЯ
Построено новое замкнутое решение связанной нестационарной задачи термоэлектроупругости для длинного пьезокерамического радиально поляризованного цилиндра. Рассматривается случай действия на его внутренней цилиндрической поверхности нестационарной нагрузки в виде функции изменения температуры при заданном законе конвекционного теплообмена на внешней лицевой стенке (граничные условия теплопроводности 1-го и 3-го родов). Электродированные поверхности цилиндра подключены к измерительному прибору с большим входным сопротивлением (электрический холостой ход).
Исследуется задача, в которой скорость изменения температурной нагрузки не оказывает влияние на инерционные характеристики упругой системы, что позволяет включить в исходные линейные расчетные соотношения уравнения равновесия, электростатики и теплопроводности относительно радиальной компоненты вектора перемещений, электрического потенциала, а также функции изменения температурного поля. В расчетах используется гиперболическая LS-теория теплопроводности.
Решение задачи осуществляется с помощью обобщенного метода биортогонального конечного интегрального преобразования, основанного на многокомпонентном соотношении собственных вектор-функций двух однородных краевых задач. Структурный алгоритм данного подхода позволяет выделить сопряженный оператор, без которого невозможно осуществить решение несамосопряженных линейных задач математической физики.
Построенные расчетные соотношения дают возможность определить напряженно-деформированное состояние, термоэлектрические поля, индуцируемые в пьезокерамиче-ском элементе при произвольном температурном внешнем воздействии. Подключение электроупругой системы к измерительному прибору позволяет измерить электрическое напряжение.
Анализ численных результатов позволяет, во-первых, установить скорость изменения температурной нагрузки, при которой необходимо использовать гиперболическую теорию теплопроводности и, во-вторых, определить физические характеристики пьезокерамиче-ского материала для случая, когда скорость изменения объема тела приводит к перераспределению температурного поля.
Разработанный алгоритм расчета может быть использован при проектировании нерезонансных пьезоэлектрических датчиков температуры.
©ПНИПУ
Получена: 14 октября 2020 г. Принята: 22 апреля 2021 г. Опубликована: 12 июля 2021 г.
Ключевые слова:
длинный пьезокерамический цилиндр, связанная задача термоэлектроупругости, нестационарное температурное воздействие, конечные интегральные преобразования.
© Шляхин Дмитрий Аверкиевич - д.т.н., доц., зав. каф., e-mail: d-612-mit2009@yandex.ru, : 0000-0003-0926-7388. Кальмова Мария Александровна - асп., e-mail: kalmova@inbox.ru, : 0000-0001 -9212-8631.
Dmitry A. Shlyakhin - Doctor of Technical Sciences, Associate Professor, Head of Department, e-mail: d-612-mit2009@yandex.ru, : 0000-0003-0926-7388.
Maria A. Kalmova - PhD student, e-mail: kalmova@inbox.ru, ¡D: 0000-0001-9212-8631.
Эта статья доступна в соответствии с условиями лицензии Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)
THE NONSTATIONARY THERMOELECTRIC ELASTICITY PROBLEM FOR A LONG PIEZOCERAMIC CYLINDER
D.A. Shlyakhin, M.A. Kalmova
Samara State Technical University, Samara, Russian Federation
ARTICLE INFO
ABSTRACT
Received: 14 October 2020 Accepted: 22 April 2021 Published: 12 July 2021
Keywords:
long piezoceramic cylinder, coupled problem of thermoelectric elasticity, non-stationary temperature effect, finite integral transformations.
A new closed-loop solution for the coupled nonstationary problem of thermoelectric elasticity is designed for a long piezoceramic radially polarized cylinder. The case of the nonstationary load acting on its inner cylindrical surface is considered as a function of temperature change at a given law of the convection heat exchange on the outer face wall (boundary conditions of heat conductivity of the 1st and 3rd types). Electrodynamic cylinder surfaces are connected to a measuring device with a high input resistance (electric idling).
We investigate the problem where the rate of the temperature load changes does not affect the inertial characteristics of the elastic system. It makes it possible to expand the initial linear computational relations with the equilibrium, electrostatics and heat conductivity equations with respect to the radial component of the displacement vector, electric potential as well as the function of temperature field changes. Hyperbolic LS-theory of the thermal conductivity is used in the computations.
The problem is solved with a generalized method of biorthogonal finite integral transformation based on a multicomponent ratio of eigen functions of two homogeneous boundary value problems. The structural algorithm of this approach allows identifying a conjugated operator, without which it is impossible to solve non-self-conjugated linear problems in mathematical physics.
The resulted computational relations make it possible to determine the stress-strain state, temperature and electric fields induced in the piezoceramic element under an arbitrary external temperature effect. By connecting the electroelastic system to the measuring tool, we can find voltage.
Firstly, the analysis of the numerical results allows identifying the rate of the temperature load changes, at which it is necessary to use the hyperbolic theory of thermal conductivity. Secondly, it allows determining the physical characteristics of the piezoceramic material for the case when the rate of changing the body volume leads to a redistribution of the temperature field.
The developed computational algorithm can be used to design non-resonant piezoelectric temperature sensors.
©PNRPU
Введение
В настоящее время пьезокерамические материалы используются при разработке датчиков температуры, работа которых основана на зависимости электрического поля, индуцируемого в пьезокерамическом элементе, от величины внешнего температурного воздействия [1-7]. Для расширения функциональных возможностей измерительных приборов данного типа возникает необходимость углубленного анализа нестационарных процессов, позволяющего понять эффект взаимодействия механических, температурных и электрических полей.
Математическая формулировка начально-краевых задач термоэлектроупругости включает систему несамосопряженных дифференциальных уравнений, исследование которых в последнее время, как правило, проводится при использовании численных методов [8-12]. Однако достаточно слабые эффекты взаимодействия полей различной физической природы удается проанализировать только с помощью замкнутых аналитических решений. При этом проблема интегрирования исходных расчетных соотношений и построение общего решения приводит к проведению расчетов в упрощенной постановке, а именно: исследуются несвязанные задачи [13,14] или анализируются бесконечно длинные тела [15-17].
Замкнутые решения динамических задач термо-электроупругости представлены в немногих исследованиях [5, 14-17]. В работе [5] на основании известных характеристик вынужденных стационарных электроупругих колебаний исследовалась плотность распределения температуры по длине конструкции. Статья [14] посвящена анализу напряженно-деформированного состояния длинного полого цилиндра в рамках Ь8-теории теплопроводности [18] при тепловом ударе без учета влияния электрического потенциала на термоупругие поля. В [15, 16] рассматриваются связанные задачи классической СТЕ-теории [19] для однородного и неоднородного пьезокерамических неограниченных слоев. Использование преобразования Лапласа позволило сформулировать в пространстве изображений интегральное уравнение Фредгольма, которое реализовыва-лось численным методом. В работе [17] предложена модель электротермоупругого полупространства с неоднородным покрытием, где для построения решения использовался численно-аналитический метод.
Можно отметить также работы [20-22], в которых используется нелинейная теория термоэлектроупруго-сти, позволяющая описать отклик пьезокерамических конструкций с учетом неоднородных начальных условий.
Целью настоящей работы является решение связанной нестационарной задачи термоэлектроупругости для длинного полого пьезокерамического цилиндра при действии на его поверхностях температурной нагрузки и учете конвекционного теплообмена с окружающей средой [19] с использованием гиперболической LS-теории теплопроводности [18]. Рассматривается случай, когда скорость изменения нагрузки существенно меньше скорости распространения упругих волн, что позволяет не учитывать инерционные свойства конструкции и использовать в расчетах уравнения равновесия [23, 24].
1. Постановка задачи. Пусть полый длинный незакрепленный пьезокерамический цилиндр занимает в цилиндрической системе координат (г, 6, г) область
О : {а < г < Ь,0 < 6 < 2л, < г < «>}. Рассматривается случай действия на внутренней (г = а) цилиндрической поверхности нестационарной нагрузки в виде функции изменения температуры ю* (г,) (граничное условие 1-го рода), а на внешней (г = Ь) лицевой поверхности задан закон конвекционного теплообмена (граничное условие 3-го рода) и известна температура окружающей среды д*.
Внутренняя электродированная поверхность заземлена, а внешняя подключена к измерительному прибору с большим входным сопротивлением, что соответствует режиму электрического холостого хода.
В общем случае дифференциальные уравнения равновесия, электростатики и теплового баланса на основании гиперболической зависимости Лорда - Шульмана имеют вид [18, 25-27]
ди „ и
дф*
да
гг
дт,
а„ — а0.
= 0, V Д. = 0 .
(1)
дя
д2 ^
д£
дг2
= ЛУ
д& дк
где агг(г*,и),а66(т*,г») - компоненты тензора механических напряжений; Д (г, К ) - радиальная составляющая вектора индукции электрического поля; s (г, 4) — объемная плотность энтропии; ©* = Т — Г0; ©*, Т, Г0 — соответственно приращение, текущая температура, а также температура первоначального состояния тела, при котором отсутствуют механические напряжения; $ге1 — время релаксации; Л — коэффициент теплопроводности материала.
Уравнения осесимметричного состояния электроупругой анизотропного среды при радиальной поляризации пьезокерамического материала с гексагональной кристаллической решеткой класса 6 mm записываются следующим образом [17, 27]:
г ди „ и дф
агг = С33 ^-+ С13~ + б33^--Узз© ,
дг г дг
абб = С13^" + С11 — + ^--Ун©
дг
дг
дф ди и
Дг =—833 ^ + + + &© '
дг дг г а объемная плотность энтропии я (г, 4) при разложении в ряд Тейлора, с учетом условия — •< 1, определяется зависимостью [15]
~ , © 5ф Тп дк
(3)
В равенствах (2),(3) и * (г, и) - радиальная составляющая векторов перемещений; ф* (г, и) — потенциал электрического поля; С^, ^, е33 — модули упругости, пьезомодули и коэффициент диэлектрической проницаемости электроупругого материала (т, я = 1,3); к — коэффициент объемной теплоемкости материала; уп, у33 — компоненты тензора температурных напряжений; g3 — компонента тензора пирокоэффициентов; оператор
дг, г,
В результате подстановки (2), (3) в (1) получаем следующую систему дифференциальных уравнений термоэлектроупругости и краевые условия рассматриваемой задачи в безразмерной форме:
,ди и „ дф 1 дф
©
V--а ——а—1 -у©+а3 — = о, (4)
дг
дг г дг
„дф „ди 1 ди ^ А
-V — + а. V-+ а--+ а, V© = 0,
4 ^ 5 ^ 6 5
дг дг г дг
+ II © + а7VU-а8^1 = 0; дг I дг дг Л 7 8 дг 1
г = Я,1 ^ + а9 и + ^ —© = 0, фк=Д = 0 , ©к=Д =®1, (5)
дг г дг
ди и ---+ а —+ а —+ а © I = 0,
дг дг г ) г=1
д©
--+ а10© | = а,„д;
дг
г = 0 и = ф = © = 0, ди = дф = д© = 0; (6)
дг дг дг
где
{и, г, Я} = {и *, г, а} / Ь,
{©, ю,, = {©*, (ю*- Т,), (д*- Т,)},
г
г
0
Shlyakhin D.A., Kalmova М.А. /РЫЯРЦ Мескапс Би1Ши 2 (2021) 181-190
ф = ф*, г = -Л К, а, = ^ , а2 = ^,
^^ з ь кь с33 е ^
Ун езз
аз =— > а4 =—-
узз ^зз ез:
ез1езз С р зз зз _ Яз^зз , аб ' рзз узз а =
а8 = узз _ е А: ' а9 С С1з С ' Сзз
а10 = а Ь, Р = Рге, , а - коэффициент теплоотдачи.
Л кЬ
Расчетные соотношения (5) учитывают отсутствие механических напряжений на цилиндрических поверхностях, заземление внутренней (г = Я) и подключение к измерительному прибору внешней (г = 1) электроди-рованных поверхностей, а также граничные условия теплопроводности.
С учетом заземления металлической подложки, напряжение холостого хода V (г,) определяется потенциалом электрического поля на внешней поверхности цилиндра:
V (г.) = ф(1, г,).
(7)
2. Построение общего решения. На первом этапе решения выполняется процедура приведения расчетных соотношений (4)-(6) к виду, позволяющему в дальнейшем использовать метод конечных биортогональных интегральных преобразований [28]. Для этого вводятся новые функции и (г, г), х(г, г), Ь (г, г) связанные с
и (г, г), ф( г, г), 0( г, г) следующими соотношениями: и(г,г) = Их (г,г)+и(г,г), ф(г,г) = Н2 (г,г)+х(г,г), (8) ©( г, г ) = Н3 (г, г)+Ь (г, г),
где
{Н , Н2, Нз}={ / ( г ) , /2 ( г ) , / ( г )} ^ (г) +
+ {/4 (г), / (г) , / (г)}3,
Подстановка (8) в (4)-(6) при удовлетворении условий
0 , дН1 Н дН
г = я,1 —1+а—+—2 - Н = о,
дг г дг (9)
Н2|г=Я = 0 Нз\г = Я =Ю1>
дН2 дН, Я1 —- - а, —1 - а —1 - а,Нъ\ = 0,
4 ^ 5 6 з 1 4
дг дг г ,, =
+ а1оНз I = ^о^
дг 1=1
позволяет получить краевую задачу относительно функций и(г,г), х(г,г), Ь(г,г) с однородными граничными условиями:
у^- " + Уд*_а 1д*_уь+а3Ь = Я, (10) дг г дг г дг г
„ди 1 ди
-у—+ а, у— + а —+ауь = я ,
/-ч 4 ^ 5 ^ 6 2 ^
дг дг г дг
удЬ -
дг
— + Р—т II Ь + а Уи - а„— 1= Я ; дг дг2 11 7 8 " 1 з
дг
г = Я,1 ^ + а9 и + дХ-Ь = 0, Х\г=я = 0, Ь\г=я = 0 , (11) дг г дг
дх ди и Л п (дЬ Т \ —- + а — + а5— + а6Ь | = 0, | — + а10Ь I = 0: дг дг г I ^ дг ),=1
г = 0 {и, х, Ь} = -{Н1, Н2, Нз},
д{и, х, Ь} _ д{Н1, Н2, Нз}
д г
д г
(12)
где
„дЯ, Н „дН, 1 дЯ, _ Я, Я = -У—1 + аН -У—2- + а--2 + УН - а—1
дг
дг
дг
„ „дН „дН 1 дН,
Я =у—2 - а У—1 - а —1 - а УН,
дг дг г дг
дН2 дг
Начально-краевую задачу (10)-(12) решаем, используя структурный алгоритм биортогонального конечного интегрального преобразования (КИП) [28]. Для этого вводим на сегменте [Я, 1] КИП с неизвестными компонентами собственных вектор-функций ядер преобразований К ,г)...К3 ,г), N (ц,г)...^3 (ц.,г):
0(Х., г) = | ^ ь (г, г)+а у и (г, г) - а8 Мр^ ^ (я,., гу-сТ, ^13) {и (г, г), х( г, г), ь (г, г )} =
да
= X0(Х1,г){Ж1 (ц.,г),N2 (ц.,г),Жз (ц,г)}\\К,\\2,
1=1
1
||К,.| I2 =| Кз(Я, г) Ыз(Х,, г)гСг,
Я
где Я, ц _ собственные значения соответствующих однородных линейных краевых задач относительно сопряженных К (Я, г) и инвариантных N (ц, г) компонент вектор-функций ядер КИП (к = 1,2, з).
В результате использования алгоритма КИП [22] получаем задачу для трансформанты 0(к1, г):
2
0
а =
г
г
г
р += ян
-г2 -г 11 н
(14)
-О
г = 0 Ог=0 = &0, = -о ; (15)
йг | г=0 —г
а также две системы дифференциальных уравнений, граничные условия относительно неизвестных компонент преобразований К (X,-,г)..К3 (X,,г):
-К К К 1 К ,2 К л V —+а V —а —т2—X,2 а7^ = 0, (16)
-г г
-г г -г
-г
,-к2 к 1 -К ,2
—V-2 + V—1 + а--1 + ^2 а<УК = 0,
7 7 27 г 83"
-г -г г -г
Задача для трансформанты О (14), (15) и сопряженная однородная задача (16), (17) относительно компонент ядра К (X,-,г)...К3 (X,,г) получены в результате
применения вырожденного преобразования (13), а соотношения (18),(19) построены путем применения к полученной (сопряженной) задаче (16),(17) аналогичного (13) КИП с компонентами ядра N (ц,г)...^3 (ц,г).
Общий интеграл уравнения (14) с учетом начальных условий (15) имеет вид
О (X,, г ) = [р( т11 — т2,. )]-1 х
-Оп
х\ р I — ОN2,. I ехр(т1,.г) — РI ^ — О0N1,. I ехр(т2,.г) +
-О
„ -К\ -К К -К „ -> V —3+—1+а,-1—а—2+Х2 К = 0;
-г -г г -г
+|Ян (х){ехрт (г — т)] — ехр[щ, (г — т)]} -т [, (20)
-К К -К-, _ ? „ г = я ,1 —1+а9—+а—2—X: аК = 0,
-г г -г
К2|г=Я = К3г=Я = 0 ,
-К2 -К К 2 . —2+—1+а2 —+х2 аК I = 0,
-г -г г
(17)
К -г
■ + ат К3 I = 0;
N (ц,,г)..Д (ц,г):
V —1 - +v —2 - а —2 ^N3 + = 0, (18)
N
-г г -г г -г
^-М, ^-К 1 -К тг Л -V —2+а V —1+а —1+а VNЪ = 0, -г -г г -г
„Ж 2 Кг -N0 , „
V + ц21 N3 + а7VNl — а8 I = 0 ;
ж N N ,г Л = я,1 —1 + ^ —+—2 - N = 0, = 0
-г
-г
^г = Я = 0,
" 2|г=Я
(19)
д^ дЖ N —2 + а^—1 + а5 1 + аб^ = 0. дг дг г / 1
-N3
-г
■ + аш N3 | = 0 ;
ще Яя = —1(ЯК + ЯК + ЯК)г-г ,
Я
1
О0 =-|(н К + Н2К + Н3 К3),=0г-г ,
Я
1=_1 & К +дН,К +дН3К 1 г-г.
11 дг 1 дг 2 дг 3 Ь
где т, т — корни характеристического уравнения:
Рт2 + т. +Х2 = 0 .
г I I I
В уравнении теплопроводности (10) перераспределение тепла в теле при малых значениях коэффициентов а, а происходит за счет учета скоростей изменения его объема и напряженности электрического поля. Это приводит к тому, что связанность термоупругих полей в большей степени оказывает влияние на трансформанту нагрузки О(Х1, г) и незначительно - на форму функций К (X1,г)..К3 (X1,г) , N (ц,-,г^.Д (ц,г) [23, 24]. Учитывая данный факт, при решении систем (16), (18) принимаем значения а = ц = 0 и в результате получаем следующие выражения функций
К (X,.,г)...К (X,-,г), N1 (ц,гу.Л, (ц,г):
К (X., г)= Д + + В31г-л , (21)
к2 (X-, г) = ^ь-ац 1п (г) + о.цу + Dз,alзт-A4 + ,
К3 (X., г) = А.01400,0 (x,г)+А (x,г)+
+^1600,- А (X,Г ) + Д40 (V ) + ОД (V ) ,
N1 (ц, г) = ЕхгА + К,г-А + £3,^,1 (г) + (г) + £5,017,
^ (ц,.г) = Еиа18гА + Е2,а19г-А +
+ £3^3 (г) + Ел,Уил (г) + £5,0201п(г) + £б,,
N3 (ц, г) = £,Л (цг)+ £Л (цг) ,
где
V,- (г ) = Л [ гА I (г )г - А-г - г -А | Я,т (г )гА-г ]
(т = 1,2),
-V,! (г )
V,3 (г и
-г
- + а5г,1 (г ) + аб-10 (ц,г)
-г ,
0
и
г
Shlyakhin D.A., Kalmova М.А. /РЫЯРЦ Мескапс Би1Ып 2 (2021) 181-190
V» (г И
dV¡2 (г)
+ а5г V, 2 (г) + аб^0 (ц.г)
Сг
ёг,
, я =(1 + а4 )-
(1 - аб + а2 аб - аз )
¿0 (ц,г)
+ (1 - аб)
с20 (ц,г )
ёг
т = 1 (цг)= . (цг),
т = 2 (цг) = Г0 (цг),
А =(а + аа)(1+а)1, ап = -а1а-1, а12 = а2А-1 +1, а13 = а2Л- -1,
^ = -(ааа-1 + а), а15 = а2а6 - а + л (а -1),
а16 = а2аб - «з - Л( аб -1) ,
^ = а [(1 - л2 )(1+а)] , а18 = аЛ-1 + а,
а19 = а^1 - а4, а20 = 1 + а5а17.
В равенствах (21) . (...),^ (...) - обыкновенные функции Бесселя 1-го и 2-го родов нулевого порядка, (...) - неэлементарные функции Ломмеля (Р = -Л,0,Л) [29], Ви...Бы, Е, —Еы - постоянные интегрирования.
Подстановка К (Я,г)..К3 (Я,г), N (ц,г)... N (ц, г) в соответствующие граничные условия
(17),(19) формирует две системы алгебраических уравнений, решение которых позволяет определить постоянные интегрирования Дг..~06г, Е1г...Е6г и собственные
значения Я, ц •
Окончательные выражения функций и (г, г),
ф(г,г), ©(г,г) получим, применяя к трансформанте
(20) формулы обращения (13). В результате, с учетом (8), имеем:
{и(г,г),ф(г,г),©(г,г)} = {Н (г,г),Н2 (г,г),Нз (г,г)}+
да
,(ц.,г),N2 (ц,г),N (ц,г)}|К||2. (22)
I=1
Полученные выражения являются сходящимися в силу полноты систем функций [Кт (Я., г), Мт (ц, г)] |да=1 на интервале г е [ Я,1].
Функции Н (г,гз (г,г) определяются при решении следующих дифференциальных уравнений:
дН Н дН 1 дН
Н
-у—1+а^-у—2+а--2+УНз - а — = 0, (23)
дг
дг
дг
у^ - а ^у^ - а1 ^ - а УН = 0, уН = 0.
г
.дН
„4 у "^'з
дг дг г дг
дг
что позволяет существенно упростить правые части расчетных соотношений (10). В результате подстановки выражений Н\...Н3 в (23) формируются две системы
уравнений относительно / (г).../ (г) и / (г).../ (г). Их решение, при удовлетворении условий (9), позволяет определить функции Н\...Н3.
3. Численный анализ результатов. В качестве примера рассматривается радиально поляризованный пьезокерамический цилиндр (Я = 0.8) состава ПКР-35, который для электроупругих материалов имеет относительно большой коэффициент линейного температурного расширения (а = 1-1* 10-5 К-1) [30]. В расчетах использовались следующие физические характеристики материала [30, 31, 32]: {Сп, С33, С13 } =
{1Э.9, 11.8, 6.2}х1010 Н/м2, е33 = 7.87х10-9 Ф/м,
к = 4 х106 Дж/(м3 К), я = 2 * 10-4 Кл/(м2 К), {е31, е33 } =
= {-7.з,16.2} Кл/м2, {уп,Уз} = {2.9з,2.4Э}х106 Н/(м2 К),
Л = з. 14 Вт/(м К), а = 5.6 Вт/(м2 К), Ры = 5х10-6 с.
Рассматривается случай действия на внутренней поверхности (г = а) температурной нагрузки в виде
(г,) = Т^ [(ег, )Н(Си - К) + Н (г, - г;„)], (24)
где //(?)- единичная функция Хэвисайда (//(/) = I при / > О, Я(?) = 0 при / < 0), Т^=Т^-Т0, Т^,, - максимальное значение внешнего температурного воздействия и соответствующее ему время в
я
размерной форме (= з7з К , Т0 = 29Э К), е = —;— .
2гтах
На рис. 1 представлены графики изменения функций ©*(г, г), ф(г, г) по радиальной координате г в различные моменты времени г (Ь = 0.02 м, = 1 с). Цифрами 1-3 соответственно обозначены результаты,
полученные при значениях г = гшах,10гшах,100гшах .
Анализ представленных графиков позволяет сделать следующие выводы:
1. При достижении температурной нагрузки максимальных значений г = гшах (рис. 1, а, кривая 1) температурное поле изменяется в области, близкой с лицевой нагреваемой поверхности (г = Я). В дальнейшем температурное поле цилиндра растет (рис. 1, а, кривая 2), и полный прогрев пьезокерамической конструкции наблюдается при г = 100гшах (рис. 1, а, кривая 3).
2. Наибольшие значения индуцируемого электрического поля ф(г, г) наблюдаются при достижении температурной нагрузки максимальных значений г = /шах (рис. 1, б, кривая 1). Далее, при постоянном значении
х
г
г
температурной нагрузки происходит уменьшение численных значений ф(г, г) (рис. 1, б, кривые 2, 3).
©*( г, г) ,0 С
U (1, t )х103
/ 3
\ 2 /
\ 1 ^-
.8 0.9 а r, t )х105
1
\
3 \
Рис. 1. Графики изменения функций ©*( r,t), ф(г,г) по радиальной координате r в различные моменты времени: а - ©*(r, t W r ; б - ф(r, tW r (1 - t = t ; 2 - t = 10t ;
V ' > ' t\ ' / V max 5 max 5
* - t = lOO^ma, , 'max
= A • ^ = kb2 tmax'
Fig. 1. Graphs of changes in functions ©*(r,t),r,t) along radial coordinate r at different moments: a - ©*( r, tr ;
6 - <fr( r, ^ r (1 - t = ^ ; 2 - i = 10tm
2
A k?
3 - t = 100i
^max ; j 2 ^max )
На рис. 2 показаны графики изменения перемещений и(1, г) во времени г в случае достижения нагрузки максимальных значений за период время г*ж = 10-4 (с) (гшах = 1,33 х 107). Сплошной и пунктирной линией соответственно обозначены результаты, полученные при использовании гиперболической и параболической (Р = 0) теорий теплопроводности.
Результаты расчета показывают, что при исследовании пьезокерамических радиально поляризованных цилиндров уточненную теорию Лоренца - Шульмана необходимо использовать при очень быстром изменении температурной нагрузки (г*ж < 10-4 с). Это приводит увеличению перемещений до 10 %.
Рис. 2. Графики изменения функции U (1, t) по времени t Fig. 2 . Graphs of function change U (l, t) by time t
При этом необходимо дать оценку влияния сил инерции упругой системы на напряженно -деформированное состояние рассматриваемой конструкции при высокоскоростном температурном воздействии. Численные результаты расчета задачи электроупругости для длинного цилиндра [33] позволяют сделать вывод, что при гармоническом воздействии силы инерции необходимо учитывать при следующем соотношении частоты вынужденных колебаний ю и
первой частоты собственных колебаний у : Ю > 0.5 .
Принимая ю = 8 =
2t;
2
max
У
и учитывая численные резуль-
таты расчета, приведенные в работе [33], получаем, что при < 10-5 данные характеристики оказывают влияние на деформированное состояние электроупругой системы. Таким образом, построенный алгоритм расчета, позволяющий учесть инерцию теплового потока, справедлив при 10-5 < < 10-4 (с).
На рис. 3, 4 представлены графики изменения температуры ©*(г, г) по радиальной координате г и разности потенциалов V(г) во времени г. Сплошной
и пунктирной линиями соответственно обозначены результаты, полученные с учетом и без учета скорости
изменение объема электроупругого тела (— VU).
дг
В настоящем примере при вычислении температурного поля ©* (г,г) и разности потенциалов V(г) наибольший эффект связанности наблюдается соответственно при г = 2гтах иг = гтах . В дальнейшем данный эффект снижается.
Заключение. Полученные численные результаты позволяют сделать вывод, что при исследовании работы электроупругого цилиндра, выполненного из пьезоке-рамического материала с коэффициентом линейного температурного расширения а, > 10—5 К 1, необходимо учитывать связанность температурных и упругих полей. При этом вследствие различных скоростей распростра-
4
3
80
2
60
0
40
1х10
20
r
0
4
2
Л
б
V (t )x109
нения электромагнитных и термоупругих полей можно не принимать во внимание влияние скорости изменения напряженности на температурное поле.
©( r, t) ,0 С
2
1
/
Рис. 3. Графики изменения © *( r ,t) по радиальной координате r в различные моменты времени i1 - t = 2tmax; 2 - t = 10^)
Fig. 3. Graphs of change © *( r, t) along radial coordinate r at different moments (1 - t = 2t , 2 - t = lOt )
V max ' max '
Рис . 4. Графики изменения разности потенциалов V(t) во времени t
Fig. 4. Graphs of difference in potentials V(t) within time t
Построенный алгоритм решения позволяет также учесть передачу тепла в виде теплового излучения. В этом случае нелинейное граничное условие теплопроводности 3 -го рода удовлетворяется путем организации итерационного процесса.
10
4
2
t
1x10
r
Библиографический список
1. Козлов В.Л. Оптоэлектронные датчики. - Минск.: Изд-во Белорус. гос. ун-та, 2005. -116 с.
2. Кульчин Ю.Н. Распределенные волоконно-оптические измерительные системы. - М.: Физматлит, 2001. - 272 с.
3. Дмитриев С.А., Слепов Н.Н. Волоконно-оптическая техника: современное состояние и новые перспективы: сб. ст. - 3-е изд. - М.: Техносфера, 2010. - 608 с.
4. Паньков А.А. Математическое моделирование пьезо-электролюминесцентного эффекта и диагностика распределения давления по длине оптоволоконного датчика// Вестник Пермского национального исследовательского политехнического университета. Механика. - 2016. - № 4. - С. 259-272.
5.Паньков А.А. Резонансная диагностика распределения температуры пьезоэлектролюминесцентным оптоволоконным датчиком по решению интегрального уравнения Фредгольма // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2018. - № 2. - С. 72-82.
6. Ларионов В.А. Резистивный датчик температуры с метрологическим контролем // Датчики и системы. - 2015. -№ 9-10. - С. 76-78.
7. Казарян А.А. Тонкопленочный датчик давления и температуры // Датчики и системы. - 2016. - № 3. - С. 50-56.
8. Abbas I.A., Youssef H.M. Finite element analysis of two-termoperature generalized magneto-thermoelasticity // Arch Appl Mech. - 2009. - № 79. - P. 917-925.
9. T. He [et al.] A generalized electromagneto-thermoelastic problem for an infinitely long solid cylinder // European Journal of Mechanics A-Solids. - 2005. - Vol. 24. - P. 349-359.
10. Youssef H.M. Theory of two-temperature generalized thermoelasticity // IMA J. Appl. Math. - 2006. - Vol. 71(3). -P. 383-390.
11. Kulikov G.M., Mamontov A.A., Plotnikov S.V. Coupled thermoelectroelastic stress analysis of piezoelectric shells // Composite Structures. - 2015. - Vol. 124. - P. 65-76.
12. Куликов Г.М., Плотникова С.В. Решение трехмерных задач термоупругости для слоистых оболочек из функцио-
нальных материалов // Вестник ТГТУ. - 2015. - Т. 21, № 1. -С. 185-190.
13. Фирсанов В.В., Нгуен, Ле Хунг. Напряженно-деформированное состояние произвольных оболочек с учетом термоэлектрического воздействия на основе уточненной теории // Тепловые процессы в технике. - 2010. - № 3. - С. 110-117.
14. Abbas I.A., Zenkour A.M. LS model on electro-magneto-thermoelastic response of an infinite functionally graded cylinder // Composite Structures. - 2013. - Vol. 96. - P. 89-96.
15. Ватульян А.О., Кирютенко А.Ю., Наседкин А.В. Плоские волны и фундаментальные решения в линейной тер-моэлектроупругости // ПМТФ. - 1996. - Т. 37, № 5. - С. 135142.
16. Ватульян А.О., Нестеров С.А. Динамическая задача термоэлектроупругости для функционально-градиентного слоя // Вычислительная механика сплошных сред. - 2017. -Т. 10, № 2. - С. 117-126.
17. Белянкова Т.И., Калинчук В.В. К моделированию преднапряженного термоэлектроупругого полупространства с покрытием // Изв. РАН. МТТ. - 2017. - № 1. - С. 117-135.
18. Lord H., Shulman Y. A generalized dynamical theory of thermoelasticity // Elasticity. -1967. - P. 299-309.
19. Коваленко А.Д. Введение в термоупругость. - Киев: Наук. думка, 1 965. - 204 с.
20. Montanaro A. Some theorems of incremental thermo-electroelasticity // Arch. Mech. - 2010. - Vol. 62. - P. 49-72.
21. Ueda S. Thermally induced fracture of a functionally graded piezoelectric layer/ Journal of Thermal Stresses. - 2004. -Vol. 27(4). - P. 291-309.
22. Yang J.S., Equations for Small Fields Superposed on Finite Biasing Fields in a Thermoelectroelastic Body // IEEE Transactions on Ultrasonics, Ferroelectricts, and Frequency Control. -2003. - Vol. 50/2. - P. 187-192.
23. Шляхин Д.А., Даулетмуратова Ж.М. Нестационарная осесимметричная задача термоупругости для жесткозакреплен-
ной круглой пластины // Инженерный журнал: наука и инновации. -2018. - Вып. 5(77). DOI.1018698/2308-6033-2018-5
24. Шляхин Д.А., Даулетмуратова Ж.М. Нестационарная связанная осесимметричная задача термоупругости для жесткозакреп-ленной круглой пластины // Вестник Пермского национального исследовательского политехнического университета. Механика. -2019. - № 4. - С. 191-200. DOI: 10.15593/peim.mech/2019.4.18
25. Лычев С.А., Манжиров А.В., Юбер С.В. Замкнутые решения краевых задач связанной термоупругости // Изв. РАН. МТТ. - 2010. - № 4. - С. 138-154.
26. Радаев Ю.Н., Таранова М.В. Волновые числа термоупругих волн в волноводе с теплообменом на боковой стенке // Вестник Самарского государственного технического университета. Серия: Физико-математические науки. - 2011. -№ 2(23). - С. 53-61.
27. Гринченко В.Т., Улитко А.Ф., Шульга Н.А. Механика связанных полей в элементах конструкций. - Киев: Наук. думка, 1989. - 279 с.
28. Сеницкий Ю.Э. Биортогональное многокомпонентное конечное интегральное преобразование и его приложение к
References
1. Kozlov V.L. Optoelektronnyie datchiki [Optoelectronic sensors]. Minsk, Publishing house of Belarus. state un-that, 2005, 116 p.
2. Kulchin Yu.N. Raspredelennyie volokonno-opticheskie izmeritelnyie sistemyi [Distributed fiber-optic measuring systems]. Moscov, Fizmatlit, 2001, 272 p.
3. Dmitriev S.A., Slepov N.N. Volokonno-opticheskaya tehnika. Sovremennoe sostoyanie i novyie perspektivyi [Fiber optic technology. Current state and new perspectives]. Moscow, Technosphere. Sat. articles, 3rd ed., 2010, 608 p.
4. Pankov A.A. Mathematical modeling of the piezoelectro-luminescent effect and diagnostics of pressure distribution along the length of the fiber optic sensorPNRPU Mechanics Bulletin, 2016, no. 4, pp. 259-272.
5. Pankov A.A. Resonant diagnostics of temperature distribution with a piezoelectric-luminescent fiber-optic sensor by solving the Fredholm integral equation. PNRPUMechanics Bulletin, 2018, no. 2, pp. 72-82.
6. Larionov V.A. Rezistivnyiy datchik temperaturyi s metrologi-cheskim kontrolem [Resistive temperature sensor with metrological control]. Sensors and systems, 2015, no. 9-10, pp. 76-78.
7. Kazaryan A.A. Tonkoplenochnyiy datchik davleniya i tem-peraturyi [Thin-film pressure and temperature sensor]. Sensors and systems, 2016, no. 3, pp. 50-56.
8. Abbas I.A., Youssef H.M. Finite element analysis of two-temperature generalized mag-neto-thermoelasticity. Arch Appl Mech, 2009, vol. 79, pp. 917-925.
9. He T. et al. A generalized electromagneto - thermoelastic problem for an infinitely long solid cylinder // European Journal of Mechanics A - Solids, 2005, vol. 24, pp. 349-359.
10.Youssef H.M. Theory of two - temperature generalized ther-moelasticity. IMA J. Appl. Math, 2006, vol. 71 (3), pp. 383-390.
11. Kulikov G.M., Mamontov A.A., Plotnikov S.V. Coupled thermoelectroelastic stress anal-ysis of piezoelectric shells. Composite Structures, 2015, vol. 124, pp. 65-76.
12. Kulikov GM, Plotnikova SV Reshenie trehmernyih zadach termouprugosti dlya sloistyih obolochek iz funktsionalnyih materialov [Solution of three-dimensional problems of thermoelas-ticity for layered shells from functional materials]. Vestnik TSTU, 2015, vol. 21, no. 1, pp. 185-190.
краевым задачам механики // Известия вузов. Математика. -1996. - № 8. - С. 71-81.
29. Камке Э. Справочник по обыкновенным дифференциальным уравнениям. - М.: Наука,1965. - 703 с.
30. ГОСТ Р 8.945-2018. Теплофизические характеристики пьезокерамик на основе ниобата лития в диапазоне температур от 300 К до 900 К. - М.: Стандартинформ, 2018. (www.docs.cntd.ru).
31. Панич А.А. Мараховский М.А., Мотини Д.В. Кристаллические и керамические пьезоэлектрики // Инженерный вестник Дона. - 2011. - № 1 (www.indon.ru).
32. Бабенков М.Б. Анализ распространения гармонических возмущений в термоупругой среде с релаксацией теплового потока // ПМиТФ. - 2013. - № 2(54). - С. 126-137.
33. Шляхин Д.А. Динамическая осесимметричная задача прямого пьезоэффекта для анизотропного пьезокерамического радиально поляризованного цилиндра // ПМиТФ. -2010. -№ 1(51). - С. 153-161.
13. Firsanov V.V., Nguyen, Le Hung Napryazhenno-deformirovannoe sostoyanie proizvolnyih obolochek s uchetom termoelektricheskogo vozdeystviya na osnove utochnennoy teorii [Stress-strain state of arbitrary shells taking into account thermoelectric effect on the basis of a refined theory]. Thermal processes in technology, 2010, no. 3, pp. 110-117.
14. Abbas I.A., Zenkour A.M. LS model on electro-magneto-thermoelastic response of an in-finite functionally graded cylinder. Composite Structures, 2013, vol. 96, pp. 89-96.
15. Vatulyan A.O., Kiryutenko A.Yu., Nasedkin A.V. Ploskie volnyi i fundamentalnyie resheniya v lineynoy termoelektro-uprugosti [Plane waves and fundamental solutions in linear ther-moelectroelasticity]. PMTF, 1996, vol. 37, no. 5, pp. 135-142.
16. Vatulyan A.O., Nesterov S.A. Dinamicheskaya zadacha termoelektrouprugosti dlya funktsionalno-gradientnogo sloya [Dynamic problem of thermoelectroelasticity for a functional-gradient layer]. Computational mechanics of continuous media, 2017, vol. 10, no. 2, pp. 117-126.
17. Belyankova T.I., Kalinchuk V.V. K modelirovaniyu pred-napryazhennogo termoelektrouprugogo poluprostranstva s pokryit-iem [On modeling a prestressed thermoelectric-elastic half-space with a coating]. Izv. RAS. MIT, 2017, no. 1, pp. 117-135.
18. Lord H., Shulman Y. A generalized dynamical theory of thermoelasticity. Elasticity, 1967, pp. 299-309.
19. Kovalenko A.D. Vvedenie v termouprugost [Introduction to thermoelasticity]. Kiev, Nauk. Dumka, 1965, 204 p.
20. Montanaro A. Some theorems of incremental thermoelec-troelasticity. Arch. Mech., 2010, vol. 62, pp. 49-72.
21. Ueda S. Thermally induced fracture of a functionally graded piezoelectric layer. Journal of Thermal Stresses, 2004, vol. 27 (4), pp. 291-309.
22. Yang J.S., Equations for Small Fields Superposed on Finite Biasing Fields in a Thermoe-lectroelastic Body. IEEE Transactions on Ultrasonics, Ferroelectricts and Frequency Control, 2003, vol. 50/2, pp. 187-192.
23. Shlyakhin D.A., Dauletmuratova Zh.M. Nestatsionarnaya osesimmetrichnaya zadacha termouprugosti dlya zhestkozakrep-lennoy krugloy plastinyi [Non-stationary axisymmetric thermoelasticity problem for a rigidly fixed round plate]. 2018, Engineer-
ing journal: science and innovations, Issue. 5 (77), DOI:1018698/2308-6033-2018-5.
24. Shlyakhin D.A., Dauletmuratova Zh.M. Unsteady coupled axisymmetric thermoelasticity problem for a rigidly fixed round plate. 2019, PNRPU Mechanics Bulletin, no. 4, pp.191-200, DOI: 10.15593 / perm.mech / 2019.4.18.
25. Lychev S.A., Manzhirov A.V., Hubert S.V. Zamknutyie resheniya kraevyih zadach svyazannoy termouprugosti [Closed solutions of boundary value problems of coupled thermoelasticity]. Izv. RAS. MTT, 2010, no. 4, pp. 138-154.
26. Radaev Yu.N., Taranova M.V. Volnovyie chisla termou-prugih voln v volnovode s teploobmenom na bokovoy stenke [Wavenumbers of thermoelastic waves in a waveguide with heat transfer on the side wall]. Vestn. Himself. state tech. un-that. Ser. Phys. - mat. Science, 2011, no. 2 (23), pp. 53-61.
27. Grinchenko V.T., Ulitko A.F., Shulga N.A. Mehanika svyazannyih poley v elementah konstruktsiy [The mechanics of related fields in structural members]. Kiev, Nauk. Dumka, 1989, 279 p.
28. Senitsky Yu.E. Biortogonalnoe mnogokomponentnoe konechnoe integralnoe preobrazovanie i ego prilozhenie k kraevy-im zadacham mehaniki [Biorthogonal multicomponent finite inte-
gral transformation and its application to boundary value problems in mechanics]. Izv. universities. Maths, 1996, no. 8, pp. 71-81.
29. Kamke E. Spravochnik po obyiknovennyim differentsial-nyim uravneniyam [Handbook of ordinary differential equations]. Moscow, Science, 1965, 703 p.
30. GOST R 8.945-2018. Teplofizicheskie harakteristiki pe-zokeramik na osnove niobata litiya v diapazone temperatur ot 300 K do 900 K [Thermophysical characteristics of piezoceramics based on lithium niobate in the temperature range from 300 K to 900 K]. Moscow, Standartinform, 2018. (www.docs.cntd.ru).
31. Panich A.A. Marakhovsky M.A., Motini D.V. Kristalli-cheskie i keramicheskie pezoelektriki [Crystalline and ceramic piezoelectrics]. Engineering Bulletin of the Don, 2011, no. 1 (www.indon.ru).
32. Babenkov M. B. Propagation of Harmonic Perturbations in a Thermoelastic Medium with Heat Relaxation. Journal of Applied Mechanics and Technical Physics, 2013, no. 2, pp. 277-286.
33. Shlyakhin D.A. Dynamic axisymmetric direct piezoelectric effect problem for an anisotropic piezoceramic radially polarized cylinder. Journal of Applied Mechanics and Technical Physics, 2010, no. 1, pp. 253-161.