基于颗粒流方法的堆积层滑坡运动过程模拟

以颗粒流离散元为研究方法对勉县杨家湾十组堆积层滑坡破坏方式与运动过程进行数值模拟研究。通过PFC^2D双轴模拟试验所标定的岩土体宏观模拟参数与室内试验所获取的宏观实测参数进行对比,确定堆积层滑坡所需的颗粒细观参数,然后将标定的细观参数代入堆积层滑坡模型,对滑坡破坏方式及运动过程进行模拟研究。结果表明：滑坡破坏在初始阶段蠕滑变形累积,滑坡体挤压坡脚,直至坡脚产生剪切破坏,并向上牵引发展,使得滑坡整体顺接触面破坏下滑并堆积于坡脚,表现为典型的牵引式渐进破坏,结果与实际情况基本吻合。研究认为,采用颗粒流方法对堆积层滑坡破坏与运动过程的模拟研究具有较高的适用性,对该类滑坡防治具有一定的参考意义。     

基于颗粒流模拟的堆积层滑坡过程分析

基于颗粒流离散元程序模拟了堆积层滑坡变形破坏及滑动的全过程。结合现场地质调查资料,总结了拜殿乡Ⅰ号滑坡的成因机制及影响因素,分析表明该滑坡是在内外因共同作用下的结果,其主要原因是滑体物质松散,雨水入渗形成软弱滑动带,导致坡体的变形滑动。采用数值试验确定土体细观参数与宏观力学性质的关系,建立拜殿乡Ⅰ号滑坡数值模型,考虑降雨对滑带土强度的软作用化,对滑坡变形破坏及滑动全过程进行了模拟。模拟结果表明：滑带土强度降低后坡脚出现剪切破坏集中区,后缘处出现拉张裂破坏集中区,随着坡体缓慢的蠕动变形,坡脚的裂隙范围逐渐扩大,坡体中后缘出现多条拉张裂隙;根据滑体的运动速度将滑动过程分为蠕动变形、滑动加速、滑动减速和缓慢堆积4个阶段;滑坡滑动情况和最终堆积状态表明拜殿乡Ⅰ号滑坡严重威胁龙潭村小学及乡村公路的安全,模拟结果对该滑坡的工程治理具有重要的指导价值。     

基于PFC2D的刘涧滑坡破坏运动过程模拟

以陕西省洛南县刘涧滑坡为研究对象,采用颗粒流离散元法对其破坏运动过程进行数值模拟。首先通过双轴数值试验对滑坡饱和土体进行细观参数标定,并与室内试验中饱和土体宏观力学参数进行对比,经验证该细观参数能应用到滑坡的破坏运动分析中,进而引入颗粒流（PFC2D）程序中平行黏结模型,采用ball-wall建模方法建立滑坡模型,对滑坡不同关键部位颗粒进行位移、速度监测,阐明其破坏运动特征。模拟结果表明,降雨为刘涧滑坡的直接诱发因素,斜坡变形破坏模式为由坡脚开挖引起的自前缘向后部牵引-孔隙水压力诱发的后部向前缘推移式滑塌。总体特征为上部推移,中部剪切,下部牵引;滑坡滑动最高时速13.4 m/s,最大滑移170 m,滑动阶段持续25 s。利用颗粒流法对滑坡的破坏运动过程模拟具有较好的适用性,可为工程决策提供依据。      

基于PFC3D的鱼鳅坡滑坡运动过程分析

以贵州省开阳县鱼鳅坡滑坡为研究对象,采用颗粒流离散元（PFC3D）对其破坏运动过程进行数值模拟。采用Ball-Wall建模方法建立滑坡模型,对滑坡不同关键部位颗粒进行位移、速度监测,阐明其破坏运动特征。结果表明,降雨为鱼鳅坡滑坡的直接诱发因素。该滑坡在破坏初始阶段以蠕滑变形为主,随着变形量的增加,滑坡体不断挤压坡脚,滑坡岩土体到达应力平衡极限,坡脚产生剪切破坏,并向上牵引发展,滑坡发生整体滑动,斜坡变形破坏模式为蠕滑-拉裂,按照力学条件为牵引式破坏。滑坡滑动最高时速12.4 m/s,最大滑移80 m,滑动阶段持续50 s。研究成果可为对该类滑坡影响范围预测,以及工程措施的制定具有一定的参考意义。     

不同墙体位移方式下被动土压力的颗粒流模拟

用二维颗粒流程序（PFC2D）对各种移动模式的刚性挡墙被动土压力进行了模拟。再根据已有的模型试验建立颗粒流模型,分别让刚性墙体朝土体一侧平移、绕墙体顶部转动和绕墙体底部转动,给出了这三种情况下墙后被动土压力的分布图形以及相应的竖向应力的分布图形,同时给出了墙后被动土压力随墙体位移的变化规律。最后,将计算结果与实验结果进行了对比表明,颗粒流数值模拟方法不仅使用方便,而且能较好地模拟墙后被动土压力的分布规律。     

金沙江白格滑坡运动过程特征数值模拟与危险性预测研究

