一、资料获取
1. 地形地貌
该坡体地势东高西低，由滑坡后缘山顶高程105m ，前缘高程50m ，相对高差达
55m 。
滑坡体上部与地形较陡的危岩土接壤坡度为4045:度，以下部分地形坡度为2225:度，两侧为南北向的小冲沟，滑坡滑移的长轴方向为340度，第四系土体西薄
东厚，西部岗顶层厚12m :，沟谷34m :；东部土层厚10m 左右，最厚处可达16m 。
此滑坡体南北纵长110m 左右，东西纵长180m 左右，总面积约为2
1.98104m ⨯，总体积约为3
9.1104m ⨯。
2. 岩土类型及其工程性质
该滑坡属上为松散堆积物下为基岩的土质滑坡。
上覆土体颜色自西向东由棕黄变为黑色，根据钻孔以及探槽等勘探资料表明，二中滑坡体的土体结构可以大致划分为七层，分别为煤矿渣堆积层，含角砾粉质粘土层，含碎石粘土层，含碎石粉质粘土层，含碎石红粘土层，粘土夹碎石层，如图1。
图1 二中滑坡体土体结构
3. 水文地质条件
在二中南侧，标高100m 以上，一般为透水不含水层；在标高100m 左右为地形急缓过渡地带，一般为弱含水层。
二中滑坡前缘基岩的稳定地下水位在标高4647m :左右。
含碎石粉质粘土为弱含水层，其下不稳定地下水位标高在5559m :之间。
地下水的渗流压力和润滑作用是该土质滑坡的主要形成原因之一。
二、滑坡稳定性分析计算
1．应用软件
计算机技术在岩土工程计算中已经得到了普遍的应用，虽然商业的专业软件层出不穷，但由于这些软件的不可见性、不可修改性，给使用者带来了诸多的不便。
Microsoft Excel 是美国微软公司研制的电子表格软件，它具有强大的计算、制图、制表和数据库操作功能，使用较为简单、方便，它在各个领域均得到了广泛的应用。
文中以滑坡剩余推力法计算为例，应用 Excel 对黄石市板岩山二中滑坡稳定性进行了分析计算。
2．剩余推力法的基本理论
剩余推力法假定条块间作用力的方向，其重要前提就是假设当前条块在分界面处对下一块体的推力的方向平行于当前条块的底滑面，然后根据平行于底滑面和垂直于底滑面两个方向的合力为零以及最前缘一块的剩余推力为零进行迭代求解。
将整个滑体进行求解后，取其中的第i 号条块，假定第1i -号条块传来的力的方向平行于第1i -号条块的底滑面，而第i 号条块传给第1i +号条块的力的方向平行于第1i +号条块的底滑面(受力图见图2、3)。
图2 折线形滑动面的滑坡计算坡面 图3 第i 号条块受力图 2.1考虑因素如下所示： a ．地震力i P
a
i=g
i c i P w K w •=•
式中：a →为地震加速度， 2m s ；
g →为重力加速度， 2m s ；
i w →为第i 条块单宽自重， KN m ； c K →水平地震系数
b ．地下水压力，若第i 条块上、下两侧面的水柱高度分别为1i h 、2i h (单位为m )，则两侧作用面上的水压力为:
21112i w i u h γ= 2221
2
i w i u h γ=
滑面上静水压力产生的扬压力为：
()3121
2
i w i i i u h h l γ=+
式中：i l →为第i 条块滑面长度， m ；
w γ→为水的容重， 3KN m ；
据此，第i 条块在自重力、地震力和静水压力作用下的平衡方程为：
11i i i i i E T S E --=-+ψ
12sin cos ()cos i i i i i i i i T w P u u ααα=++-
(){}
321cos sin sin i i i i i i i i i i i i S w P u u u tg c l k αααϕ=----+⎡⎤⎣⎦
111cos()sin()i i i i i tg k ααααϕ---ψ=---
i E ，1i E -→为第i 和1i -条块的剩余下滑力，2KN m ；
i α，1i α-→为第i 和1i -条块底部的倾角； i T 、i S 、1i ϕ-→为传递函数；
K →为滑动面安全系数
2.2 滑坡稳定性计算系数
该滑坡所在地区的地震烈度为6度，故计算方案采用6度地震的水平地震系数0.05c K =考虑。
2.3 计算
建立计算工作表录入设定的参数。
建立如表一所示格式的工作表，工作表中需要输入的项目有滑块编号i 、地震加速度a 、重力加速度g 、条块单宽自重i w 、水容重
w γ、水柱高1i h 、2i h 、条块长度i l 、条块底部倾角i α、内聚力i c 、滑动面摩擦角ϕ等。
参数工作表
利用Excel中的公式编辑功能对各计算项目录入公式。
Excel中提供了方便的编辑公式功能,可以根据计算项目的计算式录入公式。
各计算公式用的数学表达式和Excel 表达式如下表所示。
计算公式表
利用Excel的复制功能对各条块的计算项目进行求解。
通过利用Excel提供的句柄复制功能可以对各计算项目进行快速的求解。
当鼠标指向选定的单元格出现句柄标志时,向下拖动该标志完成对该列进行公式的复制,则出现每一行对应的该项计算结果,对以上各列进行复制,则出现全部结果。
计算结果表
对于任一计算剖面，首先假定一个安全系数K 值代入上式中，通过不断调整K 值的大小，直到使最后一条块的下滑力0n E =，此时的K 值即为所求该剖面抗滑稳定安
全系数(当一滑块剩余下滑力小于0Kpa 时，考虑到岩石的不抗拉性质，传递推力去取
0Kpa ，其函数表达式为=IF(B15<0,C16-C17,C16-C17+B15*B18))。
K 值确定后，每一
条块的剩余下滑力，可同样用上述公式计算。
条块剩余下滑力变化曲线
当代入的 1.362209K =时，最后一条块的剩余下滑力0.00011i E =-，且此时，所有条块的剩余下滑力也求得。
三、总结
1．在取岩体处于天然状态下，其力学参数：水容重w r 为1.83g cm ，内摩擦角ϕ为18o ，内聚力i c 为35Kpa ，地震加速度α为20.5m s ；重力加速度g 为29.8m s ，
最终计算此滑坡安全性系数为1.362209。
由计算结果可知在设定的参数下，此滑坡相对稳定。
2．此计算方法的计算过程是透明的,其原始数据、中间结果及最终结果都显示在一张表格中,一目了然。
3．Excel 可以完成重复性计算,对需要反复试算的问题或同类工程计算,只需改变一下有关参数,结果即可产生。
4．利用Excel 的模板功能,将平时工作中计算的工作簿制成模板,以利于以后反复使用。
利用Excel 进行岩土工程计算可以大大的提高工作效率。
