欢迎访问《计算力学学报》编辑部网站！

欢迎光临《计算力学学报》官方网站！

王计真,陈海波,陈红永.基于分区加速和总体共轭梯度法的耦合界面数据传递问题研究[J].计算力学学报,2015,32(2):280~286

码上扫一扫！

基于分区加速和总体共轭梯度法的耦合界面数据传递问题研究

Research of data transfer on coupling interface based on partition acceleration and global CG algorithm

投稿时间：2014-02-26 修订日期：2014-04-10

DOI：10.7511/jslx201502023

中文关键词: 数据传递耦合动力学紧支径向基函数分区加速总体共轭梯度迭代方法

英文关键词:data exchanges coupling dynamics compactly supported radial basis function partition acceleration global conjugate iterative method

基金项目:

作者	单位	E-mail
王计真	中国科学技术大学近代力学系中国科学院材料力学行为与设计重点实验室, 合肥 230027
陈海波	中国科学技术大学近代力学系中国科学院材料力学行为与设计重点实验室, 合肥 230027	hbchen@ustc.edu.cn
陈红永	中国科学技术大学近代力学系中国科学院材料力学行为与设计重点实验室, 合肥 230027

摘要点击次数: 2137

全文下载次数: 1795

中文摘要:

对于耦合动力学问题的分析过程,在界面上需频繁进行数据交换。为此,基于紧支径向基函数和多项式基函数推导了界面数据传递的插值算法,给出了传递矩阵的具体形式。通过分析时间复杂度,找出该算法在大节点量时效率不高的原因在于径向基矩阵的构造和传递矩阵的计算。为加快径向基矩阵的构造速度,提出分区加速处理以提高相关节点的搜索效率;为避免传递矩阵求解过程中的求逆运算,将其转化为多右端项的大型稀疏对称线性方程组问题,引入多右端项的总体共轭梯度迭代方法求解,并讨论了初始估计矩阵的选取方法。数值算例结果表明,结合使用分区加速原理和总体共轭梯度迭代方法,可在不损失插值精度的前提下显著提高求解效率。

英文摘要:

Frequent data exchanges take place on the interface in the solution of coupled problems.For this problem, a kind of interpolation algorithm was derived from the compactly supported radial basis functions and the polynomial basis functions;and the specific form of the transfer matrix was given as well.By analyzing the complexity of the algorithm, it was found that constructing the radial basis matrix and computing the transfer matrix consumed most CPU time, causing computational inefficiency when the number of nodes on the coupling interface was large.The partition acceleration principle was proposed to accelerate the construction of the radial basis matrix by improving the efficiency of searching the related nodes.To avoid matrix inversion, the process of computing the transfer matrix is converted to solve a multi-right-hand sparse symmetric linear system of equations, and the global CG iterative method was introduced to deal with it.Finally how to select the initial guess matrix was discussed.Numerical examples show that a combination of the partition acceleration treatment and the global conjugate gradient iterative method can improve numerial efficiency greatly without accuracy loss.

查看全文查看/发表评论下载PDF阅读器