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

книги / Механика сплошной среды. Т.4 Основы механики твёрдых сред

.pdf
Скачиваний:
25
Добавлен:
19.11.2023
Размер:
41.7 Mб
Скачать

§ 2.13. Упругие среды с электромагнитными свойствами

261

Б. Закон Ома и закон Фурье

Ксистеме (2.13.17) присоединяют обобщенный закон Ома (т. 2, (3.14.55))

иобобщенный закон Фурье (т. 2, (3.14.56)), которые для случая малых деформаций имеют следующий вид:

j = L 12V0 + RT1

е,

—q = Л • V0 + 6*L 12 • е,

(2.13.18)

где R — тензор электропроводности; Л — тензор теплопроводности (оба

являются симметричными и

положительно-определенными); L 12

тензор

Томсона, который также полагаем симметричным. Часто эффектом Томсона пренебрегают, тогда полагают, что L 12 = 0. Далее будем применять это допущение.

В. Основное термодинамическое тождество для электромагнитных сред с малыми деформациями

Кроме (2.13.17) и (2.13.18) необходимо рассмотреть определяющие соот­ ношения, связывающие тензоры напряжений сг и деформации е с векторами b, d, е, h и температурой в. Общий вид этих соотношений, вытекающий из законов термодинамики, для сред с произвольными конечными деформациями был установлен в т. 2, п. 3.14.8 и определяется формулами (т. 2, (3.14.78)), а линейные модели упругих электромагнитных сред описываются соотно­

шениями (т. 2, (3.14.84)). Для сред с

малыми деформациями можно

было

бы просто линеаризовать соотношения

(т.2,

(3.14.78) или (3.14.84)), однако

удобнее поступить иначе — так, как это

былосделано в п. 2.1.5, —

выве­

сти определяющие соотношения из основного термодинамического тождества (т. 2, (3.14.25)), которое в случае малых деформаций имеет следующий вид:

ре* — рву — Э * ё b rh —е р —е • j + + ш** = 0.

(2.13.19)

Здесь учтено, что при малых деформациях

 

(п)

_

Н

~

о

о

0

о

(2.13.20)

Т

= <т,

С = в ,

b = b,

е = е,m = m,

р = р,

j = j ,

(п)

 

 

 

 

 

 

 

_

где Т

электромагнитный тензор напряжений Коши (т. 2, (3.14.12)); Э —

симметричный электромагнитный тензор напряжений, имеющий, согласно

(т. 2, (3.14.17а)), следующий

вид:

 

^ = cr + ^(e(g)p

+ p(g)e + b(g)m + m(g)b).

(2.13.21)

Преобразуем входящие в (2.13.19) электромагнитные величины с учетом (2.13.16):

b*m = b ' b - b - h = ( ^ ) - b h,

(2.13.22)

е • р = е • d —be - ё = (e • d —

—d • ё,

(2.13.23)

262

Глава 2. Упругие среды с малыми деформациями

тогда, принимая во внимание, что е определяется по формуле (т. 2, (3.14.3))

ре = °ре + m • b = °ре + |b|2 —b • h, (2.13.24)

после подстановки (2.13.23) и (2.13.24) в (2.13.19) приводим ОТТ к следую­ щему виду:

|2

12 , ,

\ • о л .

_

.

 

V

+ 2 +

2-----b • h

- е • d j

рву -

а

е+

 

 

 

 

+ b • h + d • ё -

е • j + w* + ш** = 0.

(2.13.25)

Введем обозначения

 

 

 

 

 

 

р ё =

°ре +

---- b • h —е • d,

 

ф = ё вг],

(2.13.26)

тогда формула (2.13.25) принимает вид

 

 

 

 

рф* + руб —< r - e

+ b h

+ d e —е • j + гс* + гс** = 0.

(2.13.27)

Г. Общий вид определяющих соотношений для упругих электромаг­ нитных сред с малыми деформациями

Рассмотрим идеально-упругие электромагнитные среды, для которых ак­

