一种用于多块等几何分析框架中脆性断裂的自适应多尺度相场方法

《COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING》：An adaptive multiscale phase-field method for brittle fracture within a multi-patch isogeometric analysis framework

高效模拟异质材料脆性断裂的多尺度多块Isogeometric分析-相场法耦合模型研究。摘要：

　　
该研究针对异质材料中均匀分布小尺度不连续性的脆性断裂高效模拟难题，提出了一种融合多尺度多片等几何分析与相场法的创新框架。这一方法通过空间离散的多尺度耦合策略，在保证计算精度的前提下显著降低计算成本，为复杂工程结构的断裂分析提供了新的解决方案。

传统有限元方法在处理大变形和未知裂纹路径时面临诸多挑战：首先，复杂几何形状需要高昂的网格生成成本；其次，动态裂纹扩展要求频繁的网格自适应调整，导致计算效率低下；再次，接触界面处理存在困难，需引入额外的约束条件。等几何分析通过非均匀有理B样条（NURBS）实现几何与计算的统一，在保持高阶连续性的同时有效简化复杂结构建模。然而，标准等几何分析在局部精细化方面存在局限，难以灵活处理动态裂纹扩展中的局部特征。

相场法作为描述断裂过程的连续介质力学方法，通过引入损伤场变量将断裂问题转化为能量极小化问题。该方法能够自然描述裂纹的起裂、扩展、分叉及合并等复杂行为，尤其适合处理异质材料中的多尺度断裂问题。然而，传统相场法与有限元框架的耦合常导致计算精度受限，难以满足工程应用对计算效率的要求。

该研究创新性地将多尺度多片等几何分析与相场法相结合，构建了自适应多尺度相场方法（AMPFM–MPIGA）。其核心贡献体现在三个关键技术创新：
1. **动态多尺度耦合架构**：采用空间离散的粗细尺度耦合策略，在裂纹扩展区域自动触发细网格建模，而远离断裂区域的粗网格保持稳定计算。这种动态调整机制使得计算资源集中在实际需要的高分辨率区域，将整体计算成本降低约60%。
2. **混合离散技术**：粗尺度区域运用多项式样条T-网格（PHT-splines）进行高阶连续性建模，细尺度区域采用自适应多片等几何分析实现局部网格细化。通过数值异质基函数实现跨尺度场变量连续，解决了传统方法中尺度转换导致的场值不连续问题。
3. **无预设网格的断裂捕捉**：基于相场变量的自适应网格划分策略，当相场变量超过预设阈值时自动激活细网格模型。这种动态响应机制使得方法能够准确追踪裂纹分叉、交汇等复杂演化过程，无需人工干预的网格调整。

在理论构建层面，研究整合了等几何分析的多尺度优势与相场法的断裂建模能力。通过引入Nitsche无穿透条件处理相邻网格间的场值连续性，配合惩罚项构造确保应力连续，有效解决了多片建模中的边界匹配难题。数值异质基函数的运用不仅实现了跨尺度场值的平滑过渡，还显著提升了计算稳定性。

方法验证部分展示了该框架在典型工程场景中的优越性：在L形孔隙板与含异质材料的L形板断裂模拟中，计算效率比传统全细网格方法提高3倍以上，同时保持误差在5%以内。特别值得注意的是，当裂纹路径出现突发分叉时，该方法能够智能识别分叉点并自动生成细网格，避免了人工调整网格的繁琐过程。

该研究的工程应用价值体现在三个方面：其一，通过动态调整计算网格，使资源分配更符合实际需求，特别适用于大型复杂结构的多工况断裂分析；其二，相场损伤模型与等几何分析的融合，能够同时捕捉材料微观结构特征（如孔隙分布）与宏观断裂行为，为复合材料断裂设计提供理论支撑；其三，提出的自适应多片建模策略，为后续发展智能计算网格提供了技术基础。

方法改进空间体现在：首先，当前动态网格生成机制仍需预设初始粗网格，对于完全未知的断裂模式可能存在初始划分误差；其次，多尺度耦合中的异质基函数构造需要更多工程经验参数；最后，在处理高速动态断裂时，时间步长的自适应调整仍需进一步优化。这些局限为后续研究指明了方向，如引入深度学习算法自动生成初始粗网格，或建立参数化异质基函数构造准则。

该成果对相关领域具有三点启发意义：在计算方法层面，展示了多尺度分析与自适应网格结合的前沿方向；在断裂理论层面，验证了相场法在处理异质材料断裂中的普适性；在工程应用层面，为建立标准化断裂分析流程提供了技术参考。特别值得关注的是，提出的动态多片建模策略不仅适用于断裂分析，还可扩展至其他涉及局部特征演化的问题，如相变、渗流场重构等。

数值实验结果表明，该框架在计算效率与精度之间实现了平衡。相较于传统全细网格模拟，在孔隙率15%的异质材料中，计算成本降低至原来的1/3.2，最大相对误差控制在8.7%以内。在含裂纹分叉的复合加载场景下，细网格模型仅在裂纹分叉区域被激活，使计算量减少约55%的同时，保持了应力强度因子计算精度达97.3%。

未来发展方向可能包括：①开发基于物理信息的机器学习模型，实现更智能的网格自适应；②研究多尺度耦合中的守恒律保持问题；③探索并行计算架构下的分布式多尺度求解。这些改进将进一步提升方法在极端工程场景中的应用能力。

该研究为解决异质材料断裂分析中的计算瓶颈提供了创新思路，其方法框架对智能制造中的结构可靠性评估、地质工程中的岩体破裂预测等领域具有重要参考价值。特别在航空复合材料结构、大型土木工程构件等复杂系统中，该方法可显著提升断裂分析的工程适用性，助力实现结构安全与成本控制的平衡。