基于μCT图像的球谐随机场标定以生成喷射土三维骨料

《Probabilistic Engineering Mechanics》：Calibration of spherical harmonic-based random fields from μCT images to generate shot-earth 3-D aggregates

　　

喷射土（shot-earth）是一种环保材料，近年来在建筑领域的许多应用中重新受到关注。然而，由于其复杂的内部结构——其特征是存在形状高度多变的骨料——其力学表征面临诸多挑战。研究人员旨在开发并验证一个用于喷射土骨料几何概率建模的计算框架。研究人员利用基于球谐函数（spherical harmonic）的随机场（random field）结合蒙特卡洛（Monte Carlo）模拟生成真实的三维（3-D）骨料。尽管这些随机和数学工具先前已应用于类混凝土材料的细观尺度建模，但其在喷射土中的应用尚未被探索。生成过程基于从真实立方体喷射土样本的微计算机断层扫描（μCT）图像中获得的尺寸数据。对这些数据的统计分析用于设定构建基于球谐函数的随机场的假设，并推断所需系数的概率分布。将生成的虚拟骨料与原始骨料在体积、转动惯量和整体形状方面进行比较。定量比较表明与实验数据具有良好一致性，证明了所提出框架的有效性。该方法能够实现更准确的喷射土骨料细观尺度模拟，并为研究其力学行为的有限元建模的未来发展提供支持，最终促进喷射土在可持续建筑中的更广泛应用。此外，详细的工作流程提供了一种可重复的方法论，可推广到其他颗粒材料，从μCT数据到虚拟样本的生成。

论文解读：基于μCT图像和球谐随机场生成喷射土三维骨料

土基建筑材料历史悠久，但曾被更高性能的材料取代。近年来，随着环保意识增强，建筑行业重新关注低碳、天然的土基材料。其中，喷射土（shot-earth）作为一种新型夯土材料，通过高速喷射技术将稳定土与骨料结合，在可持续建筑领域显示出潜力。然而，其内部结构复杂，骨料形状高度不规则，导致力学特性难以准确表征。在细观尺度模拟中，骨料形态对颗粒间的摩擦、互锁及力链传递起决定性作用，因此精确构建骨料几何模型是关键。传统方法常将骨料简化为球体、椭球或多面体，虽计算高效但无法捕捉真实骨料的随机性和表面粗糙度。近年来，基于数字图像处理和三维（3-D）重建技术，球谐函数（spherical harmonic, SH）分析已成功用于生成多种颗粒材料的虚拟骨料，并能与随机场（random field）结合以引入形状随机性。然而，现有研究普遍缺乏对模拟精度的深入分析与实验校准，尤其对于研究较少的喷射土骨料，其形态特性尚未得到系统建模。为此，研究人员在《Probabilistic Engineering Mechanics》期刊上发表了这项研究，旨在建立并验证一个从微计算机断层扫描（μCT）数据到虚拟骨料生成的概率框架，从而为喷射土细观尺度模拟提供更真实的几何输入。

研究人员采用的主要技术方法包括：利用现有文献提供的立方体喷射土试样的μCT扫描图像，通过图像预处理提取完整骨料；采用谱聚类（spectral clustering）基于无量纲惯性矩比k₁、k₂将约1600个骨料分为5组，以降低计算规模；对每组骨料应用截断阶数n=20的球谐展开，获得441个系数并构建系数矩阵C；统计系数后发现其分布偏离高斯分布，且尺寸系数a_0,0与形状系数相关性弱，因此按尺寸归一化以分离尺寸与形状；构建乘积型随机场G=V·G?，其中V描述尺寸随机变量，G?描述单位尺寸缩放几何，通过平移模型将标准高斯随机场映射为具有Beta边缘分布的非高斯场；最后采用蒙特卡洛（Monte Carlo）模拟生成虚拟骨料。整个流程基于真实实验数据校准，可复现并推广至其他颗粒材料。

研究结果按论文主体结构总结如下。

1. 引言（Introduction）——论文指出现有球谐分析结合随机场的研究缺乏对模拟精度的系统校准和验证，尤其是喷射土骨料形态尚未被深入探究，因此需要开展专门研究。

2. 球谐分析用于三维形状近似（Spherical harmonic analysis for 3-D shape approximation）——通过球坐标变换将骨料表面表示为径向距离函数g(θ,φ)，并利用实系数球谐级数展开，其中零阶项代表等效球半径（即骨料尺寸），高阶项代表形状扰动。该理论为后续骨料重建奠定了数学基础。

3. 喷射土实验数据集（Shot-earth experimental dataset）——包括4个子节：3.1 预处理：从μCT图像中识别骨料，剔除被试样边界截断的颗粒以及体积小于1 mm³的颗粒，最终获得约1600个完整骨料。3.2 聚类：基于骨料主惯性矩比k₁、k₂作为形状描述符，使用谱聚类将骨料分为5组，从而在保留形态差异的同时提高计算效率。3.3 表面点检测与球坐标参数化：提取每个骨料的边界表面点，通过插值得到连续函数g(θ,φ)，为球谐分析提供输入。3.4 球谐分析应用：对每组骨料在截断阶n=20下计算441个球谐系数，统计结果显示系数具有明显的偏度和峰度，表明其非高斯特性；同时发现尺寸系数与形状系数弱相关，因此将所有系数按尺寸归一化，得到单位尺寸的缩放骨料。

4. 基于球谐的随机场模型（Spherical harmonics-based random field model）——基于上述弱相关性，研究人员建立随机场G(θ,φ)=V·G?(θ,φ)，其中V为骨料尺寸随机变量，G?为描述单位尺寸骨料几何的非高斯随机场。G?通过平移模型由标准高斯随机场N_n(θ,φ)变换得到，其边缘分布采用三参数Beta分布拟合；V也采用三参数Beta分布，并利用实际骨料尺寸数据校准。

5. 蒙特卡洛模拟（Monte Carlo simulations）——按照所提出的五步算法生成虚拟骨料，并与真实骨料进行比较。结果表明：虚拟骨料的惯性比k₁、k₂与真实骨料分布一致；对第3组和第4组，缩放体积均值误差为11%–17%，绝对体积均值误差分别为4.30%和3.23%，一致性良好。

6. 结论（Conclusions）——研究人员提出的概率建模框架成功地从μCT图像生成了与实验数据高度吻合的喷射土虚拟骨料，验证了球谐分析与随机场结合的有效性，并展示了模型在分辨率与计算成本之间的灵活性。该方法不仅适用于喷射土，还可推广至其他颗粒材料，为后续细观尺度有限元模拟和材料性能预测提供了可靠工具。

综上所述，研究结论表明：本工作提出的基于球谐随机场的概率模型，通过系统校准和验证，能够准确再现喷射土骨料的几何统计特征，显著提升了细观尺度模拟的真实性，为可持续建筑材料的研究与应用提供了重要支撑。