基于TRIP钢应力-应变曲线的应力分配分析的反问题方法

《INTERNATIONAL JOURNAL OF MECHANICAL SCIENCES》：Inverse-problem approach for stress-partitioning analysis using stress–strain curve of TRIP steel

　　

研究人员提出了一种反问题方法，将双相复合材料的应力-应变（SS）曲线的计算模型拟合到亚稳态奥氏体不锈钢的SS曲线实验数据上，以估算面心立方γ（fcc-γ）相和变形诱导体心立方α’（bcc-α’）相各自的SS曲线，以及α’相体积分数（f(α′)）随应变增加的变化。该方法被应用于两种TRIP钢的SS曲线：Fe–18Cr–8Ni–0.1C（质量%）合金（0.1C钢）和Fe–18Cr–8Ni–0.1N（质量%）合金（0.1N钢）。研究人员估算了个相（γ和α’）的SS曲线，以再现文献中报道的TRIP钢整体SS曲线。对于0.1C钢，估算的α’相流变应力显著高于0.1N钢，表明碳添加在强化变形诱导α’相方面比氮添加更有效。此外，该方法使研究人员能够在拉伸变形过程中分析相应力（γ和α’相之间的应力分配）。进一步地，研究人员估算了拉伸变形过程中f(α′)的变化，结果与文献中的实验数据吻合良好。0.1C钢中f(α′)的增加率低于0.1N钢，这对于优化TRIP钢的强度-延性平衡至关重要。

论文解读文章

**研究背景与问题**
相变诱发塑性（TRIP）钢，尤其是亚稳态奥氏体不锈钢（如304钢），在变形过程中会经历变形诱导马氏体相变（DIMT），即面心立方γ（fcc-γ）相转变为体心立方α’（bcc-α’）相。DIMT不仅通过高强度α’相增强材料强度，还通过TRIP效应改善塑性，赋予材料优异的强度-延性平衡。然而，深入了解合金元素（如碳、氮）对TRIP钢力学性能的影响，需要定量表征变形过程中α’相体积分数（f(α′)）的变化以及各组成相（γ和α’）的变形行为。现有的实验方法，如X射线衍射、饱和磁化测量和中子衍射，虽能提供f(α′)和相应力分配信息，但依赖专用设备，成本高且设施有限。因此，研究人员提出结合实验应力-应变（SS）曲线和计算建模的实用方法，以系统研究合金元素效应，为优化TRIP钢性能提供更便捷的途径。

**本研究内容与结论**
本研究提出一种反问题方法，将基于割线法（secant method）的双相复合材料SS曲线计算模型（CM）拟合到TRIP钢实验SS曲线上，同时估算γ和α’相的SS曲线以及f(α′)随应变的变化。该方法应用于文献中报道的两种TRIP钢：Fe–18Cr–8Ni–0.1C（质量%）（0.1C钢）和Fe–18Cr–8Ni–0.1N（质量%）（0.1N钢）。估算结果表明：0.1C钢中α’相的流变应力显著高于0.1N钢，证实碳添加对强化变形诱导α’相的效果优于氮添加；0.1C钢中f(α′)的增加率低于0.1N钢，这对优化TRIP钢的强度-延性平衡至关重要。此外，该方法还能分析拉伸变形过程中的相应力（γ和α’相之间的应力分配），揭示0.1C钢中α’相应力远高于γ相，而0.1N钢未呈现此趋势。研究成果发表在《INTERNATIONAL JOURNAL OF MECHANICAL SCIENCES》。该反问题方法为从TRIP钢SS曲线中提取DIMT速率和相应力分配信息提供了实用途径，有助于深入理解合金元素对力学性能的影响。

**关键技术与方法**
本研究主要采用以下关键技术方法：
1. **基于割线法的正向计算模型**：假设γ相为基体相、α’相为球形夹杂相，通过该模型计算γ+α’双相复合材料的整体SS曲线（详见附录A.1–A.3）。
2. **Swift方程**：用于近似各相（γ和α’）的SS曲线，形式为σ^(r) = a^(r) (b^(r) + ε_p^(r))^N(r)，其中r=0（基体）或r=1（夹杂）。
3. **Matsumura方程**：描述α’相体积分数f(α′)随整体塑性应变ε_p^(C)的变化，形式为f(α′) = 1 – (1 + k q [ε_p^(C)]^q)^–1。
4. **反问题优化方法**：通过Nelder–Mead算法最小化成本函数J(x) = (1/M) Σ_i=1^M (σ_expt,i – σ_calc,i(x))²，以估算参数向量x = (a⁽⁰⁾, b⁽⁰⁾, N⁽⁰⁾, a⁽¹⁾, b⁽¹⁾, N⁽¹⁾, k, q)^T。实验数据来源于文献[5]中0.1C钢和0.1N钢的SS曲线（应变间隔0.01）。

