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

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

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

 

1.2. Электромагнитное поле

11

Уравнение (10) преобразуем к виду

 

(12)

JL ^ U e dv = - $ n - h d a - jEJdo.

 

ВдВ В

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

В следующем параграфе мы получим обобщение уравне­ ния баланса энергии, учитывающее деформации тела. Со­ гласно равенству (3), В есть солеиоидальный вектор; поэтому его можно представить как ротор некоторого вектора Ао:

(13) B = rotA0.

Однако уравнение (13) не определяет вектор А0 однозначным образом. Поэтому имеем

(14)

В =

rot А,

А =

А0— gradij).

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

(16)

rot (Е + А0) =

0 ,

rot(E + А) = 0 ,

откуда

 

 

 

 

(16)

Е = — А0gradtpo,

Е = — A gradq>.

Функции ф, фо, ф связаны между собой зависимостью

(17)

 

Ф =

фо “Ь dty/dt.

С

учетом

(14) -и

(16) представим уравнения Максвелла

(1)

— (3)

при помощи потенциалов — векторного А и скаляр­

ного ф:

 

 

 

 

(18)

 

rotH =

D +

J,

(19)

 

В =

rot А,

(20)

 

Е = — grad ф — А.

К этим равенствам добавим уравнение Гаусса (4) и опре­ деляющие соотношения (6), (7). Рассмотрим выражение divh = div(EX Н); с учетом соотношения (16) получаем

div h =

— div (grad ф X H) — div (A X H)

=

 

= grad ф • rot H — div (A X H) =

div (ф rot H — A X H).

Принимая во внимание уравнение (1), имеем

(21)

div h = div [ф (D + J) — A X Н].

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

Подставляя (21) в (10) и учитывая (9), имеем

(22)

- Jnfo(D +

J ) - A X H ] d a = $ ( E b +

HB)<*0 + $ E J d o .

 

 

 

в

 

Поэтому

 

 

 

(28)

 

h =

< p ( D - f J ) - A X H .

 

Величина h = Е X

Н

является величиной

однозначной, то­

гда как, согласно соотношению (23), она таковой не является ввиду неоднозначности функций А и ф. Однако выражение

\ п • h da, где h задается формулой (23), однозначно, если

J,dB

только поверхность дВ замкнута!).

Рассмотрим пьезоэлектрическое тело (которое является диэлектриком). В принципе это тело электрически нейтраль­ но: оно содержит одинаковое число положительных и отри­ цательных электрических зарядов и не проводит токов. Вве­ дение диэлектрика в электромагнитное поле приводит к из­ менению последнего. Вследствие этого векторы Е и D ие па­ раллельны и отличаются на вектор поляризации Р. Для пье­ зоэлектриков примем те же упрощения, что и в случае ненамагничивающихся диэлектриков, а именно

(24) J = 0, М = 0, ре= 0.

При этих предположениях уравнения (18) — (20) примут вид

(25)

rot Н =

Ь,

(26)

В =

rot А,

(27)

Е =

— grad ср — А.

Определяющее соотношение (6) остается неизменным, а со­ отношение (7) значительно упростится:

(28)

В = р0Н.

Ввиду того что рв =

0, уравнение (4) становится однородным:

(29)

div D — 0,

а вследствие условия J = 0 вектор Пойнтиига принимает про­

стой вид:

 

(30)

ь = < р Ь - А х н .

Выражение (21) ие изменится, если к h из (23) прибавить функцию rot G, где G — произвольный вектор. Заметим, что

1.3. Уравнение баланса энергии

13

Следующим упрощением будет пренебрежение магнитным членом А в векторе Пойнтинга, так что

(31)

h « фЬ.

Обоснование приведенного выше упрощения (подтвержден­ ного экспериментально) было дано в интересной работе Тир-

стена [52]. Ввиду допущения А « О мы приходим к соот­ ношению

(32)

Е = — grad ф.

Таким образом, мы имеем дело с квазистатическим элек­ трическим полем, описываемым уравнением

(33)

div D =

О

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

 

(34)

D = э 0Е +

Р.

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

Вернемся к уравнению баланса энергии (12). Для пьезо­ электриков при В = О, J = 0 и с учетом (31) имеем

(35)

-Jj J Ue dv =

— $ п • Ьф da

 

В

дВ

ИЛИ

 

 

(36)

- h 5 и , dv =

- \ ф. A dv.

 

в

в

При переходе от выражения (35) к (36) использовалось преобразование Гаусса и принято во внимание уравнение (29). Окончательно, учитывая (32), получили

(37)

± \ u . d o = \ E , D , d v , £, = -Ф .,-

 

В

в

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

1.3. УРАВНЕНИЕ БАЛАНСА ЭНЕРГИИ

