金属学报(中文版) 2018 , 54 (8): 1204-1214 https://doi.org/10.11900/0412.1961.2017.00478

Orginal Article

基于序参量梯度的改进相场模型及对大尺度体系马氏体相变的模拟研究

魏铖¹^,², 柯常波², 马海涛¹^,³, 张新平²^,

1 华南理工大学土木与交通学院广州 510640
2 华南理工大学材料科学与工程学院广州 510640
3 广州大学工程抗震研究中心广州 510405

A Modified Phase Field Model Based on Order Parameter Gradient and Simulation of Martensitic Transformation in Large Scale System

WEI Cheng¹^,², KE Changbo², MA Haitao¹^,³, ZHANG Xinping²^,

1 School of Civil Engineering and Transportation, South China University of Technology, Guangzhou 510640, China
2 School of Materials Science and Engineering, South China University of Technology, Guangzhou 510640, China
3 Earthquake Engineering Research & Test Center, Guangzhou University, Guangzhou 510405, China;

中图分类号: TF777.1

文章编号: 0412-1961(2018)08-1204-11

通讯作者: 通讯作者张新平,mexzhang@scut.edu.cn,主要从事电子封装及可靠性、形状记忆合金、计算材料学与加工模拟方面的研究

收稿日期: 2017-11-15

网络出版日期: 2018-08-11

基金资助: 国家自然科学基金项目No.51205135和广东省自然科学基金重点项目No.S2013020012805资助

作者简介:

作者简介魏铖,男,1987年生,博士生

展开

摘要

基于序参量变化在界面区域远大于非界面区域的特性,构造了全局修正函数,建立了适用于可调尺度体系马氏体相变模拟的相场模型。在不改变界面能密度的情况下,通过调整界面区域的体积自由能密度差与梯度能系数,有效增大了原相场模型中的界面宽度,实现了大尺度下的高效模拟,并能很好地表征马氏体相变。结果表明,改进后的相场模型能很好地解决原相场模型在大尺度体系模拟时存在的如生长速率过快、位向关系不合理及组织形貌杂乱无序等问题,模拟结果与实验结果符合较好。

关键词： 马氏体相变 ; 序参量梯度 ; 界面宽度 ; 相场法

Abstract

The materials design and fabrication based on predicting microstructure have been drawn increasing attention from scientists and engineers. Martensite microstructure, which is well observed in many materials, has significant influence on physical and mechanical properties of the materials. Some experimental studies have been launched to understand the featured microstructure and its evolution in martensitic transformations (MT). Meantime, numerical approaches are often employed to assist the experimental studies due to the complex and nonlinear nature of MT. The phase field method is one of the most powerful tools in predicting microstructure. Due to the diffuse-interface description, phase field method can be used to simulate arbitrary morphologies without tracking the interface. As a consequence, the interface must contain enough elements to obtain reasonable results by using finite element method. On the other hand, the width of the interface is several orders smaller than its real counterpart. More computational resources are required to resolve the phase field variables at the interface with the system size increased. Therefore, the simulation is restricted in smaller system even with state-of-the-art computer power. For arbitrary model formulations, the interfacial energy depends on the interfacial width and other specific properties of materials. However, the phase field models of martensitic transformation do not have enough degrees of freedom to increase the interfacial width without changing the interfacial energy. In the present study, a scalable phase field model by introducing a global modified function is constructed to study MT, the modified function takes into account the inhomogeneous nature of order parameter gradient across the interfacial region. Through adjusting the free energy density and gradient coefficient, meanwhile keeping the interfacial energy density unchanged, the interfacial width and system size are increased, yet the MT feature can be fully characterized. The simulation results show that the modified phase field model can well solve the drawbacks such as fast growth rate of martensite, artificial orientation relationship between the variants of martensite, and disordered martensite microstructure in large scale system.

Keywords： martensitic transformation ; order parameter gradient ; interfacial width ; phase field method

PDF (8188KB) 元数据多维度评价相关文章收藏文章

本文引用格式导出 EndNote Ris Bibtex

魏铖, 柯常波, 马海涛, 张新平. 基于序参量梯度的改进相场模型及对大尺度体系马氏体相变的模拟研究[J]. 金属学报(中文版), 2018, 54(8): 1204-1214 https://doi.org/10.11900/0412.1961.2017.00478

