用点接触模拟滑面的滑坡稳定性分析方法

结合有限元强度折减法,用点接触模拟滑面,提出了适用于已知滑裂面的滑坡稳定性分析及滑坡推力计算的数值计算方法。以滑坡体为分析对象,用切向和法向两组弹簧的弹簧力来模拟滑裂面上的剪切力和法向力;结合强度折减法和摩尔-库仑强度准则对滑面弹簧刚度进行迭代计算,求解得到了满足静力平衡的数值解。该方法避免了极限平衡条分法求解边坡稳定性时需对条间力进行诸多假设的局限,也避免了有限元强度折减法求滑坡推力时需加桩计算的不便及收敛困难等问题。新方法能计算出极限状态下滑坡体上的应力、滑裂面上的边界力和稳定性系数以及给定安全系数下滑坡体上任意截面的滑坡推力大小、作用方向和位置,具有较严格的理论基础,可为实际工程中滑坡的治理提供一定的参考。     

基于反算原理的滑坡推力简易估算

滑坡滑面的抗剪强度指标往往用反算法确定。分析表明,对已变形位移的滑坡,其推力可不反算滑面抗剪强度指标而直接按现状稳定系数和设计安全系数简易地估算,一般地,滑坡推力等于下滑分力乘以设计安全系数相对于现状稳定系数的提高比例值。进而讨论了其对不同设计工况、各种典型形状滑面的适用条件。这是一条估算滑坡推力的捷径。     

基于应变软化的多阶边坡稳定分析

多阶边坡除了最危险滑面需要关注,每级台阶的潜在滑面及其他不满足规范要求的滑面亦需要得到重视。以双台阶边坡和复杂多台阶边坡为例,基于应变软化理论,利用FLAC3D软件分析边坡的剪应变发展规律,根据剪应变的集中带获取边坡的多个潜在滑面。分析滑面的剪入、剪出和滑面中部点的强度参数变化规律,并结合矢量和法分析滑面在应变软化过程中安全系数的演化规律。研究结果表明：采用该方法可获取边坡的最危险滑面及次级滑动面;随着软化的发展,滑面上点的强度参数从峰值强度逐渐降低至残余强度,然而演化的速度不同;随着应变软化过程的发展,所有滑面的安全系数均逐渐降低,最后趋于一个定值。     

基于局部强度折减法的边坡多滑面分析方法及应用研究

边坡最危险滑动面的搜索方法一直是研究的热点,但边坡内部次级滑动面也可能不满足安全设计要求,因此,考虑边坡多条滑动面的分析方法亦应得到关注。在传统强度折减法中,对边坡整个区域的抗剪强度参数进行折减,此方法仅可得到一个临界滑动面和最小安全系数。提出了一种基于局部强度折减法的多滑面分析方法,即首先定义单元安全系数的概念,并且计算边坡每个单元的安全系数,然后自动搜索出单元安全系数处于不同范围内的单元集合,对各个单元集合的强度参数进行折减计算,即可得到不同安全系数对应的滑动面。通过单台阶和双台阶边坡算例验证了该方法的可行性,结果表明,随着安全系数的增大,潜在滑面的深度和潜在滑动区域亦增大。最后把该方法应用到某隧道进口仰坡的稳定性评价中, 通过该方法得到的多级滑动面与现场监测数据吻合较好。     

滑带土抗剪强度参数的三维反分析

以滑面正应力修正模式的三维边坡稳定性分析方法为基础,建立了确定滑带土抗剪强度参数的三维反分析模型。先对滑坡体的稳定性进行调查评估,确定滑坡的安全系数与三维滑面的几何尺寸;然后,通过对滑面正应力进行合理假定与修正,建立了平衡方程,以强度参数为未知数,利用数学方法并结合试验、经验等手段求解出抗剪强度参数。该模型克服了二维反分析将会高估滑带土的抗剪强度参数的缺陷,而且计算过程简单,易于编程实现。两个简单算例表明,三维反分析模型可以考虑三维滑体的空间效应,得到的抗剪强度参数更真实、更符合实际情况,为滑坡治理提供参考。     

