基于多尺度扩展有限元法(XFEM)建模揭示显微组织调控的MoSiBTiC合金裂纹扩展与断裂韧性(Fracture Resistance)

《Advanced Engineering Materials》：Microstructure-Controlled Crack Propagation and Fracture Resistance in MoSiBTiC Alloy Revealed by Multiscale Extended Finite Element Method Modeling

　　

鉴于显微组织( microstructure )对MoSiBTiC合金断裂韧性(fracture resistance)的影响，研究人员采用扩展有限元法(Extended Finite Element Method, XFEM)研究该合金的裂纹扩展行为与局部断裂响应，进而评估显微组织形貌对断裂韧性的影响。研究人员建立有限元分析(Finite Element Analysis, FEA)模型表征均匀化(homogenized)合金的宏观断裂响应，并标定有效的断裂相关参数。结合纳米压痕(nanoindentation)与反分析(inverse analysis)，确定各组成相的弹塑性性质，并利用简化局部断裂模型进一步标定其平均有效内聚(cohesive)参数。此外，基于四点弯曲(four-point bending)试验建立多尺度FEA模型，将微裂纹路径与试验结果进行对比。研究发现，该合金的断裂过程涉及多种机制的交互作用。基于界面理想结合的局部断裂模拟结果表明，断裂韧性的变化不能仅用延性相体积分数来解释。相反，促使应力分布更对称的颗粒排布、圆度(circularity)更高及长宽比(aspect ratio)更低的规则颗粒，往往与合金较高的断裂韧性相关联。

本文对发表于《Advanced Engineering Materials》的论文《Microstructure-Controlled Crack Propagation and Fracture Resistance in MoSiBTiC Alloy Revealed by Multiscale Extended Finite Element Method Modeling》进行解读。

研究背景与意义

钼硅硼钛碳(MoSiBTiC)合金是有望替代镍基高温合金、用于下一代航空航天推进系统超高温结构材料的有力候选者，但其室温断裂韧性低是制约工程应用的关键障碍。虽然添加TiC可提高断裂抗力，单纯增加延性Mo固溶体(Mo_ss)相体积分数虽能提高韧性却会牺牲抗氧化性与蠕变强度。已有研究表明显微组织(microstructure)几何形貌对断裂韧性有重要影响但尚未被定量阐明。由于多相材料各相间力学性能差异大，传统有限元法(FEM)难以模拟复杂微观组织中不预设路径的裂纹萌生与扩展。扩展有限元法(XFEM, Extended Finite Element Method)允许裂纹独立于网格拓扑发起和扩展，适合研究多相显微组织中的裂纹演化。因此，研究人员开展本研究，通过将实验观测的显微组织纳入连续多尺度XFEM框架，量化硬质相颗粒几何特征(角度、圆度、长宽比、尺寸)对裂纹偏转、能量耗散及整体断裂韧性的影响，为超高温结构材料的显微组织导向设计提供依据。

主要关键技术方法

研究人员以名义成分为65.0Mo–5.0Si–10.0B–10.0Ti–10.0C(at%)的第一代MoSiBTiC合金为对象，经电弧熔炼和1800°C/24h氩气氛热处理后开展三类实验与模拟：(1)采用带Chevron缺口试样进行室温四点弯曲(four-point bending)试验测定合金宏观条件断裂韧性(K_Q)及载荷–位移曲线；(2)对各组成相(Mo_ss、Mo₅SiB₂即T₂相、TiC、Mo₂C)进行纳米压痕(nanoindentation)试验并结合Oliver–Pharr方法及反分析(inverse analysis)确定弹性模量与Mo_ss相的双线性硬化本构参数；(3)使用LS–DYNA软件分别建立三点模型进行XFEM模拟——均匀化四点弯曲模型标定整体有效内聚(cohesive)参数；含真实BSE–SEM提取显微组织的简化局部断裂(SLF, Simplified Local Fracture)模型在校准各相内聚参数(T₂、TiC、Mo₂C基于Griffith–Irwin关系由断裂韧性得峰值牵引力σ_max，Mo_ss由SLF早期响应拟合)后分析不同微区裂纹扩展；嵌入直缺口自由表面显微组织的多尺度四点弯曲模型模拟裂纹扩展并与BSE观测对比；采用ImageJ定量统计硬质相颗粒的角度(orientation angle)、圆度(circularity = 4πA/P²)、长宽比(aspect ratio = 长轴/短轴)及面积，分析其与临界能量释放率(G_c)的关系。

研究结果

5.1 Macroscopic Crack Propagation Behavior(宏观裂纹扩展行为)

