Моделирование структуры бинарных упорядочивающихся сплавов в приближении жесткой кристаллической решетки тема автореферата и диссертации по физике, 01.04.07 ВАК РФ

Работа выполнена в ФГБОУ ВПО «Уфимский государственный авиационный технический университет» и в Сибирском физико-технического института им. акад. В.Д. Кузнецова при ФГБОУ ВПО «Национальный исследовательский Томский государственный университет».

доктор физико-математических наук, Дмитриев Сергей Владимирович доктор физико-математических наук, профессор Потекаев Александр Иванович.

Демьянов Борис Федорович доктор физико-математических наук, профессор, ФГБОУ ВПО «Алтайский государственный технический университет им. И.И. Ползунова»; Маркидонов Артем Владимирович, кандидат физико-математических наук, доцент, ФГБОУ ВПО «Кузбасский государственный технический университет им. Т.Ф. Горбачева».

ФГБОУ ВПО «Томский государственный архитектурно-строительный университет», Томск.

Защита состоится «21» ноября 2013 г. в 13:00 часов на заседании диссертационного совета Д 212.004.04 при ФГБОУ ВПО «Алтайский государственный технический университет им. И.И. Ползунова» по адресу: 656049 г. Барнаул, пр. Ленина, 46, e-mail: Veronika 65@mail.ru.

С диссертацией можно ознакомиться в библиотеке Алтайского государственного технического университета им. И.И. Ползунова.

Автореферат разослан «16» октября 2013 г.

Ученый секретарь диссертационного совета, кандидат физико-математических наук, доцент

Отзывы на автореферат с печатью в 2- экземплярах просим присылать на e-mail и адрес диссертационного совета АлтГТУ

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность и степень разработанности темы исследования. Упорядочивающиеся сплавы и интерметаллиды могут обладать уникальными физико-механическими свойствами (жаропрочностью, сопротивляемостью окислению, кавитации, коррозии, эффектом памяти формы и др.) [1 - 8]. По целому ряду признаков их можно выделить в отдельный класс материалов, свойства которых можно варьировать за счет изменения состава и степени упорядоченности. В упорядоченном состоянии компоненты сплава занимают узлы кристаллической решетки не случайным образом, создавая определенную сверхструктуру [1]. Процесс упорядочения сплава происходит в определенном температурном интервале. С одной стороны, температура должна быть достаточно высокой, чтобы обеспечить заметный коэффициент диффузии, а с другой стороны, она должна быть ниже температуры Курнакова, определяющей точку фазового перехода порядок-беспорядок [3]. В интерметаллидных соединениях точка Курнакова близка к точке плавления. Скольжение дислокаций в упорядоченном сплаве затруднено, поскольку оно может приводить к расщеплению дислокаций с образованием особого типа планарных дефектов - сдвиговых антифазных границ [3 - 9]. Расщепление сверхдислокаций может привести в определенных случаях к образованию сидячих конфигураций и запиранию некоторых систем скольжения, что и составляет эффект так называемого сверхструктурного упрочнения, который может приводить к парадоксальному росту предела текучести сплава с ростом температуры, как, например, в сплавах со сверхструктурой 1Л2.

Вопросам упорядочения сплавов, создания структур с заранее заданными физическими и механическими свойствами, прогнозу фазовых диаграмм состояния уделяется большое внимание [1 -9], поскольку они представляют не только научный, но и прикладной интерес.

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

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

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

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

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

• приближение парных межатомных связей;

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

Кроме того, в данной работе мы ограничиваемся рассмотрением лишь бинарных сплавов состава АпВт •

Хорошо известно, что межатомные взаимодействия можно представить в виде разложения на парные, трехчастичные и т.д. взаимодействия [10]. Для многих металлов и сплавов непарные взаимодействия оказываются значительными, однако, степень непарности различна для разных материалов и её всегда можно рассматривать как поправку к двухчастичным взаимодействиям.

Пренебрежение энергией релаксации атомов кажется естественным после того как оказались отброшенными силы непарного взаимодействия.

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

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

Для достижения поставленной цели решали следующие задачи:

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

