板料冲压成形数值模拟中有限元方程求解算法的研究.pdfVIP

第 13 卷 第 4 期 塑性工程学报 Vol 13 　No4 2006 年 8 月 J OU RN AL O F PL A STICI T Y EN GIN EER IN G A u g 　2006 板料冲压成形数值模拟中有限元方程求解算法的研究 (北京航空航天大学 汽车工程系 , 北京 　100083) 　刘艳芳 施法中 徐向阳 摘 　要 : 在板料冲压成形的有限元数值模拟过程中 , 有限元方程的求解是有限元分析中的一个重要步骤 。通过分析 非线性方程组的常用求解方法发现 , 无论采用何种方法 , 最终都回归到线性方程组的求解方法上 。提出采用迭代法 取代 自主开发系统 Sheet Form 中的直接法 , 提高了计算效率 。在提出的广义相邻节点以及节点的相邻关系的两个定 义上 , 分析了整体刚度矩阵的形成原理 , 提出一种新的适合迭代法的二维变带宽存储方法 , 最大程度的节省了内 存 , 而且避免了节点编号的优化难题 。最后 , 通过实例验证了迭代算法在求解大型有限元模型时具有较高的计算 效率 。 关键词 : 有限元方程 ; 非线性方程组 ; 线性方程组 ; 迭代法 ; 存储方法 中图分类号 : T G382 　　　文献标识码 : A 　　文章编号 : 100720 12 (2006) 0400 1505 本文在研究有限元方程的常用非线性求解方法 　引　言 的基础上 , 分析了原系统 Sheet For m[4～5 ] 无法快速 分析大型零件成形的原因。本文采用 PC G 迭代法取 采用迭代法求解大型线性方程组是科学计算中 代原系统中的三角分解法 , 提高计算效率 。在提出 的热点之一 , 已经有很多非定常迭代算法[ 1 ] , 例如 的广义相邻节点和节点的相邻关系的两个定义的基 广义极小残余算法 ( Generalized Minimum Re sidual 础上 , 分析了总体刚度矩阵的生成原理以及其中非 Met ho d , 简称 GM R ES) 、拟最小残量法 ( Qua si 零子矩阵的分布规律 。结合迭代法过程中刚度矩阵 Minimum Re sidual Met ho d , 简称 QM R) 、共轭梯 恒定的特点 , 提出一种新的针对迭代法的二维变带 度法 ( Conj ugat e Gradient Met ho d , 简称 C G) 以及 宽存储方法 。该方法节省了大量内存 , 而且避免了 其它基于共轭梯度法 的变种算法 ( PC G , C GS , 节点编号优化的难题 。最后 , 通过实例验证了改进 ) B IC G , L SQ R 等 。其 中 , 预 处 理 共 轭 梯 度 法 算法在有限元方程求解的有效性 。 ( Precon ditio ned Conj ugat e Gradient Met ho d , 简称 PC G) 是求解对称正定系数矩阵最受欢迎的方法之 1 增量型弹塑性有限元方程的求解 一 , 主要思想就是首先将系数矩阵进行预处理 , 然 方法 后用 C G 法求解 , 从而提高了收敛速度 。PC G 法的 成功关键在于预优矩阵的选择 。目前获得预优矩阵 　　根据虚功原理 , 采用修正拉格朗日法描述的增 的方法主要有 3 种[2 ] : 不完全 Chole skey 分解 、取 量型有限元方程的率形式为 系数矩 阵的对角元素或利用对称逐次超松 弛法