2.5D C/SiC复合材料表面CVD-SiC涂层的热震失效微-介-宏观跨尺度建模

《Journal of Materials Research and Technology》：Micro–Meso–Macro Modelling of Thermal-Shock Failure in CVD-SiC Coatings on 2.5D C/SiC Composites

　　

化学气相沉积（Chemical Vapor Deposition, CVD）制备的碳化硅（SiC）涂层在2.5维编织碳/碳化硅（C/SiC）复合材料表面表现出约一个数量级的热震寿命统计离散性，传统经验设计准则无法对此进行有效预测。研究人员提出了一种微-介-宏观集成计算材料工程（Integrated Computational Materials Engineering, ICME）框架，该框架耦合了原子尺度晶界（Grain Boundary, GB）退化模型、相场法生成的晶粒代表性体积单元（Representative Volume Element, RVE）以及包含双重内聚区（分别位于涂层内部及涂层/基体界面）的全局-局部介观模型，可在单次热震循环条件下定量预测涂层脱层行为。研究采用两层尺度桥接策略，将内聚性能分解为受缺陷控制的绝对基准值（以拉脱实验为锚定）和受化学键控制的归一化温度趋势（从分子动力学（Molecular Dynamics, MD）模拟传递），实现了物理严谨且可实验校准的映射关系。核心发现如下：（i）MD模拟显示，在1573 K温度下，晶界的黏聚强度保持率仅为32–38%，而晶粒内部的保持率约为65%；（ii）2.5D RVE分析揭示纤维束顶部存在约1.5倍的面内应力放大效应；（iii）在1400 °C至1500 °C之间发生了温度驱动的失效模式分岔，其机理符合He–Hutchinson裂纹偏转准则；（iv）对30个试样进行的扫描电子显微镜-能谱（Scanning Electron Microscopy–Energy Dispersive Spectroscopy, SEM–EDS）断口分析证实93.3%的病例为混合模式失效，交叉验证了预测的Γ≈0.25–0.5区间。该框架预测了界面粗糙度（Ra≈6–10 μm）、晶粒尺寸（dg≈1–10 μm）的非单调设计窗口以及最小涂层厚度（～70 μm），并在1300 °C下的单循环热震实验中得到了验证——中等粗糙度（Ra≈11.6 μm）试样的强度保持率为82%，远高于近抛光表面的47%。该框架为CVD-SiC及类似非均质涂层体系提供了可迁移的设计图谱。

研究背景与意义

碳纤维增强碳化硅（C/SiC）复合材料因其优异的比强度、比刚度和极端热机械载荷下的尺寸稳定性，已成为下一代航空航天热端结构（如航空发动机燃烧室衬套、喷管组件及液体火箭推力室）不可或缺的关键材料。然而，在氧化服役环境中，碳纤维相在约450 °C以上会发生快速氧化降解，严重限制了裸C/SiC复合材料的耐久性。为此，在1000–1500 °C的服役窗口内，必须施加环境障涂层（Environmental Barrier Coating, EBC），其中化学气相沉积（CVD）SiC外层因与SiC基体化学相容性好、热膨胀系数匹配（α ≈ 5 × 10^-6K^-1）且能形成保护性SiO₂钝化膜而被广泛采用。尽管具备这些优势，但在快速升降温循环中，热震诱导的脱层和剥落仍是限制可靠性的主要瓶颈。更为棘手的是，即使在名义相同的测试条件下，涂层寿命仍表现出约一个数量级的统计离散性，这暴露了现有经验试错设计准则预测能力的局限性。虽然现有文献将此归因于界面粗糙度、晶粒尺寸和初始缺陷分布的微观结构变异性，但目前缺乏能够将这些微观特征与宏观寿命定量关联的多尺度预测模型，这正是本研究开展的动因。此外，传统的涂层热震失效模型（如热障涂层TGO驱动模型或水汽诱导退化模型）不适用于该系统，因为CVD-SiC/C/SiC的失效机制既不涉及体积膨胀的热生长氧化物，也不以水汽侵蚀为主要驱动力，而是受热弹性失配应力、温度依赖的晶界弱化以及2.5D编织结构引起的应力非均匀性三者协同控制。

关键技术方法

