《边坡稳定性的整体分析法及其应用》由会员分享,可在线阅读,更多相关《边坡稳定性的整体分析法及其应用(44页珍藏版)》请在金锄头文库上搜索。
1、边坡稳定性的整体分析法及其应用报告提纲 工作概述及3D分析的意义; 3D分析的几种常见方法; 现有方法的局限性; 滑面法向应力分布的自然形式; 总体分析法; = 0的情形; 广义特征值问题; 算例工作概述 利用整体分析法实现了边坡三维严格极限 平衡法,证明了边坡稳定性分析实际上可 归结为一个代数特征值问题,而安全系数 则是该问题的最大特征值,从而解决了极 限平衡法中的解的不收敛问题 带侧向摩擦边界的三维滑体的分析方法 官地重力坝深层抗滑稳定性分析3D分析的意义3D分析的几种常见方法 瑞典法:忽略所有的条间力,仅满足力平衡条件 Janbu法:忽略垂直于滑动方向的条块间的作用 力,假定平行于滑动方
2、向的条块间的推力作用点, 满足3个力平衡和一个力矩平衡 Spencer法:忽略垂直于滑动方向的条块间的作 用力,假定平行于滑动方向的条块间的推力都相 互平行,满足3个力平衡和一个力矩平衡现有方法的局限性 平衡条件未被全部满足; 存在大量的缺乏物理基础的假定; 数值特性差; 条分化过程过于严格总体分析法 经典的严格条分法都是通过条块的力平衡 方程消去条块底部的法向力而将不确定性 放在条间力上,所以称其为局部法,这类 方法的好处是能充分发挥人们的经验和直 觉,但数值特性差:1)所有方法都是局部 收敛的;2)所有方法都有不收敛的算例。 Bell于1968年提出了实现严格极限平衡法的 另一途径:假定滑
3、面的应力分布。Bell法的 优点是无需对条间力做出假定,因平衡条 件是针对整个滑体而不是单个条块而列出 的,所以称之为整体分析法。整体分析法续1 整体分析法的优点是其良好的数值特性且易于实 现三维严格极限平衡法 但30年来整体法一直倍受冷遇,直到杨洪杰等 (2001)、朱大勇等(2002)重新发现了这一方法Bell, J M. General slope stability analysis J. Journal of the SoilMechanics and Foundations Division, ASCE, 94(SM6): 1253-1270 , 1968.Yang, H J, W
4、ang J H, Liu Y Q. A new approach for the slope stabilityanalysis. Mechanics Research Communications, 28(6): 653-669, 2001.Zhu D Y, Lee C F. Explicit limit equilibrium solution for slope stability.InternationalJournalforNumericalandAnalyticalMethodsinGeomechanics, 26(15):15731590 , 2002.滑面法向应力分布的自然形式
5、 xxx0hn2 300是对应与瑞典法的滑面法向应力ndS dS O kdw d dS 滑面法向应力矫正的其它方式Bell方式aba xxxxba2sin0杨洪杰方式cbxax2朱大勇方式 xblxalba0几种矫正方式的对比221: Bell6& Zhu et9(F = 1.573) 2: The proposed (F = 1.580) MP(half-sine): F = 1.588 Spencer: F = 1.58911S1S2S3(30, 25)(50, 35)(50, 35) (40, 27)(52, 24)31 & 23211:Bell6F = 1.376 2: Zhu et
6、al9F = 1.376 3: The proposed F = 1.390F = 1.287F = 1.334cut off滑面法向应力的逼近普通滑体 alxT51,llTl51,aaTa xxx034125滑体的平面投影带侧向摩擦边界的滑体的插值盒 alxT51,llTl51,aaTa xxx0平衡积分方程fcFufcFewee110extSdffF0extSCdmmm nxnnc sxssc extext mmff0SwmSedScFdSfFsfsn代 alxxT0入平衡方程组0,dcQaPaagFFFSTdSlnP56ST edSflsQ56SmdS nfc016SewdSfcsd01
7、6 = 0不固结、不排水工况 有许多3D方法是专门针对这一工况的,如 Ugai (1985), Gens (1987), 以及Duncan 1996 年综述性论文中的表3所列举的其他方法, 等等。对于此种工况,这里所建议的方法 能够给出安全系数的显示表达式 因此,这一结果目前是最好的 PcPd ,det,detF0,dcQaPaagFFF011 acPadQF广义特征值问题BxAxdQAcPBF)qzeig(,BADV广义特征值问题的Matlab求解:算例1-带侧向摩擦边界的简单算例zxyABCDFEHhWH = 10 m h = 5 m = 30 W is variableF v.s. W1
8、.721.741.761.781.801.821.841.861.881.901.921.94050100150200250300 W/HFX. Zhang的椭球体算例The profile of an example from Zhang (1988)(36.6, 27.4)24.4(42.7, 6.1)weak layer: ce=0,e=10(18.3, 18.3)ce=29.3kPa e=20 19.2kN/m34.55zx椭球体的表面网格剖分令节点5在 1234内变化,F的变化范围:1.7071.71712345 1234椭球体问题的所有特征解1.70660.03470.00000.
9、03720.0093+0.3441i0.00930.3441ix1.11220.37950.01571.00000.4806+0.0329i0.48060.0329i2.04951.00000.28470.36420.90340.0074i0.9034+0.0074i2.05830.99490.28710.35980.90320.0082i0.9032+0.0082i1.11330.37900.01550.99920.4802+0.0334i0.48020.0334i2.00570.75501.00000.39440.97090.0291i0.9709+0.0291i1.00000.51080
10、.24160.47720.0269+0.0038i0.02690.0038i楔形体算例OABC xyz10041532节点5在 1234内变化,F无变化变化:1.636(1.640) 其它5个特征值:2= 0.0355,3=4=5=6= 0.0.雅砻江官地水电站概况 官地水电站是雅砻江卡拉至江口河段水电 规划五级开发方式的第三个梯级电站。位 于四川省凉山彝族自治州西昌市与盐源县 交界。 拦河坝采用碾压混凝土重力坝,最大坝高 168m,装机容量2400MW,最大坝底宽度 153.2m。自左至右依次布置左岸挡水坝段 (9个 )、河床溢流坝段(6个 ) 、右岸挡水坝 段(9个 )。基本地质条件 官地
11、水电站坝基岩体内中、缓倾角结构面 较发育,分布及其产状随机性及性状差异 较大; 左右岸坝肩中上部岩体中次生夹泥分布普 遍,尤其部分中缓倾结构面中有较厚的次 生泥,均一性差。 这些特定的地质条件使坝体坝基系统的深 浅层抗滑稳定问题很突出 。坐标系的建立坐标系的建立 对所有的坝段采用统一 的坐标系。其中: X轴坐标与坝轴线坐标 一致,原点位移0+0.0 处; Y轴垂直于X轴指向下 游,原点位移(坝)0+0.0 处; Z轴竖直向上,坐标与 高程一致。 右图显示了该坐标系的 水平坐标轴。水平坐标系载荷组合载荷组合 分析中考虑了以下三组载荷工况: 1)工况1:自重+上下游静水压力+泥沙压力+滑 移面上的
12、扬压力; 2)工况2:工况1载荷+水平载荷系数为0.284的 设计地震载荷; 3)工况3:工况1载荷+水平载荷系数为0.3262的 校核地震载荷; 其中,上下游水面高程分别为1330m和1203.7m, 上游迎水面的泥沙作用高程为1212.9m。16#坝段的两种滑移模式坝段的两种滑移模式 沿上下游方向的滑移模式: 是以N2045W/NE1525(中缓倾角裂隙)为 底滑面,沿着下游的N35W/SW30中缓倾角裂隙剪出。沿下游方向的滑移模式滑向河床的滑移模式16#坝段所用材料及结构面参数坝段所用材料及结构面参数16#坝三维极限平衡法计算结果坝三维极限平衡法计算结果工况模式16Y模式16X工况14.
13、6393.495工况24.6002.688工况34.5492.5971)沿河流方向的滑移模式16Y的安全度高于滑向河床方向滑移模式16X 的安全度;但两个方向的安全度都较大; 2)地震对模式16Y的影响明显低于16X,这是与我们的经验相符合的;24#坝段坝段 建立了两种河床滑移模式的模型。 第1种滑移模式fx11N6080W/NE3550333N6080W/NE3550N7080E/SE7085近SN/W7085N6070E/SE7085近SN/W7085近SN/W7085fN2030W/SW6070p25-2坝0+000地质横剖面图124# 坝 块f23-1f23-2IV类岩体坝体III3类
14、岩体24#坝段坝段 第2种滑移模式IV类岩体坝体III3类岩体III3类岩体fx11N6080W/NE3550333fx12N6080W/NE3550N7080E/SE7085近SN/W7085N6070E/SE7085近SN/W7085近SN/W7085fN2030W/SW6070N60W/SW3040p25-2坝0+000地质横剖面图124# 坝 块fx13f23-1f23-224#坝段滑面特性坝段滑面特性24#坝段所用材料及结构面参数坝段所用材料及结构面参数24#坝段计算结果坝段计算结果表中列出了右坝肩连同表中列出了右坝肩连同24坝段滑向河床方向的三维安坝段滑向河床方向的三维安 全系数和
15、全系数和y = 8剖面的二维安全系数,其中二维安全系剖面的二维安全系数,其中二维安全系 数也是满足全部三个平衡条件的严格极限平衡解。数也是满足全部三个平衡条件的严格极限平衡解。 可见,可见, 24#坝段坝肩的稳定性是良好的,无需再添加坝段坝肩的稳定性是良好的,无需再添加 工程加固措施。工程加固措施。工况模型1模型2 2D3D2D3D工况12.2833.0202.1402.896工况21.9672.482D2.618U1.8772.488D2.602U工况31.9262.390D2.564U1.8422.436D2.562U三维和y = 8剖面的安全系数(上标D表示地震力的分量朝下,U朝上)水电工程之一-大岗山