WEI Cheng, KE Changbo, MA Haitao, ZHANG Xinping. A Modified Phase Field Model Based on Order Parameter Gradient and Simulation of Martensitic Transformation in Large Scale System[J]. Acta Metallurgica Sinica, 2018, 54(8): 1204-1214 https://doi.org/10.11900/0412.1961.2017.00478

预测和调控材料的微观组织进而设计和制备不同工程应用所需的材料已成为材料科学与工程研究的一个热点。马氏体相变作为多种材料中一种重要的结构转变,对其微观组织的实验研究已获得许多进展^[1,2,3,4,5]。但由于马氏体相变速率快,实验难以准确捕捉随时间而变化的组织学和动力学信息;加之马氏体相变受多种复杂因素的影响,包括温度场、应力场和磁场等,这对实验研究(如原位观察)提出了严苛的要求,且实验难以准确反映单因素的影响。目前计算机模拟已成为一种有效的替代手段。构建合理高效的数值模型,对马氏体相变组织学和动力学进行模拟研究,是对实验研究的有效补充,具有重要的理论意义和应用价值。

相场法是目前模拟微观组织演化最有效的数值分析方法之一,它能够耦合温度场、应力场、磁场以及其它外部因素对微观组织的影响,且无需直接对具有复杂形貌的界面进行追踪^[6]。Wang和Khachaturyan^[7]基于微观弹性理论^[8],开创性地建立了适用于马氏体相变的相场模型,成功预测了马氏体形核和长大过程。随后,运用相场模型对马氏体相变的模拟研究逐步发展起来,包括建立了超弹性^[9,10]、弹塑性^[11,12]的相场模型,分析了惯性力^[13]、自协调效应^[14]、孔隙^[15]、非均匀形核^[16]等因素对马氏体相变的影响。然而,目前对马氏体相变的相场模拟研究仍局限于较小尺度(模拟体系格点尺度为纳米级),主要原因是相场模型中与实际相变体系中界面能相关联的界面宽度往往较小;同时,模拟过程需要在界面区域划分足够的格点数以保证相场界面的扩散性,使得相场模型对大尺度马氏体相变的模拟需要极大的计算量。已有的增大相场模型模拟尺度的方法主要基于实际体系中的界面能密度^[17],通过调整相场模型中的体积化学能系数与梯度能系数,进而增大界面宽度并扩大相场模拟尺度。由于在马氏体相变相场模型中体积化学能系数与梯度能系数由体系属性所决定,难以直接通过调整系数的方式增大模型界面宽度。相场模型的重要特征是界面具有扩散特性,通过序参量在数值上的梯度进行实现;因此,从序参量梯度入手是一种值得尝试的方法。

本工作基于序参量梯度值在界面区域远大于非界面区域的特性,构造出与序参量梯度相关的全局修正函数,通过调整界面区域的体积化学能系数与梯度能系数,实现对大尺度体系下的马氏体相变的高效率模拟。具体研究内容主要包括：(1) 分析采用全局修正函数后马氏体新相的生长速率;(2) 研究全局修正函数对变体间位向关系形成的影响;(3) 探讨体系大小与全局修正函数对NiTi合金B2-R相变微观组织形貌的影响。

1 相场模型

1.1 模型构造

马氏体相变相场模型中决定相变演化路径的能量包括体积化学自由能、梯度能和弹性应变能。其中,体积化学自由能促进马氏体相变;梯度能与体积化学自由能在界面区域的变化量共同组成界面能,控制界面移动;弹性应变能影响马氏体的形核取向及生长方向。

体积化学自由能密度φ_chem通常由Landau多项式^[18]表示：

$φ_{chem} = Δ f [{\frac{A}{2} \overset{n}{\sum_{p = 1}} η_{p}^{2} - \frac{B}{3} \overset{n}{\sum_{p = 1}} η_{p}^{3} + \frac{C}{4} (\overset{n}{\sum_{p = 1}} η_{p}^{2})}^{2}]$ (1)

式中,Δf为母相与马氏体相的体积自由能密度差;η_p (p=1, 2, $\dots$ , n)为结构序参量,p表示新相的第p个变体,n为马氏体变体的组数;η₁=η₂= $\dots$ =η_n=0表示体系此处为母相所占用,η_p=η₀ (p=1, 2, $\dots$ , n)且η_i=0 (i=1, 2, $\dots$ , n,i≠p)表示体系此处对应为马氏体相的第p种变体;η₀为体系表征马氏体相的平衡值,通常为1;A、B、C为体系的膨胀系数,一般取正数且满足以下关系式^[19]：

