结构可靠性分析有限元法.docVIP

结构可靠性分析有限元法摘要：本文将结构可靠指标计算的优化迭代法与结构分析常用的确定性有限元法结合起来，在可靠指标迭代计算过程中，所需的参数通过确定性有限元法计算得到，从而为大型复杂结构进行可靠性分析提供了有效和实用的方法。 关键词：可靠性分析；有限元法 中图分类号：TU311.2 文献标识码：B 文章编号：1008-0422（2008）08-0168-01 1概述 可靠性方法用来得到的结构各种性能的概率，结构性能通常是用结构反应，如应变和位移，应力和力，以及累积反应如累积塑性应变等来反映的。对大型复杂结构来说，其结构分析本身就是一个非常复杂的过程，同时在进行可靠性分析时，极限状态函数往往难以显式表达，而一般确定性的有限元方法是分析大型复杂结构的有效方法，但其本身并不具备可靠性分析的功能[1][4]。本文将在确定有限元基础上，将可靠度的分析的优化法与确定有限元法结合起来，用来对大型复杂结构进行可靠性分析。 2可靠指标的计算 可靠度指标的几何意义为在标准正态空间内极限状态曲面到原点的最近距离，因此计算可靠指标实际上是一个优化问题[3][4]，即 这里y=y(x),标准正态空间随机向量，由原始空间的随机向量x变换得到，G(y)=g[x(y)] 为变换后空间中的极限状态函数； 对于式（1）线性搜索循环算法，一般都可写成 其中dk――搜索方向矢量； γ――步长参数； 不同算法的主要区别在于搜索方向dk的选择的不同。对于式（1）优化问题，目前最常用的梯度法，即在标准正态空间内沿极限状态曲面的梯度方向进行搜寻，因此上述方法的关键是求解反应量在标准正态空间内对任一随机变量的梯度，可写成 （3） 其中： u为极限状态方程中的反应量，可以是位移，内力，应力等。 Jx，y为雅克比变换矩阵； 可采用直接微分或有限差分的方法从极限状态函数表达式中直接得到； 在通过确定性有限元计算得到，具体计算方法见下节介绍。 3梯度向量计算的有限元法 由前面的计算方法可知，梯度向量 是计算可靠指标的重要计算量，它的计算可以采用一般通用有限元的计算方法得到[2]。 一般通用有限元方法可写成 对式（6）中由有限元基本理论可知， ，K为结构的刚度矩阵，设 则上式可以写成 （7） 式（7）为对任一随机参数h，反应量对其梯度的一般表达式，式的（7）左边的刚度矩阵，按确定性有限元方法由单元刚度矩阵集成。式（7）右边项两项，分别是外力对任一随机变量h的偏导数和结构抗力在位移确定时对随机变量h的偏导数。 4算例 简单框架结构，其结构形式和荷载状况如图1所示，弹性模量为E=210MPa，截面积A=0.00459m2，截面惯矩I=0.0000579m4，为表1提供了随机变量及分布类型，单元采用梁单元模拟，其失效形式为弯曲失效，计算结果如表2所示，从中可以看出，计算结果与文献值一致。 5结论 本章通过对现有确定性有限元方法进行扩展，并与可靠度分析的优化法结合起来，使一般确定性有限元方法能够进行可靠性分析，这为大型复杂结构的可靠性分析提供了一种可行的方法。 参考文献： [1] Klyoshi Ishll,Makoto Suzukl.Stochastic finite element method for slope stability analysis[J].Structural Safety，4(1987)，111~129. [2] Der Kiureghian A, Ke J-B.The stochastic finite element method in structural reliability[J].Probab Engng Mech 1988;3(2):83~91. [3] Pei-Ling Liu,Armen Der Kirureghian．Optimization Algorithms for Structural Reliability[J].Structural Safety,9(1991) 161~177. [4] 董聪.现代结构系统可靠性理论及其应用[M].科学出版社，2001. 注：“本文中所涉及到的图表、注解、公式等内容请以PDF格式阅读原文。” 1