侧面碰撞模拟中效率与精度影响因素研究.doc

侧面碰撞模拟中效率与精度影响因素研究 基于动态显式非线性有限元分析的计算机仿真，与隐式算法不同，其算法是有条件稳定的，包括由简化单元计算使用的单点高斯积分引起的沙漏模态与应力波传递带来的数值震荡均需进行有效控制[41]，同时，显式积分的稳定性还受临界时间步长的控制。为了确保仿真模拟的有效性，有必要探讨网格单元尺寸大小、类型以及积分点，时间步长，接触和摩擦等因素对仿真模拟效果的影响规律。 3.1单元尺寸影响 单元尺寸主要影响模型的真实性和极限时间步长。越小的单元尺寸对压塌变形的表达越真实，过粗的网格则会导致失真的结果。但越小的单元尺寸要求的时间步长越小，因而导致计算效率下降越多。因此，在汽车侧面碰撞仿真计算中，应在确保正确表达压塌变形的前提下，尽可能选择大的单元尺寸。 已有文献对车身用薄壁梁的研究表明，沿轴向碰撞的直梁结构在局部失稳后将产生压塌失效过程，在该过程中，直梁的边会沿着半径为r的圆弧形曲线逐渐折叠，根据Werzbicki[42]对薄壁直梁件的研究，平均的折叠半径近似为: （3.1） C是截面的宽度，t是壁厚。为了充分描述出结构在碰撞过程中的折叠变形，单元的尺寸应该小于折叠圆弧长的一半，即单元边长: （3.2） 总结起来，我们认为小的单元比较容易捕捉到结构碰撞中产生的小“皱褶”，模拟的效果也越好，而粗网格单元刚性较强，模拟变形量相对较小。但过小的单元尺寸，会增加很多计算时间，所以在进行汽车侧面碰撞仿真时要综合考虑计算精度和计算效率两方面，对不同的碰撞区域选择不同的单元尺寸。例如大的变形区选择小点的单元尺寸，而非主要变形区和非变形区则选择大的单元尺寸，这样既考虑了精度，同时在可信的精度下又保证了效率。 3.2单元类型选择 对于有限元仿真分析，在建立有限元模型时选择什么样的单元进行分析是非常重要的。选择单元的类型，主要考虑结构的几何形状，其次是分析目的，加载条件和经济性。汽车结构部件主要由薄壳板件冲压后经焊接拼装而成，因此在汽车侧面碰撞分析中需要用薄壳单元来进行空间离散，薄壳单元的选择和算法对于计算的效率与正确性有着重要的意义。在非线性动力问题的研究发展中，产生了多种薄壳单元算法，LS-DYNA970为壳单元提供了16种算法，它们之间的主要差别在于采用的积分方法不同以及能否控制沙漏。对相同单元，不同的壳单元算法导致的计算时间差别会很大。本文选择了Belytschko-Tsay(BT)单元和Belytschko-Wong-Chiang(BWC)单元以及Hughes-Liu(HL)三种单元算法来进行对比，同时，选用BT单元为理论基础，但是在刚度计算时，采用BWC的刚度计算理论，在这里称为BT-BWC，该小节中将会对上述的四种单元算法进行对比分析，采用左侧B柱有限元模型为研究对象，对四种不同单元类型的计算时间进行了比较。 表3.1 不同单元类型的计算时间 单元类型 计算时间(s) BT 860 BWC 1252 BT-BWC 980 HL 1100 表3.1为四种单元算法的计算时间表，表中显示了BT单元的计算时间是最少的，BWC的计算时间最长，其次是HL。BWC单元的计算时间比BT计算时间增加了45.6%，HL单元的计算时间比BT计算时间增加了28.8%，BT-BWC单元的计算时间比BT计算时间增加了14%，从计算的精度和计算的效率两方面的来考虑，单元类型选择为BT单元类型，但是在control-shell中关于翘曲刚度的计算选用BWC的计算法。 3.3网格密度的影响 进行网格的划分就是模型的离散化，是将连续结构体离散成由各种单元组成的有限元模型，为模拟计算分析做好准备工作。离散后单元与单元之间通过节点相互连接，单元的属性及数量等应该根据所要解算问题的具体情况进行确定。为了提高计算精度和得到正确的碰撞仿真结果，对网格的总体加密有时不能完全解决问题，网格密度的分布也是重要的影响因素。是在同样的边界条件和初始条件下，不同划分密度产生不同变形的情况。一般情况下，单元划分越细则对物体的描述就越精确，计算所得的变形结果也就越接近实际情况，但是，不能为了所谓的精确度而盲目细分模型，因为这样会导致计算模型单元数量巨大、单元时间步长小，从而使计算量急剧增加，很有可能无法正常完成计算。网格单元划分要根据所要研究的问题的规模、周期、侧重点以及所拥有的计算机硬件等条件权衡考虑，力求计算精度与计算时间的最佳搭配。 不同的网格密度对能量的吸收情况不同。粗网格的计算内能偏大，刚性较强，因此在相同的受力条件下，细网格区将先于粗网格区进入屈服压塌状态。如果网格划分不当，在受力梯度不大时就会得出错