应用改进算法对称逐步超松弛预处理共轭梯度法进行大体积混凝土仿真计算.pdfVIP

第 38 卷第 3 期 河 海 大 学 学 报 ( 自 然 科 学 版 ) Vol . 38 No . 3 ( ) 2010 年 5 月 Journal of Hohai University Natural Sciences May 2010 DOI :10 . 3876/ j . issn .. 2010 . 03 . 008 应用改进算法的对称逐步超松弛预处理共轭梯度法 进行大体积混凝土仿真计算 韩 　林 ,张子明 ,倪志强 (河海大学力学与材料学院 ,江苏 南京 　210098) 摘要 : 针对大体积混凝土有限元仿真计算中存在计算规模大 、计算时间长的问题 ,为克服使用直接 解法导致占用内存多、耗时长和精度得不到保证的缺点 ,充分考虑大体积混凝土结构有限元仿真分 析中温度场和位移场求解的特点 ,用改进算法的对称逐步超松弛预处理共轭梯度法取代直接解法 , 进行温度场和位移场计算中大型稀疏对称正定线性方程组的求解. 采用算例计算对比了使用 2 种 方法计算的内存 占用 、耗时和计算精度等情况 ,结果表明 ,该方法可完全消除整体劲度矩阵带宽对 存储空间需求的影响 ,减少内存 占用和计算耗时 ,且使得计算精度可灵活控制. 关键词 : 大体积混凝土 ;对称逐步超松弛 ;预处理共轭梯度法 ;有限元方法 ;仿真计算 ( ) 中图分类号 :TU311. 4 　　　文献标志码 :A 　　　文章编号 :1000 1980 2010 030278 06 当前 ,采用有限元方法进行仿真计算是大体积混凝土结构温控分析方法的主流. 大体积混凝土结构的施 工期和运行期历时较长 ,在仿真计算中数据存储量及计算量都很大 ,计算耗时很长. 为解决此问题 ,可以从两 方面考虑 :一是从改进计算模型入手 ,即在理论指导下 ,减少网格和计算时段的划分数 ,在一定程度上减少计 [ 1] [2 ] 算规模 ,从而减少计算时间. 如朱伯芳等提出的并层算法 、弹性徐变体分区异步长算法 等. 但无论计算模 型如何改进 ,对于大体积混凝土结构仿真而言 ,计算规模依然很大. 二是通过改进计算方法 ,针对计算过程中 的主要耗时环节 ,使用高效率算法 ,在不减少问题规模并保证计算精度的情况下 ,减少计算时间. 本文即通过 改进计算方法来提高计算效率. 有限元计算可归结为大型稀疏对称正定线性方程组的求解 ,而此求解过程在整个计算中最为耗时. 通常 ( ) 使用的基于系数矩阵分解的直接解法 以下简称直接法 虽然在理论上可精确解决任何非奇异问题 ,但由于 该算法计算中操作元素过多 ,导致计算量和内存需求较大 ,在大型问题求解时多要进行分块求解 ;再者 ,在该 算法的计算机实现时 ,大量舍入误差的积累也使得计算精度得不到保证. 因此 ,可以采用其他算法代替直接 法 , 以减少计算时间 ,提高计算精度. 针对有限元方法导出的整体劲度矩阵的特点 ,采用某种预处理形式的共 ( ) [3 ] 轭梯度法 conj ugate gradient method ,CGM 进行迭代求解 ,既可充分利用系数矩阵的稀疏性 ,减少内存 占用 , 又可快速收敛至期望精度的解.J acobi 预处理共轭梯度法和多搜索方向预处理共轭梯度法已在混凝土温控