2018年10月、11月于金沙江川藏交界处江达县波罗乡白格村先后发生两次体积约2400×104 m3和850×104 m3的滑坡,两次滑坡平均运动距离1400 m,堵塞金沙江形成堰塞湖。首次形成的堵江滑坡坝天然溃决,未造成人员伤亡;然而第2次滑坡堵塞第1次滑坡自然溃口,导致堰塞湖库容迅速增加到3.85×108 m3。政府部门立即开展抢险工作,通过人工修建溢洪道的方法成功泄洪,极大程度上降低洪水风险。本文利用PFC3D颗粒流软件模拟两次滑坡的发生、运动、堆积过程,并在反演结果的基础上对白格滑坡滑源区残留潜在不稳定部分未来失稳的运动路径和堆积范围进行预测,对其危险性进行科学评价。结果表明:（1）滑坡在重力作用下失稳,除了受初始势能的影响外,微地貌也是决定滑坡运动路径与距离的关键因素之一;（2）PFC3D颗粒流数值模拟方法适用于类似于白格滑坡这类碎屑流类型的滑坡,两次滑坡反演得到的堆积厚度、堆积范围均与真实结果相近;（3）利用两次事件反演所得参数,可以预测若滑源区潜在不稳定部分同时失稳,则形成约70 m高的滑坡坝,可能再次堵塞金沙江。     

颗粒流模拟方法在岩土工程中的应用

归纳总结了颗粒流模拟方法与其它数值模拟方法的异同之处，介绍了颗料流模拟方法解题的步骤及其在岩土工程领域的应用实例。     

土坡稳定分析的颗粒流模拟

边坡的失稳破坏运动是一个存在岩土体的滑动、平移、转动的复杂过程,具有宏观上的不连续性和单个块体运动的随机性。采用颗粒流模拟土坡的变形破坏全过程,不需要假定滑移面的位置和形状,颗粒根据所受到的接触力调整其位置,最终从抗剪强度最弱的面发生剪切破坏,因此离散元法是模拟边坡变形破坏力学行为的比较理想的途径。运用PFC（颗粒流）程序对砂性土坡和黏性土坡分别进行数值模拟,分析细观参数对土坡破坏型式的影响。结果显示,土性对边坡的破坏型式有很大的影响,随着颗粒黏性的增大,边坡破坏类型从塑性破坏向脆性破坏过渡。     

颗粒流方法及PFC2D程序

归纳总结了颗粒流模拟方程产生的背景,比较了与其他模拟方法的异同之处,介绍了颗粒流模拟方法解题的步骤及其在岩土工程领域的应用实例。     

青藏高原察达高速远程滑坡运动过程与形成机理

为了了解青藏高原察达高速远程滑坡的运动过程与形成机理,运用遥感测绘、无人机地形测绘和现场勘查资料对滑坡进行分区,对滑坡形成机理进行研究,并利用PFC2D数值模拟对地震工况下滑坡运动过程进行模拟.将察达高速远程滑坡分为源区,流通区和堆积区;数值模拟结果得到滑坡平均运动速度为15~20 m/s,运动时间150 s,最大运动距离为2 800 m.察达滑坡为地震条件下诱发的高速远程滑坡,源区砾岩对上部堆积体后缘铲刮推移,使得上部堆积体产生整体变形,其运动过程可分为崩滑→铲刮→滑移→堆积4个阶段.      

考虑破碎的堆石料颗粒流数值模拟

以双江口300 m级高土石坝为工程背景,根据有关资料的试验数据,利用二维颗粒流计算机仿真技术对堆石料在一定围压下的颗粒破碎情况进行了数值模拟,数值试样由2 574个簇颗粒单元组成。数值模拟与室内平面应变试验结果对比表明,用选定的细观参数进行数值模拟与平面应变试验结果基本一致,说明用颗粒流方法进行数值模拟能够较好的再现堆石料的颗粒破碎特性。与室内试验相比,数值仿真模拟还能得到颗粒破碎的细观演化规律和破碎带的分布范围,从而为进一步建立考虑颗粒破碎的堆石坝本构模型,研究大坝填筑、水库蓄水过程以及地震情况下的颗粒破碎过程以及颗粒破碎所引起的坝体变形提供了一条途径。     

泥石流启动过程PFC数值模拟

泥石流松散碎屑物质具有散粒体的基本特征,在暴雨激发下容易形成泥石流,整个过程具有散粒体大变形的特征。颗粒流理论是基于离散单元法模拟圆形颗粒介质的运动及其相互作用,在模拟颗粒相互作用和大变形问题研究方面具有显著的优越性。为分析泥石流松散碎屑物质启动形成泥石流的过程及其与土体含水率的关系,采用二维颗粒流程序（PFC2D）分析降雨作用下松散碎屑物质启动形成泥石流的过程,揭示崩滑堆积体在降雨作用下含水率超过临界值后质点运动速度和位移增加,松散碎屑物质启动并加速而导致滑坡泥石流连锁式破坏的过程和机制。     

砂土中单桩静载室内模型试验及颗粒流数值模拟

