材料数据的验证和不确定性

分类: 分析 | 统一版 2026-04-06
CAE visualization for material data validation theory - technical simulation diagram
材料数据的验证和不确定性

材料数据验证和不确定性的理论基础

概览

🧑🎓

老师!今天是关于材料数据验证和不确定性的讲座吗?是什么内容?


🎓

材料试验数据的质量确认。温度应变速率依赖性、数据外推风险、分散性考虑。




控制方程




$$ E = \frac{\sigma}{\varepsilon} \text{ (来自拉伸试验数据)} $$
$$ \sigma_Y = \sigma_{0.2\%} \pm u_{material} $$



🧑🎓

听到这里,我终于理解了为什么材料数据验证和不确定性这么重要!


离散化方法

🧑🎓

这个方程在计算机上具体如何求解?


🎓

使用有限元法(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)| $$

    基础方程的离散形式


    🎓

    用式子表示就是这样。


    $$ E = \frac{\sigma}{\varepsilon} \text{ (来自拉伸试验数据)} $$
    $$ \sigma_Y = \sigma_{0.2\%} \pm u_{material} $$

    🧑🎓

    嗯,只有式子的话很难直观理解…这表示什么呢?


    🎓

    将连续体的控制方程离散化后,得到以下代数方程组:



    $$ [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尽可能采用

        低阶单元

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

        高阶单元

        同一网格下达到更高精度。计算成本增加,但必需的单元数往往会减少。

        牛顿拉夫逊法

        非线性问题的标准方法。收敛半径内具有2次收敛。$||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精度下降
        雅可比行列式比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)

            🧑🎓

            老师的讲解很好懂!工具名的疑惑消散了。


            常见失败和应对措施

            🧑🎓

            初学者容易犯什么失败?想预先知道!


            症状原因对策
            计算不收敛网格质量不足、边界条件不当网格改善、约束条件重新审视
            应力异常大应力奇点、网格依赖奇点回避、局部网格细化
            变位不现实材料常数错误、单位系统不一致确认输入数据
            计算时间过长不必要细化、求解方法效率低网格优化、并列计算

            质量保证检查清单

            🧑🎓

            有没有教科书里不会讲的"现场智慧"那样的东西?


            🎓
            • 用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 Fluent"


              🎓

              1970年由Swanson Analysis Systems Inc. (SASI)开发。基于APDL(Ansys参数化设计语言)。

              目前所属: Ansys Inc.



              Ansys Fluent

              🧑🎓

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


              🎓

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

              目前所属: Ansys Inc.


              🧑🎓

              哦,结构分析的话题,超级有意思!请给我更多讲解。


              功能对比矩阵

              🧑🎓

              预算和时间都有限,成本效益最强的是哪个?


              功能NastranAbaqusAnsys MechanicalFluent
              基础功能