非线性摄动随机有限元递归方程组分析及其求解.pdfVIP

·438· ‘工程力学增刊 1999年 非线性摄动随机有限元 递归方程组分析及求解 张华赵惠麟胡孔国 (南京末南大学土木t程学院．南京2l∞96) 摘要本文在理论分析基础上，将非线性问题有限元平衡方程中的各矩阵按随机变量TayIor级数展开．从 而得到关于这些随机变量的非线性方程．并利用中心二阶摄动法，推导出非线性问题摄动随机有限元的递 归方程组，运用Ncw-Raphs仰选代法求解位移的二阶摄动量。进一步得出非线性问题与线性问题摄动随机 有限元的递归方程组的统一形式及位移二阶统计特性． 关键词摄动随机有限元，非线性，1’aylor级数，递归方程组，New-RaphSon迭代法 1 引言 实际的结构工程存在各方面的不确定因素(随机波动)，随机有限元正是在考虑了这 一随机波动对结构的影响中发展起来的。80年代以来，随机有限元法已在工程领域中得到 应用，且应用范围日益扩大。摄动随机有限元法用于结构的安全、可靠性分析中，是应用 最为广泛的一种随机有限元方法。非线性问题的摄动随机有限元法比线性问题复杂得多， 但结构力学行为多为非线性，所以研究非线性问题的摄动随机有限元具有重要意义，应用 前景极其广阔。文献【l】讨论的是线性问题的摄动随机有限元递归方程组，文献[2】虽提到1F 线性问题，但没有给出递归方程组和位移的二阶摄动量的求解方法。 本文针对非线性问题。推导非线性问题摄动随机有限元递归方程组，并介绍如何运用 New．Raphson迭代法求解位移的二阶摄动量。从而得出非线性问题与线性问题摄动随机有 限元递归方程组的统一形式。最后给出了位移的均值和二阶统计量。 2非线性问题摄动随机有限元递归方程组及求解 在几何非线性与材料非线性问题中，结构的刚度矩阵是几何变形或结点位移的函数。 因此，结构系统力学行为的控制方程(con仃01equation)都可以写为： 板y)矿；F (1) 其中刚度矩阵K(v)依赖于材料的力学特性和结点的几何特性，如果材料特性或几何特 性有随机扰动，则此扰动可视为空间的随机场或随机过程。这样刚度矩阵K(V)就具有不确 定性。另外，系统承受的荷载F也可能是随机扰动的。由(1)式可以看出，结点位移列阵必 张华1975年8月博士研究生 ‘工程力学)增刊 1999年 ·439- 定也具有随机扰动的性质。 建立随机场模型后，则结构有限元平衡方程(1)中K(v)、v、F的扰动量可以用一个随 机小参数a来表示，即可表示为确定部分和随机部分之和。 足(P)=Ko(矿)(1+口) (2) y=％(1+口) (3) ，=昂(1+口) (4) 式中，a=o对应原系统。‰(v)、vo、Fo分别表示原系统的刚度矩阵、位移列阵和荷 载列阵，即为K(V)、v、F对应的均值矩阵，它们是确定性的量。a是均值为零的随机场， 它反映了K(V)、V、F的随机性。对随机场Ⅸ离散后，q可以化为随机向量。【。当有限元 离散网格确定后，荷载列阵F、位移列阵v和刚度矩阵Kfv)就可以在Ⅱ的均值处Taylor级 数展开，并略去二阶以上的项。刚度矩阵K(v)是结点位移列阵v的函数，而v又是随机 向量q的函数，所以K(v)为复合函数。 月 ． 月 H K(r)=Ko(矿)+∑(K)。口，+{∑∑(K)F口．口， (5) r=％+∑‘a，+{∑∑％叩， (6) t=I 一『．I胆l F=矗+∑Ea．+{∑∑昂叩，