CC BY
24
А. С. Михайлов
ОБ ОПРЕДЕЛЕНИИ КОЭФФИЦИЕНТА ПОВЕРХНОСТНОГО НАТЯЖЕНИЯ В ДВУХФАЗОВЫХ ЗАДАЧАХ ТЕОРИИ УПРУГОСТИ ПРИ УСЛОВИИ НЕСЖИМАЕМОСТИ ИЛИ НАЛИЧИИ ЖЕСТКИХ ВКЛЮЧЕНИЙ*
1. Введение. Рассматривается модель фазовых переходов в механике сплошных сред [1], в которой функционал энергии состоит из суммы энергии деформации и поверхностной энергии. Существует несколько способов учета поверхностной энергии [2, 3], в том числе прямой способ, при котором поверхностная энергия считается пропорциональной площади границы раздела фаз с коэффициентом пропорциональности, называемым коэффициентом поверхностного натяжения. Эта физическая характеристика зависит от рассматриваемой упругой среды, мала, но ее точная величина неизвестна. Возникает вопрос, можно ли оценить коэффициент поверхностного натяжения через другие характеристики упругой среды? Подобное исследование проводилось в работе [4].
Мы рассмотрим однородную анизотропную двухфазовую упругую среду с нулевым граничным полем смещений и дополнительным условием несжимаемости или абсолютной жесткости одной из фаз. Следуя схеме [4-7], мы исследуем соответствующий функционал энергии с нулевым и отличным от нуля коэффициентом поверхностного натяжения, рассмотрим вопрос зависимости температур фазовых переходов от размера области и, в заключение, ответим на основной вопрос — дадим метод определения коэффициента поверхностного натяжения.
2. Постановка задачи. Пусть О С К" — ограниченная область с липшицевой границей, в которой расположена двухфазовая упругая среда. Пусть функция и €
К") задает поле смещения каждой точки х € О упругой среды, а характеристические функции х = (х+ ,Х-), Х± € БУ (О) своими носителями задают распределения фаз упругой среды по области О.
Функционал энергии двухфазовой упругой среды имеет вид
I[и,х^,а] = !(Р+ (Уи,г)х+(х)+ Р- (Уи,г)х- (х)) dx + аJ\Dx\- (1)
п п
Здесь Р± — плотности энергии деформации каждой из фаз, параметр £ играет роль температуры, а = (а+ ,а-) — коэффициенты поверхностного натяжения и
аJ \Dх\ = а+ У \Dх+\ + а- ! ^х-\,
п п п
где
/ \Df \ = вир / /divhdx.
■! кеС1(П,Жп),\к\<1 ■!
п 0 п
* Работа выполнена при финансовой поддержке РФФИ (грант №05-01-01063) и Совета по грантам Президента РФ для государственной поддержки молодых российских ученых и ведущих научных школ (грант №НШ-8336.2006.1).
© А.С.Михайлов, 2006
Будем предполагать, что фазы не перемешиваются, т. е. выполняется условие
Х± —характеристические функции, х+ (х) + X-(х) = 1 при п.в. х € О. (2)
Определим плотности энергии деформации ¥± (М,А) равенствами
¥ + (м,а) = а+ м (& - с+ )(Ы - Ф - а,
¥-(м,а) = а-м- С- )(Ы - С-)+А, (3)
£=±(М + М*), *ем,
где коэффициенты а± ы € К удовлетворяют условиям симметрии
а±,к1 = а±1Ц = а±г,Ы = а±,1к
и условию положительной определенности
а1],к&3 Сы ^ Сг] , V > 0,
справедливому для всех симметричных матриц С, а симметричные матрицы в (3)
играют роль тензоров остаточного напряжения.
В дальнейшем будем рассматривать две задачи на минимум функционала (1) с плотностями энергии деформации (3) при различных дополнительных ограничениях на поле смещений.
Первое ограничение говорит о несжимаемости фазы с индексом «+». Итак, в качестве первого множества определения функционала (1) рассмотрим
А = {(и,х): Х+ и = 0}. (4)
Второе ограничение говорит о том, что фаза с индексом «+» состоит из абсолютно жесткого вещества, т. е. не испытывает деформаций:
/ иг —I— и \
02 = {{и,х) : Х+( 2 Хг) =0, 7/'../ (5)
Такие вариационные задачи рассматривались в монографии [2] и статье [8], где были доказаны теоремы существования при положительных коэффициентах поверхностного натяжения и получены необходимые условия минимума. Отличие рассматриваемой ситуации от изученной состоит в том, что плотности энергии деформации ¥± не зависят от точки х € О. Это позволяет исследовать задачи при нулевом коэффициенте поверхностного натяжения, а также пронаблюдать зависимость характера состояния равновесия от параметров задачи, чисел А и а.
3. Задача с нулевым коэффициентом поверхностного натяжения. При а =
0 функционал (1) может быть переписан следующим образом
I [u,х,t, °\ = ! а+м Сю Ы Х+ Лх - 2а+м С+1 ! йэ Х+ Лх +[а%м С+ С+ - А J Х+ Лх+ п п п
+ ! а-,ы Сгэ Сы х- Лх - 2а-э м С- ! Сэ х- Лх +[а-,к1 С- С- + А ! X- dx,
п
пгх , +u3xг
где = —3——Привлекая условие (2), перепишем последнее выражение в виде
I [u,x,t, 0] = J aij, ki (x)&j Cki dx - 2Aij JCij x+ dx + A(t), (6)
где
ац,ы(х) = %,ы Х+(х) + %,ы X (х)
Аг3 = а+),Ы С+ — ач,Ы скI,
А(Ь) = I [0^,0].
Заметим, что если минимизировать функционал (1) на множестве (5), то второе слагаемое в правой части (6) обнулится и мы придем к выводу, что
I[и, х,Ь, 0] = J Огз,кI (х)&зСкI Лх + I[0 , X , Ь, 0] ^ I[0 , X , Ь, 0].
0
Таким образом, минимум функционала (1) на множестве (5) достигается на паре (и, X), где и = 0. Подставим и = 0 в (1) при а = 0:
[0,х,^,0] = [а+з}ыС+ С + — Ь] Iх+ Лх + [ау,ыСг-Сы + ь]1(1 — х+) Лх.
о о
Правая часть последнего равенства линейно зависит от § х+ Лх, следовательно, функ-
о
ционал достигает минимума в крайних точках х+ = 1 или х+ = 0, или если множители при / х+ Лх и § (1 — х+) Лх уравновешены, при произвольной характеристической функ-
оо
ции х+. Результат вычислений сформулируем в виде теоремы.
Теорема 1. Пусть минимум функционала (1) при а = 0 на множестве (5) реализуется на паре (и, X)- Тогда и = 0 и X+ = 1, если Ь > Ь*; X+ = 0, если Ь < Ь*; X+ — произвольная характеристическая функция, если Ь = Ь*, где
1 = 2- ац,к£ц^ы)- (^)
Таким образом, состояния равновесия функционала (1) при а = 0 однофазовы всегда, за исключением одного выделенного значения температуры Ь*. Если пронаблюдать динамику изменения характера состояний равновесия при повышении температуры Ь, то можно увидеть, что среда из однофазового состояния « —» перерождается в однофазовое состояние «+» при критическом значении температуры Ь = Ь*.
Изучим вопрос о характере состояний равновесия функционала (1) при а = 0 на множестве (4). Для этого запишем (6) в виде
I[u, x, t* + е, 0] = / aij,ki (x)€ijСki dx - 2Aij f Cijx+ dx + I[0, x, t* + e, 0] =
Q Q
= / aij,ki (x)Cij Cki dx - 2 Aij J Cij x+ dx + e J (x- - x+) dx + I [0,x,t*, 0].
Проанализируем полученное соотношение. Первое слагаемое правой части (9) квадратично зависит от Ум, второе — линейно, а последние два не зависят от Ум.
Пусть Агу = с5гу, где 5гу = 0 при г = у, 5гу = 1 при г = ]. Тогда второе слагаемое в
правой части (9) перепишется в виде
2А{^ х+ Лх = 2с5{у ! х+ Лх = 2ёгу их+ Лх = 0.
п п п
Последнее равенство выполнено, т. к. пара (и, х) принадлежит множеству определения (4). Следовательно, в (9) отсутствует линейное по Чи слагаемое и минимум достигается на паре (й,Х) при и = 0. Таким образом, (и,Х) € ^2, и можно утверждать, что мы находимся в условиях теоремы 1.
Пусть Агу = с5гу, тогда можно подобрать функции (и, Х) € так, что ё1у и = 0,
Г ~ ~ г~ их. + и{.
04э,ы(х)&зЫ Лх - 2Ац ^х (1х = -Л = —Л(м, х) < 0, ^ = —3—------------------1- (Ю)
пп С учетом последнего тождества, (9) предстает в виде
1 [и,Х,Ъ* + £: 0] = —Л + £J(Х — Х+) ЛХ + 1 [0,Х,Ъ*, 0] =
я
= —Л + £J(Х — Х+) Лх + 1 [0,Х,Ъ* , 0] =
я
= —Л + £J(Х — Х+ ) Лх +1 [0,Х,Ъ* + £1 0] — £ !(Х — Х+) Лх =
я п
= —Л + £J(Х — Х+) Лх +1 [0,Х,Ъ* + £1 0] + ИМ. (11)
п
В последней цепочке Х+ = 1 при £ > 0, Х- = 1 при £ < 0, а при £ = 0 Х+ - произвольная характеристическая функция. Правую часть (11) можно оценить следующим образом:
1 [ит Х,1* + £т 0] ^ —Л + 2ИМ + 1 [0, Хт Ъ* + £т 0]
откуда видно, что при малых £ состояния равновесия функционала I[и, х, Ъ* + £, 0] двухфазовые. А именно, верна следующая теорема.
Теорема 2. Пусть минимум функционала (1) при а = 0, Ъ = Ъ* + £ на множестве
(4) реализуется на паре (и, Х) и пусть а+У ы£+ — а-у ыС— = с5гу, где 5гу = 0 при г = ],
5гу = 1 при г = ], с € К1. Тогда и = 0 и 0 ф Х+ Ф 1, если
N < ± (12)
где Ъ* задается тождеством (8), а число Л определяется в (10).
Если а+у ыС+1 — а-у ы^ы = с5гу, то верно утверждение теоремы 1.
Сделаем замечание к последней теореме. Число Л в (10) определяется неоднозначно. Оно зависит от функций (и, Х), и чем лучше мы сможем ими распорядиться, тем лучше оценим Л. В статье [7] указан метод построения оптимальной оценки Л, к сожалению, в нашем случае этот метод не работает из-за наличия дополнительного ограничения на
поле смещений (условия несжимаемости). Метод состоит в явном построении функций (й, х) таких, что Л в (10) становится большим. Для этого выписывается вспомогательная задача в шаре и строятся ее симметричные решения. Специфика нашей ситуации заключается в том, что все симметричные решения из-за наличия условия несжимаемости (4) автоматически становятся нулевыми и, как следствие, Л = 0. Построение несимметричных решений в явном виде — это отдельная задача, выходящая за рамки данной статьи. Поэтому мы ограничимся только лишь существованием таких решений.
Теорема 2 характеризует состояния равновесия функционала (1) при а = 0 на множестве (4) при температурах, мало отличающихся от критической Ь = Ь*. Возникает вопрос — что будет с состояниями равновесия при больших температурах? Ответ на аналогичный вопрос можно найти в работе [7], приведем его точную формулировку и упрощенное доказательство.
Теорема 3. Пусть минимум функционала (1) при а = 0, Ь = Ь* + е на множестве (4) реализуется на паре (й,х). Тогда й = 0 и х+ = 1 при е > е+, х— = 1 при е < е—, где
е+ 2 (.ац,к1^ц^к1 ~ ‘^‘а'1],к1^1]^к1 (а )kl,staij,kl^ijapq,stXpq)^ (13)
—2 (ац,к&Ц^к1 ~ ^аЧ],к1^г]^к1 (а ')kl,sta,ij,kl^ija,pq,stCpq)^
число Ь* задается тождеством (8), а (а±)—/^ — тензоры, обратные к а± ы.
Доказательство. Возьмем тождество (11), подставив в него Л из (10), считая, что (й, х) —пара, на которой реализуется минимум функционала (1) на множестве (4), Х+ = 1 при е > 0, х— = 1 при е < 0, а при е = 0 х+ — произвольная характеристическая функция:
1 [й, х,Ь + е, °] У агз:Ы (х)£г0 Ск1 Лх 2Аго [ х Лх+
+ е ! (х - х+) ЛХ + I[0,х,Ь* + е, 0] + |е||П| = J аго,к1 (х)£о£к1 Лх-
П П
- 2Аг^ х+ Лх + J ((|е| + е)х + (|е|-е)х+) Лх + 1 [0,х,Ь* + е, 0]- (14)
ПП
При е < 0 правую часть (14) можно упростить и, используя условие положительной определенности коэффициентов а—о ы, получить неравенство
1 [й,х,Ь* + е 0] — 1 [0,х,Ь* + е 0] > ( (а+,к1 (х)&3 £к1 - 2Аг] + (-2е))Х+ Лх =
П
= 2(е — — е) J х+ Лх + J(а+,к1 (х)6го £к1 - 2АЦбу + (—2е — ) )х+ Лх.
ПП
Исходя из определения числа е—, подынтегральное выражение последнего интеграла содержит полный квадрат, поэтому мы можем сказать, что
+
1 +е,0] -1 [0, х, Ь + е,0] > 2<е- - ;)/ х+ <ь
Следовательно, при е— — е > 0 х+ = 0, иначе пара (и, х) не доставляла бы минимум функционалу (1).
При е > 0 доказательство проводится аналогично, с заменой индекса «+» на « —», ввиду того, что
Проанализируем последние две теоремы. Нетрудно увидеть, что если выполняется условие а+ к1 С+1 — а— ы = с5^, то при Ь = Ь* у функционала (1) на множестве
(4) наблюдаются только однофазовые состояния равновесия. Если же условие не выполняется, то из теоремы 2 следует, что при температурах, близких к критической, у функционала (1) на множестве (4) наблюдаются многофазовые состояния равновесия, а из теоремы 3 следует, что при температурах, далеких от критической, все состояния равновесия однофазовые.
4. Задача с положительным коэффициентом поверхностного натяжения.
Из условия несмешиваемости (2) можно сделать вывод, что
Поэтому функционал (1) зависит от \а\ = а+ + а—. Исследуем эту зависимость. Подобная задача рассматривалась в [2, 5, 8-10], и, как отмечалось, верен следующий простой факт. Так как исходный функционал (1) линейно зависит от параметров Ь, \а\, если мы минимизируем его по (и,х) € , 2, то такой минимум будет являться вогнутой
функцией по переменным (Ь, \а\). Вогнутость функции означает выпуклость ее подгра-фика. Как простейшее следствие, мы получаем, что множества значений параметров ^, \а\) € М+, соответствующих каждому из однофазовых состояний равновесия, являются выпуклыми.
Для задачи с жесткими включениями из сказанного и теоремы 1 можно сделать следующий вывод.
Теорема 4. Пусть минимум функционала (1) при \а\ > 0 на множестве (5) реализуется на паре (и, х). Тогда и = 0 и х+ = 1, если Ь > Ь*; х+ = 0, если Ь < Ь*; х+ — либо 0, либо 1, если Ь = Ь*, где Ь* определяется равенством (8).
Аналогичное утверждение верно, если мы рассматриваем задачу на минимизацию функционала (1) на множестве (4) при условии, что а+ к1 (+ — а— к1 С,— = сдц.
Нам осталось ответить на вопрос, что происходит с состояниями равновесия функционала (1) на множестве (4), если последнее условие не выполняется. Из теоремы 2, теоремы 3 и свойства выпуклости однофазовых зон множества М+ Э (Ь, \а\) следует, что существует функция ао = ао(Ь), выпуклая на каждом из интервалов (—ж,Ь*), (Ь*, ж) такая, что график функции ао и (может быть) прямая Ь = Ь* делят М+ на три
□
п
п
п
п
зоны.
М1 = {(Ь, \а\) : Ь>Ь*, \а\ > ао(Ь)},
М2 = {(Ь, \а\) : Ь <Ь*, \а\ > ао(Ь)},
Мз = {(Ь, \а\) : Ь € М1, \а\ < ао(Ь)}.
(15)
В первой зоне, множестве М1 , все состояния равновесия однофазовые с чистой фазой «+». Во второй зоне, множестве М2, все состояния равновесия однофазовые с чистой фазой « —». В третей зоне, множестве Мз, все состояния равновесия двухфазовые. Открытым остается вопрос, как правая и левая ветки выпуклой на каждом из интервалов (—ж, Ь*), (Ь*, ж) функции ао(Ь) ведут себя при приближении к прямой Ь = Ь*. Ответ на него дают две следующие леммы.
Лемма 1. Функция ао(Ь), отделяющая зоны однофазовых состояний равновесия от двухфазовых на множестве параметров М2 Э (Ь, \а\) функционала (1), ограничена, и верна оценка
а0(Г)<^)(^7^)(^)1/П, (16)
где задаются равенством (7), а к(П) — константа в изопериметрическом неравенстве (21).
Доказательство. Метод доказательства аналогичен [5], здесь приведен из-за своей простоты и для полноты картины. Пусть пара (и, х) доставляет минимум функционалу (1) на множестве (4) при Ь = Ь*, \а\ > 0 и пусть состояние равновесия двухфазово, т. е.
0 ф х+ Ф 1. Преобразуя функционал (1) как в (9), получим
I [й,х,Ь*,а] = ! а^уМ(х)£гу £ы йх — 2ЛЦ ^ х+ йх + I [0,х,Ь*, 0] + ^\Бх\ = я я п
= J агз,ы(х)&з£ы йх — 2Л^ I ^х+ йх + а J\Пх\ + I[0,х,Ь*,а], (17)
п
где х+ ф 1. Перенесем последнее слагаемое из правой части (17) налево и учтем, что пара (й, х) доставляет минимум функционалу (1). В результате придем к неравенству
У а^,ы(х)£гз£ы йх + \Бх\ < 2Л^^ х+ йх = —2ЛЦ^ ^х— йх. (18)
п п п п
Пусть теперь хТ — та из функций х+ или х—, для которой
1 [ 1
(19)
п
Тогда из неравенства (18) следует соотношение
У а^м (х)£гз£ы йх + \а\ !\Бхт \ < 2\Л^ ^ \£^\хТ йх < п п п
<!'! 4'4'Г Лх + \М\М IГ (20)
пп Далее, пользуясь изопериметрическим неравенством
( /х^) " < к(П) ( |-Ох|, ^х : щ I Х&х < (21)
\П\
пп
и неравенством Гёльдера для функции хт
//» П — 1 л
Хт г1х < (j хт <1х^ " (у хт " ,
п п п
получим из (20) следующее:
|а| I |£Г| <; \ММк{П) I |^Г|( I Хт Лх) " •
п п п
Последний множитель можно оценить сверху, используя (19), после чего остается заметить, что в левой части неравенства число |а| не может быть большим, иначе неравенство перестанет быть верным. Таким образом, доказательство леммы завершено.
□
Лемма 2. Функция ао(Ь), отделяющая зоны однофазовых состояний равновесия от двухфазовых на множестве параметров М2 Э (Ь, \а\) функционала (1), при а+ ыСы — а— к1 С— = cSij не равна нулю и верна оценка
МП > (22)
где числа Л, Р заедаются ра,венства,ми (10), (24).
Доказательство. Перепишем (17) в виде
I[й, х, Ь*, а] — I[0, х, Ь*, а] = / а-,ы (х)£-йх — 2Л- / £-х+ йх + а $ \Вх\, (23)
п п п
где (и,х) —произвольный элемент из множества определения (4) функционала (1), хы ф 1. Положим в (23) \а\ = ао(Ь*). При таком критическом значении параметра а пара (0, х) является состоянием равновесия функционала, следовательно левая часть равенства (23) неотрицательна, откуда следует оценка
I J ач,ы(х)&з£к1 <1х 2 Ац J £^Х+<їх 1 — § \Вх+\'
П \П П/п
Далее, возьмем функции (й, Х) такие, что <ііуй = 0, Л(й, Х) > 0. Это возможно ввиду выполнения условия Аіу = ебіу. А так как пространство БУ (П) плотно вкладывается в Ь1(П), можно подобрать такую функцию Х, что
Р = J\Ех+\ < ж. (24)
п
Доказательство леммы завершено. □
Сформулируем результат о зависимости состояний равновесия функционала (1) на множестве (4) от параметров (Ь, \а\) в виде теоремы.
Теорема 5. Пусть минимум функционала (1) при \а\ = а+ + а— > 0, Ь = Ь* + є на множестве (4) реализуется на паре (й, Х) и пусть а+ ы £+ — а— к1 С— = сдц. Тогда существует финитная функция ао = ао(Ь), выпуклая на каждом из интервалов
(—ж,Ь*), (Ь*, ж), такая, что на множествах Mi, г = 1, 2, 3, задаваемых (15), выполняется следующее: х+ ф 1 если (Ь, \а\) € М1, х— ф 1 если (Ь, \а\) € М2, 0 ф х+ Ф 1 если (Ь, \а\) € М3. Кроме того, ао(Ь) ^ ао(Ь*), где для ао(Ь*) верны оценки (16), (22), а неравенства (12), (13) дают двустороннюю оценку для носителя функции ао:
(“2Щ’ Щ1 с 8ирра° -** с (е-’е+)' (25)
5. Зависимость температуры фазовых переходов от размера области и определение коэффициента поверхностного натяжения. В этом параграфе будем полагать, что вместе с областью П Э 0 задано семейство пропорциональных областей Пд:
Пд = {у € Мп : у = Их, х € П}, (26)
в каждой из которых можно рассматривать задачу на нахождение минимума функционала
Пд] = /"(¥+(Уй,Ь)х+ (у)+¥ —(Уи,Ь)х—(у)) йу + а / \Вх\ (27)
пл пл
на множестве
Бд = {(и, х): и € ’2(Пд), х.ф. х± € БУ(Пд), х+ и = 0, у € Пд}, (28)
где плотности энергии деформации задаются равенствами (3), а = (а+,а—), х = (х+,х—) определяют местоположение фаз и выполняется условие несмешиваемости
х+ + х— = 1 при п.в. у € Пд.
Сведем вариационную задачу в Пд к изученной задаче в области П. Пусть
ид(у) = Еи(х), х±(у) = х±(х), где у = Ех, х € П, у € Пд.
Тогда
I[ин, ХдЛ ^д] = ■ (29)
Кроме того, если (ид, хд) € Вд, то (и, х) € ^1. В силу этого обстоятельства, задача по определению минимума функционала энергии в множестве Пд при температуре Ь и коэффициенте поверхностного натяжения а сводится к задаче по определению минимума этого функционала в области П при той же температуре и коэффициенте поверхностного натяжения а/Е. Такую задачу мы изучали в предыдущем параграфе и проследили зависимость характера состояния равновесия от параметров. Соотношение (29) показывает, что изменение размеров области Пд сводится к изменению коэффициента а. Переформулируем теорему 5 в терминах изменения коэффициента поверхностного натяжения.
Теорема 6. Пусть минимум функционала (27) при Е = 1, \а\ = а+ + а— > 0, Ь = Ь* + е на множестве (28) реализуется на паре (и,х) и пусть Л— = а+ ыСк1 — а— ык = сб—. Тогда существуют выпуклые функции Ь+ (\а\) и — Ь— (\а\) такие, что х+ ф 1 при Ь > Ь+(\а\), х— ф 1 при Ь < Ь—(\а\), 0 ф х+ Ф 1 при Ь— (\а\) <Ь < Ь+(\а\).
Кроме того, |Ь±(0)| < то и
Ь+(0) < Ь* + е+, Ь (0) > Ь* + е—, (30)
где е± задаются равенствами (13).
Функции Ь±(\а\) = Ь* при \а\ ^ а*, где Ь* определяется тождеством (8), а на а* есть двустороняя оценка
А 4 4 *(П) (!ДгМи!) (шу'", (31)
в которой к(П) —константа в изопериметрическом неравенстве (21), числа Л, Р, задаются в (10) и (24).
Функции Ь±(\а\) определяют температуру фазовых переходов для состояний равновесия функционала (27) при Е =1. Если же мы захотим рассмотреть температуры
фазовых переходов Ь±(\а\) для функционала (27) при произвольных Е, то из соотно-
шения (29) будет следовать, что
_+± (\а\
На основании последнего равенства можно построить алгоритм по определению (оценке) коэффициента поверхностного натяжения. Будем полагать, что ІаI = ао нам неизвестен. Предположим также, что мы можем экспериментальным путем определить температуры фазовых переходов, т. е. считаем, что Ь±(ао) —наблюдаемые величины. Тождество (32) и теорема 6 говорят, что при малых значениях К (т. е. в малых областях) при изменении температуры от +то до —то состояние вещества всегда остается однофазовым (ЬД(ао) = Ь±(ао/К) = Ь* при ао/К > а*). Начинаем увеличивать параметр К, продолжая проводить эксперименты в области Пд. При постепенном увеличении параметра К наступит момент, когда в состоянии равновесия будут присутствовать обе фазы (Ь+(ао) = (ао/К) > Ь-(ао/К) = ЬД(ао) при ао/К < а*). Этот момент мы фик-
сируем. Таким образом нам удалось определить число К* такое, что при К < К* при любых температурах все состояния равновесия функционала энергии однофазовые, а при К > К* будут наблюдаться двухфазовые состояния равновесия. Следовательно ао/К* = а* и, вспоминая оценку (31) мы делаем вывод, что
*^<ВД(]ММ)(И)-. (33)
Кроме того, заметим, что эксперимент по определению температуры фазовых переходов можно ставить не при всех температурах Ь Є М, а только на интервале
Ь + < Ь < Ь + £+•
Автор признателен В. Г. Осмоловскому за постановку задачи.
A. S. Mikhailov. On determination of the surface tension coefficient in the two-phase elasticity problems under the assumption of incompressibility or with inflexible inclusions.
The homogeneous anisotropic two-phase elastic medium with one incompressible or absolutely inflexible phase and a shear field vanishing at the boundary is considered. The detailed description of the dependence of equilibrium states on the temperature and the surface tension coefficient is obtained. The method of determining the surface tension coefficient in the problem of phase transitions is suggested.
Литература
1. Гринфельд М. А. Методы механики сплошных сред в теории фазовых превращений. Москва, 1990. 180 с.
2. Осмоловский В. Г. Вариационная задача о фазовых переходах в механике сплошной среды. СПб., Изд-во СПбГУ, 2000. 260 с.
3. Михайлов А. С. Фазовые переходы в многофазовых средах при учете поверхностной энергии границы раздела фаз с помощью интеграла от старших производных поля смещений // Проблемы мат. анализа, 2001. Т. 22. С. 106-126.
4. Осмоловский В. Г. Об определении коэффициента поверхностного натяжения в механике двухфазовых упругих сред // Проблемы мат. анализа, 2005. Т. 30. C. 61-68.
5. Осмоловский В. Г. Зависимость состояний равновесия двухфазовой упругой среды от температуры при положительном коэффициенте поверхностного натяжения // Записки научных семинаров ПОМИ (в печати).
6. Осмоловский В. Г. Зависимость темпереатуры фазовых переходов от размеров области // Зап. научных семинаров ПОМИ, 2004. Т. 310. С. 98-114.
7. Осмоловский В. Г. Зависимость состояний равновесия двухфазовой упругой среды от температуры при нулевом коэффициенте поверхностного натяжения // Проблемы мат. анализа, 2004. Т. 28. С. 95-104.
8. Михайлов А. С., Михайлов В. С. Фазовые переходы в многофазовых средах // Проблемы мат. анализа, 2000. Т. 20. C. 120-170.
9. Осмоловский В. Г. Фазовые переходы для модели упругой среды с остаточным напряжением // Проблемы мат. анализа, 1997. Т. 17. С. 153-191.
10. Осмоловский В. Г. Зависимость характера состояний равновесия двухфазовой упругой среды от параметров задачи // Записки научных семинаров ПОМИ, 2000. Т. 271. С. 156-175.
Статья поступила в редакцию 10 апреля 2006 г.
