Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Электромагнитные эффекты в твердых телах

..pdf
Скачиваний:
7
Добавлен:
20.11.2023
Размер:
11.93 Mб
Скачать

1.17. Пространственные задачи термоупругости

71

через д2/д г 2, получим уравнения для функций ср, ф в цилиндри­ ческой системе координат. Перемещение u = (иг, uv, w) и тем­ пература 0 определяются в цилиндрической системе коорди­ нат следующими формулами:

 

u r =

Рс4,, W

+

ад1) Ф “

С4 4 *%Ф>

мф =

О,

 

(13)

w =

р с44~

(bS7; +

ид£) ф+

с44 (/V? + dj) ф,

 

 

где

G=

c33c44WvH^)(M-3Vr +

^ ) ^

r = (*? + ^ )1/2’ *з = г>

 

 

 

 

 

 

 

 

 

a =

Ti — х(1 + yt\), b =

 

— \ ц — 1,

х = р7р,

5 = 1 + у л »

 

f =

С\\

 

PI +

Из + 2у

Сзз

V

с13

 

с44

2

2

2

'

 

 

f1! Из

У

с44

 

С33

Функции ф и ф можно использовать для определения темпе­ ратурных напряжений в простых системах: в неограниченном пространстве, упругом полупространстве и упругом слое. В других случаях, например для определения напряженного состояния в толстой круговой или прямоугольной плите, сле­ дует ввести функции %i, %2, что позволит удовлетворить всем краевым условиям.

С помощью функций ф и ф решено несколько практически важных случаев. Найдено решение в замкнутом виде для случая действия источника тепла в упругом пространстве или полупространстве. Аналогичное решение получено для случая действия термоупругого ядра деформации. Рассмотрена так­ же задача установившегося нагрева поверхности, ограничи­ вающей упругое полупространство. При этом показано, что напряжения в направлении,, перпендикулярном плоскости, ограничивающей упругое полупространство, не исчезают. Эти напряжения исчезают только в предельном случае, при пе­ реходе от поперечной изотропии к полной.

До сих пор решено только небольшое число, главным об­ разом одномерных, задач термоупругой динамики. Заслужи­ вает внимания работа Чедвика и Сита [9], относящаяся к тем­ пературным напряжениям в телах с поперечной изотропией. Остановимся на уравнениях в перемещениях для этой дина­ мической трехмерной задачи. Отправным пунктом анализа те­

перь послужат уравнения, которые, согласно

(7),

имеют вид

V2(сц

г12)

+ V2(сц + с12) Мр ар+ с44иа 33-|- (с13Н-с41) иэ а3

 

 

 

 

рЙа=

Р0, а»

(14)

c44V -ия -}- С33М3 ?з *4* (Гр И- г44) Мр 3р

Р#зз=

Р з*

 

Се0 +

^ о ( Р « а . а + Р Ч . з ) = ^ ; 0 + Х ' 0 з з + Г ,

«,

Р = 1 ,

2.

72

Гл. 1. Основы линейной теории пьезоэлектричества

Систему этих уравнений удается частично разделить пу­ тем введения трех скалярных функций ф, %, ф, связанных с пе­ ремещениями следующим образом:

(15)«i = f i + 3C,2. »2 = Ф,2 Хл. «з = ,Ф.з*

Подставляя (15) в уравнения (14), получим

 

 

+ С44^\ 33 + ( С13 + С44) Н>. 33 - Р ^ =

Р 0 ’

(16)

(Cjg +

С44) У^ф

С44^Ф “Ь С33Ф зз

РФ =

Р 0»

 

*е0 +

