单元刚度矩阵分解法在随机平面应力问题中的应用及效率分析.pdfVIP

单元刚度矩阵分解法在随机平面应力问题中的 应用及效率分析1) 2) 鄂国康 汪明昌 蓝双文 姚伟彬 高冠鹏 （澳门大学土木及环境工程系，澳门特别行政区） 摘要 本文进一步研究和分析了一种全新的随机有限元分析法，称为单元刚度矩阵分解(ESMD)法，并将其 应用到随机平面应力问题的响应矩计算中。为了考察该方法与传统随机有限元分析法（如摄动随机有限 元法）比较后所体现的优势以及更好地评估该方法的性能，针对四节点线性四边形单元的平面应力问题 给出了数值结果，并通过所编写的计算机程序进行分析。数值分析得到的结果表明，采用单元刚度矩阵 分解法能使结构响应矩的计算效率明显提升，同时计算机的物理内存使用量将大大减少。 关键词 随机有限元，单元刚度矩阵分解法，摄动法，响应矩，计算效率 引 言 在工程结构或系统中，潜在地存在着许多不确定因素，例如结构材料和几何特性的不确定性， 施加荷载不确定性等等。因此，用确定性有限元无法有效地分析和处理那些工程中存在的大量随机 性实际问题。于是，随机有限元法，一种由传统有限元与随机分析理论相结合的方法，被提了出来 [1]。近三十年来，该方法无疑已经成为一种最有效的数值方法用来分析带有不确定参数的复杂结构 的响应量[2-6]。众所周知，在用随机有限元法进行结构分析的时候，响应矩是主要考虑的响应量。 而在使用传统随机有限元法的时候，庞大的计算量以及计算机本身物理内存的不足使得计算和求解 变得十分的困难，例如摄动随机有限元法。摄动随机有限元法在早期便被用来估计响应矩，当结构 不确定性较小时，使用随机参数的一阶和二阶矩能有效地计算出结构的响应矩，但是当结构的随机 参数有较大的不确定性时，仅仅考虑随机参数的一阶和二阶矩便不足以精确地计算出结构的响应 矩，这时便需要考虑随机参数的更高阶矩项，然而当有限单元数很大或者考虑随机参数的高阶矩项 时，使用该随机有限元法将使得计算量变得非常大从而导致计算效率很低。因此，具有更高计算效 率的随机有限元法期望被提出并用来改善上述谈及的问题。 本文对鄂国康于 1999 年提出的一种全新随机有限元分析方法[4-6] ，即单元刚度矩阵分解 (ESMD) 法，进一步进行研究和分析并应用到随机平面应力问题中。该随机有限元分析法是基于一 些新的构想，即单元刚度矩阵的分解和一些新的求解过程，完全不同于传统随机有限元法，例如摄 动随机有限元法和 Neumann 展开随机有限元法等。不失一般性，本文主要研究了在平面应力问题 中弹性模量作为一随机域下的响应矩计算。理论分析表明，采用单元刚度矩阵分解法计算量和计算 机物理内存的使用量能大大减少。为了数值评估该方法的性能，本文针对四节点线性四边形单元的 平面应力问题编写了计算机程序进行分析和研究。数值分析得到的结果很好地符合了参考文献[4-6] 中提出的理论期望，即采用单元刚度矩阵分解法能使结构响应矩的计算效率显著提升，同时计算机 的物理内存使用量也大大减少。 1 摄动随机有限元法的响应矩计算 根据随机有限元法，随机刚度矩阵 [K ]可以表示为 1) 澳门大学研究委员会项目资助 2) E-mail: ma66545@umac.mo N e [ ] [ ] [ ] K K K + 0 ∑ i αi (1)