**研究结果**
**1. 各相SS曲线及f(α′)的估算**
通过81组初始参数值进行反问题求解，成功估算出0.1C钢和0.1N钢的Swift方程参数（见表2）。图2显示，估算的各相SS曲线（α’和γ）及计算的整体SS曲线与文献实验数据（空心圆）高度吻合。0.1C钢中α’相的流变应力（约2103 MPa）远高于0.1N钢（约1528 MPa），而γ相流变应力差异不大，表明碳添加对α’相强化作用更强。f(α′)的估算结果（图2b, d）也与文献实验数据一致，0.1C钢的f(α′)增加率较低，这与其优异的强度-延性平衡相关（文献报道）。该反问题方法仅使用整体SS曲线数据，无需f(α′)实验数据，比X射线衍射等方法更简便实用。进一步对SUS301L钢在不同温度和应变率下的验证（附录A.4）表明该方法广泛适用。

**2. 相应力分配分析**
利用正向模型中的公式（附录A.3），计算了拉伸变形过程中γ相和α’相的相应力（σ*(0)和σ*(1)）。图3显示，0.1C钢中α’相应力（实心方块）显著高于γ相，且整体应力（空心方块）位于两者之间；0.1N钢则无此明显差异。这表明碳添加改变了相应力分配，进而影响TRIP钢的加工硬化行为，与文献结论一致。

**3. 参数估计的初始值依赖性**
成功估算率（J(x) < 5 MPa²）为0.1C钢42%、0.1N钢90%。初始参数值的有效性影响成功率。即使参数估计值（x^）因初始值不同而变化，f(α′)的预测误差（J_f）最大仅0.12（即误差在10%以内）。因此，建议从多组初始值出发，取参数平均作为最终估计。

**4. 夹杂相假设的影响**
分析显示，α’相形状（球形、扁球、长球）和相选择（α’或γ作夹杂）对整体SS曲线仅在低应变区（0.001–0.003）有轻微影响。由于DIMT过程中f(α′)从零开始增加，低应变区α’相体积分数极小，因此夹杂相假设的不确定性对本文方法有效性影响不大。

**讨论与结论**
**讨论部分总结**：研究指出，初始参数值的选择显著影响参数估计成功率，采用多组初始值并取平均是提高结果可靠性的有效策略。此外，夹杂相的形状和相选择在低应变区对SS曲线影响有限，且在实际变形初期f(α′)极小，因此假设α’相为球形夹杂是合理的，不会对估算结果产生显著干扰。

**结论部分翻译**：
1. 估算的各组成相（γ和α’）的SS曲线很好地再现了0.1C钢和0.1N钢整体SS曲线的实验数据。同时，估算的拉伸变形过程中f(α′)的变化与文献[5]中的实验数据吻合良好。模型参数估计的成功率显著依赖于初始估计的有效性，因此从多组初始参数值开始估计，并将估计参数值的平均值作为最终估计，是获得再现实验数据的计算结果的有效方法。
2. 0.1C钢和0.1N钢的α’相流变应力估计值存在显著差异，表明碳添加对变形诱导α’相的强化效果大于氮添加。此外，0.1C钢的f(α′)增加率低于0.1N钢，这对优化TRIP钢的强度-延性平衡至关重要。这些结果与文献[5]的发现一致。
3. 该反问题方法还能够分析拉伸变形过程中的相应力（γ和α’相之间的应力分配）。0.1C钢中α’相的应力显著高于γ相，但0.1N钢中未确认这一趋势。
该反问题方法提供了一种实用途径，从亚稳态奥氏体不锈钢的SS曲线实验数据中提取DIMT速率和γ与α’相之间的应力分配信息。因此，该方法加深了对合金元素影响TRIP钢力学性能的理解，为优化TRIP钢的强度-延性平衡铺平了道路。