通过虚拟裂纹闭合技术研究非均匀颗粒分布下弥散核燃料元件芯体-包壳界面开裂

《Mechanics of Materials》：Meat-cladding interface cracking in dispersion nuclear fuel elements under nonuniform particle distribution via the virtual crack closure technique

　　

在弥散核燃料元件中，燃料芯体与包壳之间的界面开裂会降低局部传热效率，从而威胁反应堆堆芯安全。该研究考虑界面附近颗粒的非均匀分布以及堆内热-机械耦合（thermomechanical coupling），对芯体-包壳界面的开裂行为进行了研究。研究人员采用X射线显微术（XRM）对类比弥散板内颗粒的空间分布进行了表征。随后建立了一个高保真计算模型，用于模拟芯体-包壳界面的开裂。基于该模型，研究人员利用虚拟裂纹闭合技术（VCCT）评估了不同微观结构构型下裂纹尖端的能量释放率。基于距离度量的点模式分析技术用于表征近界面区域的颗粒分布。这些方法同时也用于验证模拟模型的保真度。研究人员分析了近界面颗粒分布特征对界面能量释放率的影响。数值结果表明，颗粒团聚会加剧芯体-包壳界面的开裂行为。且更靠近弥散芯体中心的团聚区域与更高的界面能量释放率相关。该工作表明，芯体-包壳界面附近的微观结构可用于预测界面开裂与鼓泡位置，从而为评估弥散燃料制备工艺及提升反应堆堆芯安全性提供指导。

该论文发表于《Mechanics of Materials》，围绕弥散核燃料元件中芯体-包壳界面失效这一关键工程问题展开。弥散核燃料元件由弥散燃料芯体和金属包壳组成，燃料颗粒嵌入惰性金属基体中，兼具较高燃耗、较好导热性能和较高安全性等优势。然而，在现有制备过程中，燃料颗粒往往难以实现完全均匀分布，局部团聚现象较为常见。在堆内服役条件下，裂变气体可在弥散芯体与包壳之间的间隙中积聚，诱发内压升高，进而导致界面裂纹萌生、扩展并最终形成鼓泡失效。此类界面损伤不仅会削弱局部传热能力，还会进一步威胁反应堆堆芯安全。因此，阐明非均匀颗粒分布，特别是近界面区域颗粒团聚、颗粒尺寸和体积分数等特征，对芯体-包壳界面开裂的影响机制，具有明确的科学意义与工程价值。

从已有研究基础看，颗粒增强复合材料的断裂行为与颗粒空间分布高度相关。颗粒尺寸会影响复合材料的断裂韧性与流动应力；颗粒体积分数会改变材料整体抗裂能力；颗粒团聚则常常成为局部损伤和裂纹扩展的起始区域。对于弥散核燃料元件而言，芯体-包壳界面本质上属于双相材料界面，其裂纹尖端应力场具有振荡奇异性，这使传统基于界面应力强度因子的分析与数值实施更为复杂。相比之下，能量释放率G和J积分（J-integral）基于能量守恒原理，具有路径无关性，更适合表征此类异质界面的断裂驱动力。其中，能量释放率G可借助虚拟裂纹闭合技术（VCCT，一种基于裂尖结点力与裂纹张开位移求解断裂参数的方法）直接从有限元模型中获得，避免了复杂裂尖奇异单元构造及外推处理。因此，研究人员选择以界面能量释放率G作为表征芯体-包壳界面裂纹扩展驱动力的核心参数，并重点评估非均匀颗粒分布对其的影响。

研究人员首先开展了类比弥散板的颗粒分布统计分析，再据此构建高保真数值模型，用于模拟芯体-包壳界面裂纹在堆内热-机械耦合作用下的响应。文章明确将内部裂变气体压力视为界面裂纹萌生和扩展的重要载荷来源，并围绕四类关键颗粒分布特征展开分析，即燃料颗粒体积分数、颗粒尺寸、团聚程度以及团聚位置。通过将真实颗粒空间特征引入计算域，研究人员试图克服传统解析模型难以处理多颗粒非均匀分布叠加效应的问题，并建立一种能够反映实际微观组织的界面开裂评估方法。论文的核心结论是：近界面颗粒分布的不均匀性，尤其是颗粒团聚，会显著提高界面裂纹尖端的能量释放率；同时，若团聚区域位置更接近弥散芯体中心，则界面能量释放率更高，意味着界面更易发生开裂与鼓泡。这说明界面附近微观组织特征可作为预测界面损伤位置的重要依据，也为弥散燃料制备质量控制与反应堆安全评价提供了方法支撑。

就主要技术路线而言，研究人员采用X射线显微术（XRM）获取类比弥散板样品中颗粒空间分布，并利用数字图像处理开展统计表征；随后在ABAQUS 6.14平台上建立包含实际颗粒分布特征的高保真代表性体积单元模型；采用核密度估计、Voronoi图、K函数与g函数等点模式分析方法验证模型空间分布保真度；最后基于虚拟裂纹闭合技术（VCCT）计算芯体-包壳界面裂纹尖端能量释放率，并比较颗粒尺寸、体积分数、团聚程度与团聚位置等因素的影响。样本来源为采用轧制工艺制备的类比弥散板试样，其组成为ZrO₂颗粒、Zr合金基体和Zr合金包壳。

以下结合文中各部分内容对研究结果进行归纳解读。

