CN111753440A - 一种高位滑坡冲击铲刮变量的获取方法 - Google Patents
一种高位滑坡冲击铲刮变量的获取方法 Download PDFInfo
- Publication number
- CN111753440A CN111753440A CN202010655731.5A CN202010655731A CN111753440A CN 111753440 A CN111753440 A CN 111753440A CN 202010655731 A CN202010655731 A CN 202010655731A CN 111753440 A CN111753440 A CN 111753440A
- Authority
- CN
- China
- Prior art keywords
- scraping
- impact
- stress
- contact
- radius
- 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.)
- Granted
Links
Images
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
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种高位滑坡冲击铲刮变量的获取方法,首先在高位滑坡冲击铲刮层时,将冲击力分为竖向冲击力和切向冲击力,采用Hertz弹性接触理论获得竖向冲击力,采用基底摩擦阻力理论获得切向冲击力;采用莫尔库仑剪切破坏屈服准则获得铲刮层的被动屈服;根据塑性区极限应力的分布位置得出最大滑裂线位置,划分塑性范围;根据上述获得的滑体冲击铲刮层后的动态变化,获得滑坡冲击铲刮模式。上述方法能有效获取高位滑坡冲击中的各种铲刮变量,并解决动力后效应问题。
Description
技术领域
本发明涉及滑坡冲击铲刮研究技术领域,尤其涉及一种高位滑坡冲击铲刮变量的获取方法。
背景技术
目前,国内外许多研究者对滑坡灾害的动力侵蚀相关领域开展了一定研究,为揭示滑坡的动力致灾机理奠定了基础。事实上,不同类型的滑坡动力侵蚀的作用方式也有所不同,进一步细分可分为两种:一种为流化滑坡(泥石流)的裹挟侵蚀,另一种为固相态滑坡的冲击铲刮。目前大量的研究成果聚焦于泥石流运动过程中的侵蚀裹挟效应研究,并且大多数是采用基于经验的侵蚀率方法进行评估,更多是依靠人为的定性评估。然而对于滑坡、崩塌等固相体冲击作用下的铲刮理论研究尚有不足,这种冲击、撞击作用不仅会使原有滑体速度改变、体积增加、粉碎化,甚至可能会铲动周边斜坡形成次级滑坡,会直接导致高位滑坡的滑体体积、运动距离和堆积特征发生改变,影响着滑坡危险区域预测预判的精准度,使得滑坡发生规模放大,灾害救援难度增加,然而这个动力作用的规律和机理特性却不为人们所熟知。
针对滑坡冲击铲刮的分析方法主要有:一种为流化滑坡(泥石流)的侵蚀率算法,该方法建立时间和铲刮体积的关系式,人为设定铲刮率和铲刮区域,通过体积的增加量和固定的底面积值来反演计算铲刮深度,但该方法大多基于人为经验,力学机理和计算精度较为欠缺;另一种为基于农业机械学犁耕阻力的犁切算法,该方法人为设定犁切形状和铲刮几何形态,再根据极限平衡法进行计算,但该方法缺少相应的铲刮模式和动力学基础,对速度和时间的考虑较为欠缺。
发明内容
本发明的目的是提供一种高位滑坡冲击铲刮变量的获取方法,该方法能有效获取高位滑坡冲击中的各种铲刮变量,并解决动力后效应问题。
本发明的目的是通过以下技术方案实现的:
一种高位滑坡冲击铲刮变量的获取方法,所述方法包括:
步骤1、在高位滑坡冲击铲刮层时,将冲击力分为竖向冲击力和切向冲击力,采用Hertz弹性接触理论获得竖向冲击力,采用基底摩擦阻力理论获得切向冲击力;
步骤2、采用莫尔库仑剪切破坏屈服准则获得铲刮层的被动屈服;
步骤3、根据塑性区极限应力的分布位置得出最大滑裂线位置,划分塑性范围;
步骤4、根据步骤1-3获得的滑体冲击铲刮层后的动态变化,获得滑坡冲击铲刮模式。
由上述本发明提供的技术方案可以看出,上述方法能有效获取高位滑坡冲击中的各种铲刮变量,并解决动力后效应问题。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他附图。
图1为本发明实施例提供的高位滑坡冲击铲刮变量的获取方法流程示意图。
具体实施方式
下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明的保护范围。
下面将结合附图对本发明实施例作进一步地详细描述,如图1所示为本发明实施例提供的高位滑坡冲击铲刮变量的获取方法流程示意图,所述方法包括:
步骤1、在高位滑坡冲击铲刮层时,将冲击力分为竖向冲击力和切向冲击力,采用Hertz弹性接触理论获得竖向冲击力,采用基底摩擦阻力理论获得切向冲击力;
在该步骤中,在高位滑坡冲击铲刮层时,冲击铲刮的接触也根据滑体的运动分为了两个阶段:第一阶段是以滑体前缘的冲击为主,滑体同可铲刮层接触后竖向荷载大于切向荷载;第二阶段是以滑体底部的剪切为主,竖向荷载主要是滑体自重提供,切向荷载大于竖向荷载。
将冲击力分为竖向冲击力和切向冲击力,首先采用Hertz弹性接触理论获得竖向冲击力,具体过程为:
接触压力与变形之间的关系为:
式中,E为等效弹性模量,且E1为撞击体的弹性模量,E2为铲刮层(被撞击体)的弹性模量;R为等效球体半径,且R1为撞击体半径,R2为铲刮层(被撞击体)半径,若铲刮层为地面,则半径为无限大,等效半径大小为撞击体半径的大小;δn为接触变形量;ν为等效泊松比,ν1为撞击体的泊松比;ν2为铲刮层(被撞击体)的泊松比;
接触压力P与接触半径之间的关系为:
式中,a为接触半径;r为接触半径变量,取值为0-a;P(r)为随接触半径r变化的接触压应力,当r=0时为最大接触压应力;
所述接触半径a与接触深度δn之间的关系为:
a2=Rδn
为了简化计算竖向冲击力可用最大接触应力进行分析,最大竖向冲击力表示为:
式中,E为等效弹性模量,R为等效半径,m为等效质量,vn为撞击体竖向速度。
另外,所述采用基底摩擦阻力理论获得切向冲击力的过程具体为:
对于基底摩擦阻力模型,根据竖向荷载得到基底切向阻力,阻力模型为同摩擦系数相关的摩擦模型。切向冲击力的计算需用到正压力值,正压力值为竖向Hertz接触压应力,表示为:
Pt=μPn
式中,Pt为正应力,μ为孔隙水压力,Pt为切应力。
步骤2、采用莫尔库仑剪切破坏屈服准则获得铲刮层的被动屈服;
在该步骤中,采用抗剪强度q与平均应力p之间为直线关系的屈服准则,具体为莫尔库仑剪切破坏屈服准则来获得铲刮层的被动屈服,具体为:
所述莫尔库仑剪切破坏屈服准则是描述剪切面上剪应力τ与该面上正应力σ间关系,在某个平面上τ-σ的最容易破坏的组合,公式为:
用应力不变量表示为:
步骤3、根据塑性区极限应力的分布位置得出最大滑裂线位置,划分塑性范围;
在该步骤中,具体是根据地表撞击荷载应力P求得铲刮层中半无限空间中任意点M点处的附加应力状态;
再求得该点处的主应力,根据莫尔库仑剪切破坏屈服准则作为塑性临界范围的判定条件,将大小主应力代入到极限平衡条件中,得出塑性区的轨迹方程,划分塑性范围,具体过程为:
1)竖向荷载导致M点的应力状态,根据弹性力学中Boussinesq解,求得M点附加应力状态如下:
上式中,m=x/2a,n=z/2a;Pn为地面受到的竖向荷载应力大小;x为M点横坐标大小;z为M点纵坐标大小;2a为受到竖向荷载宽度;
2)切向荷载导致M点的应力状态,根据弹性力学中Cerruti解,求得M点附加应力状态如下:
上式中,m=x/2a,n=z/2a;Pt为地面受到的切向荷载应力大小;x为M点横坐标大小;z为M点纵坐标大小;2a为受到切向荷载宽度;
3)塑性区轨迹方程计算:
竖向和切向共同导致的M点的附加应力状态为:
σz=σz,n+σz,t
σx=σx,n+σx,t
τxz=τxz,n+τxz,t
M点附加大小主应力计算公式为:
假设土层原有自重应力场的侧向土压力系数为K0=1,具有静水压力性质,不改变附加应力场的大小主应力大小和方向,M点大小主应力计算公式为:
σ1=σ′1+γz
σ3=σ′3+γz
MC莫尔库仑准则的极限平衡条件为:
再根据求得大小主应力σ1、σ3,代入到上述极限平衡条件中求解,得到塑性区的轨迹方程,划分塑性范围。
并根据竖向和切向荷载的大小关系,塑性滑移线可以分为表层滑动面和深层滑动面,即竖向冲击为主的铲刮为深层滑移线,切向裹挟为主的铲刮为浅层滑移线。
步骤4、根据步骤1-3获得的滑体冲击铲刮层后的动态变化,获得滑坡冲击铲刮模式。
在该步骤中,所述滑坡冲击铲刮模式包括嵌入铲起模式、裹挟刮带模式、平推滑移模式和冲击飞溅模式,具体来说:
1)嵌入铲起模式:滑坡下滑冲击软弱下垫层,滑体具有较大冲击力,两者初始接触状态为滑体前缘嵌入进下垫层中,塑性区以竖向发展为主,滑坡运动距离相对较近,大多数在撞击点位置处堆积;
2)裹挟刮带模式:该模式最为常见,滑体具有较大切向剪切力,滑体易对下垫层材料进行剪切运动,塑性区以切向发展为主,并刮带走大量底部被屈服破坏的材料,并汇入到运动的滑体中,使滑坡规模和运动距离放大;
3)平推滑移模式:该种类型下垫层岩土材料一般为泥页岩等性质较差的层状岩土体,且前部临空或阻挡较差,滑体下滑冲击后先是导致下垫层破坏,并平推下垫层整体向前运动;
4)冲击飞溅模式:该种类型下垫层材料一般为较为坚硬,同时滑体材料的完整性相对较好,滑体高位启动具有较大动能,下滑接触后经过能量传递,产生巨大的冲击能量,导致被撞击体松散解体后飞溅形成远距离运动。
值得注意的是,本发明实施例中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求书的保护范围为准。
Claims (6)
1.一种高位滑坡冲击铲刮变量的获取方法,其特征在于,所述方法包括:
步骤1、在高位滑坡冲击铲刮层时,将冲击力分为竖向冲击力和切向冲击力,采用Hertz弹性接触理论获得竖向冲击力,采用基底摩擦阻力理论获得切向冲击力;
步骤2、采用莫尔库仑剪切破坏屈服准则获得铲刮层的被动屈服;
步骤3、根据塑性区极限应力的分布位置得出最大滑裂线位置,划分塑性范围;
步骤4、根据步骤1-3获得的滑体冲击铲刮层后的动态变化,获得滑坡冲击铲刮模式。
2.根据权利要求1所述高位滑坡冲击铲刮变量的获取方法,其特征在于,在步骤1中,所述采用Hertz弹性接触理论获得竖向冲击力的过程为:
接触压力与变形之间的关系为:
式中,E为等效弹性模量,且E1为撞击体的弹性模量,E2为铲刮层的弹性模量;R为等效球体半径,且R1为撞击体半径,R2为铲刮层半径,若铲刮层为地面,则半径为无限大,等效半径大小为撞击体半径的大小;δn为接触变形量;ν为等效泊松比,ν1为撞击体的泊松比;ν2为铲刮层(被撞击体)的泊松比;
接触压力P与接触半径之间的关系为:
式中,a为最大接触半径;r为接触半径变量;P(r)为随接触半径r变化的接触压力;
所述接触半径a与接触深度δn之间的关系为:
a2=Rδn
为了简化计算竖向冲击力用最大接触应力进行分析,最大竖向冲击力表示为:
式中,E为等效弹性模量;R为等效半径;m为等效质量;vn为撞击体竖向速度。
3.根据权利要求1所述高位滑坡冲击铲刮变量的获取方法,其特征在于,在步骤1中,所述采用基底摩擦阻力理论获得切向冲击力的过程具体为:
所述切向冲击力的获取需用到正压力值,正压力值为竖向Hertz接触压应力,表示为:
Pt=μPn
式中,Pt为正应力;μ为孔隙水压力;Pt为切应力。
5.根据权利要求1所述高位滑坡冲击铲刮变量的获取方法,其特征在于,所述步骤3的过程具体为:
根据地表撞击荷载应力P求得铲刮层中半无限空间中任意点M点处的附加应力状态;
再求得该点处的主应力,根据莫尔库仑剪切破坏屈服准则作为塑性临界范围的判定条件,将大小主应力代入到极限平衡条件中,得出塑性区的轨迹方程,划分塑性范围。
6.根据权利要求1所述高位滑坡冲击铲刮变量的获取方法,其特征在于,在步骤4中,所述滑坡冲击铲刮模式包括嵌入铲起模式、裹挟刮带模式、平推滑移模式和冲击飞溅模式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010655731.5A CN111753440B (zh) | 2020-07-09 | 2020-07-09 | 一种高位滑坡冲击铲刮变量的获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010655731.5A CN111753440B (zh) | 2020-07-09 | 2020-07-09 | 一种高位滑坡冲击铲刮变量的获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111753440A true CN111753440A (zh) | 2020-10-09 |
CN111753440B CN111753440B (zh) | 2021-03-12 |
Family
ID=72711519
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010655731.5A Active CN111753440B (zh) | 2020-07-09 | 2020-07-09 | 一种高位滑坡冲击铲刮变量的获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111753440B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112924305A (zh) * | 2021-01-27 | 2021-06-08 | 中国地质大学(北京) | 获取被动桩板拦石墙在崩塌落石地质灾害中动态响应数据的方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011152741A1 (en) * | 2010-06-02 | 2011-12-08 | Pavlo Semenivskyi | Construction system for building industry, especially in the technology of fast assembling of skeleton type buildings |
CN106917634A (zh) * | 2017-01-13 | 2017-07-04 | 吕梁市希贤治黄研究所 | 矿山采空区泥浆回填的机械设备和方法 |
CN109325250A (zh) * | 2018-07-26 | 2019-02-12 | 四川大学 | 一种高速滑坡-碎屑流运动侵蚀效应的数值模拟方法及系统 |
CN110750860A (zh) * | 2019-09-11 | 2020-02-04 | 四川轻化工大学 | 一种土坡滑坡全过程分析方法 |
US20200056097A1 (en) * | 2018-08-15 | 2020-02-20 | Wei Ming Chang | Method of treating drilling mud with oil-based drilling fluid |
CN110836690A (zh) * | 2019-11-06 | 2020-02-25 | 中国地质调查局水文地质环境地质调查中心 | 基于变形量及含水率的高位远程滑坡预警装置及方法 |
CN111238767A (zh) * | 2020-01-21 | 2020-06-05 | 中国地质环境监测院 | 一种高位滑坡碎屑流抗冲拦挡结构模拟试验平台 |
CN111398067A (zh) * | 2020-04-01 | 2020-07-10 | 中国地质科学院地质力学研究所 | 一种高位滑坡冲击铲刮模式探索的滑槽模型试验方法 |
CN111475937A (zh) * | 2020-04-03 | 2020-07-31 | 中国地质科学院地质力学研究所 | 一种流固二相流流化滑坡的模拟仿真方法 |
-
2020
- 2020-07-09 CN CN202010655731.5A patent/CN111753440B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011152741A1 (en) * | 2010-06-02 | 2011-12-08 | Pavlo Semenivskyi | Construction system for building industry, especially in the technology of fast assembling of skeleton type buildings |
CN106917634A (zh) * | 2017-01-13 | 2017-07-04 | 吕梁市希贤治黄研究所 | 矿山采空区泥浆回填的机械设备和方法 |
CN109325250A (zh) * | 2018-07-26 | 2019-02-12 | 四川大学 | 一种高速滑坡-碎屑流运动侵蚀效应的数值模拟方法及系统 |
US20200056097A1 (en) * | 2018-08-15 | 2020-02-20 | Wei Ming Chang | Method of treating drilling mud with oil-based drilling fluid |
CN110750860A (zh) * | 2019-09-11 | 2020-02-04 | 四川轻化工大学 | 一种土坡滑坡全过程分析方法 |
CN110836690A (zh) * | 2019-11-06 | 2020-02-25 | 中国地质调查局水文地质环境地质调查中心 | 基于变形量及含水率的高位远程滑坡预警装置及方法 |
CN111238767A (zh) * | 2020-01-21 | 2020-06-05 | 中国地质环境监测院 | 一种高位滑坡碎屑流抗冲拦挡结构模拟试验平台 |
CN111398067A (zh) * | 2020-04-01 | 2020-07-10 | 中国地质科学院地质力学研究所 | 一种高位滑坡冲击铲刮模式探索的滑槽模型试验方法 |
CN111475937A (zh) * | 2020-04-03 | 2020-07-31 | 中国地质科学院地质力学研究所 | 一种流固二相流流化滑坡的模拟仿真方法 |
Non-Patent Citations (2)
Title |
---|
YANG GAO .ETC.: ""Dynamic characteristics of high-elevation and long-runout landslides in the Emeishan basalt area:a case study of the Shuicheng "7.23 "landslide in Guizhou, China"", 《LANDSLIDES》 * |
高杨: ""高速远程滑坡铲刮动力学分析"", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112924305A (zh) * | 2021-01-27 | 2021-06-08 | 中国地质大学(北京) | 获取被动桩板拦石墙在崩塌落石地质灾害中动态响应数据的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111753440B (zh) | 2021-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shen et al. | DEM analyses of rock block shape effect on the response of rockfall impact against a soil buffering layer | |
Rogers et al. | Simulation of caisson breakwater movement using 2-D SPH | |
Fielke | Finite element modelling of the interaction of the cutting edge of tillage implements with soil | |
CN110261573B (zh) | 一种高位岩质滑坡稳定性动态评价方法 | |
Lee et al. | An experimental study on earthquake fault rupture propagation through a sandy soil deposit | |
CN111753440B (zh) | 一种高位滑坡冲击铲刮变量的获取方法 | |
Hu et al. | Granular effects on depositional processes of debris avalanches | |
Mao et al. | Runout prediction and deposit characteristics investigation by the distance potential-based discrete element method: the 2018 Baige landslides, Jinsha River, China | |
Wang et al. | Comparison between the first and second wave crest amplitude generated by landslides | |
Zhang et al. | A theoretical model for the estimation of maximum impact force from a rockfall based on contact theory | |
Raines et al. | Driving characteristics of open-toe piles in dense sand | |
CN106503421A (zh) | 一种边坡渐进破坏全过程计算方法 | |
CN112241603B (zh) | 高位滑坡冲击铲刮及下垫层汇入过程的数值仿真方法 | |
CN112270022B (zh) | 基于橡胶-粗糙表面接触的摩擦力学建模方法 | |
Ferrari et al. | Towards the comprehension of rockfall motion, with the aid of in situ tests | |
CN109470581A (zh) | 露天矿山边坡岩体结构面抗剪强度分级确定方法 | |
Pan et al. | Gravitational stresses in long asymmetric ridges and valleys in anisotropic rock | |
Rossiter et al. | Assessment of ice-soil interactions: Continuum modeling in clays | |
Hammerström et al. | Designed high-friction surfaces—influence of roughness and deformation of the counter surface | |
Abe et al. | Numerical simulation of a series of flume tests with dry and wet sands by using depth averaged material point method | |
HAKUNO et al. | Application of the distinct element method to the numerical analysis of debris flows | |
Yue et al. | Numerical Simulation Analysis of the Characteristics of Torpedo Anchor Penetration in Cohesive Soil | |
Numada et al. | Material point method for run-out analysis of earthquake-induced long-traveling soil flows | |
Konagai et al. | Pseudo-three dimensional lagrangian particle finite difference method for modeling earthquake induced soil flows | |
Becker | Mechanics of the Panama Canal slides |
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 | ||
GR01 | Patent grant |