Теория термоупругости, позволяющая вычислять поверхностную

advertisement
1
ТЕОРИЯ ТЕРМОУПРУГОСТИ, ПОЗВОЛЯЮЩАЯ ВЫЧИСЛЯТЬ
ПОВЕРХНОСТНУЮ ЭНЕРГИЮ И ЭНЕРГИЮ АДГЕЗИИ
Шоркин В. С., Фроленкова Л. Ю.Якушина С.И.
ФГОУ ВПО «Госуниверситет - учебно-научно-производственный комплекс»
302020 Орел, Россия
E-mail: VShorkin@yandex.ru, LaraFrolenkova@yandex.ru
Адгезия твердых тел
Физический энциклопедический словарь, М.: Советская энциклопедия,1960.
«Адгезия (прилипание) – возникновение связи между поверхностными слоями двух разнородных
(твердых или жидких) тел (фаз), приведенных в соприкосновение»
«С термодинамической стороны адгезия характеризуется свободной энергией адгезии Fa – убылью
свободной энергии на 1 см 2 поверхности адгезионного шва в изотермическом обратимом процессе
адгезии, и выражается через поверхностные натяжения
W p ( 1) , W p ( 2 ) , W p ( 1, 2 ) »
A = − Fa = W p ( 1, 2 ) − W p ( 1) − W p ( 2 )
A = − Fa
(1)
– работа адгезионного отрыва.
Характеристика соединения покрытия и подложки, измеряется при нормальном отрыве или
касательном сдвиге в единицах силы на единицу площади (МПа)
2
Е.М. Лифшиц. Теория молекулярных сил притяжения между твердыми телами, ЖЭТФ, 29, 94,
1955.
«Наличие сил притяжения между нейтральными атомами приводит , естественно, к появлению
аналогичных сил и между двумя макроскопическими телами, поверхности которых сближены до
очень малых расстояний. Вычисление этих сил, исходя из известного взаимодействия отдельных
атомов, было бы возможно, однако, лишь для достаточно разреженных тел, т.е. газов (ГАМАКЕР) –
случая, который фактически, разумеется, не может быть осуществлен. К этому вопросу можно,
однако, подойти и чисто макроскопическим образом (поскольку расстояние между телами
предполагается большим по сравнению с междуатомными расстояниями). При этом взаимодействие
тел рассматривается как осуществляющееся через посредство флуктуационного электромагнитного
поля, которое всегда присутствует внутри всякой поглощающей среды и выходит также и за ее
пределы – частично в виде излучаемых телом бегущих волн, частично в виде стоячих волн,
экспоненциально затухающих с удалением от поверхности тела. ….
Основанный на таком рассмотрении метод вычисления сил взаимодействия обладает полной
общностью.»
«Этот метод основан на введении в уравнения Максвелла стороннего «случайного» поля. … Сила F
взаимного притяжения, действующего на единицу поверхности каждого из тел мы вычисляем как xx
– компоненту максвелловского тензора напряжений.
МАКСВЕЛЛА ТЕНЗОР НАПРЯЖЕНИЙ -:
3
где
- компоненты электрич.
и магн.
полей в вакууме,
- символ Кронекера,
Смысл понятия «поверхностная энергия»
Л.И. Седов Механика сплошной среды, т. 2.
«В классической термоэластике полная внутренняя энергия термоупругого тела представляется в виде
UT = ∫ U ( ε ij , s ) ρ dV + U 0 = U1 + U 0 ,
(2)
V
где U (ε ij , s ) – некоторая определенная функция, зависящая от удельной энтропии s и компонент
тензора деформации ε ij , U 0 – аддитивная постоянная. В задачах «чистой» теории упругости
аддитивная постоянная несущественна…
При рассмотрении развивающихся трещин – … необходимо учитывать … виды энергии, связанные с
поверхностными эффектами, проявляющимися при нарушении целостности тела. Простейший
способ учета таких эффектов можно осуществить с помощью аддитивной постоянной U 0 .»
4
«В дальнейшем будем трактовать U 0 как внутреннюю энергию, которая связана с понятием
поверхностной энергии, но всё же не сводится просто к этому понятию за счет более глубокого
изменения структуры частиц материала в тонком поверхностном слое.»
«Опыт и некоторые общие физические соображения позволяют сделать допущение, что
n
dU 0 = ∑ γ i d Σ i ,
i =1
(3.5)
где d Σ i – приращения площадей поверхностей трещин в различных частях внутри тела, а γ i –
соответствующие функции, определенные в местах образования площадок разрыва d Σ i
(поверхностная плотность энергии на разрыв). Для тел, разрушающихся хрупким образом, во многих
случаях можно считать, что γ есть просто плотность поверхностной энергии, подобная плотности
энергии поверхностного натяжения для жидкостей»
«Иногда γ нельзя отождествлять с плотностью поверхностной энергии твердого тела. Опыты
показывают, что плотность энергии на разрыв
∆ U0
γ эф = lim
.
∆Σ→ 0 ∆ Σ
Во многих случаях значительно превышает плотность поверхностной энергии (например, за счет
затраты дополнительной энергии на образование пластических деформаций в окрестности трещин)»
«При изучении проблем прочности экспериментальное и, может быть, теоретическое
исследование величины γ должно составлять главную задачу»
5
2. Крокстон К. Физика жидкого состояния. Статистическое введение. М.: Мир, 1978. 400 с.
"Следуя Кирквуду и Баффу, рассмотрим гетерогенную систему, состоящую из двух фаз, и введем
прямоугольную систему координат так, чтобы поверхность раздела Гиббса лежала в плоскости x y .
При этом ось z нормальна к поверхности жидкости и направлена из жидкой фазы в газовую.
y
z
P − P( z )