Statistical analysis of particle distribution based on XRM
在这一部分，研究人员利用X射线显微术（XRM）对类比弥散板试样进行了三维或空间分布层面的颗粒表征，并借助数字图像处理技术对扫描结果进行了统计分析。该部分工作的目的在于获得颗粒真实分布的事实特征，为后续高保真界面开裂模型提供输入基础。论文强调，弥散板由ZrO₂颗粒、Zr合金基体和Zr合金包壳构成，且采用轧制工艺制备。由此可见，研究并非基于理想化随机颗粒分布直接建模，而是从实际样品出发提取颗粒空间特征，这对于提升模型与真实材料之间的一致性至关重要。该步骤奠定了后续分析颗粒非均匀性与界面断裂驱动力之间关系的基础。

Validation of the high-fidelity model
在高保真模型验证部分，研究人员在西安交通大学高性能计算平台上，利用ABAQUS 6.14开展计算，并结合核密度估计、Voronoi图、K函数和g函数对模型中颗粒分布的空间离散性进行多角度评价。文中指出，这四种方法能够从不同层面检验代表性体积单元（RVE，representative volume element，代表性体积单元）模型对颗粒空间分布特征的再现能力。当模型同时满足这些统计判据时，即可认为其具备较高保真度。该部分结果的意义在于证明，数值模型并非只是在几何上包含颗粒，而是在统计空间结构上也与真实样品相一致，因此其后得到的界面能量释放率分析结果具有更强可信度。

关于界面开裂驱动力的总体结果，论文以能量释放率G作为核心判据。研究人员通过VCCT计算裂纹尖端能量释放率，比较不同微观结构构型下的变化趋势。由于该参数直接表征裂纹单位扩展长度所需释放的能量，因此其数值越高，通常意味着裂纹扩展驱动力越强。文章指出，采用该参数可以有效回避双相界面裂纹尖端振荡奇异应力场带来的计算困难，从而更适合用于芯体-包壳界面断裂评估。这一结果层面虽然偏方法学，但实质上为后续讨论颗粒分布特征影响提供了统一、稳定的分析框架。

关于颗粒团聚程度的影响，论文的主要发现是：颗粒团聚会加剧芯体-包壳界面的开裂行为。研究人员通过数值模拟表明，随着近界面区域颗粒分布由相对均匀转向更明显的团聚状态，界面裂纹尖端的能量释放率升高。这说明颗粒团聚会强化局部热-机械失配，放大界面区域应力集中，从而提高裂纹扩展驱动力。该结论与颗粒增强复合材料中团聚易诱发局部失效的一般规律相一致，也进一步说明弥散燃料制备过程中颗粒分散均匀性控制的重要性。

关于团聚位置的影响，论文进一步指出，团聚区域距离弥散芯体中心越近，界面能量释放率越高。换言之，不同团聚区域虽然都可能恶化界面断裂响应，但其空间位置并非等效；位于更靠近芯体中心的团聚特征与更高的开裂驱动力相关。这一发现表明，近界面颗粒分布不仅要考察是否存在团聚，还需关注团聚相对界面的空间布局。该结果对于界面开裂位置与鼓泡位置预测具有直接意义，因为它说明微观组织的空间几何关系能够转化为界面损伤风险差异。

关于颗粒尺寸和体积分数的影响，论文在研究设计中将其列为重点考察因素，并通过高保真模型对其作用进行了分析。根据引言与全文结构说明，研究人员并非孤立讨论单颗粒作用，而是将颗粒尺寸、燃料相体积分数与非均匀分布特征耦合考虑，进而评估其对界面能量释放率G的影响。本文提供的材料中未展开具体定量结果与单项数值变化幅度，因此在解读上应保持原文边界：可以确认这两类参数属于系统分析对象，并被纳入界面裂纹驱动力评价框架，但不宜超出原文片段去推断其具体增减规律。

作为方法应用展示，研究人员还利用所提出的方法，自适应确定了类比弥散板中界面能量释放率的分布曲线，以估计可能发生界面失效的位置。该结果说明，本文方法不仅能够进行参数敏感性分析，还能够面向实际板状燃料样品开展局部失效风险识别。换言之，该研究已从“解释颗粒分布如何影响开裂”进一步拓展到“利用微观结构信息预测何处更可能开裂”。这使得该方法具备面向工程评估与制备工艺优化的应用潜力。

在讨论层面，论文的核心贡献主要体现在三个方面。其一，研究将真实颗粒分布统计特征引入界面断裂数值分析，增强了模型对实际弥散燃料微观组织的代表性。其二，研究针对双相界面裂纹分析中常见的应力场振荡难题，采用VCCT计算能量释放率G，建立了较为简洁且有效的界面开裂驱动力评估框架。其三，研究明确揭示了近界面颗粒团聚及其位置对界面开裂敏感性的影响，说明微观结构控制能够服务于鼓泡与界面失效预测。论文因此为弥散燃料制造容差确定、颗粒分布质量评估以及反应堆堆芯安全增强提供了理论与方法基础。

结论部分可译述为：该研究系统考察了弥散核燃料芯体中非均匀颗粒分布对芯体-包壳界面开裂的影响。通过采用基于实际颗粒分布建立的高保真模型，研究人员模拟了界面能量释放率，并对关键分布特征的作用进行了定性分析。所提出的方法可用于确定颗粒分布特征的容差范围，从而提升反应堆堆芯安全性。