碎石土滑坡的因素敏感性计算分析

通过资料搜集整理和分析、现场工程地质调查与勘探和室内外的物理力学试验,采用弹塑性接触有限元强度折减法、不平衡推力法和数理统计分析法,结合工程实例,分析了影响碎石土滑坡各因素的敏感性,结果表明,按因素敏感度从大到小排列,主要影响碎石土滑坡稳定性的因素依次为滑面岩土体内摩擦角、地形坡度、滑体饱水面积比和滑面岩土体的内聚力。地形坡度大、具有低内摩擦角滑面的碎石土滑坡的整体稳定性较差,而连续降雨、强降雨会抬高地下水位增大滑体饱水面积比,降低滑面带岩土体抗剪强度,增加下滑力,对碎石土滑坡的稳定性将产生最不利的影响。同时,在对碎石土滑坡进行加固处理时,必须优先考虑采取提高滑面等效内摩擦角、降低地形坡度结合降低地下水位措施,其次考虑采取提高滑面岩土体内聚力的措施。     

基于岩土材料软化特性的滑坡多级滑动面分析方法研究

在实际工程中,滑坡经常呈现出多级破坏的特征,而一般计算和设计中仅关注最危险的滑动面及对应的最小安全系数,这往往遗留安全隐患。在考虑岩土材料软化特性的基础上,借助FLAC3D和Matlab软件平台构建了无需人工干预便可有效模拟与评价滑坡多级破坏的理论框架。以旬阳县党家坝廉租房小区滑坡为例,通过对塑性剪应变、塑性拉应变、剪切应变增量等特征变量的渐进发展规律分析发现：滑坡多级滑动面的形成,其时间和空间顺序并不一定相对应。各级滑动面在时间上按产生的先后顺序依次为第1级主滑面、第2级主滑面、次级滑面;而在空间上按照从前到后的顺序依次为第1级主滑面、次级滑面、第2级主滑面。从计算所得的滑坡最终破坏形态中发现,坡体内部出现的滑动面条数与现场采集的拉裂缝数基本保持一致;第1级主滑面的入口位置与滑坡前缘拉裂缝位置基本吻合;而第2级主滑面及次级滑面的位置与现场勘查到的拉裂缝位置出现偏差。从材料参数的渐进发展规律中发现,滑面上强度参数的分布及大小均随着滑面的产生、发展而逐渐变化,其变化区间为峰值强度到残余强度,这是此方法可以有效模拟多级滑动面形成过程与各级滑动面之间相互影响的核心所在。借助矢量和法成功地实现了基于滑面上强度参数渐进发展规律的安全系数演化过程的确定。通过对安全系数演化过程的分析发现,在滑坡多级滑动面的形成过程中安全系数大小顺序呈现出3种不同的状态,这很好地揭示了各级滑动面在形成过程中的主、被动关系。     

基于强度折减概念的滑坡稳定性三维分析方法(I)：滑带土抗剪强度参数反演分析

滑带土抗剪强度是滑坡稳定性分析和防治工程设计中十分重要而又难于确定的参数之一。因此,基于临界状态假定的二维反分析方法得到了广泛应用。但是,自然界发生的滑坡基本上呈三维形态,其主滑方向有时变化也较大,使得有必要研究一种滑带土抗剪强度的三维反分析方法和滑坡稳定性的三维评价方法。笔者建议了一种基于强度折减概念的滑带土抗剪强度反分析方法,即通过逐步折减滑动面的强度参数,使滑动面的塑性区完全贯通,此时,塌滑体处于极限状态,所用强度参数即为滑带土的平均抗剪强度参数。从洪家渡水电站1#塌滑体计算成果来看,反演的滑带土摩擦角较二维反演值低4.1?,反映了滑坡体的三维效应,验证了所建议方法的可行性。     

