滑坡稳定性评价及治理工程设计软件
- 格式:pdf
- 大小:678.50 KB
- 文档页数:5
© 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net滑坡稳定性评价及治理工程设计软件SlopeDesigner①
李长冬a,胡新丽b
(中国地质大学a.研究生院;b.工程学院武汉430074)
摘 要:主要介绍了SlopeDesigner1.2的设计开发思路和主要功能。该软件是基于Windows面向对象技术开发的滑坡稳定性
评价与治理工程设计的一体化工程软件,采用VisualC++进行程序设计,包括稳定性分析子模块和抗滑桩设计子模块2部分,可
对地震作用、暴雨等不同工况下的滑坡进行稳定性评价和推力计算,并能够对(锚拉)抗滑桩(弹性桩、刚性桩)的弯矩、剪力及侧应
力进行计算,最后根据内力计算结果完成配筋计算。在三峡库区兴山县二里半滑坡的应用研究表明,该软件具有很好的工程实用性。
关键词:稳定性分析;滑坡推力;抗滑桩内力计算;治理工程设计;VisualC++程序设计
中图分类号:P642.22;TP319 文献标识码:A 文章编号:100027849(2006)0120095204
滑坡是自然界中一种十分严重的自然灾害,是
联合国“国际减灾十年”计划的重要研究对象之一,
同时也是斜坡破坏型式中分布最广、危害最严重的
一种地质灾害[1]。在目前滑坡防治工程的实际工作
中,迫切需要利用计算机进行辅助计算与设计。
目前在滑坡稳定性分析及治理工程软件设计领
域,国内有许强等[2]开发的Slope2CAD和殷跃平[3]
开发的GEOHZD2.0以及理正系列软件等。但是
很多软件或各模块功能单一而缺乏整体性,或价格昂
贵,使用不便,这在一定程度上限制了软件的应用与推
广,故笔者拟开发一套Windows环境下的集稳定性分
析与治理工程设计一体化的软件。
1 软件特点及编程思路
SlopeDesigner1.2版是基于Windows视窗开
发的,采用VisualC++编程,主要包括稳定性分析
子模块和抗滑桩设计子模块2部分。该软件主要具
有以下特点。
(1)稳定性分析子模块是针对多层岩土体的滑
体而设计的,不仅可以计算单层岩土体的滑体,而且
还可以计算多层岩土体的滑体。
(2)采用滑坡治理工程中应用最广泛的剩余推
力法进行稳定性计算,不仅可以计算在暴雨、地震作
用等因素下多种工况组合的滑坡稳定性,而且还可
以计算已划分任意数目条块的滑坡的剩余推力、各条
块滑面倾角、滑面长度、质量及稳定性系数,并可方便
地设定安全系数,反算该安全系数下的剩余推力。
(3)界面交互性强,方便快捷,输入数据与输出结果都显示在主界面中,且输入、输出数据均支持文
件操作,既可在网格控件中输入数据,又可由文件
(文本格式、Excel格式)导入数据,计算结果可以以
文本格式、Word格式、Excel格式及网页格式输出。
(4)在稳定性分析子模块中可以绘制滑坡的主
滑剖面图及推力曲线图,在抗滑桩设计子模块中可
以绘制抗滑桩随深度变化的弯矩图和剪力图,所绘
图形均可根据实际情况进行实时缩放,并可导出为
AutoCAD格式或BMP图片格式。
解决滑坡防治实际问题的一般思路为:滑坡地
质模型建立→滑坡稳定性分析与评价→滑坡治理工
程设计。具体流程如图1所示。
图1 SlopeDesigner编程思路
Fig.1 MainprogramsequenceofSlopeDesigner第25卷 第1期2006年 1月 地质科技情报
GeologicalScienceandTechnologyInformation Vol.25 No.1Jan. 2006
①收稿日期:2005202213 编辑:杨 勇基金项目:国家自然科学基金资助项目(40202028);湖北省自然科学基金资助项目(2003ABA018)作者简介:李长冬(1981— ),男,现正攻读岩土工程专业硕士学位,主要从事岩土体稳定性评价与治理方面的研究工作。
© 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net2 编程原理
2.1稳定性分析编程原理
在对滑坡进行稳定性分析和治理设计中,我国
铁路部门最先提出了分块传递系数法(剩余推力
法),该方法具有简单适用的优点[4]。SlopeDesig2
ner1.2版的稳定性分析子模块在采用剩余推力法
进行计算时,视岩土体为刚体,以平面二维方式来处
理,即沿滑面走向取单位宽度的滑体进行计算,不考
虑两侧滑体的摩擦力及其之间的挤压力。在考虑重
力、静(动)水压力和地震力作用的情况下,稳定性系
数K为[5]:
K=Tk/Tx(1)
式中:Tk为滑坡体的抗滑力;Tx为滑坡体的下滑力。
对单个滑块而言,其剩余推力为:
Ei=Ei-1・ψi+Txi-Tki・1Ks(2)
Tki=[Wicosαi-Qisinαi-Ui-
(PWi-1-PWi)sinαi]tanφi+ciLi(3)
Txi=Wisinαi+(PWi-1-PWi)
cosαi+Qicosαi(4)
ψi=cos(αi-1-αi)-sin(αi-1-α)・tanφiKs(5)
式中:Tki为第i条滑块的抗滑力(kN);Txi为第i条
滑块的下滑力(kN);Ei为第i条滑块的剩余推力
(kN);ψi为第i条滑块的传递系数;Wi为第i条滑
块的自重(kN);αi为第i条滑块的滑面与水平面的
夹角(取锐角);Li为第i条滑块的滑面长度(m);
PWi为第i条滑块中动水压力产生的侧水压力
(kN);Ui为第i条滑块中动力水头产生的浮托力
(kN);Qi为第i条滑块所受到的水平地震力(kN);
Ks为安全系数;ci为第i条滑块的内聚力;φi为第i条滑块的内摩擦角。
采用此法进行计算的优点是在计算滑坡稳定性
系数的同时,可以并行求解滑坡各条块的剩余推力。
2.2抗滑桩设计编程原理
SlopeDesigner1.2版的抗滑桩设计子模块可
以计算普通抗滑桩和锚拉桩(包括弹性桩和刚性桩)
的弯矩、剪力及侧应力。
对于弹性桩采用k法进行设计计算时,桩顶承
受水平荷载时的抗滑桩挠曲微分方程为[6]:
EId4xdy4+xkBp=0(6)
式中:Bp为桩的计算宽度(m);E为桩的弹性模量(kN/m2);I为桩的截面惯性矩(m4);k为地基系数
(kN/m3);x为桩身水平位移(m);xkBp为地基作
用于桩上的水平抗力;y为桩的深度(m)。
引入变形系数β,得:
β=4kBp4EI(7)
即kBp=4EIβ4(8)
式(8)可写成:d4xdy4+4βx=0(9)
解式(9),得到滑动面以下桩身任一截面的变位
和内力的计算公式。
对于刚性桩,一般采用按角变位法计算。当桩
身埋入一种地层,滑面以下为同一m值[m为地基
系数的比例系数(kN/m4)]时,在桩底为自由端的条
件下,桩身内力计算如下:对于滑面处岩土地基系数
为A,K=A+my地层中的刚性桩,此时有:Δx=
(y0-y)Δψ,然后分别求解y≤y0和y>y0时2种
情形下抗滑桩各深度点的弯矩、剪力及侧应力[7],其
中,Δx为桩身任一截面的位移(m);Δψ为桩的转角
(弧度);y为桩的深度(m);y0为滑动面到桩旋转中
心的距离(m)。
3 应用工程实例
3.1工程概况
二里半滑坡位于三峡库区兴山县高阳镇香溪河
北岸的二里半村。三峡水库蓄水后该滑坡的前部将
被淹没,在库水位变动以及库岸再造、公路切坡
等不利因素的作用下,滑坡将会因失稳而形成地
质灾害。因此,急需对该滑坡进行稳定性评价和
治理[8]。
3.2稳定性分析
根据二里半滑坡的工程地质条件及滑坡的
特征分析,可采用剩余推力法进行稳定性计算。
考虑水库蓄水后滑坡可能遇到的各类情况(高水
位、暴雨及地震等),拟定出不同工况并进行稳定
性计算。
采用笔者开发的SlopeDesigner1.2版对二
里半滑坡进行稳定性计算。选取比较复杂的工
况,即高水位(175m)时遇到暴雨及地震,其计算
结果如图2所示。在此种工况下的稳定性系数
K=0.87583,其剩余推力Ei为2443.468kN(安
全系数Ks=1.15)。
3.3抗滑桩设计
3.3.1 推力计算 二里半滑坡在采用剩余推力法
计算推力时,考虑了水库蓄水后滑坡可能遇到的各69 地质科技情报 2006年
© 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net图2 水位+暴雨+地震力工况下的计算结果
Fig.2 Calculationresultsunderwaterlevel,rainstormandearthquakecase
类情况,特别是最危险的情况。利用SlopeDesig2
ner1.2版分别对在天然状况、设计工况以及支挡后
的校核工况下的稳定性系数及推力值进行了计算,
并根据各滑块的推力值绘制了各种工况下的推力曲
线(图3)。
3.3.2 悬臂桩设计与计算 抗滑桩的布置应垂直
于滑坡的主滑方向,该方案共设计了3种类型的抗
图3 二里半滑坡不同工况下的推力曲线
Fig.3 PushingforcecurvesofErlibanlandslideindifferentcases滑桩,均为C25的钢筋混凝土人工挖孔桩,为单排
布置,共31根,其中A型桩6根(受荷段14m+锚
固段8m),B型桩16根(受荷段12m+锚固段7m),
C型桩9根(受荷段8m+锚固段6m),抗滑桩截面
均为3.0m×3.5m,间距均为5m。根据图3综合
确定滑坡设桩处的推力设计值为3400kN/m。根
据公式计算得到3种桩型的变形系数β均相等,即β=0.191,且均有β・H2>1.0,故3种桩型的抗滑
桩均应按照弹性桩进行内力计算。
将A、B、C型桩的各参数分别输入到SlopeDe2
signer1.2版的抗滑桩设计子模块中,可以求得各
桩的最大剪力、最大弯矩以及最大侧应力,并据此进
行下一步的配筋计算。以下选取A型桩为例进行
计算。对于A型桩,其运行结果如图4所示,其最
大弯矩Mmax=91.919663(MN・m),最大剪力
Qmax=-20.139041MN,最大侧应力:σmax=
3.269463MPa。
3.4锚固深度复核
取岩层产状核算系数K1=0.5,取岩石裂隙及
风化的折减系数C=0.3,岩石单轴抗压强度R=σcw=30MPa。79第1期 李长冬等:滑坡稳定性评价及治理工程设计软件SlopeDesigner