汪儒鸿，周海清，彭国园

(1.陆军勤务学院军事设施系, 重庆 401331；2.西北核技术研究所, 陕西 西安 710024)

我国是一个地质灾害发生较为频繁的国家，据国家地质环境信息网统计，仅2013年全国就发生地质灾害15 403起，其中滑坡、崩塌、泥石流共14 703起，占地质灾害总数的95.4%，共造成669人死亡或失踪，直接经济损失102亿元。政府和社会为治理这些地质灾害投入了大量的工程经费，因此加强对边坡变形破坏机理与规律的研究，具有重要的社会经济效益和工程应用价值。在不同边坡的破坏形式中，有一类边坡滑坡破坏作用特别巨大。发生和进行得非常迅速，破坏前没有明显的变形迹象，缺乏灾变预兆。2010年9月1日发生在云南省保山市隆阳区瓦马乡河东村的一起滑坡灾害，灾前没有明显的变形迹象，因此没有引起坡下村民的警觉。滑坡启动后快速下滑，解体严重，破坏力巨大，最终导致48人死亡和失踪(图1)。

图1 云南省保山市隆阳区瓦马乡河东村滑坡灾害Fig.1 Landslide in Hedong Village Wama Township, Longyang District, Baoshan City, Yunnan Province

类似的滑坡并不鲜见，2015年11月13日浙江省丽水市莲都区里东村滑坡[1]、2013年1月11日云南省镇雄县高坡村滑坡[2]、以及2015年深圳光明新区“12·20”余泥渣土受纳场滑坡[3]都具有突发、快速、解体严重的特点。

普通的滑坡变形过程通常包括蠕动这一前期变形阶段，坡体在这一阶段滑动了一定的距离之后，通常滑动面的强度只剩下残余强度，较之蠕动变形之前强度已大大降低，即边坡的抗滑力大大减少，而下滑力则并没有减少，甚至由于蠕动变形的结果坡体还具有了一定的初速度，由此滑坡进入快速下滑的剧滑阶段。而此类滑坡则没有明显的蠕动变形阶段，或者这一阶段很短，滑坡发生后直接就进入剧滑阶段。参照典型滑坡剧滑前的变化，发生这种没有前期变形积累的突然快速下滑说明边坡失稳后其滑动面强度发生了快速的下降。根据这一特点和现有对滑坡的认识，一个合理的假设就是发生该类滑坡的边坡体存在较大的结构性强度[4]，使得边坡一旦失稳发生滑移后，其滑面土体的结构性被破坏，强度突然降低，导致边坡下滑的阻力大大减少，滑动大大加速，从而产生突然的缺少前期变形迹象的快速下滑。相关研究也从侧面验证了这一推理的合理性：林威[5]通过对花岗岩残积土微观结构及工程特性的研究，得出土体结构性的存在增强了土体抵抗外荷载破坏的能力，但一旦结构受到破坏边坡稳定性就会急剧下降的结论；刘邦安[6]通过对人工制备结构性土进行室内边坡实验发现，结构性土边坡的静力破坏具有均匀变形不明显，表现出一定脆性破坏的特征。综上，研究土体的结构性对边坡失稳突发性的影响，对边坡防灾减灾有重大意义。

1 数值分析方法

边坡滑坡变形是一个十分复杂的过程，其现象涉及深层土体滑动、地表植被破坏以及地面裂隙扩张等等。为方便统计滑坡过程中坡体的位移变形量，分析土体的结构性对边坡失稳变形规律的影响，采用离散元软件PFC2D来模拟不同土体强度的结构性土质边坡，具体方法如下。

1.1 数值模型搭建

根据研究需要搭建边坡二维数值模型，即搭建具有特定尺寸及边界形状的离散元颗粒集合。通过赋予离散元颗粒集合相应材料参数，使其力学特性与真实土体相似，从而使数值边坡整体的力学响应接近真实边坡。其中颗粒摩擦强度的大小主要由接触滑动模型参数(fric)控制；颗粒结构性强度的大小主要由平行黏结接触模型控制。平行黏结模型主要包括五个特征参数:法向和切向刚度(Pb-kn/Pb-ks)、法向和切向强度(Pb-nstren/Pb-sstren)、黏结半径(Pb-rad)，反映了接触颗粒在有限区域内填充胶合材料的本构特性。由于平行黏结接触模型的定义近似于结构性土二元介质模型[7]，也有较多文献[8-10]验证过离散元分析方法模拟结构性土的可行性，故不再展开讨论。

1.2 边坡失稳方法