抗滑刚架桩挡墙的设计与应用

对于滑面较深、推力较大的滑坡,抗滑刚架桩上设置挡墙会显著地减小抗滑刚架桩长径比,故抗滑刚架桩挡墙将会越来越广泛地应用于滑坡治理工程中。在文克尔弹性地基理论的基础上,本文采用有限差分法推导了抗滑刚架桩挡墙的计算公式。所提出的计算方法已应用于工程实际中,取得了良好效果。还明确了抗滑刚架桩挡墙的适用条件。通过与抗滑刚架桩的对比分析,得出了一些有益的结论,比如在相同的条件下,抗滑刚架桩挡墙的位移约是抗滑刚架桩的55%。     

基于强度折减法的某变电所边坡稳定性分析

为评价某山区变电所高填方场地稳定性,预测该场地布设支档结构后可能产生的滑裂面的形态,采用有限元强度折减法分析了该场地的稳定性。以现场原位测试及土工试验检测数据结合工程经验,综合确定了有限元计算参数。采用场地塑性应变增量贯通区作为场地上部失稳滑裂面,并采用极限平衡法对该滑面进行了验算。经分析计算,得到了场地支档结构墙顶水平位移与边坡安全系数的关系,发现极限平衡验算结果较有限元分析结果明显偏大,偏于不安全。     

下滑坡及防治方案

下滑坡是由采煤引发的巨型深层顺层滑坡,目前处在蠕动变形阶段。本文对滑坡稳定性进行了计算,算出不同采空区面积下的稳定系数,计算了滑坡推力,提出了地表地下排水、煤层采空区设承重阻滑键、滑坡前缘设抗滑桩综合防治方案,可供防治设计施工参考。     

对库仑土压力理论的若干修正

库仑土压力理论至今仍是计算土压力的重要方法而被人们所熟知。通过分析库仑土压力的墙后土楔体的受力特点,特别是深入研究了土楔体与墙的作用力关系,对库仑土压力理论给出了一些修正。认为土楔体和挡土墙之间的作用力（即定义的土压力）,并非一定要达到极限状态,所以不能确定土压力的作用方向,但土压力的作用方向必须在其允许的角度范围之内。所以,认为库仑主动土压力为作用方向角度变化范围内的最大值,库仑被动土压力为作用方向角度变化范围内的最小值。对于墙后土楔体,认为墙体和土楔体是两个不同物体,土楔体的形成是因为土中产生潜在破裂面,而原库仑土压力理论要求墙体与土之间也达到临界状态是不必要的。墙体对土楔体的作用力（即土压力）实质就是相当于一物（墙）施加于另一物（土楔体）的力,即使土楔体滑动了,两物之间也并非要滑动。推导了主动土压力计算公式,给出了被动土压力的近似计算方案。算例证明,计算结果与原库仑理论有明显不同。该研究对库仑土压力的修正和求解值得引起重视。     

基于有限元强度折减法的复式抗滑桩内力计算及桩土作用机理分析

在边坡、滑坡防治工程中已得到广泛应用的П形、h形等复合式抗滑桩,由于其多属于复杂的超静定结构,受力形式或荷载分配方式仍存在较大争议,桩土相互作用关系复杂,传统方法在求解其内力时往往面临较大困难。有限单元法和强度折减法的结合为解决这类问题提供了一条可行的途径,该方法求解这一类问题时不需要假定抗滑桩系统的结构类型以及岩土体作用在抗滑桩上的荷载分布形式,因此获得的计算结果往往更符合实际情况。本文以某隧道进口岸坡的П形抗滑桩为例,在详细分析岸坡稳定性状况以及可能失稳模式的基础上,采用有限元强度折减法对桩身内力、桩土作用关系以及治理效果进行计算和分析,结果表明利用有限元强度折减法求解这类具有复杂超静定结构抗滑桩系统的内力是一条有效的途径,计算结果对该治理工程的设计具有一定的应用参考意义。     

