РАСЧЕТ ТЕРМОДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК СИСТЕМЫ Ре-Р МЕТОДОМ МОЛЕКУЛЯРНОЙ ДИНАМИКИ

摘要

Проблема дефосфоризации железоуглеродистых сплавов является актуальной для металлургической отрасли, так как высокая концентрация фосфора способствует появлению целого ряда крайне нежелательных явлений. Решению этой задачи посвящено множество экспериментальных работ, но справиться с ней окончательно до сих пор не удалось. Любые натурные эксперименты, направленные на исследование процесса удаления фосфора, требуют немалых материальных и временных затрат, но при этом не гарантируют получение нужного результата. Поэтому для поиска новых подходов к разрешению обозначенной проблемы гораздо рациональнее использование методов численного моделирования с привлечением вычислительных возможностей современных компьютеров. В настоящее время компьютерные эксперименты являются таким же признанным методом исследования, как теоретическое исследование и реальный эксперимент. Для изучения поведения атомов фосфора в железе с помощью численного эксперимента необходимо построить расчетную модель и апробировать ее путем вычисления различных характеристик, значения которых заранее известны. В качестве метода компьютерного моделирования в данной работе был выбран метод молекулярной динамики. С его помощью можно проводить эксперименты с заданными скоростями атомов и описывать динамику исследуемых процессов. Для описания межчастичного взаимодействия использовался потенциал, рассчитанный в рамках метода погруженного атома. Исследование проводилось на расчетной ячейке, имитирующей кристалл а-железа с атомами замещенггя фосфора. Построенная модель продемонстрировала удовлетворительные результаты при вычислении известных характеристик моделируемой системы. Были установлены зависимости изменения таких характеристик, как температурный коэффициент линейного расширения, температура плавления, скрытая теплота плавления и теплоемкость от концентрации атомов замещения фосфора, а также в отдельных случаях от величины приложенного внешнего давления. Расчеты показали, что, например, концентрация фосфора 0,5 % приводит к росту среднего теплового коэффициента линейного расширения на 9 %, снижению температуры и скрытой теплоты плавления на 5 %, а теплоемкости на 7 %.