导航：安全管理网>> 安全论文>> 建筑安全论文>>正文

路堑边坡锚固防护参数的响应面优化设计

文档作者：王宇。魏献忠，邵莲芬文档来源： 1．中国地质大学(武汉)工程学院，武汉430074；2．黄淮学院		点击数：	更新时间： 2024年12月27日
下载地址：点击这里	文件大小： 465.52 KB 共6页	文档格式： PDF	下载点数： 1 点（VIP免费）

全屏查看

内容预览 [文件共6页]

本文件共6页，如需编辑使用，请下载。

注：预览效果可能会出现部分文字乱码（如口口口）、内容显示不全等问题，下载是正常的。

	文件大小：465.52 KB 共6页文件格式：PDF 下载点数：1 点（VIP会员免费）

下一篇：边坡主动防护系统加固力简化计算方法研究

上一篇： PCCP输水管线爆管隐患分析及预防措施

文本预览

仅提取页面文字内容，供快速阅读使用。

第2 8卷第7期
2 0 1 1年7月
长 江科 学院 院 报
Journal of Yangtze River Scientific Research Institute
V0l-28 No．7
Ju1． 2 0 l l
文章编号：1001—5485(2011)07—0019—05
路堑边坡锚固防护参数的响应面优化设计
王宇 。魏献忠 ，邵莲芬 ·
(1．中国地质大学(武汉)工程学院，武汉430074；2．黄淮学院，河南驻马店463000)
摘要：边坡锚固防护结构参数的选取直接影响防治工程设计，进而影响到公路边坡的稳定性，科学合理地进行参数
选取意义十分重大。提出了响应面试验设计法，并将其应用于锚固防护结构的参数选取上，研究了锚索支护参数
对稳定系数的影响，建立了稳定系数与各影响因子的定量关系模型，给出了稳定系数的残差分布以及不同操作变
量之间的稳定系数等值线图和响应面三维图，并对试验条件进行了优化，得出了各锚索支护参数的最佳取值方案。
同时，将响应面法优化方案与极限平衡法稳定性计算结果进行了对比分析，对优化设计后边坡的稳定性进行了评
价。结果表明，响应面优化设计方案合理可靠，从而为公路边坡锚索支护参数设计和优化提供了科学依据，为公路
路堑边坡的稳定性研究提供了有用的方法。
关键词：路堑边坡；锚固防护；响应面分析；优化设计；稳定系数
中图分类号：P642 文献标识码：A
1 概述
高速公路边坡的稳定性是确保公路安全运行的
重要条件，开展高速公路边坡的开挖设计、加固措施
以及稳定性维护理论与技术的研究，对于当前快速
发展的高速公路建设、道路的安全运行是必不可少
的，在边坡工程治理设计中，寻找一种更为安全、经
济、可靠的方法尤为重要，因而进行支护结构的最优
化设计显得很有必要。
目前关于边坡稳定性最优化设计研究的方法归
纳起来有单因素分析法¨I5 (灰色关联分析等)和
多因素分析法(正交设计 。 、GIS 、基于正交设
计的RBF人工神经网络¨ “ 等)。很多专家学者
大都把研究重点放在稳定敏感性分析上面，而没有
给出稳定系数与各影响因子间的定量关系模型，更
没有给出参数的最佳取值方案。本文将响应面试验
设计法应用于锚固防护边坡稳定敏感性研究上，则
可以很好地解决这一问题。响应面法可用于确定各
因素及其交互作用对指标(稳定系数)的影响，精确
地表述因素与响应值之间的关系。与以往所研究的
正交设计法不同，它是利用中心组合试验拟合出一
个完整的二次多项式模型，通过少量试验便可获得
可靠数据并估计参数，不仅可以有效地建立回归方
程，还可以建立连续变量曲面模型，对影响反应过程
的因子及其交互作用进行评价，通过偏导数求得最
佳参数，在试验设计与结果表达方面更加优良 。
本文以国家重点公路杭州至兰州线重庆巫山一
奉节段(K22+450段)天然状态下高边坡为例，采
用响应面试验设计，对锚固防护影响因子进行响应
面分析，得出了锚索长度、锚固段长度、锚索间距和
预应力的二次多项式模型，并进行响应面优化，得出
了参数取值的最佳方案。同时对锚固支护后的边坡
进行稳定性评价，证明了响应面优化设计的科学可
靠性，为公路路堑边坡锚固支护设计提供了强有力
的支持，该方法具有广泛的应用前景。
2 响应面法简介
响应面法(Response Surface Methodology，
RSM) H ，是一种试验条件寻优的方法，适宜于解
决非线性数据处理的相关问题。它囊括了试验设
计、建模、检验模型的合适性、寻求最佳组合条件等
众多试验和统计技术；通过对过程的回归拟合和响
应曲面、等高线的绘制，可方便地求出相应于各因素
水平的响应值。在各因素水平的响应值的基础上，
可以找出预测的响应最优值以及相应的试验条件。
该方法考虑了试验随机误差，同时将复杂的未
知的函数关系在小区域内用简单的一次或二次多项
式模型来拟合，计算简便。将试验得出的数据结果
收稿日期：2010—09—13；修回日期：2011—01—06
基金项目：国家“十一五”科技支撑计划项目(2008BAC47B0)；交通部西部交通建设科技项目(200906)
作者简介：王宇(1985一)，男，河北沧州人，硕士研究生，主要从事边坡稳定性分析及边坡与滑坡工程治理设计等方面研究，(电话)
13720184614(电子信箱)good541571889@ 126．com。
长江科学院院报 2011生
进行响应面分析，得到的预测模型一般是个曲面，即
所获得的预测模型是连续的。与正交试验相比有很
大的优势，其在试验条件寻优过程中，可以连续地对
试验的各个水平进行分析，而正交试验只能对一个
个孤立的试验点进行分析。目前，响应面分析的试
验设计有多种，但最常用的是Central Composite De—
sign(简称CCD)中心组合设计和Box-Behnken De．
sign(简称BBD)中心组合设计2种。本文采用Box．
Behnken Design中心组合设计进行响应面试验设
计。
3 边坡工程地质概况
k22+450段边坡位于重庆市巫山县两坪乡溪
沟村境内，起止桩号为K22+410段至K22+500
段，长190．0 m，设计为整体式道路，路面宽度
24．5 m，线路走向约255。，顺山坡展布。边坡出露
地层岩性如下。①第四系残坡积土(Q “ )含碎石
亚黏土：黄褐色，稍湿，硬塑状，不均匀地含5％ ～
40％ 的泥灰岩碎石，粒径2～9 em，个别达15 cm，土
层厚度8—11 m，其中重度 =22 kN／m ，内聚力
C=19．6 kPa，内摩擦角 =30．1。；②三迭系中统巴
东组第3段(T2b，)泥灰岩：浅灰一灰色，薄一中厚层
状，隐晶质结构，节理裂隙较发育，部分充填方解石
脉，岩心破碎。岩层产状179。／28。，与山坡形成逆
向坡，有利于边坡的稳定，厚度约6～18 m。开挖土
层主要为上层含碎石亚黏土、碎石土，开挖等级为Ⅳ
级。钻孔资料揭露无地下水。在勘察区布置了3条
剖面线，本文取Ⅱ一Ⅱ剖面进行均匀试验分析研究，
边坡治理采用预应力锚索防护。
4 锚固防护边坡优化设计
4．1 边坡稳定性计算模型
边坡稳定性分析法主要有极限平衡法和数值分
析法2种。极限平衡法是当前国内外应用最广泛的
边坡稳定性分析方法，其中以条分法最为重要。研
究人员在此方法长期的使用和发展中已经积累了丰
富的经验，分析方法成熟可靠，为边坡的稳定性分析
奠定了基础。为此，本文借助于slide边坡极限平衡
分析软件，用简化的Bishop法，视土条为刚体，求出
各种参数组合下的的边坡稳定系数，对边坡稳定性
影响因素进行均匀分析，计算模型如图l所示。
4．2 响应面试验设计
考虑锚索长度￡、锚固段长度 、锚索间距D和
图1 极限平衡计算模型
Fig．1 M odel of limit equilibrium unnlysis
预应力Ⅳ共4种因素对锚索锚固防护边坡稳定性
的影响，对K22+450段处边坡进行响应面试验设
计。参数取值范围按边坡工程手册及室内、现场试
验参数值来确定，并按照单因素试验结果，将其分为
高、中、低(即1，0，一1)3个等级水平。参数取值范
围和因素水平见表1。
表1 因素取值范围
TabIe 1 Range offactors
因素 锚索
L／
长
m
度 锚 度 锚 距 预
N／
应
k
力
N
蒸垦 = =! = 垒 ：
4．3 试验结果与分析
根据Box—Behnken中心组合设计的试验原理和
单因素试验结果，取稳定系数F 为考察指标(响应
值)，借助Design—expert软件进行4因素3水平共
27个试验点的响应面分析试验(见表2)。
表2 响应面试验的因素和水平
TabJe 2 Factors and levels of RSM
27个试验点可分为析因点(1—24)和零点(25～
27)，其中析因点为自变量，取值在4个因素所构成的
三维顶点，零点为区域的中心点，零点试验重复3次，
用以估计试验误差。采用slide边坡稳定性分析程
序，对杭兰高速公路K22+410至K22+500标段的路
堑边坡进行极限平衡分析，计算结果如表3所示。
从表4不同模型方差分析的均方及检验结果综
合来看，二次方程模型的拟合效果要好于其它模型，
Design．Expert系统在这里推荐了2个模型，即平均
模型和二次方程模型，这里应选择相对高次多项式，
即二次方程式模型。表5对能够拟合数据的各种多
第7期 王 宁等 路堑边坡锚固防护参数的响应面优化设计 21
项式模型的复相关系数及均方差和偏差平方和的结
果进行了比较，从比较结果来看，二次多项式模型为
最优。图2为拟合模型的学生化残差分布情况，从
图2中也可以看出，其残差各点的分布几乎在一直
线上，模型拟合效果较好。通过对二次多项式模型
及模型中的各影响因素的置信度分析，二次多项式
模型拟合试验数据的效果是显著的，而且锚索长度
表3 试验方案及结果
Table 3 Diferent tests and results
锚索长度锚固段长度锚索间距 预应力
A B C D
稳定系数
F
I．287 71
1．255 27
1．350 96
1．184 42
1．215 52
1．275 51
1．199 89
1．198 89
1．252 99
1．198 49
1．327 02
1．352 71
1．229 51
1．235 27
1．255 28
1．199 82
1．255 17
1 337 43
1．255 27
1．255 28
1．241 21
1．275 5l
1．241 lO
1．275 51
1．357 43
1．353 51
—1 O 一1 1．235 27
表4 多种模型方差分析比较
Table 4 Comparison of vadance analysis
0f different models
案豢 平方和 均方 F值 概率 ，
表5 R 综合分析
Table 5 R。model statistics for Box．behnken design
类型 播囊 R 校R正2值预g测2值 葬
2F1 0．012 0．963 1 0．940 0 0．877 1 7．454 ×10—3
二次方程6．527×lO一 O．991 6 O．98l 7 0．95l 4 2|建94
议
5 x
采
1
用
0-
三次方程9．497×10 o．999 9 o．999 6 0．991 4 · 矣 o
学生化残差
图2 稳定系数残差分布
Fig．2 Normal probability plot of residual
for the stability factor
与锚索间距和预应力之间交互作用的影响对稳定系
数的影响最大。
以稳定系数为响应值，经回归拟合后，去掉所有
不显著的交互项，是以增加方程的失拟程度为代价
的，因而可以尝试在此基础上增加交互项，为此，手
动优化响应面分析结果，只保留AC和AD两个交互
项(见表6)，模型的失拟值(Lack of Fit)为3．696×
10 达最小，模型最为理想。此时方程拟合的效果
最好，各因子与响应值的二次回归方程为
F = 0．830 56 +0．076 580 ×L 一2．220 8 × 10一。×
L1—0．023 506 ×D 一6．901 × 10一 × N 一
8．051× 10一 ×￡× D +4．190 × 10一 × L ×Ⅳ 一
1．904 27×10一 × +1．654 17×10一 ×L． +
0．012 775×D。+1．942×10～ × Ⅳ 。
由方差分析结果(表6)可知，回归方程描述了各
因子与响应值之间线性关系显著性，由F值检验来判
定，概率P值越小，则其相应变量的显著性越高。
图3为锚索长度和预应力交互作用对稳定系数
影响的等高线和响应面图：图(a)中，锚固段长度和
锚索间距取中间水平时，锚索长度和预应力的适当
增大，可以提高边坡的稳定性，但是当稳定系数达一
定值时，锚索长度和预应力的增长呈负相关趋势；图
(b)中，其它2个因素为中间水平时，可清楚地看出
随锚索长度和预应力的增大，稳定系数不断增大。
图4为锚索长度和锚索间距交互作用对稳定系
数影响的等高线和响应面图：图(a)中，锚固段长度
和预应力为中间水平时，当锚索长度取一定值时，稳
定系数随锚索间距的增大而减小，当锚索间距一定
时，在一定范围内，稳定系数随锚索长度的变化而变
化不大；图(d)中，可更清楚地看出二者的交互作用
较为明显，二者关系并不总是呈正相关或负相关趋
势，而是在一定范围内，稳定系数取得最大值。
o 0 ● 。0 O O o o o o 0 o o o o o
0 o ● O o o o o o o o o o o 0
o o o o o 0 0 0 O o O o 0 o 。 o
o o o o。o一一O o o 0 o 0 o o o o
跚婀～。 4 m n H ” ∞
长江科学院院报 2011年
表6 响应面分析拟合回归方程的方差分析结果
Table 6 Analysis of variances with regression
equation of RSM
( 啊厘回团
图3 F。=，(A。D)的等高线和响应面图
Fig．3 Response surface and contour of F。=，(A，D)
4．4 试验方案优化与讨论
在试验结果分析及模型拟合的基础上，对试验
参数进一步进行优化，即在获得最佳稳定系数的前
提下，锚索支护参数取值的最优方案。分析表明，经
锚固防护后边坡的稳定系数取值范围为1．134 22—
1．353 51，经优化给出了55种优化方案。表7为在
获得最大稳定系数时的锚索支护参数的优化方案，
从优化结果可以看出，锚索长度L=19．69 m、锚固
段长度L。：7．84 m、锚索间距D：3．06 m和预应力
N=721．43 kN，可以得到稳定系数 最大值为
1．353 95。
稳定系数 5 00
●DesignPo~ts
。
点
I已}
磊4．00
孥
xl=A 锚索长度
x2 c：锚索间距 3
． 50
Actual Factors
B．锚固段长度=7．00
D 预应力 70o0o 3 00
6．00 1700 18 00 19．00 2000
A：锚索长度／m
fa)等高线圈
(b)响应面图
图4 F。=，(A，C)的等高线和响应面图
Fig．4 Response surface and contour of Fl：，IA。C)
表7 响应面优化方案
Table 7 RSM optimization designs
注：F划线是指55种优化方案的最优参数组合结果
借助Design-Expert软件的优化设置，同样可以
得出以稳定系数作为目标值的1．184 22—1．353 51
之间的任一F 相对应的锚索参数取值，这些结果可
以很方便地应用于边坡的加固设计上面。
为了评价响应面优化参数的合理可靠性，建立图
1所示模型，对表7所示的优化方案进行验证，表8给
出了极限平衡法稳定性计算与最优方案的对比结果，
可以看出二者得出的结果基本相同，误差范围很小。
同时也进行了其它组方案的验证，得出相似的结论。
边坡岩土体物理力学参数取值：重度 =22 kN／m。，黏
聚力C=19．6 kPa，内摩擦角p=30．1。。
第7期 王 宇等 路堑边坡锚固防护参数的响应面优化设计
表8 结果对比分析
Table 8 Comparative analysis of the results
5 结 语
(1)响应面法通过少量实验便可获得可靠数据
并估计参数，不仅可以有效地建立回归方程，还可以
建立连续变量曲面模型，对影响反应过程的因子及
其交互作用进行评价，确定最佳水平范围。将其应
用于锚索结构支护参数的选取上面，经优化分析可
以得出合理可靠的方案，可广泛应用于其它工程中，
发展前景广大。
(2)响应面法作为一种优化方法，考虑了试验
随机误差，而传统优化是不考虑试验随机误差的。
同时，响应面法将复杂的未知的函数关系在小区域
内用简单的一次或二次多项式模型来拟合，计算比
较简便。文中是以稳定系数作为单一响应值进行了
响应面试验设计，下一步的工作可以考虑其它因变
量作为响应值，进行响应分析。例如可以考虑工程
造价的经济方面等进行试验设计。
(3)响应面法是试验设计的一个基本方法，该
方法同样可用于边坡稳定可靠性的计算。文中得出
了关于稳定系数的响应面函数，由此可以进行边坡
稳定可靠性方面的模拟分析。
参考文献：
[1] XIE Y．Application of Grey Relational Analysis to the Op—
timal Selection of Landslide Treatment Scheme[C]∥ETP／
IITA World Congress in Applied Computing， Computer
Science and Computer Engineering．Sanya，PRC，Aug．08
— 09，2009．
[2] 郑敏洲，简文彬，吴茂明．滑坡敏感性因子分析及其治理
措施研究[J]．岩土工程界，2004，7(12)：29～31．
(ZHENG Min·zhou，JIAN Wen—bin，WU Mao—ruing．Sen—
sitive Factors of Landslide and the Treatment Measles
[J]．Geol画cal Exploration for Non-ferrous Metals，2004，
7(12)：29—31．(in Chinese))
[3] 刘立鹏，于红杰，陈奇．力学、几何参数对土质边坡稳
定性的影响[J]．岩土工程技术，2008，22(3)：89—92．
(LIU Li—peng，YU Hong-jie，CHEN Qi．Sensibility Anal·
ysis of Soil Slope Stability on Mechanical an d Geometrical
Parameter[J]． Geotechnical Engineering Technique，
2008，22(3)：89—92．(in Chinese))
[4] CHUNG C J F，FABBRI A G．Probabilistic Prediction
Models for Landslide Hazard Mapping[J]．Photogrammet—
ric Engineering＆ Remote Sensing，1999，65(12)：1388
— 1399．
[5] BINAGHI E，LUZI L，MADELLA P，et a1．Slope Insta—
bility Zonation： A Comparison Between Certainty Factor
and Fuzzy Dempster—Shafer Approaches l J 1．Natural Haz—
ards，1998，l7：77—97．
[6] 马显春，王雷，赵法锁．滑坡稳定影响因子敏感性分析
及治理方案探讨[J]．地质力学学报，2008，(3)：382—
388．(MA Xian．chun，WANG Lei，ZHA0 Fa—SUO．Sensi—
tivity Analysis on Factors Influencing Izndslide Stability
and Discussion on Its Strength Design{J 1．Journal of Geo．
1ogical Mechanics，2008，(3)：382—388．(in Chinese))
[7] 葛华，刘汉超．万州草街子双堰塘滑坡稳定性影响因
素敏感性分析[J]．中国地质灾害与防治学报，2003，14
(2)：15—18．(GE Hua，HU Han—chao．Sensitivity Anal·
ysis on Influence Factors for Stability of Caojiezi—Shuangyantang
Landslide in Wanzhou f J 1．The Chinese Journal of
Geological Hazard and Control，2003，14 (2)：15一l8．
(in Chinese))
[8] 毛新生，张玉灯，郑涛．基于正交设计的边坡稳定性影
响因素敏感性分析[J]．灾害与防治工程，2008，(1)：36
— 40．(MA0 Xin—sheng，ZHANG Yu—deng，ZHENG Tao．
Sensitivity Analysis on Efect Factors of Slope Stability
Based on Orthogonal DesignI J}．Disaster and Control En．
gineering，2008，(1)：36—40．(in Chinese))
[9] 胡德勇，李京，陈云浩，等．基于GIS的热带雨林地区
滑坡敏感性分析一马来西亚金马伦高原个案研究[J]．自
然灾害学报，2008，17(6)：148—152．(HU De—yong，LI
Jing，CHEN Yun—hao，et a1．GIS—Based Susceptibilitv A．
nalysis of Landslide in Tropical Rainforest Area：A Case
Study of Cameron Highland，Malaysia[J]．Journal of Nat—
ural Disasters，2008，17(6)：148—152．(in Chinese))
[1O]夏元友，熊海丰．边坡稳定性影响因素敏感性人工神经
网络分析[J]．岩石力学与工程学报，2004，23(16)：2703
— 2707．(XIA Yuan—you，XIONG Hai—feng．Sensibility A—
nalysis of Slope Stability Based on Artificial Neural Net·
work I J I．Chinese Journal of Rock Meehanics and Engi—
neering，2004，23(16)：2703—2707．(in Chinese))
l 11 l ARORA M K。DAS GUPTA A S．GUPTA R P．An Artifi．
cial Neural Network Approach for Landslide Hazard Zonation
in the Bhagirathi(Ganga)Valley，Himalayas[J]．In．
ternational Journal of Remote Sensing，2004，25(3)：559—
572．
[12]汪仁官．试验设计与分析[M]．北京：中国统计出版社，
1998．(WANG Ren—guan．Design and Analysis of Experi．
ments[M]．Beijing：China Statistical Publishing House，
1998．(in Chinese))
l 13 j RAYMOND H M．Response Surface Methodology．Current
Status and Future Directions『J]．Journal of Qualitv Tech．
nology，1999，31(1)：30—44．
1 14 l CARTER W H，WAMPLER G L，STABLEIN D M．Re．
view of the Application of Response Surface Method ology in
the Combination Therapy of Cancer J J]．Cancer Treatment
Repo rts，1983，70：133—140．
(编辑：曾小汉)
(下转第36页)
长江科学院院报 2011垒
Research on础sk Probability M odel and Calculation M ethod
for River Closure System
HE Chang—hai 。LIU Yong—yue
(1．Experimental Research Center on Hydraulic Model，State Key Laboratory of Water Resources and
Hydropower Engineering Science，Wuhan University，Wuhan 430072 China；2．Department of Power Generation，
Heilongjiang Provincial Research Institute of Electric Power Exploration and Design，Harbin 150010 China)
Abstract；To make scientific decision on river closure for water conservancy and hydropower projects，this paper
firstly makes a bIief review on the research history of river closure risk．Based on the actual engineering．the rail—
dom factors of river closure system are identified，and the resulting factors of possible risk events are analyzed．Taking
the hydrological，hydraulic and construction uncertainties into consideration，a new mathematical model is pro—
posed and the calculation of risk probability for river closure system is improved．The mod el uses average velocity，
water depth an d average dumping intensity at the closure gap axis as risk variab les to assess the integrated risk of
the river closure system．Furthermore，by comparing the defects and merits of diferent calculation methods，the
Monte·Carlo method of calculating the risk probability based on complete hydraulic calculation is put forward．At
last，problems including the correlation of risk variables，dynamic and static risk differentiation，sampling frequen—
cy，and sampling eITor are discussed and a specific numerical example is presented．The calculation results show
that the integrated risk of the fiver closure system varies greatly with the change of design ed discharge and is closely
correlated with dumping intensity which indicates the organizing ability of the construction．
Key words：river closure system；risk probability；Monte—Carlo method；construction uncertainty；sampling error
(上接第23页)
Response Surface Optimization Design for Cutting
Slope Anchoring Parameter
WANG Yu ，WEI Xian．zhongz
， SHAO Lian—fen '
(1．School of Engineering，China University of Geosciences，Wuhan 430074，China；
2．Huanghuai University，Zhumadian 463000，China)
Abstract：Scientific and reasonab le parameter selection for anchored protection structure of cutting slope is of great
significance as it can directly affect the engineering design，and thereby afecting the stab ility of highway slope．A
Response Surface Methodology(RSM)design is proposed and is applied to the parameter selection of anchoring
protection structure．The impact of anchoring support parameters on slope stability factor is studied，and a quantita—
tive relation model between the stability factor and influencing factors is established．The residual distribution of
stability factor and the stab ility coefficient contour and response surface 3-D graph of diferent operation variables
are also given． Meanwhile．the test conditions are optimized and the optimal value of each parameter is obtained．
Moreover，the response surface optimization is compared with the computation result of limit equilibrium method，
an d the optimized slope stability is evaluated．It is concluded that the optimal response surface design is reasonable
and reliable．As an applicable method for the study of cutting slope stab ility，it provides a scientific basis for the
design and optimization of cutting slope anchored support param eters．
Key words：cutting slope；anchored protection；RSM analysis；optimized design；stab ility factor

网友评论 more

建筑安全论文最新内容

12-27