研究人员通过均匀化合金的四点弯曲XFEM模拟发现，最大主应力(Max. Prin. S)在预制裂纹尖端集中形成对称蝶形应力场，峰值荷载时应力达约700 MPa(即均匀化合金内聚峰值牵引力σ_max)，超过后引发损伤及裂纹扩展。模拟宏观裂纹路径与实验试样表面观测的总体扩展趋势吻合，证实宏观均质模型可捕捉裂纹启裂及早期扩展的主要响应，但因无法描述微观多机制损伤导致荷载跌落段存在差异，需从微观角度进一步分析。

5.2 Microscopic Crack Propagation Behavior(微观裂纹扩展行为)

在多尺度模型中嵌入真实显微组织后发现，启裂前脆性相沿裂纹扩展方向形成拉应力集中带。模拟显示初始阶段裂纹以穿晶(transgranular)方式穿过T₂和TiC等脆性相，与BSE观测相符；扩展后期观察到Mo_ss/T₂界面脱粘(interfacial debonding)、裂纹在Mo_ss中止(crack arrest)并在前方再萌生(renucleation)及二次裂纹，Mo_ss形成的韧带桥接(ligament bridging)是增韧因素之一。由于模型未引入界面单元故未完全再现界面脱粘过程，但表明Mo_ss/T₂界面易成应力集中弱区。实际三维裂纹演化导致模拟与表面BSE图像中裂纹轻微偏转存在细微差别。

5.3 Comparison of Crack Propagation and Fracture Properties among Different Microstructures(不同显微组织间裂纹扩展与断裂性能对比)

研究人员选取6个100×100 μm代表性显微区域(Microstr 1–6)进行SLF模拟，计算单位厚度临界能量释放率G_c(反映断裂抗力)与裂纹路径长度。结果显示Microstr 4–6的G_c(138–164 mJ·m^?2)高于Microstr 1–3(78–108 mJ·m^?2)，且前者裂纹尖端最大主应力分布保持较对称，后者随拉伸位移增加应力场变得不对称。较对称应力场使裂尖前方的载荷被更大范围材料分担，利于稳定扩展；不对称应力场促使裂纹沿微观弱路径偏转加速失稳。Mo_ss相体积分数与G_c无显著相关性，说明在本研究范围内断裂抗力变化不能仅用延性相含量解释，显微组织形貌影响更直接，其中硬质T₂相颗粒取向排列与裂尖应力对称性有关。

5.4 Characteristics and Quantitative Analysis of Microstructural Morphology(显微组织形貌特征与定量分析)

用ImageJ定量统计硬质相(T₂、TiC、Mo₂C)颗粒形貌参数发现：发生开裂的颗粒较全体颗粒总体呈现更低圆度(circularity)、更高长宽比(aspect ratio)及更大面积，说明裂纹倾向于穿过形状不规则(低圆度、高长宽比)和尺寸较大的颗粒。对沿裂纹路径开裂硬质颗粒的参数与G_c线性回归显示，G_c随圆度增大而升高(R²=0.62)，随长宽比减小而升高(R²=0.57)，即在理想结合界面假设下更规则(高圆度、低长宽比)硬质颗粒有利于提升断裂抗力。T₂相角度统计显示较高断裂抗力微区中T₂颗粒排布更易产生对称应力场(如沿裂纹扩展方向排列或更随机分布保持对称)。需注意上述定量关系基于各相各向同性及界面理想结合假定，实际中断裂还受界面结合强度、取向及晶体学各向异性影响。

讨论与结论翻译(Conclusions)

研究人员得出以下结论：

1. 1.
   基于四点弯曲试验与二维均匀化FEA分析，初步表征并标定MoSiBTiC合金宏观断裂响应与峰值强度；结合纳米压痕试验与反分析确定各组成相弹塑性参数，并利用简化局部断裂(SLF)模型标定各组成相的平均有效断裂相关参数，为多相显微组织构型的断裂模拟提供基础输入数据。
2. 2.
   多尺度XFEM模拟与BSE观测对比表明，MoSiBTiC合金复杂断裂行为无法仅用单一二维裂纹路径充分描述；裂纹扩展涉及多种机制交互作用，包括脆性相穿晶断裂、界面脱粘、二次裂纹、裂纹在Mo_ss中止及再萌生，这些复杂断裂机制的协同作用促成相对平稳的能量耗散过程与荷载渐降，延迟裂纹启裂后的失效。
3. 3.
   不同显微组织的局部断裂模拟表明，在本研究考察范围内，MoSiBTiC合金断裂韧性的变化不能仅用Mo_ss相体积分数解释；相比之下，显微组织形貌对裂纹扩展行为影响更直接。具较高断裂抗力的显微组织通常伴随裂尖附近较对称应力场，这与硬质颗粒特别是T₂相的取向排布有关。
4. 4.
   硬质相颗粒定量分析揭示，形貌更规则(较高圆度与较低长宽比)的颗粒有助于提升断裂抗力。需注意形态–断裂抗力定量关系是在各组成相各向同性及界面理想结合假设下获得，界面脱粘、界面取向及晶体学各向异性对裂纹扩展的影响有待进一步研究。