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

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

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

 

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) — заданные поверхностные заряды на дВ. Там, где заданы перемещения, = 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

, а также а, и являются

материальными

постоянными.

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

принимают упрощенный вид:

 

 

(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+

^44677) 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 — ЗД>.