抗滑桩全桩计算的有限差分法

有限差分法求解分析抗滑桩(或预应力锚索抗滑桩)常采用Winkle弹性地基梁法。在对有限差分法分析的基础上,文章改进了有限差分法求解抗滑桩全桩内力和变形的计算方法。对全桩整体采用有限差分方法求解其内力,避免了将桩身分为受荷段和锚固段及利用滑面处连续条件求解的复杂过程。该方法对滑坡推力与地基系数取值,可依据实际情况假定为广义形式来模拟土体对抗滑桩的作用工况。能很好地求解预应力锚索作用于桩身或桩身设多道锚索工况下的锚索拉力值。并利用Excel对该方法编程,能便捷地计算抗滑桩内力。分析计算结果,当差分节点单元长度适当时,计算结果能够较好地满足工程应用。对比有限元计算分析结果,相同节点单元长度下两者计算相差较小,且节点长度越小,越接近有限元结果。     

返包式土工格栅加筋土高挡墙现场试验研究

为了研究返包式土工格栅加筋土高挡墙结构的受力、变形状态,分析其作用机理,进行了包括加筋土墙体基底应力、墙背侧向土压力、拉筋拉力和墙面水平变形等内容的现场试验,研究了加筋土墙体基底垂直应力、不同层位的拉筋拉力沿筋长的分布规律,加筋土挡墙潜在的破裂面位置,墙背侧向土压力沿墙高的分布规律以及墙面水平变形规律。测试结果表明,加筋土挡墙基底垂直土压力沿土工格栅拉筋长度方向呈非线性分布,最大值发生在拉筋中部附近,向拉筋两端方向逐渐减少;实测墙背侧向土压力沿墙高呈非线性形式分布,其值小于主动土压力;上部墙体拉筋应变沿筋长呈单峰值分布,下部墙体拉筋应变沿筋长呈双峰值分布;上部墙体潜在破裂面形状与“0.3H法”接近,而下部墙体潜在的破裂面形状与朗肯主动土压力理论接近;施工期墙面最大水平变形位置在墙高的下部,竣工后墙面最大水平变形发生在墙顶处等结论。     

一种基于强度折减法的次级滑动面分析方法研究

传统强度折减法通常对整个边坡区域进行折减,因而只能得到最小安全系数及对应的最危险临界滑动面。但在工程实践中,边坡的治理范围不能局限于最危险滑动面所包围的区域,不满足规范要求的次级滑动面范围内的坡体亦应得到治理。为了克服传统强度折减法在搜索次级滑动面方面的缺陷,提出了一种新的强度折减法。一般而言,边坡的屈服区域主要集中在滑动面两侧附近,形成一个沿滑动面分布的剪切破坏带,因此,在用强度折减法进行边坡稳定性分析时,只需对剪切破坏带范围内的局部坡体进行折减,则可得到该滑动面的安全系数和滑动面。基于以上认识,假设剪切破坏带沿对数螺旋曲线分布,不断调整对数螺旋曲线的形状以得到不同范围的剪切破坏带,对各种不同的剪切破坏带范围内的坡体进行局部折减计算,即可得到不同的安全系数及对应的滑动面。通过两个边坡算例（单台阶、双台阶）验证了该方法的可行性,并最后讨论了剪切破坏带的宽度对结果的影响。两个算例的计算结果表明,该方法不仅可以得到最危险临界滑动面,亦可得到任意安全系数对应的滑动面。     

北川地区擂鼓断裂在汶川地震中的地表破裂及其意义

2008年5月12日汶川MS8.0级特大地震发生后,在北川县擂鼓地区出露了一条擂鼓同震地表破裂带,该破裂带呈近南北向展布,位于映秀－北川断裂的中北段,其北东起于擂鼓镇柳林村北部,南西止于石岩村南部,以脆性破裂为特征,分别由3条呈北北东、北西西、北北东走向的地表破裂组成,延伸长度约4～5km,并切割了多种地貌单元,其平均垂直断距为1.5m,平均水平断距为1.4m,垂直与水平断距之比为1.07∶1。通过对该地表破裂带野外测量数据、几何展布结构及其成因机制的初步分析表明： 擂鼓断裂是出露于擂鼓地区的捩断层,具有捩断层的基本特性;其主要特征包括：1)擂鼓断裂的形成是由于在汶川地震中其东西两侧逆冲块体之间的差异性运动而引起; 2)断裂呈近南北向展布,与映秀－北川主干断裂近垂直相交;3)断面倾角较陡,为高角度断面的逆断层,具有逆冲兼走滑特征。     

Numerical Analysis of Deformation Control of Deep Foundation Pit in Ulanqab City

The deformation of a deep foundation pit in Ulanqab city was analyzed by in situ monitoring and numerical simulation data. First, a limit equilibrium method was used to evaluate the stability of foundation pit under non-supporting and ribbed plate anchored retaining wall support. The results show that the ribbed plate anchored retaining wall support could control the soil deformation effectively in deep foundation pit. Additionally, a three-dimensional finite difference model of excavation zone and reinforcement structure was established using FLAC3D software to simulate the excavation of deep foundation pit. The simulation results were verified by monitoring field data. The simulation results show that all the bolts are in the tension state, and the maximum axial force is located at the end of the bolt, with a value of 46.4 kN. Foundation pit excavation is a continuous stress release process, making the rock and soil of side wall of foundation pit slide towards the free surface. The largest deformation is 4.0 mm, near the side wall of foundation pit, consistent with the deformation monitoring results.

     

滑坡治理中抗滑桩桩位分析

采用有限元强度折减法对单排桩、桩间距为4 m的抗滑桩在不同设桩位置加固滑坡进行数值模拟,研究发现不同的桩位影响滑坡的稳定安全系数、滑动面的位置和形状。桩位不同决定了是桩后土体还是桩前土体滑落。桩位选择的一个重要原则是加固后滑坡体的稳定安全系数必须大于设计要求的安全系数,否则就会出现桩端越顶破坏或桩前土体滑落。斜坡在同一稳定安全系数下,桩位不同时抗滑桩所受的滑坡推力、桩的内力（剪力、弯矩）和桩的挠度不同。抗滑桩位于斜坡中部时,桩身长度较长,推力和内力大,是不合理的桩位;但桩身挠度大,提供的桩前抗力也大。按常规方法与按有限元法计算桩前无土体的悬臂桩桩上推力是相近的。把桩视作埋入滑坡体中的梁单元,按有限元法算出的推力比较小,因为它已经充分考虑了桩前土体的抗力。     

Folding of a detachment and fault – Modified detachment folding along a lateral ramp,southwestern Montana,USA

The inversion of the Middle Proterozoic Belt sedimentary basin during Late Cretaceous thrusting in Montana produced a large eastwardly-convex salient, the southern boundary of which is a 200 km-long oblique to lateral ramp subtended by a detachment between the Belt rocks and Archean basement. A 10 km-long lateral ramp segment exposes the upper levels of the detachment where hanging wall Belt rocks have moved out over the Paleozoic and Mesozoic section. The hanging wall structure consists of a train of high amplitude, faulted, asymmetrical detachment folds. Initial west-east shortening produced layer parallel shortening fabrics and dominantly strike slip faulting followed by symmetrical detachment folding. “Lock-up” of movement on the detachment surface produced regional simple shear and caused the detachment folds to become asymmetrical and faulted. Folding of the detachment surface after lock-up modified the easternmost detachment folds further into a southeast-verging, overturned fold pair with a ramp-related fault along the base of the stretched mutual limb.