Рассмотрим пьезоэлектрик, подвергнутый действию внеш­ них нагрузок и находящийся в электромагнитном поле. К внешним нагрузкам отнесем массовые силы X и поверх­

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

ностные нагрузки р. Эти причины, вызывающие движение тела, являются функциями положения х и времени t.

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

(1)

7Г(Ж +

+

Здесь Ж — кинетическая энергия, выраженная интегралом

(2)

Т = |- J VtVidv, v t =

 

в

— внутренняя энергия, которую представим в виде

(3)

< U = $ U d v ,

 

в

где U — внутренняя энергия (механическая и электромаг­ нитная), отнесенная к единице объема. Через 3 обозначена мощность внешних сил:

(4)

Х(о{ dv - f J PiVida.

 

зв

Наконец, 3) — поток электромагнитной энергии через по­ верхность дВ, ограничивающую область В тела:

(б)

3> = — J

Ditii d a = ^ E iD id v .

 

дВ

В

Принцип сохранения энергии (1) означает, что прирост во времени кинетической и внутренней энергий равен мощности внешних сил и притоку электромагнитной энергии. Прини­ мая во внимание равенства (2) — (5), представим принцип сохранения энергии в следующем виде:

(б)

j viv i + £/) dv = J Xfli dv +

J p p i da -f

J

dv.

В

В

дВ

В

 

Сначала преобразуем поверхностный интеграл в правой ча­ сти равенства (6). Учитывая, что р, = а/,«/, где Л/ — состав­ ляющие нормали п к поверхности дВ, имеем

\ plvi d a = J <уцП/Vi da — ^

/ dv.

ОВ

дВ

В

 

1.3. Уравнение баланса энергии

15

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

 

(7)

J Ф ~ [fa/*. / + Xt — рvt) vt +

j + EiDi)} dv =

0.

 

в

 

 

Так

как соотношение (7) должно быть справедливым для

произвольной области В, то

 

 

(8)

U = (сг^, / + Xi pV{) Vi +

cijiVi, j + EiDi.

 

Таким образом, мы получили локальный принцип сохранения энергии.

Потребуем, чтобы выражение (8) оставалось неизменным при перемещениях тела как жесткого целого [59]. Сначала рассмотрим жесткое равномерное движение тела

(9)

Vi -> Vi +

bt,

где b — произвольный

постоянный

вектор. Подставляя (9) в

(8), приходим к равенству

 

(10)U — (оц, j-\-Xi — рг)/) (Vi 4* b{) 4- OjiVi,} 4- EiDt.

Предположим теперь, что величины р, U, X-h оц и £,• не пре­ терпевают изменений. Вычитая соответствующие стороны уравнений (10) и (8), получим равенство

(11)

(<ty./ + * i - P * i ) * i = 0,

которое должно выполняться при произвольных значениях bi. В результате получаем уравнение движения — первое урав­ нение Коши

(12)

оц, j 4- Xi = piii, x e f i .

Уравнение (12) значительно упрощает равенство (8) ба­ ланса энергии. Тогда имеем

(13)

0 = OjiVi, j E i D i .

Потребуем, чтобы выражение (13) оставалось неизменным при вращении тела как жесткого целого с постоянной угло­ вой скоростью Q. Поэтому полагаем

(14)

v -► v 4- Q X г,

откуда

 

 

(15)

V i->v{ BikiQiXk,

Vi,j-+Vi,j — в/y/Q/.

Вектор

Q — постоянный вектор

вращения. Подставляя (15)

в (13),

получим

 

(16)

О = оц {Vi, I — e^Q/) 4- Е^Ьц

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

Вычитая равенство (13) из (16) и полагая, что величины сг//, Ei, U не претерпевают изменений, приходим к уравнений

(17)

 

==

Поэтому для Qi ф 0

 

(18)

ъцьр7fc = 0

или Оц = ац-

В результате получаем, что тензор напряжений является симметричным. Это — содержание второго уравнения Коши.

Для дальнейших рассуждений определим тензор дефор­ мации 8,-у и тензор вращения со//:

(19)

е// == 7* (u i, / + и ], /)>

 

UJ.i)'

Сложение этих тензоров дает

 

 

(20)

м/. / =

е//+

Qi/-

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

 

 

(21)

vi, j =

&ij +

®i/-

Заметим, что a//d>//= 0, где тензор сг// симметричный, а тен­ зор ©// кососимметричный. Приведенные выше утверждения позволяют представить уравнение баланса энергии с учетом (17) и (21) в виде

(22)

0 = <7//ё// + EiDi.

Отсюда видно, что внутренняя энергия есть функция дефор­ мации 8// и электрического смещения £)/. Внутренняя энергия является функцией состояния, а ее производная — полной производной; следовательно,

(23)

0

Сравнивая уравнения (22) и (23), имеем

(2 4 )

Эти равенства выполняются для произвольных ёц и £>/, если

(25)

_

dU

„ _ dU

° Ч ~

дгц '

t l ~~ dDt ’

 

Соотношения (25) пригодятся нам для введения определяю­ щих уравнений в предположении, что U = f/(e//, £),-), т. е. что независимыми функциями являются е// и Z)/, а зависимы­ ми — функции ail и Ei.

Для дальнейших рассуждений выгодней использовать та­ кие определяющие уравнения, в которых величины ац, Di будут зависеть от е// и Е(. С этой целью введем электриче­

 

1.4. Определяющие соотношения

17

скую энтальпию, определенную следующим образом:

 

(26)

Н = U — EtD{.

 

Исключая из

уравнения

(22) и

(26)' функцию Я,

получим

(27)

Я =

ч„ё„ -

О Д .

 

Здесь П а й

(в/,-, Я,). Следовательно,

 

да)

Сравнивая уравнения (27) и (28), получаем соотношение

(29)

+

=

которое должно выполняться при произвольных величинах в//, Ё{. Отсюда

/ол\

ая

п

ая

^

аЧ ~ de.f

‘Di==

дЯ, *

1.4.ОПРЕДЕЛЯЮЩИЕ СООТНОШЕНИЯ

Разложим электрическую энтальпию Я (ец, В/) в ряд Маклорена в окрестности естественного состояния, т. е. состоя­ ния, для которого е(/- = 0, Ei = 0:

(1) Я (в„, Ei) = Я (0, 0) + -дНд^ °) е„ +

Et +

+ * [ % № * * + ^

+ •

Значение Я (0,0) есть электрическая энтальпия для есте­ ственного состояния; без потери общности можно принять ее равной нулю. Введем обозначения

 

 

аЯ(0,0)

~

ф

ая (0, Q) _

Ь{,

 

(2)

 

двЦ

дЕ1

~

 

 

а2я (о, о)

_

а2я (о, о)

_

 

а2я (о, о)

Эц.

 

dBu dBki

~ ClIkh

dEkdBn

6кф

dEt dEj

 

 

Таким образом, выражение (1) запишется в виде

 

(3)

Я {etj, Ei)=* ацВц +

btEt +

ll2ciikietj4 i ~~ekn^kBn ~

 

 

 

 

 

 

 

 

~

....

Желая получить линейные определяющие соотношения, огра­ ничимся в ряде Маклорена лишь линейными и квадратич­ ными членами, Из соотношений оц = дЯ/де,у, D* = —dH/dEf

18

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

получим линейные определяющие соотношения

(4)

Оц = Cijkfiki ~~ ещЕ&

(5)

D{= eikizki + BikEk.

В этих зависимостях принято, что ац = 0, 6,- = 0, так как мы предполагаем, что в теле отсутствуют начальные деформации

и начальное электрическое поле. В соотношениях

(4), (5)

Cijki — упругие

жесткости (при

постоянных

£,•),

виц — пьезо­

электрические

постоянные,

а э

— диэлектрические

постоян­

ные (при постоянных е*/). Заметим, что из

уравнений

(3) —

(5) можно вывести следующие соотношения:

 

 

 

 

д 2Н

д 2Н

 

ИЛИ

дО ц

d ° k i

 

 

 

д е ц d ekl

d h l д г Ц

d ekl -

д в 1}

 

 

 

 

 

 

