FLAC3D在滑坡数值分析及监测预警设计中的应用

１４７８地下空间与工程学报第２卷

料达到屈服极限时的塑性流动现象。由于ＦＬＡＣ３Ｄ采用了显示差分格式和混合离散技术，可以精确的模拟材料的屈服和塑性流动。它的三维网格可以由用户根据所模拟物体的形状进行调整，网格由若干六面体单元组成，每个单元在指定的力和边界条件下遵循指定的线性或非线性本构关系，如果单元应力使得材料发生屈服或是塑性流动，单元网格可以在大变形的计算模式下发生很大的变形‘３｜。

２八字门滑坡概况

八字门滑坡为一土质古滑坡，滑坡分布高程１３９～２８０ｍ，面积１２万ｍ２，滑坡体积４００万ｍ３。滑坡体位于侵蚀构造中低山山前斜坡，地形上前凸中平后凹，以上陡下缓的连续斜坡形态兼具三级堆积平台地貌。滑坡体内冲沟较发育，滑坡后壁高程２５０ｍ，滑坡璧较陡，坡度４０。～６０。，其剖面形态呈凹槽状。滑坡前缘为人长江河口附近的香溪河，河道较直，受三峡库区蓄水影响，目前水位已达１３９ｍ，形成了水面宽阔的库区。

滑坡体簸箕状展布于岸坡坡脚，长３８０ｍ，宽１００～５００ｍ，厚１０—３５ｍ。

滑体主要为结构松散的崩坡积、坡积物。滑带在上部松散堆积层与下部基岩面间连续分布，滑带物质以粉质粘土含角砾为主，颜色较杂。滑床基岩在１３９ｍ蓄水位以上主要为侏罗系下统（Ｊ，）岩层，水位以下主要为三叠系中统（Ｔ２）岩层。

目前，八字门滑坡具有发生失稳的可能性。

滑坡变形影响因素：

２．１三峡水库蓄水

在三峡水库蓄水１７５ｍ后，滑坡的滑体大部分将处于库水位以下，地下水位强烈上升；同时在库水位１４５～１７５ｍ变化大幅度涨落期间，地下水产生的渗透动水压力，会加剧滑坡的整体失稳。三峡水库蓄水后即库水位调节是影响八字门滑坡稳定的最重要因素。

２．２降雨

连续的降雨易使滑坡体处于饱和状态，自重不断增加，地下水的下渗，由于滑面相对隔水，上层地下水在此集中。使滑面的强度大大降低，有可能诱发滑坡失稳。据调查近几年来，在连续暴雨之后，在滑体中前部多次发生小型的塌滑，在滑坡后缘和边界产生裂缝。从地表ＧＰＳ的监测点水平位移～时间曲线来看，２００４年５月一８月变形速率增大，而５月～８月正是八字门地区暴雨频繁的季节。说明降雨是影响八字门滑坡稳定性的重要因素。２．３人类工程活动

移民兴建住房，增大了上部荷载；复建公路的修建、移民开垦耕地，破坏植被，给地表水的下渗创造了条件。

３滑坡数值模拟

３．１计算模型

计算模型的长（ｘ轴方向，沿河流方向为正）为４２０ｍ，宽（Ｙ轴方向，垂直河流向上方向为正）为４３０．５ｍ，厚度（Ｚ轴方向，垂直向上为正）６０ｍ，见图ｌ。仅考虑自重力的作用，不考虑构造应力场，滑坡四周设法向约束，底部设固定约束。

图ｌ八字门滑坡ＦＬＡＣ３Ｄ计算模型

Ｆｉｇ．１ＴｈｅＦＬＡＣ３ＤｍｏｄｅｌｏｆＢａｚｉｍｅｎｌａｎｄｓｌｉｄｅ３．２参数选取

在边坡的数值分析过程中，滑坡物质物理力学参数的选取对计算的结果准确与否至关重要。滑坡合理参数的选用，需要通过多种手段联合实现。一般是通过实验室内土工试验结合以往的工程实践通过反算分析取得的。更加准确的参数要通过现场大型原位剪切试验取得，而这种试验一般会受到时间、资金以及现场试验条件的限制而无法进行。

关于合理参数的选取，田斌”ｕ教授提出的方法是：综合考虑滑坡体稳定现状及发展趋势、滑坡体物质组成、启动成因、工程地质与水文地质条件、室内实验结果、反算结果、人工神经网络预测等方法得出的结果，加以修正后使用。

ＦＬＡＣ３Ｄ计算模型范围内各岩土体的物理力

学参数取值如表ｌ。

２００６年第８期王振勉，等：ＦＬＡｃ３Ｄ在滑坡数值分析及监测预警设计中的应用１４７９

表１计算模型参数

Ｔａｂｌｅ１Ｍｅｃｈａｎｉ∞ｌｐａｒ锄ｅｔｅ璐ｏｆｍａｔｅｒｉａｌｓ

ｉＩｌｃａｌｃｌｌｌａ廿ｏｎｍＯｄｅｌ

粼称ｐ。篡，，鬻鬻泊警鬻篙掌基岩２８００１．０８ｅ１０１．５６ｅ１０Ｏ．２２２．４５ｅ６４２滑带１９９０３．３７９ｅ６３．２６７ｅ７０．２１．４７ｅ４１６

滑体２０００５．２５ｅ６２．４５ｅ７０．２５９．８ｌｅ３２２３．３模拟结果及分析

以上工作完成以后，开始对天然状态滑坡进行三维数值模拟。首先对滑坡进行初始状态的应力平衡，这个过程相当于给边坡施加初始地应力。计算过程中使所有部分服从弹性本构关系，可以使得计算的速度大大提高。再次将本构关系调整为摩尔一库仑模型，恢复滑坡各部分原来的物理力学参数，进行天然状态的计算。对滑坡的天然状态进行模拟可以得到天然状态下滑坡的应力、变形、塑性区分布。整个计算过程历时６０００时步。３．３．１位移分析

由数值分析得出滑坡区位移分布如图２所示。可以看出，滑坡体位移集中在浅表层的中部，位移方向平行坡面向下，其值在滑坡体中心部位最大，并沿滑坡中心往周围逐渐减小。

对比图３（八字门滑坡ＧＰｓ监测点水平位移～时间曲线）和图２可知，运用ＦＬＡＣ３Ｄ得到的位移变化范围为ｏ．１２５—０．３５ｍ，滑坡体上的４个监测点所得到的累计位移量为Ｏ．２６ｍ～Ｏ．３ｍ，二者的数据并不十分吻合，这是因为，Ｆｌａｃ３ｄ计算出的变形结果是在自重条件、天然状态下得出的，而实际监测结果是在自重平衡后综合降雨、人类活动等因素引起的变形。

图２滑坡体位移等高线