(p v# + РЧ зз)= *v?0 +

м

зз +

w,

(17)

 

v2(C„ -

Cl2) Vfy + c44x, 33-

PX =

0-

 

Заметим, что только потенциалы ф и <р связаны с температур­ ным полем. Волна SH описывается потенциалом %. На эту волну поле температур не действует. Очевидно, что среде, ха­ рактеризуемой поперечной изотропией, можно сообщить три типа волн в каждом направлении: квазипродольную, квазипоперечную и чисто поперечную. Квазипродольная и квазипоперечная волны модифицированы посредством связи с полем температуры; эти волны затухают и подвергаются дисперсии.

Проведенный обзор показывает, что осталось еще обшир­ ное поле научной деятельности. Мало развиты исследования несвязанной динамической термоупругости. Ожидают разра­ ботки связанные задачи, как квазистатические, так и динами­ ческие. Приведенный здесь результат является обобщением задачи Стернберга и Мак-Доуэлла [47].

Для случая произвольной анизотропии тела до сих пор не найдено ни одного решения, даже в случае неограниченного пространства. Имеется только формальное решение, получен­ ное путем применения четырехкратного преобразования Фурье [6].

1.18. ПЛОСКИЕ ЗАДАЧИ АНИЗОТРОПНОЙ ТЕРМОУПРУГОСТИ

Плоские задачи теории температурных напряжений для анизотропных тел исследованы в научной литературе значи­ тельно шире, нежели пространственные. Так, Моссаковский [37] рассмотрел действие источника тепла в полубесконечной пластине, проявляющей изогональную анизотропию. В своей работе он использовал метод функций комплексной перемен­ ной и ввел обобщенную функцию Эйри и Маргерра.

В случае ортотропной пластины обозначим через Е\ и Е%

модули

Юнга в

направлении осей х\

и х2> через G

модуль

сдвига

и через

vi, v2 коэффициенты

Пуассона [49].

Пусть

о&1, а2— коэффициенты теплопроводности в направлениях осей Х\ и Хч. Уравнение теплопроводности для ортотропной пла­

 

1.18. Плоские задачи анизотропной термоупругости

73

стины примет вид

 

(1)

М + 1 2д 1 - с гд()В=-----W.

 

Определяющие соотношения, соотношения между напряже­ ниями, деформациями и приращением температуры, имеют следующий вид:

8 ц = ЯцСТц Ч~ #12^22 Ч " а 10» ®22 = = ^21 4 “ ^22*^22 4 “ а 20> ®12 ^ ^66^12»

(2)

______ 1

______1

Y I

___

Y 2

______I

aw

, а22

, ai2 ~ Ei

°г1

Е2

а&е> 2Q

 

 

Elv2 = E2v l.

 

 

Подставляя эти соотношения в уравнения неразрывности

(3)

^ е22Ч-^2е11= 2а1д2е12

 

и выражая напряжения через функцию Эйри

(4)

°ар = - W

+

«, Р =

1, 2,

приходим к дифференциальному уравнению

 

(В)

%*d\F + 2w*d\d\F + d\F =

- Е х (ахд\ + а2д|) 0,

где

* * = ! ; .

2пи2 = В

, ( 4 - - ^

- ) .

Решение уравнения (5) составим из двух частей — из частного

интеграла F неоднородного уравнения (5)

и общего интеграла

F однородного квазибигармонического уравнения:

(6)

%*d\F + 2r\y?d]df -f d\F =

0.

Представленный здесь способ действия, использующий функ­ цию напряжений Эйри, удобен тогда, когда краевые условия заданы в напряжениях.

На рассматриваемую задачу об ортотропной пластине лег­ ко перенести «плиточную аналогию», рассмотренную Дюба и Треммелем для изотропных пластин [16, 56]. Если решить соотношения (2) относительно напряжений и внести эти по­ следние в уравнения равновесия, то вместе с равенством (1) получим систему трех уравнений

(7)

 

 

L u i i ^ - W b , з,

«3 = 0,

/ = 1 , 2 , 3 ,

 

где

■^п =

c\fi\ Ч- свбЩ*

 

L22=

C(jg<9^ 4~ ^22^2*

 

 

 

 

 

Ll2=

L2l =

,2Ч" Cjg) 5j52,

Ln =

\ d ] +

K d \ - c tdt,

.CJ

 

со II

1 TD Cb

 

 

