Towards Automatic Stress Analysis using Scaled Boundary Finite Element Method with Quadtree Mesh of High-order Elements

本文提出了一种结合高阶单元四叉树网格与缩放边界有限元法（SBFEM）的自动应力与断裂分析技术，该方法无需特殊处理悬挂节点或渐近富集即可精确模拟裂纹尖端奇异性并有效处理曲线边界。

要点

这篇论文介绍了一种让计算机自动分析材料受力（比如哪里会断、哪里压力最大）的**“智能网格技术”**。

为了让你更容易理解，我们可以把这项技术想象成**“用乐高积木自动拼出复杂的地图，并精准预测哪里会裂开”**。

以下是用通俗语言和比喻对这篇论文核心内容的解读：

1. 以前的痛点：拼图太麻烦

在传统的工程分析（有限元分析 FEA）中，要把一个物体（比如一个有圆孔的金属板）画成网格（像渔网一样）来算受力，非常耗时且容易出错。

• 比喻：想象你要用正方形的小瓷砖去铺一个有圆形游泳池的地面。
  • 问题 1（挂点）：如果大瓷砖旁边是小瓷砖，它们边缘对不上，就像瓷砖中间多出了个“悬空”的角，导致地面不平，计算就会出错。
  • 问题 2（曲线）：正方形瓷砖很难完美贴合圆形的游泳池边缘，要么留缝隙，要么要把瓷砖切得粉碎，导致工作量巨大。
  • 问题 3（裂缝）：如果游泳池边缘有裂缝，传统方法需要在裂缝尖端把瓷砖切得极小极密，像显微镜一样，非常费时间。

2. 新方案：Scaled Boundary + 四叉树（SBFEM + Quadtree）

作者提出了一种新组合拳，把**“四叉树网格”（一种自动分层细分的网格技术）和“标度边界有限元法”**（SBFEM）结合起来。

核心魔法一：把“不规则”变成“多边形”

• 比喻：以前我们强迫所有瓷砖必须是正方形。现在，作者发明了一种**“万能变形胶”**。
• 原理：不管你的网格边缘是直的还是弯的，也不管旁边是大瓷砖还是小瓷砖（解决了“挂点”问题），作者直接把每个网格单元看作一个**“多边形”**。
• 效果：
  • 自动适应：如果大瓷砖旁边是小瓷砖，不需要切分，直接把它们当作一个多边形处理，就像把不同大小的积木块直接粘在一起，中间没有缝隙，也没有“悬空”的角。
  • 无需特殊处理：不需要像以前那样专门去修补那些“悬空”的节点，系统自动就能算好。

核心魔法二：用“高阶元素”画圆

• 比喻：以前用正方形瓷砖拼圆，需要切很多块（h-细化）。现在，作者允许每个瓷砖的边是**“弯曲的”**（高阶元素，p-细化）。
• 效果：就像用几块弯曲的乐高积木就能拼出一个完美的圆，而不需要把圆切成几百块小方块。这样既保留了网格的整齐，又完美贴合了复杂的曲线边界。

核心魔法三：裂缝的“透视眼”

• 比喻：这是最厉害的地方。当材料出现裂缝时，裂缝尖端是受力最集中、最容易断裂的地方（应力奇异点）。
• 原理：传统的 SBFEM 有一个特性，它把裂缝尖端当作一个“中心点”，像阳光从中心向四周辐射一样。
• 效果：
  • 无需加密：在裂缝尖端，不需要把网格切得极小。系统直接通过数学公式“看穿”了裂缝尖端的应力变化规律。
  • 自动计算：就像你不需要把放大镜放在裂缝上，系统直接告诉你裂缝尖端的应力有多大。这让分析带裂缝的结构变得超级简单。

3. 整个流程有多自动化？

作者设计了一个**“全自动流水线”**：

1. 输入：你只需要告诉电脑“这里有个圆孔”、“这里有裂缝”、“我想让这里的网格密一点（通过几个种子点控制）”。
2. 生成：电脑自动画出一个像树状图一样的网格（四叉树）。
3. 修剪：遇到圆孔或裂缝，电脑自动把边缘的网格“修剪”成多边形，并自动把裂缝尖端周围的网格合并成一个大块。
4. 计算：因为大部分网格都是标准的，电脑只需要算几种“标准模板”，剩下的直接套用，速度极快。

