Neumann随机有限元的一种推广形式.pdfVIP

维普资讯 第22卷第6期 计 算 力 学 学 报 Vo1．22，No．6 2005年12月 ChineseJournalofComputationalM echanics December2005 文章编号：1007—4708(2005)06—0681—04 Neumann随机有限元的一种推广形式 杨 杰 ， 陈 虬 (西南交通大学 应用力学与工程系，成都 610031) 摘 要：针对Monte—Carlo随机有限元，研究在多样本环境下如何提高单样本计算效率的问题。证明了关于线性问 题Neumann随机有限元法的算法与线性系统的等刚度迭代法是一致的；在此基础上，将迭代效率更高的算法一预 处理共轭梯度法 1入线性Monte—Carlo随机有限元系统，建立了相应的有限元列式；最后，利用两个算例比较了本 方法与Neumann随机有限元法，结果显示随着随机变量离散度的增大，Neumann法的计算效率急剧下降，而本方 法的计算效率表现稳定；对于离散度不太大的随机变量，Neumann法已经不收敛，而本方法在随机变量离散度很 大的情况下，仍然保持很好的收敛性。 关键词 ：Monte—Carlo随机有限元；Neumann随机有限元；预处理共轭梯度法 中图分类号 ：TB125 文献标识码 ：A 限元法的算法与线性系统 的等刚度迭代法是一致 1 引 言 的。在此基础上，通过对预处理阵的比较认识 ，将预 Monte—Carlo随机有 限元是对工程结构进行 处理共轭梯度法有机地 引入线性 Monte—Carlo随 随机分析的一个重要方法，该方法通过抽样计算获 机有限元系统，在理论上证明了预处理共轭梯度法 得随机结构相关量的数值特征 ，其实施的难点在于 有更好的迭代效率和收敛性 ，因而应用该方法于多 计算耗时。目前 的研究侧重于如何减少样本数 样本环境下的单样本计算，将较Neumann法有更 目L1]，如重要抽样法 、分层抽样法、响应面法以及积 大的优越性 ，本文最后两个算例也说明了这一点。 分法等n]；另一方面，关于在多样本情况下，单样 由于着眼于计算效率的研究，本文在叙述上，将不 本计算策略的研究尚未得到重视。 区别 Monte—Carlo随机有限元和多样本计算 问题 从计算角度 ，Monte—Carlo随机有限元可看作 这两个概念 。 是多样本计算问题 。在多样本情况下，由于各样本 2 Neumann法与等刚度迭代算法 之间在所研究问题的几何尺寸、物理性质 、载荷等 结构性质方面以及在结构有限元网络的划分、单元 2．1 Neumann法介绍 的选取等计算性质方面均有相 同或相近之处，可以 Neumann随机有限元法 的基本思想是 ，首 预计各个样本的求解信息也应有近似之处 ，因而这 先求解某一样本 (称之为特征样本)，获得其刚度阵 就为利用已求解 问题的信息去加速待求样本问题 K。和其逆阵K ；对于其余样本，采用Neumann展 的求解提供 了可能。在这方面，较经典成果是 式获得该样本的刚度阵逆阵K_。，K 的前 项展 Shinozuka和Yamazaki提出的基于Neumann展式 开式为 的Monte—Carlo随机有 限元法 ](简称Neumann K一 一 (K。+ AK)一 ≈ 随机有限元法或Neumann法)。该方法的实质是通 ( 一 P +