研究人员建立了一个物理信息驱动的分层多尺度框架来预测热震诱导的降解。首先，利用分子动力学（MD）模拟，采用Tersoff多体势描述Si–C相互作用，计算了无缺陷单晶和三种代表性双晶晶界模型的温度依赖性降解动力学及牵引分离行为，提取了晶界和晶粒内部的归一化温度弱化函数。其次，通过相场法（Phase Field Method, PFM）生成了考虑CVD-SiC柱状织构的二维多晶微观结构，构建了晶粒尺度代表性体积单元（RVE），并通过双重尺度均匀化提取了涂层的各向异性有效力学性能。随后，在介观尺度建立了包含双重内聚区模型（Cohesive Zone Model, CZM）的有限元模型，分别表征涂层内部断裂和涂层/基体界面脱粘。为了解决几何复杂性带来的计算量挑战，研究采用了全局–局部嵌套子模型策略，即在全周期编织RVE中识别应力热点，再将边界条件映射到具有高分辨率正弦界面粗糙度的局部子模型中。最后，研究人员进行了单循环热震实验（最高1300 °C）和拉脱测试，测量残余黏结强度以校准和验证模型参数。

研究结果

3.1 原子基础：晶界退化主导局部失效

MD模拟揭示了晶界相较于晶粒内部具有显著且更严重的退化路径。在1573 K的实验验证温度下，晶界黏聚强度保持率仅为32–38%，而缺陷自由晶粒内部保持率约为65%，两者相差约30个百分点。这种差异为连续介质框架中将晶界弱化作为损伤局部化的主要物理驱动因素提供了定量依据。此外，敏感性分析表明，在1300 °C的验证条件下，连续介质预测结果稳健，但在更高温度下对原子间势函数的不确定性变得敏感。

3.2 宏观拓扑：编织几何引发局部应力放大

通过对2.5D编织RVE的分析，研究人员量化了编织拓扑对应力分布的影响。纤维束峰顶处的面内压应力达到4.88 GPa，相比谷底区域产生了约1.5倍的应力放大效应。这种由编织结构引起的空间应力梯度导致了寄生拉伸分量的产生，成为损伤萌动的相关驱动力。虽然弹性RVE预测涂层内主要为压应力，但实验观察到的微裂纹辅助松弛现象表明实际峰值应力应低于此上限估计。

3.3 介观机制：最优参数窗口

研究预测了热震耐受性与界面粗糙度之间的非单调依赖关系。当粗糙度极高（R_a=30.5 μm）时，界面发生灾难性脱粘；当粗糙度极低（R_a→0）时，损伤转变为贯穿厚度的涂层开裂；最佳性能出现在中间粗糙度范围（R_a≈6–10 μm）。这一预测得到了实验支持：未抛光组（R_a≈11.6 μm）保留了82%的初始黏结强度，而抛光组（R_a≈0.49 μm）仅保留了47%。此外，晶粒尺寸效应也表现出非单调性：极细晶粒（d_g?0.1 μm）因晶界面积分数过高导致有效韧性急剧下降；粗晶粒（d_g?30 μm）则因裂纹偏转效率降低而易开裂；最佳窗口为d_g≈1–10 μm。涂层厚度方面，研究发现存在一个最低完整性阈值，推荐h≥70 μm以避免薄涂层（h≤40 μm）中的能量释放模式转变为界面分层主导。

3.4 温度驱动的失效模式分岔

在固定几何参数下，随着热震温度升高，系统在1400 °C至1500 °C之间发生了明显的失效模式分岔。低于1400 °C时，损伤主要表现为局限于涂层厚度内的离散分段裂纹，界面基本完好；高于1500 °C时，分段裂纹到达界面后发生偏转，触发渐进式界面脱粘。这一转变符合He–Hutchinson裂纹偏转准则，即随着温度升高，界面断裂能相对于涂层体韧性的退化速度更快，导致韧性比低于偏转临界阈值。

讨论与结论总结

本研究建立的微-介-宏观集成计算材料工程（ICME）框架成功解决了CVD-SiC涂层热震失效预测中的尺度跨越难题。研究最重要的贡献在于提出了两层尺度桥接策略，将受缺陷控制的绝对基准值与受化学键控制的归一化温度趋势解耦，解决了原子尺度与连续介质尺度之间经典的“量级鸿沟”问题，无需进行任意拟合即可实现物理保真的参数映射。

基于该框架，研究人员明确了晶界差异弱化是层间损伤局部化的根本原因，并量化了2.5D编织结构导致的应力放大效应及其与界面形貌的协同作用。论文发表在《Journal of Materials Research and Technology》，研究给出的设计图谱表明，通过控制界面粗糙度在R_a≈6–10 μm、晶粒尺寸在d_g≈1–10 μm并确保涂层厚度大于70 μm，可以在单循环热震下获得高于90%的强度保持率。此外，研究通过SEM–EDS断口分析证实了93.3%的失效为混合模式，与He–Hutchinson准则预测的Γ≈0.25–0.5区间高度吻合，为模型的物理合理性提供了直接证据。这项工作不仅为CVD-SiC涂层提供了可迁移的设计工具，也为其他涉及原子退化与介观几何共同控制失效的非均质涂层系统（如EBCs、TBCs）提供了方法论模板。