以有限元、有限差分法为主的边坡稳定分析方法通常采用强度折减[11]来使边坡失稳，但这类方法并不适用于结构性土坡的滑坡模拟。因为结构性土具有应变软化特性，对该类土体采用人工强度折减不一定能反映其土体强度的真实非线性变化。更重要的是，结构性土目前尚无公认的有限元数值本构模型。因此，采取逐级增加坡顶荷载的方法使结构性土坡滑坡失稳。首先，使边坡在自重状态下计算至应力平衡，令坡体的位移形变归零；然后，在坡顶局部范围内对边坡逐级施加竖向荷载，每施加一级荷载，使边坡颗粒自主运动并重新计算至应力平衡。当施加至某一级荷载后，边坡颗粒无法再靠自身运动达到平衡，其最大位移无法收敛，边坡彻底失稳，称该级荷载为边坡最终荷载。在整个加载过程中，利用离散元软件的可视化技术，监测颗粒之间平行黏结接触来评估边坡土体结构性的完整性。

需要说明的是：逐级增加坡顶荷载使边坡失稳的方法属于特定数值分析方法，并不对应某类真实滑坡情况。但从破坏本质上来说，该过程与真实边坡在各种不利工况影响下滑面逐渐损伤的过程相对应，符合边坡渐进破坏[12-13]的定义，其结果具有一定分析价值。

1.3 滑坡变形度量

在边坡坡顶、坡中以及坡脚等不同位置选取表面颗粒作为滑坡监测点，提取各监测点的信息来度量边坡位移变形的大小。以某监测点为例：在坡顶荷载作用下，边坡变形导致该点颗粒产生相应位移，反映在二维数值平面上就是颗粒坐标的变化。由于岩土材料和边坡滑动方向的不确定性，在各级荷载作用下该点颗粒的位移应分解在x轴、y轴两个方向上来同时进行研究。通过统计该点颗粒在滑坡过程中实时x轴、y轴方向上的累积坐标增量，即可分析边坡在该点处的滑坡变形特征。

2 算例

2.1 基本模型搭建

在图2所示的框架内搭建结构性土坡的离散元模型。选取A-D四处表面颗粒作为滑坡监测点，并每隔一定距离给颗粒上色以便于观察坡体变形。

图2 PFC2D边坡模型Fig.2 PFC2D slope model

2.2 材料参数赋予

应用PFC2D搭建模型材料的力学性质主要定义颗粒密度、粒径分布、基本颗粒形状及颗粒基本细观属性，并通过控制颗粒的接触本构模型来模拟材料的本构模型[14]。材料宏观参数与微观参数的对应建立在大量的双轴压缩实验模拟[15]及参数调试上。该过程相对复杂，又并非本文重点，这里不做深入讨论，仅展示三种不同结构性强度土坡的颗粒主要微观参数及对应宏观参数(表1)。

以高结构性强度组及低结构性强度组边坡为例，两组边坡的颗粒集合在100 kPa围压作用下的偏应力-轴向应变曲线如图3所示。可以看出高结构性强度组颗粒的峰值强度更高，其残余强度与峰值强度之比较低结构性强度组颗粒更小。说明高结构性强度组颗粒的结构性较强，具有更加明显的应变软化现象，这也进一步验证了利用平行黏结接触模型模拟土体结构性的可行性。同理，将三组颗粒集合在不同围压下进行离散元双轴压缩测试，绘出其各自的应力莫尔圆，即可得到对应宏观力学参数。

表1 颗粒宏微观参数Table 1 Microscopic macro parameters

图3 偏应力-轴向应变曲线Fig.3 Deviatoric stress-axial strain curve

2.3 施加坡顶荷载

根据边坡的坡顶长度和颗粒数量把坡顶均布荷载换算成作用在每个颗粒的集中荷载施加在颗粒上。为保证所施加的为柔性荷载，坡顶的一排颗粒也用平行黏结模型黏结起来。根据土体参数设置，三组结构性土坡的坡顶初始荷载均设为0.25 kN/m，并同样以此大小逐级递增，在这个过程中每施加一级荷载使边坡重新计算至应力平衡，再继续施加下一级别荷载。当施加至某一级荷载之后，颗粒的最大位移无法收敛，边坡失稳，则将该级荷载定为最终荷载。

2.4 结果及分析

2.4.1最终荷载与边坡稳定性

数值模拟得出三组不同结构性强度土坡的最终荷载大小如表2所示。以高结构性强度土坡为例：图4表示离散元边坡颗粒的实时状态。由于坡顶荷载是逐级增加的，图4a为边坡在坡顶荷载为16 kN/m的条件下达到应力平衡时的状态。此时颗粒最大位移以及最大不平衡力都是收敛的，说明边坡还未失稳；图4b及图4c则是对边坡施加最终荷载以后，在坡顶荷载16.25 kN/m的条件下滑坡的动态发展过程。可以看出边坡此时已经失稳，坡体裂隙大量发育，颗粒最大位移也不收敛。