тивные переменные Л = { ф , a, b, d,

е • j —гс* —гс**} являются

функциями

от реактивных переменных 1Z = {в, 0, h, е}:

 

Л =

/(7г),

(2.13.28)

в частности,

 

 

ф = ф(е, в, h, е).

(2.13.29)

Подставляя (2.13.29) в ОТТ (2.13.27), приходим к следующим определя­

ющим соотношениям:

 

 

 

 

 

дф

а = р

дф

, Ъ = —р

дф

дф

w = е • j, ш* = 0.

= —

де

,

d = р-^-

V =

~дв9

 

dh

де

(2.13.30)

Д. Линейно-упругие электромагнитные среды

 

В

модели

линейно-упругих

электромагнитных

сред потенциал ф

(2.13.29) полагают квадратичной функцией векторных и тензорных аргумен­ тов:

РФ = рФвР) + ф

■■

• • е +

■■(ЗЛ -

• /х • h -

• х

• е-

h u e /З- e + q / i h

+ qe-e —е •

Зл/г

^

и Зт

• • в, (2.13.31)

6Ше • • е h • 6Шр

 

 

в

в

 

 

 

 

 

'Фв =

Фо +

cv{ 9 ')d 9 -9

СУ ф й в ',

Р = 9 - 9 0.

 

 

 

J

V

 

 

 

 

%%

§ 2.13. Упругие среды с электромагнитными свойствами

263

Здесь фо ёо —вщначальное значение свободной энергии; cv{6 ) — тепло­ емкость; 4С — тензор модулей упругости; рь — симметричный тензор второго ранга — тензор магнитной проницаемости; к — симметричный тензор второго ранга — тензор диэлектрической проницаемости; и — тензор второго ранга, называемый тензором магнитоэлектричества; 3М е и 3М/г — тензоры третьего ранга, симметричные по вторым и третьим индексам (в силу симметрии тензора е): dM e = dMe , dM/j = , называемые тензором пьезоэлектрических коэффициентов и тензором пьезомагнитных коэффи­ циентов соответственно; /3 = • • а (где а — тензор теплового расширения); qh — вектор пиромагнитного эффекта;

а = • • в - е • 3М е - h 3M/i - /М,

 

в

РП = РПв + /3 • • £ + q/г • h + qe • е -

а • • /3$, ‘Пв = Щ + (cv/e r) de',

b = ц ■h + 3M/j • • г + v ■e + q h$,

d = 3Me ■■e + K - e + u T-h + qeP.

 

(2.13.32)

Эти соотношения совместно с (2.13.17) и (2.13.18) образуют полную систему определяющих соотношений для упругих электромагнитных сред.

Е. Уравнения термомеханики

Система уравнений термомеханики, состоящая из уравнений термоупруго­ сти и теплопроводности, для рассматриваемого типа электромагнитных сред

имеет вид

 

_

. ° р

 

° d2u

(2.13.33)

р

-^ 2

= v - c r

+ p f + p" iFem

° р 0 ^

= V - ( \ - V e ) + °pqm + e - j ,

(2.13.34)

 

£ = ^(V<X>u + V ® u T).

(2.13.35)

Эту систему дополняют выражением для пондеромоторной силы (т. 2, (2.10.16а)) с учетом малости деформаций:

pfem = рее + 1 j х b + (V 0

е) • р + (V 0

Ь) • m + 1 J^(p х Ь).

(2.13.36)

Выражение для производной от р, согласно (2.13.32), имеет вид

 

1 =

+ * - | -

- * * •

<2 1 3 -3 7 >

Подставляя (2.13.37) в уравнение баланса энтропии (2.13.34), приходим к

следующему уравнению теплопроводности:

 

 

 

о дв

_

^

де

дЪ.

де о_

.

опч

рСе ^

= V- (A - V 6 » )-/3 "

— - q h - —

- q e - — + pqm + e -),

(2.13.38)

где qe вектор пироэлектрического эффекта.

264

Глава 2. Упругие среды с малыми деформациями

Подставляя выражение (2.13.31) в (2.13.30), получаем определяющие со­ отношения для модели линейно-упругих электромагнитных сред.

Ж. Замкнутая система уравнений термомеханики и электродинамики

Система семи скалярных уравнений электродинамики (уравнения (2.13.17г) и (2.13.17д), как всегда, исключаются из общего числа), трех скалярных уравнений движения (2.13.33) и скалярного уравнения тепло­ проводности (2.13.38), в которые подставлены определяющие соотношения (2.13.18), (2.13.32), соотношения Коши (2.13.35) и выражение (2.13.36) для пондеромоторной силы, является замкнутой системой относительно семи скалярных электромагнитных функций (компонент векторов е и h, а также плотности заряда ре), трех компонент вектора перемещений и и температу­ ры в:

е, h, ре, и, в || х \ t.

(2.13.39)

Граничные и начальные условия к этой системе имеют вид (2.13.7), (2.13.8), (2.13.15); к ним также добавляют условие для электрического заряда при t = 0: ре = ре0.

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

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

уравнений:

 

 

 

 

 

 

 

 

 

дре

+ V - j

= 0,

1

dd

= V

х

1

дЪ

V х е,

dt

 

 

с

dt

 

 

с

dt

 

 

 

 

V

• b = 0,

V • d = ре,

 

 

(2.13.40)

 

 

j = R -1 • е,

b = fi • h, d = х • е,

 

состоящую из семи скалярных уравнений относительно семи скалярных неиз­

вестных:

(2.13.41)

е, h, ре || хг, t.

2.13.6. Упругие диэлектрики с поляризацией и намагниченностью

Для диэлектриков электрические заряды и электрические токи отсутству­ ют, поэтому в системе уравнений (2.13.17), (2.13.33), (2.13.38), (2.13.18), (2.13.32), (2.13.35) (2.13.36) следует принять

Ре = о, j = 0.

(2.13.42)

§ 2.13. Упругие среды с электромагнитными свойствами

265

В результате число неизвестных в (2.13.39) и число уравнений сократится на единицу.

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

1)инерционным членом в уравнении движения (2.13.33);

