侧面碰撞仿真基本理论与方法.docVIP

侧面碰撞仿真基本理论与方法 2.1引言 对于碰撞过程的研究方法，基本可以分为实验法和仿真法两种。80年代之前，当有限元应用于汽车的碰撞分析时还只能计算非常简化的车辆模型，不仅单元数量少，而且模型与实体之间差别巨大；80年代之后，随着有限元理论和方法的不断完善和计算机技术的飞速发展，汽车碰撞的精确数值分析与计算己成了现实，并且正在逐步取代部分试验室的工作。汽车碰撞分析的有限元法本身也是一个复杂的课题，它涉及到众多的学科和领域，尤其是有限元理论、计算数学、计算力学、弹塑性力学、计算机图形学等，并且这些理论本身也在不断发展和完善。 汽车侧面碰撞是一个动态的大位移和大变形过程，接触和高速冲击载荷影响着碰撞的整个过程，系统具有几何非线性、材料非线性和边界非线性等多重非线性。因此，汽车侧面碰撞的数值计算一般涉及非线性动态有限元，时间积分法，碰撞结构单元等。本文是采用了有限元仿真的办法。因此，在下面就着重介绍一下这种办法基本理论与方法。 2.2有限元简介 有限元法是大型复杂结构或多自由度体系分析的有力工具，近20年来已广泛地用于：工程结构、传热、流体运动、电磁等连续介质的分析，并在气象、地球物理、医学等领域得到应用和发展。电子计算机的出现和发展，使有限元法在许多实际问题的应用变成现实，并且有广阔的前景。有限元法是关于边界问题的、建立在待定场函数离散化基础上的一种数值方法。这种方法设想把所研究的连续区域分割成有限数量的子区域——称为有限单元，同时把待定的场函数用分别在各单元中连续的近似场函数代替。在每个单元中的近似函数完全由该单元边界上的节点的近似函数值决定。这些函数值如果通过某种平均的方法，例如变分法求得，那么就能获得这个问题的分区近似解。古典的变分法或者称直接法就是 Ritz.Galerkin 方法，它源于最小势能原理。有限元法属于 Ritz.Galerkin 方法的范畴，但又是这种古典变分法的发展和改进。这里的关键就是选取分块多项式函数类SN作为真实函数类的空间。具体说，就是先把区域Ω剖分为一系列单元，然后在每一个单元内作多项式插值，使它们在相邻单元的公共面（边）上满足一定的连续性条件，以保证用这种函数组成的有限维空间SN是实际函数类H的空间。从这个意义上讲，有限元法的基础是变分原理和分块多项式插值[]。用有限元法近似求解变分问题时，产生误差的原因很多。首先是用子空间 Sh 逼近容许函数类空间 H，即用投影逼近真解；其次是区域逼近（即用多角形代替 Ω）和边界条件的近似；再次是计算刚度矩阵和载荷向量时采用数值积分；最后还有求解代数方程组的舍入误差等等。但一般来说还是第一种原因造成的逼近误差[]。 图 2.1 物体变形过程描述 设质点 a 在t = 0时刻的位置坐标为xa（α =1，2，3）。这时由质点构成的物体图形称为初始态图形。设在此后的任意时刻 t，质点移动到另一位置，由于质点间的相互位置关系发生了变化，物体因此改变了形状。在t时刻，物体经过变形后，质点 a 的位置坐标为xi（i=1，2，3）。它是坐标xa和时间t的函数，即 （2.1） 2.3.2虚位移原理和虚应力原理 1、虚位移原理 虚位移原理实质上是平衡方程的弱形式，它是以积分的形式描述变形物体的平衡，是与平衡方程等价的。也就是说，由虚位移原理可以导出平衡方程，反之，由平衡方程也可以导出虚位移原理。 如果对于任一虚位移，有 　　 　　　　　　　　　　　　　　　 　　　　　　　　　　　　　　　　　　　　　　　　 　　　　　　　　　　　　式(2.2)　　 恒成立(其中的“：”是张量的双点乘积运算符号)，并且在力的边界上满足力的边界条件： 式(2.3) 这就是虚位移原理。 3、虚应力原理 对于任意的虚应力场和虚位移场，等式 式(2.4) 恒成立，这就是虚应力原理。 式中: E—应变; —虚应力场; 和—是较位移u高一阶的二阶张量; F, S, dA, dV的意义同前。 2.3.3能量方程 能量方程为： 式(2.5) 式中，V为现时构形的体积；为应变率张量；q为体积粘性阻力偏应力张量： 静水压力： 伽辽金法弱