МЕХАНИКА
УДК 539.3:536.2
Г. Н. Кувыркин И. С. Родикова
ОЦЕНКА ТЕРМОМЕХАНИЧЕСКИХ СВОЙСТВ МЕТАЛЛОВ И СПЛАВОВ В ЗОНЕ ФАЗОВОГО ПРЕВРАЩЕНИЯ1
Предложена схема определения термомеханических свойств металлов и сплавов в зоне фазового превращения на основе двойственной вариационной формулировки задачи.
В работе [1] приведены основные соотношения, описывающие фазовые превращения металлов и сплавов, сопровождающиеся изменением характерного размера кристаллической решетки или изменением типа ее симметрии, а также соотношения, позволяющие оценить границы возможного изменения модулей всестороннего сжатия и сдвига в процессе фазового превращения. Данная статья является продолжением работы [1] и посвящена оценке других свойств металлов и сплавов в зоне фазового превращения.
Для оценки теплопроводности металлов и сплавов в зоне фазового превращения положим, что в рассматриваемом гетерогенном материале поле температуры таково, что dTa/dxi = dT/dxi = const, i = 1, 2, 3, a = 1, 2, где T,Ta — абсолютная температура эквивалентной гомогенной среды и a-го компонента материала. Функционалы, аналогичные функционалам энергии в механике деформируемого твердого тела для неоднородной и однородной сплошной среды, при отсутствии источников объемного тепловыделения имеют вид
1 Г а (Т) dT dT f
J«[T«] = 2/ Xa\j dXf dXadV - ' qaTadS ,
V S
1 [T]=2/
jdTdTdV - IqTdS,
(1)
dxj dx
v S
a
(T)
\ \(т)
где л ^ , л - — компоненты тензоров теплопроводности а-го компонента материала и эквивалентной гомогенной среды; да,д — заданная на граничной поверхности Б тела объемом V плотность теплового потока для неоднородной и однородной среды соответственно; XI ,х2
хРабота выполнена при поддержке программы "Университеты России" (проект УР 03.01.139) и РФФИ (проект № 05-01-00596)
— объемные доли мартенсита и аустенита в сплаве в зоне фазового превращения, х1 + х2 = 1. Положим также, что
2
^ J qaTadS = J qTdS.
a=l s S
2
Тогда из условий минимума Ia[Ta] и I[T] следует неравенство
La [
а=1
а=1
Если в гетерогенном теле задано поле вектора плотности теплового
а
потока такое, что qi = qi = const, а = 1,2, то двойственные по отношению к выражениям (1) функционалы имеют вид
2
Ja[q] = -1 / Xah]qtdV - T**'qlnldS,
у я (2) 3 ^=- 2/кгз ^ - / т**стлв,
V я
а
где к^, к^ — компоненты тензоров термического сопротивления;
Т**, т** — заданные на граничной поверхности Б значения темпе-
2
EOi
Ja [q] и J [q] следуют неравенства
а=1
2 (т) (т) kij > Xakij, а так как (kj)-1 = Aj ), (kj)-1 = Aj , то оконча-
=1
тельно
1 -1 21
аЛ
2 (т) 2 (т) -1
EXaAij > j XalKj) . (3)
)
а=1 \а=1
Для изотропного гомогенного материала и гомогенных компонентов из формулы (3) следуют неравенства
£ > > (£ А^] 1 , (4)
а=1 \а=1 Аа )
а(т) (т)
в которых А^ = -а .
В неравенствах (3) и (4) верхние границы аналогичны соответствующим оценкам по Фойгту, а нижние — по Рейссу в механике. Оценки (3) и (4) компонентов тензора теплопроводности эквивалентной гомогенной среды далеки от свойств реального материала.
(т)
Для получения более узких границ значений А]) рассмотрим двухкомпонентный материал. Для средних по объему представительного элемента гетерогенного материала, а также для средних по объему компонентов материала значений дг, дТа/дхг справедливы соотношения
а а (Т) дТа 12 дТ дТ, дТ2
а = -Аг] 9X7' Чг = + Х24' дх = Х1 дХ" + Х2дХ7' (5)
где а = 1, 2; дг, дТ/дхг — средние значения компонентов вектора плотности теплового потока и градиента температуры в объеме представительного элемента. Соотношения (5) позволяют установить связи
dTa а дТ 1 2 ^ 1 2 _
ij =
(6)
qt = At]qj , дХГ = Dij д^ ' XlAij + X2Aij = Sij ' Xl+ X2Dij = Sij '
(Т) 1(Т) 1 2(Т) 2 11 2 2 А(т = Х1Агк + Х2Агк , кг] = Х1кгкАк] + Х2кгкАк] ,
— символ Кронекера. Если А(Т) и кг]- известны, то тензоры с
а(т) а
компонентами Аг]- и кг] определяются из соотношений
а (Т) (Т )\ а (т) V (Т)
Ха ( Агк - Агк ) &к] = Аг] - Аг] ,
4 (7)
/а V \ а V
Ха ( кгк - кгП Ак] = кг] - кг] , а, V =1, 2, а = V .
Умножив первое из уравнений (7) на дТ/дх], а второе — на Ч], с учетом равенств (6) получим
Ха ( Аг, - Аг, ) дТ^ = ГА(Р - ■ -
г] г] дх] г] г] дх]
] ] (8) /а V А а / v А
Ха ( кг] — кг] ) Ч] = ( кг] — кг] ) .
Если компоненты гетерогенного двухкомпонентного материала и эквивалентного гомогенного материала изотропны, то из соотношений (8) следуют равенства
дТа = А(Т) - лУТ ) дТ а = 1/А(Т) - 1/лУТ ) (9)
Ха дхг лат) - аут) дхг, ХаЧг 1/лат) - 1/лУт) Чг, (9)
где учтено ка = 1 /а7Т), к = 1 /А(т).
Аналог потенциальной энергии деформации рассматриваемого гетерогенного материала объемом V
а
U [T]=-1 i q IdV=-1
V
X1/ 1 ЁdV+X2
'2 dT2 dV+ q ^ öX2dV+
+
V
dT ÖT2
ÖXi ÖXi
dV I = U[T] + U2[T] + Ui2[T],
где и1[Т], и2[Т], и12[Т] — средние значения аналогов потенциальной энергии деформации компонентов материала и потенциальной энергии их взаимодействия. Третье слагаемое в правой части последнего равенства можно представить в виде
Ui2[T] = 2ч ^' - ^О (л(Т' - л?') 11dV,
У
(T K(T)
где Ax =
лг )л;
2
xirf) (лр - л2")2'
л>т > = х i л1Т) + Х2л2Т),
(л<Т>)- 1 = хi СлГМ + Х2 ( л
>
-1
л
(T), V,
(T)
л^ — оценки тепло-
проводности по неравенствам (4).
В том случае, когда эффективное значение А(т) фиксировано, любое изменение ^ и дТ/дх^ приводит к соответствующему изменению средних значений ^ дТа/дхг в компонентах материала. С другой стороны, можно задать вариации значений ^ и дТа/дХг, не изменяя значения переменных qi и дТ/дх^, что приводит к изменению эффективного значения А(т).
Рассуждая далее аналогично тому, как это было сделано в работе [1] при получении оценок модуля всестороннего сжатия, можно получить неравенства
л^ + л^) > л(т) >
2 л (т )л(т)
V + ля
л
(T)
(T)
(10)
Отметим, что неравенства (10) эквивалентны первой группе неравенств (24) [1]. "Вилка" Хашина-Штрикмана для определения границ А(т) аналогична соответствующим неравенствам для модуля всестороннего сжатия, т.е.
А1Т) + Х2(-2 --)А1 )(т) > А(т) > -2т)+
1 1 + Х1«1(А2Т) - А1т)) 2
+
Х1(л1Т) - л2Т))
1+ Х2«2(л1Т) - л2Т))
л1Т) > л.
(т ).
2
где аа = 1/а7Т) , а = 1, 2.
Очевидно, что оценка теплопроводности по неравенствам (10) дает более узкий диапазон возможных значений А(т) по сравнению с "вилкой" Хашина-Штрикмана.
Для оценки величины температурного коэффициента линейного расширения а(т) и удельной массовой теплоемкости при постоянной деформации се и напряжениях двухкомпонентной смеси воспользуемся выражениями для средних по объему V представительного элемента значений объемной плотности свободной энергии рА и термодинамического потенциала Гиббса рФ [2]:
V1 урА^ = 2К ^^^ Хаак^ + Ха] Ха] -
3tf,af- ^Ö2,
a=1
(11)
V /P$dV =18K: ^ + 4^ (z ** j ^ IV \a=1 / \a=1 / \a=1
1 r*
- 1 «f^ExA* - 2TrÖ2,
a=1 0
где K*, p*, a*T), c*, c* = c* + 9ToK* (a*T)) /Р — эффективные свойства смеси, находящейся в однородном поле температуры в = T — T0, To = const, | T — To |/To ^ 1.
Следуя работе [2], сначала рассмотрим на границе S объема V граничное условие
a0, n, = p0 = const. (12)
Полагая, что массовые силы отсутствуют, имеем
J [a ,в]= I [и, в], (13)
где J [a , в] = — J pФdV — J a,«¿dS — дополнительная работа тела, I [и, в] = J pAdV — J p0MidS — полная потенциальная энергия тела.
V
Здесь и Би — части поверхности Б, на которых заданы компоненты векторов плотности поверхностных сил р° и перемещений йг. При выполнении граничных условий (12) из равенств (11) следуют соотношения
1 1 2
^I [и , 0] = -2 (а?;. + 3К*а*Т ^) £ Х^ - ^,
а=1 °
V3[а , 0] = 1 (£ Х^; + а*^;) а?; + 2|02,
которые записаны с учетом закона Дюамеля-Неймана для эффективных величин
2 2 2 Е Х^а; = к* Е Х^с^Ьг; + Е Хае; - 3К*а*т)0Ь; ,
а=1 а=1 а=1
2 1 Л а , 1 2
Е Xa^ij = Е Xa^fefeÜj + Е ХаS j + «f^ü
УК* ' ^
а=1 а=1 а=1
и очевидных равенств
2
/ a0dV = a0n,U"dS = a0 E Xa£j
/л a a=1
V я
Если на границе Б элементарного объема V задано граничное условие
"г^^ г;
U = u0 = Xj , = const, (14)
то, как и ранее, при отсутствии массовых сил
2
V1 [u,= 24 ( £ Xa^ij - BKaf^j ) - 02 ,
a=1
-13 [а , 0] = 1 (с° - )0Ьг^ Е + ^02 .
а=1
Если в выражениях 3 [а ,0] и I [и, 0] истинные поля деформации и напряжений при заданных граничных условиях (12) и (14) заменить на произвольные, не обязательно удовлетворяющие уравнениям равновесия в объеме V, но обязательно граничным условиям, то должны выполняться неравенства
-3[а ,0] < -3 [а ,0] < I [и, 0],
(15)
-3 [а ,0] < I [и, 0] < I [и, 0],
1 , 1 , п лгл„,- ) a С
где 7 [а ,9], I [и, 9] — дополнительная работа тела и полная потенциальная энергия тела, вычисленные на произвольных полях деформации и напряжений. Неравенства (15) записаны с учетом соотношения (13). Воспользовавшись граничными условиями (12), запишем первое неравенство (15) в явном виде:
1 о 1 в
а2° + Т~ Sj Sj + a(T}ва°° + c) dV <
V] кк ] ] кк " 2Т0
V
1 1 92
* 1 (а0 )2 1 „0 „0 а(т)9а0 с; 9 *
* - 18ККк) - „0] - а 9а0к - 2Т0 *
* V1 У Кекк +Рег]ег] -3Ка(Т)9^кк-Се^- V1 %^ .
V V
(16)
Вторая цепочка неравенств (15) после подстановки граничных условий (14), примет вид
1 1 1 92 1
-V.) ( 18Ка2к + „г] + а(Т)9^кк + — ) ^+V?4 I ^ *
V V
< 2K(£°°°)2 + 4 - 3Kaf)e£°kk - c:^ <
1 1 92
* V1 ( 2Ке2к + ^ег]ег] - 3Ка(Т)9^кк - Се— ) ^ .
V
Примем в неравенствах (16) аг] = а]. Тогда, полагая неизвестные деформации однородными, получим
-18 (К) <"°»)2 - 4 (£) 4 -<а<Т'«к - (С)^ *
1 1 92
* -7^«)2 „0]„0] - а ;Т)9а0к - с^ * (17)
18K*X °k/ * j j * °° 17 2Т°
°
1в
2
* 2)(^кк)2 + Мег]ег] - 3(Ка(Т))9^кк - (се)— - а]^ ,
где угловые скобки обозначают операцию статистического осреднения.
Величину ег] выберем такой, чтобы верхняя граница в выражении (17) была минимальной. Это приводит к условию
(К)£кк¿г] + 2(р)ег] - 3(Ка(Т))9^ - а0] = 0.
1
С помощью этого условия правую часть неравенств (17) можно записать иначе, что приводит их к другой форме:
18« K> - K) <*>2+4«£ > - £) sjsj+
+ (<а(Т») - а<т>)0<4 + (<с») - с;) ¿т > 0. (18)
118 ( K)(* )2+4 (i^) sj +(a2--<K<f)) ^+
+ (^ + 2 (<*)) «2 > 0. (19)
Переход от <ce) к <c;) в неравенстве (19) произведен с помощью зависимости <с;) - <се) = 9T0 {K(a(T))2) . Примем теперь, что компоненты тензора напряжений aj = a0£j, a0 = const. Тогда неравенства (18) и (19) можно записать в виде
18 ((K> - K) x2 + «а(Т>( - af') x > 0, (20)
1 ( 1 1 ) x2 + - iKO^ ,x +
18 V К* (К)у у * (К)
+ (^ + 2 ((^Г) - > (2,)
где х = а°/0. Поскольку неравенства (20) и (21) должны выполняться при любых значениях х, заключаем, что
1 \ 1 11 -к) - к > 0, К-(К) * 0;
(с;)-с; > 0+ 2 ((К К ))2) - в2) > 0.
Первая система неравенств дает оценку модуля всестороннего сжатия по Фойгту и Рейссу. Вторая система дает оценку удельной массовой теплоемкости материала с;:
<С) + 9To(- (*(a(TЧ2)) < c: < <С>-
(22)
Необходимым и достаточным условием того, чтобы квадратичная форма у = ах2 + Ьх + с при а > 0 была неотрицательной, является требование неположительности ее дискриминанта Д = Ь2 - 4ас.
Поэтому из неравенств (20) и (21) следует
а
(г )\
а
*т ))2
а
(T) _
Ка(т )N
1
К
2 9
1
К
((c.) - c.)
2T°
< 0,
(23)
1
К
(k ),
- (c.) 2T°
+ 2 (<K(а'т')2> - ^ <
0.
(24)
Система неравенств (23), (24) содержит в качестве неизвестных величин а !Т) и с*. Учитывая, что эффективный модуль всестороннего сжатия может быть рассчитан существенно точнее по вилке (24) из работы [1], будем в дальнейшем считать, что его точное значение известно. Тогда система неравенств (23), (24) будет определять область допустимых значений эффективных величин а !Т) и с*.
(т)
Неравенства (23), (24) ограничивают возможные значения а( ) и с* внутренними областями парабол, уравнения которых получаются из неравенств (23), (24) заменой знака неравенства на равенство. Обе параболы имеют оси симметрии, параллельные оси с*. Очевидно, что вилка для с* определяется максимальным и минимальным значениями на границе области пересечения парабол. Эта вилка совпадает с неравенствами (22).
Для нахождения вилки величины а !Т) заменим в выражениях (23) и (24) неравенства на равенства и исключим с*. Это приводит к уравнению
а
(г )2
1 + £ - 2а*т) c + b(а(т)) + c2 - bc + Ь <а(т)>2 = 0 , (25)
где b=2(( Ю- К
9
c =2
[К(а(т-
b = 20.-J-
, П к* (К)
Ка(т )>2
c=
Ка(т)
w
(K)
Решение уравнения (25) определяет границы для эффективных значений температурного коэффициента линейного расширения:
,<тc + b (а(Т')
±
D
D=
1 + b (а(т)) 1 + b (а(т))
£ (c(b + b) - (<а(т)> - c)2
(26)
1/2
2
1
c
Так как диапазон изменения температуры, в котором происходят фазовые превращения, достаточно узок, то в зависимости от направления этого процесса или 0 = Т- Мв, или 0 = Т- А и | 0 | / Мв ^ 1, | 0 | / А ^ 1, где А, Мв — температуры начала превращения мартенсита в аустенит и аустенита в мартенсит. В этом случае в соотношении (7) [1] положим
В (Т,Х) = -1Т02 + ^0Х + В1(Х), В1(Х) = -17еХ2 (27)
и примем
с!Т) = а!Т)0, с(Х) = а(Х)Х, (28)
где «гр , «гх) — компоненты тензоров температурных коэффициентов линейного расширения и коэффициентов фазовой деформации. С учетом соотношений (27) и (28) выражение (7) [1] примет следующий вид:
РА(сЫ , Т , Х) = 1 Сг;ДО сЫсг; - )сг;0 - ^^г.;Х-
-1 Р^02 + р^0Х - 2Р^Х2 , (29)
2 Т° Т° 2
где в(Т) = Одав(х) = Ода«?, Т° = М5 или Т° =
Так как диссипативная функция Ьр ^ 0 всегда, то с точностью до линейных слагаемых по компонентам тензора деформации сг; из равенства (29) имеем
г дА / ше \
ьр = -Р^Х = Р 0 + 7еХ Х > 0 .
дХ Т°
Последнее неравенство выполняется всегда, потому что при нагревании материала происходит превращение мартенсита в аустенит и объемная доля мартенсита Х убывает, Х < 0. При обратном превращении с понижением температуры объемная доля мартенсита Х возрастает, Х > 0.
Рассмотрим представительный объем V материала. Пусть этот объем находится в однородном поле температуры. Тогда при сг; = 0 и г = 0 из уравнения (12) [1] следует
Рсе0 - ртеХ + р(те0/Т° - 7еХ)Х = фф, (30)
где = q(t)Б/V, Б — площадь поверхности представительного
элемента, q(t) — плотность подводимого к этой поверхности теплового потока, — количество теплоты, полученное единицей объема представительного элемента в единицу времени. Найдем верхнюю и нижнюю границы удельной массовой конфигурационной теплоемкости материала ше и коэффициента 7е из графика зависимости количества теплоты , переданного единице массы материала при его
Рис. 1. Зависимость количества теплоты, переданного 1 кг материала никелида титана при его нагреве на 1 И, от температуры
нагревании на 1К [ 3 ] (рис. 1), по соотношению (30) с использованием неравенств (22) и равенства ca = ce + 9T0K(a(T))2 :
pcfown(x)AT - S(AT, Ax,T,x)
paX
< m£(x) <
^ pcUP(x)AT - S(AT, Ax,T,x) " PAx :
(31)
mdoWn(x)e/(Tox) < 7e(x) < mUp(x)0/(Tox),
(32)
где cfown,
„up _
нижняя и верхняя границы удельной массовой те-
down
е ,
m
up _
нижняя и верхняя
плоемкости материала (см. (22)), т\ границы удельной массовой конфигурационной теплоемкости материала, Б(АТ, Ах, Т, х) — количество теплоты, переданное единице объема материала для изменения температуры тела на АТ и изменения доли мартенсита в материале на Ах; Т, х — исходное значение температуры и объемной доли мартенсита. Соотношения (31) и (32) получены с учетом предположения, что диссипация энергии отсутствует. Предположим, что объемная доля мартенсита линейно зависит от температуры материала и Б(АТ, Ах , Т, х) будем вычислять как
площадь трапеции Б (А Т, Ах , Т, х) = - (Ят+а.т ) А Т. Таким обра-
2
зом, в случае обратного фазового превращения, выражения (31), (32) можно представить в следующем виде:
е
c
down
(X) - 1QT + QT+ДТ') (As - A,) « m,(x) i
i I c?'(X) - 1 Qt + QT+ЛТ^ (A, - A,),
m
down
(X)< 7e(x) < mup(x). (34)
f
A,
TA
f
A,
Численный анализ границ величин се , ше , 7е по неравенствам (22), (33), (34) проводился для никелида титана (49% И) [3,4] (рис. 2-4).
'sз
550
500
450
400
кг К
350
\
\ /
ч /
0
0,2
0,4
0,6
0.8
Рис. 2. Зависимость нижней (сплошная линия) и верхней (штриховая линия) оценок теплоемкости никелида титана от объемной доли мартенсита
т.
■е>
кг
1,5-10
1.0-10"
0,5-10
А
\ н
У V
0
0,2 0,4
0,6
0.8
Рис. 3. Зависимость нижней (сплошная линия) и верхней (штриховая линия) оценок конфигурационной теплоемкости никелида титана от объемной доли мартенсита
7 s
Дж
кг
8 ■ LO
б-LO
4-10"
2-LO
W
1
// №
0
0.2
0,4
0,6
0,8
Рис. 4. Зависимость нижней (сплошная линия) и верхней (штриховая линия) величины 7е никелида титана от объемной доли мартенсита
а.
В том случае, если охлаждение (нагревание) тела проходит при = 0, то деформация тела будет состоять из температурной и фа-
зовой:
£(Т,Х) = а(Т),
ij ij
в ± - Хо),
(35)
Выбор знака перед вторым слагаемым в правой части равенства (35) определяется свойствами кристаллической структуры и составом рассматриваемого материала. Так, при возрастании температуры от значения Л3 длина I цилиндрических образцов из одних (изотропных) материалов в зоне фазового превращения может увеличиваться, а из других (также изотропных) — уменьшаться [3-5]. При понижении температуры от значения Ы3 до окончания обратного превращения эффект противоположен.
На рис. 5 и 6 приведены расчетные границы температурного коэффициента линейного расширения а(т) и коэффициента фазовой деформации а(х) по неравенствам (25) и (35) Рис. 5. Зависимость коэффициента температурно-для никелида титана го линейного расширения никелида титана от
объемной доли мартенсита
€,05 0,04 0,03 0,02 0,01
0
0,2 0,4 0,6 0,8 /.
Рис. 6. Зависимость коэффициента фазового расширения никелида титана от объемной доли мартенсита
(53,8% Ti) по экспериментальным данным из работы [3]. Отметим, что верхние и нижние границы коэффициентов а(Т) и а(х) на соответствующих рисунках практически неразличимы.
В заключение можно отметить, что предложенный подход к определению термомеханических свойств металлов и сплавов в зоне фазового превращения позволяет получить оценки необходимых для расчетов значений параметров материалов.
СПИСОК ЛИТЕРАТУРЫ
1.Кувыркин Г.Н., Родикова И. С. Термомеханическая модель поведения металлов и сплавов в зоне фазового превращения // Вестник МГТУ. Сер. "Естественные науки". - 2006. - № 1. - С. 65-76.
2. Шермегор Т. Д. Теория упругости микронеоднородных сред. - М.: Наука, 1977. -400 с.
3. Сплавы с эффектом памяти формы / Ооцука К., Симидзу К., Судзуки Ю. и др.: Пер. с японск. - М.: Металлургия, 1990. - 224 с.
4. Лихачев В. А., Кузьмин С. Л., Каменцева З. П. Эффект памяти формы. - Л.: Изд-во ЛГУ, 1987. - 216 с.
5. Абдрахманов С. Деформация материалов с памятью формы при термосиловом воздействии. - Бишкек: Илим, 1991. - 115 с.
Статья поступила в редакцию 19.09.2005 Георгий Николаевич Кувыркин родился в 1946 г, окончил МВТУ им. Н.Э. Баумана в 1970 г. Д-р. техн. наук, профессор кафедры "Прикладная математика" МГТУ им. Н.Э Баумана. Автор более 100 научных работ в области термомеханики и структурной механики конструкционных элементов.
G. N. Kuvyrkin (b. 1946) graduated from the Bauman Moscow Higher Technical School in 1970. D. Sc. (Eng), professor of "Applied Mathematics" department of the Bauman Moscow State Technical University. Author of more 100 publications in the field of thermal mechanics and thermal structural mechanics of construction elements.