Ｆｉｇ．２Ｄｉｓｐｌａｃｅｍｅｎｔｃｏｎｔｏｕｒｓｏｆｌａｎｄｓｌｉｄｅ

ＯＦ碉

ｌ十ｚＧＬ１（｜

ｌ十ｚＧｌｌ』

１２１１鱼！！Ｉ

图３八字门滑坡ＧＰＳ监测点水平位移～时间曲线

Ｆｉｇ．３

Ｖ撕ａｔｉｏｎｏｆ

ｄｉｓｐｌａｃｅｍｅｎｔｏｆＧＰＳｍｏｎｉｔ０Ｉｉｎｇｐｏｉｎｔｓ证ｔｈｔｉｍｅ

３．３．２应力及塑性区分析

天然状况下ｙｙ方向应力分布图见图４。由图４可知，滑坡体浅表面大部分存在着拉应力分布区，特别在中下部拉应力相对较大；垂直方向上压应力值随深度的增加而同步递增。数值分析结果表明滑坡体应力场分布符合斜坡应力场普遍分布规律。

图５为剪切塑性区分布图。可以看出，剪切塑性区范围较广，主要分布在滑坡体后缘、前缘以及滑坡体右侧区域。

４结论

由以上分析可以得出以下结论。

（１）滑坡体变形较大的区域集中在表层中部，这说明滑坡体中部是变形的关键部位，在监测预警设计工作中，应把位移监测点重点布设在这一区域，以便真实反应滑坡体变形情况。

（２）滑坡体中下部拉应力相对较大，剪切塑性

凇耋｜枷｜瑚

㈣

１４８０地下空间与工程学报第２卷

区范围较广，主要分布在滑坡体后缘、前缘以及滑

坡体右侧区域。滑坡监测预警设计中推力孔和测

斜孔主要布置在这些区域，滑坡体的应力状态可以

得到真实反应。

（３）滑坡监测点得到合理布设，从而为下一步

的预测预报提供可靠的研究数据。

参考文献：

图４）ｒｙ方向应力等高线［２］

Ｆｉｇ．４Ｓｔｒｅｓｓｃｏｎｔｏｕｒｓａｌｏｎｇｙｙ

图５剪切塑性区分布图

Ｆｉｇ．５Ｄｉｓｔｒｉｂｕｔｉｏｎｏｆｓｈｅａｒｆａｉｌｕｒｅｚｏｎｅ［３］

［４］

冯晓等．长江三峡库区危岩滑坡监测预警体系规划

研究［Ｊ］．工程勘察，２００２，（２）：４３—４６

张雪东等．呷爬滑坡稳定性的ＦＬＡｃ一３Ｄ数值模拟

分析［Ｊ］．岩土力学，２００３，（１０）：１１３一１１６

ＦＬＡＣ，ＦａｓｔＬ删舀锄Ａｎａｌｙｓｉｓ０ｆＣｏｎｔｉｎｕａ，Ｖｅｒｓｉｏｎ４．

Ｏ，ｕｓｅｒ’ｓｍａｎｕａｌ［Ｍ］．ＩｔａｓｃａｃｏｎｓｕｌｔｉｎｇＧｒｏｕｐ，Ｉｎｃ．

ＵＳＡ，２００５，５

田斌等．滑带土结构强度特征及其强度参数取值研

究［Ｊ］．岩土力学与工程学报，２００４，２３（１７）：２８８７—

２８９２

【精品】第9章边坡稳定性分析

第9章边坡稳定性分析 学习指导：本章介绍了边坡的破坏类型,即:岩崩和岩滑;着重介绍了边坡稳定性分析与评价基本方法,包括圆弧法岩坡稳定分析、平面滑动法岩坡稳定分析、双平面滑动岩坡稳定分析、力多边形法岩坡稳定分析及近代理论计算法；介绍了岩坡处理的措施。 重点：1边坡的变形与破坏类型； 2影响边坡稳定性的因素； 3边坡稳定性分析与评价. 9。1边坡的变形与破坏类型 9。1.1概述

随着社会进步及经济发展，越来越多地在工程活动中涉及边坡工程问题，通过长期的工程实践，工程地质工作者已对边坡工程形成了比较完善的理论体系,并通过理论对人类工程活动，进行有效地指导。近年来，随着环境保护意识的增加及国际减轻自然灾害十年来的开展，人类已认识到：边坡诞生不仅仅是其本身的历史发展，而是与人类活动密切相关；人类在进行生产建设的同时，必须顾及到边坡的环境效应，并且把人类的发展置于环境之中,因而相继开展了工程活动与地质环境相互作用研究领域，在这些领域中，边坡作为地质工程的分支之一，一直是人们研究的重点课题之一。 在水电、交通、采矿等诸多的领域，边坡工程都是整体工程不可分割的部分，为保证工程运行安全及节约经费，广大学者对边坡的演化规律、边坡稳定性及滑坡预测预报等进行了广泛研究。然而，随着人类工程活动的规模扩大及经济建设的急剧发展,边坡工程中普遍出现了高陡边坡稳定性及大型灾害性滑坡预测问题。在我国，目前的露天采矿的人工边

坡已高达300—500m，而水电工程中遇到的天然边坡高度已达500—1000米，其中涉及的工程地质问题极为复杂，特别是在西南山区，边坡的变形、破坏极为普遍,滑坡灾害已成为一种常见的危害人民生命财产安全及工程正常运营的地质灾害。

滑坡稳定性定量分析法(最新)