r1

r2
x1 , y1 , z1
r
x
x2 , y 2 , z 2
l =1
d
2π ρ 2 ∞
3
∇
≡
P = ρkT −
∇ 1Φ ( r ) g ( 2 ) ( r ) r d r – давление в однородной фазе, 1
;
d r1
3 ∫0
6
P ( z1 )
γ=
1 x 22

= k T ρ ( 1 ) ( z1 ) − ∫
∇ 1 Φ ( r ) ρ ( 2 ) ( z1 , r ) d 3 r2 – текущее давление;
2 r
+∞
∫ [ P − P( z) ] d z
– поверхностное натяжение.
−∞
«Фаулер … Кирквуд и Бафф … сузили переходную область до математической поверхности разрыва
плотности, совпадающей с поверхностью раздела Гиббса. … Для поверхностного натяжения …
Кирквуд и Бафф получают следующее выражение
π ρ 2L ∞
γ=
∇ 1Φ ( r ) g(L2 ) ( r ) d r ,
∫
8 0
а для добавочной поверхностной энергии – выражение
π ρ 2L ∞ L
uS = −
g( 2 ) Φ( r ) r 3 ( r ) d r
∫
2 0
Кроме того, имеется связь
u
d (γ / T )
= − S2 , γ ( T ) = γ 0 (1 − T / Tc )µ ,
dT
T
Для КБ- модели всегда γ ( скачк ) > γ ( эксп ) , u S ( скачк ) > u S ( эксп ) .
»
µ = 1,24 »
Л.Д. Ландау, Е.М. Лифшиц. Теоретическая физика. т.V. Статистическая физика. Часть. I. ,
М.: Наука, 1976. 584 с.
7
«Работа, затрачиваемая при обратимом изменении площади Σ на величину d Σ , может быть написана
в виде
d R = αd Σ
(154.1)
… α …называется коэффициентом поверхностного натяжения.
Формула в точности соответствует формуле d R = − p d V для работы при обратимом изменении
объема тела.
Можно сказать, что α играет для поверхности такую же роль, как давление для объема»
«Рассматривается система двух изотропных фаз жидкость-пар одного и того же вещества.
Без учета поверхностных эффектов дифференциал энергии системы из двух фаз… при заданном
объеме V всей системы имеет вид dE = TdS + µdN (в равновесии температуры T и химические
потенциалы µ обеих фаз одинаковы…). При учете же наличия поверхностных эффектов к правой
стороне этого равенства надо, очевидно, прибавить еще выражение (154.1):
d E = T d S + µ dN + α d Σ .
(154.2)»
Все характеристики термодинамического состояния делятся на объемные и поверхностные. Но при
этом возникает неопределенность из-за неопределенности толщины поверхностного слоя и
количества попавших в него частиц.
«Мы сделаем разделение однозначным, подчинив его следующему естественному условию: объемы
V1 и V2 обеих фаз определяются так, чтобы наряду с равенством V1 + V2 = V имело место также и
8
равенство n1V1 + n2V2 = N . Эти два равенства фиксируют … объемные значения всех …
термодинамических величин».
Методы вычисления поверхностной энергии
Магомедов М.Н. О зависимости поверхностной энергии от размера и формы нанокристалла.
ФТТ, 2004, т. 46, № 5, С. 924 – 937.
«Рассмотрим систему из N одинаковых атомов при температуре T и давлении P . Изменение
свободной энергии такой системы при варьировании температуры, объема V , числа атомов и
площади поверхности Σ равно
d F ( T ,V , N , Σ ) = − S d T − P d V + µ d N + σ d Σ .
(1)
Здесь S , µ и σ – энтропия, химический потенциал удельная (на единицу площади) свободная
поверхностная энергия. … Согласно (1), функция σ ( T ,V , N ) определяется выражением
σ ( T ,V , N ) = ( d F / d Σ ) T ,V , N .
(2)
Основная проблема всех теоретических подходов заключалась в нахождении производной ( d F / d Σ )
именно при постоянных значениях объема и числа атомов в системе» … Автор получает выражение
свободной энергии нанокристалла, позволяющее вычислить производную (2).
Партенский М.Б. Самосогласованная теория металлической поверхносчти. УФН, 1979, т. 128,
в. 1, С. 69 – 106.
9
«Поверхностная энергия есть работа по разделению кристалла на две части вдоль некоторой
плоскости приходящаяся на единицу площади вновь образуемой поверхности. Если исходный
кристалл макроскопичен и образуемые фрагменты эквивалентны, то
{
[ ]
[ ]}
σ = ( 2Σ )− 1 2 (G [ n] + Ees [ n] ) − (G n∗ + Ees n∗ ,
(4.4)
n и n∗ – электронные плотности соответственно во фрагментах и в исходном кристалле, E es –
полная классическая кулоновская энергия положительных и отрицательных зарядов,
 
1 n( r ) n r ∗
 
G [ n] = Ψ T + V Ψ − ∫∫  ∗ ( d r ) d r ∗
2
r −r
( )
(
)
– функционал плотности, T , V описывают кинетическую энергию и кулоновское взаимодействие
электронов, Ψ – вектор основного состояния системы»
Векилов Ю.Х., Вернер В.Д., Самсонова М.Б. Электронная структура поверхностей
непереходных металлов. УФН, 1987, т. 151, в. 2, С. 341-376.
«В последнее время при исследовании свойств поверхностей металлов всё большее внимание
уделяется изучению геометрии расположения ионов вблизи поверхностного слоя, которая может
значительно отличаться от объемной и изменять физические свойства поверхности» …
«Исправить недостатки модели парных потенциалов можно введением многочастичного
взаимодействия: показано, что уже при учете трехчастичного взаимодействия ионов получается
10
значение ∆ 12 (расстояние между первым и вторым атомными слоями), соответствующее
экспенименту»
Выводы
1.
2.
3.
4.
5.
Возможность вычисления энергии и силы адгезии, а также свободной энергии твердого тела
имеется.
В начальном и конечном состояниях, поверхностную энергию, энергию адгезии предлагается
вычислять как разность ее значений, отнесенную к единице площади образовавшейся
свободной поверхности, или как производную по величине ее площади.
Температуры элементов тела до разделения и после него, или в глубине тела и вблизи
поверхности и на ней, считаются одними и теми же.
Отмечена необходимость многочастичного взаимодействия в поверхностном слое.
Обращено внимание на сложность учета объемных эффектов в поверхностном слое и
удобство гипотезы о нулевой толщине поверхностного слоя распределения поверхностной
энергии
11
Сделано следующее
1. Допускается, что материал, для которого определяется поверхностная энергия является сплошным,
однородным, изотропным, термоупругим.
2. Считается, что любое тело B получается путем мгновенного изотермического выделения из
бесконечно протяженной среды (тела) Ω . Отсчетным для него является то термодинамическое
состояние, которым оно обладало бы, находясь в составе тела Ω, отсчетной конфигурацией –
конфигурация внутри Ω . При этом внешние воздействия отсутствуют, температура T равна Т0.
Свойства тела Ω являются известными. Тело B занимает в отсчетном состоянии ( t = 0 ) область V
(объем), а в текущем – Vt = V (t ) ( t > 0 , t – момент времени, в который тело B находится в текущем
состоянии). Границами этих областей являются гладкие поверхности (площади) А и At = A( t ) .
3. Предполагается, что процесс выделения B из Ω происходит мгновенно, изотермически. После
этого часть C = Ω \ B , дополняющая B в Ω , мгновенно удаляется от B на бесконечность.
4. При вычислении поверхностной энергии σ возможны варианты.
4.1. Для полубесконечного тела B , ограниченного плоскостью A , поверхностная энергия σ при
t = 0 + 0:

1
σ = lim  ∫ ( F p∗ −F p0 ) dV  .
(1)
A→∞ A
V

∗
∗ 
0
0 
Здесь F p = F p (r ) , F p = F p ( r ) – объемные плотности распределения потенциальной энергии в теле
B сразу после его выделения ( t = 0 + 0 ) и непосредственно перед ним ( t = 0 − 0 ). Положения частиц

 

до выделения r и после него R , совпадают: R = r .
12
4.2. В случае криволинейной поверхности A :
∫ σ dA = ∫ w p dV , где w p = F p∗ − F p0 .
A
V
(2)
Тогда:
  
 
∃a = a ( r ) : σ = n ⋅ a ,
 
 
σ
dA
=
n
⋅
a
dA
→
σ
=
n
⋅ a и div a = w
∫
∫
(3)
p
A
A
 
∃ Ψ = Ψ ( r ) : a = grad Ψ
После этого поверхностная энергия определяется как решение задачи:
grad div Ψ = w p ,
(4)

Поскольку важна только нормальная составляющая вектора a считается:
 
r = rA → Ψ = const = 0 ,
(5)
После этого:

σ = n ⋅ gradΨ .
(6)
При t > 0 тело B деформируется, в нем меняется распределение температуры.
4.3. Допускается, что при t → ∞ динамический процесс затухает, если он проходит адиабатически, то


1
0
∀ t > 0 : σ = lim  ∫ ( F −F ) dV  , F = F ( R( t )) .
(7)
A→∞ A
V

5. Объемная плотность свободной энергии равна
F = F p + F f + Fe .
Для неметаллов Fe = 0 .
(8)
13
2




3
5
π
k
T
2
− hν / k T
 
Ff = ∑
d ν , Fe = ne µ 0  1 − 
(9)
∫ ν ln 1 − e
5
12
µ

j =1
 0  
0
n –среднее числу атомов и ne – число свободных электронов в единице объема; h и k соответственно
постоянные Планка и Больцмана; ν – возможное значение частоты колебаний атомов, заключенное в
пределах от нуля до частоты обрезания Дебая ν D ; ν j – частота моды с номером j , число которых в
единице объема равно 3 n ; µ 0 – энергия Ферми.

(
F
=
F
R
(t )) в любой момент времени определяется формулой
6. Величина p
p
3n
1
27 n k T
hν j +
2
ν 3D
νD
(
)
 1 ( 2)  

1
( 3)   
(10)
∫  2! Φ ( R1 , R2 ) + 3 ! ∫ Φ ( R1 , R2 , R3 ) dVt 3  dVt 2 ,

Vt 
Vt

 



R1 , R2 ) dVt 1 dVt 2 , Φ ( 3 ) ( R1 , R2 , R3 ) dVt 1 dVt 2 dVt 3 – потенциалы парного и тройного

