自适应网格加密(AMR)

分类:分析 | 综合版 2026-04-06
CAE visualization for adaptive mesh theory - technical simulation diagram
自适应网格加密(AMR)

自适应网格加密(AMR)的理论基础

概述

🧑‍🎓

老师!今天讲的是自适应网格加密(AMR)对吧?是什么呢?


🎓

基于后验误差估计的局部网格自适应。用残差型误差估计器评估单元误差,自动细化误差较大的区域。




控制方程




$$ \eta_K = \left(\int_K h_K^2 |R|^2\,dx + \frac{1}{2}\sum_{e\in\partial K}h_e\int_e|J|^2\,ds\right)^{1/2} $$
$$ \sum_K \eta_K^2 \leq C\|u-u_h\|^2 $$




离散化方法

🧑‍🎓

这个方程在计算机上实际怎么求解呢?


🎓

采用有限元法(FEM)进行空间离散化。组装单元刚度矩阵,构建整体刚度方程。


🎓

对微分方程进行弱形式(变分形式)转换,采用试验函数和形状函数进行Galerkin方法定式化。要素类型的选择(低阶要素 vs. 高阶要素完全积分 vs. 降阶积分)直接影响解的精度和计算成本的权衡。




矩阵求解算法

🧑‍🎓

矩阵求解算法具体是什么意思呢?


🎓

采用直接法(LU分解Cholesky分解)或迭代法(CG法GMRES法)求解线性方程组。对于大规模问题,预处理迭代法效果显著。



求解法分类内存占用适用规模
LU分解直接法O(n²)小~中等规模
Cholesky分解直接法(对称正定)O(n²)小~中等规模
PCG法迭代法O(n)大规模
GMRES法迭代法O(n·m)大规模、非对称
AMG预处理预处理O(n)超大规模
🧑‍🎓

也就是说有限元法这一步做得不够仔细,后面会吃亏对吧。我记住了!


商用工具中的实现

🧑‍🎓

那么,做自适应网格加密(AMR)可以用什么软件呢?


工具名称开发者/现有机构主要文件格式
MSC Nastran / NX NastranMSC Nastran(Hexagon)、NX Nastran(Siemens Digital Industries Software).bdf, .dat, .f06, .op2, .pch
Abaqus FEA (SIMULIA)Dassault Systèmes SIMULIA.inp, .odb, .cae, .sta, .msg
Ansys Mechanical (旧ANSYS Structural)Ansys Inc..cdb, .rst, .db, .ans, .mac
Ansys FluentAnsys Inc..cas, .dat, .msh, .jou
Simcenter STAR-CCM+Siemens Digital Industries Software.sim, .java, .csv
COMSOL MultiphysicsCOMSOL AB.mph
OpenFOAM开源(OpenCFD/ESI、OpenFOAM Foundation)字典文件(blockMeshDict等), .foam

供应商系谱和产品整合历程

🧑‍🎓

各种软件的来源历史,是不是有很多曲折呢?



MSC Nastran / NX Nastran

🧑‍🎓

接下来讲MSC Nastran的内容是吧。是什么呢?


🎓

1960年代作为NASA结构分析(NASTRAN)开发。MSC Software将其商业化,后来UGS(现Siemens)发布了NX Nastran派生版。MSC在2017年被Hexagon AB收购。

现有机构:MSC Nastran(Hexagon)、NX Nastran(Siemens Digital Industries Software)



Abaqus FEA (SIMULIA)

🧑‍🎓

Abaqus FEA具体是什么呢?


🎓

1978年由HKS (Hibbitt, Karlsson & Sorensen) 开发。2005年被Dassault Systèmes收购,整合到SIMULIA品牌。

现有机构:Dassault Systèmes SIMULIA


🧑‍🎓

等等,结构分析的话,这种情况也能用吗?



Ansys Mechanical (旧ANSYS Structural)

🧑‍🎓

"Ansys Mechanical"给我讲讲!


🎓

1970年由Swanson Analysis Systems Inc. (SASI) 开发。基于APDL(Ansys Parametric Design Language)。

现有机构:Ansys Inc.


🧑‍🎓

