采用响应面函数分析土坡稳定的可靠度.pdf

·422· 采用响应面函数分析土坡稳定的可靠度 薛松涛h 2，徐荣梅1， 陈镕1， 王远功1 (1．同济大学结构工程与舫灾研究所．上海200092：2．近畿大学理工学部建筑学科．日本大阪577-8802) 摘要：本文提出了用响应面法分析±坡稳定可靠度的方法。利用确定性分析中简化的毕肖普法对土坡进行建模和分析， 在适当的区间内用蒙特卡罗法确定各变量的样本点，经过若干次重复计算，由简化的毕肖普法计算的结果可以确定该 土坡的响应面函数．然后利用这一函数借助予JC法计算士坡稳定可靠度的近似值。用同样的方法可得到第二、第三个 响应面函数，并使可靠度计算结果逐步达到良好的精度要求。这一方法可以减少大蕴的计算工作盏。 关键词：土坡稳定： 可靠度： 响应面法； 蒙特卡罗抽样； 毕肖普法； 1引言 传统的土坡稳定问题分析通常采用定值分析方法，工程中用得较多的是简化的毕肖普法、扬布法、费 伦纽斯法等等。定值分析方法是以按经验确定的安全系数K为度量指标，它忽略了荷载效应(滑动力矩) 及抗力(抗滑力矩)计算中各要素实际存在的随机性质对工程结构安全性的影响，与土坡稳定的实际一】：作 状态不符。 如果在分析中能考虑各种计算参数的随机变异性和计算模型的不确定性对土坡稳定性的影响，把土坡 稳定问题作为一个随机过程看待，对土坡稳定可靠度研究，在理论上和工程实践中都有重要的价值。结构 可靠度计算通常采用近似算法，常用的有一次二阶矩法、JC法、Monte-Carlo法、随机有限元法等。由于 有时需要考虑几何、物理等非线性性质和随机特性，采用上述方法计算结构可靠度，考虑到迭代次数太多 或需要大量的抽样在提高精度或减少计算量方面存在一些实际困难。近十儿年来，国外许多学者在响应面 函数分析结构可靠度方面进行了许多研究[i-4]，该方法对大型复杂结构的可靠度计算比较适宜。国内响 应面法在许多方面已得到不少应用，但尚未见到有关响应面函数方法在土坡稳定可靠度上的研究。本文采 用响应面函数研究土坡稳定的可靠度，并由计算分析得到比较合理的结果。 2可靠度分析方法 2．1随机变量的选取 影响土坡稳定的因素较多，但各种研究表明，在影响土坡稳定的诸要素中，对土坡稳定可靠度分析成 果影响最大的主要是土体的抗剪强度指标凝聚力c和内摩擦角9。而土体的容重y因为变异系数较小，在 可靠度计算时可忽略其变异性的影响。因而一般可选取G缈作为基本随机变量进行研究。 2．2极限方程的建立 结构的极限状态实质上是结构工作状态的一个极值，若超过这一极值，则结构处于不安全、不耐久或 不适应的状态。冈此，在原有确定性分析的基础上，将已经建立的士坡稳定安全系数公式按荷载效应和抗 力的关系，进行适当的变换就可以建立基于各种分析方法的极限状态方程。 由图1所示简化的毕肖普法计算简图[7]，可建立如下的极限状态方程： 式中c。，缈，为第i条士条的有效应力强度指标；彬为第i条土条的臼重：∥，为第i条士条底面中点处的孔 基金项m围客杰出青年辩譬麓金资助埂目(5902fi8∞) 作者简介：静松涛(1％3)．熬．j：海人．工攀博士．教授．媾士生肆姊． 隙水应力；口。为第i条±条底面中点的法线与竖直线的夹角：k为土坡稳定安全系数。 ·423· 图1简化的毕肖普计算简图 2．3响应面法的基本原理 响应面是采用一个合适的梯度函数近似于未知函数，通过数值分析建立结构响应y和结构模型参数 xI，x2，?，x。之间的关系，即Y=f(xl，x2，?，x。)，其中f(xl，x2，?，x。)是未知函数。通过若干次数值 分析，可以获得结构响应y和模型参数xl，x2，?，X。的一组样本。借助于数据拟合，可以确定分析模型中 的待定参数。 本文选取二次多项式作为响应面函数，并按C．G．Bucher[3]的建议，响应面函数取随机变量的平方项， 不含混合项。设结构极限状态函数为： g(x)=g(xl，x2) 取下式为响应面函数 g’(Ⅳ)=a+∑b。x。 (2) (3) 式中的系数a，b。，d．需要由5个样本点值建立方程来确定，由于式中不含混合项，在迭代过程中g’(X)将沿 实响应 图2二个随机变量非线性响应图形 着坐标轴x，逼近结构的极限状态函数。确定响应面函数所取的样本点是有限的，不可能有效地覆盖整个样 本空间。本文采用文献[4]中提出的方法，第一次以均值点∥，为中心点，在区间(1u，一弦。，∥，+fa，)中 用蒙特卡罗法选取样本点，迭代尺度参数厂取1．0或2．O。用g(X)在样本点处的函数值可以确定响应面 函数中的系数。响应面函数得到以后，即可计算出极限状态面上验算点的第一个估计值％。第二步选新 2，Xd：∑Ⅲ + ·424· 的样本中心点xM，点x材在均值点∥，到点xD的连线上，并保持在新的中心点处满足g’(X)=o，即 x肼=∥；+(XD-,I．