4. 论文证明了什么？

作者做了 5 个实验，包括：

• 带圆孔的板子（验证曲线拟合）。
• 有很多随机孔洞的板子（验证自动处理复杂形状）。
• 从圆孔长出多条裂缝的板子（验证多裂缝处理）。
• 两条裂缝交叉的板子（验证复杂裂缝）。
• 一个真实的核反应堆模型（验证实际应用）。

结果：

• 快：生成网格只需几秒钟。
• 准：计算结果和理论公式、传统商业软件（如 ANSYS）的结果高度一致。
• 省：用更少的网格数量，达到了更高的精度（特别是用了高阶元素后）。

总结

这篇论文就像给工程师配了一副**“智能眼镜”和一套“自动拼图机器人”**。
以前工程师需要花大量时间手动画网格、修补漏洞、处理裂缝尖端；现在，只需要输入简单的几何信息，系统就能自动生成完美的网格，自动处理所有复杂的边界和裂缝，并且算得又快又准。这让复杂的工程结构分析变得像搭积木一样简单和自动化。

技术摘要

论文技术总结：基于高阶单元四叉网格的自动应力分析与扩展边界有限元法

1. 研究背景与问题 (Problem)

传统的有限元分析 (FEA) 在网格生成方面面临巨大挑战，尤其是处理复杂几何边界、几何特征间的快速过渡以及奇点（如裂纹尖端）时。

• 网格质量与效率的矛盾：为了获得高精度，需要在复杂区域进行细密网格划分，但这增加了计算成本。自适应网格技术虽然能平衡精度与效率，但高质量的网格生成通常占据分析时间的绝大部分。
• 四叉网格 (Quadtree Mesh) 的局限性：四叉网格是一种高效的自适应分层网格技术，能自动生成不同尺寸的方形单元。然而，将其应用于传统有限元存在两个主要问题：
  1. 悬挂节点 (Hanging Nodes)：相邻单元尺寸不同导致边界上出现未连接的中间节点，造成位移不协调，需要特殊的形函数或单元细分处理。
  2. 曲线边界拟合：四叉网格由水平和垂直线组成，难以精确拟合曲线边界。通常需要过度细化网格或进行复杂的边界修剪，导致网格质量下降或计算量增加。
• 裂纹分析的挑战：在裂纹尖端附近，传统方法通常需要局部网格细化或使用渐近富集技术（如 XFEM）来捕捉应力奇异性。

2. 方法论 (Methodology)

本文提出了一种结合扩展边界有限元法 (Scaled Boundary Finite Element Method, SBFEM) 与高阶单元四叉网格的自动化应力与断裂分析技术。

核心原理

• SBFEM 与多边形单元：将四叉网格中的每个单元（包括被边界切割的不规则单元）视为一个扩展边界多边形 (Scaled Boundary Polygon)。
  • SBFEM 是一种半解析方法，只需在边界上离散线单元，域内的解通过缩放中心到边界的缩放关系半解析地表达。
  • 悬挂节点处理：在 SBFEM 框架下，悬挂节点被视为多边形顶点的普通节点，无需特殊的形函数推导或单元细分，直接作为一维边界线单元的端点处理，天然保证了位移协调性。
  • 高阶单元：每个单元边界上的线单元可以直接使用高阶谱单元（如 4 阶、5 阶等），无需增加节点数量即可提高精度。

算法流程

1. 几何定义：使用符号距离函数 (Signed Distance Function) 定义几何域和边界。
2. 种子点控制：用户只需在边界和感兴趣区域定义少量的“种子点”来控制网格密度，无需手动划分区域。
3. 四叉网格生成：
   • 从根单元开始递归细分，直到满足种子点密度限制和 2:1 平衡规则。
   • 边界修剪：利用符号距离函数识别与边界相交的单元，将其修剪为多边形。对于靠近边界的顶点，若距离过近则投影到边界上，避免产生极短边。
   • 曲线拟合：利用高阶线单元（p-细化）精确拟合曲线边界，避免了传统 h-细化带来的过度网格加密。
