книги / Электромагнитные эффекты в твердых телах
..pdf
|
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) |
Е = — А0— gradtpo, |
Е = — 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>. Пред положим еще, что приращения Ьщ согласованы с условиями,