Исследование фазовых равновесий "жидкость–твердое" трехкомпонентных сплавов легкоплавких металлов расчетным методом

Размещено на http://www.allbest.ru/

ИССЛЕДОВАНИЕ ФАЗОВЫХ РАВНОВЕСИЙ ЖИДКОСТЬ - ТВЕРДОЕ ТРЕХКОМПОНЕНТНЫХ СПЛАВОВ ЛЕГКОПЛАВКИХ МЕТАЛЛОВ РАСЧЕТНЫМ МЕТОДОМ

RESEARCH OF PHASE BALANCES LIQUID-SOLID OF TERNARY ALLOY MATERIALS OF LOW-MELT METALS BY CALCULATING METHOD

Доценко Сергей Павлович

Dotsenko Sergei Pavlovich

Dr. Sci. Agr., assistant professor

Арустамова Ирина Сергеевна

Arustamova Irina Sergeevna

Cand. Agr. Sci., senior research worker

Арутюнян Маргарина Мкртычевна

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

Article is devoted to the research and working out of calculation method of thermodynamic data and phase diagrams of ternary alloy materials of low-melt metals with the use of data bases known in literature.

Ключевые слова: РАСЧЕТНЫЙ МЕТОД, ТЕРМОДИНАМИЧЕСКИЕ ДАННЫЕ,

ФАЗОВЫЕ ДИАГРАММЫ, ТРЕХКОМПОНЕНТНЫЕ СПЛАВЫ, ЛЕГКОПЛАВКИЕ МЕТАЛЛЫ, БАЗЫ ДАННЫХ.

Key words: CALCULATION METHOD, THERMODYNAMIC DATA, PHASE DIAGRAMS, TERNARY ALLOY MATERIALS, LOW-MELT METALS, DATA BASES.

Условия, отвечающие фазовому равновесию многокомпонентных систем в случае отсутствия или с учетом приемлемой ошибки допущения отсутствия растворимости в твердом состоянии, содержат химические потенциалы всех компонентов, которые должны быть выражены через энергию Гиббса к-компонентной системы, выраженную в рамках принятого подхода через соответствующие энергии двойных систем.

Согласно работе [1], парциальную величину L_i экстенсивного свойства L можно представить следующим образом.

Пусть L относится к одному молю смеси, состоящей из n₁, n₂,...n_k молей компонентов. Тогда, по определению парциального свойства [1], можно записать:

(1)

Пусть L выражено в мольных долях компонентов, т.е.

(2)

Тогда L является функцией от k-1 независимых аргументов, поскольку один из них определяется из условия . Будем считать, что таким зависимым аргументом является мольная доля k-го компонента. Вычислим производную как производную от сложной функции:

(3)

С учетом (2) получаем

(4)

Подставляя (3) в (1), с учетом (4) получаем:

(5)

Запишем это выражение для k-го компонента:

(6)

В уравнении член отсутствует в силу принятой независимости L от x₁. Вычитая уравнение (5) из (6), получаем:

(7)

Полагая L свободной энергией Гиббса, уравнение (7) перепишем в виде:

(8)

что совпадает с полученным с помощью уравнения Гиббса - Дюгема уравнением в работе [1].

В нашем конкретном случае парциальные избыточные энергии Гиббса компонентов в тройной системе рассчитываем по уравнениям [2]:

(9)

Для расчета избыточных термодинамических свойств многокомпонентных систем обычно используют уравнение Колера [3], которое дает хорошее согласование с экспериментальными значениями, поэтому оно было использовано в данной работе.

Например, избыточную интегральную энергию Гиббса тройной системы можно представить уравнением:

, (10)

где - ДG_ij значения интегральных избыточных энергий Гиббса двойных систем при у_i = N_i/(N_i + N_j) в системе i-j;

в принятой нами системе обозначений их можно представить [2]:

, (11)

где y - атомная доля второго компонента в двойных системах.

Подставим соответствующие значения, с учетом того, что сумма атомных долей компонентов в тройной системе N₁+N₂+N₃=1, раскроем скобки и получим:

(12)

В эвтектических системах линию ликвидус находят из равенства химических потенциалов жидкой и твердой фаз [1]. Уравнение линии ликвидус запишем в виде:

(13)

где N_i - атомная доля, ДG_i(N_i, Т) - избыточная парциальная энергия Гиббса, ДH_i^пл, Т_i^пл - энтальпия и температура плавления i-го компонента.

Например, для тройной системы линии ликвидус находят из равенства химических потенциалов жидкой и твердой фаз [2]. Систему уравнений линий ликвидус в поле каждого компонента представим в виде:

(14)

Вычисление координат эвтектики реальных растворов сводится к решению системы уравнений (14) относительно Т и N_i [1].

Расчет диаграмм состояния тройных систем легкоплавких металлов является этапом разработки многокомпонентных эвтектических сплавов, которые могут быть использованы как ТАМ. Экспериментальное построение диаграммы состояния тройной системы требует проведения значительного количества измерений методом ДТА, которые можно с удовлетворительной точностью заменить расчетными данными.

