Молекулярно-динамическое моделирование образования кластеров углерода при отпуске мартенсита Fe-C

Abstract

Дальнейшее улучшение свойств сталей невозможно без понимания процессов, происходящих на различных стадиях термообработки, на атомистическом уровне. В данной работе методом молекулярной динамики с использованием межчастичных потенциалов на основе модели погруженного атома (ЕАМ) проводилось моделирование мартенсита системы железоуглерод. Наблюдаемая диффузия углерода по октаэдрическим междоузлиям при повышенных температурах (несколько сотен °C) приводит к формированию атомами С ближнего упорядочения в виде периодических плоских скоплений, разделенных между собой областями решетки, практически свободными от углерода. Определено, что области кластеризации имеют ориентацию относительно решетки железа (102), что согласуется с результатами экспериментов, где исследовались структуры, получаемые при отпуске мартенсита на стадии двухфазного распада. Результаты атомистического моделирования показывают, что при кластеризации углерода происходит увеличение отношения параметров решетки c/a как в областях кристаллической решетки, в которых происходит формирование скоплений, так и в зонах, не содержащих атомы С. Толщина кластеров оказалась равной 17 А, а областей, не занятых углеродом - 30 А. Уменьшение энергии системы с течением наблюдаемого процесса можно рассматривать как движущую силу реакции с величиной 453,61 Дж/моль, что находится в качественном согласии с другими работами. Further improvement of the properties of steel is impossible without understanding of the atomic- level processes that take place at the different stages of heat treatment. In this work a simulation of iron- carbon martensite was performed using the method of molecular dynamics with interatomic potential based on the embedded atom model (EAM). The observed diffusion of carbon on octahedral interstices at high temperatures (at several hundreds of °C) causes the formation of short-range ordering of C atoms by way of periodical plain clusters divided by lattice regions, which almost do not contain carbon. We found that the cluster regions are orientated relative to iron lattice with (102) indices, what is consistent with the results of the experimental studies of the structures produced during martensite tempering at the stage of two-phase decomposition. The atomistic simulations results show that carbon clusterization causes the increasing of the lattice parameters relation c/a, both in the lattice regions where clusters are formed, and in the zones which do not contain any carbon atoms. The last fact is explained due to necessity of crystallographic coupling of these two zones. The thickness of the clusters turned out to equal 17 A, and that of the regions not filled with carbon - 30 A. During the simulation the total energy of modeling system decreases, and that can be considered as the reaction driving force with the value of 453,6 J/mole, which shows a qualitative agreement with other works.