2)скоростями изменения электромагнитных величин в уравнениях Макс­ велла (2.13.40);

3)вкладом электромагнитных и механических величин в энтропию;

4)квадратами электромагнитных величин в уравнении движения (2.13.33) после подстановки в него (2.13.36) и (2.13.21), поскольку полагают справедливыми следующие неравенства:

|| р

||< ||

V • а

||,

|| е ® р

||< || а ||,

|| b ® m ||< ||

а ||,

м 1 д ,

ч и

и

и

и 1

и и

и

и 1 ab и и

и

H - S ( p x b ) | |« | | V . < T | |,

| | - ^ | | « | | V x h | | .

| | - _ | | « | | V x e | | ,

 

 

 

 

II v - т IKII

т | | .

 

(2.13.42)

Всилу отсутствия электрических токов и электрических зарядов (2.13.42),

атакже в силу допущения 2, уравнения Максвелла (2.13.17) сводятся к

следующей системе:

V х h = 0,

V

х е = 0,

(2.13.43)

V • d = 0,

V - b = 0.

(2.13.44)

Первые два уравнения этой системы допускают простые интегралы:

h = -V<£, e = - W ,

(2.13.45)

где Ф и U — скалярные поля, называемые магнитным и электрическим потенциалами соответственно. Очевидно, что выполняются соотношения

V x h = - V x V $ E 0 , V х е = - V х V U = 0.

(2.13.46)

В результате остаются только два скалярных уравнения электродинамики (2.13.44) относительно двух скалярных функций U и Ф. Эти два уравнения решают совместно с системой уравнений термомеханики (2.13.32), (2.13.33),

(2.13.35), (2.13.38), которая с учетом допущений

1, 3 и 4 принимает следую­

щий вид:

 

 

_

 

 

 

 

V • а + pi = 0,

 

 

0

дв

_

