Excel数据分析在土坡稳定性计算中的应用

://OI10.15923.cnki.cn22-1382t.2018.2.17 Dj

Excel数据分析在土坡稳定性计算中的应用

()安徽工业大学工商学院,安徽马鞍山 243100

周 靖

摘 要:以一均质黏性土坡为例,基于整体圆弧滑动法,采用E得到xcel软件进行数据分析,安全系数值。

关键词:黏性土坡;稳定性分析;整体圆弧滑动法;Excel

()中图分类号:TU432 文献标志码:A 文章编号:1674-1374201802-0200-04

AlicationofExcelbaseddataanalsisppy

insoilsloestabilitalculationpyc

(,A,M)Industrial&CommercialColleenhuiUniversitfTechnoloaanshan243100,Chinagyogy

ZHOUJing

theinteralarcslidinethodtogetthesafetoefficients.ggmyc

:;;;KeordsstickoilsloestabilitnalsisinteralarcslidinethodExcel.yspyayggmyw

:,AbstractTakinomoeneousclaloeasanexamlewealxcelfordataanalsisbasedongahgysppppyEy

0 引 言

1 计算方法

土坡稳定性一般用土坡稳定性安全系数来表

示。相对于无黏性土坡而言,黏性土坡的安全系数计算比较复杂。目前黏性土坡稳定性计算的方法有很多种,常用的方法有整体圆弧滑动法(包、瑞典条分法(括稳定数法)包括总应力法和有效应力法)和折线滑动法等。几种计算方法都是基于整体圆弧滑动法,假设滑动面为圆筒面。对于均质简单土坡可以转化为平面问题,如图1所示。

土坡是指具有倾斜坡面的土体。通常可分为

天然土坡(由于地质作用自然形成的土坡,如山坡、江河岸坡等)和人工土坡(经人工挖、填的土工。由于建筑物边坡,如基坑、渠道、土坝、路堤等)土坡表面倾斜,土体在自重及外荷作用下,将出现

]1

。滑坡将危及到滑坡体及自上而下的滑动趋势[

附近人的生命和财产安全。此外,河岸的滑坡还会造成很大的波浪,使长距离内产生灾难;土坝、河堤的滑坡还会引起垮坝,乃至发生大洪水,其损

]2

。失不堪设想,因此,边坡的稳定计算意义重大[

收稿日期:2018-02-19

:)基金项目安徽省 2018年度高校优秀拔尖人才培育项目(x2018156gyq

,:女,汉族,安徽石台人,安徽工业大学讲师,硕士,主要从事数值模拟方向研究, 作者简介:周 靖(1986-)E-mail307235739@q.q

com.

第2期 周 靖:Excel数据分析在土坡稳定性计算中的应用

201

圆弧法对整个下滑土体进行受力分析,采用Ex-图 C1 均质土坡整体圆弧滑动

点为坡脚,CB为坡面,BA为坡顶。假定

滑动面为圆弧面AC为刚性体(滑弧上的法向反力,圆心为ON,滑通过动土圆心体ABC)

。利用力矩平衡条件进行安全系数K的评价[

3

]:K=

抗滑力矩τfAC⌒

aR滑动力矩=

G(1

)其中,τf为土的抗剪强度,

可运用库伦公式τf=c+σtanφ计算求得,对于土坡而言,式(()A⌒

C2)中σ即为重力的2余弦分量;为滑弧弧长,a为重力偏心距,即重心离滑弧圆心的水平距离。

这种计算方法的难点在于滑动土体的重量和重心的位置很难确定。由于土体为均质土体,因此转化几何问题即下滑土体的面积和弧长很难确定。因此演化出多种计算方法。条分法分析土坡稳定性迄今已近百年,目前常用的计算方法还是

条分法[4]

视为平行四边形。条分是将土体分成若干,重心即为平行四边形的形心土条,将土,条

测量出土条底面中点的法线与竖直线的交角,单独分析每一个土条的受力及力矩的平衡。但由于计算量大、测量结果不精确,导致手动计算难度很大,计算结果的准确性也不高。计算土坡稳定性安全系数的方法通常有两种:一是对构成土体的土条进行受力分析;二是对土坡圆弧滑动体进行

整体稳定性分析[

5]

处。目前有很多。学但这两种方法均存在不足之者运用[[,][][0E]

x等cel6-7

]软件,进Auto-析A,D28提高了计算速度与准确性,MATLAB9,FLAC1行分,但计算方法都是基于条分法,对土条进行受力分析,对土条进行受力分析存在静不定问题。为解决此问题,往往将土条所受的某些应力当零处理。因此,由此法计算的土坡稳定性系数必然存在误差比较大的问

题[8]

。文中采用的是第二类计算方法,直接基于

cel软件建立土坡及滑裂面的几何模型,

然后将土体根据步长划分成近似三角形和平行四边形土条,计算出下滑土体的总面积即重力G矩平衡条件找出重心。最后代入式(1),得到安全再根据力

系数。

2 计算模型

E用Micxrcoesl中大量的公式函数可以应用选择,使oftExcel可以执行计算,分析信息并管理电子表格或网页中的数据信息列表与数据资料图表制作,可以实现许多方便的功能,使用者可获得更多便利。文中即利用Excel强大的数据处理、便捷的公式函数和所支持的面向对象化的程

序设计语言5

]面列举一个计算实例VBA,

