[0037] 下面结合实施例对本发明作进一步说明,但不应该理解为本发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下,根据本领域普通技术知识和惯用手段,做出各种替换和变更,均应包括在本发明的保护范围内。
[0038] 实施例1:
[0039] 参见图1,本实施例公开基于传递系数法和表面位移的滑坡稳定性评价方法,包括以下步骤:
[0040] 1)参见图2,将待评价滑坡划分为n个垂直滑块,并由上到下依次标记为1、2、……、i、……、n。
[0041] 参见图3,对于任意一个滑块i,在不考虑外部作用时,
[0042] 滑块i的法向方向受力平衡关系式为:
[0043] Ni=Wicosαi+Pi-1sin(αi-1-αi) (1)
[0044] 滑块i的下滑力总和为:
[0045] Ti=Wisinαi+Pi-1cos(αi-1-αi) (2)
[0046] 滑块i的抗滑力为:
[0047]
[0048] 式中,Ni为滑块所受正压力,kN;αi为第i个滑块与水平面的夹角,°;Wi为滑块的重力,kN;Bi为滑块的宽度,m;Pi-1为第i块滑块受到上一滑块的剩余推力,即该滑块受到的不平衡力,kN;Ti为第i块滑块的下滑力,kN;ci为滑面的黏聚力,kPa; 为滑面的内摩擦角;Sfi为滑块抗滑力,kN。
[0049] 在传统的传递系数法中,通常通过对抗滑力进行折减或对下滑力进行扩大的方法,进而求得稳定性系数。但这种方式无法得到变形数据。本实施例为了建立滑坡变形与稳定性之间的关系,对传统传递系数法的基本假定进行以下改进:
[0050] a.假设滑块能够发生变形,认为滑块底部(即滑带)的变形满足理想弹塑性剪应力—剪位移本构模型,其极限剪应力满足莫尔库伦破坏准则。
[0051] b.认为各滑块在上部推力作用下发生平动的基础上,还发生了侧向和竖向的变形,上一滑块对下一滑块的剩余推力并不是集中力作用,而是分布力。
[0052] c.对抗滑力不进行折减,也不对下滑力进行扩大。
[0053] 因此,剩余推力Pi可表示为:
[0054] Pi=Ti-Sfi (4)
[0055] 2)根据滑块表面位移实测数据和MINDLIN解,计算各个滑块底部滑带的剪切位移。
[0056] 对于任一滑块,其表面水平位移△x等于滑带的水平位移sx加上不平衡力对滑块侧边单独引起的水平位移ux。同理,其表面竖向位移△y等于滑带的竖向位移sy加上不平衡力对滑块侧边单独引起的竖向位移uy。
[0057] 参见图4,明德林(mindlin)理论是在布西奈斯克(Boussinesq)解答的基础上推导出竖向或水平向荷载作用在半无限弹性体内部时体内任一点的应力与应变的数学解。
[0058] 在实际工程中,滑块所受不平衡力引起滑块变形。对于传递系数法中的任意滑块,其受到的不平衡力就是上部滑块对其的剩余推力与作用于该滑块上其它外力的合力。对于传递系数法中划分的条块通常除第一个滑块和最后一个滑块是三角形外,其余滑块均为四边形(矩形或梯形)且侧边为竖直状。因此,可引入MINDLIN解,得到半无限弹性体内深度c处作用水平向和垂直向集中力P时,体内任一点的应力与应变的数学解。
[0059] 当弹性半无限体内部深度c处作用水平集中力Q时,则离地面深度z处的任意一点M的侧向、竖向变形分别为:
[0060]
[0061]
[0062] 式中,Q为水平集中力,G为剪切模量,ν为泊松比,其余长度参数含义详见图4所示。
[0063] 滑坡中的除第一个三角形滑块外的任意一个滑块,认为其侧边受到上一滑块的推力作用,但推力方向并不是水平方向或竖直方向,因此需要进行力的分解。参见图5,滑块i的合力方向与滑块底部滑面平行,Pi-1为滑块所受不平衡力,将Pi-1分解成水平方向的分量Px和竖直方向的分量Py。此时,可根据MINDLIN解计算得到水平和竖直方向集中力作用下滑块侧边所产生的水平位移及竖向位移。
[0064] 2.1)通过对MINDLIN解进行简化,在二维情况下,水平集中力作用点所处竖直面上各点的变形计算公式即可表示为:
[0065]
[0066] uz=0 (8)
[0067] 从上式可知,在水平集中力的作用下,滑块侧面上的竖向位移为0,但水平位移不可忽略。
[0068] 将滑块上部推力考虑成三角形分布,根据式7,得到了水平集中力作用下滑块侧边变形的解析解。参见图6,按照滑块上部推力的三角形分布,需采用数值积分的方式进行计算。
[0069] 参见图7,7a表示滑块不平衡推力分布,7b表示各段受力分布,7c表示等效集中力。对于任意滑块i,其受到上部推力的水平分力为Px,将其分布为三角形荷载后,用m个节点将侧边等分为m-1段,那么对于第j个节点上的荷载大小即为:
[0070]
[0071] 式中,Px为该滑块所受上部推力的水平分力;m为滑块侧边节点划分数量,m>1;j为节点号;Hi为滑块i的侧边高度。
[0072] 将j-1节点和j节点的荷载(kPa)等效成集中力(kN),其值即为图7b中梯形的面积,如下式:
[0073]
[0074] 式中,hj为第j个节点与第j-1个节点之间的距离。
[0075] 将公式10代入公式9,可得:
[0076]
[0077] 根据上述分解,可获得侧边的m-1个等效集中力,再将式11代入式7,即可得到每个等效集中力作用范围下侧边1~m个节点上的水平位移。由于地表位移是关注的重点,则只需求得各段等效集中力对滑块侧边顶部节点(节点1)产生的水平位移,然后进行累加,进而得到该点的水平位移。因此,在采用公式7计算时,对于顶点的深度z取为0,而各等效集中力作用点深度cj则应根据分的段数进行取值:
[0078]
[0079] 对于顶点(节点1),任意一等效集中力对其产生的水平位移即为:
[0080]
[0081] 从上式可发现,任意一等效集中力对顶点处的水平位移是相同的,将m-1个等效集中力累加,得到滑块侧边在上部推力作用下顶点(节点1)的水平位移,其值仅与荷载大小、侧边高度以及土体参数有关,与节点数无关。
[0082]
[0083] 2.2)通过对MINDLIN解进行简化,在二维情况下,竖向集中力作用点所处竖直面上各点的变形计算公式可表示为:
[0084] ux=0 (15)
[0085]
[0086] 根据步骤2.1),可得到任意一等效竖向集中力对滑块侧边顶点产生的竖向位移:
[0087]
[0088] 从上式可发现,任意一等效集中力对顶点处的竖向位移也是相同的,将m-1个等效集中力累加,得到滑块侧边在上部推力作用下顶点(节点1)的竖直位移,
[0089]
[0090] 2.3)根据滑坡监测数据,可以得到滑块的表面位移实测数据。根据表面位移实测数据,计算得到滑块底部滑带的剪切位移。其中,滑块底部滑带剪切位移的水平分量为:
[0091] sx=Δx-ux (19)
[0092] 滑块底部滑带剪切位移的竖向分量为:
[0093] sy=Δy-uy (20)
[0094] 式中,△x为滑块总的水平位移,由滑块表面位移实测数据获取。△y为滑块总的竖向位移,由滑块表面位移实测数据获取。ux为不平衡力对滑块侧边单独引起的水平位移。uy为不平衡力对滑块侧边单独引起的竖向位移。
[0095] 3)根据步骤2)计算得到的滑带剪切位移,计算滑坡稳定性系数。采用稳定性系数评价滑坡的稳定性。
[0096] 引入滑带的理想弹塑性剪应力—剪位移关系曲线。当滑块所受下滑力等于抗滑力时,该滑块底部滑带处于B点状态,定义该滑块底部滑带处于临界剪位移状态。当滑块所受下滑力小于抗滑力时,该滑块底部滑带处于O-B线性阶段,定义该滑块处于未达到临界剪位移状态。当滑块所受下滑力大于抗滑力时,该滑块底部滑带处于B点以后状态,定义该滑块发生了剪切破坏。参见图8,对于任意滑块i,其底部滑带的剪应力—剪位移关系为:
[0097]
[0098] 式中,ai为剪切刚度,kPa/mm;si为剪切位移,mm;sif为滑块的临界剪切位移,mm。
[0099] 3.1)根据滑坡实际发生的位移与其临界位移,得到任意滑块i的位移稳定性系数为:
[0100]
[0101] 式中,sfi为滑块i底部滑带的临界剪位移,根据滑带土的剪切试验得到。si为滑块i底部滑带的现状剪位移。Li为滑块i底部滑带的长度。
[0102] 3.2)对各滑块的位移稳定性系数加权计算,得到滑坡整体的实时稳定性系数:
[0103]
[0104] 式中,αi为滑块i底边与水平面的夹角。uxi为不平衡力作用下滑块i侧边顶点的水平位移。uyi为不平衡力作用下滑块i侧边顶点的竖向位移。△xi为任意滑块侧边顶部总的水平位移,由滑块表面位移实测数据获取。△yi为任意滑块侧边顶部总的竖向位移,由滑块表面位移实测数据获取。
[0105] 3.3)采用稳定性系数评价滑坡的稳定性。