面向服务的有限元并行计算网格系统设计研究.pdf

第23卷第3期 湖 北 工 业 大 学 学报 2008年06月 of V01．23No．3 JournalHubei of Jun．2008 UniversityTechnology [文章编号]1003--4684(2008)03—0001—03 面向服务的有限元并行计算网格系统设计 宋庭新，李 慧 (湖北工业大学机械工程学院，湖北武汉430068) [摘要]为了提高有限元法计算的速度和共享网络资源，结合目前主流的网格系统体系结构和有限元并行 求解算法，提出了面向服务的有限元并行计算网格系统的总体框架，讨论了有限元并行计算基本理论和基于 Web Services技术的计算网格系统设计方法． [关键词]Web服务；有限元法；并行计算；网格计算 [中圈分类号]TP312 [文献标识码]：A 随着科学研究和工程技术的不断发展，出现了 过程主要有以下几个阶段[4]：1)数据准备阶段，包括 诸如航天飞机、空间站等大型或超大型结构，以及复 结构离散化、单元信息、节点坐标生成以及边界条件 杂非线性力学分析计算等问题，这些都需要借助有 和初始条件的确定；2)单元分析阶段，包括计算单元 限元模型进行大量计算．然而传统的串行计算方法 刚度矩阵、单元载荷向量等，对于动力分析，还需要 因运算速度的限制，已不能满足大型工程的复杂计 计算单元质量矩阵；3)生成总体刚度矩阵和总载荷 算，这导致人们开始致力于有限元并行计算方法的 向量；4)约束处理；5)方程组求解，得到全部节点位 研究‘11． 移；6)计算应变和应力． 网格计算技术的发展，使人们认识到可以充分 并行计算主要是针对2)，3)，5)这3个步骤进 利用网格系统提供的环境进行有限元并行计算．“网 行的．对于2)，串行单元分析是按照单元顺序进行 格”是一种面向虚拟组织的资源共享技术[2]．在网格 的，每分析完成一个单元，就将其叠加到总体刚度矩 计算环境中实现有限元并行计算的实质就是将分布 阵中．当采用精度较高的高阶单元时，逐个单元分析 于网络中的计算机通过局城网或广域网连接起来， 将花费大量的时间．由于各单元刚度矩阵的计算是 将大规模的计算任务分配给各个计算结点，从而大 独立的，所以可以同时计算所有单元的单元刚度矩 大提高有限元的计算速度．网格计算在有限元并行 阵，从而达到完全并行的目的．3)中的总体刚度矩阵 计算中的应用是一项新兴的实践研究，在计算力学 和总载荷量的集成，可以采用EBE策略．对于5)中 领域具有重要的应用价值． 大型方程组的求解，传统的直接解法的分解、消去会 占用较大的存储空间，因此，在并行计算中，一般都 1 有限元并行求解方法 采用迭代解法，例如并行预处理共轭梯度法、并行极 小余量法等Is]． 常用的有限元并行计算方法包括区域分解法和 在程序设计中，主要针对单元刚度矩阵计算和 EBE方法等[3]．区域分解法遵循子结构计算思想，线性方程组求解采用并行算法，计算流程如图1所 对要求解的结构采用某种剖分策略，使得各个计算 示．图1中，方程组并行求解主要是对预处理共轭梯 机所承担的部分刚度／质量矩阵可以进行内部自由