对土坡进行稳定性分析[,计算模型如图2所示。

。下图2 计算模型示意图

土坡坡高H为均质黏性土,黏聚力20cm为,1坡角0kPβ为a

,内45摩°

,角土体为20°,重度为φ为,圆心为18kON/点m3,半径,假定滑动面R为40m为。

AC弧,

通过坡脚3 数据分析

根据计算模型通过Excel建立几何模型,

将坡脚C点设为原点,水平方向为x方向,竖直方向为y方向,坡面CB直线的几何方程表示、坡顶BA都可以用折线和出来,坡面tanβ标x·x,坡顶CB为y=

BA为y=20。滑弧AC圆心横坐0为将可调节变量步长11.69487,圆心纵坐标y0为38.25219。

为若干个土条,将两端土条视为三角形Δx设为0.5,将滑动土体划分,中间土条视为平行四边形,通过坐标换算得到相应的面积,

根据力矩的平衡将每一个土条对x轴和y轴取C202

长春工业大学学报 第39卷

滑土体重心横坐标x纵坐标y3.8891,c为2c为9.620118。整体通过图形拟合得到Excel计算模型,如图3所示。

矩,将总的弯矩之和除以面积即得到重心坐标,下

力,不需要利用CAD绘图建立几何模型,Excel可以直接建立几何模型,将土坡及滑裂面表示出来,并为后续计算提供依据。

,与滑裂面圆心的夹角θ为2求得滑动力矩M3°

⌒

在计算出重心坐标之后,可以推导得到重心

力矩Mf=(GcosθtanC·R为φ+c)·A代入式(最终得到安全系数K=248200.1,1)完成了该土坡的稳定性分析。此计算1.883924,

在设定好公式函数之后可以重复利用。

裂面半径R、滑裂面圆心与原点连线后与x轴的

图3 Excel计算模型示意图

为1利用式(得到抗滑=G(xc-x0)31742.2,2)

由上可知,该算例中首先设定设计参数,如滑

夹角α、坡角β、内摩角φ、坡高H、滑裂面圆心横坐标x0,纵坐标y0,计算步长Δ土的重度γx、(。等,见表1该参数都为可调节变量)

利用E 由图3可知,xcel强大的数据处理能

表1 设计参数列表

R40

α73

β45

φ20

H20

x011.69487

y038.25219

Δx0.5

γ18

计算参数 通过力矩平衡公式找出重心坐标,

如土条分割线横坐标x、与土坡交点纵坐标y1、与滑裂面交点纵坐标y2、土条面积S、单位土条

水平方向偏心距e竖直方向偏心距e土条xx、y、方向面积矩Mx、见表2。y方向面积矩My等,

表2 重心计算参数示例

x0

y10

y2012

y2-y1

0

SexeyMxMy0.51.01.52.045.045.546.046.046.547.0…

-0.14931-0.29155-0.42680-0.5551516.0989016.8696517.6813917.6813918.5391819.44933

…

0.51.5…

0.6493061.2915481.9268052.5551483.9010983.1303512.3186132.3186131.4608210.550673

…

0.1623270.4852140.8045881.1204882.1340241.7578621.3622411.3622410.9448580.502873

…

0.250.751.251.7544.7545.2545.7545.7546.2546.75…

0.0876730.2647860.4454120.62951217.8659818.2421418.6377618.6377619.0551419.49713

…

0.0405821.0057351.96085595.4975679.5432762.3225262.3225243.6997023.50933

…0.36391

0.0142320.1284780.3583730.70536038.1264232.0671725.3891225.3891218.004419.804587

…

202020202020

表2数据在建立好公式函数之后可以根据步长下拉至滑裂面端点即坡顶B再根据A的位置,圆弧法进行一系列的换算,最终换算参数如土条

总面积S、土条x方向面积矩总和、土条y方向重心与圆心连线后与y方向的夹角θ、滑动力矩

面积矩总和、划裂土体重心横坐标x纵坐标yc、c、

第2期 周 靖:Excel数据分析在土坡稳定性计算中的应用

203

抗滑力矩Mf等,见表3。M、

表3 换算参数列表

ΣS600.202

ΣMx14338.3

ΣMy5774.01

xc23.8891

cyθ0.402628

M131742.2

Mf248200.1

9.62012

所有的换算参数都是利用Excel提供的公式函数换算得到的,最终计算得到安全系数值。 在整个数据分析过程中可变量是可以调节的土坡稳定性计算软件,只需要利用Microsoft

Excel执行计算。也不需要使用AutoCAD软件导入几何模型,直接利用Excel即可完成几何模的,也就是说在建立好计算模型之后,可以应用于不同物理性质指标的土体,取不同的滑裂面半径,求得相应的安全系数值。在整体圆弧法计算中,滑裂面的位置是假定的,滑裂面有无数个,也就是滑动土体的圆心可以任意指定,找出安全系数的最小值才是所要求的土坡的稳定安全系数Kmin这是手动计算无法完成的,但是通过这个计算模,

型就可以轻松实现,只需要改变圆心坐标跟滑弧半径即可。与此同时,在同一个算例中步长也可以调节,步长越短,滑动土体的重力以及滑裂面的弧长计算越精确,重心坐标值也越精确,计算精度也会随之提高。

结 语

在土坡的稳定性分析中,利用,基于圆弧法可以准确E、x快捷地得到cel强大的

数据处理能力计算结果。该算例不需要利用,直接利用AutoCAD建立几何模型,方便计算模型的反复使用Excel就可以完成几何模型的导入,可以针对不同坡脚、坡高的土坡进行数据分析。再通过,在计算模型确立以Excel强大的公式函数建立计算模型后,可以输入不同的参数值,得到不同的计算结果,这样就可以通过改变物理性质指标计算不同性质的土体。还可以改变滑弧圆心坐标及半径取不同的滑裂面,得到不同的计算结果,最终统计分析出最危险的滑裂面,得出相应的最小安全系数值。通过步长的调节改变计算的精度。

该模型应用起来非常简单,不需要购买昂贵

型的建立。利用Excel强大的数据分析能力,输入参数值可以轻松地得到计算结果。但这个模型也有几个需要改进的地方,本身整体圆弧法分析土坡的稳定性就有很多不足,假定的滑动面与实际滑动面有偏差,需要与实例进行对比分析考证。最小安全系数值需要进行无数次指定不同的滑裂面找出最小值,实际操作起来有一些困难。因此计算模型有待于进一步发展。参考文献:

[1] 刘松玉社,2016.土力学[.

M].4版.北京:中国建筑工业出版[2] 程文华,刘爱国,高建勇行土坡稳定计算[J].山西建筑.利用A,2u0t1o2C(A2D和0):28E4x-c2e8l进[3] 李伟昌,

程紫华5.

稳定性的评价[[4] 陈书申,陈小平J.扩大分条数目的瑞典分条法对土坡.土力学与].山西冶金,地2基01基6础(5[)M:42]-.544.版.武汉:武汉理工大学出版社,[5] 寇海磊2010,29.土坡稳定性分20析15方.法综述[(13):83-83.

J].价值工程,[6] 邓东平,李亮,赵炼恒应用[.Excel在边坡稳定性分析中的[7] 戴自航J],.岩土工程技术,沈蒲生.土坡稳2定00分9(析6)简:3化04-B3i0sh8.值解[J].岩土力学,2002,23(6):760-764.op法的数[8] 李军,王星华,秦曦青析研究[J].岩土力学,2.W00in6d,o27ws下的土坡稳定性分[9] 胡辉.基于Matlab的边坡稳定分析(8):1(3瑞典条分法65-1368.)的解析计算[[10] 陈忠源剪J]胀.城市道桥与防洪,角对土质边坡稳2定017安(9全)系:19数8-的199影.响[J].长春工业大学学报.

,2015,36(2):235-240.4