弹塑性摩擦接触多极边界元法的规划-迭代型算法(英文)  

Program-iteration pattern algorithm for the elasto-plastic frictional contact FM-BEM

作　　者：于春肖[1] 申光宪[2] 

机构地区：[1]燕山大学理学院,河北秦皇岛066004 [2]燕山大学机械工程学院,河北秦皇岛066004

出　　处：《中国科学技术大学学报》2008年第1期70-76,共7页JUSTC

基　　金：the National Natural Science Foundation of China(No.70671088);the Hebei Provincial Natural Science Foundation of China(No.E2007000381)

摘　　要：提出一种基于多极边界元法(FM-BEM)的规划-迭代型不完全广义极小残值法(简称IGMRES(m)),并建立其收敛性理论.新求解算法采用截断技术,在迭代时仅使用前面计算出的部分向量构造新的递推式以计算后面的向量,矩阵和向量的乘积采用多极展开法(FMM)计算,使得计算量和存储量大为减少.通过数值试验证明,新算法可有效地处理弹塑性摩擦接触迭代的繁杂和费时问题,在确保数值计算精度的前提下,可大大减少迭代次数,显著提高计算效率.A new program-iteration pattern algorithm, incomplete generalized minimal residual method （IGMRES（m）） based on the fast multipole boundary element method （FM-BEM）, was proposed for the solution of highly nonlinear equations and its convergence theory was established. With help of truncation technology, a new recursion formulae with the proposed method using only some of the calculated vectors to compute the following vectors, which could greatly reduce the computation and memory requirement. The fast multipole method （FMM） was used to compute the matrix-vector products. Numerical experiments proved that the new algorithm is highly efficient for computing elasto-plastic frictional contact problems, especially for complicated iteration and time-consuming calculation. And it can greatly reduce the iteration times and improve computational efficiency with ensured numerical accuracy.

关 键 词：FM-BEM IGMRES(m)算法 FMM 弹塑性摩擦接触 截断指标 

分 类 号：O241.6[理学—计算数学]