稀疏PCE与LAR

分类:分析 | 综合版 2026-04-06
CAE visualization for sparse pce theory - technical simulation diagram
稀疏PCE与LAR

稀疏PCE与LAR的理论基础

概述

🧑‍🎓

老师!今天是讲稀疏PCE与LAR的对吧?这是什么东西呢?


🎓

高维问题中PCE的维度诅咒对策。通过最小角回归等稀疏回归方法仅选择有效项。



🧑‍🎓

我明白了。那么如果高维问题处理得当的话,基本上就没问题了对吧?


控制方程




$$ Y \approx \sum_{\boldsymbol{\alpha}\in\mathcal{A}} c_{\boldsymbol{\alpha}}\Psi_{\boldsymbol{\alpha}}(\mathbf{X}) $$
$$ \hat{\boldsymbol{c}} = \arg\min\|\boldsymbol{\Psi c} - \mathbf{y}\|_2^2 + \lambda\|\boldsymbol{c}\|_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)超大规模
🧑‍🎓

也就是说在有限元法的阶段中做不好的话,后面就会很吃亏啊。铭记在心!


商业工具中的实现

🧑‍🎓

那么做稀疏PCE与LAR需要什么软件呢?


工具名称开发商/现在主要文件格式
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数据交换格式。支持形状+产品制造信息。
IGES.igs/.iges中立CAD早期CAD数据交换标准。曲面数据互操作性存在问题。正向STEP迁移。
VTK.vtk/.vtu可视化可视化工具包格式。用于ParaView等。
🎓

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


🧑‍🎓

原来格式看起来简单,其实超级复杂啊。


实务上的注意事项

🧑‍🎓

有没有教科书里讲不到的「现场智慧」呢?


🎓

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


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


  • 🧑‍🎓

    哇,稀疏PCE与LAR原来这么深奥啊……不过在老师的讲解下终于理顺了思路!


    🎓

    嗯,不错!真正的学习就是要动手实践。有疑问随时来问。


    验证数据的可视化

    定量展示理论值与计算值的对比。误差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%: 需检查

    稀疏PCE与LAR的数值计算方法

    数值方法的详细说明

    🧑‍🎓

    具体怎样的算法来求解稀疏PCE与LAR?




    离散化的数学表述



    🎓

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

    基础方程的离散形式


    🎓

    用式子表示就是这样:


    $$ Y \approx \sum_{\boldsymbol{\alpha}\in\mathcal{A}} c_{\boldsymbol{\alpha}}\Psi_{\boldsymbol{\alpha}}(\mathbf{X}) $$
    $$ \hat{\boldsymbol{c}} = \arg\min\|\boldsymbol{\Psi c} - \mathbf{y}\|_2^2 + \lambda\|\boldsymbol{c}\|_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)$ 的数量级误差递减(光滑解的情况下)


        🧑‍🎓

        原来网格细分看似简单,其实深度很大啊。


        求解器设置推荐

        🧑‍🎓

        具体怎样的算法来求解稀疏PCE与LAR?


        参数推荐值说明
        迭代法收敛判定$10^{-6}$残差范数标准
        预处理方法ILU(0) 或 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%: 需检查

        稀疏PCE与LAR的实际应用

        实践指南

        🧑‍🎓

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


        🎓

        讲解稀疏PCE与LAR的实务分析流程和注意事项。



        分析流程

        🧑‍🎓

        从第一步开始教我!怎么开始?


        🎓

        1. 预处理 (前处理)

        • CAD数据导入和几何简化
        • 材料特性定义
        • 网格划分(单元类型、尺寸的决定)
        • 边界条件和荷载条件的设置

        🎓

        2. 求解 (求解)

        • 求解器设置(求解方法、收敛准则、输出控制)
        • 作业投入和计算执行
        • 收敛性监控

        🎓

        3. 后处理 (后处理)

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


        网格划分的最佳实践

        🧑‍🎓

        网格好坏怎样判断?



        单元质量指标

        🧑‍🎓

        「单元质量指标」请给我讲讲!


        指标理想值可接受范围影响
        长宽比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个以上等级吗
            • 验证了力平衡(反力合计)吗
            • 检查过结果在物理合理范围内吗
            • 和已知理论解或基准问题比对了吗


            • 🧑‍🎓

              哇,稀疏PCE与LAR原来这么深奥啊……不过在老师的讲解下终于理顺了思路!


              🎓

              嗯,不错!真正的学习就是要动手实践。有疑问随时来问。


              验证数据的可视化

              定量展示理论值与计算值的对比。误差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%: 需检查

              稀疏PCE与LAR软件比较

              商业工具比较

              🧑‍🎓

              有各种软件吧?各自的特点给我讲讲!


              🎓

              比较支持稀疏PCE与LAR的主要商业CAE工具的功能,详述各产品的历史背景。



              支持工具列表

              🧑‍🎓

              那么做稀疏PCE与LAR需要什么软件呢?


              工具名称开发商/现在主要文件格式
              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

              🧑‍🎓

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


              🎓

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

              现在的所属:Ansys Inc.


              🧑‍🎓

              哇,结构分析的内容超有意思!请再给我讲讲。