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

254

.pdf
Скачиваний:
4
Добавлен:
15.11.2022
Размер:
10.82 Mб
Скачать

Методвычисленияповерхностнойэнергиииэнергииадгезииупругихтел

взаимодействия, способные совершать работу на определенных на граничной поверхности тензорных характеристиках неравномерности распределения перемещений в ее окрестности: nu(k) A(k) , …. .

Система (1.4)–(1.7) дополняется условием инвариантности по-

тенциальной

энергии по отношению к

жестким движениям

тела B, которое выражается в условии

симметрии тензора

(

)

 

P(k ) =P(k(1))

P(k(2)) (P(k(3)) ...) .

 

Замыкание системы (1.4)–(1.7) осуществляется путем использования соотношения

m раз

 

}

 

P(n) =P0(n) −Β(n)Θ+∑( mur) ... C(m,n).

(1.8)

m=1

При получении этого соотношения, а также аналогичного соотношение для энтропии

m раз

 

 

}

 

 

S =S0 +∑( mur) ... Β(m) +cε

Θ

(1.9)

m=1

 

T0

 

( сε – удельная теплоемкость материала при отсутствии деформаций)

учтено, что дифференциал функции F(k ) =F(k ) (T(k ) ,{ ju(k ) }) является

полным. При этом для ее приращения справедливо выражение, совпадающее по форме с (1.3). Так что коэффициенты разложений (1.8), (1.9) определяются равенствами

0(n)

 

F

 

(m,n)

 

 

2 F

 

 

P

=

 

 

, C

 

=

 

,

(1.10)

( nu)

 

( mu)( nu)

 

Β(n) =−

 

2 F

 

 

=−P(n) .

 

(1.11)

 

T ( nu)

 

 

 

 

∂Θ

 

 

Вычисления в (1.10) и (1.11) должны проводиться для случая, ко-

гда T(k) =T0(k ) , nu(k ) =0 .

Используется также уравнение притока тепла

 

∂Θ

m

 

λ ( 2Θ)cε

T0

(

u)

m...разΒ(m) = ϑ.

(1.12)

t

t

 

 