打造最便宜 滑坡稳定性定量分析方法 目前，滑坡稳定性分析和工程治理主要是依据工程地质类比、自然历史分析、工程地质力学分析、极限平衡力学计算、弹塑性有限元计算等进行的，且在一定的程度上都有一定的实效性和可靠性。滑坡是一个复杂的、非线性的动态系统，且大型滑坡规模大、机制复杂、破坏性强，不仅失稳影响范围广，而且防治难度高、治理措施复杂。采用工程地质类比、历史反演和地质力学分析，需弄清地层结构、地质构造、地壳演化历史等问题。通过对滑坡形成的地质环境条件、影响因素、变形破坏及形成机制等特征的综合性分析,滑坡堆积体在天然状态下处于稳定状态, 在连续降雨、暴雨影响下处于基本稳定状态。在连续降雨、暴雨及地震等影响下处于欠稳定状态。 一、传统的稳定系数法。 稳定系数预测法是最早的滑坡空间预测方法，它是基于极限平衡法理论提出来的，是将有滑动趋势范围内的边坡土体沿某一滑动面切成若干竖条或斜条,在分析条块受力的基础上建立整个滑动土体的力 或力矩平衡方程,并以此为基础确定边坡的稳定安全系数。这些方法均假设土体沿着一个潜在的滑动面发生刚性滑动或转动。简化的极限平衡法有瑞典法，Bishop法、Spencer法,Janbu法, Sarma法等。通过计算滑坡体的安全系数Fs，来预测边坡的稳定性。 Fs=F抗滑力/F下滑力 当Fs<1.0,不稳定状态; 当Fs=1.0,临界状态; 当Fs>1.0,稳定状态。 二、数值分析方法。 ①有限单元法 有限元法是目前使用最广泛的一种数值分析方法。优点是部分地考虑了边坡岩体的非均质和不连续性,可以给出岩体的应力、应变大小与分布；避免了极限平衡分析法中将滑体视为刚体而过于简化的缺点；能近似地从应力应变去分析边坡的变形破坏机制，分析最先、最容易发生屈服破坏的部位和需要首先进行加固的部位等。但是对于大的变形和位移不连续问题的求解还不理想。 ②离散单元法 离散单元法是处理结构控制型岩体工程问题较成熟方法。该程序不但允许有限位移和离散体的转动及脱离,而且在计算过程中可以自动判别块体之间可能出现新的接触关系,因此它可以方便地实现对复杂结构体变形破坏的模拟，可以将所研究的区域划分为一个个多边块体单元，单元之间通过接触关系,建立位移和力的相互作用规律,通过迭代使得每一个块体都达到平衡状态。在稳定分析中,它的功能在于反映岩块之间接触的滑移、分离和倾翻等大位移的同时,又能计算岩块内部的变形与应力,该法的另一个优点是利用显式时间差分解求解动力平衡方程,可方便地求解非线性大位移和动力稳定。 ③统计分析方法。 这是目前国内外研究人员研究滑坡稳定性使用较多的一类方法。统计分析方法建立在对滑坡影响因子和滑坡分布关系的分析之上,因此,它能最大程度反映滑坡分布与致灾因子之间的关系,使地质灾害危险性评价更加趋近于客观现实。包括信息量法、多元统计方法、聚类分析方法等。 三、瑞典法的基本理论 瑞典圆弧滑动法是条分法中最古老而又最简单的方法。除了假定滑裂面是个圆柱面外, 在求条底反力时忽略了条间力的作用, 且在求安全系数时仅考虑对同一点的力矩平衡。其安全系数方程为:

数值分析课后题答案

数值分析 第二章 2．当1,1,2x =-时，()0,3,4f x =-,求()f x 的二次插值多项式。 解： 0120121200102021101201220211,1,2, ()0,()3,()4;()()1 ()(1)(2)()()2()()1 ()(1)(2) ()()6 ()()1 ()(1)(1) ()()3 x x x f x f x f x x x x x l x x x x x x x x x x x l x x x x x x x x x x x l x x x x x x x ==-===-=--==-+-----==------= =-+-- 则二次拉格朗日插值多项式为 2 20 ()()k k k L x y l x ==∑ 0223()4() 14 (1)(2)(1)(1)23 537623 l x l x x x x x x x =-+=---+ -+= +- 6．设,0,1,,j x j n =L 为互异节点，求证： （1） 0()n k k j j j x l x x =≡∑ (0,1,,);k n =L （2）0 ()()0n k j j j x x l x =-≡∑ (0,1,,);k n =L 证明 （1） 令()k f x x = 若插值节点为,0,1,,j x j n =L ，则函数()f x 的n 次插值多项式为0 ()()n k n j j j L x x l x == ∑。 插值余项为(1)1() ()()()()(1)! n n n n f R x f x L x x n ξω++=-= + 又,k n ≤Q

(1)()0 ()0 n n f R x ξ+∴=∴= 0()n k k j j j x l x x =∴=∑ (0,1,,);k n =L 0 000 (2)()() (())()()(()) n k j j j n n j i k i k j j j i n n i k i i k j j i j x x l x C x x l x C x x l x =-==-==-=-=-∑∑∑∑∑ 0i n ≤≤Q 又 由上题结论可知 ()n k i j j j x l x x ==∑ ()()0 n i k i i k i k C x x x x -=∴=-=-=∑原式 ∴得证。 7设[]2 (),f x C a b ∈且()()0,f a f b ==求证： 21 max ()()max ().8 a x b a x b f x b a f x ≤≤≤≤''≤- 解：令01,x a x b ==，以此为插值节点，则线性插值多项式为 10 101010 ()() ()x x x x L x f x f x x x x x --=+-- =() () x b x a f a f b a b x a --=+-- 1()()0()0 f a f b L x ==∴=Q 又 插值余项为1011 ()()()()()()2 R x f x L x f x x x x x ''=-= -- 011 ()()()()2 f x f x x x x x ''∴= --

滑坡稳定性分析知识讲解

滑坡稳定性分析

习题一岩村滑坡稳定性评价 一、目的 学会滑坡机理分析、稳定性定价和定量计算的基本方法，了解滑带土抗剪强度指标选择的基本途径，掌握滑坡防治工程要点。 二、滑坡概况 l、自然地理 岩村滑坡位于四川盆地某城市市中区，地处长江和佳江的交汇地带，呈半岛状，土地资源十分紧张。在经济建设迅速发展的80年代，市中区斜坡土地得到了大量的利用，交通线路不断改进，高层建筑逐渐增多。但与此同时滑坡灾害事件也日趋严重，岩村滑坡就是灾害之一。 该地区属于亚热带气候，温暖潮湿，雨量充沛，多年平均降雨量在1200mm以上，并常有暴雨出现。长江和嘉陵江是市中区两大地表水系，水位年平均变化幅度达20m以上，平均低水位158m，高水位181m，1981年为百年一遇的特大洪水，水位达193m。三峡工程按175m高程修建大坝，使该地区最高洪水位达205m左右。 2、地质概况 滑坡区基岩地质构造属川东隔档式褶皱中的一复向斜内部，岩层产状平缓，倾角10°以下，倾向在SW200°～270°范围变化。无明显的断裂构造，优势节理产状：75°∠82°；346°∠81°，263°∠85°。 基岩地层为侏罗系泥岩砂岩互层，为内陆河潮沉积，呈紫红色。相对坚硬的砂岩组成了滑坡区的上部平台状地形，泥岩及崩积物则组成斜坡主体。崩积物主要由砂岩块石及泥岩风化粘土组成，厚度分布特点是斜坡上部薄，中前部相对较厚。人工堆石为近期在砂岩体中开挖地下洞室而堆弃于斜坡后部的基岩大块石。

滑坡区属河流侵蚀、剥蚀的低山丘陵地貌，斜坡顶部为平台，河谷岸坡的坡度由上至下逐渐变缓，在纵剖面上呈内凹的地形。 下伏基岩相对不透水，为弱含水层。据洞室调查，基岩洞室绝大多数为干洞，偶见裂隙有渗水现象。斜坡地带入渗的地表水则汇集于基岩顶面，形成崩积层中的上层滞水。 该地区新构造运动不强烈，属受活断裂包围的稳定地块，地震基本烈度为Ⅵ度。 3、滑坡特征 滑坡主滑方向为NW方向，后缘有一系列NE-SW方向的拉张裂缝，居民建筑物受到严重影响。据调查，人工洞室开挖于1970-1980年之间，地面裂缝最早发现在1981年。1981年四川盆地普降暴雨，江河水位达百年一遇特大水位。滑坡的活动已严重威胁经由滑坡区的主干公路的正常通车。滑坡现处于蠕滑阶段，且在每年的雨季，位移明显增大。 表1-1钻孔地质描述