2) Изучение допустимых изменений параметров ближнего порядка для двумерных и трехмерных структур стехиометрических составов А„Вт при произвольных значениях энергий парных межатомных взаимодействий.

3) Исследование кинетики упорядочения сплавов состава А„Вт как для модельных сплавов на двумерных решетках, так и для сплавов на основе ГЦК и ОЦК решеток.

4) Исследование энергии сублимации бинарного сплава с учетом влияния состава сплава на межатомные взаимодействия.

1) Разработана математическая модель и её программная реализация, описывающая процесс упорядочения сплава состава А„Вт, заданного на произ-

вольной решетке, по вакансионному механизму диффузии в приближении жесткой кристаллической решетки и парных межатомных взаимодействий.

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

3) Показано, что кинетика упорядочения сплавов по вакансионному механизму зависит от параметров парного взаимодействия ^, , р^ даже при неизменных энергиях упорядочения ц =<?% +<Рвв гДе '

номер координационной сферы. В случае обменного механизма диффузии кинетика упорядочения зависит только от энергий упорядочения сплава , по не от конкретных значений <р%.

4) Показано, что при учете линейного изменения энергий парных взаимодействий от состава сплава становится возможным моделирование некоторых элементов фазовых диаграмм состояния, включающих промежуточные фазы стехиометрических составов Лп 1!т.

Теоретическая и практическая значимость работы.

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

Методология и методы исследования.

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

Положения, выносимые на защиту:

Математическая модель процесса упорядочения сплавов состава А„Вт по вакансионному механизму диффузии в приближении жесткой кристаллической решетки.

Алгоритм расчета диапазонов изменения параметров ближнего порядка для структур стехиометрических составов Л„В,„, заданных на различных решетках, при произвольных значениях энергий парных межатомных взаимодействий в нескольких первых координационных сферах.

