CN111931369B - 降雨型滑坡稳定性分析及运动距离测算方法、设备及介质 - Google Patents
降雨型滑坡稳定性分析及运动距离测算方法、设备及介质 Download PDFInfo
- Publication number
- CN111931369B CN111931369B CN202010779873.2A CN202010779873A CN111931369B CN 111931369 B CN111931369 B CN 111931369B CN 202010779873 A CN202010779873 A CN 202010779873A CN 111931369 B CN111931369 B CN 111931369B
- Authority
- CN
- China
- Prior art keywords
- boundary
- unit
- slope
- boundary unit
- landslide
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 16
- 238000000691 measurement method Methods 0.000 title claims abstract description 11
- 239000002689 soil Substances 0.000 claims abstract description 49
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 30
- 238000000034 method Methods 0.000 claims abstract description 28
- 238000006073 displacement reaction Methods 0.000 claims abstract description 27
- 238000004364 calculation method Methods 0.000 claims abstract description 25
- 230000008569 process Effects 0.000 claims abstract description 9
- 239000000758 substrate Substances 0.000 claims abstract description 6
- 238000012876 topography Methods 0.000 claims abstract description 5
- 239000011148 porous material Substances 0.000 claims description 7
- 230000001133 acceleration Effects 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 6
- 230000035699 permeability Effects 0.000 claims description 5
- 238000003860 storage Methods 0.000 claims description 5
- 210000003888 boundary cell Anatomy 0.000 claims description 4
- 210000004027 cell Anatomy 0.000 claims description 4
- 238000005520 cutting process Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 17
- 238000004590 computer program Methods 0.000 description 7
- 230000008595 infiltration Effects 0.000 description 6
- 238000001764 infiltration Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 230000006355 external stress Effects 0.000 description 3
- 230000005012 migration Effects 0.000 description 3
- 238000013508 migration Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000002265 prevention Effects 0.000 description 2
- 239000002352 surface water Substances 0.000 description 2
- 101710154918 Trigger factor Proteins 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 230000001687 destabilization Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000009736 wetting Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Abstract
本发明公开了一种降雨型滑坡稳定性分析及运动距离测算方法、设备及介质,获取待预测区域的基础地形、地质参数、降雨强度参数,滑坡点的位置参数;建立刚性块体模型,当边坡安全系数Fs<1时边坡失稳;将判断的已失稳边坡划分边界单元和滑体单元,计算边界单元所受土压力P;选用摩擦型流变基底阻力模型,获得摩阻力T;确定边界单元所受合力F;根据边界单元所受合力F计算滑体、边界单元的速度、位移、高度及土压力系数梯度;将计算得到的边界单元土压力系数梯度、速度、位移、高度代入边界单元所受土压力P,依次迭代计算得到滑坡运动过程中的速度和距离。本发明计算精确,广泛适用于降雨地区的路基、路堑、水利水电及建筑基坑边坡工程。
Description
技术领域
本发明属于地质灾害防治安全技术领域,具体涉及一种降雨型滑坡稳定性分析及运动距离测算方法、设备及介质,广泛适用于降雨型滑坡稳定性分析和运动范围评估。
背景技术
滑坡是一种严重的地质灾害。据统计,我过已受到和可能受到滑坡灾害威胁的地区占全国陆地面积的1/5~1/4,而西南地区发生的滑坡约占全国滑坡次数的一半以上。随着国民经济建设和大型工程的修建,滑坡灾害有逐年增加的趋势。
降雨引起滑坡称为降雨型滑坡,滑坡的发生固然离不开不良的地质条件和岩土体结构条件等内因,但另一方面离不开各种外部应力的作用,引起滑坡的外部应力有很多种,在诸多外部应力中降雨,尤其是大量的降雨或暴雨无疑是触发滑坡的最主要的因素之一。降雨对滑坡的触发作用主要体现在:在降雨过程中,地表水进入岩土体増加了坡体的自重,增大了孔隙水压力,使原本处于极限平衡状态的土体滑动;降雨由地表进入地下,转变成为地下水,浸泡软化滑动面,降低了坡体的抗剪强度;坡体多次的于湿交替变化导致岩土体开裂,产生大量的节理裂隙。研究表明,在具备了滑坡的地势和地质条件时,降雨在很大程度上是滑坡发生的最主要激发因素。
现有滑坡运移距离测算方法多采用经验统计法,该方法较为粗糙、无法准确获得滑坡运移距离。近年来得到较快发展的滑坡力学计算方法,如连续介质模型、条分模型、离散元模型等仍处于研究阶段,其中条分模型具有计算效率高、适用范围广的特点,逐渐被广大工程技术人员所采用,然而,采用条分模型的滑体内部相互作用是当前研究的重难点,这对于滑坡运移距离影响较大。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足,提供一种降雨型滑坡稳定性分析及运动距离测算方法、设备及介质,能够在降雨地区内根据降雨量的变化对降雨型滑坡的稳定性进行快速评价,对滑坡运动的范围进行测算。
本发明采用以下技术方案:
一种降雨型滑坡稳定性分析及运动距离测算方法,包括以下步骤:
S1、获取待预测区域的基础地形、地质参数、降雨强度参数,滑坡点的位置参数;
S2、根据步骤S1获得的参数建立刚性块体模型,根据降雨强度计算边坡安全系数Fs,确定降雨强度p与边坡安全系数Fs之间的关系,当Fs<1时,边坡失稳;
S3、将步骤S2判断的已失稳边坡划分边界单元和滑体单元,计算边界单元所受土压力P;选用摩擦型流变基底阻力模型,获得摩阻力T;确定边界单元所受合力F;根据边界单元所受合力F计算滑体、边界单元的速度、位移、高度及土压力系数梯度;将计算得到的边界单元土压力系数梯度、速度、位移、高度代入边界单元所受土压力P,依次迭代计算得到滑坡运动过程中的速度和距离。
具体的,步骤S2中,降雨强度p与边坡安全系数Fs之间关系如下:
其中,α为边坡倾角,W为滑体模型重量,γw为水的重度,为内摩擦角,Ks为饱和渗透系数,m=1-1/n,n为曲线形状参数。
具体的,步骤S3中,边界单元所受土压力P具体为:
其中,为第i边界单元土压力系数梯度;ac为第i边界单元离心加速度;g为重力加速度;Hi为第i边界单元的高度,α为边坡倾角。
具体的,步骤S3中,边界单元基底摩阻力T具体为:
其中,ru为孔隙水压力系数;φ为外摩擦角;
边界单元所受合力F具体为:
F=γHisinα+P-T
其中,γ为水的重度,Hi为第i边界单元的高度,α为边坡倾角。
具体的,步骤S3中,边界单元速度V为:
其中,Vi、Vi'为第i边界单元新、旧速度;Δt为计算时间步长,Hi为第i边界单元的高度,γ为水的重度,F为边界单元所受合力。
具体的,步骤S3中,边界、滑体单元位移S、s为:
其中,Si、Si'为第i边界单元新、旧位移;Si、Si+1为第i、i+1边界单元位移;sj为第j滑体单元位移,Vi、Vi'为第i边界单元新、旧速度。
具体的,步骤S3中,滑体、边界单元的高度h、H为:
其中,Dj为第j滑体单元体积,hj为第j滑体单元的高度,Si为第i边界单元新位移,Hi为第i边界单元的高度。
具体的,步骤S3中,第i边界单元土压力系数梯度为:
其中,kj为第j滑体单元的新土压力系数,hj为第j滑体单元的高度,Hi为第i边界单元的高度,sj为第j滑体单元位移,Si为第i边界单元新位移。
本发明的另一个技术方案是,一种计算设备,包括:
一个或多个处理器、存储器及一个或多个程序,其中一个或多个程序存储在所述存储器中并被配置为所述一个或多个处理器执行,所述一个或多个程序包括用于执行所述的方法中的任一方法的指令。
本发明的另一个技术方案是,一种存储一个或多个程序的计算机可读存储介质,所述一个或多个程序包括指令,所述指令当由计算设备执行时,使得所述计算设备执行所述的方法中的任一方法。
与现有技术相比,本发明至少具有以下有益效果:
本发明一种降雨型滑坡稳定性分析及运动距离测算方法,将降雨强度p代入计算得到降雨下边坡稳定性系数,根据稳定性系数值来判断边坡是否失稳,能够根据降雨量的变化对降雨型滑坡的稳定性进行快速评价,对滑坡运动的范围进行测算,其计算方法较为简洁、明了,预测得到滑坡运动距离更为精确。
进一步的,边坡稳定性系数Fs<1时,边坡失稳,滑坡体可分为多个边界单元及条块单元,
边界单元受后部条块单元推动,计算边界单元所受土压力P。
进一步的,选用摩擦型模型计算边界单元基底阻力T,以便进行下一步边界单元合力计算,结合边界单元所受土压力P、基底阻力T,计算边界单元所受合力F,以便下一步计算边界单元运动速度V。
进一步的,根据边界单元合力F计算得到其运动速度V,以便下一步计算边界及滑体单元位移S、s。
进一步的,根据边界单元运动速度V计算得到边界及滑体单元位移S、s,以便下一步计算边界及滑体单元高度h、H,根据边界及滑体单元位移计算边界及滑体单元高度h、H,为下一步计算边界单元土压力系数梯度做准备。
进一步的,根据边界及滑体单元的位移及高度,计算滑体单元土压力系数k及边界单元土压力系数梯度将计算得到的土压力系数梯度、边界单元速度及边界单元高度带入权利1中的公式计算下一时步边界单元所受土压力P。
综上所述,本发明方法广泛适用于降雨地区的路基、路堑、水利水电及建筑基坑边坡工程,能够快速计算边坡稳定性及失稳后滑坡运动距离,简洁明了,能够被广大工程人员理解运用。同时,计算精确,对于降雨地区边坡施工、滑坡灾害治理及边坡防灾减灾工作的开展具有十分重要的现实意义。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
图1为本发明流程示意图;
图2为降雨型滑坡稳定性分析计算流程示意图;
图3为滑坡运动距离计算流程示意图;
图4为边坡降雨入渗计算简图;
图5为降雨型滑坡稳定性分析的块体受力计算简图;
图6为运动距离测算分析的边界、滑体单元计算模型图;
图7为图6中边界、滑体单元刚度计算示意图;
图8为实例中滑坡运动距离计算中滑坡运动示意图;
图9为实例中滑坡运动距离计算中滑坡运动示意图。
具体实施方式
请参阅图1,本发明一种降雨型滑坡稳定性分析及运动距离测算方法,包括边坡稳定性计算及滑坡运动距离计算,包括以下步骤:
S1、获取待预测区域的基础地形、地质参数、降雨强度参数,滑坡点的位置参数;
通过获取待预测区域的基础地形、地质参数、降雨强度参数,滑坡点的位置参数建立刚性块体模型,根据降雨强度计算稳定系数判断边坡是否失稳,当边坡失稳后对已失稳的滑坡划分边界单元和滑体单元,计算滑坡运动过程中的速度和距离。
S2、建立刚性块体模型,根据降雨强度计算稳定系数,当Fs<1时,滑坡失稳,请参阅图2,具体包括:
S201、建立刚性块体模型;
刚性块体模型,即将边坡简化为一个块体,不考虑其变形特性。
S202、计算降雨工况下滑坡模型的力平衡公式,请参阅图5,具体计算如下:
其中,Pw为孔隙水压力;W为滑体模型重量;为内摩擦角;α为边坡倾角。
土体水力学特性
Van Genuchten关于土-水特征曲线的函数为:
其中,θ为体积含水率(%);h为压力水头(m);|hw|为h的绝对值;θs为饱和含水率(%);θr为残余含水率(%);β和n为曲线形状参数,且n>1,β的单位是(m-1);m=1-1/n。
相应的水力传导方程的函数表达式为:
其中,K为渗透系数(m/s);Ks为饱和渗透系数(m/s)。
根据假定,假设传导区含水率均匀分布,其计算模型如图4所示;图中x*、y*分别为平行与垂直坡面的坐标方向;θw为降雨入渗情况下传导区土体的体积含水率;θi为边坡的初始体积含水率。
设坡面上z*向的入渗率为i,并假定坡面为参考面,根据达西定律可得
其中,θw为降雨入渗情况下传导区土体的体积含水率;θi为边坡的初始体积含水率;
hw为含水率θw对应的压力水头;z为竖直坐标方向;z*为垂直坡面的坐标方向。
根据坐标变换可得
z=x*sinα+z*cosα
联合可得
坡面入渗率由两部分组成,一部分由重力梯度所引起,另一部分由水压梯度所引起。
由于降雨强度小于饱和渗透系数,因此,不会出现坡面积水,则坡面上z*方向的入渗率为i=pcosα
联立可得
求得湿润锋以上土体的压力水头hw,进一步求得此时的孔隙水压力为:
pw=γwhw
联立得到降雨强度p与孔隙水压力Pw之间关系:
根据边坡安全系数Fs为:
联立得到降雨强度p与边坡安全系数Fs之间关系,通过降雨强度判断边坡是否处于失稳状态,计算公式如下:
其中,γw为水的重度。
S3、对已失稳的滑坡划分边界单元和滑体单元,计算滑坡运动过程中的速度和距离,请参阅图3,具体包括:
S301、将块体滑坡模型划分成边界单元和滑体单元,其中,边界单元为无厚度单元;
S302、按下式计算边界单元所受土压力P;
其中,为第i边界单元土压力系数梯度,与切向应变增量Δε相关,见步骤S308;ac为第i边界单元离心加速度,计算公式/>取决于单元速度vi及滑动路径R的垂直曲率半径;g为重力加速度,取g=9.8m/s2;Hi为第i边界单元的高度;
S303、选用摩擦型流变基底阻力模型,获得摩阻力T;
其中,ru为孔隙水压力系数;φ为摩擦角,与滑体基底层下部岩土材料相关。
S304、计算边界单元所受合力F;
F=γHisinα+P-T
S305、计算边界单元速度V;
其中,Vi、Vi'为第i边界单元新、旧速度;Δt为计算时间步长;
S306、计算边界、滑体单元位移S、s;
其中,Si、Si'为第i边界单元新、旧位移;Si、Si+1为第i、i+1边界单元位移;sj为第j滑体单元位移;
S307、计算滑体、边界单元的高度h、H;
其中,Dj为第j滑体单元体积。
S308、请参阅图6和图7计算滑体单元土压力系数k及边界单元土压力系数梯度
S3081、计算滑体单元切向应变增量Δε,如下式:
S3082、运用计算机对滑体单元主、被动土压力系数ka、kp计算如下:
其中,δ为第j边界单元底部摩擦角;
S3083、计算刚度系数Sc;
当Δεj<0时,即压缩状态,Sc=(ka-kp)/0.05;
当Δεj≥0时,即卸载状态,Sc=(ka-kp)/0.025。
S3084、计算第j滑体单元土压力系数k为:
kj=k'j+ksΔεj
其中,kj、kj'为第j滑体单元的新、旧土压力系数;
S3085、计算第i边界单元土压力系数梯度
将计算得到的土压力系数梯度边界单元速度V及边界单元高度H代入步骤S301计算边界单元所受土压力P。
将相关参数更新后,代入步骤S302,依次迭代计算得到滑坡运动过程中的速度和距离,当迭代计算至各边界单元速度为0时,滑坡运动停止,将此时各边界及滑体单元位移相加得到滑坡总运动距离。
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。通常在此处附图中的描述和所示的本发明实施例的组件可以通过各种不同的配置来布置和设计。因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例
本发明一种降雨型滑坡稳定性分析及运动距离测算方法,选取某降雨型滑坡,测定滑体的密度、坡度、内摩擦角、降雨强度、含水量等参数;将边坡简化为块体计算模型,计算降雨状态下滑体的安全系数,计算对应于降雨强度p的边坡安全系数Fs,当Fs<1时,将滑坡体划分为边界单元及滑体单元,计算边界单元受到土压力及摩阻力,进而确定边界单元所受合力;根据边界单元所受合力计算滑体、边界单元的速度、位移、高度及土压力系数梯度等;将计算得到的边界单元土压力系数梯度、速度、高度代入下一时步骤边界单元土压力计算公式,依次迭代计算。
本发明已在滑坡动力学软件中实现,某边坡高度55m,降雨强度678mm,软件模拟滑体长度33.4m,如图8所示,共划分为32段滑体单元,33段边界单元,从平均坡角33°边坡由上而下滑动,内摩擦角10°,滑块底部摩擦角20°,边坡滑行状态如图8和图9所示,边坡滑动水平距离约54.1m,图8为滑坡运动发生前滑体示意图,原图中深色区域为滑坡体;图9为滑坡运动结束后滑体示意图,原图中深色区域为滑坡体。
综上所述,本发明一种降雨型滑坡稳定性分析及运动距离测算方法,通过边坡坡度、降雨强度等参数,计算降雨状态下滑体的安全系数,当边坡安全系数小于1时,计算失稳后滑体运动距离。本发明专利相比于其他降雨型滑坡稳定性分析方法更为简洁、明了,计算得到滑坡运动距离相比其他方法也更为精确。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上内容仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明权利要求书的保护范围之内。
Claims (4)
1.一种降雨型滑坡稳定性分析及运动距离测算方法,其特征在于,包括以下步骤:
S1、获取待预测区域的基础地形、地质参数、降雨强度参数,滑坡点的位置参数;
S2、根据步骤S1获得的参数建立刚性块体模型,根据降雨强度计算边坡安全系数Fs,确定降雨强度p与边坡安全系数Fs之间的关系,当Fs<1时,边坡失稳;刚性块体模型,即将边坡简化为一个块体,不考虑其变形特性,降雨强度p与边坡安全系数Fs之间关系如下:
其中,α为边坡倾角,W为滑体模型重量,γw为水的重度,为内摩擦角,Ks为饱和渗透系数,m=1-1/n,n为曲线形状参数;
S3、将步骤S2判断的已失稳边坡划分边界单元和滑体单元,计算边界单元所受土压力P;选用摩擦型流变基底阻力模型,获得摩阻力T;确定边界单元所受合力F;根据边界单元所受合力F计算滑体、边界单元的速度、位移、高度及土压力系数梯度;将计算得到的边界单元土压力系数梯度、速度、位移、高度代入边界单元所受土压力P,依次迭代计算得到滑坡运动过程中的速度和距离;
边界单元所受土压力P具体为:
其中,为第i边界单元土压力系数梯度;ac为第i边界单元离心加速度;g为重力加速度;Hi为第i边界单元的高度,α为边坡倾角;
边界单元速度V为:
其中,Vi、Vi'为第i边界单元新、旧速度;Δt为计算时间步长,Hi为第i边界单元的高度,γ为水的重度,F为边界单元所受合力;
边界、滑体单元位移S、s为:
其中,Si、Si'为第i边界单元新、旧位移;Si、Si+1为第i、i+1边界单元位移;sj为第j滑体单元位移,Vi、Vi'为第i边界单元新、旧速度;
滑体、边界单元的高度h、H为:
其中,Dj为第j滑体单元体积,hj为第j滑体单元的高度,Si为第i边界单元新位移,Hi为第i边界单元的高度;
第i边界单元土压力系数梯度为:
其中,kj为第j滑体单元的新土压力系数,hj为第j滑体单元的高度,Hi为第i边界单元的高度,sj为第j滑体单元位移,Si为第i边界单元新位移。
2.根据权利要求1所述的方法,其特征在于,步骤S3中,选用摩擦型流变基底阻力模型,边界单元的摩阻力T计算如下:
其中,ru为孔隙水压力系数;φ为外摩擦角;
边界单元所受合力F具体为:
F=γHi sinα+P-T
其中,γ为水的重度,Hi为第i边界单元的高度,α为边坡倾角。
3.一种计算设备,其特征在于,包括:
一个或多个处理器、存储器及一个或多个程序,其中一个或多个程序存储在所述存储器中并被配置为所述一个或多个处理器执行,所述一个或多个程序包括用于执行根据权利要求1或2所述的方法中的任一方法的指令。
4.一种存储一个或多个程序的计算机可读存储介质,其特征在于,所述一个或多个程序包括指令,所述指令当由计算设备执行时,使得所述计算设备执行根据权利要求1或2所述的方法中的任一方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010779873.2A CN111931369B (zh) | 2020-08-05 | 2020-08-05 | 降雨型滑坡稳定性分析及运动距离测算方法、设备及介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010779873.2A CN111931369B (zh) | 2020-08-05 | 2020-08-05 | 降雨型滑坡稳定性分析及运动距离测算方法、设备及介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111931369A CN111931369A (zh) | 2020-11-13 |
CN111931369B true CN111931369B (zh) | 2024-04-09 |
Family
ID=73307518
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010779873.2A Active CN111931369B (zh) | 2020-08-05 | 2020-08-05 | 降雨型滑坡稳定性分析及运动距离测算方法、设备及介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111931369B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112395768B (zh) * | 2020-11-25 | 2022-04-19 | 成都理工大学 | 一种地震边坡永久位移计算方法 |
CN112733412B (zh) * | 2020-12-28 | 2023-03-24 | 河海大学 | 一种水动力作用滑坡运动机制研究的速度等效表征方法 |
CN113065185B (zh) * | 2021-03-31 | 2022-10-25 | 上海交通大学 | 用于膨胀土边坡防护的土工编织袋结构单元及设计方法 |
CN116245283A (zh) * | 2023-03-08 | 2023-06-09 | 北京七兆科技有限公司 | 弃渣场次生灾害风险评估方法、装置、设备及存储介质 |
CN116739184B (zh) * | 2023-08-08 | 2023-11-07 | 四川川核地质工程有限公司 | 一种滑坡预测方法及系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102306233A (zh) * | 2011-06-15 | 2012-01-04 | 浙江大学 | 降雨作用下流域滑坡时空预测方法 |
KR101580062B1 (ko) * | 2014-12-05 | 2015-12-28 | 연세대학교 산학협력단 | 실시간 산사태 예경보 방법 및 시스템 |
CN107066771A (zh) * | 2017-06-13 | 2017-08-18 | 成都理工大学 | 一种平推式滑坡运动距离计算方法及应用 |
CN107688895A (zh) * | 2017-08-09 | 2018-02-13 | 中国地质环境监测院 | 一种冻融型滑坡安全性分析及运动距离测算方法 |
WO2018119880A1 (zh) * | 2016-12-29 | 2018-07-05 | 柳成荫 | 一种基于降雨量及土壤湿度的降雨型滑坡预警方法及装置 |
CN108984821A (zh) * | 2018-05-25 | 2018-12-11 | 重庆交通大学 | 基于降雨诱发变形的结构面控制边坡稳定性评价方法 |
CN111210073A (zh) * | 2020-01-06 | 2020-05-29 | 四川省公路规划勘察设计研究院有限公司 | 一种滑坡灾害预测方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20180292299A1 (en) * | 2014-01-13 | 2018-10-11 | Hubei University Of Technology | Method of critical displacement forecast based on the deformation failure mechanism of slope |
-
2020
- 2020-08-05 CN CN202010779873.2A patent/CN111931369B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102306233A (zh) * | 2011-06-15 | 2012-01-04 | 浙江大学 | 降雨作用下流域滑坡时空预测方法 |
KR101580062B1 (ko) * | 2014-12-05 | 2015-12-28 | 연세대학교 산학협력단 | 실시간 산사태 예경보 방법 및 시스템 |
WO2018119880A1 (zh) * | 2016-12-29 | 2018-07-05 | 柳成荫 | 一种基于降雨量及土壤湿度的降雨型滑坡预警方法及装置 |
CN107066771A (zh) * | 2017-06-13 | 2017-08-18 | 成都理工大学 | 一种平推式滑坡运动距离计算方法及应用 |
CN107688895A (zh) * | 2017-08-09 | 2018-02-13 | 中国地质环境监测院 | 一种冻融型滑坡安全性分析及运动距离测算方法 |
CN108984821A (zh) * | 2018-05-25 | 2018-12-11 | 重庆交通大学 | 基于降雨诱发变形的结构面控制边坡稳定性评价方法 |
CN111210073A (zh) * | 2020-01-06 | 2020-05-29 | 四川省公路规划勘察设计研究院有限公司 | 一种滑坡灾害预测方法及装置 |
Non-Patent Citations (2)
Title |
---|
强降雨作用下浅层滑坡的入渗及稳定性;王进;冷先伦;阮航;刘世伟;崔臻;;东南大学学报(自然科学版)(S1);全文 * |
连续强降雨工况土质边坡非饱和渗流及稳定性分析;邹文华;刘辉;邓小钊;刘耀坤;;中外公路(06);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111931369A (zh) | 2020-11-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111931369B (zh) | 降雨型滑坡稳定性分析及运动距离测算方法、设备及介质 | |
Tang et al. | A novel approach for determining landslide pushing force based on landslide-pile interactions | |
CN105224704B (zh) | 一种海底不稳定性的评价方法 | |
Chen et al. | A dynamic model for rainfall-induced landslides on natural slopes | |
Chai et al. | Coefficient of consolidation from non-standard piezocone dissipation curves | |
CN103135128A (zh) | 一种地震荷载作用下三维边坡稳定性预测方法 | |
Zhang et al. | New stability calculation method for rock slopes subject to flexural toppling failure | |
Abe et al. | Numerical study on dynamic behavior of slope models including weak layers from deformation to failure using material point method | |
CN110261573A (zh) | 一种高位岩质滑坡稳定性动态数值评价方法 | |
CN103469780A (zh) | 一种重力坝坝基深层抗滑稳定临界滑动面的计算方法 | |
CN103163563A (zh) | 一种三维边坡稳定性预测方法 | |
CN111553101A (zh) | 一种隧道开挖上覆岩层开裂预报方法及围岩支护方法 | |
Li et al. | One-dimensional large-strain consolidation of soft clay with non-Darcian flow and nonlinear compression and permeability of soil | |
Cao et al. | Analysis of large-strain consolidation behavior of soil with high water content in consideration of self-weight | |
Nguyen et al. | Postfailure characterization of shallow landslides using the material point method | |
Takaine et al. | Predictive simulation of deformation and failure of peat-calcareous soil layered ground due to multistage test embankment loading | |
Di et al. | An operator‐split ALE model for large deformation analysis of geomaterials | |
CN113720992B (zh) | 一种利用雨滴下落法模拟降雨作用对岩土体影响的方法 | |
Badakhshan et al. | A simplified method for prediction of ultimate bearing capacity of eccentrically loaded foundation on geogrid reinforced sand bed | |
Smith et al. | Modified Boussinesq's equations for nonuniform tire loading | |
Talbot et al. | Modeling large-deformation features of the Lower San Fernando Dam failure with the Material Point Method | |
Lu et al. | Numerical and experimental analyses for bearing capacity of rigid strip footing subjected to eccentric load | |
CN102930148B (zh) | 一种基于随机起动的管涌渗透系数确定方法 | |
CN108319761A (zh) | 确定群桩位移的方法和装置 | |
Singh et al. | Evaluating the slope instability of the Amiya Slide |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant |