基于HDMR响应面的结构可靠度计算方法.pdfVIP

：妻 耋 釜 蔽术 预《左技末》2017年第1期总第12o 基于HDMRI]~应面的结构可靠度计算方法 官 华 张守龙 方继伟 (宁波市交通建设工程试验检测中心有限公司 浙江宁波 315124) 摘 要：介绍了响应面法在可靠度分析中的应用进展，提出了一种基于高维模型表示(HDMR)扩展的迭代响应 面法。在详述采用HDMRI~应面法计算结构可靠度的基本原理后 ，采用Matlab语言对若干数值算例进行计 算 ，同时给出了其它计算方法的结果以对比分析。结果表明：本文方法具有较高的效率和可靠的精度，在 结构可靠度分析中具有一定的实用性。 关键词：结构可靠度分析 HDMR 扩展 迭代 响应面法 Matlab DOI：10．13211 cnki．pstech．2017．01．004 引言 Faravelli131则系统地阐述了以空间变量均值为抽样 对于大型复杂结构的可靠度问题，由于较难 中心的包含交叉项的二次多项式作为响应面函 得到其极限状态方程的准确表达式，通常采用有 数 ，并通过方差分析考虑设计变量变异性的影 限元结合Monte—Carlo(M—c)抽样法、响应面法 响。作者同样采用了2水平因子设计和最小二乘 或随机有限元法进行分析。有限元结合M—c法是 估计来拟合响应面，并对一个轻水反应堆压力容 一 种较为直接的计算方法，但要得到精确解通常 器进行了可靠度分析。然而，当采用2水平因子 要耗费巨大的计算成本。响应面法 (RSM)则是 试验设计时，要得到一个二阶多项式响应面函数 一 种近似分析方法，其思想是通过少量的试验 需要进行2次抽样 (N为随机变量个数)，当N较 (抽样)来拟合一个响应面函数作为真实极限状 大时必然耗费巨大的计算成本，而且以均值为抽 态函数的近似表达式，再基于近似响应面函数进 样中心所得的响应面可能与真实的极限状态表面 行可靠度求解。该方法作为一种近似方法，有着 不相符。Bucher等 J贝0提出一种改进的响应面 较高的效率和可靠的精度，但对于非线性程度较 法，即建议采用不含交叉项的二次多项式作为响 高的极限状态函数，其求解精度还有待商榷。随 应面函数，并选取位于标准正态随机变量空间坐 机有限元法 (SFEM)则是一种更直接的计算方 标轴上的 与 ±厂 (i=1，2，…，N)为抽样点 (其 法，但需要对确定性有限元程序进行改造，而 目 中 与叮分别为随机变量 的均值与标准差， 前要形成一个通用的随机有限元程序尚有一定困 并假定厂：3)，再通过回归分析得到响应面函数 难。因此，在实际工程中，响应面法相较于其它 g(x)，并依此求解初始设计点 。继而用线性插 两种方法有着更高的效率和实用性，这也使得该 值函数得到一个新的抽样中心点 并以此获得新 方法得以广泛的应用于结构可靠度分析中。 的响应面函数，最终根据新的响应面函数求解失 1响应面法的发展 效概率。此方法只需要进行4N+3次抽样，相比 响应面法首次I~tBox等 J作为一种数理统计工 Wong与Faravelli的方法在效率上有明显提高。然 具提出，用于搜寻某种化学过程的优化反应所需 而 ，Rajashekhar等 在研究中发现Bucher与 要的条件。随后该方法逐渐发展成为结构可靠度 Bourgund的方法仅对响应面作一次更新是不够 分析的常用方法之一。Wongt21首次将响应面法应 的，而且在非线性程度较高时令厂为常数3不一定 用