m=1`

 

 

 

 

 

 

 

 

 

241

Л.Ю. Фроленкова, В.С. Шоркин

Здесь ϑ – с объемной плотностью скорости возникновения тепла. Начальным условием для него является требование Θ=0 . В любой последующий момент на границах тела задаются значения температуры или нормальной составляющей теплового потока.

Как отмечено в [11], вычисления можно провести, если известна зависимость объемной плотности свободной энергии от температуры и конфигурации деформированного состояния. В работе [10] предложено считать

F =Fr +F f +Fe .

(1.13)

Здесь второе и третье слагаемые – свободную энергию фононного и электронного газов соответственно – предложено вычислять по известным из физики выражениям [12]. Первое слагаемое в (1.13) – потенциальная энергия, определяется выражением

n1

j раз

}

Fr =Φ( j+1) (T ,L ,..., L )dV ...

(k ) (k , p) 0(k ) 1(k ,k ) n 1(k , p) 1(k ) j=1 V

Потенциал парного взаимодействия

dVj ( p)

Φ(2)(k ,k )

V

UV

, k p,

, V = (k )

( p)

 

V(k ) ,

k = p.

можно выбрать, на-

пример, в форме, аналогичной форме потенциала Морзе [9, 13] в физике твердого тела,

Φ(2)(k ,p) (2)0(k ,p) (e( k , p ) L1( p ,k ) 2eβ( k , p ) L1( p ,k ) ).

Для n – частичного взаимодействия предполагается [14]

Φ((nk ),p) (0(nk) ,p) Πn1 (e( k , p ) Lj ( p ,k ) 2eβ( k , p ) Lj ( p ,k ) ). j=1

Если k = p (рассматривается взаимодействие частиц из одного материала), то параметр β(k ,k ) предложено определить методами физи-

ки твердого тела [9, 15] через характеристики потенциального взаимодействия атомов и среднее межатомное расстояние l0(k ,k ) , определенное

при температуре T0(k ,k ) .

Параметры Φ(2)0(k ,k ) и Φ(3)0(k ,k ) можно выразить через модуль Юнга E(k ) и коэффициент Пуассона ν(k ) [9].

242

Методвычисленияповерхностнойэнергиииэнергииадгезииупругихтел

При этом предполагается, что переходной группы физических параметров между P(1) и P(2) нет.

2. Поверхностная энергия

Поверхностная энергия Wp равна изменению свободной энергии

рассматриваемого тела в изотермическом обратимом процессе образования единицы площади его свободной поверхности А [16]. Суммарное изменение свободной энергии Ηp при образовании всей поверхности А

(суммарная поверхностная энергия) определяется поверхностным интегралом

Ηp =Wp dA.

A

С другой стороны [17–22], это изменение свободной энергии оказывается сосредоточенным в трехмерном поверхностном слое. Ввиду неопределенности его границ [16] это изменение свободной энергии можно считать распределенным по всему объему V изучаемого тела B. Если Η – суммарная свободная энергия тела B, а ∆Η – ее изменение, возникшее при образовании A и распределившееся по области V с объемной плотностью w ( w =w(R) , R V , R – радиус-вектор точки об-

ласти V в конечном состоянии тела B, ∆Η=wdV ), то, приравняв H

V

и Hp , можно получить интегральное уравнение для определения Wp (при условии известности распределения w по области V):

Wp dA=wdV.

(2.1)

A

V

 

Обычно считается, что вдоль поверхности A величина Wp не ме-

няется и является физической характеристикой материала (например, [23, 24]). В физических расчетах [1, 22] область V в прямоугольных декартовых координатах (x1, x2 , x3 ) задается в виде

V: 0< x x1 <+∞ , −∞< x2 , x3 <+∞; A: x x1 =0 .

(2.2)

В данной работе сформулируем задачу по расчету распределения Wp вдоль криволинейной поверхности A.

243

Л.Ю. Фроленкова, В.С. Шоркин

Функция Wp =Wp (R) определена на двухмерной гладкой замкну-

той области A. Для этой функции всегда можно найти такую векторную функцию a =a(R) , которая определена на трехмерной области V, непрерывна и дифференцируема в каждой внутренней ее точке, а на границе A удовлетворяет условиям

Wp =n a , an (n a) =0 .

(2.3)

Подставляя первое из равенств в левый интеграл равенства (2.1), можно получить

n a dA=wdV ,

diva dV =wdV.

(2.4)

A

V

V

V

 

Для векторного поля a =a(R)

всегда можно найти такое скаляр-

ное поле Ψ =Ψ(R) , R V , [25], что

 

 

 

a =grad Ψ,

Ψ(R)R A =0 .

(2.5)

Равенство (2.1) считается справедливым для сплошных тел B произвольных размеров с произвольной формой гладкой границы A. Поэтому учет (2.5) для второго из интегралов (2.4) приводит к классической краевой задаче для функции Ψ =Ψ(R) :

divgrad Ψ=w , Ψ(R)R A =0 .

(2.6)

Задача (2.6) имеет решение и притом единственное. Значит, существует и единственно решение задачи (2.3)–(2.6) по определению распределения поверхностной энергии вдоль свободной от внешних воздействий поверхности A.

В частном случае, когда область V определена выражением (2.2) (полубесконечное тело), а значение w зависит лишь от расстояния x до плоской границы A, решение задачи (2.3)–(2.6) имеет вид

Wp =w(x)d x .

(2.7)

0

 

Это выражение совпадает с тем, которое используется для расчета поверхностной энергии методами физики твердого тела [1, 22].

244

Методвычисленияповерхностнойэнергиииэнергииадгезииупругихтел

Вычисление поверхностной энергии путем решения задачи (2.3)– (2.6) или с помощью (2.7) можно осуществить, если известно распределение w =w(R) , R V .

Сразу после мгновенного выделения B(k ) из (k ) взаимное распо-

ложение его частиц и температура остались теми, которыми они обладали в составе (k ) . Поэтому различие свободных энергий тела B(k )

в состояниях перед реальным выделением и после него обусловлено только разницей потенциальных энергий, возникшей в теле B(k ) за счет

удаления частиц,

принадлежащих

дополнению

B(k )

в составе

(k ) .

С учетом этого

 

 

 

 

 

 

( j+1)

,...,ln1(k ,k ) )dV1(k )

 

 

 

n1

Φ(k ,k ) (T0(k ) ,l1(k ,k )

... dVn1(k )

 

w(k ) =( k )

 

 

.

(2.8)

j=1

( j+1)

 

 

 

 

Φ(k ,k ) (T0(k ) ,L1(k ,k ) ,...,Ln(k ,k ) )dV1(k ) ... dVn1(k )

 

V( k )

 

 

 

 

Момент выделения B(k ) из (k ) служит началом переходных про-

цессов. Эти процессы, а также распределение w(k ) (r)

в каждый момент

времени можно рассчитать после решения начально-краевой задачи

(1.4)–(1.12).

Предполагается, что решение системы (1.4)–(1.12) имеет предел при t →∞ , совпадающий с решением системы, получаемой из (1.4)– (1.12) путем исключения из нее слагаемых, содержащих производные по времени.

Если обмена энергии у тела B(k ) после его образования не было,

то суммарная поверхностная энергия в любой момент времени является одной и той же.

3. Адгезия твердых тел

Процесс вступления в адгезионный контакт тел B(1) и B(2) можно

условно разбить на этапы.

Первый этап – образование участков свободных поверхностей: A(1,2) у тела B(1) и A(2,1) у тела B(2) , вдоль которых происходит адгезион-

ный контакт. В процессе первого этапа у каждого из тел B(1) и B(2) формируются их поверхностные энергии Wp(1) и Wp(2) .

245

Л.Ю. Фроленкова, В.С. Шоркин

Второй этап – «жизнь» тел B(1) и B(2) с «намеченными» участками контакта A(1,2) и A(2,1) без взаимодействия между собой. В процессе

второго этапа тела могут подвергаться внешним воздействиям. Из-за этого могут меняться их термоупругое состояние, физические свойства, свободная энергия. А следовательно, и поверхностная энергия Wp(k) .

Третий этап начинается, когда реального контакта поверхностей еще не произошло, но тела B(1) и B(2) приблизились друг к другу бес-

конечно близко. Тела B(1) и B(2) начинают воздействовать друг на друга

посредством силовых полей.

Четвертый этап заключается в формировании адгезионного контакта. Он начинается в момент касания поверхностей A(1,2) и A(2,1) и за-

канчивается формированием трехмерного переходного слоя B(1,2) между телами B(1) и B(2) , не имеющего четких границ. Далее предполагает-

ся следующее.

Считается, что поверхности A(1,2) и A(2,1) конгруэнтны. Случай не-

конгруэнтности анализируется в работах [2, 4, 5], далее не рассматривается.

Переходный слой является трехмерным, в специальную область с особыми термоупругими свойствами не выделяется. Случай, когда его толщина h фиксируется с постулированием распределения свойств материала по нему или когда h =0 , а свойства материала при пересечении слоя меняются скачком, не рассматривается. Большое внимание рассмотрению таких слоев уделено в работах [7, 8, 26]. Не рассматривается также случай, когда одно из тел «выращивается» на поверхности другого [27].

Допускается, что тела B(1) и B(2) выделяются из бесконечных сред

(1) и (2) одновременно в момент времени t =0 . В момент t =0+0 они

сначала мгновенно сблизились на бесконечно близкое расстояние, а затем вступили в адгезионный контакт сразу вдоль обеих поверхно-

стей A(1,2) и A(2,1) .

Мгновенность выделения B(1) и B(2) из (1) и (2) с мгновенным их сближением и соединением обеспечивает адиабатичность процесса

246

Методвычисленияповерхностнойэнергиииэнергииадгезииупругихтел

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

нечных средах до вступления в адгезионный контакт, можно связывать исключительно с процессами отделения B(k) от (k) и образования из-

за адгезии единого тела B =B(1) U B(2) .

У частиц тел B(1) и B(2) в этот момент вновь (после выделения)

мгновенно меняется число окружающих их соседей, при этом деформации и изменения температуры отсутствуют. Каждая частица dB(k) B(k), k =1,2 , мгновенно, наряду с взаимодействием с другими

частицами dB(k)1 B(k)

и парами частиц dB(k)1, dB(k)2 B(k),

начинает

взаимодействовать с частицами

dB(p)1 B(p),

p =1,2 , p k

и парами

частиц dB(k)1 B(k),

dB(p)2 B(p),

а также

dB(p)1 B(p),

dB(k)2 B(k)

и dB(p)1,dB(p)2 B(p) . Потенциалы этих взаимодействий предполагаются

известными.

Изменение их свободной энергии, происходящее в ходе этого процесса, определяется выражением, аналогичным выражению (2.8), но более общим, учитывающим взаимное влияние частиц разных тел:

 

 

 

 

 

 

=

2

 

 

 

F(1,2)

Fa d A(1,2) =∑ ∫ wa(k ) dV(k ) =

 

 

 

 

 

 

 

( )

( k )

 

 

 

 

 

 

 

A 1,2

k=1 V

 

2

n1

Φ((kj+,k1)) (T0(k ) ,l1(k ,k ) ,...,ln1(k ,k ) )dV1(k ) ... dVn1(k )

=∑ ∑

 

k=1

j=1

( k )

 

 

 

2

n1

 

 

1

 

 

 

 

∑ ∑

 

 

 

Φ((kj+,p1)) (T0(k ) ,L1(k ,p) ,...,Ln(k ,p) )dV1(k ) ... dVn1( p) .

 

( j +1)!

k , p=1

j=1

 

 

 

 

 

 

 

 

V

 

 

 

 

 

 

 

 

( k )

 

 

После мгновенного присоединения тел B(1) и B(2) друг к другу

в следующий за этим момент времени при их адгезии начинает формироваться переходный слой B(1,2) . При вступлении в состояние ад-

гезии на поверхности контакта, по предположению, мгновенно начинают выполняться условия непрерывности термодинамических параметров.

247

Л.Ю. Фроленкова, В.С. Шоркин

Далее рассматривается объединенное тело B =B(1) U B(2) , занимающее в отсчетном состоянии область V =V(1) UV(2) , ограниченную со стороны B(1) поверхностью Aa (1) , а со стороны B(2) поверхностью Aa (2) . Тела граничат между собой вдоль поверхности Aa (1,2) .

Пусть wa(k ) =FaV (k ) F(k ) – объемная плотность изменения свободной энергии тела B(k ) к моменту времени t, отсчитываемому от момента его выделения из (k ) , а Wap(k) – его поверхностная энергия (вдоль поверхности Aa (1,2) ). Из-за соединения тел B(k) в единое целое

wa(k ) w(k ) и Wap(k) Wp(k) как на участках Aa (k) , так и, тем более, на участке Aa (1,2) . Несмотря на это, для каждой из частей B(k) , k =1,2, объединенного тела B =B(1) U B(2) справедливы соотношения (2.1)–(2.7) с той

лишь разницей, что входящие в них величины относятся к состоянию адгезии. При этом выражение (2.7) относится только к части Aa(k) по-

верхности тела B(k) . На оставшейся его части A(1,2) должно выполняться условие сопряжения для полей Ψa(k ) (t,r) . Если контакт есть, но адгезии нет, то сопряжения полей Ψa(k ) (t,r) нет. Таким образом, задача о распределении поверхностной энергии тел B(k ) в случае, когда они

находятся в состоянии адгезии, приобретает вид (2.4), (2.5).

Условия сопряжения, накладываемые на распределения Ψa(k ) (t,r) , предлагается принять в виде

 

Ψa(1) (t,r) a(2) (t,r) , r A(1,2) ;

(3.1)

Ψ

a(1)

(t,r) =

Ψ

a(2)

(t,r) r A .

(3.2)

n

n

 

 

 

(1,2)

 

Условие (3.1) соответствует условию непрерывности изменения характеристик материала и его термодинамического состояния при пересечении переходного слоя. Условие (3.2) соответствует равенству поверхностных энергий тел B(k ) по обе стороны поверхности A(1,2) .

В соответствии с методом сжатия переходного слоя [28] это является

248

Методвычисленияповерхностнойэнергиииэнергииадгезииупругихтел

условием определения поверхности A(1,2) в объединенном теле B =B(1) U B(2) . В общем случае она не совпадает с исходным положени-

ем. В рассматриваемой ситуации это совпадение допускается.

Если сложить найденные величины Wap(1) и Wap(2) для точек поверхности A(1,2) , то можно получить приходящееся на единицу ее площади изменение свободной энергии Wp(1,2) системы B =B(1) U B(2) по отношению к свободной энергии совокупности несоединенных тел B(k) .

При этом энергия адгезии Fa =Wp(1,2) Wp(1) Wp(2) .

Для проведения расчетов поверхностной энергии Wap(k ) (t,r) необходимо произвести расчет распределения wa(k ) (t,r) . Для этого реша-

ется сопряженная задача следующего вида.

В каждой из областей V(k ) выполняются балансовые соотношения

для импульса (1.4) и энергии (1.12) с учетом выражений (1.8)–(1.11). В выражениях (1.8)–(1.11), в отличие от ситуации с одиночными телами B(k ) , необходимо учесть их взаимное влияние. Предполагается, что

в любой момент времени после вступления тел в адгезионный контакт взаимное влияние осуществляется только за счет потенциальных полей, создаваемых частицами тел B(k ) . Таким образом, в выражениях

(1.10), (1.11) необходимо учитывать, что

2 n1

F =Fe +F f +∑ ∑ Φ( j+1) (T ,L ,...,L )dV ... dV .

(k ) (k ,p) 0(k ) 1(k , p) n(k , p) 1(k ) n 1( p) p=1 j=1 V( k )

Здесь, как и ранее, величины Fe и F f определяются из [12].

На свободных от взаимного контакта участках границ Aa(k ) вы-

полняются краевые условия для внутренних напряжений (1.6), (1.7) или перемещений, а также условия, накладываемые на поле темпера-

тур T

=T

или тепловых потоков λ

dT(k )

=q

(T

, q

должны

 

(k )

A(k )

(k ) d n

A(k )

A(k )

A(k )

 

 

 

 

(k )

 

 

 

 

быть заданы). На поверхности A(1,2) задаются условия сопряжения.

249

Л.Ю. Фроленкова, В.С. Шоркин

 

 

 

 

 

 

∂Θ

T

 

 

 

∂Θ

T

 

 

Θ

+T

+T

; k

 

(1)

+

0(1)

 

=k

 

(2)

+

0(2)

 

; (3.3)

 

 

 

 

(1)

0(1)

(2)

0(2)

(1)

 

n

n

 

(2)

 

n

n

 

 

 

 

 

 

 

 

 

 

 

 

 

u(1) (r)u(2) (r)=τuτA (r)+nun A (r), r A(1,2) ;

 

n (P(1)(1) P(1)(2) )A (n P(1)(2) )=n (P(2)(1) P(2)(2) )A (n P(2)(2) );

(3.4)

u

(rr)u

2)

(r)=С%(r), С%(r)=с(r)(nn)+ε(r)(ττ) ; n P(2) =n P(2) ;

(3.5)

(1)

(

 

(1)

(2)

 

 

 

 

S(1) =S(2) ; F(1) =F(2) .

 

(3.6)

 

Варьируя функционал, определяющий свободную энергию объе-

диненного тела B =B(1) UB(2) ,

а также его частей B(1) и B(2) ,

можно по-

лучить выражение

 

 

 

 

{(n(1,2) q(1) )δT(1) +(n(2,1) q(2) )δT(2)} + { n(1,2) (P(1)(1) P(1)(2) )A ×

 

 

×(n(1,2) P(1)(2) ) δu(1) + n(2,1)

(P(2)(1) P(2)(2) )A (n(2,1) P(2)(2) ) δu(2) } +

(3.7)

 

+{(n(1,2) P(1)(2) ) δ( nru(1) )+(n(2,1) P(2)(2) ) δ( nru(2) )} +...=0.

 

Допускается, что адгезионный контакт наступает мгновенно. При этом на границе контакта мгновенно начинают выполняться условия адгезии – непрерывность изменения характеристик термодинамического состояния при переходе через поверхность контакта. Это означает, что, во-первых, на поверхности контакта выполняются равенства (3.6). Во-вторых, что вариация свободной энергии равна нулю. В-третьих, что вариации параметров термодинамического состояния материала по обе стороны от поверхности контакта равны между собой как соответствующие одному и тому же состоянию элементарной частицы, разделенной поверхностью контакта. Равенство вариаций

δT(1) =δT(2) ; δu(1) =δu(2) ; δ( n u(1) ) =δ( n u(2) )

(3.8)

означает, что варьируемые величины могут отличаться друг от друга на заранее известную функцию. Это приводит к появлению первых равенств в каждой из пар (3.3)–(3.5). Вместе с тем произвол каждой из вариаций в совмещении с равенствами (3.7) и (3.8) приводит к появлению вторых равенств в каждой из пар (3.3)–(3.5).

250

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]