книги / Электромагнитные эффекты в твердых телах
..pdf
|
2.1. Введение |
81 |
мильтона (7). Согласно равенствам |
|
|
|
и |
|
( И ) |
^ dt ^ piiibUidv, |
в |
получим |
и в |
|
|
|
|
(1б) $ [(a/<f , + |
Xt — Р«,) Ьи. + (£{■ — ф , —£5) 6Pt + |
|
в |
|
|
+ (— э0ф. и + Pi. i) бф] dv — J э0ф. а бфdv + |
||
|
В' |
|
+ |
\ [(Pi — °iinJ)bul + (90\q>,i\ — Pi)ni b(p]da = 0. |
|
|
дВ |
|
Здесь |<р,,|— скачок функции <p,t- на поверхности дВ. Заметим, что в объемный интеграл входят виртуальные приращения Ьш, 6Pi, бф, а в поверхностный интеграл — только виртуаль ные члены бщ и бф. Ввиду произвольности виртуальных при ращений из (15) получаем следующие уравнения Эйлера;
(16) |
— Pdi |
О, |
-х е В, |
|
(17) |
|
|
О, |
|
(18) |
Эоф. и + |
£/,< = |
0, J |
|
(19) |
Ф, и = 0 |
для |
х е |
В'. |
Здесь уравнение (16) представляет собой уравнение движе
ния, а |
(18) идентично |
уравнению |
электромагнитного поля |
|
Di, i = |
0. Подстановка в Д , ,• = 0 соотношений Д = |
эоД -f- Pi |
||
и Ei = |
—ф,; непосредственно приводит к уравнению |
(18). |
||
Введенный Тупином «баланс межмолекулярных сил» пред |
||||
ставлен уравнением (17). Заметим, |
что для Е] — 0 и £ f = |
|||
= —Ei уравнение (17) |
переходит в уравнение £,• == —ф. |
|||
К уравнениям поля |
(16) — (19) присоединяются так назы |
ваемые естественные краевые условия, получаемые из поверх ностного интеграла в выражении (15):
(2°) |
чу1я /= р (, |
(21) |
щ [ - э 0|<р.,Ц-Р,]=О, XSOB- |
Условие (20) представляет заданные нагрузки на дВ, а уело* вне (21) — заданные поверхностные заряды на дВ. Там, где заданы перемещения, 6ш = 0, а там, где задан электрический потенциал, бф = 0. Заметим, что краевое условие не связано с уравнением (17). Это следует из того, что в поверхностном интеграле, входящем в уравнение баланса (15), отсутствует член бPi.
82 Гл. 2. Теория пьезоэлектричества Тупина и Миндлина
Зададим внутреннюю энергию |
UL = |
Pi) в виде |
|
квадратичной формы |
|
|
|
(22) |
№ = ‘/2<f(/iIe„ew + |
+ |
fu N Pk. |
Здесь cfjkl — коэффициенты упругой податливости при по стоянных Pf, aztj — коэффициенты поляризации при постоян
ной деформации в,у. Используя соотношения (9), приходим к следующим определяющим уравнениям:
(23) |
|
0"jy === |
-{- ffcijPfc, |
|
(24) |
- |
ELi = f m H i + |
aliP r |
|
Здесь |
|
|
|
|
^ijk l~ c]ikb |
cifkl ~ |
cijlk> |
ciikl = |
Ckllf* fkij = fk]i> a] l ==aiJ‘ |
В этих соотношениях напряжение а»/ и локальная электриче ская сила Ef зависят от деформации и вектора поляризации.
Дифференциальные уравнения (16) — (19) вместе с краевыми условиями (20), (21) и определяющими соотношениями (23),
(24) |
составляют линейный |
вариант теории |
пьезоэлектриче |
ства Тупина. |
|
|
|
Миндлин [32] получил соотношения между постоянными |
|||
cfjki> |
ekij> э*/ и постоянными |
tfJkr fkij, а?ц |
и показал, что |
дифференциальные уравнения пьезоэлектричества Тупина пе реходят в уравнения классической теории Фойгта.
Рассмотрим еще основные энергетические теоремы, кото рые представляют собой различные формы принципа сохра нения энергии. Отправным пунктом рассуждений послужит принцип виртуальных работ Лагранжа
(25') |
^ {Xi — pui)6uidv-\- ^ pibuida — ^oifbzijdv. |
||
|
в |
дВ |
в |
Подставляя в (25') определяющее соотношение |
|||
(26) |
Оц = |
Cijkieki + |
fki}Pk> |
приходим к уравнению |
|
|
|
(25") |
[ (X* — piii)buidv-\- |
\ Pt Ьщ da = |
|
|
в |
дВ |
|
|
|
|
~ ^ (cijki4i + fkijPk) bei{dv. |
|
|
|
в |
В дальнейшем предположим, что |
виртуальные перемещения |
||
образуются действительными, |
|
||
(27) |
bui — ijfdi = vt dt, |
dt =*= ёц dt. |
2.1. Введение |
83 |
Подставляя эти соотношения в равенство (25"), получим
(25'") |
± (Ж + |
Ж) + ful \ ё„Р* dv = |
J p,v, da + $ X,v, dv, |
|
|
|
|
d B |
В . |
где |
Ж = |
у J WiVi dv, Ж = у |
J ст гИък1 dv |
|
|
|
в |
в |
|
{Ж — кинетическая энергия, Ж — работа |
деформации). По |
множим определяющее соотношение (24) на Р,-, полагая, что
£9 = |
0, Е^ = — Ег |
Интегрируя по области тела, получим |
||
(28) |
\ fkiPifik dv + \ auPtPi dv = |
J EiPi d v . |
|
|
|
в |
в |
в |
|
Складывая уравнения (25'") и (28), находим |
|
|||
(29) |
+ |
+ |
|
|
|
|
= ( ptv t da + |
^ Xflidv + |
\ EtP{dv. |
|
|
д В |
В |
В |
Здесь введено обозначение
^ = Т Sanp ip i dv-
В
Заметим, что, согласно представлению (22),
(30) |
Ж + |
9 + |
fkit J et,Pk.dv = <UL = \ U L dv. |
||
|
|
|
|
в |
в |
Учитывая (30), представим выражение |
(29) в виде |
||||
(31) |
- ± W + |
<UL) = |
\ p tV td a + ^ X f lid v + ^ E iP id v . |
||
|
|
|
д В |
В |
В |
Воспользуемся, наконец, уравнением электрического поля
(32) |
Dh. * = 0. |
Помножим это уравнение на ф и проинтегрируем по области тела. После простых преобразований, используя формулу Гаусса, имеем
(33) |
^ b tynt da — ^ £>i<pt {dv = 0. |
|
|
дВ |
В |
Наконец, |
учитывая |
зависимости £>«• = Эо£/ + Pi, £« = —Ф,;* |
приведем уравнение (33) к виду |
||
(34) |
^ P tEi dv — — J b {<pnt da —э0 J Ф. <Ф.i dv. |
в |
дВ |
в |
84 Гл. 2. Теория пьезоэлектричества Тупика и Миндлина
Исключая из |
уравнений |
(31) |
и (34) |
общий интеграл |
||
EiPidv, находим, что |
|
|
|
|
||
(35) |
(Ж + |
<U) = $ |
PiVi da + |
J Xtvt dv - |
$ D m i da, |
|
|
|
д В |
|
|
В |
д В |
где |
|
<U= |
<UL+ |
у J э0ф, ,ф. I dv |
|
|
|
|
|
|
в |
|
|
— внутренняя интегральная энергия. Основная энергетическая теорема (35), полученная здесь в рамках теории Тупина, на ходит себе аналог в форме теоремы (9) из § 1.6 в теории Фойгта.
Выведем теперь теорему взаимности работ для рассматри ваемой теории Тупина. Отправным пунктом рассуждений бу дут уравнения двух систем причин и следствий
(36) |
ад, / + |
Xi = |
О, |
(37) |
а'„ , + |
*; = |
0. |
Ограничимся статическими задачами. Умножая уравнение (36) на и\, а уравнение (37) на щ и вычитая затем одно из другого, проинтегрируем результат по области тела. После простых преобразований получаем известные уравнения
(38) |
$ (р.и\ ~ p\ut) da + |
fa.u'i—X,lui) d v = ^ ( o iJe.'lf— <r'tJsi/) dv. |
|
|
д В |
В |
В |
Теперь воспользуемся определяющим соотношением для оц:
|
|
Gli = |
Cllkfikl + fkljPk |
|
и аналогичным выражением для |
Подставляя эти соотно |
|||
шения в (38), получим |
|
|
||
(39) |
^ |
{Piu'i — p'iui) da + |
^ {Хьи\ — Х\иг) dv = |
|
|
ЗВ |
в |
|
|
|
|
|
= |
$ htj i^ijPk ~~ вцР'ь) dv - |
|
|
|
|
в |
Исподьзуем теперь равенство (17) |
(в предположении, что |
|||
E°t = |
6) |
для двух систем |
|
|
(40) |
|
Я ^ -Ф ^ -О , |
|
|
(41) |
|
B |
f - Ф ^ - 0 . |
|
Умножим уравнение (40) на Р\, а уравнение (41) на Pi. За тем вычтем одно из другого и проинтегрируем по области
2.1. Введение |
85 |
тела: |
|
|
|
|
|
(42) |
|
{ (£ f Р\ - E’Lp,) dv - 5 (Ч> |
- ф;tP,) dv = |
0. |
|
|
|
в |
в |
|
|
В первый интеграл подставим определяющее соотношение |
|||||
|
|
Щ = |
ikie kl "Ь a ijP] |
|
|
и аналогичное соотношение для E'tL. |
Таким образом, прихо |
||||
дим к равенству |
|
|
|
||
(43) |
|
$ /*„ («;,£* - |
dv - 5 (Ф. ,Р' - <?',Р() dv. |
||
|
|
В |
в |
|
|
Из уравнений (39) и (43) |
получаем следующий вид теоремы |
||||
взаимности работ: |
|
|
|
||
(44) |
$ |
{рр\ — pfa) da + |
|
dv= $ (ф |
tPt) dv. |
|
дВ |
В |
|
В |
|
Наконец, используем уравнение электромагнитного поля (18) для обеих рассматриваемых систем причин и следствий
(45) |
э0Ф. ц + Pi, t — |
|
(46) |
эоФ;а |
p 'i,i~ О* |
Первое уравнение умножим |
на ф', второе — на ф. После ин |
|
тегрирования по области тела имеем |
(47) |
J (p i. У ~ |
К У d v |
— |
эо \ (Ф. нФ' — ф' ^£ф)d v |
** 0. |
|
|
В |
|
|
В |
|
|
После |
простых преобразований уравнение (47) приведем квиду |
|||||
(48) |
$ (/> ; , - Р > .,) dv - |
J [(э ^ |
f - Pt) Ф' - |
|
||
|
В |
|
дВ |
- ( э иф'г-р ;)ф ]/г <dar= 0, |
||
или |
|
|
|
|
||
|
|
|
|
|
|
|
(49) |
J (PfФ,'i — Р У i) сГо — J |
(£>У — Z>y) nt da. |
||||
|
В |
|
|
дВ |
|
|
Подставляя (49) |
в (44), приходим к окончательному виду тео |
|||||
ремы взаимности работ: |
|
|
|
|
||
(60) |
^ р У da + |
^ ХЖ dv -f |
^ Dlq>'nt da ■= |
|
||
ов |
в |
|
ов |
|
|
|
|
|
= |
jj |
p'tu,i da + J X\ut dv + |
J D'i<fni da. |
|
|
|
|
OB |
|
В |
dB |
Эта теорема идентична теореме, выведенной в рамках теории Оойгта (ср. с уравнением (21) из § 1.8 при со->0).
86 |
Гл. 2. Теория пьезоэлектричества Тупина и Миндлина |
2.2.УЧЕТ ГРАДИЕНТА ПОЛЯРИЗАЦИИИ
ВЛИНЕЙНОЙ ТЕОРИИ ПЬЕЗОЭЛЕКТРИЧЕСТВА
Экспериментальные исследования показывают, что пьезо электрический эффект может проявиться и в центрально-сим метричных кристаллах. Проявляются, кстати, и другие ано малии. В тонких слоях диэлектрика, помещенного между металлическими электродами, поведение электрического по тенциала и поляризации отличается от поведения, следующего из классической теории пьезоэлектричества. Если предполо жить, что ф = ±const на границах слоя, то из классической теории получаем ненапряженное состояние, линейный харак тер по толщине слоя электрического потенциала и . постоян ную поляризацию. Однако эксперименты, произведенные Ми дом [29, 30], обнаруживают нелинейное поведение функций ф и Pi и появление напряжений в слое диэлектрика. Эти недо статки классической теории Фойгта побудили изыскивать ее обобщение. Названную аномалию хорошо объясняет теория Миндлина [31, 34], в которой дополнительно учитывается влияние грааиента поляризации на функции а.., ф, Ejr. При
обобщении теории пьезоэлектричества Миндлин исходил из линейной теории диэлектрика Тупина, изложенной в § 2.1. Развитие этой теории основано на введении в выражение эн тальпии (уравнение (6) § 2.1) градиента поляризации:
(1) |
|
H = UL( &ц, P h Pt, /) — 72э0ф. i<P. i + |
Ф. iPi- |
|
|
|||
Следовательно, |
|
|
|
|
|
|
||
(2) |
ЬН = |
dUL |
dUL |
6Pi + |
dUl |
|
|
|
дец |
дР, |
дР |
|
|
|
|||
|
|
|
|
t- i |
|
|
|
|
|
|
|
|
|
ЭоФ. i ^Ф , * - f - Ф. г |
+ |
Pi бф, i. |
|
Вводя определение |
тензора |
Ец = |
dUL/dPji /, |
получим |
из (2) |
|||
(4 . |
ЬН = <F„ 6в„ + |
(ф ,-£{■) ЬР,+Е„ ЬР, ,- э „ ф , бф ,+ |
/>, бф |
Дальнейший ход действий аналогичен предложенному в пре дыдущем параграфе. Принцип Гамильтона
ti |
k |
(б ) 6 ^dt |
\ ( K - H ) d v + |
tt |
в* |
. bui + Е\ЬР,) civ ~f-
“Г ^ Pi but da |
0 |
dB |
|
2.2. Учет градиента поляризации в линейной теории |
87 |
приводит к уравнению
и
иВ
идв
Заметим, что здесь появляются два новых члена: £/«. /б-Р* в объемном интеграле и £/,/г/бЛ в поверхностном интеграле. Этот последний приводит в дальнейшем к новому краевому условию. С учетом произвольности виртуальных приращений dm, бPi, бф получим следующие уравнения:
(7) |
|
°7л / ~Ь %i= |
P^f ' |
||
(8) |
£ / и + |
Е \ — ф t. + |
£9 = |
0 I для х е й , |
|
(9) |
~ |
э0Ф, и + |
Pi, i = |
0 |
|
( 10) |
|
Ф, а = |
0 |
для |
х е й ' |
сестественными краевыми условиями
( И ) |
oiln] = |
pt |
* для ХЕЙЙ. |
( 12) |
Ejitlj==0 |
||
(13) |
(— э0|ф, t\ + Pi)nt = |
0 |
J |
По отношению к теории, не учитывающей градиента поляри зации, изменению подверглось втррое уравнение Эйлера, ко торое дополнено членом Ец, /. Добавилось новое краевое усло вие (12). Возможны также и другие краевые условия на дБ. Там* где на границе заданы перемещения, следует принять
Ьщ = 0; там |
же, |
где задана поляризация, следует принять |
|
dPi = 0; |
наконец, |
там, где задан электрический заряд, пола |
|
гаем бф = |
0. |
f/L(e,7, Pit Pij) в виде |
|
Зададим |
Разложим |
теперь UL в окрестности естественного состояния, |
в котором |
деформация и поляризация равны нулю, а несим- |
11 Обозначения в, Р, Q относятся к фиксированной деформации, поля ризации и градиенту поляризации. В дальнейшем эти обозначения будем опускать.
88 Г л. 2. Теория пьезоэлектричества Ту пина и М инд лина __________
метричный тензор Ец принимает значение Ъц. Учитывая со отношения
п Ч ч |
„ _ dU |
p L _____dU^_ |
Е |
|
dUL |
|
(15) |
|
t i ~ |
дР{ |
’ |
‘7 |
dPjti ’ |
получим следующие определяющие соотношения:
(16) |
= |
Ciikfikl + fkijPk + dkiijPi, k, |
(17) |
E f — |
f jk iBki “1 a jk P k i j k i p i, k » |
(18)E„ = b%+ d m 4 t + jstlPk + bm P
Анализ этих соотношений приводит к утверждению, что пьезоэлектрический эффект проявляется и в центрально-сим метричном теле. Действительно, в этом случае равны нулю тензоры fkij, jjki, поскольку тензоры нечетного порядка в цен трально-симметричных телах должны отсутствовать, но от личны от нуля оставшиеся тензоры четного порядка йьщ и bijki. Заметим далее, что тензор Ец несимметричен. В есте ственном состоянии В// = О, Р,*= О, Р/,/ = 0, и только тензор Ец не исчезает в естественном состоянии — его значение равно
Эти последние величины можно трактовать как началь ные данные, введенные в тело в естественном состоянии.
Подстановка соотношений (16) — (18) в уравнения (7) —
(9)приводит к системе семи дифференциальных уравнений
относительно |
Р,, ср. |
Заметим, что величины bjj входят в решение через крае вое условие (12). Из системы уравнений (7) — (9) видно, что введение градиента поляризации не повышает порядок диф ференциальных уравнений.
2.3.ИЗОТРОПНЫЕ ДИЭЛЕКТРИКИ
Рассмотрим центрально-симметричный и изотропный ди электрик. В этом случае энергия UL принимает вид
(1) UL = bQP]t f + Va^P/Py+'^ioPi. iP/. /+7а (&44+M Pi. iPf. <+
+Va (*44—b-п) Pj, iPi, / + 1/2Ci28ii8//-fc44ef/efy+ cf12pf <8у/+ 2 d44P*, ,ew,
где |
Ci2, C44, d.12, |
612, ^44» 677 |
, а также а, и b° являются |
|||
материальными |
постоянными. |
Определяющие соотношения |
||||
принимают упрощенный вид: |
|
|
||||
(2) |
аЦ = |
cl2uk, к6ц + |
с44 (Щ, i+uj, t) + d l2Pk' kbtl+ d AA(PL,+P ,. ,), |
|||
(3) |
Ец = |
d^k, кЬц "j~ ^44(ui, 1 ~b Uj. y) + |
bi2Pk кЬц -j- ' |
|||
„ |
|
D |
+ |
(Л. ( + |
P,,t) + |
b„ <PM - P t, ,) + № ,„ |
|
2.8. Изотропные диэлектрики |
89 |
|||
Подстановка |
соотношений |
(2) — (4) |
в дифференциальные |
||
уравнения |
|
|
|
|
|
(6) |
U/i, / + Xi = |
piii, |
|
|
|
(6) |
+ |
f + |
2;° = |
0, |
|
(7) |
— э0Ф, a + |
Pi,i = Pe |
|
|
приводит к следующей системе трех связанных уравнений:
(8) |
c44V2u + |
(ci2+ |
сАА) grad div и + |
d44V2P + |
|
|
|
|
|
+ (rf|2 + ^44) grad div P + |
X = |
pii, |
|
(9) |
d44V2u + |
(^12+ |
^44) grad div u 4* (^44"I" ^77) V2P + |
|
|
|
|
+ |
(&12+ |
^44— 677) grad div P —aP — grad <p + |
E° = |
0, |
|
(10) |
|
— 30V2<p + divP = pe |
для x e B , |
|
|
|
(11) |
|
|
У2ф = 0 для |
x e B ', |
|
|
в которые в качестве неизвестных входят функции: переме щение и, поляризация Р и электрический потенциал <р. К диф ференциальным уравнениям (8) — (10) добавим краевые усло вия
(12) О ц П ^ р р |
Ejtnf = |
0, (э0 (ф+, — Ф^) + Р{) tit = 0 |
на дВ. |
Вернемся к |
системе |
(8) — (10). Это очень сложная |
систе |
ма связанных уравнений, которая в таком виде трудна для
решения. Ясно, что действие источников — массовых сил X, |
|
внешнего электрического поля Е° и электрических |
зарядов |
ре — вызывает в теле деформированное состояние, |
которому |
сопутствуют поляризация и электрический потенциал. Поля эти вызваны также краевыми и начальными условиями.
Структура дифференциальных уравнений (8) — (10) ука зывает на то, что на границе можно задавать как поляриза цию, так и электрический потенциал, что невозможно в клас сической теории пьезоэлектричества. Значительное упрощение уравнений мы получим путем разложения представленных в
них векторов на потенциальную и соленоидальную части: |
|
||||||
(13) |
u = |
grad ф + rot Н, |
Р = |
grad %-j- rot К, |
|
|
|
|
|
|
|
div Н = |
0, |
divK = |
0, |
(14) |
X « |
р (grad ft + rot r\), |
E° = |
grad т + rot £, |
|
|
|
|
|
|
|
div л = |
0, |
div£ = |
0. |
Эти действия аналогичны тем, которые применяются в клас сической динамической теории упругости. Подстановка разло жений векторов (13) и (14) в уравнения (8) — (10) приводит
90 |
Гл. 2. Теория пьезоэлектричества Тупина |
и Миндлина |
|
к системе независимых уравнений |
|
|
|
(15) |
(cuV2 — рд?)Ч> + diy-% = |
— pft, |
|
(16) |
+ (buV2 — а)% — ф = |
— т, |
|
(17) |
V2 (X — э0ф) = |
рв, |
|
(18) |
(c44V2 - р5?) Н + d44V2K = |
- |
pri, |
(19) |
d4iV2H + ((b44+ b71) V2 - а) К = |
- |
g. |
Здесь |
введены обозначения: с\ \ = |
с\2 + 2с44, d\ \ — dn + 2d44, |
bn = |
Ьц -f 2644. Следует еще дать |
интегральные представле |
ния для источников Ь, т, ц, £. Подвергнем выражение (14) i операции дивергенции и ротора. Получим уравнения Пуассона
(20) pV2^ = Xisi, PV2T); = — ei[k^k. /••
Решение уравнений (20) приводит к интегралам [41, с. 181]
,21) #(x) = - J r ^ |
i ( x ') ^ 7 T O L T dl/(x'), |
V |
1 |
(22) Т]г(х) = 4^ - J 8ijkXli* )-g^ Х(хГхТ dV
V |
Л |
R (х, х') = |
| х — х' |. |
Аналогичные формулы получим для источников т, £:
(23)x ( * W —
|
|
V |
1 |
(24) |
{,(«).— |
« « |
д а -i r |
|
|
у |
л |
|
|
R (х, х') = | х — х' |. |
|
Решим |
систему |
уравнений |
(15) — (17), формально исключая |
последовательно функции ф и %. В результате получим урав нения
(25) |
[ > ' - У |
2( 1 |
+ ^ |
) + |
^ ] ф = |
|
|
|
|
|
|
: КЬ<У- - й) о* - |
+ э0-р )], |
||
(26) |
[ ^ - * |
( i + |
j £ a s ) + a j - £ ] x - |
|
|||
|
|
|
= |
- |
a t |
[с" п ? + |
э» ,р) " d"pv2*]- |
(27) |
|
|
V2P = |
p„ |
Р = |
X — ЗД>. |
|