[论文解读] Towards Automatic Stress Analysis using Scaled Boundary Finite Element Method with Quadtree Mesh of High-order Elements
本文提出了一种基于四叉树高阶单元的缩放边界有限元法(SBFEM)的自动应力与断裂分析框架。通过将四叉树单元建模为高阶多边形单元,该方法消除了悬挂节点问题,并在无需增强的情况下精确捕捉应力奇异性,从而实现完全自动化的网格生成,且用户输入极少,同时在裂纹和曲边界建模中表现出高收敛性与高精度。
This paper presents a technique for stress and fracture analysis by using the scaled boundary finite element method (SBFEM) with quadtree mesh of high-order elements. The cells of the quadtree mesh are modelled as scaled boundary polygons that can have any number of edges, be of any high orders and represent the stress singularity around a crack tip accurately without asymptotic enrichment or other special techniques. Owing to these features, a simple and automatic meshing algorithm is devised. No special treatment is required for the hanging nodes and no displacement incompatibility occurs. Curved boundaries and cracks are modelled without excessive local refinement. Five numerical examples are presented to demonstrate the simplicity and applicability of the proposed technique.
研究动机与目标
- 开发一种用于结构力学中应力与断裂问题的完全自动化的网格生成与分析技术。
- 解决基于四叉树的有限元分析中悬挂节点与曲边界拟合的挑战。
- 在模拟裂纹尖端时,无需使用特殊形函数、子三角剖分或渐近增强。
- 通过极少的用户定义参数,实现在应力与断裂分析中的高阶精度与快速收敛。
- 在具有裂纹与曲边界的复杂几何体上,验证该方法的鲁棒性与高效性。
提出的方法
- 该方法采用缩放边界有限元法(SBFEM),将四叉树单元建模为任意边数的高阶多边形单元。
- 通过将悬挂节点视为多边形单元的顶点,自然处理了该问题,避免了位移不相容性。
- 曲边界与裂纹尖端通过缩放边界多边形表示,无需过度细化网格。
- 直接在每个四叉树单元上应用高阶形函数,实现精确的应力表示与高收敛速率。
- 网格生成完全自动化,仅需通过种子点控制细化密度,裂纹尖端无需特殊处理。
- 通过利用四叉树的层次结构,仅计算必要单元,高效地组装刚度矩阵。
实验结果
研究问题
- RQ1四叉树网格是否可在有限元分析中有效使用,而不会引入悬挂节点问题?
- RQ2SBFEM中的高阶多边形单元是否能在无需渐近增强的情况下,准确模拟裂纹尖端的应力奇异性?
- RQ3是否可在用户输入极少的前提下实现自动网格生成,同时保持对曲边界与奇异性模拟的高精度?
- RQ4所提出方法的收敛速率与相同网格密度下的标准FEM相比如何?
- RQ5该方法是否能以极低计算成本高效处理复杂工程问题,如带裂纹的核反应堆?
主要发现
- 对于双裂纹板模型,2阶与5阶单元之间的差异小于0.04%,表明收敛速度极快。
- 在带裂纹的核反应堆示例中,2阶与5阶单元的裂纹张开位移差异小于0.1%,证实了高阶单元的高精度。
- 对于包含160个单元的反应堆象限四叉树网格,网格生成时间低于0.8秒,使用5阶单元的求解时间低于0.7秒。
- 该方法无需对悬挂节点进行特殊处理,因为多边形单元通过基于顶点的公式自然容纳了它们。
- 高阶单元的使用使得曲边界可精确表示,而无需过度局部细化,从而减少了所需单元数量。
- 刚度矩阵计算高效,仅需计算160个单元中的44个,展示了良好的计算可扩展性。
更好的研究,从现在开始
从阅读论文到最终审阅,大幅缩短您的研究时间。
无需绑定信用卡
本解读由 AI 生成,并经人工编辑审核。