X , _ х,Х о

_

 

Рсе

—V • (Л • V0) + р qm,

 

V • d = 0, V - b

= 0,

сг

= 4С • • в —е • 3М е —h • 3M/j —/3 г?,

b = ц, ■h + 3М/( • • е + и ■е + q/г

(2.13.47)

266 Глава 2. Упругие среды с малыми деформациями

d = к • е + 3М е • • е + и т • h + qe $,

e = i( V ( g ) u + V (g > u T), Ь = ^ Ф , e = - W .

Граничные условия к системе (2.13.47) для электромагнитных величин обычно записывают в виде заданных значений потенциалов Ue и Фе или потоков электрической и магнитной индукции dne, Ъпе\

и \Т:и = ие, Ф|Еф = Фе или п • d |E^ = dne, n • b |Efc = bne. (2.13.48)

Механические и тепловые граничные условия имеют вид (2.13.15):

n - a \ ^ = t n e , u|Eu =

Ue,

= Qe, -n -A -V(9| E? = Qe,

(2.13.49)

где E = Eu U E,, = ECTU E? = Et/ U Ed = ЕФU S b.

 

Присоединяя к системе

(2.13.47)—(2.13.49) единственное начальное усло­

вие

£ = 0:

в = в0,

(2.13.50)

 

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

и, 0, U, Ф ||

х \ t.

(2.13.51)

Эта задача является линейной относительно неизвестных функций.

Более того, задача теплопроводности

(для

нахождения температуры в)

не связана с задачей электроупругости и может быть решена независимо. Уравнения электродинамики и теории упругости в системе (2.13.47) за счет наличия тензоров пьезоэлектрических и пьезомагнитных коэффициентов 3М е и 3М/г остаются связанными и решаются совместно.

Если же рассматривают сплошные среды, являющиеся диэлектриками, но не обладающие пьезоэлектрическими и пьезомагнитными эффектами, т. е. для которых 3М е = 0 и 3М/г = 0 (например, для изотропных сред эти тензоры всегда нулевые, см. [18]), то задачи электродинамики и теории упругости в системе (2.13.47)—(2.13.49) также становятся несвязанными.

В частном случае, когда отсутствуют еще и пироэлектрические, пиро-

магнитные и магнитоэлектрические эффекты, т. е. и

= 0,

= 0, qe = О,

из (2.13.47) и (2.13.48) получаем два независимых

скалярных уравнения

эллиптического типа с соответствующими граничными условиями для U и Ф:

V - ( * - V U )

= 0,

u\llu = Ue,

- n - x - W | Sd = dne,

(2.13.52)

У - ( ^ ^ Ф )

= 0,

Ф|Еф = Фе,

- п - / х ^ Ф | ъ ь = Ъпе,

(2.13.53)

которые образуют постановку задачи электростатики (2.13.52) и магнито­ статики (2.13.53).

§ 2.13. Упругие среды с электромагнитными свойствами

267

Несложно заметить, что постановки этих задач совпадают формально с постановкой стационарной задачи теплопроводности, если температурное поле в(хг) в теле не зависит от t:

V- ( A- V0) = o, 0|Efl=0e, - n - A - V 0 | Sq = qe-

(2.13.54)

Методы решения всех трех задач (2.13.52)—(2.13.54) одинаковые.

В частности, если тело изотропное, то тензоры х, ц и Л шаровые: х = хЕ,

ц

= рьЕ, А = АЕ, где х — коэффициент диэлектрической проницаемости;

р

— коэффициент магнитной проницаемости; А — коэффициент тепло­

проводности.

 

 

 

 

Если х, р и А не зависят от х \ то задачи (2.13.52)—(2.13.54)

становятся

краевыми задачами для уравнения Лапласа:

 

 

 

AU = 0 в V,

и к , = ие,

- х п • VC/|Sd = dne,

 

АФ = 0 в У,

ф||£Ф= Ф е ,

-ЦП

= Ъпе,

(2.13.55)

Ав = 0 в V,

в \^е = ое,

A n-V #L

= qe.

 

Г л а в а 3

У П Р У Г И Е С Р Е Д Ы С К О Н Е Ч Н Ы М И Д Е Ф О Р М А Ц И Я М И

§ 3 .1 . З а м к н у т ы е с и с т е м ы у р а в н е н и й д л я у п р у г и х ср е д в п р о с т р а н с т в е н н о м о п и с а н и и

3.1.1. вR U V F -система уравнений термоупругости

Перейдем к формулировке задач для общего случая упругих, т. е. идеаль­ ных твердых сред с произвольными конечными деформациями, определяющие соотношения для которых были установлены в т. 2, §§ 3.5-3.9. Для твердых сред с конечными деформациями эта формулировка включает в себя две группы уравнений:

законы сохранения (т. 2, (2.12.5));

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

Общая система законов сохранения (т. 2, (2.12.5)) в пространственном

(эйлеровом) описании состоит из уравнения неразрывности ( а =

1), уравне­

ния движения ( а =

2), уравнения энергии ( а

= 3), динамического уравнения

совместности деформаций ( а

= 5) и кинематического уравнения ( а = 6):

 

 

 

!

+ V . p v = 0,

 

(3.1.1а)

 

^ot

+ V

pv <g>V = V

• Т + p i,

(3.1.16)

^ + V

(pve +

q) =

V • (T v) + pi v +

pqm,

(ЗЛЛв)

 

afdt T + V - p ( v ® F T

Fcx)v) =

0,

(3.1.1г)

 

+ V (pv ® u) = pv.

 

(3.1.1д)

Вообще система законов сохранения (т. 2, (2.12.5)), как отмечалось, со­ стоит из шести групп уравнений ( а = 1, ..., 6), суммарно образующих 18 скалярных уравнений: 1 + 3 + 1 + 1+ 9 + 3 = 18. Однако после добавления к этой системе определяющих соотношений необходимо исключить из системы (т. 2, (2.12.5)) одно из уравнений термодинамики. Действительно, вспомним (см. т. 2, § 2.5), что часть определяющих соотношений эквивалентна ОТТ, которое было получено суммированием двух уравнений термодинамики (см.

§ 3.1. Замкнутые системы в пространственном описании

269

т. 2, § 3.3), поэтому ОТТ заменяет одно из них. Чаще всего исключают уравнение баланса энтропии (т. 2, (2.12.5) при а = 4). В этом случае система (3.1.1 а)—(3.1.1 д) состоит из 17 скалярных уравнений и содержит 27 скаляр­ ных неизвестных:

р, v, u, Т, е, в, q, F || х, t.

(3.1.1е)

Каждое векторное поле v, u, q эквивалентно трем скалярным функциям, тензорное поле Т — шести скалярным функциям, а поле F — девяти скалярным функциям.

Для замыкания этой системы к ней следует присоединить определяющие соотношения.

Часть этих соотношений имеет одинаковый вид для всех идеальных твер­ дых сред — это закон Фурье (т. 2, (3.12.7)) и выражения для плотностей

полной энергии е и внутренней энергии е (т. 2, (3.3.12) и (3.5.4)):

 

 

 

 

 

q = -Л

• V0,

 

 

 

(3.1.2а)

 

 

 

 

 

|

lv l2

 

/

I

у2

 

 

(3.1.26)

 

О с т а л ь н ы е

 

£ = е + ^ Г = ф - в 1ю + ^ -

 

д л я и д е а л ь н ы х

т в е р д ы

о п р е д е л я ю щ и е

 

с о о т н о ш е н и я

в и с я т

о т в ы б о р а

т о й

и л и

и н о й

м о д е л и

с р е д ы .

В и д э т и х

с о о т н

п о л у ч е н

в т. 2 ,

§ 3 .8 .

 

 

 

в т . 2 ,

п . 3 . 8 . 1 0

н а и б о л е е

о б щ и м

В о с п о л ь з у е м с я

у с т а н о в л е н н ы м

л е н и е м

1

о п р е д е л я ю щ и х

с о о т н о ш е н и й

и д е а л ь н ы х

т в е р д ы х с р е д в

( 3 . 8 . 1 1

) )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T =

£

G

( F ,

0 ) .

 

 

 

 

 

 

 

ф = ^ Ф 3 ( C G(F)), в) = ф(¥,

в),

 