除直接观察以外，还可以结合颗粒位移矢量图来分析边坡稳定性。图5a为边坡在坡顶荷载16 kN/m的条件下达到应力平衡时的颗粒位移矢量图，此时边坡还未形成完整的滑面；而在最终荷载施加以后，颗粒不断滑动，逐渐形成了完整的滑面，边坡最终失稳，其动态发展过程如图5b及图5c所示。以上分析可知：各边坡在最终荷载施加之前，都能够通过颗粒的自主运动达到应力平衡，边坡尚未失稳。而一旦施加最终荷载以后，由于外界荷载超出了边坡所能承担的限度，导致边坡失稳，产生滑坡。

表2 各边坡最终荷载Table 2 Final load on the slope

2.4.2边坡失稳与土体结构性的损伤

在整个加载过程中，边坡土体结构性的完整情况也在不断发生变化，同样以高结构性强度土坡为例：图6中黑线表示颗粒之间的平行黏结接触。图6a为边坡施加坡顶荷载之后，于坡顶荷载为16 kN/m的条件下达到应力平衡时的颗粒平行黏结接触情况。可见颗粒之间的黏结接触已经出现损伤，特别是坡顶范围内有明显断裂现象，但从整个边坡来看并没有形成贯通的破坏面；图6b及图6c为边坡在施加最终荷载之后的颗粒平行黏结接触破坏情况。此时颗粒的平行黏结接触出现大量断裂破坏，并形成贯通的破坏面。

图4 离散元边坡颗粒Fig.4 Discrete element slope particles

图5 颗粒位移矢量图Fig.5 Particle displacement vector

图6 颗粒平行黏结接触示意图Fig.6 Particle parallel adhesion contact diagram

将其与图4和图5进行对比：土体结构性的损伤贯穿于结构性土坡整个滑坡过程，且同样动态反映出边坡的实时稳定性。在施加坡顶荷载以后，边坡土体的结构性处在不断损伤的情况下，但一直到最终荷载施加之前都没有形成贯通的破坏面，说明局部平行黏结接触的断裂对边坡整体稳定性的影响并不大，或者说此时潜在滑面土体的残余强度仍然满足承担外界荷载的要求；而在施加最终荷载以后，与边坡失稳相对应的就是颗粒平行黏结接触大量断裂，并形成贯通的破坏面。图6所示平行黏结接触的破坏贯通面与图5所示边坡滑面位置基本重合，说明潜在滑面的土体结构性受损最大，反之也说明由于土体结构性的损伤主要集中于边坡潜在滑面位置，导致该部分土体强度迅速下降，边坡沿该面加速下滑，最终失稳。因此边坡潜在滑面土体结构性的损伤是结构性土坡失稳的主要原因之一。

2.4.3失稳变形比较

潜在滑面土体结构性的损伤是导致结构性土坡失稳的主要原因之一，其破坏面的贯通是该类边坡失稳的重要标志。因此可通过统计各边坡在土体结构性的破坏面刚好贯通时的边坡位移变形大小，即各监测点的累积方向坐标增量，来评估不同结构性土坡的失稳突发性。

各边坡监测点的统计数据如表3所示。可以看出土体的初始结构性强度越高，边坡失稳前累积的滑动变形越小，失稳的突发性越大。为验证该统计结果的合理性，并直观展示不同结构性土坡在滑面损伤过程中累积位移形变的不同。这里选取高、低两种不同结构性强度土坡在各自潜在滑面土体的结构性完全破坏、形成贯通破坏面时的颗粒位移矢量图进行比较(图7)。可以看出高结构性强度土坡失稳时颗粒位移矢量较小，说明边坡滑动变形小，其失稳风险性较大。

表3 各监测点累积方向坐标增量Table 3 Cumulative direction coordinate increment of each monitoring point

“-”表示数据过小，不予统计。

图7 颗粒位移矢量图Fig.7 Particle displacement vector

3 结 论

以边坡渐进破坏理论为基础，对土体的结构性对边坡失稳突发性的影响进行了探讨。通过离散元数值模拟方法比较了不同结构性强度土坡的失稳变形规律，所得结论如下：

(1)土体结构性的损伤贯穿于结构性土坡滑坡的整个过程。该损伤集中于边坡潜在滑面的位置，由于滑面土体结构性的损伤导致该区域土体强度下降，进一步加剧了坡体的下滑，最终导致边坡失稳。

(2)结构性土坡在滑坡过程中会形成与滑面类似的土体结构性损伤贯通面，是结构性土坡失稳的重要标志之一。

(3)土体初始结构性强度越高的边坡，在边坡滑坡过程中累积的位移变形越小，失稳突发性越高。