EXCEL在岩土工程计算中的应用

EXCEL 在岩土工程计算中的应用摘要 Microsoft Excel 是Office 中的一个电子表格软件，它在自动化办公中得到了广泛应用，其强大的计算功能和丰富的函数，可以方便地应用于工程计算中去。

本文以滑坡剩余推力法计算和地基沉降量计算为例，介绍Excel 在岩土工程计算中的的一些应用。

关键词 岩土工程 EXCEL 剩余推力法 地基沉降量1. 前言计算机技术在岩土工程计算中已经得到了普遍的应用，虽然商业的专业软件层出不穷，但由于这些软件的不可见性、不可修改性，这给使用者带来了诸多的不便。

Microsoft Excel 是美国微软公司研制的电子表格软件，它具有强大的计算、制图、制表和数据库操作功能，使用较为简单、方便，它在各个领域得到了广泛的应用［1］［2］。

笔者就EXECL 以岩土工程勘察中的滑坡剩余推力法计算和地基沉降量计算为例，介绍一些常见计算问题的解决方法和实践体会，以达到抛砖引玉的目的。

２．EXCEL 的一些主要特点［3］（1）具有强大的函数计算功能：它的内部函数包括对数函数、三角函数、工程函数、字符串函数及逻辑函数等等，它支持公式的编辑、复制、粘贴；同时还支持Visual Basic 编程，通过宏和Visual Basic 可以定义用户自定义函数。

（2）具有强大的数据库功能：可以对数据进行修改、插入、删除、查询、替换、排序、筛选、链接等操作。

（3）计算结果自动更新：更改原始数据后，计算结果自动更新。

３．用EXECL 计算滑坡剩余推力 在《岩土工程勘察规范》（GB50021-94）中用剩余推力法计算滑坡稳定系数Fs 和各滑块的剩余下滑力Pi 计算公式如下［4］：∑∏∑∏-=---=--+⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫⎝⎛=111111n i n n i j j i n i n n i j j i S T T R R F ψψP i = P i-1·ψ+ Ｆst ·T i - R i()()ii i i i i i i i i i n i i i n i j j i i i i i i j Q T Q N L C N R T θθϕψψψψψϕθθθθψsin cos tan .......tan sin cos 1211111==+==---=-++-=+++∏下面以山西大运高速公路清徐县附近一滑坡的稳定性计算为例，介绍利用EXCEL 进行计算的步骤。

根据野外钻探揭露，滑体的主滑断面如图1所示，断面几何参数及滑面的力学参数见图2。

（1） 建立计算工作表录入原始数据建立如图2、图3格式的工作表，工作表中需要输入的项目有分块编号，分块边界的横坐标X ，顶部纵坐标Y 1，底部纵坐标Y 2，滑体容重γi 、滑动面粘聚力C i 、摩擦角φi 及设计安全系数F s 。

录入数据后，需要计算的项目有滑体几何特征参数滑面长L i 、滑面倾角αi ，各滑块重力Q i 、滑动分力T i 、抗滑分力R i 、传递系数ψi ，传递系数Πψj ，传递抗滑力R i Πψj ，累计传递抗滑力∑R i Πψj ，传递滑动分力T i Πψj ，累计传递滑动力∑T i Πψj ，剩余下滑力P i 和滑坡稳定系数F s 。

Ｅxcel中提供了方便的编辑公式功能，我们可以根据计算项目的计算式录入公式。

有时为了引用某一固定的单元格，我们可以使用Excel中的绝对引用（＂$＂）功能。

各计算公式用的数学表达式和Ｅxcel 表达式见表1：剩余推力法计算法的数学表达式及其相应Ｅxcel表达式表1（3）利用Ｅxcel的复制功能对各条块的计算项目进行求解通过利用Excel提供的句柄复制功能可以对各计算项目进行快速的求解。

当鼠标指向选定的单元格出现句柄标志时，向下拖动该标志完成对该列进行公式的复制，则出现每一行对应的该项计算结果，对以上列的每一行对应每一条块的剩余下滑推力，各列进行复制，则出现全部结果，如图2、图3所示。

其中Pi其余各列为中间结果。

（4）计算滑坡的稳定系数Fs的计算式根据滑坡稳定系数计算公式，引用计算工作表中相关单元格数据，进行计算，稳定系数Fs为：＂B1＂=(P17+L17)/(R17+K17)通过以上步骤，我们即可以完成以前繁琐的滑坡稳定性评价中的计算，而且概念清晰，计算过程一目了然。

4. 地基沉降量计算《建筑地基基础设计规范》（GBJ7-89）规定的最终沉降量计算式如下［5］：()∑=---==ni i i i i sis s z z E p s s 1110'ααψψ式中： s —地基最终沉降量（mm ）；s ’—按分层总和法计算出的地基沉降量（mm ）；ψs —沉降计算经验系数；n —地基沉降计算深度范围内所划分的土层数；p 0—对应荷载标准值时的基础底面处的附加压力（kPa ）； E si —基础底面下第i 层土的压缩模量（Mpa ）；z i 、z i-1—基础底面至第i 层土、第i+1层土地面的距离（m ）；图2. 剩余推力法计算工作表图图3. 剩余推力法计算工作表图iα、1-i α—基础底面计算点至第i 层土、第i ＋1层土底面范围内平均附加应力系数。

