Numerical simulation of earthquake-induced loess landslides based on particle flow method
-
摘要: 基于颗粒流理论研究土质边坡动力稳定性及其滑动过程是近年来滑坡研究的一个新热点。在野外调查和室内试验的基础上,通过标定土体细观参数、模型建立、动力输入、动态监测等过程,利用PFC2D程序模拟了西吉县兴平乡堡湾村下马达子滑坡的失稳破坏运动过程,得到了该滑坡的破坏运动机理。得到如下结论:① 下马达子滑坡的失稳机制是在地震作用下斜坡前缘牵引、后缘推挤,使得坡肩受拉发生破坏,失稳后坡肩位置较大的速度和位移是地震滑坡破坏力强、致灾范围大的主要原因;② 黄土地震滑坡的滑坡后壁相对平缓,这是区别于重力滑坡的重要特征之一;③ 颗粒流模拟得到的滑坡前后相对高差和长度与实际情况较为吻合,因此,颗粒流方法可以用于地震滑坡滑距的预测。Abstract: Research on the dynamic stability and sliding process of soil slopes based on particle flow theory is a new hot spot in landslide research in recent years. On the basis of field investigation and indoor experiment, the PFC2D program was used to simulate the instability failure movement process of the Xiamadazi landslide in Baowan village, Xingping township, Xiji county through the process of calibrating soil parameters, model establishment, power input, dynamic monitoring, etc., and the failure movement mechanism of the landslide is obtained. The following conclusions are obtained: ① The instability mechanism of the Xiamadazi landslide is that under the action of the earthquake, the leading edge of the slope is pulled and the trailing edge is pushed, causing the shoulder to be pulled and damaged. The larger velocity and displacement at the shoulder position after the instability is the main reason for the strong destructive force and the large disaster range of the earthquake landslide; ② The back wall of the earthquake-induced loess landslides is relatively flat, which is one of the important characteristics different from gravity landslides; ③ The relative altitude difference and length before and after the landslide obtained by particle flow simulation is more consistent with the actual situation. Therefore, the particle flow method can be used to predict the slip distance of earthquake landslides.
-
引言
受黄土动力易损性和南北地震带的双重影响,黄土地区发生地震滑坡的风险很高(徐张建等,2007)。历史研究表明,黄土地区一旦发生中强地震,往往诱发数量众多的滑坡,这些滑坡单体规模大、影响范围广、冲击速度快,对人民生命财产安全和城市建设的危害程度极高(常晁瑜等,2020)。实践证明,只有掌握地震滑坡发生机理,了解滑坡发生的过程和运动特征,才能采取有效的措施,避免或减轻地震滑坡的灾害(Zhou et al,2002)。
随着数值模拟方法的发展和计算机运算能力的提升,基于有限单元法、有限差分法和块体离散单元法等数值模拟方法,地震作用下斜坡的安全系数计算问题取得了极大进步(张江伟,2016;薄景山等,2019)。然而,强震作用下的滑坡是一个存在滑动、平移、转动的复杂过程,具有宏观上的不连续性和单个块体运动的随机性,当涉及到岩土体的大变形及岩土破坏问题时,上述数值模拟方法均具有一定的局限。而离散元颗粒流方法不受变形量限制(Hadjigeorgiou et al,2009;Tang et al,2009),可方便地处理非连续介质力学问题,有效地模拟岩土体的开裂、分离等非连续现象,对于探索地震滑坡的反映机理、运动过程和致灾范围具有一定优势(周健等,2000;曹文等,2017;石崇等,2018)。许多学者采用颗粒流离散元程序对边坡的滑移破坏机理开展了一系列的研究工作,例如:贺续文等(2011)讨论了节理连通率对边坡破坏形式的影响;周喻等(2016)从细观力学角度深入研究了顺层断续节理岩质边坡模型破坏过程的力学机制;李新坡和何思明(2010)对节理岩质边坡的破坏和运动过程进行了研究,分析了不同关键参数对破坏后的堆积形状和运动距离的影响;Scholtès和Donzé (2012)研究了含裂隙岩体的渐进性破坏。
通过以上分析,静力作用下滑坡运动失稳过程的研究进展较快,但针对地震滑坡失稳破坏运动过程的研究依然不成体系。鉴于PFC在模拟边坡失稳破坏方面的优势,本文拟在野外调查和室内试验的基础上,通过标定土体细观参数、模型建立、动力输入等过程,利用PFC2D程序模拟西吉县兴平乡堡湾村下马达子滑坡动力荷载下的失稳破坏运动过程,得到该地震滑坡的破坏运动机理,以期为黄土地区滑坡防治提供一定参考。
1. 颗粒流方法基本理论
1971年,伦敦大学帝国学院Cundall博士在分析准静力或动力条件下岩石边坡的运动时,借鉴分子动力学理论方法,首次提出了离散单元法(discrete element method,缩写为DEM)的概念(Cundall,Strack,1979)。PFC (particle flow code)程序又称颗粒流方法(王光谦,倪晋仁,1992),是基于离散单元模型框架,由计算机引擎和图形用户界面构成的细观分析软件,主要用于模拟有限尺寸颗粒的运动和相互作用。
PFC的计算循环使用时间-步长显示算法,需要对每个颗粒反复应用牛顿第二定律,并对每个颗粒间的接触单元反复应用力-位移定律。在接触模型中,颗粒接触点的接触力和相对位移可以分解为沿法线和切线方向的两个分量,通过力-位移定律可以由法向刚度和切向刚度把接触力的法向、切向分量分别与法向、切向相对位移联系起来。第i个接触点的力
$ {F_i} $ 可以沿法向和切向进行分解,分解为法向力$ {{\boldsymbol{F}}^{\rm{n}}} $ 和剪切力$ {{\boldsymbol{F}}^{\rm{s}}} $ ,即$$ {{\boldsymbol{F}}_i} = {{\boldsymbol{F}}^{\rm{n}}}{{\boldsymbol{n}}_i} + {{\boldsymbol{F}}^{\rm{s}}}{{\boldsymbol{t}}_i} \text{,} $$ (1) 式中,
${{\boldsymbol{n}}_i}$ ,${{\boldsymbol{t}}_i}$ 分别为接触平面的单位矢量。两个接触实体之间的法向“重叠”量表示为
${{\boldsymbol{U}}^{\rm{n}}}$ ,则法向力${{\boldsymbol{F}}^{\rm{n}}}$ 可以表示为$$ {{\boldsymbol{F}}^{\rm{n}}} = {K^{\boldsymbol{n}}}{{\boldsymbol{U}}^{\rm{n}}} \text{,} $$ (2) 式中,
${K^{\boldsymbol{n}}}$ 为接触点的法向刚度,${\boldsymbol{n}}$ 为接触面的单位向量。对于集合体来说,颗粒受到的剪切力与颗粒的相对运动、加载历史以及应力路径相关。因此,在颗粒流模型中,剪切力的计算以增量表示。切向力增量为
$$ \Delta {{\boldsymbol{F}}^{\rm{s}}} = - {k^{\rm{s}}}\Delta {{\boldsymbol{U}}^{\rm{s}}} \text{,} $$ (3) $$ {{\boldsymbol{F}}^{\rm{s}}} \leftarrow {{\boldsymbol{F}}^{\rm{s}}} + \Delta {{\boldsymbol{F}}^{\rm{s}}} {\text{≤}} \mu {{\boldsymbol{F}}^{\rm{n}}} \text{,} $$ (4) 式中,
$ {k^{\rm{s}}} $ 为接触剪切刚度,$ \mu $ 为颗粒摩擦系数,$\Delta {{\boldsymbol{U}}^{\rm{s}}}$ 表示在计算时步内接触位移的切向增量,由下式得到:$$ \Delta {{\boldsymbol{U}}^{\rm{s}}} = {v^{\rm{s}}}\Delta t \text{,} $$ (5) 式中,
$ \Delta t $ 为时间步长增量,${v^{\rm{s}}}$ 为速度。2. 参数标定
利用颗粒流进行数值模拟,本质上就是从其细观力学特征出发,将材料的力学响应问题从物理域映射到数学域内进行数值求解(Wang et al,2003)。土的结构特征除土颗粒的大小、形状、表面特性及粒度级配特征外,还包括颗粒间的排列与集合关系,孔隙的大小,颗粒间联结的特点等,这些土体细观结构的变化很大程度上影响到土体宏观力学特性的变化(周健,池永,2003;Härtl,Ooi,2008;陈达等,2018)。利用PFC2D进行数值模拟之前,需要先假定土体颗粒间的本构特性,再用颗粒流方法模拟双轴压缩或直剪等土工试验过程(Park,Song,2009),对其力学性质进行数值试验,以此得到细观力学参数与宏观力学参数之间的对应关系。
通过试验,得到黄土的土工试验参数列于表1。通过细观直剪数值试验模拟宏观大型直接剪切试验,首先赋予土体颗粒微观力学参数(表2),剪切盒模型尺寸为500 m×250 m,生成颗粒总数为1 706个(图1)。颗粒间的接触模型选用平行黏结模型,在表2的细观力学参数条件下,对直剪模型分别施加100,200,300和400 kPa的竖向应力,得到剪切应力-位移关系曲线,然后拟合出土体颗粒抗剪强度曲线(图2),从而得到数值模拟所用岩土体颗粒的黏聚力c为19.25 kPa,内摩擦角φ为13.5°。经过直接剪切数值试验后,得到的土体颗粒内摩擦角和黏聚力大小与表1中宏观参数基本对应,故此土体颗粒细观参数可用于数值模拟。
岩性 饱和密度/(g·cm−3) 孔隙比 含水率 饱和度 液限 塑限 压缩模量/MPa 黏聚力/kPa 内摩擦角/° 黄土 2.00 0.62 19.6% 86% 29.4% 18.5% 9.59 18.9 14 颗粒半径
/m颗粒密度
/(g·cm−3)阻尼比 摩擦系数 法向刚度
/MPa切向刚度
/MPa平行黏结参数 抗拉强度/Pa 黏聚力/Pa 摩擦角/° 摩擦系数 0.01—0.05 2 400 0.5 0.5 50 70 4×105 2×105 45 0.5 3. 模型建立
3.1 滑坡基本情况
1920年海原MS8.5特大地震诱发了数以百计的黄土滑坡,造成了惨重的人员伤亡和经济损失,由于特殊的地理位置和气候条件,尽管经历了近百年时间,这些黄土地震滑坡的外部形态多数依然保存完好,具有极高的学术研究价值(许冲等,2018;常晁瑜等,2019)。
本文作者参与了对海原特大地震诱发滑坡的再调查,实地调查1920年海原特大地震诱发黄土滑坡620个,卫星影像补充调查滑坡605个,总计1 225个。调查内容包括地理位置、滑动类型、地层岩性、坡形、斜坡结构类型、滑坡基本特征以及影响因素等。调查所得滑坡分布如图3所示。滑坡集中分布在三个区域:① 宁夏西吉与甘肃静宁、会宁三县交会地带,区内滑坡分布密集,发育特征明显,代表滑坡为西吉县兴平乡堡湾村下马达子滑坡;② 宁夏海原九彩乡境内,区内滑坡主要沿深切河流两侧发育,代表滑坡为海原县李俊滑坡;③ 固原市区东部河川乡山地区域,区内滑坡分布较为零散,代表滑坡为固原市石碑塬滑坡。
下马达子滑坡野外编号为XJ06007,位于西吉县兴平乡堡湾村喜家湾组葫芦河支流滥泥河左岸黄土斜坡上,地理位置为(105.621°E,35.893°N),是1920年海原地震诱发的大型黄土地震滑坡。下马达子滑坡整体形状呈长舌状,滑坡周界清晰,主滑方向为310°,滑坡体上陡下缓,上部平均坡度约15°,下部平均坡度小于10°,长度达到1 280多米。由于上部土体下滑,导致上部较薄,厚度只有8—13 m,而下部滑体较厚,约15—25 m;后缘较陡,坡度约20°,该段长约190 m,滑体上发育大量落水洞。中部向下到居民区长约640 m,该段坡体平缓,滑体已被改造为梯田。下部滑坡堆积区长约450 m。滑坡上部宽约 140 m,中部较宽,最宽约270 m。滑坡总体积约360万 m3。滑坡区相对高差为220 m (图4)。
根据野外勘探(图5),下马达子滑坡地层岩性主要为上更新统马兰黄土、中更新统离石黄土和新近系砖红色泥岩,滑坡发育在黄土与泥岩的不整合接触上。由于黄土垂直节理极其发育,且其结构松散,当土体含水量大时,其抗拉、抗剪强度大大降低,其它条件具备时极其易发生变形。而且此处坡体为阴坡,受北部西伯利亚冷空气影响,黄土堆积厚度大,这也是影响坡体失稳的一个重要因素。
3.2 模型建立
基于下马达子滑坡坡面图(图5),根据等体积法进行地形复原,所得滑坡发生前斜坡剖面图示于图6,而后通过颗粒流软件PFC2D进行建模,对泥岩颗粒施加重力后,颗粒做自由落体运动,循环一定次数后形成矩形,而后进行削坡操作,利用二次函数将泥岩层削出凹面,再对斜坡开挖形成基本泥岩层形状(图7),此时模型长约1 040 m,高度约200 m。对颗粒间不平衡力进行设置,多次运算后滑床颗粒静止平衡,此时锁定泥岩层颗粒速度及位移。接着采用落球法生成黄土颗粒,采用平行黏结模型,分三次下落到泥岩层指定位置,每次下落后黄土颗粒进行一次自重平衡,通过不断运算调整颗粒间不平衡力大小,直至不平衡力下降至允许范围,多次试算平衡后生成滑坡发生前基本模型,如图7所示。
3.3 动力输入
采用拟静力法对地震进行输入,拟静力法(刘红帅等,2007;Zhou et al,2015)是一种用静力学方法近似解决动力学问题的简易方法,它发展较早,迄今仍然被广泛使用。其基本思想是在静力计算的基础上,将地震作用简化为一个惯性力系附加在研究对象上,其核心是设计地震加速度的确定。由现场勘察及后期地震灾害分析报告可知,该滑坡所在地在1920年海原地震中的烈度为Ⅸ,故此次模拟采用Ⅸ度区的设计加速度即0.4g施加于模型上,使颗粒在水平方向上产生相应的惯性力。
4. 结果分析
为了获得在外力输入下滑坡的运动情况,在关键位置设置6个测点(图8),以便于观察模拟中滑坡的发生机理和运动特征。
4.1 滑坡过程
滑坡模型运行60 s后停止,为了解滑坡的发生过程,截取不同时间下滑坡的位移图,示于图9。从图中可以清晰地看出:外力输入后,模型开始失稳,在20 s时,边坡坡度发生改变,上部滑体在地震波作用下产生明显位移,并挤压坡面颗粒,推动坡面颗粒向前运动,滑坡前缘颗粒在地震波作用下逐渐脱离滑床,开始运动;30 s时,坡面颗粒在上部滑体的重力和水平力的双重作用下,抗滑能力进一步降低,此时已可辨别出滑动面,滑面颗粒沿滑面不断向下运动,同时牵引后缘颗粒下滑;40 s时,随着滑坡后缘颗粒的逐渐下滑,滑坡后缘已经破坏,后缘颗粒已全部沿滑动面下滑,滑面颗粒在后缘颗粒的推移和前缘颗粒的牵引的共同作用下,速度进一步增大,位移增长迅速;50 s时,前缘颗粒在地面摩擦力作用下逐渐静止,在坡角处产生堆积,前缘颗粒位移达到40 m左右,此时坡面抗滑力逐渐增大,滑体与滑床接连处颗粒速度减小,逐渐趋于静止;60 s时,滑面颗粒运动缓慢且处于减速阶段,逐渐趋于稳定,此时滑坡后缘颗粒在滑面位移达到最大值100 m左右,边坡再次趋于稳定。
整体上,滑坡的发生过程中,滑坡后缘颗粒的推挤作用和前缘颗粒的牵引作用均发挥了重要作用,滑体在推挤和牵引的双重作用下,逐渐发生失稳破坏,形成贯通滑动面,与滑床间摩擦力降低,直至发生滑动破坏。
4.2 速度变化
为了监测滑床内部与滑体的速度特征,检测了6个测点颗粒的速度变化,得到速度曲线示于图10。由图10可以看出:在整个滑坡过程中,滑床内部颗粒产生的速度变化较小,趋近于零;而滑坡后缘测点(测点1,2,3)速度变化较为剧烈。说明在地震波输入后滑坡后缘运动剧烈程度大于滑坡前缘(测点5),这与滑坡后缘失稳后重力势能转化为动能有关,测点3处速度高于其它五个测点,表明其首先被破坏并持续滑动。
4.3 位移变化
通过对比滑坡发生前后的变化情况(图11),可以看出模型中颗粒发生了较大的位移(Behbahani et al,2013)。为详细研究不同测点位移情况,通过对五个测点(测点1—5)的水平位移和垂直位移进行监测,得到如图12所示的位移曲线。通过对比可知,测点位置不同,水平和垂直方向的位移大小关系也不同,但大体上呈现出水平位移约为垂直位移的1—2倍。测点3处颗粒在水平和竖直方向上位移均最大,可见此处颗粒运动距离最大。测点5在运行50 s左右达到垂直位移最大值,但与其它测点不同的是,该测点处颗粒垂直位移未进一步明显发展,即该处颗粒在竖直方向不再发生明显运动,这是由于滑体在坡脚的堆积挤压作用,在竖直方向上,使得该测点颗粒克服重力,达到平衡状态,不再产生加速度,从而出现位移无明显变化的现象。在水平方向上,通过滑坡过程位移图可知,滑坡后缘处(测点1,2,3)水平位移基本较滑坡前缘处(测点5)大,这是由于滑坡前缘颗粒在破坏后期脱离滑面后产生堆积,受到地面摩擦力作用趋于稳定;而测点1处颗粒位移较小则是因为此处颗粒靠近滑坡后缘,滑坡后缘颗粒经滑面滑下后产生堆积再次趋于稳定,滑体与滑面间抗滑力逐渐恢复,此处颗粒未运动至坡脚。
4.4 滑动面情况
结合运动终止时刻的颗粒位移图像,可以确定下马达子滑坡的滑动面如图13所示,这与现场勘察结果显示滑坡体上陡下缓(图5)基本吻合,表明颗粒流方法可以用于地震滑坡滑动面的确定。
5. 讨论与结论
本文通过标定土体细观参数、模型建立、动力输入等过程,利用PFC2D程序模拟了西吉县兴平乡堡湾村下马达子滑坡的失稳破坏运动过程,得到如下结论:
1) 下马达子滑坡的失稳机制是地震作用下的斜坡前缘牵引、后缘推挤、坡肩受拉破坏,由于坡肩高度高且具有临空面,失稳后坡肩位置速度和位移均较大,是地震滑坡破坏力强、致灾范围大的主要原因;
2) 由于黄土特有的垂直节理,黄土高原区常常存在高陡的斜坡,且重力滑坡通常具有高陡的后壁,而地震作用下滑坡的后壁相对平缓,这是地震滑坡区别于重力滑坡的重要特征之一;
3) 通过颗粒流模拟得到滑坡前后相对高差为250 m左右,滑坡后总长度达到1 080 m,滑坡后坡角约13.0°;与实际情况相对高差220 m,滑坡总长度1 040 m,滑坡后坡角约11.94°相比较为符合,颗粒流方法可以用于地震滑坡滑距的预测。
值得指出,拟静力法是将地震期间最大惯性力施加在土体上,认为土体中各点的最大加速度同时出现,且没有考虑土体随时间变化的非线性,定量分析还不够精准。今后的研究中,需要利用动力时程分析法考虑地震动的时间过程对滑坡失稳破坏过程的影响,计算滑坡的演化过程。
-
表 1 土工试验参数
Table 1 Geotechnical test parameters
岩性 饱和密度/(g·cm−3) 孔隙比 含水率 饱和度 液限 塑限 压缩模量/MPa 黏聚力/kPa 内摩擦角/° 黄土 2.00 0.62 19.6% 86% 29.4% 18.5% 9.59 18.9 14 表 2 细观试验参数
Table 2 Micro-scale test parameters
颗粒半径
/m颗粒密度
/(g·cm−3)阻尼比 摩擦系数 法向刚度
/MPa切向刚度
/MPa平行黏结参数 抗拉强度/Pa 黏聚力/Pa 摩擦角/° 摩擦系数 0.01—0.05 2 400 0.5 0.5 50 70 4×105 2×105 45 0.5 -
薄景山,段玉石,常晁瑜,李孝波,彭达,闫东晗. 2019. 斜坡地震稳定性研究的若干问题[J]. 自然灾害学报,28(1):1–8. Bo J S,Duan Y S,Chang C Y,Li X B,Peng D,Yan D H. 2019. Some problems of study on slope stability under earthquake[J]. Journal of Natural Disasters,28(1):1–8 (in Chinese).
曹文,李维朝,唐斌,邓刚,李俊峰. 2017. PFC滑坡模拟二、三维建模方法研究[J]. 工程地质学报,25(2):455–462. Cao W,Li W C,Tang B,Deng G,Li J F. 2017. PFC study on building of 2D and 3D landslide models[J]. Journal of Engineering Geology,25(2):455–462 (in Chinese).
常晁瑜,杨顺,焦淙湃,彭达. 2019. 1920年海原特大地震诱发黄土滑坡滑距统计与预测分析[J]. 防灾科技学院学报,21(2):36–43. doi: 10.3969/j.issn.1673-8047.2019.02.006 Chang C Y,Yang S,Jiao C P,Peng D. 2019. Sliding distance statistics and prediction of the loess landslide triggered by the 1920 Haiyuan earthquake[J]. Journal of Institute of Disaster Prevention,21(2):36–43 (in Chinese).
常晁瑜,薄景山,李孝波,乔峰,闫东晗. 2020. 地震黄土滑坡滑距预测的BP神经网络模型[J]. 地震工程学报,42(6):1609–1614. doi: 10.3969/j.issn.1000-0844.2020.06.1609 Chang C Y,Bo J S,Li X B,Qiao F,Yan D H. 2020. A BP neural network model for forecasting sliding distance of seismic loess landslides[J]. China Earthquake Engineering Journal,42(6):1609–1614 (in Chinese).
陈达,薛喜成,魏江波. 2018. 基于PFC2D的刘涧滑坡破坏运动过程模拟[J]. 煤田地质与勘探,46(4):115–121. doi: 10.3969/j.issn.1001-1986.2018.04.019 Chen D,Xue X C,Wei J B. 2018. Simulation of failure process of Liujian landslide based on PFC2D[J]. Coal Geology and Exploration,46(4):115–121 (in Chinese).
贺续文,刘忠,廖彪,王翠翠. 2011. 基于离散元法的节理岩体边坡稳定性分析[J]. 岩土力学,32(7):2199–2204. doi: 10.3969/j.issn.1000-7598.2011.07.046 He X W,Liu Z,Liao B,Wang C C. 2011. Stability analysis of jointed rock slopes based on discrete element method[J]. Rock and Soil Mechanics,32(7):2199–2204 (in Chinese).
李新坡,何思明. 2010. 节理岩质边坡破坏过程的PFC2D数值模拟分析[J]. 四川大学学报(工程科学版),42(增刊):70–75. Li X P,He S M. 2010. Numerical analysis of the failure of heavily jointed rock slopes using PFC2D[J]. Journal of Sichuan University (Engineering Science Edition)
,42(S1):70–75 (in Chinese). 刘红帅,薄景山,刘德东. 2007. 岩土边坡地震稳定性评价方法研究进展[J]. 防灾科技学院学报,9(3):20–27. doi: 10.3969/j.issn.1673-8047.2007.03.006 Liu H S,Bo J S,Liu D D. 2007. Development on study of seismic stability evaluation methods of rock-soil slopes[J]. Journal of Institute of Disaster-Prevention Science and Technology,9(3):20–27 (in Chinese).
石崇,张强,王盛年. 2018. 颗粒流(PFC5.0)数值模拟技术及应用[J]. 岩土力学,39(增刊):36. Shi C,Zhang Q,Wang S N. 2018. Numerical simulation technology and application of particle flow (PFC5.0)[J]. Rock and Soil Mechanics,39(S2):36 (in Chinese).
王光谦,倪晋仁. 1992. 颗粒流研究评述[J]. 力学与实践,14(1):7–19. Wang G Q,Ni J R. 1992. A review of research on particle flow[J]. Mechanics and Practice,14(1):7–19 (in Chinese).
许冲,田颖颖,马思远,徐锡伟,周本刚,吴熙彦,庄建琦,高玉欣,吴玮莹,黄学强. 2018. 1920年海原8.5级地震高烈度区滑坡编录与分布规律[J]. 工程地质学报,26(5):1188–1195. Xu C,Tian Y Y,Ma S Y,Xu X W,Zhou B G,Wu X Y,Zhuang J Q,Gao Y X,Wu W Y,Huang X Q. 2018. Inventory and spatial distribution of landslides in Ⅸ-Ⅺ high intensity areas of 1920 Haiyuan (China) M8.5 earthquake[J]. Journal of Engineering Geology,26(5):1188–1195 (in Chinese).
徐张建,林在贯,张茂省. 2007. 中国黄土与黄土滑坡[J]. 岩石力学与工程学报,26(7):1297–1312. doi: 10.3321/j.issn:1000-6915.2007.07.001 Xu Z J,Lin Z G,Zhang M S. 2007. Loess in China and Loess landslides[J]. Chinese Journal of Rock Mechanics and Engineering,26(7):1297–1312 (in Chinese).
张江伟. 2016. 土质边坡地震动力响应特性分析及稳定性评价[J]. 国际地震动态,(1):44–45. Zhang J W. 2016. Seismic dynamic response characteristics and stability evaluation of soil slopes[J]. Recent Developments in International Seismology,(1):44–45 (in Chinese).
周健,池永. 2003. 砂土力学性质的细观模拟[J]. 岩土力学,24(6):901–906. doi: 10.3969/j.issn.1000-7598.2003.06.006 Zhou J,Chi Y. 2003. Mesomechanical simulation of sand mechanical properties[J]. Rock and Soil Mechanics,24(6):901–906 (in Chinese).
周健,池永,池毓蔚,徐建平. 2000. 颗粒流方法及PFC2D程序[J]. 岩土力学,21(3):271–274. doi: 10.3969/j.issn.1000-7598.2000.03.020 Zhou J,Chi Y,Chi Y W,Xu J P. 2000. The method of particle flow and PFC2D Code[J]. Rock and Soil Mechanics,21(3):271–274 (in Chinese).
周喻,韩光,吴顺川,胡乃联. 2016. 断续节理岩体及岩质边坡破坏的细观机制[J]. 岩石力学与工程学报,35(增刊):3878–3889. Zhou Y,Han G,Wu S C,Hu N L. 2016. Meso mechanical failure mechanism of rock mass and slope with intermittent joints[J]. Chinese Journal of Rock Mechanics and Engineering,35(S2):3878–3889 (in Chinese).
Behbahani S S,Moarefvand P,Ahangari K,Goshtasbi K. 2013. Pomiary przemieszczeń i sił kontaktu pomiędzy cząstkami materialnymi w trakcie wybierania wyrobiska pochyłego przy pomocy programu PFC2D[J]. Arch Min Sci,58(2):495–504.
Cundall P A,Strack O D L. 1979. A discrete numerical model for granular assemblies[J]. Geotechnique,29(1):47–65.
Hadjigeorgiou J,Esmaieli K,Grenon M. 2009. Stability analysis of vertical excavations in hard rock by integrating a fracture system into a PFC model[J]. Tunn Undergr Space Technol,24(3):296–308. doi: 10.1016/j.tust.2008.10.002
Härtl J,Ooi J Y. 2008. Experiments and simulations of direct shear tests:Porosity,contact friction and bulk friction[J]. Granular Matter,10(4):263–271. doi: 10.1007/s10035-008-0085-3
Park J W,Song J J. 2009. Numerical simulation of a direct shear test on a rock joint using a bonded-particle model[J]. Int J Rock Mech Min Sci,46(8):1315–1328. doi: 10.1016/j.ijrmms.2009.03.007
Scholtès L,Donzé F V. 2012. Modelling progressive failure in fractured rock masses using a 3D discrete element method[J]. Int J Rock Mech Min Sci,52:18–30. doi: 10.1016/j.ijrmms.2012.02.009
Tang C L,Hu J C,Lin M L,Angelier J,Lu C Y,Chan Y C,Chu H T. 2009. The Tsaoling landslide triggered by the Chi-Chi earthquake,Taiwan:Insights from a discrete element simulation[J]. Eng Geol,106(1/2):1–19.
Wang C,Tannant D D,Lilly P A. 2003. Numerical analysis of the stability of heavily jointed rock slopes using PFC2D[J]. Int J Rock Mech Min Sci,40(3):415–424. doi: 10.1016/S1365-1609(03)00004-2
Zhou J X,Zhu C Y,Zheng J M,Wang X H,Liu Z H. 2002. Landslide disaster in the loess area of China[J]. J Forestry Res,13(2):157–161. doi: 10.1007/BF02857244
Zhou X P,Qian Q H,Cheng H,Zhang H P. 2015. Stability analysis of two-dimensional landslides subjected to seismic loads[J]. Acta Mech Solida Sin,28(3):262–276. doi: 10.1016/S0894-9166(15)30013-6
-
期刊类型引用(3)
1. 付亚伟. 地震和降雨共同作用下植生混凝土生态护坡稳定性分析. 长沙理工大学学报(自然科学版). 2024(01): 100-109 . 百度学术 2. 杜国梁,刘俊江,丛璐,谷丽莹,乐婧瑜,王闫超. 数值模拟交通荷载作用下黄土边坡动力响应研究. 中国地质灾害与防治学报. 2024(05): 59-66 . 百度学术 3. 常晁瑜,徐久欢,薄景山,杨济源,孙雪晨,顾骏. 基于颗粒流的地震液化型滑坡运动学特征分析. 地震工程与工程振动. 2022(06): 153-161 . 百度学术 其他类型引用(3)