$\{\begin{array}{l} \frac{A}{2} - \frac{B}{3} + \frac{C}{4} = - 1 \\ A - B + C = 0 \end{array}$ (2)

体系的梯度能密度φ_grad通常用梯度项表示：

$φ_{grad} = \overset{n}{\sum_{p = 1}} \frac{1}{2} k_{η} {(\nabla η_{p})}^{2}$ (3)

式中,k_η为序参量场的梯度能系数,由界面能密度确定。

界面能由体积化学能的变化量和体系梯度能组成,界面能密度γ可表示为^[20]：

$γ = \int_{l_{1}}^{l_{2}} \{Δ φ_{chem} + φ_{grad}\} d l =$

$\int_{l_{1}}^{l_{2}} \{Δ f ∙ θ (η_{1}, η_{2}, \dots, η_{n}) + \frac{k_{η}}{2} \overset{n}{\sum_{p = 1}} {(\nabla η_{p})}^{2}\} d l$ (4)

式中,l₁和l₂分别为界面两侧的位置;l为界面的法向坐标;Δφ_chem为界面区域体积化学能密度的变化量,其值为体积化学能密度φ_chem与马氏体相具有的体积化学能密度之差;因此,Landau多项式θ (η₁, η₂, $\dots$ , η_n)为：

$θ (η_{1}, η_{2}, \dots, η_{n} {) = \frac{A}{2} \overset{n}{\sum_{p = 1}} η_{p}^{2} - \frac{B}{3} \overset{n}{\sum_{p = 1}} η_{p}^{3} + \frac{C}{4} (\overset{n}{\sum_{p = 1}} η_{p}^{2})}^{2} + 1$ (5)

当界面宽度趋于稳定时,梯度能密度与体积化学能密度的变化量相等,界面能密度可采用下式^[17]计算：

式中,a_p和b_p分别为界面两侧序参量η_p的值,通常取0或1。参数a_p、b_p、Δf和θ (η₁, η₂, $\dots$ , η_n)均为常量,由实际体系中的界面能密度γ可直接确定k_η。

在马氏体相变中,弹性应变能G_elast主要包括新相和母相间由于晶格畸变产生的弹性应变能及马氏体变体之间自协调而产生的能量,可表示为：

$G_{elast} = \int_{V} \frac{1}{2} C_{ijkl} ε_{ij}^{el} ε_{kl}^{el} d V$ (7)

式中,V为所研究体系的体积;C_ijkl表示弹性模量张量(i, j, k, l为空间向量); $ε_{ij}^{el}$ 和 $ε_{kl}^{el}$ 是弹性应变张量,为总应变与相变应变之差。其中,相变应变 $ε_{ij}^{t}$ 通过下式^[18]求得：

$ε_{ij}^{t} = \overset{n}{\sum_{p = 1}} η_{p} ε_{ij}^{0} (p)$ (8)

式中,应变张量 $ε_{ij}^{0} (p)$ 描述的是母相转变为第p个马氏体变体时对应的无应力相变应变。

将马氏体相变中的体积化学能、梯度能和弹性应变能构成的Ginzburg-Landau函数代入非守恒序参量场的Allen-Cahn方程^[21]‍,得到以下控制方程：

$\frac{\partial η_{p}}{\partial t} = - L [Δ f (A η_{p} - B η_{p}^{2} + C η_{p} \overset{n}{\sum_{q = 1}} η_{q}^{2}) -$

$k_{η} Δ η_{p} - σ_{ij} ε_{ij}^{0} (p)] + ξ_{p} (r, t)$ (9)

式中,t为体系组织演化的时间;L为相场动力学系数;Δ为Laplace算子; $ξ_{p} (r, t)$ 为序参量对应的Langevin噪声,用以模拟形核阶段的热起伏,一般假定Langevin噪声项与时间和空间不相关,且满足以下条件：

$< ξ_{p} (r, t) > = 0$ (10)

$〈ξ_{p} (r, t) ξ_{q} (r', t')〉 = - 2 k_{B} TLδ (r - r') δ (t - t')$ (11)

式中, $〈〉$ 表示取均值, $〈ξ_{p} (r, t)〉$ 和 $〈ξ_{q} (r', t')〉$ 分别为序参量场热起伏对空间和时间的平均;r、r′和t、t′分别为体系中对应的不同位置和不同时间;k_B为Boltzmann常数,T为热力学温度,δ为Kronecker Delta函数。Langevin噪声可在每一时间步采用满足标准Gauss分布的随机数来实现。

1.2 模型改进

相场模型中的界面宽度λ、母相与马氏体相的体积自由能密度差Δf及梯度能系数k_η关系可表示为^[17]：

其中,b和a为界面两侧序参量值,通常分别取1和0;Δf和k_η均由体系属性所决定,难以直接调整^[22]。界面宽度λ内通常需要3~5个格点以保证模拟的准确性。若采用的格点间距d大于λ,则界面宽度将因格点限制被增大至λ′,致使原有的界面平衡被破坏。根据式(4),若λ′=Nλ时,积分长度相应地也增大N倍,体积化学能密度变化量的积分 $\int_{l_{1}}^{l_{2}} Δ φ_{chem} d l$ 相应地增至 $N \int_{l_{1}}^{l_{2}} Δ φ_{chem} d l$ 。同时,由于界面宽度被增大,序参量梯度 $\nabla η_{p}$ 降至 $\frac{\nabla η_{p}}{N}$ ,梯度能密度 $\frac{k_{η}}{2} \overset{n}{\sum_{p = 1}} {(\nabla η_{p})}^{2}$ 减小至 $\frac{k_{η}}{2 N^{2}} \overset{n}{\sum_{p = 1}} {(\nabla η_{p})}^{2}$ ,梯度能密度的积分 $\int_{l_{1}}^{l_{2}} \frac{k_{η}}{2 N^{2}} \overset{n}{\sum_{p = 1}} {(\nabla η_{p})}^{2} d l$ 变为 $\frac{1}{N} \int_{l_{1}}^{l_{2}} \frac{k_{η}}{2} \overset{n}{\sum_{p = 1}} {(\nabla η_{p})}^{2} d l$ 。原有平衡中,梯度能密度等于体积化学能密度的变化量。因此,若将原宽度的界面能密度记为γ,则计算的界面能密度(γ $'$ )为：

$γ' = (N + \frac{1}{N}) \frac{γ}{2} > γ, N \in (1, + \infty)$ (13)

由此可见,当格点间距d大于λ时会高估界面能密度,影响相变演化路径。由于在大尺度体系下加密格点需要极大的计算量,只有选取合适的方法增大界面宽度,才能实现高效模拟。

相场模型的模拟区域可分为界面区域与非界面区域(马氏体相或母相区域)。非界面区域序参量为0或者1,其梯度为0;而在界面区域中序参量沿界面法向有显著变化,其梯度较大。因此,根据序参量的梯度值,能够很好区分界面区域与非界面区域。基于此特性构造序参量梯度的全局修正函数g_p;在不改变界面能密度的前提下,仅针对界面区域调整参数Δf和k_η,在增大界面宽度进而扩大模拟体系尺度的同时,尽量避免了参数调整影响到基于界面能密度的体系属性。全局修正函数g_p的表达式为：

$g_{p} (υ_{p}) = (\frac{λ_{m}}{λ} - 1) ∙ υ_{p}^{2} (3 υ_{p}^{2} - 8 υ_{p} + 6) + 1, \frac{λ_{m}}{λ} \in (1, + \infty)$ (14)

其中,λ为理论界面宽度,λ_m为调整后的界面宽度;ν_p为等效梯度：

$υ_{p} = λ_{m} |\frac{\partial η_{p}}{\partial l}|$ (15)

式中, $\frac{\partial η_{p}}{\partial l}$ 为第p个变体的序参量在界面法向的梯度。g_p满足以下关系式：

$\{\begin{array}{l} g_{p} (0) = 1; g_{p} (1) = \frac{λ_{m}}{λ} \\ g'_{p} (0) = 0; g'_{p} (1) = 0 \end{array}$ (16)

图1为λ_m=10、λ=1和λ_m=8、λ=1时,g_p分别随ν_p变化的曲线。由图1和式(16)可见,当ν_p=0时,g_p=1,表示在非界面区域不进行调整;当ν_p≈1时,g_p≈λ_m/λ(如图1所示,当ν_p为1时,点A和点B对应的g_p分别为10和8),表示在界面区域对参数Δf和k_η进行相应调整可以增大界面宽度。

显示原图| 下载原图ZIP| 生成PPT

图1 全局修正函数g_p随等效梯度ν_p的变化曲线(λ_m/λ=8和λ_m/λ=10)

Fig.1 Global modification function (g_p) vs equivalent gradient (ν_p) at different width ratios (λ_m/λ), in which g_p does not adjust any coefficients in the non-interface region (ν_p=0, g_p=1), but modifies the coefficients in the interfacial region (see points A and B)

利用g_p对Δf和k_η进行调整：

$\{\begin{array}{l} Δ f_{p} = \frac{Δ f}{g_{p}} \\ k_{p} = k_{η} ∙ g_{p} \end{array}$ (17)

式中,Δf_p为第p个变体调整后的体积自由能密度差;k_p为第p个变体调整后的梯度能系数。根据式(12)和(17),同时调整Δf和k_η后,界面能密度γ保持不变,界面宽度λ则增至g_pλ。需要指出的是,若调整后的Δf_p/k_p过小,由于数值计算的不稳定性,可能无法形成界面,界面宽度的增大在不同实际体系中有相应的限度。

选取Ti-50.5Ni (原子分数,%)作为模拟研究对象,研究其中B2到R相的相变过程。B2相属于高对称立方晶系;R相属于菱方晶系,有4个变体;在B2相参考坐标系下4个变体的相变应变张量分别为^[23]：

$\{\begin{array}{l} ε_{ij}^{0} (1) = [\begin{matrix} 0 & 0.0048 & 0.0048 \\ 0.0048 & 0 & 0.0048 \\ 0.0048 & 0.0048 & 0 \end{matrix}] \\ ε_{ij}^{0} (2) = [\begin{matrix} 0 & - 0.0048 & - 0.0048 \\ - 0.0048 & 0 & 0.0048 \\ - 0.0048 & 0.0048 & 0 \end{matrix}] \\ ε_{ij}^{0} (3) = [\begin{matrix} 0 & 0.0048 & - 0.0048 \\ 0.0048 & 0 & - 0.0048 \\ - 0.0048 & - 0.0048 & 0 \end{matrix}] \\ ε_{ij}^{0} (4) = [\begin{matrix} 0 & - 0.0048 & 0.0048 \\ - 0.0048 & 0 & - 0.0048 \\ 0.0048 & - 0.0048 & 0 \end{matrix}] \end{array}$ (18)

模拟采用的温度为245 K;假设B2相与R相具有相同的弹性系数,C₁₁=162 GPa,C₁₂=129 GPa,C₄₄=34 GPa^[24]。相场方程中的各参数选取如下^[25]：L=5.4774 m²/(Ns);A=0.14,B=12.42,C=12.28;Δf=1.3×10⁷ J/m³;k_η=1.1575×10^-11 N。未采用全局修正函数所计算得到的相场模型界面宽度λ=0.67 nm。

2 模拟结果与分析

2.1 生长速率

采用3种相同网格数(800 $\times$ 800)而不同尺度的二维体系,模拟体系大小通过增大格点间距进行改变(在模拟计算数值稳定的范围内调节),分别为：134 nm $\times$ 134 nm、1072 nm $\times$ 1072 nm和4288 nm $\times$ 4288 nm,对应的格点间距分别为0.17 nm (0.25λ)、1.34 nm (2λ)和5.36 nm (8λ)。初始状态均为在模拟区域中心放置相同大小(10.72 nm $\times$ 10.72 nm)的R相核心,R相的相变应变由式(18)中的 $ε_{ij}^{0} (1)$ 给出。

图2和3分别为修正前和修正后新相在不同尺度体系中的演化过程及新相的等效半径随相变时间变化(等效半径r =$\sqrt{s/\pi}$,s为新相面积)的模拟结果,其中红色和蓝色分别代表R相和B2相。由于各尺度下相同的物理属性,新相应呈现相同的生长速率;但由图2可见,新相生长速率却随体系尺度的增大而增快。这是由于界面宽度随着格点间距增大,原有的平衡遭到破坏,由式(13)可知,界面能密度增大致使界面区域序参量变化速率加快,新相的演化速率显著增快。相比加密格点需增加极大的计算量,图3模拟时采用相同格点数,引入全局修正函数g_p,参数λ_m分别选取5.36 nm和21.44 nm,采用式(17)同时调整k_η和Δf,由式(6)和(12)可知界面宽度增大后,界面能密度γ保持不变;从图3所示引入全局修正函数后2个较大尺度体系的模拟结果看,新相在2个不同尺度体系中呈现基本相同的演化速率。

显示原图| 下载原图ZIP| 生成PPT

图2 新相在不同尺度体系中的演化过程和等效半径变化

Fig.2 Morphological evolution of new phase (red color) simulated in systems with scales of 134 nm×134 nm (a1~a3), 1072 nm×1072 nm (b1~b3) and 4288 nm×4288 nm (c1~c3) at time of t=0.00 s (a1~c1), t=0.15 s (a2~c2) and t=0.45 s (a3~c3), and change of the equivalent radius with time (In which each simulation system owns the same initial square nucleus size of 10.72 nm×10.72 nm) (d)

显示原图| 下载原图ZIP| 生成PPT

图3 引入全局修正函数后新相在不同尺度体系中的演化过程和等效半径变化

Fig.3 Morphological evolution of new phase (red color) simulated with employing the global modification function in systems with scales of 1072 nm×1072 nm (a1~a3) and 4288 nm×4288 nm (b1~b3) at time of t=0.00 s (a1, b1), t=0.15 s (a2, b2) and t=0.45 s (a3, b3), and change of the equivalent radius with time (In which each simulation system owns the same initial square nucleus size of 10.72 nm×10.72 nm) (c)

2.2 位向关系

采用3种具有相同网格数(600 $\times$ 400)而不同尺度的二维体系(即101 nm $\times$ 67 nm、2010 nm $\times$ 1340 nm和20100 nm $\times$ 13400 nm)进行模拟,对应的格点间距分别为0.17 nm (0.25λ)、3.35 nm (5λ)和33.5 nm (50λ);其中对后2种较大尺度体系在引入全局修正函数前后进行对比。模拟选取B2-R相变作为研究对象。图4和5分别是引入全局修正函数前后的模拟结果。其中,所有体系均等分成左侧R相的变体1(蓝色)和右侧变体2 (红色),2变体的相变应变分别如式(18)中 $ε_{ij}^{0} (1)$ 和 $ε_{ij}^{0} (3)$ 所示,在二维面上的位向关系^[26]为(1 1)。对比图4b~e可见,演化初期弹性应变能逐渐降低,且可促进变体协调配置。随着演化的进行,小尺寸体系中弹性应变能持续降低,不断促使变体协调配置,最终弹性应变能接近于0,体系形成合理的位向关系。但对大尺寸体系,由式(13)可知,由于体系原有的平衡被破坏而使变体间界面区域限制界面移动的体积化学能增大,促进界面移动的梯度能减小,界面能对界面移动起限制作用。当界面能驱动力大于应变能驱动力时,受到界面能的限制,应变能难以促进变体协调配置,变体间无法形成合理的位向关系。

显示原图| 下载原图ZIP| 生成PPT

图4 2个R相变体在不同尺度体系中的位向关系及应变能密度变化

Fig.4 Morphological evolution of two R-phase variants in their initial states (variant 1 in blue and variant 2 in red) (a) simulated in systems with different scales of 101 nm×67 nm (b), 2010 nm×1340 nm (c) and 20100 nm×13400 nm (d), and change of strain energy density with time (e)

从图5所示引入全局修正函数(λ_m分别为13.4和134 nm)后2个大尺度体系的模拟结果可见,采用式(17)同时调整梯度能系数k_η和体积自由能密度差Δf后,变体间界面区域的体积化学能与梯度能重新处于平衡状态,界面能不再限制变体间的界面移动,弹性应变能得以在演化过程中持续降低,促使变体协调配置,弹性应变能最终接近于0,体系呈现合理的位向关系。

显示原图| 下载原图ZIP| 生成PPT

图5 修正函数后2个R相变体在不同尺度体系中的位向关系及应变能密度变化

Fig.5 Morphologicl evolution of two R-phase variants (variant 1 in blue and variant 2 in red) simulated with employing the global modification function in systems with scales of 2010 nm×1340 nm (a) and 20100 nm×13400 nm (b), and change of strain energy density with time (c)

2.3 B2-R相变的微观组织形貌

采用具有相同网格数(600 $\times$ 600)的3种不同尺度体系(即101 nm $\times$ 101 nm、2010 nm $\times$ 2010 nm和20100 nm $\times$ 20100 nm)进行模拟,对应的格点间距分别为0.17 nm (0.25λ)、3.35 nm (5λ)和33.5 nm (50λ)。其中,对后2种较大尺度体系在引入全局修正函数前后进行对比。R相的形核过程采用符合Gauss分布的随机噪声项模拟,位移场和相场均采用周期性边界条件。

理论上,R相4个变体在B2相所在参考坐标系的二维模拟中所形成的组织^[27]具有如图6所示特征。

显示原图| 下载原图ZIP| 生成PPT

图6 R相4个变体分别在(001)面上形成2种界面示意图

Fig.6 Schematic configuration of four R-phase variants in the (001) plane with variants 1 and 4 (a), and variants 1 and 2 (b)

图7和8分别是引入全局修正函数前后B2-R相变在不同尺度体系下的演化过程,其中蓝色代表母相B2相,浅蓝色代表R相变体1,绿色代表R相变体2,黄色代表R相变体3,红色代表R相变体4。由图7a1~a3和d可见,在随机噪声的形核阶段R相各变体都以均匀方式形核,由于变体相互挤压,体系的弹性应变能增大;随相变进行,小尺度体系中R相变体以逐渐清晰且有序的结构形成“鱼骨”状的二维组织,体系的弹性应变能降低,模拟结果与实验观察结果^[27]符合得很好。当体系尺度增大后,如图7b1~b3和c1~c3所示,由于界面宽度被强制增大,致使体积化学能与梯度能在界面区域分布不合理,阻碍弹性应变能促进变体协调配置,结果R相变体无法进一步形成具有清晰位向关系的自协调组织。

显示原图| 下载原图ZIP| 生成PPT

图7 B2-R相变在不同尺寸体系中的演化过程及应变能密度变化

Fig.7 Morphological evolution of B2-R phase transformation simulated in systems with scales of 101 nm×101 nm (a1~a3), 2010 nm×2010 nm (b1~b3) and 20100 nm×20100 nm (c1~c3) at time of t=0.05 s (a1~c1), t=0.50 s (a2~c2), t=40.10 s (a3), t=16.52 s (b3) and t=16.13 s (c3), and change of strain energy density with time (d) (In which the blue region stands for B2-phase, the light blue region stands for variant 1, the green region stands for variant 2, the yellow region stands for variant 3, and the red region stands for variant 4)

由图8所示2种较大尺度体系中引入全局修正函数(λ_m分别取13.4和134 nm)后的模拟结果可见,采用式(17)同时调整梯度能系数k_η和体积自由能密度差Δf对界面区域的体积化学能与梯度能进行修正后,R相得以形成具有清晰位向关系的结构,与实验观察到的组织形貌特征符合较好。

显示原图| 下载原图ZIP| 生成PPT

图8 修正后B2-R相变在不同尺寸体系中的演化过程及应变能密度变化

Fig.8 Morphological evolution of B2-R phase transformation simulated with employing the global modification function in systems with scales of 2010 nm×2010 nm (a1~a3) and 20100 nm×20100 nm (b1~b3) at time of t=0.05 s (a1, b1), t=0.50 s (a2, b2), t=67.71 s (a3), t=83.76 s (b3), and change of strain energy density with time (c) (In which the blue region stands for B2-phase, the light blue region stands for variant 1, the green region stands for variant 2, the yellow region stands for variant 3, and the red region stands for variant 4)

3 结论

(1) 根据相场模型中界面与非界面区域序参量变化的差异,构造了能有效区分2种区域的全局修正函数,并以此改进了马氏体相变的相场模型;改进的模型可在界面区域调整体积自由能密度差Δf和梯度能系数k_η,避免了影响体系整体属性的同时相当于增大了原相场模型中的界面宽度,进而扩大了马氏体相变的模拟尺度。

(2) 相比原有的马氏体相变相场模型,运用改进后的模型在相同的网格数下能够准确分析大尺度体系的界面能密度,可避免相同物理属性的新相随模拟尺度的增大而生长速率加快的问题。

(3) 应用改进的相场模型在不同尺度下对B2-R相变进行模拟,R相变体都得以形成合理的位向关系并呈现清晰有序的“鱼骨状”微观组织,与实验结果相符。

The authors have declared that no competing interests exist.

〈

〉