若为矩形基础上的均布荷载，平均附加应力系数i α的计算公式如下［6］：⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡+++⎪⎭⎫ ⎝⎛-+⎪⎭⎫ ⎝⎛+++⎪⎭⎫⎝⎛++⎪⎭⎫ ⎝⎛-+++⎪⎭⎫ ⎝⎛-+⎪⎭⎫ ⎝⎛+++⎪⎭⎫ ⎝⎛++⎪⎭⎫ ⎝⎛-++=22222222222222111111111ln 1111ln21n m z L arctg m n m m n m z L m m m n m m m m n m z B i i i i πα其中m=L/B ，n=z i /B 。

式中： L —矩形基础的长边长度（m ）；B —矩形基础的短边长度（m ）；z i —意义同前。

通过上式可以方便地计算出矩形面积上均布荷载作用下角点平均附加应力系数α，由于计算公式较复杂，我们可以用EXCEL 中的Visual Basic 语言编辑器创建一个用户自定义函数来实现。

Visual Basic 编辑后的函数在插入函数中的“用户定义”显示，用户可以像使用内部函数一样直接使用该函数。

求矩形面积上均布荷载作用下角点平均附加应力系数α函数的程序如下：Function alfa(zi, l, b, m, n) If zi = 0 Thenalfa = 0.2500 Elsex1 = Sqr(1 + m * m + n * n) x2 = Sqr(1 + m * m) x3 = l / x1 / zix = (x1 - m) * (x2 + m) / (x1 + m) / (x2 - m) y = (x1 - 1) * (x2 + 1) / (x1 + 1) / (x2 - 1)alfa = 0.5 * (b * Log(x) / zi + l * Log(y) / zi + Atn(x3)) / 3.14159 End IfEnd Function以陈仲颐主编的，清华大学出版社出版的《土力学》P134页的例题（如图4所示）为例［7］，介绍用Excel 求矩形面积上均布荷载作用下角点沉降量∑△S i'的计算方法：（1） 设计如图5的工作表格工作表中需要录入的数据有计算深度Z （m ），基础长度L （m ），基础宽度B （m ），基底附加压力P 0（kPa ），第i 层土的压缩模量E si （kPa ），需要计算的有m ，n ，平均附加应力系数αi ，αi Z i ，αi Z i －αi-1Z i-1，△S i '及累计沉降量∑△S i '。

（2） 对工作表各计算栏插入公式和函数对于需要计算的m ，n ，αi Z i ，αi Z i －αi-1Z i-1，△S i '及累计沉降量∑△S i '，它们的计算公式简单，可以直接输入：m ： E3=B3/C3， n ： F3=A3/C3，αi ： 引用定义函数；αi Z i ： H3=G3*A3， αi Z i －αi-1Z i-1： I4=H4-H3， △ S i '： J4=1000*$B$1/D4*I4， ∑△S i '：K4=K3+J4。

平均附加应力系数αi 直接利用上面定义的函数来求，选择“菜单栏”中的“函数”，在弹出的“粘贴函数”菜单中函数分类“用户定义”的函数名一栏，会发现有一个“地基沉降.xls!alfa.alfa ”函数，根据函数中各参数的含义引用相关的单元格即可完成平均附加应力系数的计算。

（3） 对各计算栏复制公式求结果同样利用句柄的复制功能进行公式复制，即可求出各栏的计算结果，单元格K11即为∑△S i'的计算结果。

然后根据规范的规定确定沉降计算经验系数ψs ，两数相乘就可以求出矩形面积上均布荷载作用下角点的最终沉降量。

最后按叠加原理计算出建筑物地基的最终沉降量。

6．结束语从以上两例可以看到，尽管通过编程可以完成各类岩土工程计算，但EXCEL 有其独特的优点： 1.其计算过程是透明的，其原始数据，中间结果及最终结果都显示在一张表格中，一目了然。

2. EXCEL 可以完成重复性计算，对对需要反复试算的问题或同类工程计算，只需改变一下有关参数，结果即可产生。

3.我们可以利用EXCEL 的模板功能，将平时工作中计算的工作簿制成模板，以利于以后反复使用。

可见利用EXCEL 进行岩土工程计算可以大大的提高我们的工作效率。

参考文献：1． 梁秀娟 赵兴民 Excel 在水利水电工程地质勘察中的应用 《工程地质计算机应用》 2000；3；21～242．何仲伟王辉用Excel求解土坝安全鉴定中的土工试验数据统计问题《工程地质计算机应用》2000；3；25～283．［美］MARK DODGE等著《Microsoft EXCEL for Windows 95 使用指南》清华大学出版社1996；35～484．中华人民共和国国家标准《岩土工程勘察规范》（GB50021—94）中国建筑工业出版社1995；69～71 5．中华人民共和国国家标准《建筑地基基础设计规范》（GBJ7—89）中国建筑工业出版社1996；22～24 6．《工程地质手册》编写委员会《工程地质手册》（第三版）中国建筑工业出版社1994；972～978 7．陈仲颐主编《土力学》清华大学出版社 1996；134～136。