有限元法分析结果的误差分析

     有限元法分析一般包括四个步骤：物理模型的简化、数学模型的程序化、计算模型的数值化和计算结果的分析。每一个步骤在操作过程中都或多或少地引入了误差，这些误差的累积最终可能会对计算结果造成灾难性的影响，进而蒙蔽我们的认识和判断。

     第一步，物理模型的简化，主要有几何实体、连接 / 装配关系、环境边界条件和材料特性的简化，进而构建数学模型。这些简化或者说假设，是必要的，也是必须的，但是也由此在模型中引入了理想化误差（ idealization error ）。有些理想化误差是非良性奇异的，比如几何实体简化时细节部位上忽略小的圆 / 倒角，连接 / 装配关系简化时忽略焊缝和螺栓连接等，往往导致模型发生结构方面（诸如 L 形截面的角点）的奇异，即结构奇异（奇异的数学定义是在某一点处导数无穷）；有些理想化误差是良性奇异的，比如边界条件简化时添加集中载荷和孤立点约束，导致模型发生边界条件的奇异，即边界奇异；其它理想化误差，比如几何实体简化时三维壳 / 面体简化为二维壳 / 面、三维梁简化为一维梁，边界条件简化时非均匀温度场和压力场简化为均匀温度场和压力场等，只会影响计算结果的准确度，不会引发计算结果方面的数值奇异，即应力奇异和位移奇异等。理想化误差是在有限元法分析开始之前引入的，因此我们不可能通过改进有限元分析技术来达到消除其的目的，而只能通过修改数学模型本身来实现消除其的目的。
   

  第二步，数学模型的程序化，主要有几何实体的单元离散、单元网格的装配连接、模型环境边界条件的添加，进而构建计算模型。几何实体的离散，和单元类型（形状和精度）、单元尺寸以及分网方式的选择有关，不可避免地会引入离散化误差（ discretization error ）。离散化误差，是根植于有限元法分析本身的，因此只能通过改进有限元分析技术或者技巧来尽力消除 / 减小这方面的误差，比如采用规则化的单元形状避免单元在形状上产生奇异（即单元奇异）、提高单元精度和增加网格密度减小计算方面的误差等方法。单元网格的装配连接一般采用 MPC 多点约束法，因而会引入人为误差（ artificial error ），这方面误差的消除更多是需要长期计算经验的积累。模型环境边界条件的添加，其误差影响依赖于第一步的理想化简化。
    

 第三步，计算模型的数值化，主要是用数值计算方法（程序求解器）求解、逼近真实的解析值，因而必然存在数值化误差（ numerical error ）。数值计算方法的精度（非人为可控）越高，计算结果的误差就越小，但计算的工作量也越大。实际考虑到计算精度和计算资源的利用，必然要做一个适当的统一。
     

第四步，计算结果的分析，主要是利用数值计算结果来分析、评判，或预知真实的物理模型，由此也存在着认知误差（ recognized error ）。认知误差的消除，一方面需要真实物理试验的指导，另一方面依赖于分析人员的工程经验和认知能力。同时，不要忘记了我们的前提假设，即第一步物理模型的简化，或假设。

转载请注明精川材料检测地址：www.jctest.vip