公路边坡稳定性评价方法及滑坡防治措施实用版

YF-ED-J2674 可按资料类型定义编号 公路边坡稳定性评价方法及滑坡防治措施实用版 In Order To Ensure The Effective And Safe Operation Of The Department Work Or Production, Relevant Personnel Shall Follow The Procedures In Handling Business Or Operating Equipment. （示范文稿） 二零XX年XX月XX日

公路边坡稳定性评价方法及滑坡 防治措施实用版 提示：该解决方案文档适合使用于从目的、要求、方式、方法、进度等都部署具体、周密，并有很强可操作性的计划，在进行中紧扣进度，实现最大程度完成与接近最初目标。下载后可以对文件进行定制修改，请根据实际需要调整使用。 引言 近年来，随着国民经济的飞速发展，“村 村通公路”工程的进一步实施，在地形困难路 段修建的公路越来越多。受各种条件的限制， 大填、大挖方路段频繁出现，相伴而来出现了 较多的路堤边坡失稳，边坡及路堑边坡坍塌等 地质灾难现象，给公路建设、运营带来巨大的 经济损失。因此在公路建设中需要选用合理的 方法评价其边坡稳定性，根据评价结果确定合 理的边坡治理措施进而做到既保证公路运营的

安全，又节约投资。由此看来，稳定性评价的方法显得至关重要。本文对边坡稳定性评价方法和滑坡防治措施进行研究，为二程技术人员在实际工程中选用合理的评价方法和防治措施提供参考。 1、公路边坡病害的分类 边坡病害可分为以下3类。 1、1滑坡 滑坡是路基山坡土体或岩体由于长期受地下水、地表水活动的影响使其结构逐渐失去支撑力，在自重的作用下，整体沿着一定软弱面向下滑动。滑坡按其引起滑动的力学特性来区分，可分为牵引式和推移式滑坡。牵引式滑坡是下部先滑动，使上部失去支撑而变形滑动，一般速度较慢，可延续相当长时间，横向张性

滑坡稳定性分析

习题一岩村滑坡稳定性评价 一、目的 学会滑坡机理分析、稳定性定价和定量计算的基本方法，了解滑带土抗剪强度指标选择的基本途径，掌握滑坡防治工程要点。 二、滑坡概况 l、自然地理 岩村滑坡位于四川盆地某城市市中区，地处长江和佳江的交汇地带，呈半岛状，土地资源十分紧张。在经济建设迅速发展的80年代，市中区斜坡土地得到了大量的利用，交通线路不断改进，高层建筑逐渐增多。但与此同时滑坡灾害事件也日趋严重，岩村滑坡就是灾害之一。 该地区属于亚热带气候，温暖潮湿，雨量充沛，多年平均降雨量在1200mm以上，并常有暴雨出现。长江和嘉陵江是市中区两大地表水系，水位年平均变化幅度达20m以上，平均低水位158m，高水位181m，1981年为百年一遇的特大洪水，水位达193m。三峡工程按175m高程修建大坝，使该地区最高洪水位达205m左右。 2、地质概况 滑坡区基岩地质构造属川东隔档式褶皱中的一复向斜内部，岩层产状平缓，倾角10°以下，倾向在SW200°～270°范围变化。无明显的断裂构造，优势节理产状：75°∠82°；346°∠81°，263°∠85°。 基岩地层为侏罗系泥岩砂岩互层，为内陆河潮沉积，呈紫红色。相对坚硬的砂岩组成了滑坡区的上部平台状地形，泥岩及崩积物则组成斜坡主体。崩积物主要由砂岩块石及泥岩风化粘土组成，厚度分布特点是斜坡上部薄，中前部相对较厚。人工堆石为近期在砂岩体中开挖地下洞室而堆弃于斜坡后部的基岩大块石。 滑坡区属河流侵蚀、剥蚀的低山丘陵地貌，斜坡顶部为平台，河谷岸坡的坡度由上至下逐渐变缓，在纵剖面上呈内凹的地形。

下伏基岩相对不透水，为弱含水层。据洞室调查，基岩洞室绝大多数为干洞，偶见裂隙有渗水现象。斜坡地带入渗的地表水则汇集于基岩顶面，形成崩积层中的上层滞水。 该地区新构造运动不强烈，属受活断裂包围的稳定地块，地震基本烈度为Ⅵ度。 3、滑坡特征 滑坡主滑方向为NW方向，后缘有一系列NE-SW方向的拉张裂缝，居民建筑物受到严重影响。据调查，人工洞室开挖于1970-1980年之间，地面裂缝最早发现在1981年。1981年四川盆地普降暴雨，江河水位达百年一遇特大水位。滑坡的活动已严重威胁经由滑坡区的主干公路的正常通车。滑坡现处于蠕滑阶段，且在每年的雨季，位移明显增大。 表1-1钻孔地质描述 表1-2岩土体物理力学性质指标 表1-3滑带土抗剪强度指标实验值

数值分析习题集及答案