4. 裂纹处理：
   • 裂纹尖端：将裂纹尖端选为缩放中心。包含裂纹尖端的多个四叉单元被合并为一个大的多边形单元。
   • 裂纹路径：裂纹路径上的单元被裂纹线分割。
   • 奇异性：SBFEM 天然能解析裂纹尖端的应力奇异性（通过特征值分析），无需局部网格细化或渐近富集。
5. 刚度矩阵组装：利用四叉网格的层级结构，仅计算少数几种“主单元”（Master Cells）的刚度矩阵，其余单元通过旋转或提取获得，极大提高了计算效率。

3. 主要贡献 (Key Contributions)

1. 悬挂节点的无缝处理：提出了一种无需单元细分或特殊形函数的方法，将悬挂节点直接作为多边形顶点处理，彻底解决了四叉网格在传统 FEA 中的位移不协调问题。
2. 完全自动化的网格生成：仅需定义几何边界和少量种子点，即可自动生成包含高阶单元、适应复杂边界和裂纹的平衡四叉网格。
3. 高精度曲线边界拟合：通过结合 SBFEM 和高阶单元，能够用较少的单元精确拟合曲线边界，避免了传统方法中因边界拟合不良导致的应力集中误差。
4. 高效的裂纹分析：利用 SBFEM 的半解析特性，直接在裂纹尖端单元内解析应力奇异性，无需局部网格细化，显著降低了裂纹扩展模拟的网格重构成本。
5. 计算效率优化：基于四叉树结构的刚度矩阵组装策略，大幅减少了重复计算，使得大规模高阶单元分析成为可能。

4. 数值结果 (Results)

论文通过五个算例验证了该方法的有效性和精度：

1. 无限板圆孔拉伸：
   • 与解析解对比，位移和应力结果高度吻合。
   • 展示了高阶单元（p-细化）相比低阶单元具有更快的收敛速度。
   • 通过改变板尺寸比 ($L/a$)，验证了该方法模拟无限域的能力，且网格模式保持一致。
2. 多孔正方形板：
   • 处理了 9 个随机分布的圆孔，展示了在几何特征过渡区的自动化网格生成能力。
   • 与商业软件 ANSYS 对比，在相同或更少自由度下，高阶 SBFEM 获得了更高的精度和更快的收敛率。
3. 孔边多裂纹：
   • 模拟了从中心孔发出的多条裂纹（2、4、8 条）。
   • 无需裂纹尖端局部细化，计算得到的应力强度因子 (SIF) 与参考文献高度一致。
4. 交叉裂纹：
   • 处理了两条相互交叉的裂纹，展示了算法处理复杂裂纹构型的能力。
   • 网格生成完全自动，无需用户定义裂纹尖端区域或进行试错性细化。
5. 带裂纹的核反应堆压力容器：
   • 模拟了具有复杂几何形状的实际工程结构。
   • 仅需极短的时间（<0.8s）生成网格，且高阶单元在少量单元下即可达到高精度。

5. 意义与结论 (Significance & Conclusion)

• 工程应用价值：该技术为复杂几何结构（特别是包含裂纹和曲线边界的结构）的应力和断裂分析提供了一种全自动、高效且高精度的解决方案。
• 消除人工干预：极大地减少了网格生成所需的人工操作和试错过程，特别适用于需要频繁进行网格重构的裂纹扩展模拟。
• 理论突破：成功将 SBFEM 的半解析优势与四叉网格的自动化优势相结合，克服了传统有限元在处理悬挂节点和曲线边界时的固有缺陷。
• 未来展望：该方法为大规模工程问题的自动化仿真设计奠定了基础，特别是在需要高精度应力分析和断裂力学评估的领域具有广阔的应用前景。

总结：本文提出了一种基于 SBFEM 和高阶四叉网格的自动化分析框架，通过巧妙处理悬挂节点和曲线边界，实现了无需局部细化即可精确模拟裂纹奇异性，显著提升了计算效率和网格生成的自动化水平。