变尺度混沌优化算法在二密滑坡锚固方案优化设计中的-吉林大学学报
( )
第 卷 第 期 吉林大学学报 地球科学版
43 1 Vol.43 No.1
年 月 ( )
2013 01 JournalofJilinUniversit EarthScienceEdition Jan.2013
y
变尺度混沌优化算法在二密滑坡
锚固方案优化设计中的应用
1 1 2 1 1
, , , ,
周福军 陈剑平 栾 海 徐黎明 牛岑岑
,
吉林大学建设工程学院 长春
1. 130026
,
吉林省公路勘测设计院 长春
2. 130021
: , , 。
摘要 滑坡治理是一个复杂的非线性问题 影响因素众多 各个因素之间相互制约 传统滑坡支护治
, 。 ,
理方案存在不同程度的浪费 有必要进行方案细部结构的优化设计 基于退火策略变尺度混沌优化算法
。 , 、
对二密滑坡锚固支护方案的锚索结构进行细部优化设计 针对锚固方案 选取锚索布设道数 锚索水平间
距、 、 、 , ,
锚索钢绞线根数 锚索锚固角 锚固孔径 个因素作为优化设计变量 以工程造价作为优化目标 同时
5
, 。
保证滑坡稳定安全 进行方案优化分析 经过优化后的设计方案造价比规范法设计方案的造价减少
, , 。
33.17% 并且采用优化的支护方案锚固后的滑坡体稳定系数提高 潜在滑动面稳定系数提高 3.31% 经
, 、 、 、 。
过比较分析 混沌优化设计方案更加科学 合理 经济 安全
: ; ; ; ; ;
关键词 滑坡 支护方案 混沌理论 优化 结构设计 护坡
中图分类号: 文献标志码: 文章编号: ( )
P642.22 A 16715888201301019207
- - -
MutativeScaleChaosO timizationA liedinAnchor
p pp
StructureDesinofErmiLandslide
g
1 1 2 1 1
, , , ,
ZhouFuun ChenJian in LuanHai XuLimin NiuCencen
j p g g
, , ,
1.De artmento ConstructionEn ineerin JilinUniversit Chan chun 130026China
p f g g y g
, ,
2.JilinHi hwa Surve andDesi nInstitute Chan chun 130021China
g y y g g
: ,
Abstract The controllin desin oflandslidesisacom lexnonlinear roblem asman factors
g g p p y
affectthedesin rocessandthe influenceeachother.Thereismore orlessawaste oftraditional
g p y
,
landslidesu ortin scheme itisnecessar too timizethedesin ofdifferentsu ortin schemes.
pp g y p g pp g
,
Basedontheannealin im rovedmutative scalechaoso timizationalorithm asu ortin schemefor
g p p g pp g
Ermilandslideissu estedthrouh com uter rorammin technolo .Five factorssuchasanchor
gg g p p g g gy
channel , , ,
number anchorhorizontalsacin anchorsnumbers anchoranchorin an leand anchor
p g g g
anchorin drillin diameterarereardedastheo timizin desinvariables.Theenineerin costand
g g g p g g g g
safet ofsloearethe ofo timizationscheme.Thecalculationresultshowsthatthecostofo timal
oal
y p g p p
su ortin ,
schemewillreduce33.17% thanthedesinbasedonsecifications andsafet factorof
pp g g p y
, ,
risesabout3.31%.Thechaoso timizationdesinschemeismuchmorescientific reasonable
landslide
p g
econom andsafe.
y
: ; ; ; ; ;
Ke words landslides su ortin scheme chaostheor o timization structuredesin sloe
y pp g y p g p
收稿日期:
- -
: ( , ); ( )
基金项目 国家自然科学基金项目 4087217040472136 高等学校博士学科点专项科研基金项目 20090061110054
: ( ), , , 、 , :
作者简介 周福军 1985 男 博士研究生 主要从事岩体工程 地质灾害研究 Email zhoufuun0349 126.com
- j @
: ( ), , , , 、 , : 。
通信作者 陈剑平 1957 男 教授 博士生导师 主要从事岩体工程 地质灾害研究 Email chen w 126.com
- @
jp q
, :
第 期 周福军 等 变尺度混沌优化算法在二密滑坡锚固方案优化设计中的应用
1 193
rotection
p
混沌优化算法的基本出发点就是利用混沌的遍
0 引言
,
历性特点 使混沌算法在计算过程中可以不重复地
[ ]
9
, 。
我国正处在现代化建设的高速发展阶段 在水 经历一定范围内的所有状态 李兵等 采用
、 、 、 , ,
利 铁路 公路 矿山和城镇建设中 滑坡造成的危害 映射模型系统产生混沌变量 并将产生的
Lo istic
g
[ ]
1
、 。 ,
之严重 治理费用之昂贵是众所周知的 混沌变量映射到优化变量的取值区间 然后进行混
, 、 。
滑坡防治是一项复杂的系统工程 滑坡的类型 沌搜索 这种基于载波以及二次载波的混沌优化算
、 、 、 、 , 、 ,
性质 规模 机制 发育阶段 稳定状态各不相同且复 法 思路很清晰 易于编程 比较适合连续变量的优
。 。 : ,
杂多变 针对具体滑坡所采用的防治措施适用性不 化问题 但是也存在明显的缺点 首先 搜索遍历性
, , 。 , ,
同 防治要求也不同 更使得滑坡防治复杂多变 效果不好 搜索过程中轨道点分布不均匀 区间两端
, ; ,
在滑坡防治中为了避免出现多次治理 以及防 点数较多 其次 由于 映射产生的混沌变
Lo istic
g
, , , ; ,
止古滑坡复活 个别设计人员的设计过于保守 盲目 量 导致变量只是在单侧领域内搜索 再次 当搜索
。 ,
增加设计安全储备 这一做法的弊端是治理过程中 起始点选择不合适或者优化变量区间很大时 搜索
[ ]
11
。 , 结果很难达到或是接近最优解 。
增大投资成本 从而 滑坡治理方案的优化设计就
, , 退火策略的变尺度混沌优化算法步骤
显得非常必要 即在保证滑坡稳定的前提下 以造价 1.2
。 利用退火策略的变尺度混沌优化方法进行目标
最低的方案作为最优设计方案 传统滑坡设计方案
,
优化只是凭借专家经验对有限几个可行性方案进行 函数的最优化计算 搜索过程中在变量区间内不断
, 。
。 缩小搜索空间 并改变第二阶段搜索的调节参数
优劣比较 近年来研究发展出很多新的滑坡治理定
[ ]
2
, ,
: 、 , 基于退火的思想 随着搜索循环次数的增加 对调节
量优化方法 复合型法 单纯型法 拉格朗日松驰
[ ] [ ] [ ]
3 4 56
-
、 、 、 参数进行修正。具体的算法实现步骤如下:
变量 复形耦合法 模拟退火法 遗传算法
-
[ ] *
7 ) 。 ,
、 。 设定初值 为混沌变量初值 记录优化
灰色模型 神经网络方法 等 1 x f
i
*
, ,
滑坡治理方案优化问题是高度复杂的非线性问 目标值 记录优化目标解 基于 映射产
Lo istic
p g
i
k
生混沌变量 , :
, : 并线性映射到优化变量区间内
题 传统的优化方法普遍存在着以下问题 对输入参 x
i
数敏感, , 。 k ( ); ()
迭代收敛速度慢 容易陷入局部极小值等 x x 1 x 1
= -
i μ i i
m m m k
: ; ; ( ) 。 ()
混沌算法具有独特的性质 随机性 遍历性 p = a + b -a x 2
① ② ③ i i i i i
[ ]
89
-
: ; ; ;
规律性 。混沌算法可以在一定范围内按其自身 式中 为优化搜索次数 为迭代次数 为
=4 m k p
μ i
m
, , 变量区间内的混沌优化变量; 为变量区间的下
的规律不重复地经历所有状态 是一种全局性算法 a
i
m
; 。
可以作为搜索过程中避免陷入局部极小值的一种优 界 为变量区间的上界 用优化变量进行优化搜
b
i
。 : ( ) * , * , * ( ); ,
化机制 笔者采用混沌优化算法对滑坡治理的锚固 索 若f p <f 则p =p f =f p 否则 混
i i i i
k+1 k k
, , ( ), 。
措施进行优化计算 以工程造价为优化目标 并与传 沌变量继续迭代x =4x 1-x 返回进行计算
i i i
, 。 ) ,
统方案进行对比 分析优化方案的科学性 2 迭代计算 k+1 次后 当前的目标函数为
*k+1 , *k+1 。 ),
f 目标解为p 多次进行步骤 1 直至满
1 混沌优化算法基本原理与数学模型 i
*k+1 *k * *
, ,
足 - 时 目标值 为最优值 目
| |<ε
f f f f
o
* *
1.1 混沌优化算法基本原理 标解 为最优解 。 为优化可以接受的极小
p p ε
i io
, , 。
从数学角度讲 工程优化设计就是建立一个以 值 本文取 =0.001
ε
、 ) ,
工程造价为目标函数 以工程设计的安全可靠性能 3缩小各个变量的搜索范围 生成新的混沌变
。 :
指标等为约束条件的数学模型 工程优化设计过程 量
m 1 m m
, + * * ( );
可以转化为在满足所有约束条件下 寻求使目标函 a = p - b -a
i i β i i i
[ ]
10 m 1 m m
( ) 。 + * * ( )。 ()
数 工程造价 值最小的设计方案 b = p + b -a 3
i i β i i i
m 1 m
( , ,…, ), : + / 。
设目标函数为f x x x 最优化问题就 进一步迭代 m
1 2 n
β = β
* *
是在一定的等式约束条件和不等式约束条件下, : ( , )。
求 其中 00.5 将目标解 p 还原到混
β ∈ i
( ) 。 :
解目标函数f x 的极小值 沌区间变量内
( )
吉林大学学报 地球科学版 第 卷
194 43
m 1
* +
p -a
* i i 。 ()
x = 4
i m 1 m 1 二密滑坡稳定性分析及锚固优化
+ + 2
b -a
i i
从而新的混沌变量表达为 2.1 滑坡概况及稳定性评价
( ) * k, 。 ()
1 x x 0.5 5
y = -α +α α=
i i i 二密滑坡位于鹤大高速公路通化路段二密镇一
)
4 利用新生成的混沌变量在最优值附近进行第
, , , ,
带 发育于二密河右岸 坡体长 283 m 宽 86 m 滑
, )、 ) 。 2
二阶段的细搜索 返回第 1 2 步骤进行迭代计算
, ;
坡平面形态呈不规则形状 面积为 10808 m 坡体
, ,
当满足精度要求时 结束搜索迭代 保存当前的函数 3
, , , 。
后部厚 前部浅 平均厚度 12m 体积 126096m
* *
最优值f 和最优解p 。
o io 年 月中旬在路段施工开挖过程中部分
2010 5
) , ,m+1 m / ,
5减小 值 m =m +1 = m 重复第
α α α
,
路段边坡发生塌方 并在山体上部发现较多的张开
)、) 。 ,
3 4 步骤 直到 小于某一值 输出优化设计变量
α
。 ,
裂隙 随着时间的推移 滑坡所在区域地表裂缝数
和优化函数值。
, 、 ,
量增多 裂缝规模增大变长 变宽 滑坡区地表裂缝
) ,
重复以上所有步骤 次 找出 次目标函
6 M M
, 。 ,
最长可达50m 最宽可达20 30cm 同时 为了监
~
* *
,
数的最优函数值 和最优化变量 进行对比 从
f p
o io
,
测滑坡区地表位移情况 在滑坡区内布置 个位移
20
中找到最优的结果。
, ,
监测点 为期一个月的监测发现 平行于主滑线方向
,
通过以上思路 笔者用VC6.0编制了退火策略
,
产生的水平位移最大 累积水平位移最大值可达
, 。
的变尺度混沌优化程序 程序流程见图 1
。 ,
29.8cm 滑坡继续发展将造成严重的破坏后果 严
重影响高速公路的正常施工和人们的安全。
、
滑坡区第四系松散堆积物为崩积 坡积和塌滑
, , 。
堆积物 岩性为碎石混合土 厚度 1.00 8.18 m
~
( )
滑坡区基岩为上白垩统三棵榆组 K sh 一套火山
2
, ,
沉积岩系 由于火山喷发和沉积环境不同 岩相变化
, ,
复杂 滑体发育较多的软弱夹层 成分主要为粉泥质
。
物质 滑坡区内主要岩性及其物理力学性质见表
。
1
表 滑坡岩体物理力学参数
1
Table1 Phsicalandmechanical arametersofrockmasses
y p
3 / /()
岩土名称 重度/( / ) C kPa °
kN m ψ
碎石混合土 19.0 5.00 35.00
灰绿色角砾凝灰岩 24.0 7.80 41.26
紫灰色凝灰质粉砂岩 24.0 3.29 42.51
软弱夹层 19.4 8.00 11.00
: ; 。
注 为内聚力 为内摩擦角
C
ψ
,
滑坡区内典型滑坡纵剖面为 剖面 岩性特
22
-
。 ,
征详见图 从图 可知 滑坡体整体岩层软硬相
2 2
间, ,
两层主要软弱夹层控制滑坡的发生 滑动面在此
。 , 。
层产生 粉砂岩岩层倾向坡外 滑坡为顺层滑坡
,
图 混沌优化算法程序流程图 针对滑坡区内 剖面 采用极限平衡法进行
22
1 -
。
Fi.1 Flow chartofthe roram com utin ofchaos 稳定性计算 从图 剖面岩性分布特征确定 剖
g p g
p g 2 22
-
otimizationalorithm
p g ,
面主滑动面发育在下部软弱夹层底 且主滑面非圆
, 。
弧形 故笔者所研究的滑坡属于非圆弧形滑坡 根
据软层发育情况,
上部软层为潜在滑动面可能发育
的位置。采用简布法( )、 法分别进行稳
Janbu Sarma
, :
第 期 周福军 等 变尺度混沌优化算法在二密滑坡锚固方案优化设计中的应用
1 195
图 滑坡 剖面图
2 22
-
Fi.2 Landslide22sectiondrawin
g - g
, 。 ,
定性计算分析 结果见表 可以看出 滑坡体稳定 / 。
2 kN m
, ,
性系数小于 并监测显示滑坡已处于滑动状态 需 ) 。 ,
1 2 单根锚固力计算 参考规范要求 并考虑防
, 。
及时采取治理措施 防止产生更加严重的后果 止产生群锚效应, ,
锚索水平间距取 6m 锚索道数为
表 滑坡 剖面稳定性计算结果
2 22
- , 。 ( )
道 锚索垂直间距 采用公式 计算锚索单
4 4m 7
Table2 Stabilit resultsof22section
y -
根锚固力:
滑面名称 稳定系数 稳定系数
Janbu Sarma / 。 ()
T = TJ N 7
l anchor
主滑动面 0.811 0.835
: ( );
式中 为单根锚索锚固力 为锚索水平间
T kN J
l
潜在滑动面 0.826 0.884
( ); 。 () ,
距 为锚索道数 根据公式 计算 可得
m N 7
anchor
单根锚索锚固力为 。
1034kN
) 。
2.2 滑坡锚固治理及混沌优化 3锚固段长度计算 防止锚索体从胶结体中拔
基于规范设计滑坡锚固治理方案 , :
2.2.1 出 按下式计算
, [ ] / 。 ()
根据工程治理要求 查阅文献 12 采用不平衡 L KT n dC 8
m1 = l π 1
,
推力法进行滑坡剩余下滑力计算 计算选取软弱夹 , :
防止胶结体与锚索体一起沿孔壁滑移 按下式计算
,
层底层面为滑坡主滑动面 滑坡治理要求达到的安 / 。 ()
L KT DC 9
m2 = l π 2
全系数为 。滑坡 剖面剩余下滑力为 : , ;
1.3 2-2 597 式中 为安全系数 取值 为钢绞线根
K 2.0 4.0 n
~
/ 。 , ,
kN m 滑动面埋深较深大 且倾角很小 采用锚索 ; ( );
数 为钢绞线直径 为锚固钻孔孔径
d mm D
,
支护方案对滑坡进行防治 锚索选取国标 丝钢绞 ( ); ( );
7 mm C 为砂浆与钢绞线允许的黏结强度 MPa
1
, [ ]。
线 钢绞线参数见文献 12 以下为锚索设计主要 ( ),
C 为砂浆与岩石的胶结系数 MPa 为砂浆强度的
2
步骤: / 除以安全系数 。
110 1.75 3.0
~
) 。 ,
1锚固力计算 锚索计算过程中 不考虑锚索 ) 。
4 锚索钢筋截面面积要求
。 ( )
垂直滑动面产生的法向阻滑力 采用公式 计算
6 rN
0 a 。 ( )
A 10
s ≥
所需锚固力: f
2 y
ξ
2
/ 。 () 式中: 为钢绞线截面面积( ); 为锚索轴向拉
A m N
T P cos 6 s a
= θ
( ); ,
: ( / ); 力设计值 锚筋抗拉工作条件系数 本工程
式中 为设计锚固力 为滑坡剩余下滑 kN s
T kN m P ξ
力( / ); 为锚索倾角()。计算得锚固力为 为永久性锚索取 ; 为边坡工程重要性系数,
kN m ° 689 0.69r
θ 0
( )
吉林大学学报 地球科学版 第 卷
196
文档评论(0)