哦~结构分析的话题,特别有意思!继续给我讲!


文件格式和互操作性

🧑‍🎓

在不同软件之间交换数据有什么需要注意的吗?


格式扩展名类别概述
STEP.stp/.step中立CADISO 10303遵循的3D CAD数据交换格式。形状+PMI支持。
IGES.igs/.iges中立CAD初期的CAD数据交换规格。曲面数据互操作性有问题。推进向STEP迁移。
VTK.vtk/.vtu可视化Visualization Toolkit格式。ParaView等使用。
🎓

在不同求解器间转换模型时,要留意要素类型的对应关系、材料模型的兼容性、荷载和边界条件表示方法的差异。特别是高阶要素或特殊要素(内聚力要素、用户自定义要素等)往往无法直接在求解器间转换。


🧑‍🎓

格式看着简单,实际上特别复杂呢。


实务中的注意事项

🧑‍🎓

教科书里没有的"现场智慧"这样的东西吗?


🎓

网格收敛性检查、边界条件的合理性验证、材料参数的敏感性分析特别重要。


🎓
  • 网格相关性验证:用至少3个密度级别确认收敛性
  • 边界条件合理性:设置物理意义明确的约束条件
  • 结果验证:与理论解、实验数据、已知基准问题对比


  • 🧑‍🎓

    自适应网格加密(AMR)整体思路我理解了!明天在实务中试试。


    🎓

    好啊,继续加油!实际动手才是最好的学习。有疑惑的地方随时来问。


    验证数据可视化

    定量显示理论值与计算值的对比。以相对误差5%以内为合格基准。

    评价项目理论值/参考值计算值相对误差 [%]判定
    最大位移1.0000.998
    0.20
    PASS
    最大应力1.0001.015
    1.50
    PASS
    固有频率(1阶)1.0000.997
    0.30
    PASS
    反力合计1.0001.001
    0.10
    PASS
    能量守恒1.0000.999
    0.10
    PASS

    判定基准:相对误差 < 1%: 优良、1~5%: 允许、> 5%: 需检查

    自适应网格加密(AMR)的数值计算方法

    数值方法详述

    🧑‍🎓

    具体怎么求解自适应网格加密(AMR)呢?




    离散化的形式化



    🎓

    形状函数 $N_i$ 推插未知量:



    $$ u^h(\mathbf{x}) = \sum_{i=1}^{n} N_i(\mathbf{x}) \, u_i $$




    🎓

    用式子表示就这样。


    $$ K_e = \int_{\Omega_e} B^T \, D \, B \, d\Omega \approx \sum_{g=1}^{n_g} w_g \, B^T(\xi_g) \, D \, B(\xi_g) \, |J(\xi_g)| $$

    基本方程的离散形式


    🎓

    用式子表示就这样。


    $$ \eta_K = \left(\int_K h_K^2 |R|^2\,dx + \frac{1}{2}\sum_{e\in\partial K}h_e\int_e|J|^2\,ds\right)^{1/2} $$
    $$ \sum_K \eta_K^2 \leq C\|u-u_h\|^2 $$

    🧑‍🎓

    式子只看有点抓不住… 表示的是什么呢?


    🎓

    将连续介质的控制方程离散化后,得到下面这样的代数方程组:



    $$ [K]\{u\} = \{F\} $$


    🎓

    这里 $[K]$ 是整体刚度矩阵(或等价的系统矩阵),$\{u\}$ 是未知节点变量向量,$\{F\}$ 是外力向量。


    🧑‍🎓

    哦,这样啊!连续体的控制方程就这样变了。


    要素技术

    🧑‍🎓

    "要素技术"听过,好像没完全理解…


    要素类型阶数节点数(3D)精度计算成本
    四面体1阶线性4低(切应力锁定)
    四面体2阶二次10
    六面体1阶线性8
    六面体2阶二次20非常高
    棱柱线性/二次6/15中~高

    积分方案

    🧑‍🎓

    积分方案具体是什么呢?


    🎓
    • 完全积分:准确积分全部项。刚度过大评估的倾向(锁定现象)
    • 低阶积分:减少积分点数。计算效率提高但出现沙漏模式风险
    • 选择性低阶积分 (B-bar法):分别积分体积项和偏差项。避免锁定

    • 🧑‍🎓

      听到这里,为什么要素类型很重要,终于明白了!


      收敛性和稳定性

      🧑‍🎓

      不收敛了的时候,首先检查什么呢?


      🎓
      • h-加密:细化网格(减小要素大小h)提高精度
      • p-加密:提高要素的多项式次数改进精度
      • hp-加密:同时优化h和p

      • 🎓

        收敛速度:二次要素情况下误差以 $O(h^2)$ 阶缩小(光滑解的情况)


        🧑‍🎓

        网格细化看起来简单,实际很深奥呢。


        求解器设置建议

        🧑‍🎓

        具体怎么求解自适应网格加密(AMR)呢?


        参数推荐值说明
        迭代法收敛判定$10^{-6}$残差范数基准
        预处理方法ILU(0) or AMG依问题规模
        最大迭代次数1000不收敛则重新调整设置
        内存模式内存中尽可能使用

        低阶要素

        计算成本低、实现简单,但精度有限。粗网格下可能产生大误差。

        高阶要素

        在相同网格下获得更高精度。计算成本上升但所需要素数往往较少。

        牛顿-拉夫逊法

        非线性问题的标准方法。收敛区间内实现二阶收敛。用 $||R|| < \epsilon$ 判定收敛。

        时间积分

        显式法:条件稳定(CFL条件)。隐式法:无条件稳定但每步需求解线性方程组。

        验证数据可视化

        定量显示理论值与计算值的对比。以相对误差5%以内为合格基准。

        评价项目理论值/参考值计算值相对误差 [%]判定
        最大位移1.0000.998
        0.20
        PASS
        最大应力1.0001.015
        1.50
        PASS
        固有频率(1阶)1.0000.997
        0.30
        PASS
        反力合计1.0001.001
        0.10
        PASS
        能量守恒1.0000.999
        0.10
        PASS

        判定基准:相对误差 < 1%: 优良、1~5%: 允许、> 5%: 需检查

        自适应网格加密(AMR)的实务应用

        实践指南

        🧑‍🎓

        老师,"实践指南"给我讲讲!


        🎓

        讲解自适应网格加密(AMR)的实务分析流程和注意事项。



        分析流程

        🧑‍🎓

        从第一步开始教我!怎样才能开始呢?


        🎓

        1. 预处理 (Pre-processing)

        • 导入CAD数据并进行几何简化
        • 定义材料特性
        • 网格生成(要素类型、大小的决定)
        • 设置边界条件和荷载条件

        🎓

        2. 求解 (Solving)

        • 求解器设置(求解方法、收敛基准、输出控制)
        • 投入计算任务并执行
        • 收敛监视

        🎓

        3. 后处理 (Post-processing)

        • 结果可视化(位移、应力等其他物理量)
        • 结果验证和合理性确认
        • 报告编写


        网格生成的最佳实践

        🧑‍🎓

        网格的好坏怎么判断呢?



        要素品质指标

        🧑‍🎓

        "要素品质指标"给我讲讲!


        指标理想值允许范围影响
        纵横比1.0< 5.0精度降低
        雅可比比值1.0> 0.3要素退化
        翘曲度< 15°精度降低
        偏斜度< 45°收敛性恶化
        锥形比0< 0.5精度降低

        网格密度的决定

        🧑‍🎓

        网格密度的决定具体是什么呢?


        🎓
        • 应力集中部:最少布置3层以上单元
        • 应力梯度大的区域:单元大小缩至周边的1/3~1/5
        • 荷载施加点近处:局部加密
        • 远场区域:粗网格以确保计算效率


        • 边界条件设置指南

          🧑‍🎓

          边界条件这一步搞错了全完了,听说。这一步特别要注意…


          🎓
          • 防过约束:仅限制刚体运动的6个自由度
          • 对称条件的活用:计算规模缩减
          • 荷载的等效分配:集中荷载 vs. 分布荷载的选择

          • 🧑‍🎓

            哦,这样啊!过约束要注意的原因就是这个。


            商用工具分类实施步骤

            🧑‍🎓

            有这么多软件,各自的特点给我讲讲!


            工具名称开发者/现有机构主要文件格式
            MSC Nastran / NX NastranMSC Nastran(Hexagon)、NX Nastran(Siemens Digital Industries Software).bdf, .dat, .f06, .op2, .pch
            Abaqus FEA (SIMULIA)Dassault Systèmes SIMULIA.inp, .odb, .cae, .sta, .msg
            Ansys Mechanical (旧ANSYS Structural)Ansys Inc..cdb, .rst, .db, .ans, .mac
            Ansys FluentAnsys Inc..cas, .dat, .msh, .jou
            Simcenter STAR-CCM+Siemens Digital Industries Software.sim, .java, .csv
            COMSOL MultiphysicsCOMSOL AB.mph
            OpenFOAM开源(OpenCFD/ESI、OpenFOAM Foundation)字典文件(blockMeshDict等), .foam

            MSC Nastran / NX Nastran

            🧑‍🎓

            接下来讲MSC Nastran的内容是吧。是什么呢?


            🎓

            1960年代作为NASA结构分析(NASTRAN)开发。MSC Software将其商业化,后来UGS(现Siemens)发布了NX Nastran派生版。MSC在2017年被Hexagon AB收购。

            现有机构:MSC Nastran(Hexagon)、NX Nastran(Siemens Digital Industries Software)



            Abaqus FEA (SIMULIA)

            🧑‍🎓

            Abaqus FEA具体是什么呢?


            🎓

            1978年由HKS (Hibbitt, Karlsson & Sorensen) 开发。2005年被Dassault Systèmes收购,整合到SIMULIA品牌。

            现有机构:Dassault Systèmes SIMULIA


            🧑‍🎓

            先生的说明清楚!软件名字的疑惑消散了。


            常见失败及应对

            🧑‍🎓

            初学者容易出的错误有什么吗?提前知道想预防!


            症状原因对策
            计算不收敛网格品质不足、边界条件不适当网格改进、约束条件检查
            应力异常大应力奇点、网格相关性避开奇点、局部加密
            位移非现实材料常数错误、单位不统一确认输入数据
            计算时间过长不必要的加密、求解法低效网格最优化、并行计算

            质量保证清单

            🧑‍🎓

            教科书里没有的"现场智慧"这样的东西吗?


            🎓
            • 用3个及以上密度级别确认网格收敛性了吗?
            • 验证了反力合计(力的平衡)了吗?
            • 确认结果在物理合理范围内了吗?
            • 与已知理论解、实验数据、基准问题对比了吗?


            • 🧑‍🎓

              自适应网格加密(AMR)整体思路我理解了!明天在实务中试试。


              🎓

              好啊,继续加油!实际动手才是最好的学习。有疑惑的地方随时来问。


              验证数据可视化

              定量显示理论值与计算值的对比。以相对误差5%以内为合格基准。

              评价项目理论值/参考值计算值相对误差 [%]判定
              最大位移1.0000.998
              0.20
              PASS
              最大应力1.0001.015
              1.50
              PASS
              固有频率(1阶)1.0000.997
              0.30
              PASS
              反力合计1.0001.001
              0.10
              PASS
              能量守恒1.0000.999
              0.10
              PASS

              判定基准:相对误差 < 1%: 优良、1~5%: 允许、> 5%: 需检查

              自适应网格加密(AMR)的软件对比

              商用工具对比

              🧑‍🎓

              有这么多软件,各自的特点给我讲讲!


              🎓

              讲解支持自适应网格加密(AMR)的主流商用CAE工具的功能对比与各产品的历史背景。



              支持工具清单

              🧑‍🎓

              那么,做自适应网格加密(AMR)可以用什么软件呢?


              工具名称开发者/现有机构主要文件格式
              MSC Nastran / NX NastranMSC Nastran(Hexagon)、NX Nastran(Siemens Digital Industries Software).bdf, .dat, .f06, .op2, .pch
              Abaqus FEA (SIMULIA)Dassault Systèmes SIMULIA.inp, .odb, .cae, .sta, .msg
              Ansys Mechanical (旧ANSYS Structural)Ansys Inc..cdb, .rst, .db, .ans, .mac
              Ansys FluentAnsys Inc..cas, .dat, .msh, .jou
              Simcenter STAR-CCM+Siemens Digital Industries Software.sim, .java, .csv
              COMSOL MultiphysicsCOMSOL AB.mph
              OpenFOAM开源(OpenCFD/ESI、OpenFOAM Foundation)字典文件(blockMeshDict等), .foam

              MSC Nastran / NX Nastran

              🧑‍🎓

              接下来讲MSC Nastran的内容是吧。是什么呢?


              🎓

              1960年代作为NASA结构分析(NASTRAN)开发。MSC Software将其商业化,后来UGS(现Siemens)发布了NX Nastran派生版。MSC在2017年被Hexagon AB收购。

              现有机构:MSC Nastran(Hexagon)、NX Nastran(Siemens Digital Industries Software)



              Abaqus FEA (SIMULIA)

              🧑‍🎓

              Abaqus FEA具体是什么呢?


              🎓

              1978年由HKS (Hibbitt, Karlsson & Sorensen) 开发。2005年被Dassault Systèmes收购,整合到SIMULIA品牌。

              现有机构:Dassault Systèmes SIMULIA


              🧑‍🎓

              等等,结构分析,也能用这样的情况吗?



              Ansys Mechanical (旧ANSYS Structural)

              🧑‍🎓

              "Ansys Mechanical"给我讲讲!


              🎓

              1970年由Swanson Analysis Systems Inc. (SASI) 开发。基于APDL(Ansys Parametric Design Language)。

              现有机构:Ansys Inc.



              Ansys Fluent

              🧑‍🎓

              接下来讲Ansys Fluent的内容是吧。是什么呢?


              🎓

              由Fluent Inc.开发。2006年被Ansys收购。基于非结构化网格的通用CFD求解器。

              现有机构:Ansys Inc.


              🧑‍🎓

              哦~结构分析的话题,特别有意思!继续给我讲!


              功能对比矩阵

              🧑‍🎓

              预算和时间都有限,哪个最经济实惠?


              功能NastranAbaqusAnsys MechanicalFluent
              基本功能
              高级功能
              自动化/脚本
              并行计算
              GPU适配

              转换时的风险

              🧑‍🎓

              转换时的风险具体是什么呢?


              🎓
              • 要素类型不兼容:求解器专有要素无法在中立格式中表示
              • 材料模型差异:同名的材料模型内部实现可能不同
              • 边界条件重新定义:很多情况下需要手工重新设置
              • 结果数据比较:输出变量定义差异(节点值 vs. 单元值、积分点值)

              • 🧑‍🎓

                哦,这样啊!不同工具间模型转换有这样的事。


                许可模式

                🧑‍🎓

                "许可模式"听过,好像没完全理解…


                工具许可特点
                商用FEA节点锁定/浮动费用高但有官方支持
                OpenFOAMGPL免费但支持需付费
                COMSOL节点锁定/浮动按模块购买
                Code_AsterGPLEDF开发的开源求解器

                选择指引

                🧑‍🎓

                结局怎么选择,判断标准给我讲讲?


                🎓

                选择自适应网格加密(AMR)工具需考虑以下几点:


                🎓
                • 分析规模:可扩展至数万~数亿自由度
                • 物理模型:所需构成关系、要素类型的支持情况
                • 工作流:与CAD的联动、自动化的容易程度
                • 成本:初期投资 + 年度维护 + 教育成本
                • 支持:技术支持的质量和响应速度


                • 🧑‍🎓

                  自适应网格加密(AMR)整体思路我理解了!明天在实务中试试。


                  🎓

                  好啊,继续加油!实际动手才是最好的学习。有疑惑的地方随时来问。


                  验证数据可视化

                  定量显示理论值与计算值的对比。以相对误差5%以内为合格基准。

                  评价项目理论值/参考值计算值相对误差 [%]判定
                  最大位移1.0000.998
                  0.20
                  PASS
                  最大应力1.0001.015
                  1.50
                  PASS
                  固有频率(1阶)1.0000.997
                  0.30
                  PASS
                  反力合计1.0001.001