基于缩减样本窗口的响应面重构方法.pdfVIP

第28卷第5期 计 算 力 学 学 报 Vo1．28，No．5 2011年 1O月 ChineseJournalofComputationalMechanics October2011 文章编号 ：1007—4708(2011)05-0699—06 基于缩减样本窗 口的响应面重构方法 程 晔 “ ， 周翠英 ， 黄林冲 ， 文建华 。 (1．中山大学 工学院，广州 510275；2．中山大学 岩土工程与信息技术中心 ，广州 510275) 摘 要 ：具有隐式功能函数的复杂结构可靠度分析常采用响应面方法。对于功能函数非线性程度较高的结构可 靠度分析 的响应面迭代重构 问题 ，由于一般响应面方法样本窗 口过大，导致迭代分析精度难以提高。本文提 出了 响应面重构中样本窗口的确定原则，并建立了结构可靠度分析的缩减窗 口序列响应面重构方法。多个数值试验 表明：该方法迭代收敛迅速，计算精度较高，应用过程简单，特别是对于非线性程度较高的隐式可靠度分析具有较 高的效率和精度 。 关键词 ：结构可靠度；响应面方法 ；重构；功能函数；缩减 ；样本窗 口 中图分类号 ：TU311．4 文献标志码 ：A 序列响应面迭代过程 中，通过不断缩减样本窗 口大 1 引 言 小 ，提高在接近实际极 限状态曲面设计验算点附近 大型复杂结构或岩土工程可靠度分析的困难 响应面的拟合精度 ，从而提高响应面方法的可靠度 在于极 限状态曲面的非线性程度较高且显式表达 分析精度 。 难 以获得 ，目前 已有方法主要有 MonteCarlo方 2 一般响应面重构方法 法 、随机有限元方法和响应面方法 RSM(Response SurfaceMethod)[1]。MenteCarlo方法虽具有较 响应面方法最早 由数学家 Box和 Wilson于 好的计算精度 ，但计算工作量太大，目前还难以在 1951年提出，该方法通过一系列确定性 “试验”构 工程中推广应用；随机有 限元方法不 易编程实现 ， 造一个曲面来模拟真实曲面，该曲面即称之为响应 对非线性问题和基本．变量变异系数大的情况 ，误差 面。将响应面方法应用于可靠度分析中，一般先设 仍较难控制 。响应面方法是 目前最具生命 力的复 计一系列变量值 ，每一组变量值构成一个试验点也 杂结构可靠度分析方法之一 ，具有思路简 明，计算 称样本点，利用确定性分析获得每个样本点的功能 量小和精度较高的特点。 函数值 ，通过变量组和功能函数值组构造一个明确 自从 1985年Wong[4等将响应面方法首先运 的函数关系，用 以近似代替难 以直接表达 的真实功 用于边坡可靠度分析 以来 ，围绕响应面构造方法 ， 能函数 ，再利用常用 的可靠度分析方法 ，如一次二 响应面函数设计 、响应面函数验算点解算方法等方 阶矩或二次二 阶矩方法求解。本文以Bucher[5于 面取得了大量研究成果∞ ，然而已有部分研究表 1990年提 出了二次响应面法为例说明其分析过 明，响应面方法迭代计算过程中，当迭代验算点接 程。近似响应面函数一般取如下形式 的不含交叉 近实际极 限状态曲面设计验算点时，由于响应面样 项的二次曲面 本点分布范围较大，易出现收敛速度减慢甚至难以 Y一g(z，z：，… )一口。+∑bz+∑ 收敛 的现象；即使迭代收敛 ，计算验算点附近响应