过在普通碳素钢中加入合金元素改变钢铁的零界点从而改变热处理工艺[4]。

文玉华, 朱如曾, 周富信等人介绍了分子动力学的基本原理，有关的有限差分技 术，势函数的发展，初始条件和边界条件的选取，平衡态系综及其调控，感兴趣量的提 取，分子动力学的特别用途以及与其他计算方法的结合，最后还指出了分子动力学模拟 方法本身进一步研究的方向。同时指出分子动力学模拟最重要的两大要素是初始位形的 给定和原子间作用势[5]。崔守鑫, 胡海泉, 肖效光等人通过原子或分子运动的角度，介 绍了分子动力学方法的基本理论以及对一些感兴趣量的提取.并且还简要地说明了分子动 力学的用途以及与其它方法的结合[6]。进行分子动力学模拟的基础是选用合适势函数来 描述系统内原子间的相互作用势。陈强, 曹红红, 黄海波回顾了历史上出现的各种原子 间互作用势函数。对各种势模型的理论背景、具体形式、应用范围、特点进行了分析评 述[7]。李昶运等人运用分子动力学的理论,结合嵌入原子模型(EAM)势函数对金属铝、 铜以及铝-铜界面偶的扩散过程进行了理论模拟与计算。对于体扩散,在温度不太高时,由 于原子的跳跃频率很低,想要产生具有统计意义的原子位移需要经过很长的时间且超出了 分子模拟所能允许的时间尺度,所以对体扩散的分子进行模拟具有一定的困难[8]。徐云 波等人通过了解热变形对相变孕育期的影响，计算确定形变诱导相变的发生。温度取决 于变形条件和合金成分，随着碳含量的增大而减小。当变形诱导铁素体相变发生时，位 错密度的流动曲线会呈现明显的软化特性。通过高温大变形，铁素体转变提高，得到更 高的容积率。模拟结果与实验结果吻合良好[9]。张洪涛等人通过引入铁素体形成温度的 差异和多元体系的二元 Fe-C 参数，并参照热力学模型对 Fe-C 二元系统，多元结构钢先 共析铁素体形成一个简化的热力学模型提出了建议。得出结论，简化在多组分结构钢先 共析铁素体形成的热力学模型是可行的，计算出的温度的温度与实验的相吻合[10]。

然而，采用分子动力学方法直接研究铁碳合金的相变过程，直接观察原子运动的文 献并不多见。大多数的研究还停留在对于纯铁相变的模拟。Engin 和 Urbassek 模拟了分 析了常用于描述铁元素相变的势，得出了只有 Meyer-Entel 势才能双相描述奥氏体与马氏 体相变的结论[11],并采用该势成功模拟了纯铁中的相变过程并分析了相变热力学及动力 学[12]。Tateyama 等人采用分子动力学方法，在体系中插入符合西山关系及 K-S 关系的 马氏体/奥氏体相界面，并比较了不同初始相界条件下相变的异同[13]。这些作者们又改 良了 FS 势，第一次采用分子模拟学方法实现了模拟可逆的马氏体相变[14]。对于铁-碳 合金，王等人利用 Meyer-Entel 势描述铁原子间相互作用，Johnson 势描述铁-碳间的相互 作用，Tersoff 势描述碳-碳间的相互作用，得出了与实验结果一致的模拟结论[15]。