Утверждение о том, что вакансионный механизм диффузии приводит к зависимости кинетики упорядочения от соотношения значений парных межатомных потенциалов , (fi^, <,->%, в то время как при обменном механизме диффузии кинетика упорядочения зависит только от энергии упорядочения <4 = <$л + <Рев

2<Рлв > гДе ' _ номер координационной сферы.

Вывод о существенном увеличении числа структур, реализуемых в сплавах состава , при учете линейной зависимости энергий парных межатомных взаимодействий от состава сплава.

Диссертационная работа выполнена в соответствии с Планом НИР Института проблем сверхпластичности металлов РАН по теме на 2008-2011 гг. «Исследование нелинейных явлений в твердых телах, далеких от термодинамического равновесия» (№ гос. per. 01200951797). Работа была частично поддержана Российским фондом фундаментальных исследований, гранты РФФИ 09-08-00695_а, 11-08-97057-р_поволжье_а. Получена финансовая поддержка по результатам конкурса УМНИК в 2011г.

Обоснованность и достоверность результатов

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

Апробация работы. Результаты диссертационной работы доложены и обсуждены на следующих научных форумах: V Всероссийской научно-практической конференции с международным участием «Инновационные технологии и экономика в машиностроении», Юрга (2007); Всероссийской молодежной научной конференции «Мавлютовские чтения» посвященная 75-летию УГАТУ, Уфа (2008); Международной молодежной научной конференции «XXXIV Гагаринские чтения» «XXXVI Гагаринские чтения», Москва (2008, 2010); 3-ей, 4-ой, 5-ой, 7-ой Всероссийской зимней школы-семинара аспирантов и молодых ученых (с международным участием) «Актуальные проблемы науки и техники», Уфа (2008, 2009, 2010, 2012); Первые московские чтения по проблемам прочности материалов посвященные 85-летию со дня рождения профессора В.Л. Лебедева и 90-летию со дня рождения профессора Л.М. Утевско-го, Москва (2009); X Международной научно-практической конференции «Проблемы и перспективы развития литейного, сварочного и кузнечно-штамповочного производств», Барнаул (2009); IV Международной научной конференции студентов, аспирантов и молодых ученых «Научный потенциал студенчества в XXI веке», Ставрополь (2010); 16-ой Всероссийской научной конференции студентов-физиков и молодых ученых (ВНКСФ-16), Волгоград (2010); XIX Петербургских чтениях по проблемам прочности посвященных

130-летию со дня рождения академика АН УССР H.H. Давиденкова, Санкт-петербург (2010); XII Всероссийской молодежной школе-семинаре по проблемам физики конденсированного состояния вещества (СПФКС-12), Екатеринбург (2011); Международной школы-конференции для студентов, аспирантов и молодых ученых «Фундаментальная физика и ее приложения в естествознании», Уфа (2011); XV Международной научной конференции, посвященной памяти генерального конструктора ракетно-космических систем академика М.Ф. Решетнева «Решетневские чтения», Красноярск (2011).

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

Публикации. Основное содержание диссертационной работы отражено в 36 научно-технических публикациях, из них - 9 публикаций в изданиях, включенных в перечень рекомендованных ВАК, 5 свидетельств о государственной регистрации программ для ЭВМ и 3 свидетельства о государственной регистрации баз данных.

Структура и объем диссертации. Диссертационная работа состоит из введения, четырёх глав, выводов и списка литературы. Работа изложена на 122 страницах, содержит 53 рисунков и 16 таблиц. Список литературы включает 139 наименований.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

Во введении обоснована актуальность выбранной темы исследования, сформулированы цель, научная новизна и практическая значимость диссертационной работы. Приведено краткое описание работы по главам.

ГЛАВА 1. УПОРЯДОЧИВАЮЩИЕСЯ СПЛАВЫ: ОБЛАСТИ ПРИМЕНЕНИЯ, СТРУКТУРА, СВОЙСТВА И МЕТОДЫ МОДЕЛИРОВАНИЯ. В данной главе, на основе анализа литературных данных, прежде всего, обсуждаются особые свойства упорядочивающихся сплавов и интерметаллидов, используемые в ряде современных технологий. Во втором разделе рассмотрены структура и свойства упорядочивающихся сплавов. Далее обсуждаются особенности фазовых переходов порядок-беспорядок в твердых телах. Четвертый раздел посвящен описанию численных подходов, используемых при моделировании упорядочивающихся сплавов и интерметаллидов. В частности, упоминаются такие методы как первопринципные расчеты, метод молекулярной динамики, метод Монте-Карло, метод кристаллогеометрического анализа упорядоченных сплавов, модели Изиига для описания фазовых переходов порядок-беспорядок. Описаны основные предположения, используемые в данной работе при построении компьютерных моделей, основными из которых, как уже отмечалось, являются учет только парных составляющих межатомных взаимодействий и принятие модели жесткой кристаллической решетки.

На основании сделанного литературного обзора сформулированы частные задачи исследования и описаны подходы к их решению, применяемые в данной диссертационной работе.

ГЛАВА 2. ДОПУСТИМЫЕ ОБЛАСТИ ИЗМЕНЕНИЯ ПАРАМЕТРОВ БЛИЖНЕГО ПОРЯДКА СПЛАВОВ НА ОСНОВЕ РЯДА ДВУМЕРНЫХ И ТРЕХМЕРНЫХ РЕШЕТОК. В данной главе предложен алгоритм определения границ изменения параметров порядка для структуры стехиометрии АпВт, заданной на произвольной решетке, при учете парных межатомных взаимодействий в I первых координационных сферах. При решении данной задачи учитываются только структурные характеристики сплава, а величины энергий межатомных взаимодействий не вовлекаются в рассмотрение. Алгоритм применен для сплавов АпВт стехиометрического состава при учете взаимодействий атомов в двух первых координационных сферах. Рассмотрены двумерные квадратная и гексагональная решетки, а также ГЦК и ОЦК решетки.

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

Атомы сортов А и В в сплаве стехиометрического состава А„Вт располагаются в узлах решетки, у которой на I -й координационной сфере имеется N1 атомов. Концентрации атомов сортов А и В равны

Через р% обозначим вероятность того, что на /-й координационной сфере атома сорта А находится атом сорта В. В бинарной структуре АпВт существуют следующие связи между вероятностями р\\, р%, р% и р\[]л:

Обозначим через энергию связи пары атомов сортов А и В, расположенных на расстоянии равном радиусу г'-й координационной сферы.

Потенциальная энергия структуры в расчете на один атом, при учете взаимодействия атомов в I первых координационных сферах, запишется в виде

Энергия полностью разупорядоченного состояния структуры определяется выражением (3) для = р'^ =сл, р^ = р% = св, что дает

Выберем энергию £Л50п1. в качестве точки отсчета, и будем характеризовать энергию любой структуры разностью

где введены параметры ближнего порядка

и энергии упорядочения

При выводе соотношения (5) использовались соотношения (1, 2) и тот факт, что сА + с,. = I.

Как видим из (5), энергия рассматриваемой структуры однозначно определяется координационными числами .V,, энергиями упорядочения о>. , а также параметрами а,, которые, по сути, являются параметрами ближнего порядка Каули [11].

Для определения диапазона изменения параметров порядка а1 , где / = 1,2. / - номер координационной сферы, использовался метод концентрационных волн [1], адаптированный для целей данного исследования.

Для трехмерного кристалла рассматривается расчетная ячейка с наложенными периодическими граничными условиями, содержащая N:txNyxNг атомов. Далее вычисляются значения

мх лгу лгг ! - - - ' ч • (8) 2л, к? 27! 2я1к?> 1 ^ +-— +--— + +С,

где 0<1<МХ, 0<7<Л'11, 0<1<Ыг. Векторы кт,кт,к0) в выражении (8) имеют целочисленные значения 0<4°,*®.*® <МХ , О^'Д^Д-™ ,

()<к'.[>,к',2',к[3> <№,. Сдвиги фаз ех,Ег,е3 введены для того, чтобы избежать обращения чисел С¡¡, в ноль. Константа С определяет стехиометрию структуры. Для структуры состава ДА используется

С помощью чисел Су( производится заполнение расчетной ячейки атомами сортов А и В таким образом, что узлу кристаллической решетки с номером (»',Л0 приписывается сорт А если С„, >0 и сорт В если С^ <0. После заполнения ячейки атомами проверяется соблюдение стехиометрии и, если она нарушена, данная структура не рассматривается. Для структур со стехиометрией рассчитываются параметры ближнего порядка а,, ; = 1 . 1 и отображаются точкой в 1 -мерном пространстве (аВ5 в приближении жесткой кристаллической решетки. / А.Р. Халиков, C.B. Дмитриев, А.И. Потекаев // Письма о материалах. - 2013. - Т.З - С. 225-229.

Программы для ЭВМ:

10. Свидетельство о государственной регистрации программы для ЭВМ №2012613171 «Моделирование процесса упорядочения сплавов по вакансиониому механизму в модели твердых сфер» / А.Р. Халиков, A.M. Искандаров, C.B. Дмитриев //М.: РосПатент, 03.04.2012.

11. Свидетельство о государственной регистрации программ для ЭВМ №2008613834 «Программа генерирования элементов эвтектических сплавов» / А.Р. Халиков, A.A. Танеев, К.В. Мавродиев, A.A. Сягаев // М.: РосПатент, 11.08.2008.

12. Свидетельство о государственной регистрации программы для ЭВМ №2008613833 «Расчет эвтектических концентраций и температур» / А.Р. Халиков, A.A. Танеев, К.В. Мавродиев, A.A. Сягаев // М.: РосПатент, 11.08.2008.

13. Свидетельство о государственной регистрации программы для ЭВМ №2008615428 «Расчет эвтектических концентраций бинарных систем» / А.Р. Халиков// М.: РосПатент, 14.11.2008.

14. Свидетельство о государственной регистрации программы дм ЭВМ №2009616486 «Моделирование эвтектических концентраций многокомпонентных систем» / Халиков А.Р // М.: РосПатент, 23.11.2009.

Свидетельства о регистрации баз данных:

15. Свидетельство о государственной регистрации базы данных №2008620302 «База данных по эвтектическим сплавам двухкомпонент-ных систем» / А.Р. Халиков, A.A. Танеев, А.А Сягаев // М.: РосПатент, 11.08.2008.

16. Свидетельство о государственной регистрации базы данных №2009620550 «Эвтектические составляющие» / А.Р. Халиков // М.: РосПатент, 23.11.2009.

17. Свидетельство о государственной регистрации базы данных №2010620086 «Элементы образующие многокомпонентные диаграммы состояния эвтектического типа» / Халиков А.Р. // М.: РосПатент, 07.12.2009.

Статьи в сборниках научных трудов:

18. Халиков, А.Р. Базы данных для создания композиционных материалов / А.Р. Халиков, A.A. Танеев // Ползуновский альманах. -2006. - №3 -С.157-160.

19. Халиков, А.Р. Метод расчета эвтектических концентраций бинарных систем / А.Р. Халиков, A.A. Танеев // Ползуновский альманах. -2008. - №3 -С.47-48.

20. Халиков, А.Р. Моделирование эвтектических концентраций бинарных и тройных диаграмм состояния металлических систем / А.Р. Халиков // Инженерная физика. - 2010. №2,- С.25-35.

21. Халиков, А.Р. Энергия сублимации бинарного сплава с учетом влияния состава сплава на межатомные взаимодействия / А.Р. Халиков, А.М. Искандаров // Письма о материалах. - 2011. - Т1 вып.:4. - С. 226- 230.

Публикации в материалах международных, всероссийских и региональных

22. Танеев, A.A. Расчет температур и упрочняющей фазы для получения композиционных материалов методом пропитки. Сб. статей V Всероссийской научно-практической конференции с международным участием «Инновационные технологии и экономика в машиностроении» / A.A. Танеев, Р.Р. Кабиров, А.Р. Халиков // Юрга: Изд-во: «ТПУ», -2007. - С.291-297.

23. Халиков, А.Р. Разработка методики расчета эвтектических концентраций и температур Сб. трудов Всероссийской молодежной научной конференции, посвященной 75-летию УГАТУ «Мавлютовские чтения» / Халиков А.Р.. A.A. Танеев // Уфа: Изд-во «Уфимск. Гос. Авиац. техн. Ун-т», -2007.-Т. 2-С. 167- 168.

24. Халиков, А.Р. Автоматизированная система расчета эвтектических концентраций и температур диаграмм состояния. Сб. статей III всероссийской зимней школы-семинара аспирантов и молодых ученых «Актуаль-

ные проблемы в науке и технике» / Халиков А.Р., К.В. Мавродиев, A.A. Танеев // Уфа: Изд-во «Диалог», - 2008. Т. 2. - С. 18- 24.

25. Танеев, A.A. Определение элементов образующих многокомпонентные эвтектические сплавы. Научные труды Международной молодежной научной конференции в 8 томах / Танеев A.A., А.Р. Халиков // Москва: МАТИ Ответственный редактор Н.И. Сердюк, - 2008. Т.1. - С.210- 212.

26. Халиков, А.Р. Определение эвтектических концентраций бинарных и многокомпонентных систем. Сб. статей IV всероссийской зимней школы-семинара аспирантов и молодых ученых, «Актуальные проблемы в науке и технике» / А.Р. Халиков, A.A. Танеев // Уфа: Изд-во «Диалог», - 2009. -С.390- 393.

27. Халиков, А.Р. Прогноз эвтектических концентраций бинарных и тройных диаграмм состояния металлических систем. Сб. трудов X Международной научно-практической конференции «Проблемы и перспективы развития литейного, сварочного и кузнечно-штамповочного производства» /

A.Р. Халиков // Барнаул: Изд-во «АлтГТУ» Под общ. ред. А.М. Гурьева и

B.А. Маркова, - 2009. Вып.5. - С.208- 210.

28. Халиков, А.Р. Моделирование эвтектических концентраций бинарных и многокомпонентных диаграмм состояния. Сб. трудов «Первые московские чтения по проблемам прочности материалов посвященные 85-летию со дня рождения профессора B.JI. Лебедева и 90-летию со дня рождения профессора Л.М. Утевского» / А.Р. Халиков // Москва. - 2009. - С. 160.

29. Халиков, А.Р. Эвтектические составляющие элементов. Сб. трудов Международной молодежной научной конференции в 8 томах «XXXVI Гага-ринские чтения» / А.Р. Халиков, A.A. Танеев // Москва: МАТИ. — 2010. — Т.1.-С. 124-126.

30. Халиков, А.Р. Расчет эвтектических концентраций трехкомпонентных систем. Сб. трудов Международной научной конференции по естественно-научным и техническим дисциплинам «Научному прогрессу - творчество молодых» / А.Р. Халиков, A.A. Танеев // Йошкар-Ола: Изд-во «МГТУ»,-2010.-4.1 -С. 312-313.

31. Халиков, А.Р. Математическое моделирование физического процесса образования эвтектики. Сб. трудов IV Международной научной конференции студентов, аспирантов и молодых ученых «Научный потенциал студенчества в XXI веке» / А.Р. Халиков // Ставрополь: Изд-во «СевКав-ГТУ», - 2010. - Т.1. - С.104- 105.

32. Халиков, А.Р. Эвтектические составляющие элементов и соединений и их связь различными физико-химическими величинами. Сб. трудов XVI Всероссийской научной конференции студентов-физиков и молодых ученых (ВНКСФ-16) / А.Р. Халиков, A.A. Танеев // Екатеринбург, Волгоград: Изд-во «АСФ России», - 2010. - С.769- 770.

33. Халиков, А.Р. Расчет эвтектических концентраций диаграммы состояния Pb-Sn-Bi. Сб. трудов XVI Всероссийской научной конференции студентов-физиков и молодых ученых (ВНКСФ-16) / А.Р. Халиков, A.A. Танеев // Екатеринбург, Волгоград: Изд-во «АСФ России», - 2010. - С.769- 770.

34. Халиков, А.Р. Влияние стехиометрии на энергии бинарных сплавов в ра-зупорядоченном состоянии и после отжига. Сб. трудов XV Международной научной конференции, посвященной памяти генерального конструктора ракет-космич. систем акад. М.Ф. Решетнева «Решетневские чтения» / А.Р. Халиков, A.M. Искандаров, C.B. Дмитриев // Краснярск: под общ. ред. Ю.Ю. Логинова, Сиб. гос. аэрокосмич. ун-т, - 2011. - Ч. 2. - С.517-518.

35. Халиков, А.Р. Учет влияния состава на энергию сублимации бинарного сплава. Фундаментальная физика и ее приложения в естествознании. Сб. трудов Международной школы-конференции для студентов, аспирантов и молодых ученых «Фундаментальная математика и её приложения в естествознании» / А.Р. Халиков, C.B. Дмитриев // Уфа: РИЦ БашГУ отв. ред. Альмухаметов Р.Ф.,-2011. - С.66.

36. Халиков, А.Р. Зависимости энергий парных взаимодействий от состава для структур АпВга ГЦК решетки. Актуальные проблемы в науке и технике. ТЗ. Управление в социально-экономических системах. Естественные науки. Сб. трудов 7 Всероссийской зимней школы-семинара аспирантов и молодых ученых / А.Р. Халиков, A.M. Искандаров, Ф.Р. Кувандыков, C.B. Дмитриев // Уфимск. гос. авиац. тех. ун-т. - Уфа: УГАТУ, - 2012. - С. 350-353.

Подписано в печать 14/10/13. Формат 60x84 1/16. Бумага офсетная. Печать ризографическая. Тираж 100 экз. Заказ 1158 Гарнитура «TimesNewRoman». Отпечатано в типографии «ПЕЧАТНЫЙ ДОМЪ» ИП ВЕРКО. Объем 1,4 п.л. Уфа, Карла Маркса 12 корп. 5. т/ф: 8(347) 27-27-600, 27-29-123

ФГБОУ ВПО УФИМСКИЙ ГОСУДАРСВЕННЫЙ АВИАЦИОННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ФГБОУ ВПО НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ ТОМСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

На правах рукописи

Халиков Альберт Рашитович

МОДЕЛИРОВАНИЕ СТРУКТУРЫ БИНАРНЫХ УПОРЯДОЧИВАЮЩИХСЯ СПЛАВОВ В ПРИБЛИЖЕНИИ ЖЕСТКОЙ КРИСТАЛЛИЧЕСКОЙ РЕШЕТКИ

01.04.07 - Физика конденсированного состояния

Диссертация на соискание ученой степени кандидата физико-математических наук

Доктор физико-математических наук Дмитриев Сергей Владимирович

Доктор физико-математических наук Потекаев Александр Иванович

ГЛАВА 1. УПОРЯДОЧИВАЮЩИЕСЯ СПЛАВЫ: ОБЛАСТИ ПРИМЕНЕНИЯ, СТРУКТУРА, СВОЙСТВА И МЕТОДЫ

1.1. Значение упорядочивающихся сплавов и интерметаллидов в технике. 14

1.2. Структура и свойства упорядочивающихся сплавов. 14

1.3. Параметры ближнего и дальнего порядка. 19

1.4. Фазовые переходы. Переход порядок-беспорядок. Точка Курнакова. 21

1.5. Численные подходы к изучению упорядочивающихся сплавов и интерметаллидов. 23

1.5.1. Первопринципные расчеты. 23

1.5.2. Метод молекулярной динамики. 25

1.5.3. Метод Монте-Карло. 27

1.5.4. Кристаллогеометрический анализ упорядоченных сплавов 28

1.5.5. Модель Изинга для описания фазового перехода порядок-беспорядок . 29

1.5.6. Методы, применяемые в данной работе. 30

1.6. Учет зависимости энергий связи пары атомов от концентрации. 31

1.7. Постановка задач исследования. 35

ГЛАВА 2 ДОПУСТИМЫЕ ОБЛАСТИ ИЗМЕНЕНИЯ ПАРАМЕТРОВ БЛИЖНЕГО ПОРЯДКА СПЛАВОВ НА ОСНОВЕ РЯДА ДВУМЕРНЫХ

И ТРЕХМЕРНЫХ РЕШЕТОК. 37

2.1. Принимаемые допущения. 37

2.2. Потенциальная энергия структуры. 38

2.3. Сопоставление вакансионного и обменного механизма диффузии 40 2.3.1. Обменный механизм диффузии. 40

2.3.2. Вакансионный механизм диффузии. 42

2.4. Допустимые области изменения параметров дальнего порядка . 44

2.4.1. Определение границ изменения параметров порядка для структур стехиометрии АпВт. 45

2.4.2. Сплавы стехиометрии АпВт на квадратной и гексагональной решетке. 46

2.4.2.1. Сплав АВ на квадратной решетке. 46

2.4.2.2. Сплав АгВ на квадратной решетке. 48

2.4.2.3. Сплав АЪВ на квадратной решетке. 49

2.4.2.4. Сплав А3В5 на квадратной решетке. 50

2.4.2.5. Сплав АВ на гексагональной решетке. 50

2.4.3. Сплавы стехиометрии АВ на ОЦК и ГЦК решетке. 50

2.4.3.1. Сплав АВ на ГЦК решетке. 53

2.4.3.2. Сплав АВ на ОЦК решетке. 53

2.5. Выводы по главе. 54

ГЛАВА 3. МОДЕЛИРОВАНИЕ КИНЕТИКИ УПОРЯДОЧЕНИЯ БИНАРНОГО СПЛАВА ПО ВАКАНСИОННОМУ МЕХАНИЗМУ ДИФФУЗИИ В ПРИБЛИЖЕНИИ ЖЕСТКОЙ КРИСТАЛЛИЧЕСКОЙ

3.1. Описание модели. 56

3.1.1. Температура Курнакова. 58

3.2. Влияние энергий парных связей на кинетику упорядочения модельных сплавов. 58

3.2.1. Результаты моделирования кинетики упорядочения на квадратной и гексагональной решетке. 59

3.2.2. Результаты моделирования кинетики упорядочения сплавов АВ и АЪВ на основе ГЦК решетки. 63

3.2.3. Результаты моделирования кинетики упорядочения сплавов АВ на основе ОЦК решетки. 67

3.3. Выводы по главе. 70

ГЛАВА 4. ЭНЕРГИЯ СУБЛИМАЦИИ БИНАРНОГО СПЛАВА С УЧЕТОМ ВЛИЯНИЯ СОСТАВА СПЛАВА НА МЕЖАТОМНЫЕ ВЗАИМОДЕЙСТВИЯ. 73

4.1. Влияние стехиометрии на энергии бинарных сплавов состава АпВт на основе ГЦК решетки в разупорядоченном состоянии и после отжига. 73

4.2. Влияние линейной зависимости энергий парных связей от концентрации для сплавов состава АпВт, на основе ГЦК решетки в разупорядоченном состоянии и после отжига. 77

4.3 Моделирование термических кривых охлаждения в процессе

упорядочения сплавов стехиометрии АВ, А3В и А2В5 95

4.4. Выводы по главе. 101

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ. 103

СПИСОК ЛИТЕРАТУРЫ. 105

Актуальность и степень разработанности темы исследования.

Упорядочивающиеся сплавы и иитерметаллиды могут обладать уникальными физико-механическими свойствами (жаропрочностью, сопротивляемостью окислению, кавитации, коррозии, эффектом памяти формы и др.) [1 - 8]. По целому ряду признаков их можно выделить в отдельный класс материалов, свойства которых можно варьировать за счет изменения состава и степени упорядоченности. В упорядоченном состоянии компоненты сплава занимают узлы кристаллической решетки не случайным образом, создавая определенную сверхструктуру [1]. Процесс упорядочения сплава происходит в определенном температурном интервале. С одной стороны, температура должна быть достаточно высокой, чтобы обеспечить заметный коэффициент диффузии, а с другой стороны, она должна быть ниже температуры Курнакова, определяющей точку фазового перехода порядок-беспорядок [3]. В интерметалл идных соединениях точка Курнакова близка к точке плавления. Скольжение дислокаций в упорядоченном сплаве затруднено, поскольку оно может приводить к расщеплению дислокаций с образованием особого типа планарных дефектов -сдвиговых антифазных границ [3 - 9]. Расщепление сверхдислокаций может привести в определенных случаях к образованию сидячих конфигураций и запиранию некоторых систем скольжения, что и составляет эффект так называемого сверхструктурного упрочнения, который может приводить к парадоксальному росту предела текучести сплава с ростом температуры, как, например, в сплавах со сверхструктурой Ь\2.

Вопросам упорядочения сплавов, создания структур с заранее заданными физическими и механическими свойствами, прогнозу фазовых диаграмм состояния уделяется большое внимание [1 —9], поскольку они представляют не только научный, но и прикладной интерес.

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

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

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

• приближение парных межатомных связей;

Кроме того, в данной работе мы ограничиваемся рассмотрением лишь бинарных сплавов состава АпВт.

Цель работы: установление связи между энергиями парных межатомных взаимодействий и структурой, а также кинетикой упорядочения сплавов состава А„Вт по вакансионному механизму диффузии в приближении жесткой кристаллической решетки.

Для достижения поставленной цели решали следующие задачи:

2) Изучение допустимых изменений параметров ближнего порядка для двумерных и трехмерных структур стехиометрических составов АпВт при произвольных значениях энергий парных межатомных взаимодействий.

3) Исследование кинетики упорядочения сплавов состава АпВт как для модельных сплавов на двумерных решетках, так и для сплавов на основе ГЦК и О ЦК решеток.

1) Разработана математическая модель и её программная реализация, описывающая процесс упорядочения сплава состава АпВт , заданного на произвольной решетке, по вакансионному механизму диффузии в приближении жесткой кристаллической решетки и парных межатомных взаимодействий.

3) Показано, что кинетика упорядочения сплавов по вакансионному механизму зависит от параметров парного взаимодействия <р^]в даже при

неизменных энергиях упорядочения ^,=<Раа+^вв

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

конкретных значений ср^А, ср(в]в, (р

📎📎📎📎📎📎📎📎📎📎