并行PCG算法及其在水坝应力分析中的应用

来源 :清华大学 | 被引量 : 0次 | 上传用户:jiebaidexue
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
考虑横缝开合变化的拱坝地震非线性响应计算分析在拱坝动力分析中至关重要的作用.其中较多使用的有限元程序:ADAP88(Arch Dam AnalysisProgram---1988),随着横缝条数的增加,计算时间急剧增加.其中方程组求解的计算时间占整个问题求解时间的70﹪以上甚至更高,当所计算的时间步数较多时,方程组求解的时间是决定性的.因此要提高程序的效率应该从改进方程组求解的效率入手.原程序采用的线性方程组解法为直接Cholesky分解后求解,此方法适合于中小规模,当问题规模较大时计算效率较低.为提高计算效率,我们将改用带预处理的共轭梯度法(PCG)对程序中的线性方程组进行求解,并将该PCG法并行实现.选择合适的预处理方法是PCG法的核心.通过对预处理矩阵分解因子的分析我们发现,要达到好的预处理效果,原矩阵某些位置的元素是应当保留的.在满足这个条件的前提下,并结合有限单元划分的信息,以及负载平衡方面的考虑,我们对矩阵进行了分块.预处理矩阵的生成以及预处理方程组的求解都可以在各个分块上分别进行.同样出于负载平衡的考虑,在矩阵向量相乘时我们还对原矩阵进行了重新分组.在采取了以上措施后,使用该PCG法对23856阶的实例进行了并行数值计算试验.在清华大学计算机系高性能计算技术研究所"可扩展高性能集群计算机系统"上分别使用1、3、5、9、17、36个CPU进行计算,并对得到的加速比和效率进行了比较和分析.
其他文献
该文分为六章:第一章为综述,介绍了单种群,两种群,三种群及多种群模型,以及二维多项式系统中的Darboux方法,并引入了该言语要讨论的几个开问题.第二章将Darboux方法推广至高
要不要现场指证权据悉,武汉市将赋予交警现场指证权,即不管司机是否认账,交警可现场认定其按了喇叭,并给予处罚。将依《道路交通安全法》有关法规处以100元罚款。相当多的司
该文首先讨论了直角梯形上的Poincaré不等式,使式中常数具体化;并推广到三维空间中,进一步研究了三维空间中单位正方体及直角梯形台上的Poincaré不等式,同样给出了不等式中
陈世英是武都县汉王镇陈家坡村的党支部书记。他是一个憨厚的人,做事有板有眼,村民们很信任他。所以,尽管当初他不大愿意当村干部,但在乡亲们一致推选下,挑起了引领全村人脱
该文主要研究连续的和离散的2+1维孤子方程的分解及其拟周期解的构造. 文中发展了一个分离技术,由此可将连续的和离散的2+1维孤子方程分解为相容的常微分方程或相容的常微分
该文对von Neumann代数中套子代数的相似与摄动、张量积、局部映射问题、套代数的Jordan结构和Lie结构问题进行系统深入的研究.全文共分五章:第一章研究因子von Neumann代数中
该文应用经济学中供求均衡的分析方法,通过对国际上现行资产证券化融资运行机制、特点及实例的分析,从四个方面对中国现行体制下开展资产证券化业务的可能性、必要性进行了分
该文以武汉钢铁集团公司两个主体厂——冷轧厂、二炼钢厂的实际生产问题作为背景,来研究一类不允许等待的混合流水车间排序问题的数学模型、算法和应用.从实际中归纳出理论,
该文分别得到Runge-Kutta方法、单支方法和线性多步方法关于时滞奇异摄动问题的整体误差估计.此外,研究人员还获得用一般线性方法求解多刚性奇异摄动问题的整体误差估计.数据
随着并行计算机的飞速发展,并行计算已成为数值求解刚性微分方程的十分重要的手段之一,为此迫切需要研究刚性微分方程的高效并行算法.文献中对Runge-Kutta法及块方法的并行计