车辆薄板有限元分析中的多因子不完全分解预处理解法.pdfVIP

车辆薄板有限元分析中的多因子不完全分解 预处理解法 姚 松 田红旗 1 中南大学轨道交通安全教育部重点实验室，湖南长沙，410075 dynacn@126.com 摘要：薄板是轨道车辆结构的主要形式，本文基于离散Kirchhoff 假设的DKT 弯曲板单元推 导了四边形弯曲板单元DKQ 的构造过程，并进一步阐述了用于一般薄板问题分析的平板单 元的构造。提出了一种“多因子不完全分解” 的预处理方法，与共轭梯度迭代法结合能够 大大加快薄板问题大型稀疏方程组的收敛速度，经过数值试验，说明该方法是稳定可靠的。 该方法避免了常规不完全分解不适用于薄板这样的 “病态”结构的情况。在此基础上，编 写了一般薄板问题分析的有限元程序，程序对结构刚度矩阵采用压缩存贮的方法，节约了大 量内存空间。本文还对分解算法中的可选参数进行了优化研究。通过一个数值试验，本程序 计算结果与商业有限元软件ANSYS5.7 的结果完全一致。 关键词：薄板结构，DKQ 单元，预处理，不完全分解，共轭梯度法 1 概 述 有限元单元法已经成为结构分析的重要方法，薄板结构是轨道车辆的主要结构形式，因 此薄板结构有限元分析已成为车辆结构分析中的重大课题。早期的弯曲板单元大多基于经典 的薄板理论，在以该理论为基础的板单元的能量泛函中，包含位移的二阶偏导数，要求位移 为 1 C 类连续。这给构造板单元带来了困难，由此研究人员将注意力转向了中厚板单元，大 多采用Mindlin −Re issner 中厚板理论，其能量泛函仅包含位移的一阶导数，只要求位移 是C0 类连续，但是用厚板理论建立的单元仅对中厚板有效，当板逐渐变薄时，单元刚度矩 阵中的剪切项占主导地位，计算出的弯曲变形远小于实际变形；当板非常薄时，求得的位移 趋向于零，从而产生了“剪切闭锁”现象。 基于离散的Kirchhoff 假设，通过挠度和转角分别独立插值，然后在若干个离散点上强 迫挠度与转角满足薄板经典理论中的 Kirchhoff 约束，构造出三角形（DKT ）和四边形 （DKQ ）薄板弯曲单元，其泛函的表达式又回复为经典薄板理论的泛函表达式，又自然解 决了“剪切闭锁现象”问题。多个文献表明 DKT 元与 DKQ 元在求解薄板弯曲问题时都显 示出良好的性能，具有较高的精度。在对实际车辆结构进行有限元分析时，由于结构受力复 杂，在承受板平面内的载荷的同时，也有可能板平面外的载荷，因此在进行分析时所采用的 平板单元是平面应力单元与 DKT 弯曲单元的组合而成。由于三角形平面应力单元为常应变 单元，为了提高分析的精度，在本文中我们讨论由四边形膜单元和 DKQ 单元组合而成的平 板单元。 1 教育部博士点基金（20020533007 ）项目资助 1 采用有限元法求解薄板弯曲问题最终归结为求解一组稀疏对称正定的线性方程组， K ⋅u b ，K 为整体刚度矩阵，u 为待求解的位移向量，b 为载荷向量。有限元求解主要 分为直接求解器和迭代求解两大类[1-2]，直接求解是当前应用最为广泛的求解技术，其存贮 方案多采用一维变带宽，通过对总体刚度矩阵直接进行LLT 或LDLT 分解，然后再回代求解， 采用这类技术经过长期使用比较成熟，但是该方法的劣势在于：分解后L 不再是稀疏矩阵， 在分解的过程中会产生大量的“填入”元，因此对于大型稀疏矩阵的 Cholesky 分解不仅耗 费时间，而且占用内存。在求解大型结构问题时速度比较缓慢，而且所需存贮空间和计算量 随结构规模增大而急剧增加，以致于限制了求解规模。 对于象整车结构分析问题，有限元离散方程组的阶数可以达到几十万阶，采用迭代算法 可以仅仅保存刚度矩阵中的非零元素，由给定初值通过若干迭代步骤获得满足一定精度的近 似解。传统的迭代法包括：Jacobi ，Gauss −seidel ，SOR 等等，这些方法收敛速度过慢