克里金法(高斯过程回归)代理模型

分类:分析 | 整合版 2026-04-06
CAE visualization for kriging surrogate theory - technical simulation diagram
克里金法(高斯过程回归)代理模型

克里金法(高斯过程回归)代理模型的理论基础

概述

🧑‍🎓

老师!今天是克里金法(高斯过程回归)代理模型的讨论是吗?它是什么样的?


🎓

同时推估预测值和预测不确定性的代理模型。作为计算成本高昂的仿真的替代品用于不确定性分析。




支配方程




$$ \hat{Y}(\mathbf{x}) = \boldsymbol{\mu} + \mathbf{r}^T\mathbf{R}^{-1}(\mathbf{y}-\boldsymbol{\mu}) $$
$$ \hat{\sigma}^2(\mathbf{x}) = \sigma^2\left(1-\mathbf{r}^T\mathbf{R}^{-1}\mathbf{r}\right) $$



🧑‍🎓

等等,克里金法,也就是说在这种情况下也能用吗?


离散化方法

🧑‍🎓

这个方程在计算机上具体怎么求解?


🎓

使用有限元法(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 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的话题吧。内容是什么?


🎓

作为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
    通过
    最大应力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%: 需检讨

    克里金法(高斯过程回归)代理模型的数值计算方法

    数值方法详解

    🧑‍🎓

    克里金法(高斯过程回归)代理模型具体用什么算法求解?




    离散化的形式化



    🎓

    形状函数 $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{Y}(\mathbf{x}) = \boldsymbol{\mu} + \mathbf{r}^T\mathbf{R}^{-1}(\mathbf{y}-\boldsymbol{\mu}) $$
    $$ \hat{\sigma}^2(\mathbf{x}) = \sigma^2\left(1-\mathbf{r}^T\mathbf{R}^{-1}\mathbf{r}\right) $$

    🧑‍🎓

    唔…仅看公式没那么直观…这是什么意思?


    🎓

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



    $$ [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不收敛时需重新调整设置
        内存模式内核尽可能采用

        低阶单元

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

        高阶单元

        用相同网格达到更高精度。计算成本增加,但所需单元数常减少。

        牛顿-拉夫逊法

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

        克里金法(高斯过程回归)代理模型的实务应用

        实践指南

        🧑‍🎓

        老师,请告诉我"实践指南"!


        🎓

        讲解克里金法(高斯过程回归)代理模型的实务分析流程与注意点。



        分析流程

        🧑‍🎓

        最开始该怎么做?先从哪儿开始?


        🎓

        1. 预处理(Pre-processing)

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

        🎓

        2. 求解(Solving)

        • 求解器设置(求解法、收敛基准、输出控制)
        • 作业投入与计算执行
        • 收敛监视

        🎓

        3. 后处理(Post-processing)

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


        网格生成最佳实践

        🧑‍🎓

        怎样判断网格的质量好坏?



        单元品质指标

        🧑‍🎓

        请介绍"单元品质指标"!


        指标理想值许可范围影响
        纵横比1.0< 5.0精度降低
        Jacobian比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的话题吧。内容是什么?


            🎓

            作为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
              通过
              最大应力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%: 需检讨

              克里金法(高斯过程回归)代理模型的软件比较

              商用工具对比

              🧑‍🎓

              有各种各样的软件吧?请介绍一下各自的特点!


              🎓

              讲解克里金法(高斯过程回归)代理模型对应的主要商用CAE工具的功能对比和各产品历史背景。



              支持工具清单

              🧑‍🎓

              那么做克里金法(高斯过程回归)代理模型需要什么软件?


              工具名称开发方/现在主要文件格式
              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的话题吧。内容是什么?


              🎓

              作为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

              🧑‍🎓

              接下来是Ansys Fluent的话题吧。内容是什么?


              🎓

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

              现在所属:Ansys Inc.


              🧑‍🎓

              哇,结构解析的故事太有趣了!请告诉我更多!


              功能对比矩阵

              🧑‍🎓

              预算和时间都有限,最划算的是哪个?


              功能NastranAbaqusAnsys MechanicalFluent
              基本功能
              高级功能
              自动化/脚本