挡土墙电子表格开发报告 - 岩土工程研究所.pdfVIP

挡土墙电子表格开发报告 1 开发背景 在国家“十二五”期间，在科技支撑项目“重大水利工程建设信息化标准化专项技术开 发研究”的资助下，研发了挡土墙电子表格。挡土墙电子表格是在Excel 界面上，利用其强 大的计算功能并结合VBA 进行二次开发，研制的计算挡土墙后主动土压力的电子表格。该 程序主要功能是完成指定滑面的库仑主动土压力计算，并自动搜索出土压力极大值及对应的 临界滑裂面位置，可计算墙后填土为任意形式 （包括分层土、粘性土、有地下水等）的挡土 墙的库仑主动土压力。其基本原理是：将作用于支挡结构上土压力问题看作是一个具有垂直 表面并且在此表面作用有外荷载的边坡稳定问题，结合库仑土压力理论，在通用条分法基础 上，对墙后滑动土体划分土条，并计算各土条间的土压力，最终得到挡土墙后库仑主动土压 力；计算三个滑裂面的土压力值，做出抛物线，找到极大值 Pa 及对应的滑面位置。该计算 方法适用范围广，可以为水利及道路桥梁工程设计人员提供设计依据。 挡土墙电子表格是在边坡稳定分析电子表格 Lossap 基础上开发的，借用Excel 界面， 主要的计算都是在Excel 中输入公式进行的，方便用户进行核算；对于较为复杂的判断和计 算，利用Excel 的VBA 进行程序的二次开发完成。另外，针对该程序，创建了“程序菜单” 及相关操作按钮，界面友好，方便用户使用。 2 基础理论和规范、文献依据 2.1 分层土库仑主动土压力计算基本原理 在利用库仑土压力理论确定挡土墙后主动土压力时，通常假定墙后填土破坏时将沿着土 体内某一直线滑裂面滑动，通过分析由墙土结合面、滑裂面和土体顶部直线构成的三角形土 体的力的平衡，计算出土体作用在墙上的主动土压力P 。通过不断变动直线滑裂面位置，可 以找到使P 获得极大值Pa 的临界滑裂面，如图 1 所示。 1 挡土墙电子表格开发报告 图1 应用库仑理论计算主动土压力的方法 （a ）作用于一滑动土体上的力；（b ）静力平衡条件；（c ）确定临界滑裂面 这种在静力平衡条件下建立力的三角形的计算方法（如图 1(b) ），只适用于墙后填土为 无粘聚力的单层土的情况。此时，滑裂面下固定土体对滑动土体的支承反力R 与滑裂面的  夹角即为土体的内摩擦角 （如图 1(a) ）；而对于分层土，该支承反力R 应为各土层的支承 反力R 的合力，此时R 值及与滑裂面夹角难以确定，因此，无法利用库仑土压力理论解决 i 分层土主动土压力的计算问题。另外，该方法只在墙后填土为无粘性土时提供了数值解，对 于粘性土，只能通过图解法进行计算，计算过程复杂，因此实际工程少用。 Chen Li (1998)[1] 曾提出一个建立在通用条分法基础上计算主动土压力的数值分析方 法。通过力平衡方程求得土压力值，以力矩平衡方程确定主动土压力作用点位置，再通过改 变滑裂面位置来确定极大值 Pa 及对应的临界滑裂面位置。该方法可以解决任意形式滑裂面 的分层土(可以包含粘性土)的主动土压力计算问题。 由此，如果将通用条分法引入到库仑土压力理论中，将墙后土体划分土条，引入条间力 倾角的假定，分条计算土压力值，就可以解决分层土的库仑主动土压力的计算问题。 2.1.1 公式推导 设想某一边坡的滑动土体沿滑裂面y = y (x)下滑，如图 2 所示。对滑动土体进行垂直条 分，取出其中一个土条，分析作用在其上的力，主要有以下几种： 1) 土条重量W，浸润线上为天然容重，浸润线下为饱和容重。 2) 坡表面垂直荷重qx 。 3) 地震力，水平地震力Q =W，其作用点与土条底距离为h 。 e 4) 作用在土条垂直边上的总作用力G （即土骨架间的法向有效作用力和水压力之和）， 它与水平线的夹角为，其作用点的纵坐标值为y t 。 对土条建立x 和y 方向的静力平衡方程，得到静力平衡的微分方程：