层级妥当性确认

分类:解析 | 综合版 2026-04-06
CAE visualization for hierarchical validation theory - technical simulation diagram
层级妥当性确认

层级妥当性确认的理论基础

概要

🧑‍🎓

老师!今天是关于层级妥当性确认的话题吧?那是什么呢?


🎓

在单体→组件→子系统→系统的层级中进行验证与妥当性确认。




支配方程




$$ E_{system} \leq \sum_{i} E_{subsystem,i} + E_{coupling} $$
$$ \text{CL}(E) = P(|S-D| < \delta) $$




离散化手法

🧑‍🎓

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


🎓

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


🎓

进行弱形式(变分形式)的转换,使用试验函数和形状函数,采用Galerkin法进行定式化。单元类型的选择(低阶单元 vs. 高阶单元完全积分 vs. 降阶积分)直接关系到解的精度和计算成本的折衷。




矩阵求解算法

🧑‍🎓

矩阵求解算法,具体是指什么呢?


🎓

通过直接法(LU分解Cholesky分解)或迭代法(CG法GMRES法)求解线性方程组。对于大规模问题,预处理迭代法最有效。



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

也就是说,在有限元法这一步如果不认真,之后会吃大亏啊。我铭记在心!


商用工具中的实现

🧑‍🎓

那么,进行层级妥当性确认需要什么样的软件呢?


工具名称开发方/现在主要文件格式
MSC Nastran / NX NastranMSC Nastran(Hexagon)、NX NastranSiemens 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的故事吧。什么内容呢?


🎓

作为NASA结构分析(NASTRAN)在1960年代开发。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参数化设计语言)。

现在的所属:Ansys Inc.


🧑‍🎓

哦~,结构解析的故事,超级有趣!请多讲点。


文件格式和互操作性

🧑‍🎓

不同软件间传递数据时有什么注意点吗?


格式扩展名类型概要
STEP.stp/.step中立CADISO 10303遵循的3D CAD数据交换格式。形状+PMI对应。
IGES.igs/.iges中立CAD早期的CAD数据交换规范。曲面数据的互通性存在问题。向STEP过渡进行。
VTK.vtk/.vtu可视化Visualization Toolkit格式。用于ParaView等。
🎓

在不同求解器间转换模型时,需要注意单元类型的对应关系、材料模型的相容性、荷载和边界条件的表示差异。尤其是高阶单元或特殊单元(粘聚单元、用户定义单元等),在求解器间通常无法直接转换。


🧑‍🎓

原来如此……格式看似简单,其实深藏玄机呢。


实务上的注意事项

🧑‍🎓

教科书上没有的"现场智慧"有没有呢?


🎓

网格收敛性的确认、边界条件的妥当性验证、材料参数的敏感性分析都超级重要。


🎓
  • 网格依存性的验证:至少用3个网格密度水平确认收敛性
  • 边界条件的妥当性:设置物理上有意义的约束条件
  • 结果的验证:与理论解、实验数据、已知基准问题的比较


  • 🧑‍🎓

    层级妥当性确认的整体图像掌握了!明天开始在实务中注意。


    🎓

    嗯,你的状态不错!亲手动手做是最好的学习。有不明白的地方随时来问。


    验证数据的可视化

    定量显示理论值与计算值的比较。合格标准为误差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%: 需检讨

    层级妥当性确认的数值计算手法

    数值手法的详细

    🧑‍🎓

    具体用什么算法求解层级妥当性确认呢?



    🧑‍🎓

    等等等等,层级妥当性确认的话,也可以用于这种情况吗?


    离散化的定式化



    🎓

    使用形状函数 $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)| $$

    基础方程的离散形


    🎓

    用数式表示就是这样。


    $$ E_{system} \leq \sum_{i} E_{subsystem,i} + E_{coupling} $$
    $$ \text{CL}(E) = P(|S-D| < \delta) $$

    🧑‍🎓

    嗯,只有式子的话摸不着……表示什么呀?


    🎓

    将连续体的支配方程离散化,得到以下代数方程系:



    $$ [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)$ 的数量级递减


        🧑‍🎓

        原来……细分网格看似简单,但其实学问很深呢。


        求解器设置的推荐

        🧑‍🎓

        具体用什么算法求解层级妥当性确认呢?


        参数推荐值备注
        迭代法的收敛判定$10^{-6}$残差模数基准
        预处理手法ILU(0) or AMG根据问题规模
        最大迭代次数1000不收敛时重新检查设置
        内存模式In-core尽可能

        低阶单元

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

        高阶单元

        在相同网格下实现更高精度。计算成本虽然增加,但必要的单元数通常会减少。

        牛顿·拉夫逊法

        非线性问题的标准手法。在收敛半径内二阶收敛。以 $||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%: 需检讨

        层级妥当性确认的实务应用

        实践指南

        🧑‍🎓

        老师,请讲讲"实践指南"!


        🎓

        讲解层级妥当性确认的实务解析流程和注意点。


        🧑‍🎓

        师兄说"层级妥当性确认的实这一步一定要好好做",现在我明白意思了。


        解析流程

        🧑‍🎓

        从头开始讲讲!首先做什么?


        🎓

        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 NastranSiemens 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的故事吧。什么内容呢?


            🎓

            作为NASA结构分析(NASTRAN)在1960年代开发。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个以上网格密度水平确认了收敛性吗
            • 验证了力的平衡(反力合计)吗
            • 确认结果处于物理上合理的范围内吗
            • 与已知理论解或基准问题进行了比较吗


            • 🧑‍🎓

              层级妥当性确认的整体图像掌握了!明天开始在实务中注意。


              🎓

              嗯,你的状态不错!亲手动手做是最好的学习。有不明白的地方随时来问。


              验证数据的可视化

              定量显示理论值与计算值的比较。合格标准为误差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%: 需检讨

              层级妥当性确认的软件比较

              商用工具比较

              🧑‍🎓

              有各种各样的软件吧?每个的特点讲讲!


              🎓

              支持层级妥当性确认的主要商用CAE工具的功能比较,以及各产品的历史背景。


              🧑‍🎓

              等等等等,层级妥当性确认的话,也可以用于这种情况吗?


              对应工具列表

              🧑‍🎓

              那么,进行层级妥当性确认需要什么样的软件呢?


              工具名称开发方/现在主要文件格式
              MSC Nastran / NX NastranMSC Nastran(Hexagon)、NX NastranSiemens 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的故事吧。什么内容呢?


              🎓

              作为NASA结构分析(NASTRAN)在1960年代开发。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参数化设计语言)。

              现在的所属:Ansys Inc.



              Ansys Fluent

              🧑‍🎓