基于随机有限元方法的结构可靠性预计.pdfVIP

磊季．；；叁 的应力分布 ，然后结合应力一强度干涉模型，进行 引言 结构的可靠性预计。 航天工程结构一般要求很高的可靠性 ．在这类结 1 随机场 的离散 构的静强度和断裂可靠性分析中．结构设计上的圆 角、缺口等引起的应力集中对可靠度有显著的影响。 即使在常幅载荷下．由于个体的差异．工程结构 结构可靠性预计广泛采用应力一强度干涉模 的应力响应也是不确定的．其中随机因素包括弹性模 型 工程上进行结构的可靠性预计时 ．往往将结构 量、泊松比、抗拉强度 、缺 口半径、几何尺寸等。 在外载荷作用下的应力响应视为常值 实际上 ．由 统计资料和理论研究表明．弹性模量 E较好 于受微观组织不均匀性等内在分散性和载荷历程、 地服从正态分布 ：泊松 比的离散程度较小 ，认为服 工作环境等外在分散性的影响．实际工程结构的应 从均匀分布 ：材料 的抗拉强度较好地服从正态分 力响应表现出很强的随机性 因此 。考虑这些随机 布 ：结构的几何尺寸的离散程度一般 以公差的形式 因素影响的结构可靠性预计很有必要 给出 ，加工尺寸服从正态分布 ，分布参数采用 近年来．国内外很多研究者相继把随机有限元 “3 ”法则确定 引人到结构可靠性领域中．形成 了随机有限元法模 由于参数的空间变异性 ．如果用单一的随机变 型 如将积分法应用于随机结构的有限元分析而形 量来描述 ．可能会低估实际的可靠度 ．因此用随机 成的积分随机有限元法 ．以及用 K—L级数展开表 场来描述结构的几何尺寸和材料属性等随空间坐标 示输入随机场并用非正交多项式混沌表示输出随机 的随机变化情况 由于有限元中使用离散的基本方 响应 ．从而建立一系列确定性递推方程来对问题进 程 。因此首先要对随机场进行离散。为方便起见 ． 行求解的谱随机有限元法等。以上方法只能响应各 随机场离散与结构离散共用一套网格 ．对结构进行 阶统计量 ．无法得到响应量具体 的概率分布函数 ． 有限元网格划分 的同时也进行了随机场的离散 ．这 因此应用受到很大的局限 样随机场就离散成随机变量 本文基于随机有限元方法．将影响结构应力分 在每个单元里 。随机变量可分解为均值和因变 布的载荷 、材料属性等参数考虑为随机变量 ．采用 异引起的摄动量两部分： 蒙特卡罗随机有限元 (MCSFEM)方法 ．计算结构 (e) ) ) 叼 (，Y)= (1+r／) (1) 厂—] ZH-IL lANG YU KEKAOXlNC-, I8 L— — — — — — L 2011年第 3期 ◆总第 153期 季 罗基 E(叼 )：0 (2) E(77“叼 )； ， ， = l，2，…，n (3) 式中：叼 为随机场 叼 (，y)的均值 ；叼 为第 个单 元随机变量相对于均值 的摄动量 ；11,为单元数 ；6 为 Kronecker函数 2 基于MCSFEM 的结构可靠性预计 结构的可靠度就是 “强度大于应力的概率”。 用 来表示安全余量 ．则可以写出安全余量方程 为 = 尺一S (4) 图l 基于MCSFEM的结构可靠性预计流