解检验(Solution Verification)

分类:分析 | 综合版 2026-04-06
CAE visualization for solution verification theory - technical simulation diagram
解检验(Solution Verification)

解检验(Solution Verification)的理论基础

概要

🧑‍🎓

老师!今天是讲解检验(Solution Verification)的话题吧?它是什么东西呢?


🎓

特定问题的离散化误差估算。Richardson外推法和GCI进行不确定性定量化。




支配方程




$$ \hat{f}_h = f_h + \frac{f_1 - f_2}{r^p - 1} $$
$$ \text{GCI}_{fine} = \frac{F_s|\varepsilon_{12}|}{r_{12}^p - 1} $$




离散化手法

🧑‍🎓

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


🎓

使用有限元法(FEM)进行空间离散化。组建要素刚度矩阵,构造全局刚度方程。


🎓

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




矩阵求解算法

🧑‍🎓

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


🎓

通过直接法(LU分解Cholesky分解)或迭代法(CG法GMRES法)求解线性方程组。大规模问题中前处理迭代法效果显著。



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

也就是说有限元法那块搞不好,后面会吃大亏对吧。我记住了!


商用工具中的实现

🧑‍🎓

那进行解检验(Solution Verification)能用什么软件呢?


工具名称开发者/现属主要文件格式
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。2017年MSC被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中立CAD符合ISO 10303的3D CAD数据交换格式。支持形状+PMI。
IGES.igs/.iges中立CAD早期CAD数据交换规范。曲面数据兼容性有问题。正逐步被STEP取代。
VTK.vtk/.vtu可视化Visualization Toolkit格式。用于ParaView等。
🎓

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


🧑‍🎓

啊……看似简单的格式,实际上内涵很深呢。


实务注意事项

🧑‍🎓

有没有"教科书上没写"但工程中很重要的经验呢?


🎓

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