д 2Н

д 2Н

 

и л и

dD t _ d D j

 

 

 

д Е { dEj

д Е } d E t

~ dE J ~

d E t

*

 

 

 

 

 

Эти соотношения приводят к следующим результатам

(сим­

метрия тензора сцы усматривается уже в формуле

(4)):

(7)

 

c ijki = c ki ijt

э,7 = Эц.

 

 

 

 

Из того факта, что тензоры ац и вц являются тензорами

симметричными, вытекают

следующие соотношения:

 

(8)

Cijki ~ CjiM =

Cijik, в^ц = e^ji.

 

 

 

Заметим, что тензор эсимметричен, а трехиндексный тен­ зор ещ симметричен только относительно индексов i, /. Со­ отношения симметрии (7)i и (8)ь (8)2 уменьшают число постоянных от 81 до 21. В случае наиболее сложного триклинного кристалла имеем 21 постоянную с-цм, 18 пьезоэлектри­ ческих постоянных виц и 6 диэлектрических постоянных э«•/. Рассмотрим отдельно случай изотропного тела: Величины Сцы, Skij, эц должны быть в этом случае изотропными тен­ зорами

^

Сцк1 = И- (bik^ji +

+ кЬцЬк1,

 

Cklj == &kijC> Эц =

эбц-

Здесь р, %— постоянные Ламе, гци — кососимметричный тен­ зор Риччи. Подставляя (9) в (4), имеем

(10)

 

оц =

2цец +

"kbn&kk,

 

так

как член

BkijeEk = 0

ввиду

симметрии тензора

ац. Из

соотношения

(5) получаем

 

 

 

 

(11)

 

0 / =

б/лэ£й =

э£'/,

 

так

как член

ешеА/ = cziklzki =

0

ввиду симметрии

тензора

 

В результате получаем, что в изотропном тел$ связности

1.5. Дифференциальные уравнения теории пьезоэлектричества

19

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

Возвращаясь к выражению электрической энтальпии Н —

— U — EiDi, запишем его в виде

(12) £/ = # + £*£,.

Подставляя выражения (3) и (5) в (12), приходим к вы­ ражению для внутренней энергии V:

(13) U — ll2CijklQijeltl + 'l&ijEiEj-

Так как U — неотрицательная скалярная величина, правая часть выражения (13) должна быть положительно опреде­ ленной квадратичной формой, что обусловливает устойчи­ вость решения.

1.5. ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ ТЕОРИИ ПЬЕЗОЭЛЕКТРИЧЕСТВА

Приведем сводку основных соотношений и уравнений тео­ рии пьезоэлектричества. Имеем два основных уравнения — уравнение движения и уравнение электрического поля:

( 1)

°7(, / 4*

= P^i,

0.

(2)

Di. t =

x e S , t >

О,

 

К этим уравнениям добавляются определяющие соотношения

(3)

вц — cijkl4l ~~ekijEk>

 

 

(4)

Di =

Giki&ki Н~

E/t

ф, ft,

и определение тензора деформаций

 

 

(5)

 

Бц = Уо(^/. / ~f~ ^/. *)•

Подставляя

(3) и

(4) в равенства

(1)

и (2), приходим к си­

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

(6)

cijkiuk,ti + eki]4>k, i +

— Plli>

(7)

eikiuk, и ~~ э;аФ.ы 0-

Эти уравнения связаны между собой. К ним следует доба­ вить граничные условия. Пусть поверхность дВ, ограничи­ вающая тело, состоит из частей Гi и Гг, так что дВ = Ti U Гг, причем Ti П Гг = 0. Пусть на части поверхности Ti заданы перемещения (У/, а на части Г2— нагрузки В этом случае будем иметь дело с граничными условиями

(8)

Ut = U{ (x, t) на Г,, p{ = o}i(xt /)п/(х) на Г2.

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

Пусть теперь на поверхности Гз задан электрический потен­ циал ф, а на Г4— поверхностный заряд а, причем дБ = = Гз U Г4, Гз П Г4= 0. Тогда будем иметь дело с граничными условиями

(9)

ф= ф (х, t) на. Г3> Dknk = — а на Г4.

Равенства (6) являются уравнениями движения; уравне­ ние (7) описывает квазистатическое электрическое поле. К граничным условиям (8) следует еще добавить начальные условия для перемещений:

(10) и, (х, 0) = /*(*). щ{х, 0) = gi(x), X G B, t 0.

Зная решения системы уравнений (б) и (7) — функции щ, ф, определим последовательно поле Ek по формуле £*. = —ф, л и деформации ъц по формуле (5). Напряжения о,-/ и электри­ ческое смещение £>,• определим из определяющих уравнений

(3) и (4). Знание функций £,• и D,- позволяет тогда опреде­ лить вектор поляризации Р, по формуле

(И)

Pt = Dl - 9 0Ei.

Здесь Эо — электрическая проницаемость в вакууме. Рассмотрим частные случаи равенств (6) и (7). Если пье­

зоэлектрический эффект в теле не проявляется, то уравне­ ние (6) принимает вид

( 1 2 )

c t j k i u k, ц + X i —

С другой стороны, рассматривая тело как недеформируемое, имеем

(13)

э^ф, ij = 0.

Уравнение (12) следует дополнить граничными (8) и началь­ ными (10) условиями. Уравнение (13) дополняется гранич­ ным условием (9).

1.6. ПРИНЦИП ВИРТУАЛЬНЫХ РАБОТ. ЕДИНСТВЕННОСТЬ РЕШЕНИЙ

Отправным пунктом наших рассуждений будет принцип виртуальных работ Лагранжа для деформируемого тела

( 1 )

J {Xi piii)buidv +

J р^Ьщ d a =

t ацЬъцйу.

 

в

зв

в

Здесь Ьщ— виртуальное приращение перемещения щ. Пред­ положим, что виртуальные приращения Ьщ являются произ­ вольными дифференцируемыми функциями класса С<1>. Пред­ положим еще, что приращения Ьщ согласованы с условиями,