F p ( R) =
( 2)
где Φ (
взаимодействия.


=
= F p ( R( t = 0 + 0 ) = F p ( r ) ,
(11)
 
  
  
( 2)
( 2)  
( 3)
когда Φ ( R1 , R2 ) = Φ ( r1 , r2 ) , Φ ( R1 , R2 , R3 ) = Φ ( 3 ) ( r1 , r2 , r3 ) , Vt = V .
(12)



F p0 = F p0 ( r ) = F p ( R( t = 0 − 0) = F p ( r ) ,
(13)
 
  
  
( 2)
( 2)  
( 3)
когда Φ ( R1 , R2 ) = Φ ( r1 , r2 ) , Φ ( R1 , R2 , R3 ) = Φ ( 3 ) ( r1 , r2 , r3 ) , Vt = Ω .
(14)
7. Выражения для потенциалов Φ ( k ) , k = 2, 3 , постулируются. Для изотропной однородной среды
принято
F p∗

F p∗ ( r )
Φ
( 2)
(
(
)
 
R1 , R2 ) = Φ ( 2 ) L1 2 ,
Φ
( 3)
(
)
  
( R1 , R2 , R3 ) = Φ( 3) L12 , L13 .
14
(15)
Вид функций при деформировании не меняется.
7.1. По аналогии с учетом ванн-дер-ваальсова взаимодействия частиц радиуса r0 = 0,5 l 0 с объемом
3
dVk ≈ ∆Vk = (4 / 3) π r0 дискретной среды принято:
Φ
( 2)
=
H ( 2)
( l ij + l0 )
6
, Φ
( 3)
=
(
)
H ( 3 ) 3 cos ϑ i cos ϑ j cos ϑ k + 1
( l ij + l 0 )
9
.
Здесь i, j, k – номера взаимодействующих частиц, i, j, k = 1, 2, 3; l 0 = 3 n = L0


( 2)
( 3)
векторами l ij и l ik ; H , H – константы.
(16)
 
R= r ,
ϑ i – угол между
2
 4π r0 3  ( 2 ) 3  ∞ 2


 3  H = π ∫ α ( i ω) d ω
0


3
,
(17)
 4π r0  ( 3 ) 3  ∞ 3


 3  H = π ∫ α ( i ω) d ω
0


Зависимость поляризуемости материала α( ω ) от частоты ω флуктуационных колебаний внутреннего
электрического поля известны.
( 2)
( 3)
7.2. Потенциалы Φ и Φ определяются формулами:
3
(
)
(
)(
)
Φ ( 2 ) = D ( 2 ) e − 2ξ l1 − 2e −ξ l1 , Φ ( 3 ) = D ( 3 ) e − 2ξ l1 − 2e −ξ l1 e − 2 ξ l2 − 2e −ξ l2 .
15
(18)
( 2)
( 3)
Параметры D , D , ξ определены, опираясь на данные механических экспериментов.
За время t → 0 после момента t = 0 + 0 тело B испытает изменение температуры


Θ ( r , t ) = T ( r , t ) − T0 и относительное смещение частиц
  
  
∆ u = u2 − u1 = u( r + l ) − u( r ) =
n раз
( −1) n n    n ,
∑ n! ∇ u l =0 ⋅ ... ⋅ l
n =1
(
∞
)
(19)
 
    
l ≡ l12 ≡ l 2 = r2 − r1 , r1 ≡ r – радиус-вектор центра инерции произвольной частицы dB2 относительно
частицы dB1 ≡ dB .

Получаем, что F = F T , ∇ n u . Полагая эту функцию полиномом второй степени от Θ = T − T0 , и

∇ n u , можно получить, что объемная плотность приращения свободной энергии при переходе из
промежуточного состояния в текущее определяется выражением
∞
n раз
1 ∞ m  m раз ( n,m ) n раз n  T 1 c ε 2
n
0( n )
nT
w = w Θ, ∇ u = S 0 Θ + ∑ P
⋅ ... ⋅ ∇ u + ∑ ∇ u ⋅...⋅C
⋅...⋅ ∇ u + 2 T Θ +
2
n =1
m , n =1
0
{
( {
}
( {
})
})
(
)
(
)
(
n раз
∞
1 ∞

( n)
nT
+  ∑ Θ B ⋅ ...⋅ ∇ u + ∑ ∇ m u
2  m , n =1
n =1
(
)
(
) ⋅...⋅B
m раз
( n) T

Θ .

)
(20)
16
При этом тензоры напряжений P ( n ) и энтропия S (в отсчетном состоянии считается равной нулю)
определяются формулами
∞
∂F
( n)
(
)
(
)
(
)
( m ,n )
n
0
n
n
m  m раз
P =
P
=
P
−
Β
Θ
+
∇
u
⋅

⋅
С
,
,
(21)
∑

n
∂∇ u
m =1
∞
Θ

S = ∑ ∇ m u ⋅  ⋅ Β ( m ) + cε .
(22)
T
m =1
При этом уравнения притока тепла и равновесия принимают вид:
m  m раз
∞
d
Θ
d
∇
u
λ ⋅ ⋅ ∇ 2 Θ − cε
− T0 ∑
⋅  ⋅ Β ( m ) = −ϑ ;
dt
dt
m =1
(23)


dv

( 1)
( 2)
( 3)
∇ ⋅ P − ∇ ⋅ P ⋅∇ P − + f = ρ
r ∈V ,
;
(24)
dt



 ;
( 1)
( 2)
( 2)
( 0)
(25)
r
∈ A,
n⋅ P − ∇ ⋅ P − − ∇A ⋅ n⋅ P − = Π ≡ σ




( 2)
( 3)
( 3)
( 1)
;
(26)
r
∈ A.
n⋅ P − ∇ ⋅ P − − ∇ A ⋅ n⋅ P − = Π ≡ Ν
Изменение плотности ρ учитываются законом сохранения массы (Седов Л. И. Механика сплошной
среды. – СПб.: Издательство "Лань", 2004. – Т. 1. – 528 с.)
∂ρ

+ ∇ ⋅ ( ρ v ) = 0 , ρ = const .
(27)
∂t
Для выполнения закона сохранения момента импульса необходима симметрия тензора
{
[
[
(
)
(
)
[
(
(
(
(
)
)
(
)
)]}
(
)]
)]
[ (
[ (
)]
)]
(
(
))
P = P ( 1) − ∇ ⋅ P ( 2 ) − ∇ ⋅ P ( 3 ) −  .
17
(28)
18
Для функций Φ ( 2 ) , Φ ( 3 ) допускаются представление рядами Тейлора, в которых сохранены лишь

слагаемые со степенью ( ∆u1k ) n не выше второй. В результате:
1 
T 1 
T 
w p = ∫  ∇ 2 Φ ( 2 )  ∆ u12
+ ∆ u12  ∇ 2 2 Φ ( 2 )  ∆u12
dV2 +

2! V 
2!

(
(
)
)

1   3
1 3 
T
T 
( 3)
( 3)
+ ∫   ∑ ∇ p Φ  ∆ u1 p + ∑ ∆ u1i  ∇ i ∇ j Φ ∆ u1 j  dV3  dV2 ,
3! V   p= 2
2! i , j = 2


d
n +1 d
∇ n =  = ( − 1)
 , n = 2, 3 .
dr
d l1n
С учетом этого:
∂P 0 ( n )
( n)
( n ,m )
0( n )
0( n )
0( n )
( n ,m )
0( n )
Β =−
= Pp
+ Pe
+ Pf
= Cp
,P
,C
,
∂T T =T
(
)
(
)
(29)
(30)
0

1 1
1 3  1
( 2 ) n
( 3 ) n
P
= ∫ ∇ 2Φ
l12 dV2 + ∑ ∫  ∫ ∇ p Φ l1 p dV3  dV2 ,
2! V n!
3! p= 2 V V n!

2
2nе µ 0 
5  π kT  
9 n k θ D γ xD x 3d x
 , pf =
1 + 
pe =
∫ e x − 1,
5  12  µ 0  
xD
0


 1 n

1
1 n 2 ( 2 ) m
1 2
( n ,m )
( 3 ) m
Cp
= ∫
l 1 ∇ 1 Φ l1 dV1 +
l
∇
∇
Φ
l
dV
∑
∫
p
p q
q
1  dV2 .
∫
2! V m! n!
3! p ,q =1 V V m! n!

(
r 0( n )
(
)
)
(
)
(
)
(31)
(32)
(33)
(
19
)
  
  
Классический тензор упругих постоянных C cl = ei e j ek em [ C cl ] ijrm , ( e1 , e 2 , e 3 ) – ортонормированный
(
)[
  
базис, i , j , k , m = 1,2,3 , тензору C ( 1,1) = ei e j ek em C ( 1,1)
[
Матрицы этих тензоров [ C cl ] ijrm и C ( 1,1)
[C ( ) ]
1,1
kλ
[C ( ) ]
1,1
ijkm
ijkm
]
ijkm
]
ijkm
, должен совпадать с тензором (33).
имеют вид:
λ
λ
λ + 2µ
 λ
λ + 2µ
λ

 λ
λ
λ + 2µ
=
0
0
 0
 0
0
0

0
0
 0
0 0 0
0 0 0

0 0 0
,
µ 0 0
0 µ 0

0 0 µ
η + χ + ε ( η + χ ) 3 + ε ( η + χ ) 3 + ε 0
( η + χ ) 3 + ε η + χ + ε ( η + χ ) 3 + ε 0

( η + χ ) 3 + ε ( η + χ ) 3 + ε η + χ + ε 0
=
( η + χ) 3 0
0
0
0
0
0
0
0 ( η + χ)

0
0
0
0
0
Здесь: λ, µ – коэффициенты Ламе.
(34)




.
0


3 0
( η + χ ) 3
0
0
0
0
0
0
(35)
20
2 ( 2)
∞
( l12 )
2π
16 π
4 d Φ
16π 2
d l12 ,
η=
J ,χ=
J 11 , ε =
J 12 , J = − ∫ l12
2
5
15
d l12
27
0
∞ ∞
J 11 = − ∫  ∫ l12
0 0
4
∂ 2 Φ ( 2 ) ( l12 )
∂ l12
2
2 ( 2)
∞ ∞
 2

 3
∂
Φ ( l12 )
3
d l12  l13 d l13 , J 12 = − ∫  ∫ l12
d l12  l13 d l13 ,
∂ l12 ∂ l12
0 0


(36)
(37)
2
2π
16 π 2
16 π 2
2
π
16
π
(38)
J+
J 11 +
J 12 = λ ,
J+
J 11 = µ .
15
45
27
15
45
Параметр ξ определен на основании описания нелинейного дисперсионного закона для плоских
продольных упругих волн высокой частоты с точки зрения предложенной модели и с точки зрения
физики твердого тела.
β≈
45 3
.
2 l0
(39)
8. При вычислении энергии адгезии используются потенциалы, описывающие взаимодействие частиц
разных материалов (α ) и (β ) (аналогия с Бараш Ю. С. Силы Ван-дер-Ваальса. – М.: Наука, 1988. – 344 с.):
2)
Φ (( αβ
) =
2)
H (( αβ
)
,
( l( α ) ij + l0( α ) ) ( l( β ) ji + l0( β ) )
( 3)
H
( αβγ) ( 3 cos ϑ i cos ϑ j cos ϑ k + 1)
( 3)
Φ ( αβγ) =
( l( α ) ij + l0( α ) ) 3 ( l( β ) jk + l0( β ) ) 3 ( l( γ ) kj + l0( γ ) ) 3 .
3
3
(40)
(41)
21
2
∞
 4π r03( α ,β ) 
3

(
)
2

 H
=
∫ α ( α ) ( i ω) α ( β ) ( i ω) d ω ,

 ( αβ)
3
π
0


(42)
3
∞
 4π r03( α ,β )  ( 3 )
3


 H
(43)
=
∫ α ( α ) ( i ω) α ( β ) ( i ω) α ( γ ) ( i ω) d ω .

 ( αβγ)
3
π
0


Здесь α , β , γ = 1,2 – номер материала.
Допускается также, что вступающие в адгезионный контакт тела B( 1) и B( 2 ) образовались
одновременно в момент времени t = 0 , а в момент t = 0 + 0 , непосредственно следующий за
моментом t = 0 , образовалось тело B( 1)  B( 2 ) . При этом деформации и изменения температуры
отсутствуют. Каждая частица dB( α ) ⊂ B( α ) , α = 1, 2 , мгновенно, наряду с взаимодействием с другими
частицами dB( α ) 1 ⊂ B( α ) и парами частиц dB( α ) 1 , dB( α ) 2 ⊂ B( α ) , начинает взаимодействовать с
частицами dB( β ) 1 ⊂ B( β ) , β = 1, 2 , β ≠ α , и парами частиц dB( α ) 1 ⊂ B( α ) , dB( β ) 2 ⊂ B( β ) , а также
dB( β ) 1 ⊂ B( β ) , dB( α ) 2 ⊂ B( α ) и dB( β ) 1 , dB( β ) 2 ⊂ B( β ) . Изменение их свободной энергии, происходящее в
этом процессе, определяется выражением:
1 2
1
( 2)
∆F( 1, 2 ) = ∫ Fa dA =
Φ
d
V
d
V
+
Φ (( α3 ),β , γ ) d V( α ) d V( β ) 1 d V( γ ) 2 −
∑
( α)
( β) 1
∫
∫
∫
∫
∫
( αβ)
2 ! α ,β =1V( α ) V( β )
3 ! V( α ) V( β ) V( γ )
Aa



(
(
2)
3)

− ∑  ∫ d V( α ) ∫ Φ ( α , α ) d V ( α ) 1 + ∫ d V( α ) 2 ∫ Φ ( α , α , α ) d V( α ) 1   .
V

α =1 V( α )
V
V
(
)
(
)
(
)
α
α
α



2
(44)
22
В случае, когда тела B( 1) и B( 2 ) полубесконечны, а поверхность Aa адгезионного контакта –
плоскость, распределение величины Fa вдоль нее можно считать равномерным.
σ, Дж/м2
5
4,5
4
3,5
3
2,5
2
1,5
1
0,5
0
Si
Al
- справочные данные [*];
W
Ni
V
Ti
Mo
- вычисленные значения [**];
Cu
Cr
Fe
- вычисленные значения [***]
*Свойства элементов. Физические свойства: Справочник. / Е. Б. Самсонов.– М.: Металлургия, 1976. – 600 с.
** Магомедов М. Н.
*** Фроленкова Л. Ю., Шоркин В. С.
23
Таблица 1 – Результаты расчета энергии адгезии для некоторых пар материалов
Состав
Энергия адгезии Fa
по литературным данным [*]
рассчитанная
системы
Cu-Al
3,10
2,978
Fe-Cu
4,00
3,763
Cr-Fe
5,00
5,286
*Вакилов А. Н., Мамонова М. В., Прудников В.В. Адгезия металлов и полупроводников в рамках
диэлектрического формализма // Физика твердого тела. – 1997. – т. 39. – № 6. – С. 964 – 967.
24
Спасибо за внимание
Download