基于龙格库塔法的弹塑性有限元并行计算-计算机工程与应用.PDFVIP

52 2011 ，47 （27 ） Computer Engineering and Applications 计算机工程与应用 基于龙格库塔法的弹塑性有限元并行计算 付朝江 FU Chaojiang 福建工程学院 土木工程系，福州 350108 Department of Civil Engineering ，Fujian University of Technology ，Fuzhou 350108 ，China FU Chaojiang.Elastic-plastic finite element parallel computing based on Runge-Kutta method.Computer Engineering and Applications ，2011 ，47 （27 ）：52-54. Abstract ：An elastic-plastic parallel finite element method based on domain decomposition is studied under the environment of MPI cluster.An algorithm for integrating stress-strain relation based on the third and the fourth order Runge-Kutta method is presented.This substepping scheme controls the errors in the integration process by adjusting the substep size automatically. A parallel substructure preconditioned conjugate gradient algorithm combined with minimal residual smoothing method is de- veloped.The resulting algorithms are implemented on a parallel environment defined by a cluster of workstation.The results show high speedup and efficiency.The algorithm is efficient for parallel computing. Key words ：elastic-plastic finite element analysis ；preconditioned conjugate gradient algorithm ；Message Passing Interface （MPI ）； parallel computing 摘 要：基于MPI 集群环境对弹塑性区域分解有限元并行计算进行研究。提出了基于三阶和四阶的龙格库塔（Runge-Kutta ）方法 对应力-应变关系进行积分的算法。积分过程中自动调整子步大小来控制积分过程中的误差。研制了采用最小残余平滑法的子 结构预处理共轭梯度并行求解算法。算法在基于工作站机群的并行环境下实现。计算结果表明：该算法具有良好的并行加速比 和效率，是一种有效的并行求解算法。 关键词：弹塑性有限元分析；预处理共轭梯度法；消息传递接口（MPI ）；并行计算 DOI ：10.3778/j.issn. 1002-8331.2011.27.015 文章编号：1002-8331（2011 ）27-0052-03 文献标识码：A 中图分类号：TP31 1 引言 对弹塑性问题，如果在积分点的应力引起塑性屈服和采 网络并行计算具有良好的性价比和可扩展性，以及灵活 用各向同性硬化法则，应力可通过求解系统微分方程得到，即 的体系结构等优点。随着基于消息传递的并行程序设计软件 dσ t Î[0 1] dt D ep (σ)Dε ， （1） 平台的出现，极大地方便了并行编程，同时保证了程序的可移