固态相变的计算机模拟.pptVIP

用元胞自动机方法进行组织模拟 元胞自动机—Cellular automaton (CA) Von Neumann创立 定义：在空间，时间以及系统状态都是离散化的动力学系统。空间被一定形式的网格分割为许多单元，在这规则的空间格子中的每个单元，称为元胞（Cell)，它所具有的物理状态是系统有限数目状态中的一种状态。 本工作首先讨论不同形核位置对再结晶演变的影响。分别模拟了均匀形核、晶界形核、混合形核模型三种情况。 1.均匀连续形核 形核不考虑具体位置，再结晶晶核在体系中均匀分布而且连续形核。形核率采用下式计算： 其中， 为形变储能，用 为形核发生的临界储能，用 再结晶晶核的长大过程被视为热激活的过程，驱动力来自变形区域和再结晶区域之间的自由能差，生长过程中，晶粒的生长情况由位移来表示： 再结晶晶粒界面迁移速率用下式表示： 其中： M—界面迁移率 P—驱动力 其中M用下式来表示: 而界面迁移的驱动力来自变形产生的位错密度，用下式来表示： 其中： —位错密度 —剪切模量 —柏氏矢量的大小 用局部应变作为形核判据 晶体塑性有限元的计算结果 本文将晶体塑性有限元与元胞自动机方法进行耦合，建立了介观尺度的再结晶形核生长动力学模型，模拟了不同形核条件下再结晶过程。在力的作用下，由于变形的不均匀性，形核既可在晶界发生也可在晶内发生。 本工作建立起了非均匀变形再结晶形核的判据，分别以位错密度和局部应变作为判断依据。该判据尝试以最大位错密度的60%，70%，80%作为再结晶的临界形核判据，其准确性需要结合实验验证。 冷速1.0 ℃ /s的连续冷却相变过程中, 温度为716 ℃时的相变组织和碳浓度沿直线分布 奥氏体－铁素体连续冷却相变模拟 界面碳浓度不等且介于初始浓度和该温度下的平衡浓度之间，γ/α界面为非平衡界面。 716 ℃相变组织 碳浓度分布 变形态奥氏体－铁素体相变模拟 变形态奥氏体的特点 变形储能使奥氏体相能量增加—表现为奥氏体相中空位和位错密度大为增加，奥氏体晶粒内部产生变形带 变形使奥氏体晶粒被压扁或拉长，奥氏体晶界面积增加 在不同热变形条件下，Fe-0.19%C-0.445%Si-1.46%Mn钢真应力-真应变曲线 变形态奥氏体－铁素体相变模拟 变形结束时，三种变形条件下奥氏体相中的变形储能分别为： 11.9, 9.8, 6.6 J/mol。 三个试样分别在850℃以10/s的应变速率,900℃以10/s的应变速率,900℃以1/s的应变速率变形到真应变0.6停止变形,变形结束后立刻以15℃/s的冷速冷却，直到相变结束。 经历不同热变形后的奥氏体在冷速为15.0℃/s时的连续冷却相变动力学曲线 奥氏体中的变形存储能对相变后铁素体晶粒尺寸的影响 变形态奥氏体－铁素体相变模拟 (a)6.6J/mol (b)9.8J/mol (c)11.9J/mol 具有不同变形储能的奥氏体在连续冷却（15.0℃/s）铁素体相变结束时的相变组织及对应的碳浓度场 变形态奥氏体－铁素体相变模拟 (a) 未变形 (b)20% (c)50% (d)60% 压缩变形引起的奥氏体晶粒拓扑结构的变化(压扁拉长，晶界面积增加) 变形态奥氏体－铁素体相变模拟 模拟的不同冷速下的最终铁素体平均晶粒尺寸和实验结果的比较 3.塑性变形 ? 不同塑性变形条件下奥氏体晶粒结构 (a)压下量0%,(b) 压下量20%,(c) 压下量50%,(d) 压下量60% a b c d 4.变形态奥氏体分解为铁素体相变模拟结果 变形储能对形变驱动力的影响 a b c d 模拟的变形态奥氏体连续冷却过程中不同温度时的组织 (a)705℃,(b)690℃,(c)667℃, (d)633℃ 铁素体相变结束时体系中的碳浓度分布 形核模型： 长大模型： 其中， 驱动力P通过 计算 应变分布情况 位错密度分布情况 变形晶粒分布情况 形核位置分布情况 局部应变分布情况 在局部应变较大的区域形核密度较大，局部应变较小的区域形核密度较小甚至根本不形核。 以局部应变作为形核判据的再结晶情况下晶粒长大模拟 小结 * 固态相变的计算机 模拟 固态相变的计算机模拟 相变指的是当外界约束（温度或者压强）作连