（适合课程《数值方法A 》和《数值方法B 》） 第一章 绪 论 1. 设x >0,x 的相对误差为δ,求ln x 的误差. 2. 设x 的相对误差为2％,求n x 的相对误差. 3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位 有效数字: ***** 123451.1021,0.031,385.6,56.430,7 1.0.x x x x x =====? 4. 利用公式(3.3)求下列各近似值的误差限: * * * * * * * * 12412324(),(),()/,i x x x ii x x x iii x x ++其中* * * * 1234,,,x x x x 均为第3题所给的数. 5. 计算球体积要使相对误差限为1％,问度量半径R 时允许的相对误差限是多少? 6. 设028,Y =按递推公式 11783 100 n n Y Y -=- ( n=1,2,…) 计算到100Y .若取783≈27.982(五位有效数字),试问计算100Y 将有多大误差? 7. 求方程2 5610x x -+=的两个根,使它至少具有四位有效数字(783≈27.982). 8. 当N 充分大时,怎样求 2 11N dx x +∞+?? 9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2 ? 10. 设2 12S gt = 假定g 是准确的,而对t 的测量有±0.1秒的误差,证明当t 增加时S 的绝对误差增加, 而相对误差却减小. 11. 序列{}n y 满足递推关系1101 n n y y -=-(n=1,2,…),若02 1.41y =≈(三位有效数字),计算到10 y 时误差有多大?这个计算过程稳定吗? 12. 计算6 (21)f =-,取 2 1.4≈,利用下列等式计算,哪一个得到的结果最好? 3 6 3 11,(322), ,9970 2. (21) (322) --++ 13. 2 ()ln(1)f x x x =- -,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一等 价公式 2 2 ln(1)ln(1)x x x x - -=-+ + 计算,求对数时误差有多大? 14. 试用消元法解方程组{ 10 10 12121010; 2. x x x x +=+=假定只用三位数计算,问结果是否可靠? 15. 已知三角形面积 1sin , 2 s ab c = 其中c 为弧度, 02c π << ,且测量a ,b ,c 的误差分别为,,.a b c ???证 明面积的误差s ?满足 . s a b c s a b c ????≤ ++ 第二章 插值法 1. 根据( 2.2)定义的范德蒙行列式,令

工程地质数值法-FLAC3d模拟

《工程地质数值法》课程论文 论文题目：水窖开挖工程模拟 专业方向：建筑与土木工程 姓名： 学号： 2015年9月1日

1. 工程概况 在地面上开挖一个尺寸为2m ×4m ×5m 的存水地窖，其横截面图如图1。要求分析开挖后土体受力情况。 图1 睡觉开挖横截面图 2. 分析目的 利用有限单元法分析得出梁以下各个单元不同受力状况，计算出梁下应力分布状况。 3. 计算参数 将地下6m ×8m ×8m 的土体作为计算范围，水窖尺寸为2m ×4m ×5m ，G=100MPa ，3m kg 2000-?=ρ，35.0=ν。 4. 有限元模型 计算模型（见图）

建立一个6m×8m×8m，的立方体模型，将整个水窖包含其中，然后取其中一半进行开挖模拟，分析开挖后应力应变及位移变化情况。施工过程中在底部位置监测开挖时的各种变化。具体命令流入下： new gen zone brick p0=(0,0,0) p1=(6,0,0) p2=(0,4,0) p3=(0,0,8) size 6,4,8 plo blo group model mohr ini dens 1000 fix x range x -0.1 0.1 fix x range x 5.9 6.1 fix y range y -0.1 0.1 fix y range y 3.9 4.1 fix z range z -0.1 0.1 prop bulk 1e8 shear 0.3e8 fric 35 prop coh 1e10 tens 1e10 set grav 0 0 -9.81 solve save 初始应力.sav model null ran z 3 8 x 2 4 y 0 2 prop bulk 3e8 shear 1e8

滑坡稳定性分析计算

对最不利滑移横断面进行各种工况稳定性分析计算，计算过程如下： 一、天然工况 滑坡剩余下滑力计算 计算项目：滑坡推力计算 1 ===================================================================== 原始条件: 滑动体重度= 19.000(kN/m3) 滑动体饱和重度= 25.000(kN/m3) 安全系数= 1.250 不考虑动水压力和浮托力 不考虑承压水的浮托力 不考虑坡面外的静水压力的作用 不考虑地震力 坡面线段数: 6, 起始点标高 4.000(m) 段号投影Dx(m) 投影Dy(m) 附加力数 1 13.600 0.700 0 2 12.250 7.000 0 3 2.000 0.000 0 4 12.000 8.000 0 5 24.500 0.500 0 6 127.000 27.000 0 水面线段数: 1, 起始点标高 0.000(m) 段号投影Dx(m) 投影Dy(m) 1 0.000 0.000 滑动面线段数: 5, 起始点标高 0.000(m) 段号投影Dx(m) 投影Dy(m) 粘聚力(kPa) 摩擦角(度) 1 12.000 0.600 10.000 14.500 2 9.900 1.300 10.000 14.500 3 28.000 9.000 10.000 14.500 4 8.400 2.800 10.000 14.500 5 117.000 29.000 10.000 14.500 计算目标：按指定滑面计算推力 -------------------------------------------------------------- 第 1 块滑体

数值分析课后题答案

数值分析 2?当x=1,—1,2时，f(x)=O, 一3,4,求f(x)的二次插值多项式。解： X 0 =1,x j = — 1,x 2 = 2, f(X。)= 0, f (xj = -3, f (x2)= 4; l o(x)=(x-xi^~x2\=-1(x 1)(x-2) (x o -X/X o _x2) 2 (x -x0)(x -x2) 1 l i(x) 0 2(x-1)(x-2) (x i ~x0)(x i ~x2) 6 (x—x0)(x—x,) 1 l2(x) 0 1(x-1)(x 1) (X2 -X°)(X2 - X i) 3 则二次拉格朗日插值多项式为 2 L 2(X)= ' y k 1 k ( x) kz0 = -3l°(x) 4l2(x) 1 4 =(x_1)(x—2) 4 (x-1)(x 1) 2 3 5 2 3 7 x x - 6 2 3 6?设Xj, j =0,1,||(,n 为互异节点，求证： n (1 )7 x：l j(x) =x k(k =0,1川,n); j=0 n (2 )7 (X j -x)k l j(x)三0 (k =0,1川,n); j ￡ 证明 (1)令f(x)=x k

n 若插值节点为X j, j =0,1,|l(, n，则函数f (x)的n次插值多项式为L n(x)八x k l j(x)。 j=0 f (n 十)(?) 插值余项为R n(X)二f(X)-L n(X) n1(X) (n +1)!

.f(n1)( ^0 R n(X)=O n 二瓦x k l j(x) =x k(k =0,1川,n); j ：o n ⑵、(X j -x)k l j(x) j卫 n n =為(' C?x j(—x)k_L)l j(x) j =0 i =0 n n i k i i =為C k( -x) (、X j l j(x)) i =0 j=0 又70 _i _n 由上题结论可知 n .原式二''C k(-x)k_L x' i=0 =(X -X)k =0 -得证。 7设f (x) c2 la,b 1且f (a) =f (b)二0,求证: max f(x)兰一(b-a) max a $至小一*丘f (x). 解：令x^a,x^b，以此为插值节点，则线性插值多项式为 L i(x^ f(x o) x x f (xj X o —人x -X o X —X o x-b x-a ==f(a) f(b)- a - b x -a 又T f (a) = f (b)二0 L i(x) = 0 1 插值余项为R(x)二f (x) - L,(x) f (x)(x - X Q)(X - xj 1 f(x) = 2 f (x)(x -X g)(X -xj

FLAC3D滑坡模拟

FLAC 3D滑坡模拟 一、源程序 ; Create Material Zones gen zone brick size 5 5 5 & p0 (0,0,0) p1 (3,0,0) p2 (0,3,0) p3 (0,0,5) & p4 (3,3,0) p5 (0,5,5) p6 (5,0,5) p7 (5,5,5) gen zone brick size 5 5 5 p0 (0,0,5) edge 5.0 group Material ; Create Bin Zones gen zone brick size 1 5 5 & p0 (4,1,0) p1 add (3,0,0) p2 add (0,3,0) & p3 add (2,0,5) p4 add (3,6,0) p5 add (2,5,5) & p6 add (3,0,5) p7 add (3,6,5) gen zone brick size 1 5 5 & p0 (6,1,5) p1 add (1,0,0) p2 add (0,5,0) & p3 add (0,0,5) p4 add (1,6,0) p5 add (0,5,5) & p6 add (1,0,5) p7 add (1,6,5) gen zone brick size 5 1 5 & p0 (1,4,0) p1 add (3,0,0) p2 add (0,3,0) & p3 add (0,2,5) p4 add (6,3,0) p5 add (0,3,5) & p6 add (5,2,5) p7 add (6,3,5) gen zone brick size 5 1 5 & p0 (1,6,5) p1 add (5,0,0) p2 add (0,1,0) & p3 add (0,0,5) p4 add (6,1,0) p5 add (0,1,5) & p6 add (5,0,5) p7 add (6,1,5) group Bin range group Material not ; Create named range synonyms range name=Bin group Bin range name=Material group Material ; Assign models to groups model mohr range Material model elas range Bin ; Create interface elements int 1 face ran plane ori (4,0,0) nor (-5,0,2) dist 0.01 z (0,5) y (1,6) int 2 face ran plane ori (0,4,0) nor (0,-5,2) dist 0.01 z (0,5) x (1,6) int 1 face ran x 5.9 6.1 y 1 6 z 5 10 int 2 face ran x 1 6 y 5.9 6.1 z 5 10 int 1 maxedge 0.55 int 2 maxedge 0.55 ; Move bin toward material ini x add -1.0 range Bin ini y add -1.0 range Bin ; Assign properties

滑坡稳定性影响因素及分析

滑坡稳定性影响因素及分析 滑坡是在一定的内因、外因等地质环境条件和其它因素综合作用下产生的，影响因素包括：地质条件、地形地貌、人类活动、气候及迳流条件、其它因素。就本滑坡隐患体而言，各因素对其的影响如下： ①地质条件 岩土体的本身特性是影响边坡稳定性的主要因素；对岩质边坡来说主要包括软弱结构面存在与否及其强度、结构面特别是主要结构面的产状、结构面的组合关系、结构面的结合情况、渗透性、与临空面的相对关系；对土质边坡来说主要包括土体强度、软硬接触面的渗透性。滑坡隐患体及边坡出露的地层为泥盆系佘田桥组，岩性为砂岩，受地形地貌、构造侵蚀、剥蚀及风化作用影响，第四系及土状风化物厚度变化较大；原始地形较平缓的人工切坡坡面及坡顶局部地段第四系及土状风化物厚度大。第四系坡残积土其孔隙性大且含较多碎石，抗剪强度较低，坡度较陡时其自稳性差；中上部基岩埋藏多较浅且表部风化较强烈；整个山体岩体裂隙发育，地层及裂隙产状较杂乱(图2-1)，地层产状多近坡向或与坡向小角度斜交，岩体呈碎裂结构、电阻较高，结构面结合多数差~较差，易产生松动变形。 ②地形地貌因素 勘查区属中低山地貌，高差较大，山脊地形坡度较陡（坡度25~30°），两侧地形陡峻(坡度40~45°)，但从调查情况来看，沟谷处及外围天然斜坡未见有滑坡现象，天然条件下斜坡是稳定的；但切坡以后，山体前缘产生高陡临空面，所形成的上缓下陡地形不利于斜坡的稳定。 ③人类活动因素 人类工程活动破坏原有的地形地貌，使在自然条件下已经达到平衡状态的岩土体应力进行重新分布，斜坡产生变形，当岩土体中应力无法平衡时，边坡将发生失稳破坏。就本区而言，切坡产生高陡地形，

滑坡勘查中滑坡稳定性分析实例

滑坡勘查中滑坡稳定性分析评价实例 中国建筑材料工业地质勘查中心河南总队吴德运 关键词：滑坡稳定性安全系数稳定状态 滑坡地质灾害每年均会给社会造成较大的人员伤亡和财产损失，滑坡的产生受多种引发因素影响，往往也是多种因素叠加的结果。如何准确分析滑坡的稳定性是治理滑坡的关键。本文是以一个滑坡实例，评价滑坡稳定性的分析过程。 1 滑坡区自然条件及地质环境条件 1.1 自然条件 该滑坡处于中纬度带，属亚热带季风气候区，多年平均降雨量1100mm，最大年降雨量1522.4mm，最小年降雨量694.8mm。5～9月为雨季，其降雨量占全年降雨量的70%以上。一小时最大降雨量达75.2mm，一日最大降雨量达193.3mm。 1.2 地质环境 1.2.1 地形地貌 滑坡区属鄂西中低山地貌单元。由于地壳长期间歇性抬升，形成山高坡陡、河谷深切的地貌特征。 1.2.2 地层岩性 滑坡区分布的地层有： 第四系：残坡积碎石土、残坡积堆积土。 三叠系中统：中厚至厚层微晶白云质灰岩、泥灰岩、中厚层泥质条带灰岩、肉红色中厚层亮晶鲕状灰岩及灰绿色泥岩。岩层产状总体向北东向倾，倾角为35o－70°之间。 1.2.3 水文地质条件 受地层岩性结构和地质构造影响，滑坡区内地下水主要以三叠系中统岩溶裂隙水和第四系松散岩孔隙水的形式赋存。 2．滑坡基本特征及类型 2.1 滑坡地形地貌 滑坡区地形南高北低，地形总坡度15o－20o，为侵蚀构造低山区。滑坡区最低点标高330m，最高点滑坡后缘，标高364m，相对高差34m。

2.2 滑坡空间形态 该滑坡为覆盖层滑坡，平面形态呈舌形，地形上为围椅状，滑坡两边周界清晰。滑坡体北低南高，主滑坡轴线长86m，前缘宽98m，标高330m ，后缘宽66m，标高364m。滑坡的面积为0.732×104m2，总体上是前厚后薄，中间厚两侧薄的态势，滑体平均厚度为5m，体积约3.66×104m3。 滑坡主滑方向为311度，滑体坡度15～30度，中部滑坡平台呈舒缓波状，中部靠后缘出现陡坎。 2.3 滑坡物质组成及结构特征 (1)滑体 滑体物质组成主要为第四系崩坡积碎块石夹粉质粘土，黄褐－黄灰色，稍密－中密，碎块石直径一般为0.4－0.8m，最大达1.2m，成分主要为泥灰岩、灰岩，其含量约占70%。滑体厚度一般为2.3-6.7m。 (2)滑带 滑带主要成分为粉质粘土夹砾石，灰黄－褐黄色，粉质粘土呈可塑状，含量约70%，具有挤压条纹状构造，砾石成份为泥灰岩、灰岩，呈次棱角状－次圆状，直径2～20mm。部分砾石表面见擦痕，表面具滑感。 (3)滑床 滑床为三叠系中统泥灰岩，强～中风化程度，浅灰－黄灰色，中厚层～厚层状构造，岩石较为破碎，地层倾向为19～40度，倾角41～75度，岩石节理裂隙发育，裂隙面倾角为60～75度，裂隙面均较平直，略具起伏，稍粗糙，多为泥质、铁质充填，部分为钙质充填。 2.4 滑坡水文地质 本滑坡地下水主要为第四系覆盖层松散岩类孔隙水和基岩裂隙水。 覆盖层孔隙水水量贫乏，赋水性弱，主要接受大气降水次为农作物灌溉渗入补给。地下水沿基岩面排泄，或渗入下伏基岩裂隙中。基岩浅部裂隙发育，含裂隙水，赋水性弱，动态变化大。补给主要靠覆盖层地下水渗入，排泄主要受微地貌控制，流量小。 2.5 滑坡岩土物理力学性质 2.5.1滑体岩土物理力学性质 滑体主要由第四系崩坡积碎块石夹粘性土组成，碎石含量达70％以上，受取样条件限制，滑体中采取的原状样土工试验所作的物理力学指标仅能代表碎石土中所夹粉

数值分析第四版习题及答案

第四版 数值分析习题 第一章绪论 1.设x>0,x得相对误差为δ,求得误差、 2.设x得相对误差为2％,求得相对误差、 3.下列各数都就是经过四舍五入得到得近似数,即误差限不超过最后一位得半个单位,试指 出它们就是几位有效数字: 4.利用公式(3、3)求下列各近似值得误差限: 其中均为第3题所给得数、 5.计算球体积要使相对误差限为1％,问度量半径R时允许得相对误差限就是多少? 6.设按递推公式 ( n=1,2,…) 计算到、若取≈27、982(五位有效数字),试问计算将有多大误差? 7.求方程得两个根,使它至少具有四位有效数字(≈27、982)、 8.当N充分大时,怎样求? 9.正方形得边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝? 10.设假定g就是准确得,而对t得测量有±0、1秒得误差,证明当t增加时S得绝对误差增 加,而相对误差却减小、 11.序列满足递推关系(n=1,2,…),若(三位有效数字),计算到时误差有多大?这个计算过程 稳定吗? 12.计算,取,利用下列等式计算,哪一个得到得结果最好? 13.,求f(30)得值、若开平方用六位函数表,问求对数时误差有多大?若改用另一等价公式 计算,求对数时误差有多大? 14.试用消元法解方程组假定只用三位数计算,问结果就是否可靠? 15.已知三角形面积其中c为弧度,,且测量a ,b ,c得误差分别为证明面积得误差满足 第二章插值法 1.根据(2、2)定义得范德蒙行列式,令 证明就是n次多项式,它得根就是,且 、 2.当x= 1 , -1 , 2 时, f(x)= 0 , -3 , 4 ,求f(x)得二次插值多项式、 3. 4., 研究用线性插值求cos x 近似值时得总误差界、

土质滑坡稳定性分析

土质滑坡稳定性分析 影响滑坡稳定性的因素有很多，其中对滑坡稳定性影响较大的因素有降雨和地震，不同条件下滑坡的稳定性是不同的。文章以圆弧条分法分析了汶川地震灾区某滑坡的稳定性，结合现场的工程地质勘察，计算了滑坡的安全系数，分析不同条件下滑坡的稳定性，并给出相应的处理意见。 标签：滑坡稳定性；地震；降雨；稳定性分析 引言 5.12汶川地震发生后，诱发了为数众多的崩塌、滑坡、泥石流等次生地质灾害，这些重大地质灾害隐患点险情紧迫、危害巨大、危险程度高，严重危及着城区居民生命财产安全。文章结合地震区的某土质滑坡，运用圆弧条分法，分析了在自重、降雨、地震不同的情况下滑坡的稳定性[1-5]。 1 地质环境条件 1.1 地形地貌 勘查区位于白龙江南侧，属河谷地貌，位于白龙江一级阶地上。微地貌位于凸出的五山岭山脊两侧，总体地势中部高，东西两侧低，西侧（左侧）地形较平缓，东侧（右侧）地形起伏大。该滑坡前缘位于一冲沟的丘间梯田，沟底部分地段基岩出露，地面高程为611.50～618.00m；滑坡后缘为五山岭山脊的平坝边缘，地面高程为631.70～631.90m，相对高差约为13.00～20.00m，地势较为平缓。整体坡度角一般为20～30°。 1.2 地层岩性 勘查区基岩出露较差，仅在滑坡左侧冲沟边有出露。主要出露地层为第四系人工填土、冲洪积粉质粘土、卵石土及志留系黄坪组下段千枚岩（Shn1），现就与工程密切的地层由新至老简述如下： （1）第四系。第四系松散土层主要为冲洪积粉质粘土及卵石土层（Q4al+pl）。冲洪积粉质粘土，厚度一般约3m，最厚段可达6.50m，主要分布于五山岭山顶及两侧斜坡一带；冲洪积卵石土层，厚度较大，一般20～30m，分布于整个勘查区。 （2）基岩。工作区内基岩主要为志留系黄坪组下段（Shn1），其岩性主要为千枚岩，岩体较破碎，表层风化较严重，强度较低。 1.3 地质构造及地震

FLAC3D原理..

2、2 三维数值模拟方法及其原理 2、2、1 FLAC3D工程分析软件特点 FLAC3D就是由美国Itasca Consulting Group, Inc、为地质工程应用而开发得连续介质显式有限差分计算机软件。FLAC即Fast Lagrangian Analysis of Continua 得缩写。该软件主要适用于模拟计算岩土体材料得力学行为及岩土材料达到屈服极限后产生得塑性流动,对大变形情况应用效果更好。 FLAC3D程序在数学上采用得就是快速拉格朗日方法,基于显式差分来获得模型全部运动方程与本构方程得步长解,其本构方程由基本应力应变定义及虎克定律导出,运动平衡方程则直接应用了柯西运动方程,该方程由牛顿运动定律导出。 计算模型一般就是由若干不同形状得三维单元体组成,也即剖分得空间单元网络区,计算中又将每个单元体进一步划分成由四个节点构成得四面体,四面体得应力应变只通过四个节点向其它四面体传递,进而传递到其它单元体。当对某一节点施加荷载后,在某一个微小得时间段内,作用于该点得荷载只对周围得若干节点(相邻节点)有影响。利用运动方程,根据单元节点得速度变化与时间,可计算出单元之间得相对位移,进而求出单元应变,再利用单元模型得本构方程,可求出单元应力。在计算应变过程中,利用高斯积分理论,将三维问题转化为二维问题而使其简单化。在运动方程中,还充分考虑了岩土体所具有得粘滞性,将其视作阻尼附加于方程中。 FLAC3D具有一个功能强大得网格生成器,有12种基本形状得单元体可供选择,利用这12种基本单元体,几乎可以构成任何形状得空间立体模型。 FLAC3D主要就是为地质工程应用而开发得岩土体力学数值评价计算程序,自身设计有九种材料本构模型: (1)空模型(Null Model) (2)弹性各向同性材料模型(Elastic, Isotropic Model) (3)弹性各向异性材料模型(Elastic, anisotropic Model) (4)德拉克-普拉格弹塑性材料模型(Drucker-Prager Model) (5)莫尔-库伦弹塑性材料模型(Mohr-Coulomb Model)

渝黔高速向家坡滑坡特征与稳定性分析

渝黔高速向家坡滑坡特征与稳定性分析 林道刚吴汉辉杨转运刘会 摘要：总结了以往对向家坡滑坡治理工作经验教训，并通过定性的工程地质分析和定量计算评价了向家坡滑坡稳定性的变化趋势，提出改建公路超挖深切坡角、持续长时间强降雨形成的空隙水压力以及滑坡内部的物质组成是导致古滑波复活进而发生滑动破坏的主要因素，结论指出该滑坡治理中需考虑膨胀岩的膨胀性和加强内外排水 向家坡古滑坡位于重庆一黔江高速公路K13+500-K13+960段．距重庆市南岸区四公里(四公里是重庆市的一个地名)以东3 km，有机耕道直接通往该古滑坡。交通较为便利。该古滑坡形成于1998-05-18．受当时当地连续暴雨的影响，该滑坡出现大面积跨塌．在其后缘出现圆弧形张拉裂缝，经有关专家现场踏勘后，确定该滑坡为浅层牵引式土层滑坡。体积为8 910 m3，提出了在山顶修建截排水沟(3号、4号排水沟)、地下盲沟及在K13+560-K13+840段左侧设置抗滑桩(般桩长8-10 m，最大桩长11 m，总共55根)的综合治理方案。由于原设计桩的深度不够．1999年6月，发现K13+570-K13+720段路堑边坡顺路线方向又出现长达140 m的多条裂缝，原边坡出现了较大的滑移变形，山顶部分桩(K13+600-K13+850段)有明显的位移，再次补充评价后认为该滑体范围大，土层较厚，一般为5-12 m，最大厚度可达19 m，滑坡的体积(100-180)×104 m3，为一大型滑坡。根据该勘察报告，施工方案改为：在山顶加一排抗滑桩共17根，深度一般18-20 m，中部按1：1-1：1．25的坡比分两级削坡，并采用格子梁加锚索锚固。在坡脚设48根锚拉（Kl3+593~K13+922段)，两端用抗滑挡墙加固，在桩后作钢筋砼挡板，同时在前面施工片石砼挡墙。并在坡脚再增设15根锚索桩，在半坡设12根抗滑桩和在锚固桩上设4根锚杆加固锚固桩，锚索桩间用C20钢筋砼现浇挡墙，工程于2002年11月通过验收。但2004年6~8月，由于渝黔高速公路施工开挖，在坡顶又发现新的裂缝。前排桩以向公路倾斜为主，半坡桩出现向坡外倾斜、下沉的现象，滑坡区仍然在缓慢变形阶段，直接危及临近公路及立交桥正常营运。因此，对向家坡滑坡成因机制、滑坡特征与稳定性进行研究，并对其实施工程治理，这将对维护渝黔高速公路的正常营运和保持社会稳定起到积极作用。同时也是实施科学合理加固的重要基础。 1滑坡区工程地质环境条件 1．1 地形地貌 滑坡区内地势东高西低，为阶梯状斜坡地貌，上陡下缓，下部斜坡自然坡度20°～30°：上部斜坡坡度大于40°，局部为陡崖；由于公路的修建及滑坡的前期治理，现状地形可明显分为3级台阶。滑坡上部缓坡平台地形坡度18°，高程385~400 m；中部呈阶梯状，前缘为高12~14 m的桩锚挡土墙，公路为勘测区的最低点，高程一般342~346 m。滑坡后壁为高约40 m的陡壁，南侧和北侧为冲沟。滑坡前缘为走向为SN的路堑边坡。 1．2地层岩性 滑坡区分布的地层为第四系全新统、侏罗系下统珍珠冲组(J1z)、自流井组（J1-2Z)和三叠系上统须家河组(T3xj)。滑坡岩土层见工程地质综合剖面图1。

某滑坡稳定性分析

清平水库瓦窑堡滑坡稳定性分析 杨荣科,辜明清 (四川省水利水电勘测设计研究院勘察分院，四川郫县611731) 摘要：瓦窑堡滑坡是位于清平水库坝址上游左岸的一个大型古滑坡，水库蓄水后的滑坡稳定性评价是水库区重大工程地质问题之一。根据大量的勘察试验资料，分析了滑坡的成因和形成机制，利用反演计算进行了滑坡稳定分析评价。 关键词：滑坡；抗剪强度；稳定性；清平水库 1引言 瓦窑堡滑坡是位于清平水库枇杷岩坝址上游1.4 km左岸的大型滑坡，水库蓄水后该滑坡的稳定性是近坝库岸的主要工程地质问题。分析滑坡的工程地质条件，针对滑坡形成机制，采用反演进行稳定性评价，是对古滑坡稳定性评价较适用的方法。 2滑坡基本特征和工程地质条件 瓦窑堡滑坡地面高程884～1 155 m。据地表地质测绘，滑坡体长约450 m，宽290～450 m，厚30～65 m，体积约364×104 m3。滑坡体在平面上呈“板斧”形，两侧以冲沟为界，下游侧缘冲沟切割至滑坡床基岩，沟深3～10 m，沿滑面无地下水点出露，见图1。滑坡后缘地形坡度30°～45°，并见张开5～15 cm 的拉裂缝；中部地形平缓，坡度12°～30°，呈阶梯形；前缘剪切口明显，与2al）接触。滑坡体总体地形坡向N60°～70°W。 河床砂卵石（Q 4 ）中厚层灰岩，下部滑坡区出露的地层主要有：二叠系上统长兴组（P 2c 常夹碳质页岩；龙潭组（P ）上部为碳质页岩夹煤层，下部为厚3.4～6.5 m的 2l ）为中厚层灰岩夹泥质灰岩；地表分布第四系坡粘土岩；二叠系上统茅口组（P 1m 积层（Q 4 dl ）。滑坡地段在构造上位于照壁山倒转向斜核部附近，有近 ）从滑坡后缘一带通过。瓦窑堡断裂走向北东，倾南北向断裂之瓦窑堡断裂（F 5 向北西，倾角54°左右，延伸约24 km，上、下盘均为灰岩，滑坡一带下盘为龙潭组之碳质页岩。断层破碎带一般厚10～40 cm，由断层角砾、挤压破碎透镜体等组成。