Тройные диаграммы
.pdfвидуса и солидуса), или после гомогенизационного отжига неоднородного литого сплава.
Нами выполнены расчеты возможного переохлаждения, соответст вующего неравновесным условиям кристаллизации в тройных сплавах на ос нове палладия. Для разных составов это переохлаждение составляет 2060 ºС, что соответствует нашему предположению о возможных неравновесных ус ловиях кристаллизации сплавов, данные по которым приводятся в литературе
(см. таблицу 6.3.1).
Лекция 17
План лекции
Общие представления о методах расчета диаграмм состояния тройных и многокомпонентных систем с использованием термодинамических моделей растворов.
Применение термодинамической модели идеальных растворов замещения.
Применение термодинамической модели регулярных растворов замещения.
6.4. Методы расчета диаграмм состояния тройных и многокомпонентных систем с использованием термодинамических моделей растворов
6.4.1. Общие представления о методах расчета диаграмм состояния
Расчетные методы для анализа диаграмм состояния применяют доста точно давно1,2. Только появление современной компьютерной техники и раз личных автоматизированных систем, например, MathCad, сделало реальным выполнение чрезвычайно трудоемких вычислений и представление результа тов расчетов в графической форме.
Одним из методов, позволяющим описать фазовые равновесия в трой ных и многокомпонентных системах является метод, использующий понятия:
__________________
1Кауфман Л., Бернстейн Х. Расчет диаграмм состояния с помощью ЭВМ. М.,1972 2Расчет фазовых равновесий в многокомп. системах. Под.ред.В.Н.Романенко . – М., 1987
130
«термодинамические модели идеальных, регулярных, субрегулярных раство ров», «квазихимическое приближение». В этом случае применяют статисти ческие теории, оперирующие моделями межатомного взаимодействия в рас творе1,2. Энергию раствора рассматривают как сумму энергий взаимодейст вия между парами атомов, находящихся на ближайших расстояниях. При этом полагают, что энергии взаимодействия между однородными и разно родными атомами известны и имеют постоянные значения независимо от ме стоположения связи или что их можно определить путем сопоставления с экспериментом. При этом необходимо рассматривать отдельно теорию рас творов замещения и внедрения.
В любом случае при таком подходе условием равновесия фаз при за данной температуре в тройных и многокомпонентных системах является ра венство химических потенциалов компонентов в сосуществующих фазах.
Рассмотрим применение наиболее разработанных моделей растворов замещения к описанию фазовых равновесий в тройных и многокомпонент ных системах.
6.4.2. Применение термодинамической модели идеальных растворов замещения
Идеальный раствор проще всего определить как раствор, теплота сме шения которого равна нулю, а энтропия смешения равна конфигурационной энтропии хаотического расположения атомов. Величина энергии Гиббса иде ального nкомпонентного раствора, имеющего структуру φ, равна
G |
ϕ |
= |
n |
ϕ |
0ϕ |
+ RT |
|
å |
x G |
|
|||
|
|
|
i = 1 |
i |
i |
|
n |
ϕ |
ϕ |
(6.4.1) |
å |
(x |
ln x ) , |
|
i = 1 |
i |
i |
|
0φ
где xi атомная (мольная) доля iго компонента в растворе; Gi начальное
значение энергии Гиббса чистого iго компонента φмодификации; n число компонентов в растворе; R универсальная газовая постоянная; T абсолютная температура.
Эти растворы иногда называют совершенными растворами. Если энер гия Гиббса хотя бы одного из компонентов не удовлетворяет уравнению (6.4.1), то такой раствор не будет идеальным.
__________________
1Кожеуров В.А. Статистическая термодинамика. М., 1975 2Глазов В.М., Павлова Л.М. Химическая термодинамика и фазовые равновесия. – М., 1988
131
Для трехкомпонентного раствора, например α, системы АВС выра жение (6.4.1) имеет вид:
Gα = xAαGA0α + xBαGB0α + xCαGC0α + RTxαA ln xAα + RTxαB ln xBα + RTxCα ln xCα ,
(6.4.2)
где xαA + xαB + xCα =1.
Аналогичные выражения можно записать для четверного, пятерного и более компонентного идеального раствора замещения.
Для химического потенциала (μi) любого компонента идеального рас твора замещения справедливо соотношение
μ |
i |
= μ |
0 |
+ RTln x , |
(6.4.3) |
|
|
i |
i |
|
0
где μi — химический потенциал чистого iго компонента;
xi — молярная или атомная доля iго компонента.
Для примера рассмотрим двухфазное равновесие в трехкомпонентной системе АВС с неограниченной растворимостью компонентов в твердом и жидком состоянии (рисунок 6.4.1). Условию равновесия жидкой фазы и α твердого раствора при заданной температуре Т1 соответствует равенство хи мических потенциалов (парциальных энергий Гиббса) компонентов А, В и С в сосуществующих фазах, то есть
μАж |
= μАα ü |
|
||
μВж |
|
ï |
|
|
= μВα ý |
(6.4.4) |
|||
ж |
= μ |
α ï |
||
|
||||
μС |
С þ |
|
Рисунок 6.4.1 – Иллюстрация к описанию двухфазного равновесия в трехкомпонент ной системе АВС
132
Подставив выражения для химических потенциалов компонентов А, В и С в (6.4.4) и выполнив несложные преобразования, получим
|
|
|
ì |
|
|
|
|
xж |
|
|
|
|
|
|
|
|
|
|
|
|
ïDGα ® ж |
+ RTln |
А |
|
|
= 0 |
|
|
|
|
|
||||
|
|
|
α |
|
|
|
|
|
|
|
|||||||
|
|
|
ï |
A |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ï |
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
А |
|
|
|
|
|
|
|
|||
|
|
|
ï |
|
|
|
|
xж |
|
|
|
|
|
|
|
|
|
|
|
|
ï |
|
|
|
|
|
|
|
|
|
|
|
(6.4.5) |
||
|
|
|
íDGα ® ж |
+ RTln |
В |
|
|
= 0 , |
|
|
|
||||||
|
|
|
xα |
|
|
|
|
|
|||||||||
|
|
|
ï |
В |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ï |
|
|
|
|
В |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
ж |
|
|
|
|
|
|
|
|
|
|
|
|
|
ï |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ïDGα ® ж |
+ RTln |
хС |
|
|
= 0 |
|
|
|
|
|
||||
|
|
|
α |
|
|
|
|
|
|
|
|||||||
|
|
|
ï |
С |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
î |
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
С |
|
|
|
|
|
|
|
|
|
|
æ |
μ |
0ж |
- μ |
0α ö |
= DG |
α®ж |
|
æ |
0ж |
- μ |
0α ö |
= DG |
α®ж |
и |
|||
где разности ç |
А |
÷ |
A |
, ç μ |
В |
|
В |
÷ |
В |
||||||||
è |
|
|
А ø |
|
|
è |
|
|
ø |
|
|
æμ0ж - μ0α ö = DGα®ж соответствуют изменению энергии Гиббса при темпе
ç С С ÷ С è ø
ратуре Т чистых компонентов А, В и С модификации α при плавлении. Эти величины получили название «параметры стабильности фаз».
Таким образом, для получения изотермического сечения системы АВ С (рисунок 6.4.1) расчетным путем с использованием термодинамической модели идеальных растворов замещения необходимо задаться температурой Т1, знать тип кристаллической решетки раствора и значения параметров ста бильности фаз.
Нами* была создана электронная база данных “Periodic.mdb” в формате программы Microsoft Access, содержащая в себе справочную информацию по 109 элементам Периодической системы Менделеева. Вид окна формы «Elements» электронной базы данных приведен на рисунке 6.4.2.
Готовый к распечатке пример отчета для элемента молибден показан на рисунке 6.4.3.
Кроме того, создана электронная книга для определения величины мольной постоянной энергии стабильности фазы и температуры плавления элемента.
Ларри Кауфманом были получены функциональные зависимости меж ду разностями энергий Гиббса различных фаз и температурой для многих металлов. Анализ выражений позволил ввести понятие «мольная постоянная энергии стабильности твердой фазы», ее размерность Дж×(моль–1∙К–1). Представление энергий Гиббса фаз с использованием этой величины являет ся более удобным. Полученные данные были включены в электронную кни гу. Фрагмент справочной таблицы представлен на рисунке 6.4.5.
_________________
*В работе принимал участие Коновалов П.В.
133
Рисунок 6.4.2 Вид окна электронной базы данных «Elements»
Рисунок 6.4.3 – Вид отчета «Elements» электронной базы данных
134
Элемент |
|
Мольная постоянная стабиль |
|
Температура |
|
|
|
|
|
||||
|
|
ности твердой фазы (Gkn) |
|
плавления (Tn) |
|
|
|
|
кал∙моль–1∙К–1 |
Дж∙моль–1∙К–1 |
|
К |
|
Al |
|
2.75 |
11.506 |
|
933 |
|
Be |
|
1.6 |
6.694 |
|
1562 |
|
Co |
|
2.32 |
9.707 |
|
1768 |
|
Cr |
|
2 |
8.368 |
|
2136 |
|
Cu |
|
2.3 |
9.623 |
|
1357 |
|
Hf |
|
2 |
8.368 |
|
2504 |
|
Рисунок 6.4.4 — Фрагмент справочной таблицы из электронной книги
При заданной температуре три уравнения (6.4.5) содержат четыре не известные (если учесть, что xαA + xαB + xαC =1 и xжA + xBж + xCж =1), то есть их
система является ненасыщенной. Для решения системы уравнений было вве дено четвертое «насыщающее» уравнение.
Полученная система уравнений позволяет находить температуры соли дуса и ликвидуса для тройного сплава с заданной концентрацией компонен тов. Она также дает возможность отыскивать концентрации компонентов в твердой и жидкой фазах для заданных состава сплава и температуры (этих данных достаточно для построения серии конод).
.
6.4.3.Применение модели регулярных растворов замещения
Втеории регулярных растворов предполагают, что энтропия такая же, как и в идеальном растворе, а отклонение от идеальности обусловлено не равной нулю энтальпией смешения. Регулярные растворы образуются с по глощением или выделением тепла
Энергия Гиббса однофазного регулярного раствора, имеющего струк туру φ и n сортов атомов равна:
G |
ϕ |
n |
ϕ |
0ϕ |
n |
ϕ |
|
= å xi |
Gi |
+ RT å (xi |
|||
|
|
i = |
1 |
|
i = |
1 |
ϕ |
n-1 |
n |
ϕ ϕ |
ϕ |
, |
(6.4.6) |
|
ln x |
)+ å |
å |
x |
x E |
|
||
i |
i=1 j=i+1 |
i |
j ij |
|
|
||
|
|
|
n
где åxiϕ =1, xiϕ атомная доля iго компонента в растворе со структурой φ;
i=1
Т абсолютная температура; Gi0ϕ энергия Гиббса φмодификации чистого i
го компонента; Eijϕ энергия взаимообмена или энергия взаимодействия ато мов iго и j–го компонентов в φфазе.
135
Для тройного раствора, например α, системы АВС выражение (6.4.6) имеет вид:
Gα = x |
G0α + x |
G0α + x |
G0α + RTx |
ln x |
+ RTx |
ln x |
+ RTx |
ln x |
+ |
|||
A |
A |
B |
B |
C |
C |
A |
A |
B |
B |
C |
C |
(6.4.7) |
+ x x Eα + x x Eα + x x Eα . |
|
|
|
|
|
|
||||||
A B |
AB |
A C |
AC |
B C |
BC |
|
|
|
|
|
|
В трехкомпонентной системе АВС (рисунок 6.4.1) равновесие жид кой и αфазы при температуре T1 можно описать уравнениями:
ì |
|
α ® ж |
|
xж |
æ |
|
ж ö |
2 |
ж |
æ |
ж ö |
2 |
ж |
æ |
α ö2 |
|
α |
||
ï |
|
|
А |
|
|
||||||||||||||
ï |
DG |
|
+ RTln |
|
+ ç x |
|
÷ |
E |
|
+ ç x |
|
÷ |
E |
|
- ç x ÷ |
E |
|
||
А |
α |
|
|
|
|
АВ |
|||||||||||||
|
|
è |
В |
ø |
|
АВ è |
С |
ø |
|
АС è |
В ø |
|
|||||||
ï |
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ï |
|
|
|
А |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ï + xж xж (Eж + Eж - Eж ) - xα xα (Eα + Eα - Eα ) = 0
ï |
В С АВ |
|
АС |
|
ВС |
|
В С |
АВ |
|
АС |
ВС |
|
|||||
ï |
|
|
xж |
|
|
|
|
2 |
|
|
|
2 |
|
|
|
2 |
|
ï |
α ® ж |
|
|
|
æ |
ж ö |
ж |
æ |
ж ö |
ж |
|
æ α ö |
α |
||||
|
В |
|
|
|
|||||||||||||
ïDG |
|
+ RTln |
|
|
+ |
ç x |
÷ |
E |
|
+ ç x |
÷ |
E |
|
- ç x ÷ |
E |
|
|
|
|
|
|
|
|
||||||||||||
ï |
В |
|
xα è А ø |
|
АВ è С ø |
|
ВС è А ø |
|
АВ |
||||||||
í |
|
|
|
|
|||||||||||||
ï |
|
|
В |
|
|
|
|
|
|
|
|
|
|
|
|
ï + xж xж (Eж + Eж - Eж ) - xα xα (Eα + Eα - Eα ) = 0
ï |
А С АВ |
ВС |
|
|
АС |
|
А С |
|
АВ |
|
ВС |
АС |
|
|
|||||
ï |
|
|
xж |
|
|
|
|
|
2 |
|
|
|
|
2 |
|
|
α ö2 |
|
|
ï |
α ® ж |
|
|
æ |
|
ж ö |
ж |
æ |
ж ö |
ж |
æ |
|
α |
||||||
|
С |
|
|
|
|||||||||||||||
ïDG |
|
+ RT ln |
|
+ |
ç x |
|
÷ |
E |
|
+ ç x |
÷ |
E |
|
- ç x |
÷ |
E |
|
||
|
xα |
|
|
|
|
||||||||||||||
ï |
С |
|
|
è |
В |
ø |
|
ВС è |
А |
ø |
|
АС è |
В ø |
|
ВС |
ïС
ï |
+ xж xж (Eж + Eж - Eж ) - xα xα (Eα + Eα - Eα ) = 0 |
ï |
|
î |
А В АС ВС АВ А В АС ВС АВ |
æ |
α ö2 |
E |
α |
- ç x ÷ |
АС |
||
è |
С ø |
|
æ |
α ö2 |
E |
α |
- ç x ÷ |
ВС |
||
è |
С ø |
|
æ |
α ö2 |
E |
α |
|
- ç x |
÷ |
АС |
||
è |
|
А ø |
|
+
+ (6.4.8)
+
При условии, что величины ЕαАВ,ЕαАС,ЕαВС,ЕАВж ,ЕАСж ,ЕВСж равны нулю, система уравнений будет описывать двухфазное равновесие в приближении модели идеальных растворов.
Аналогичные системы уравнений можно записать для равновесия двух твердых растворов замещения (например, α и β).
Оценку параметров взаимодействия компонентов в равновесных фазах при использовании термодинамической модели регулярных растворов заме щения позволяет осуществить методика, предложенная Кауфманом и Берн стейном1.
С этой целью была создана электронная рабочая книга «InterEnergy.xls», лист которой «ВводВывод» показан на рисунке 6.4.5.
__________________
1Кауфман Л., Бернстейн Х. Расчет диаграмм состояния с помощью ЭВМ. М.,1972
136
Для ввода информации достаточно ввести химические символы инте ресующих элементов в соответствующие ячейки. Для введенных элементов автоматически производится подбор необходимых справочных данных (но мер группы, величина энтальпии, удельный объем, тип кристаллической ре шетки) и производится расчет.
Рисунок 6.4.5 Вид листа «ВводВывод» электронной книги “InterEnergy.xls”
В работе был разработан метод расчета положения коноды, позволяю щий определять состав и долю фаз сплава, находящихся в равновесии. Дан ная методика была апробирована на системе NbWMo (рисунок 6.4.6).
Результаты расчета для этой системы в приближении регулярных рас творов приведены на рисунке 6.4.7. Сравнение теоретически рассчитанных положений изотерм солидуса с экспериментальными данными показало удовлетворительную точность модели регулярных растворов для системы
NbWMo.
Подобные системы уравнений составлены для систем с четырьмя и большим числом компонентов.
137
Рисунок 6.4.6 Построение конод
Рисунок 6.4.7 Проекция изотерм солидуса на концентра ционный треугольник системы NbWMo, слева — резуль тат теоретических расчетов, справа — построение на основе справочных экспериментальных данных. (Изотермы подпи саны в градусах Цельсия)
Анализ уравнений показал, что основой любой многокомпонентной системы служат соответствующие двойные системы, которые в большинстве своем изучены достаточно хорошо. Поэтому, первое, что необходимо опре делить: какое термодинамическое приближение позволяет удовлетворитель но описать линии фазовых равновесий в рассматриваемых двухкомпонент ных системах. Цель этой процедуры выбор термодинамической модели для дальнейшего описания поверхностей фазовых равновесий в более сложных системах.
138
Контрольные вопросы к разделу 6
1.Назовите правила о соприкасающихся пространствах состояния и числе фаз в этих пространствах.
2.Приведите пример использования этих правил на изотермическом разрезе тройной диаграммы состояния.
3. Приведите пример использования этих правил на политермическом разрезе тройной диаграммы состояния.
4.Каковы геометрические основы изображения четверных и более сложных систем?
5.Какие расчетные методы можно использовать для нахождения кри тических точек тройных, четверных и более сложных сплавов?
6.Какие расчетные методы можно применять при построении диа грамм состояния тройных и более сложных систем?
139