🎓
  • 网格依赖性验证:至少用3种网格密度确认收敛性
  • 边界条件合理性:设置物理上有意义的约束条件
  • 结果验证:与解析解、实验数据、已知基准问题的对比


  • 🧑‍🎓

    现在对解检验(Solution Verification)的全貌有把握了!明天起要在实务中意识到这些。


    🎓

    很好!实际动手实践是最好的学习方法。有不懂的地方随时来问。


    验证数据的可视化

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

    评估项目理论值/参照值计算值相对误差 [%]判定
    最大变位1.0000.998
    0.20
    合格
    最大应力1.0001.015
    1.50
    合格
    固有振动数(1阶)1.0000.997
    0.30
    合格
    反力合计1.0001.001
    0.10
    合格
    能量守恒1.0000.999
    0.10
    合格

    判定基准:相对误差 < 1%: 优秀、1~5%: 可接受、> 5%: 需检查

    解检验(Solution Verification)的数值计算方法

    数值方法的详细内容

    🧑‍🎓

    具体用什么算法来求解解检验(Solution Verification)呢?



    🧑‍🎓

    前辈说过"解检验得好好做",现在明白是怎回事了。


    离散化的定式化



    🎓

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

    基本方程的离散形式


    🎓

    用公式表达就这样。


    $$ \hat{f}_h = f_h + \frac{f_1 - f_2}{r^p - 1} $$
    $$ \text{GCI}_{fine} = \frac{F_s|\varepsilon_{12}|}{r_{12}^p - 1} $$

    🧑‍🎓

    光看公式我还是有点摸不着头脑…表示什么意思呢?


    🎓

    连续体的支配方程离散化后,得到这样的代数方程组:



    $$ [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)$ 阶降低(光滑解的情况)


        🧑‍🎓

        网格细分看似简单,实际上内涵很深呢。


        求解器配置的建议

        🧑‍🎓

        具体用什么算法来求解解检验(Solution Verification)呢?


        参数建议值说明
        迭代法收敛判定$10^{-6}$残差范数基准
        前处理方法ILU(0) or AMG问题规模决定
        最大迭代次数1000不收敛时需重审设置
        内存模式In-core尽可能采用

        低阶要素

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

        高阶要素

        同一网格下实现更高精度。计算成本增加,但所需单元数通常较少。

        牛顿-拉夫逊法

        非线性问题的标准方法。在收敛半径内二次收敛。以 $||R|| < \epsilon$ 判定收敛。

        时间积分

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

        验证数据的可视化

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

        评估项目理论值/参照值计算值相对误差 [%]判定
        最大变位1.0000.998
        0.20
        合格
        最大应力1.0001.015
        1.50
        合格
        固有振动数(1阶)1.0000.997
        0.30
        合格
        反力合计1.0001.001
        0.10
        合格
        能量守恒1.0000.999
        0.10
        合格

        判定基准:相对误差 < 1%: 优秀、1~5%: 可接受、> 5%: 需检查

        解检验(Solution Verification)的实务应用

        实践指南

        🧑‍🎓

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


        🎓

        解检验(Solution Verification)的实务分析流程和注意要点。


        🧑‍🎓

        前辈说过"解检验得好好做",现在明白是怎回事了。


        分析流程

        🧑‍🎓

        从最初一步开始教我!我该怎么开始呢?


        🎓

        1. 前处理 (Pre-processing)

        • 导入CAD数据并简化形状
        • 定义材料特性
        • 网格生成(单元类型·尺寸的决定)
        • 设定边界条件与荷载条件

        🎓

        2. 求解 (Solving)

        • 求解器设置(求解法、收敛基准、输出控制)
        • 任务提交与计算执行
        • 收敛状态监视

        🎓

        3. 后处理 (Post-processing)

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


        网格生成的最佳实践

        🧑‍🎓

        网格好坏怎样判断呢?



        要素品质指标

        🧑‍🎓

        请讲讲"要素品质指标"!


        指标理想值允许范围影响
        长宽比1.0< 5.0精度降低
        Jacobian比1.0> 0.3要素退化
        翘曲< 15°精度降低
        倾斜度< 45°收敛性恶化
        Taper比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。2017年MSC被Hexagon AB收购。

            现属单位:MSC Nastran(Hexagon)、NX Nastran(Siemens Digital Industries Software)



            Abaqus FEA (SIMULIA)

            🧑‍🎓

            前辈的解释好好理解了!工具名称的疑惑解开了。


            常见失败与对策

            🧑‍🎓

            初学者容易犯什么失败呢?想提前知道!


            症状原因对策
            计算不收敛网格品质不良、不适当的边界条件改善网格、检查约束条件
            应力异常大应力奇点、网格依赖回避奇点、局部网格细化
            位移非现实材料常数错误、单位混用确认输入数据
            计算时间过长不必要的细化、低效求解法网格优化、并列计算

            质量保证检查清单

            🧑‍🎓

            教科书上没写但现场很重要的"秘诀"吗?


            🎓
            • 在3种以上网格密度下确认了收敛性吗
            • 验证了力的平衡(反力合计)吗
            • 结果在物理的合理范围内吗
            • 与已知理论解或基准问题做过比较吗


            • 🧑‍🎓

              现在对解检验(Solution Verification)的全貌有把握了!明天起要在实务中意识到这些。


              🎓

              很好!实际动手实践是最好的学习方法。有不懂的地方随时来问。


              验证数据的可视化

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

              评估项目理论值/参照值计算值相对误差 [%]判定
              最大变位1.0000.998
              0.20
              合格
              最大应力1.0001.015
              1.50
              合格
              固有振动数(1阶)1.0000.997
              0.30
              合格
              反力合计1.0001.001
              0.10
              合格
              能量守恒1.0000.999
              0.10
              合格

              判定基准:相对误差 < 1%: 优秀、1~5%: 可接受、> 5%: 需检查

              解检验(Solution Verification)的软件比较

              商用工具比较

              🧑‍🎓

              有各种各样的软件对吧?各自特点请讲讲!


              🎓

              支持解检验(Solution Verification)的主要商用CAE工具的功能对比与各产品的历史背景。


              🧑‍🎓

              前辈说过"解检验得好好做",现在明白是怎回事了。


              支持工具一览

              🧑‍🎓

              那进行解检验(Solution Verification)能用什么软件呢?


              工具名称开发者/现属主要文件格式
              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。2017年MSC被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开发的开源求解器

                选择指导

                🧑‍🎓

                最后应该选哪个,判断基准是什么?


                🎓

                进行解检验(Solution Verification)时工具选择要考虑这些:


                🎓
                • 分析规模:对