Lat=

^-32 “

^-23~

 

 

 

 

 

 

 

 

 

 

nt ca

 

 

 

 

 

 

r ^

V

 

 

ъЕ\

С

1 — ViV2

C,J,J

1 — v,v2

,£2

C2I — , — V,Vj '

<,1 2 = | _ V,V!’

пп

о

i -

£ i(a i+ a 2V2)

g ,(a ,-a ,v L)

E

= E2Vt.

св0 =» 2G,

P

! -

VjVl .

Г2

1 — VlV,

* i

74 Гл. 1. Основы линейной теории пьезоэлектричества

Выразим

функции иа

( а = 1 ,2 ),

и3 = 0 через

три

функции

X,-

(i = 1,2, 3)

следующим образом:

 

 

 

 

 

 

Xi

L.\2

A 3

Ln

Xi

Аз

^11

^12

Xl

(9)

w ,=

X2

•^22

^23

1 ^2 L2\

X2

^23 . «3 =

^21

^-22

X2

 

 

Хз

^32

о

 

Хз

•^33

^-31

^32

Хз

 

 

CO CO

 

Функции Xi должны удовлетворять уравнению

A 'l

C- to

7 -

CO

(10)

■^21

L>22 ^ 2 3

%i = - b

i3W, i = 1,

2, 3,

 

или

^-31

L$O ^ 3 3

 

 

 

 

