Файл: Хачатурян, А. Г. Теория фазовых превращений и структура твердых растворов.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 16.10.2024

Просмотров: 188

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

§22. Энергия внутренних напряжений упруго анизотропного кристалла,

содержащего когерентные включения новой фазы [155, 156]

Как уже отмечалось в предыдущем параграфе, фазовый пере­ ход, протекающий в состоянии, свободном от внутренних напря­ жений, приводит к изменению объема и формы превратившегося кристалла. Это изменение обусловлено перестройкой кристалли­ ческой решетки и может быть выражено через тензор однородной

деформации е?,-. Деформацию е?;- удобно отсчитывать от состояния, в котором находится матрица в отсутствие внутренних напряже­

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

Пусть в бесконечной среде — матрице — находится ѵ типов включений, каждое из которых характеризуется своей деформа­

цией вij = Eij(p) (р = 1, 2, ..., ѵ), т. е. своей кристаллической ре­ шеткой, отличной от решетки матрицы. Примем, что постоянные упругости всех фаз, образующих гетерогенное состояние, одина­ ковы. В этом случае упругая энергия единицы объема может быть записана в форме

V

 

/ ( Г) = ~ 2 а°ІІ % (г) 8« + ~ Y hilm&ißlm,

(22.1)

_ _ _ _

P = 1

 

 

 

где \ j lm — тензор модулей упругости,

а?}- (р)

— постоянный тен­

зор,

характеризующий свойства фазы

р (его

физический

смысл

будет установлен ниже); Ѳр (г) — функция формы включений, рав­ ная единице, если радиус-вектор г попадает в какое-либо включе­ ние сорта р, и равная нулю в противоположном случае; индексы г, /, I, т описывают декартовы координаты. По дважды повторяю­ щимся индексам подразумевается суммирование. Деформация

&U, так же как и деформация (р), отсчитывается от состояния, отвечающего недеформированной (и ненапряженной) матрице. Выражение (22.1) представляет собой два первых неисчезающих члена разложения плотности свободной энергии по деформации еі;-. Коэффициенты этого разложения

V

2 °Ь (р) (г) и хШт р=і

являются константами материала. Разложение (22.1) ограничено квадратичными членами по деформации, так как учет членов бо­ лее высокого порядка вывел бы нас за рамки линейной теории'упругости. Присутствие в выражении (22.1) члена, линейного по де­ формации, отражает тот факт, что рассматриваемая система про­

200


странственно неоднородна и поэтому ее ненапряженное состояние не является недеформированным *). Иными словами, при Оц(г) = = О, где Oij(r) — тензор напряжения, егДг) ^ 0. Деформация Eij(r) в последнем случае оказывается равной одному из следую­

щих значений: е®,(1), Sy (2), ..., Sij(p), ..., e?j(v) в зависимости от того, во включение какого типа попадает вектор г. Если система является пространственно однородной и ненапряженное состоя­ ние является одновременно недеформированным состоянием, то

линейные члены по

в выражении для плотности упругой энер­

гии отсутствуют.

 

 

 

 

Выражение для тензора напряжений ац(т) можно получить

по обычным правилам,

дифференцируя по

плотность упругой

энергии (22.1):

 

 

 

 

(г) —

2

(р) (Г) "t*

(22.2)

 

р=1

 

 

 

Пусть радиус-вектор г находится внутри включения p-то типа, находящегося в ненапряженном состоянии. Тогда, по определе­ нию, имеем следующие соотношения:

Ой (г) = 0, 'Ѳр(г) = 1, eitn = e?m(jtj).

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

(?) — nfilm[p)i (22.3)

которое, по существу, представляет собой определение постоян­ ных материала о%(р), фигурирующих в выражении (22.1), через

характеристики кристаллогеометрии фазового превращения ец (р). Воспользовавшись выражением (22.2), можно представить уравне-

ние равновесия dOijldr) =

0

в

виде

 

^ ^ Г = І

<

( Р

) ^ Ѳ р ( г ) -

(22.3а)

І

 

р = 1

j

 

Выражая деформацию еІт (г) через вектор смещения и (г):

, V

1 ( ди1(г) ,

дит(г) \

(22.4)

S‘m(г) -

2 \ дгт

+

drt )

 

подставим (22.4) в (22.3а). При этом получим:

u 5 r = i = ! i ( 4 e P(r).

(22.5)

При получении выражения (22.5) мы воспользовались симмет­ рией kijlm‘ = hijml-

х) Ситуация здесь аналогична той, которая имеет место в теории термо­ упругих напряжений (см., например, [161]).