(3.1.36)

 

( п )

где тензорные функции Т G и потенциал ф зависят от выбора модели Ап, Вп, Сп или Dn и определяются соотношениями (т. 2, (3.8.109)—(3.8.115))

T c ( F , e ) = 4 E G .. T G = 4 E g -- Т с ( ¥ , в ) = ^ 2 ^ 1 4 F G

4 3 = d l ^ / d C c ,

7=1.

^ = p ( d < p / d l % )

4 3 = hGp s)(C G) + (1 -

( n )

M 4 S)(o T• C G ■O),

( n )

C G = n - I II

(n)3

4E G = УЗ

^ 2 ( К ~ Ш(кОРа ® Pa

+ (! - M P a ® P a )) “

E,

a=\

 

 

О

о

(3.1.4)

Eaf3Pa®Pf3® (hGpp(g)pa + (l - h G)p p ® p a).

a ,/3=1

270

Глава 3. Упругие среды с конечными деформациями

 

Здесь Ц(s)] — инварианты тензоров к J

относительно группы симметрии Gs

рассматриваемой упругой среды;

 

 

 

Ра,

Ра || F

(3.1.5)

собственные значения и собственные векторы, являющиеся

функция­

ми

тензора F; коэффициенты Еар выражаются только через Аа

(см. т. 2,

(3.2.38а)); О = р^ 0 р'3 — тензор поворота, сопровождающий деформацию, также зависит только от F; ha и ha — функции-указатели класса модели:

1,

G = А, В,

-

1,

G = А, С,

(3.1.5а)

ha

G = C,D,

ha

0,

G = B ,D .

0,

 

 

Внутреннюю энергию е можно рассматривать как функцию F и в, т. е.

 

e = V

> - ^

= e(F,

в).

(3.1.6)

Таким образом, получаем систему из 10 скалярных определяющих со­ отношений (векторное соотношение (3.1.2а) эквивалентно трем скалярным, соотношение (3.1.26), в которое подставлено (3.1.36), — одному скалярному соотношению, а (3.1.3а) — шести скалярным соотношениям), которые вместе с (3.1.1) образуют замкнутую систему из 27 скалярных уравнений относи­ тельно 27 скалярных неизвестных функций (3.1.1е).

Число уравнений и неизвестных в этой системе может быть сокращено,

если подставить

определяющие соотношения (3.1.2) и (3.1.3) в (3.1.16) и

(3.1.1 в). В результате этого

получим

систему

из

17 скалярных уравнений

относительно 17

неизвестных:

 

 

 

 

 

0,

р, z/, v,

F || х,

t.

(3.1.7)

Такую систему уравнений (3.1.1) называют вRU VF-системой уравнений термоупругости.

(П)

(п)

г

, ч

Замечание 3.1.1. Выбирая функции T Q = ^ G ( C G ,0) =

 

для мо"

7—1

делей Ап, Вп, Сп или Dn в каком-либо виде из числа рассмотренных в т. 2, § 3.8, получим частные представления тензорных функций (3.1.4). Так, если рассматривают модели Ап, то общее представление функции таково (см. т. 2, (3.8.48)):

(п)

(п) (п)

(п)

(п)

(3.1.8)

С,

в) = 4М •• С + £ С 2 + 6L --- (С ®

С),

где 4М — квазилинейный тензор упругости, £ — параметр квадратичной упругости, 6L — тензор квадратичной упругости, вообще говоря, завися­

( п )

щие от инвариантов 1 ф (С )

Для квазилинейных моделей Ап полагают (см. т. 2, (3.8.54)): £= 0, 6L = 0,