(л,а? +

х д -

 

 

а*) Х( .— т №

(II)

С(а,) (ц?а?+ ai) (ц’з*+

где

 

I U ± V°* — b

 

a >

1,

 

 

 

(12)

и ? , - » 2!

ff'

 

 

a = l ,

r w

' r y .

 

 

 

u v

a <

1.

 

 

 

 

 

 

Функции хь %2 удовлетворяют однородному, а функция хз —

неоднородному дифференциальному уравнению. Способ дей­ ствия здесь следующий. Из уравнения (11) для i = 3 опреде­

ляем частный интеграл хз» а с помощью функций хь Хз УД°В' летворяем заданным краевым условиям.

1.19СВЯЗАННОСТЬ МЕХАНИЧЕСКИХ

ИЭЛЕКТРОМАГНИТНЫХ ВОЛН

Впредыдущих рассуждениях мы имели дело со связан­ ностью квазистатического электрического поля с движением упругого тела. В такой теории уравнения движения теории упругости были связаны посредством пьезоэлектрических по­ стоянных с уравнением Гаусса divD = 0. Теперь рассмотрим более общую проблему, а именно как упругие, так и электро­ магнитные динамические задачи. Отбросим ранее сделанное предположение

(1)

rotE = 0, Е = — gradqp,

вводившее

природу квазистатического электрического поля,

и обратимся к полной системе уравнений Максвелла (однако

в предположении, что

= О, J = О, М = 0)

(2)

rot И =

D,

rot Е =

— В,

(3)

div 0 =

0,

divB =

0,

1.19. Связанность механических и электромагнитных волн

75

дополненную соотношениями

 

 

(4)

D = e 0E +

P, В = ц0Н.

 

Производя

над равенством

(2)2 операцию ротора, используя

уравнение

(2) t и соотношение (4)2, приходим к уравнению

(5)

rot rot Е = — p,0d2D/d/2.

 

Представим теперь вектор D с помощью определяющего со­ отношения для квазистатической задачи (формула (5) из § 1.4)

(6)

Di = бiki%i + 3ife£fc.

Подставляя (6) в (5), получаем систему трех уравнений, в ко­ торых в качестве неизвестных функций выступают перемеще­ ния щ и составляющие электрического поля Остальные три уравнения найдем из уравнений движения

(7)

Оц, / + Х{ = рй/,

в которые подставим определяющие соотношения

(8)

оif = ci}klult k ekijEk.

Рассмотрим простой пример, относящийся к распростране­ нию монохроматической плоской волны, упругой^ и электро­

магнитной, в кристалле, принадлежащем классу 42т тетраго­ нальной системы и характеризуемом десятью независимыми упругими, пьезоэлектрическими и диэлектрическими постоян­ ными. Полагаем, что волна распространяется в направлении оси х\. Определяющие соотношения в рассматриваемом слу­ чае имеют вид

(9)

СГи —

C\\U\t 1, <Ti3— £44^3, 1+ Si4E2, <J12— ^66^2. i

e36E3i

(10)

D i=3uE \,

D2 = — ^14^3, i "Ь Эц£2> Z)3 = — е3бМо. i +

э33£3.

Подставляя (9) и (10) в уравнения (5) и (7), получим си­ стему пяти дифференциальных уравнений

(11)

 

 

(сцд| ~ Р^<) и\ = 0,

 

 

( 12)

(

(С66^Г

Р^<) ^2 *36^3 ~

 

 

{

(д\ -

№xidt) Е3+

V«eXdPlU2=

 

 

 

(13)

|

{СНд\ ~ Р<5?) «3 +

eUdlE23 °*

 

1

Е2“Ь P’0ei4<^ l M3 ^

 

 

 

Шестое уравнение здесь не представлено,

поскольку Еt =

0,

D\ =

0. Заметим, что продольная волна ii\

не возмущена элек­

тромагнитным

полем. Волны «2,

£з и п3,

Е2 связаны. Всего

имеем

пять фазовых

скоростей.

Первая из них относится

к

76

Гл. 1. Основы линейной теории пьезоэлектричества

 

продольной

волне (11):

о = (с п /р )1/2.

Остальные

фазовые

скорости получим из системы уравнений

(12) и (13),

которые

можно преобразовать к виду

 

 

 

(14)

[(0? -

ц0э33<5?)

- р0?) - |1о4$9?] («2, В,) =

о,

(15)

[(а ? -й 0эпа*)(с<45 ; - р а г ) - м

где?](Нз, я а) = о .

В случае монохроматической волны,

распространяющейся в

направлении оси хи принимаем

 

 

 

(16)

Uj =

u.je-Uat- kXi\

*= Я/в- 4

 

/ = 2, 3.

Вводя (16) в уравнение (14), получим следующие характе­ ристические уравнения для величин %= k/со = \/v:

(17)

 

 

— £2 (*рг + y s + л) +

~^ут— 0.

 

Здесь введены следующие обозначения:

 

2

 

£«=-1

 

с — (

С°6 У /2

у

1

 

 

 

" •

 

©

\

р

/

(ЦоЭзз)1/2 ’

Сев

Решение биквадратного уравнения (17) дает

 

(18)

6?,2 — '|~ {J r + “i72 +

Tl± [(рг +

ут + .)2-^ W n

ИЛИ

 

 

 

 

 

 

 

 

 

(19)

k\ 2 =

 

 

+

+

+

т? +

л )2+ - ^ - ] 1/2}«

Из формул (18) и (19) видно, что

k\ >

0, /г |> 0 . Решение

уравнения

(14), следовательно, имеет вид

 

 

 

и2 =

е~ш {AeiklXt +

Be~ikiXi + CeikiXi - f De~iktXl}t

где

E3 — e~iat { AeiklXl Be~lklXi)%\ +

(CeiklX>De~lklXl)%2}>

 

 

 

 

 

 

 

 

 

 

 

.

V2 - k \

 

 

. v 2 - k \

 

(0

*36

 

K = l - T k T '

 

 

V = = T ’ T = = ^ 7 -

Поперечная волна w2 и электромагнитная волна Е3 распро­ страняются с одной и той же фазовой скоростью v = £-1 = = со//г. Поскольку существует два корня k\t Л2, то мы имеем дело с двумя волнами. Так как k\, k% пропорциональны со, то волна не подвержена дисперсии, а так как они вещественны, волна не затухает.

Такое же решение можно провести для волн ц3, Я2, описы­ ваемых уравнением (15). Отличие здесь состоит лишь в дру­ гих коэффициентах (с44, ец, ей). Знание перемещения и и элек­ трического поля Е позволяет определить напряжения по фор­ муле (9) и вектор D по формуле (10). Составляющие поля Н

1.19. Связанность механических и электромагнитных волн

77

определим из уравнения

(2) i, а вектор В— по формуле

(4) 2.

Предположим далее,

что проводимость тока отлична от

нуля. В этом случае имеем дело с довольно сложной систе­

(21)

rot Н =Ь+J,

rot Е =—В,

мой уравнений электромагнитного поля

(22)

div D = 0,

div В= 0.

Примем следующий вид определяющихP, В соотношений:

(23) D = e 0E + = ц0Н, J = crE.

JЗдесь предполагается пропорциональность между векторами

и Е и изотропность материала относительно электропрово­ димости. Из уравнений (21) j, 2 с учетом (23) получаем сле­ дующие волновые уравнения:

(24) rot rot Е = — Р0-^У1г— [•‘•offЩ-

Дальнейший ход действий аналогичен тому, который был

установлен ранее для случая J =

0._Для волн и2, £з и

£ 2

в случае кристалла с симметрией

42т получим следующие,

более общие уравнения:

 

 

 

 

(25)

[(c66d2 -

pdf) (д] -

|х0эзза| -

ц0<тд,) -

ц0в326Э Д ] (и2, Б3) *

0,

(2в)

[(с44а2-

р<52) (а2 -

р0эид2 -

ц0ога,) -

ц0е24а2а2] (н3, е2) =

о.

Подставим (16) в уравнения (25). Для рассматриваемой здесь плоской волны, распространяющейся в направлении оси х\, получим следующее характеристическое уравнение:

с, - г 2( ж + т ^ + ч + - ^ ) + ^ + г^ - = ° -

(27)

со

 

Т ’

откуда

л2, 2 = - y { - ^ + - p J + 'rl + i'X=fc

(28)

V ( * 2 + V 2

11

с7 ( V 2

**) } ’

Х = аро/ю.

Ясно, что корни k$ (р = 1,2) комплексные. Следовательно,

(29)

= ар +

Р = 1, 2.

Волны затухают. Фазовая скорооть и коэффициент затухания найдем по формулам

(30)

со _

со

р = Im ftp, р — 1, 2.

 

Re*ft ’

 

 

 

РР

78

Гл. 1. Основы линейной теории пьезоэлектричества

Корни k$ зависят нелинейно от частоты колебаний. Следова­ тельно, волны м2, Ез подвержены дисперсии-. Решение урав­ нений (25) имеет вид

Соотношения между постоянными А,

В,

... и А, В, ... найдем

из следующей системы уравнений:

 

 

 

(сб6^1

Р^1) И2 “Ь *36^3

О»

(33)

 

 

д () Е 3 =

°*

Р0е 36д *д 1и 2 + ( д 1 ~

»0ЭЗЗд 1 -

являющейся обобщением уравнений (12).

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

(34)

rot rot Е = — Ро-^2----- ~ д Г '

(35)

/ — рй* = О,

(36)М - и — П {Уцёц + giEt) = 0.

Последнее уравнение является связанным уравнением тепло­ проводности. В уравнения (34) и (35) подставим определяю­ щие соотношения теории термопьезоэлектричества

(37)

aij — Cijkfikl ekijEk — Yi/®>

(38)

Di — eikieki + BikEk + £ $ •

Подставляя (37)

и (38) в дифференциальные уравнения (34)

и (35), получим вместе с уравнением (36) систему семи урав­ нений, в которых неизвестными функциями являются переме­ щение и, электрическое поле Е и температура 0. Все волны будут затухающими и подвергаться дисперсии. Это следует, с одной стороны, из необходимости учитывать член dE/dt, а с другой — из связанности поля температуры с остальными полями.

Глава 2

ТЕОРИЯ ПЬЕЗОЭЛЕКТРИЧЕСТВА ТУПИНА И МИНДЛИНА

2.1. ВВЕДЕНИЕ

Тела, которые не являются проводниками, мы будем назы­ вать диэлектриками. Эти тела электрически нейтральны и со­ держат одинаковое число положительных и отрицательных

зарядов. Введение диэлектрика в электрическое поле изме­

няет это поле. Векторы э0Е и D уже не будут параллельны,

их разность дает вектор поляризации

(1)

Р= D — э0Е.

Нелинейная теория диэлектриков была разработана Тупином [54, 55]. В эту общую теорию, учитывающую эффекты фото­ упругости, эффект Фарадея и другие, входит как частный слу­ чай теория Фойгта классического пьезоэлектричества.

Представим в сильно сокращенном виде ход действий Тупина, приводящий к соотношениям и дифференциальным уравнениям пьезоэлектричества. Разложим внутреннюю энер­ гию U диэлектрика на энергию, связанную с деформацией и поляризацией, UL{ti},Pi), и энергию, связанную с электриче­ ским полем. Отправным пунктом рассуждений будет уравне­ ние баланса энергии в виде (см. формулу (22) § 1.3)

(2)

U = aijBi} + EiDt.

Это уравнение показывает, что

(4)

и = UL (вф Pi) + 7*0а д , Ef = - ф. I

или

 

(4')

U = UH*t}, Pt)~\~ 7гэоФ. /Ф. I-

Вводя электрическую энтальпию

(5)

Н — U — EiD{

и учитывая выражение (40, получаем

(6)

Н — UL {el}, Pt) — У2Э0Ф, *<р. i + Ф. iP{.

Рассмотрим тело В, занимающее объем у, ограниченный по­ верхностью дВ. Эта поверхность отделяет тело В от ва­ куума В'. В дальнейших рассуждениях используем принцип

80 Гл. 2. Теория пьезоэлектричества Тупина и Миндлина

Гамильтона, примененный Тупином к диэлектрикам:

/2 ti

(7) б \ d t

'\{ K ^ - H )d v +

$ л П (* .6 м ,. +

£°6Р.)сго +

и

в•

11 LB-

 

 

 

 

+

^

Pt btii da~\=0,

 

 

 

d B

J

Здесь К — V2Pvtvi — кинетическая энергия, отнесенная к еди­ нице объема, Е°. — внешнее электрическое поле и В* =*

= В + В'. Определим виртуальное приращение бЯ:

(8) 6Я = б8‘/ + 6Pi ~ Эо<Р>*6<Р- *' + Ф. *•bpi + p i б(Р. I-

Исходя из определения напряжения оц и локальной электри­

ческой силы Ef и из того, что

UL =

UL(eijt Pi), причем

(9)

 

atj = dUL/de.r

=

- dUL(dPt,

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

(10)

6Я =

oif be.. + (ф t Е\) ЬРг — э0ф , бсрfi - f Р { бф

Преобразуем следующие выражения:

 

 

Оц Ье1} = (oij bu{\ j Oijt 1 Ьщ,

(! О

 

Ф. * вф. i =

(ф. 16ф). i — Ф, и 6ф,

 

 

Pi бф., =

{Pi бф), i — P t,i бф

и подставим их в бЯ из уравнения

(10). Тогда

(12)

бЯ =

- cr/it} but — {Е\ — ф {) бР{ — (— э0ф н + Р, t) бф +

 

 

 

 

+ [ст/у bui + (— э0ф, i + Pt) бф], г

Следовательно,

(13') -45Hdv = \ [а„,6и, + (fif--9,,)6Р, +

Вв

+(— э0ф, и + Pi, i) бф] dv —

и

— ^ [°li

+ (— э0ф. f + Pj) щ бф] da для тела В

дВ

 

 

 

 

 

(13") — б \ н dv =

? э0ф, и бф do — ^ э0Ф, da для тела В'.

 

В'

В'

дВ

В уравнении (13")

принят во внимание тот факт, что в теле

В' и /эs 0, Р/з= 0.

Подставим (13') и (13") в принцип Га-