201


Уравнение (22.5) удобно решать с помощью преобразований Фурье:

00

00

 

f (к) = gjjd3r /( r ) e-*r,

/(г) = й ( ^ / ( к ) е ікг,

(22.6)

где f( к) есть фурье-образ функции / (г), являющейся фурье-ори- гиналом, и к — волновой вектор — параметр преобразования Фурье. Умножая правую и левую части уравнения (22.5) на ехр(—ikr) и интегрируя по бесконечному пространству, получим:

V

{hiimbjkt) йт(к) == — і 2

оу (р) bßp (к).

(22.7)

р=1

 

и (к) = ^dsr exp (— ікг) и (г), Ѳр (к) =

^d3r exp (— ikr) Ѳр (г).

(22.8)

В отличие от исходного уравнения (22.5), уравнение для фурьѳобраза и (к) является чисто алгебраическим и может быть легко решено. Уравнение (22.7) можно переписать в более компактной операторной форме:

V

 

б-* (k) I и(к)> = - і S 3° (р) I к> Ѳр(к),

(22.9)

р = і

 

где G-1 (к) и ои (р) — операторы, матричные элементы которых в

тензорном представлении есть "кілт kjkl и с\} (р ) соответственно. Решение уравнения (22.9) можно получить, умножая (22.9) на

оператор G (к), обратный оператору G-1 (к )*):

V

и (к) = - і& (к) 2 о0 (р) I к> Ѳр (к).

(22.10)

р=і

 

По определению, матричные элементы оператора G (к) явля ются фурье-образами тензора Грина анизотропной задачи теории упругости. Для того чтобы более детально определить структуру

оператора G (к), рассмотрим выражение для матричных элемен­ тов (компонент тензора) обратного оператора G_1(k):

[G-1(k)ly =

Ximnjkmkn = к2 [fi-1 (n)]„

(22.11)

или

 

(22.12)

Ü - Ң к ) = 1 с*& -Ң п ),

где

 

 

n = kJk,

[Q-1 (n)]y = Ximnjnmnn.

(22.13)

x) ö (k) G-1 (k) = 1, гдѳ 1 — единичный оператор, или в матричных обозначениях [б (к)]н [Gr1 (k)];m = 6jm.

202


Йз определения G(k)G *(к) = G (к)Лг2й -1 (п) =

І следует, что

G(k) = -^Q (n),

(22.14)

где 12 (п) — оператор, обратный к й -1(п), т. е.

12 (п) й -1 (п) = і ,

или в тензорных обозначениях й ;г(и)[й-1 (п)]і7- = бі;- Из выраже­

ния (22.14) следует, что оператор 6 (к) является однородной функцией |к|.

Воспользуемся теперь решением (22.10) уравнений теории уп­ ругости для вычисления полной упругой энергии среды. Полная упругая энергия может быть получена из выражения (22.1) для плотности упругой энергии в результате интегрирования ее по всему бесконечному пространству:

О О

V

 

 

и упр = $

d3rj~— 2 4 ІР) <5р (г) ец (г) + 4 -

(22.15)

В соответствии с определениями (22.6) фурье-оригиналов через

фурье-образы заменим функции Ѳр (г) и е;Дг) в (22.15) через их фурье-образы:

^ r ) = ! i w 0 p ( k ) e i k r ’ <2 2 Л 6 )

«*ц w =

(k) eikr = $W

‘ І г {кіЩ(k) + kiSl (k)] eikr- (22-17)

Подставляя

(22.16) и (22.17) в

(22.15) и учитывая соотношение

 

J d3r

(2я)3 б (к + к'),

где б (к + к') — дельта-функция Дирака, получим:

V

упр = $ tSf S 4 (р) % (к)*і«Г(к) + 4- (W A ) Sm(к)Щ (к)} .

'' L р=х

(22.18)

Подинтегральное выражение в (22.18) удобно переписать в опера­ торной форме:

и упр =

V

<«* (к)I (Р)I к>Ѳр(к) + і- (S'* (к)I G-і (к)IS (к))}.

2

" ' '

^ pel

 

(22.19)

Подставляя (22.10) в (22.19), получим:

 

^ упр = ^ш

{ 2

-

<к 11°0 ) 6 (к) ° 0 (?) I к> % (к) в; (к) +

 

 

+

~

<к IS0(р) G (к) G-1 (к) G (к) о° (?) I к>},

(22.20)

203