通过自行设计的可视化模型箱,进行了单桩静载室内模型试验,考虑了不同桩径、不同土体密实度等影响因素,研究了单桩的沉降模式、桩侧摩阻力和桩端阻力的发挥性状及随沉降的发展模式,并通过拍摄跟踪和图像处理技术,研究了桩端和桩周土体的孔隙率的变化规律;对二维颗粒流程序进行开发,模拟了单桩静载过程中桩端阻力和桩周土体孔隙率的变化情况。结果表明,颗粒流数值模拟能够较好地模拟单桩静载过程,研究成果进一步揭示了单桩静载过程的宏细观机制。     

堆石料单轴流变试验的颗粒流模拟

岩石断裂力学的亚临界裂缝扩展理论认为微裂缝扩展可导致岩石破碎,即岩石颗粒破碎具有时间效应。根据亚临界裂缝扩展理论,提出了考虑微裂缝扩展导致堆石颗粒破碎时间效应的数值流变模拟新方法,并进行了考虑颗粒不同典型破碎模式的单轴流变颗粒流数值试验。在对比数值与室内流变试验曲线的基础上,分析了数值流变过程中颗粒破碎情况与颗粒体内部结构发展过程等。研究成果表明,两种试验手段得到的堆石流变的发展趋势基本一致,由微裂缝扩展引起的颗粒延时破碎是堆石流变的主要原因之一,深化了对堆石料变形机制的认识。     

千将坪滑坡动力学过程仿真模拟

滑坡体的失稳破坏是一个动态过程,动力学行为起着非常重要的作用。采用DDA方法对长江三峡库区千将坪滑坡的运动全过程进行了数值模拟研究。模拟方案充分依据该滑坡的地质、地形特征,按不同岩土体和地质结构面类型进行块体单元划分,模拟了滑坡发生、发展的渐进破坏过程以及滑坡触发后的运动情况。模拟结果表明,千将坪滑坡是以斜坡坡脚的局部破坏为其运动的开始阶段,并进一步牵引上部滑体,在地下水压力作用下最终产生整体滑动。     

石灰岩细观力学特性的颗粒流模拟

岩体地区地质灾害的发生和发展取决于岩石细观组分的运动学行为。研究岩石运动学行为时通常将岩石作为整体研究对象较多,而直接以细观组分为对象的研究较少。以石灰岩为例,根据室内试验获得的岩石力学性质指标,使用基于非连续介质理论的颗粒流方法,将材料离散成刚性颗粒组成的模型,把颗粒细观变化与宏观力学特性联系起来,建立了石灰岩的细观结构模型,获得了颗粒接触力、颗粒接触模量、接触连接强度和连接刚度比等细观力学参数。由于文中直接以细观成分为研究对象、反映了岩石和岩体组成的本质特点,所得结论不仅对含裂隙岩石本构关系研究具有广阔的应用前景,而且对岩体工程性质和地质灾害机制研究也具有重要的理论意义。     

石笼挡土墙的颗粒离散元细观力学模拟研究

针对在施工过程中,石笼挡土墙墙面位移的发展、墙背填土中土应力的变化的问题,从细观层面,采用基于离散单元法的颗粒流程序,通过多组试件的模拟压缩试验,确定了填土、填石、石笼等材料合理的力学参数值,以工程中的实体石笼挡土墙为原型,建立了墙背阶梯式石笼挡土墙数值模拟模型。通过在填土表面施加均布荷载,研究了挡土墙的力学行为。模拟结果显示：填土中的水平应力随深度的增加对均布荷载的敏感性降低,填土深度越浅,应力值越早得到增加,增加幅度越大,这主要是由于石笼挡土墙上部发生了较大的位移,均布荷载的作用大部分消耗在上层填土中,仅有小部分能量传递到较深填土;石笼挡土墙变形主要发生在填土阶段,应加强石笼挡土墙填土阶段的施工管理;无论是在填土阶段,还是均布荷载作用阶段,石笼挡土墙墙面的位移值均是中上部较大,下部较小。由此可推断,墙背阶梯式石笼挡土墙的破坏形式为倾覆式。     

桩承式路堤承载特性的颗粒流模拟

现有的土拱效应计算方法中,由于采用的计算模型不同,计算结果差异很大。文中克服传统连续介质力学模型的宏观连续性假设,采用二维颗粒流程序（PFC）建立基于模型试验的细观数值分析模型,对桩承式路堤中土体接触力、应力分布、主应力方向、竖向位移进行分析,并比较计算和实测结果,研究土拱效应的荷载传递机制。同时,对不同桩帽、桩间距、填土高度、颗粒大小、摩擦角的情况进行PFC方法的参数敏感性分析。研究结果表明,桩承式路堤桩顶处局部范围可按弹性核考虑;土拱的分布型式受桩帽型式、桩净距、格栅的影响;实际土拱作用的影响范围主要集中在路堤底面以上约1倍桩净距的区域;土拱内部的竖向应力和水平应力均随深度非线性改变,桩土应力比随着荷载水平、土体内摩擦角、颗粒大小的增大而增加。