Расчет диаграммы состояния тройных систем полезен, по нашему мнению, в связи с тем, что для разработки многокомпонентных ТАМ достаточны полуколичественные данные по всей диаграмме, а область эвтектики или монотектики может быть уточнена экспериментальными методами. Таким образом, значительно уменьшается объем экспериментальных исследований многокомпонентных систем.

В данной работе вычислены термодинамические свойства и диаграммы состояния двух эвтектических систем: галлий - олово - цинк и галлий - индий - цинк. В таблицах 1, 2 представлены избыточные интегральные энергии Гиббса и температуры ликвидуса: экспериментальные [3], рассчитанные по параметрам чистых компонентов [3] и по методике, приведенной в данной работе.

Таблица 1 - Рассчитанный и экспериментальный интегральный избыточный изобарный потенциал (ДG^изб) при 763 К, кДж/моль в трех разрезах системы галлий - цинк - олово N_Ga : N_Sn = 1 : 4

Таблица 2 - Рассчитанные и экспериментальные температуры ликвидуса в трех разрезах системы галлий - цинк - индий, N_Ga : N_In = 1 : 3. Температура ликвидуса, ⁰С

Температуры ликвидуса, рассчитанные по методике, приведенной в данной работе, достаточно близки к экспериментальным.

Составы и температуры плавления эвтектических сплавов, рассчитанные с использованием данных фазовых равновесий двойных систем, находятся в околоэвтектических областях, что существенно сокращает последующее экспериментальное уточнение состава эвтектик.

Галлий образует области расслаивания в двойных системах с кадмием, висмутом, свинцом, таллием и ртутью. Для определения областей расслаивания по данным для чистых компонентов в работе [3] рассчитывали парциальные избыточные энергии Гиббса жидких сплавов. По ним были рассчитаны области расслаивания в двойных системах галлия со свинцом, таллием, висмутом, кадмием и ртутью. Полученные результаты качественно отражают фазовые равновесия жидкость - жидкость в этих системах, хотя количественные результаты не удовлетворительны.

В работе [4] рассмотрена возможность расчета формы купола расслаивания в двойных системах на основании разных представлений энтальпии и энтропии смешения: через атомные, объемные и поверхностные доли. Учет различия атомных объемов позволяет объяснить асимметрию купола расслаивания. На этой основе рассчитаны линии расслаивания в системах Gа - Сd и Fe - Sп (согласие с опытом неудовлетворительное). Для оценки уровней избыточной свободной энергии при парных взаимодействиях, необходимых для построения модели, использовали известные данные по трем двойным системам. Адекватность полученной модели проверяли реконструкцией двойной системы Сd - Zп, показавшей достаточную точность расчетных методов. Получены кривые относительного изменения свободной энергии смешения для сплавов системы РЬ - Zn при 714, 805 и 877 К. Сопоставление расчетных кривых свободной энергии смешения при указанных температурах с экспериментальными данными по тройной системе указывает на удовлетворительное их соответствие.

Для тройных систем следует различать два варианта: область расслаивания может лежать целиком в поле одного компонента или в двух полях, т.е. пересекаться с пограничной кривой. На рисунке 1 изображена плоская диаграмма для первого варианта [5].

На рисунке 1 заштрихованной является область расслаивания, l₁ -эвтектическая точка в системе А-В, l₂ - эвтектическая точка в системе В-С, l₃ - эвтектическая точка в системе В-С, Е - точка тройной эвтектики. В подобном типе диаграмм состояния монотектический процесс моновариантный, неизотермический, поэтому для целей аккумулирования тепла и термостатирования тепловыделяющих элементов приборов такие сплавы не подходят.

Рисунок 1 - Плоская диаграмма состояния тройной системы, когда область расслаивания лежит в поле одного компонента

Если область расслаивания накладывается на два поля, то через нее проходит одна из кривых вторичного выделения. В этом случае приходится иметь дело с монотектическим четырехфазным равновесием, которое может быть выражено схемой [5]:

Ж₁ - А + В + Ж₂. (15)

Причем при отнятии теплоты от системы процесс идет направо, а при поглощении теплоты системой - налево. Это равновесие нонвариантное, так как в его осуществлении участвуют четыре фазы, поэтому в ходе этого процесса температура и концентрации жидких фаз остаются постоянными, и может измениться лишь число жидких и твердых фаз. Диаграмма состояния представлена на рисунке 2 [5].

Рисунок 2 - Плоская диаграмма состояния тройной системы с простой эвтектикой, когда область расслаивания накладывается на два поля

На диаграмме (см. рисунок 2) линия l₃d вторичного выделения А и С пересекает область гетерогенных жидких состоянии bкс в точках d и f . Вдоль линии l₃d идет процесс Ж₁ - А + С, а в точке d возникает вторая жидкая фаза состава, определяемого точкой f, и начинается четырехфазный нонвариантный изотермический монотектический процесс Ж₁ - А + С + Ж₂.

Точки d и f находятся на бинодали и одновременно на линии вторичных выделений. Внутри области расслаивания линия вторичных выделений становится коннодой, связывающей точки d и f, отвечающие растворам, равновесным с твердыми компонентами А и С. При продолжающемся отнятии теплоты количество первого жидкого слоя уменьшается, а количество фаз А, С и второго жидкого слоя увеличивается. Процесс закончится, когда вся первая жидкая фаза будет израсходована. Затем происходит вторичное выделение тех же твердых фаз А и С, но уже из второго жидкого слоя (Ж₂-А + С).

Этот процесс идет по линии Ш. Как и в случае системы без расслаивания, затвердевание заканчивается кристаллизацией тройной эвтектики в точке Е: Ж₂ - А + В + С. Таким образом, в подобном типе диаграмм есть тройная монотектика в точке d, координаты которой желательно рассчитать как можно точнее, и сплавы, соответствующие монотекттической точке, могут быть использованы в качестве изотермически плавящихся тепловых аккумуляторов.

Существует вариант диаграмм состояния с тремя компонентами, когда два компонента образуют монотектики с третьим компонентом (рисунок 3).

Рисунок 3 - Плоская диаграмма состояния тройной системы, когда имеются две области расслаивания

Искомая точка d тройной монотектики будет определяться двойной системой с большим куполом расслаивания.

В нашей работе была исследована возможность расчета координат монотектической точки в многокомпонентных сплавах [5].

Для монотектической точки можно записать систему уравнений:

(16)

Монотекгическая точка расположена всегда ближе к тугоплавкому компоненту [5] (компонент А (1) на рисунке 3). Поэтому, решая систему уравнений (16) относительно N₁', N₁", N₁, можно определить точку пересечения кривых 1М и КМ, т. е. координаты монотектической точки.

В изучаемых нами сплавах следующие двойные системы образуют монотектики: свинец - галлий, висмут - галлий и кадмий - галлий. Данные по величине области расслаивания в этих системах приведены в таблице 3 [6]. сплав уравнение многокомпонентный легкоплавкий

Таблица 3 - Границы области расслаивания на монотектической горизонтали и температура критической точки области расслаивания

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

Для расчета координат монотектической точки в тройной системе, например Рb - Вi - Gа (1 - 2 - 3), можно записать следующую систему уравнений:

(17)

Или, согласно принятой нами модели изобарно-изотермического потенциала:

(18)

где ДG_i^?изб и ДG_i^??изб - парциальные изобарно-изотермические потенциалы на границах расслаивания; им соответствуют мольные доли компонентов N_i^? и N_i^?? на границах области расслаивания.

Совместное решение такой системы уравнений относительно N₁, N₂, N₃, N₁^?, N₂^?? и Т, с учетом N₁ + N₂ + N₃ = 1, дает координаты монотектической точки.

Монотектические диаграммы состояния являются более сложными, чем эвтектические, так как дополнительно к линиям равновесия жидкость - твердое на диаграмме имеются линии равновесия жидкость - жидкость. Расчетное установление границ области расслаивания способствует более экономному нахождению монотектической точки. Поэтому представляет интерес расчет линий ликвидуса в тройных диаграммах с различной величиной области расслаивания.

Данные по температурам ликвидуса разрезов диаграммы состояния галлий - цинк - свинец приведены в таблице 4.

Таблица 4 - Температуры ликвидуса (°С) экспериментальные (t_эксп) и расчетные (t_расч) сплавов системы галлий - цинк - свинец. Концентрация компонентов в массовых долях (X_i). Отклонение расчетной температуры от экспериментальной, рассчитанной для наших данных (Дt)

Экспериментальные и расчетные данные подтверждают хорошую сходимость по температуре ликвидуса.

Список литературы

1. Глазов, В.М. Химическая термодинамика и фазовые равновесия / В.М. Глазов, Л.М. Павлова. - Металлургия, 1981. - 336 с.

2. Доценко, С.П. Исследование фазовых равновесий жидкость - твердое 2- компонентных сплавов легкоплавких металлов на основе баз термодинамических данных / С.П. Доценко, И.С. Арустамова, А.В. Лепеха, А.Б. Фурсина // Эл. журнал. www/ej.kubagro.ru. № 27(03), март, 2007. Идентиф. № ИНФОРМрегистра: 042070012?0075. 1 п. л.

3. Данилин В.Н. Физико-химические основы создания галлийсодержащих тепловых аккумуляторов: Автореф. дис. … д-ра хим. наук. - Краснодар, 1982. - 34 с.

4. Gluck R., Preddel B. Thermodynamiks of the decomposition of liguid alloys//Eff. Gravity solidification Immiscible Alloys Proc. RIT/ESA/SSC Workskop, 18 - 20. - Jar. - 1984. - Paris. -1984. - P.13-23.

5. Аносов, В.Л. Основы физико-химического анализа / В.Л. Аносов, М.И. Озерова, Ю.Я. Фиалков. - М.: Наука, 1976. - 686 с.

6. Элиот Р.П. Структуры двойных сплавов. - М.: Металлургия, 1970. - Т. 1. - 455 с., - Т. 2. - 472 с.

Размещено на Allbest.ru