CN104268643B - 一种三维蠕变边坡滑动位移的预测方法 - Google Patents
一种三维蠕变边坡滑动位移的预测方法 Download PDFInfo
- Publication number
- CN104268643B CN104268643B CN201410470310.XA CN201410470310A CN104268643B CN 104268643 B CN104268643 B CN 104268643B CN 201410470310 A CN201410470310 A CN 201410470310A CN 104268643 B CN104268643 B CN 104268643B
- Authority
- CN
- China
- Prior art keywords
- cos
- shearing
- bar post
- beta
- alpha
- 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
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
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)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
Abstract
本发明公开了一种基于滑动位移分析的三维边坡稳定性预测方法,它包括选定具体待预测的三维蠕变边坡、将该三维蠕变边坡离散化、建立岩土体的西原体模型、获得条柱底面的剪力与关键点竖向位移的关系和求解不同时刻对应的关键点竖向位移等几个步骤,可以得到三维蠕变边坡的滑动位移随时间变化的关系。本发明的优点是:考虑了岩土体的蠕变特性和三维蠕变边坡的滑动位移信息,提高了计算精度,预测结果更为可靠。
Description
技术领域
本发明属于地质灾害防控技术领域,具体涉及一种三维蠕变边坡滑动位移的预测方法。
背景技术
岩土边坡的蠕变目前在国内外的研究尚少,但因对岩石特性包括岩土体流变特性研究不够而导致延误施工甚至工程失败的先例不胜枚举,意大利瓦依昂(Vajont)库岸蠕滑破坏是其中之一。自1966年第1届国际岩石力学会议以来的历次国际会议上,都有不少关于岩石流变性研究的论文。在第1届国际岩石力学会议上,Zischinsky用流变学模型描述了高边坡的变形,并指出岩土体的蠕变在高边坡变形中起重要作用。
工程实践与研究表明,岩土边坡工程的破坏与失稳,在许多情况下并不是在开挖形成以后立即发生,岩土体应力与变形是随时间变化而不断地调整,其调整的过程往往需要延续一个较长的时期。边坡蠕变是指组成边坡的岩体和土体在自重应力以及以水平应力为主的构造应力场的作用下,变形随时间而持续增加的性质。产生边坡变形的原因是多方面的,地质作用、地下水流、温度变化、植被作用等都可以产生边坡的宏观变形。但就岩土体本身而言,边坡与时间有关的变形主要是由岩土体蠕变引起的,因此,研究边坡变形时,应特别注意岩土体材料的蠕变特性与边坡蠕变特性之密切关系。
目前对蠕变边坡的稳定性分析,主要还局限于分析二维蠕变边坡问题,例如,“考虑岩土蠕变特性的边坡长期稳定性研究”,蒋海飞,胡斌,刘强,王新刚,金属矿山,2013年第12期,第131~157页,记载了采用考虑岩土蠕变特性的强度折减法对十堰某建筑场地边坡进行位移计算,得到了监测点在不同折减系数下的蠕变曲线。但是自然界的边坡工程问题通常具有三维特征,故用二维蠕变边坡的方法来分析三维蠕变边坡,与实际的情况存在明显差异,由此获得的滑动位移结果必然不可靠,得到的可靠性不能完全满足工程要求。
发明内容
针对现有技术用二维蠕变边坡的分析方法存在与实际状态的差异,本发明所要解决的技术问题就是提供一种三维蠕变边坡滑动位移的预测方法,它能够结合三维边坡滑动的位移量,计算获得该三维蠕变边坡的滑动位移随时间变化的曲线,从而提高预测的可靠性。
本发明所要解决的技术问题是通过这样的技术方案实现的,它包括有以下步骤:
步骤1、选定具体待预测的三维蠕变边坡,确定该三维蠕变边坡潜在滑裂面的几何尺寸,将潜在滑裂面和边坡表面用方程表示,确定岩土的土体指标参数;
步骤2、将三维蠕变边坡离散化,三维蠕变边坡被垂直离散为m行和n列条柱,每个条柱按所在的行号i和列号j定义为[i,j];假定行方向的条柱间作用力与水平面夹角为±α,假定列方向的条柱间作用力与水平面夹角为±β;
步骤3、根据西原体模型,确定第[i,j]条柱底面的剪力与该条柱的剪切位移之间的关系;利用离散后的各条柱的位移协调条件,得到第[i,j]条柱底面的剪力与关键点竖向位移Δ0的关系;
步骤4、根据3个力的平衡方程和3个力矩的平衡方程,并结合[i,j]条柱底面的剪力与关键点竖向位移的关系,建立一个含有未知数α、β、Δ0和蠕变时间t的非线性方程组;求解该方程组,得到不同的时刻所对应关键点的竖向位移Δ0;进而得到该三维蠕变边坡的滑动位移随时间的变化规律。
由于本发明在三维蠕变边坡离散化后,根据3个力平衡方程和3个力矩的平衡方程,并且考虑了三维蠕变边坡的滑动位移量,得到三维蠕变边坡的滑动位移随时间变化的曲线,提高了预测结果的可靠性。另外,所有建模过程都能程序化,便于操作和编程,通过计算机实现三维蠕变边坡滑动位移的预测,大大的减少了人为的计算量。所以本发明的优点是:考虑了岩土体的蠕变特性和蠕变边坡的滑动位移信息,提高了计算精度,预测结果更为可靠。
附图说明
本发明的附图说明如下:
图1为本发明的一个实施例的边坡潜在滑裂面和边坡表面的剖面图;
图2为三维蠕变边坡离散化结构图;
图3为离散化条柱的受力模型图;
图4为岩土体的西原体模型示意图;
图5为实施例的关键点竖向位移Δ0随时间变化的关系图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明:
步骤1、选定具体待预测的三维蠕变边坡,确定三维蠕变边坡潜在滑裂面的形状及边坡体几何尺寸,将边坡表面和潜在滑裂面用方程表示;边坡表面的参数有边坡斜面在水平面上的投影长度l,边坡斜面在竖直方向上的投影长度H;潜在滑裂面的参数依据实际几何形状确定,以及确定岩土的土体指标参数;
边坡表面的表达式为
潜在滑裂面的表达式为:
本实施例的三维蠕变边坡如图1所示,其边坡表面为倾斜面,取其倾斜角为45°,边坡斜面的水平方向的投影为20m,在竖直方向上的投影为20m;滑裂面为椭球体滑裂面,设y方向上滑裂面宽度的半长轴为40m;x、z方向的半轴长度为20m。该边坡的岩土体材料为各向同性均质材料,其指标参数为:粘聚力c=10kPa,内摩擦角容重γ=22.0kN/m3。
步骤2、将该三维蠕变边坡体离散化
如图2所示,将三维蠕变边坡离散化。该三维蠕变边坡被垂直离散为200行和200列条柱。每个条柱按所在的行号i和列号j定义为[i,j];假定行方向的条柱间作用力与水平面夹角为±α,假定列方向的条柱间作用力与水平面夹角为±β。第[i,j]条柱的受力模型如图3所示。
步骤3、建立岩土体的西原体模型;
(1)、确定第[i,j]条柱底面的剪力与该条柱的剪切位移之间的关系
建立岩土体的西原体模型如图4所示,三维蠕变边坡的岩土体的西原体模型参数为G1=43×106Pa,G2=37×106Pa,η1=3.739×1012Pa·s,η2=3.739×1012Pa·s。该岩土体的长期抗剪强度为τlong,长期抗剪强度为抗剪强度的70%,那么长期抗剪强度为为:
其中σi,j为第[i,j]条柱底面的正应力。
那么长期抗剪力为:
其中,Ai,j为第[i,j]条柱的横截面积,为第[i,j]条柱底面法向力的单位向量在z轴上的分量。
边坡条柱中的岩土体的西原体模型的数学表达式如下,
(a)当条柱底面的剪力Si,j小于长期抗剪力时,
bi,j是第[i,j]条柱行方向的宽度,并且有边界条件Δi,j(0)=0,可以得到第[i,j]条柱底面的剪力与该条柱的剪切位移之间的关系为:
式中Δi,j(t)为第[i,j]条柱底面的剪切位移,Si,j为第[i,j]条柱底面的剪力,G1、G2、η1、η2为西原体模型中的参数,t是发生蠕变的时间。
(b)当条柱底面的剪力Si,j大于等于长期抗剪力时,
并且有边界条件Δi,j(0)=0和可以得到第[i,j]条柱底面的剪力与该条柱的剪切位移之间的关系为:
其中Ai,j为第[i,j]条柱的横截面积,为第[i,j]条柱底面法向力的单位向量在z轴上的分量。
(2)、确定[i,j]条柱底面的剪力与关键点竖向位移Δ0(t)的关系
根据边坡体的刚性假定,离散的各条柱满足一定的协调条件,即各条柱的水平位移相等。
首先选定该三维蠕变边坡的关键点(假设为[1,j0]条柱,其竖向位移为Δ0(t)),那么该三维蠕变边坡水平位移为:
其中为[1,j0]条柱的水平位移,Δh(t)为所有条柱的水平位移(即该边坡的水平位移),为[1,j0]条柱底面剪力的单位向量在x、z方向上的分量。
那么,[i,j]条柱底面的剪切位移Δi,j(t)可以得到:
将式(8)分别代入到式(4)和(6)得到:
(a)当条柱底面的剪力Si,j小于长期抗剪力时,
(b)当条柱底面的剪力Si,j大于长期抗剪力时,
步骤4、根据3个力的平衡方程和3个力矩的平衡方程,并结合[i,j]条柱底面的剪力与关键点竖向位移的关系,建立一个含有未知数α、β、Δ0和蠕变时间t的非线性方程组;求解该方程组,得到不同的时刻所对应关键点的竖向位移Δ0;进而得到该三维蠕变边坡的滑动位移随时间的变化规律。
(1)、确定[i,j]条柱底面的剪力与该条柱的剪切位移之间的关系
[i,j]条柱的受力模型如图3所示,[i,j]条柱底面上法向力方向的单位向量为
式中,f是滑裂面的函数,为条柱底面上法向力方向的单位向量,其中分别是该向量在x、y、z方向上的分量,
条柱底面上剪力方向的单位向量为
式中,为条柱底面剪力的单位向量,其中分别是该向量在x、y、z方向上的分量,
根据三个轴向力的平衡,每个条柱上可建立如下方程:
沿x方向力的平衡方程为:
沿y方向力的平衡方程:
沿z方向力的平衡方程:
分别联立式(13)~(15)与式(9)或(10)得到第[i,j]条柱底面的法向力和剪力与关键点竖向位移的关系。
(1)当条柱底面的剪力Si,j小于长期抗剪力时,
第[i,j]条柱底面的法向力和剪力与关键点竖向位移Δ0的关系为:
(2)当条柱底面的剪力Si,j大于等于长期抗剪力时,
第[i,j]条柱底面的法向力和剪力与关键点竖向位移Δ0的关系为:
式(16)~(21)中,Ni,j和Si,j分别为第[i,j]条柱底面的法向力和剪力,分别是[i,j]条柱底面剪力的单位向量在x、y、z方向上的分量,为[1,j0]条柱底面剪力的单位向量在x、z方向上的分量,分别是[i,j]条柱底面上法向力方向的单位向量在x、y、z方向上的分量,Wi,j是[i,j]条柱的重量。
将法向力Ni,j和剪力Si,j代入到3个坐标轴的力矩平衡方程,力矩的平衡方程如下,
边坡绕x轴的力矩平衡方程:
边坡绕y轴的力矩平衡方程:
边坡绕z轴的力矩平衡方程:
可以建立一组非线性方程组,记为:
通过赋予不同的时刻值t,并求解非线性方程组(25),得到不同的时刻对应关键点竖向位移Δ0。编制计算机程序,代入相关参数,可以得到该实施例的关键点竖向位移Δ0随时间变化的关系。如图5所示,纵坐标为关键点竖向位移Δ0,横坐标为变化的时间,曲线图表示了关键点竖向位移变化量与对应时间的关系。该曲线图比二维蠕变边坡的分析方法所得的蠕变曲线更接近于实际状态的变化。
Claims (3)
1.一种三维蠕变边坡滑动位移的预测方法,包括步骤1、选定具体待预测的三维蠕变边坡,确定该三维蠕变边坡潜在滑裂面的几何尺寸,将潜在滑裂面和边坡表面用方程表示;确定岩土的土体指标参数;
步骤2、将三维蠕变边坡体离散化,三维蠕变边坡被垂直离散为m行和n列条柱,每个条柱按所在的行号i和列号j定义为[i,j];假定行方向的条柱间作用力与水平面夹角为±α,假定列方向的条柱间作用力与水平面夹角为±β;
其特征是,还包括以下步骤:
步骤3、根据西原体模型,确定第[i,j]条柱底面的剪力与该条柱的剪切位移之间的关系;利用离散后的各条柱的位移协调条件,得到第[i,j]条柱底面的剪力与关键点竖向位移Δ0的关系;
步骤4、根据3个力的平衡方程和3个力矩的平衡方程,并结合[i,j]条柱底面的剪力与关键点竖向位移的关系,建立一个含有未知数α、β、Δ0和蠕变时间t的非线性方程组;求解该方程组,得到不同的时刻所对应关键点的竖向位移Δ0;进而得到该三维蠕变边坡的滑动位移随时间的变化规律。
2.根据权利要求1所述的一种三维蠕变边坡滑动位移的预测方法,其特征是:在步骤3中,边坡条柱中的西原体模型分为两个部分:
(1)当条柱底面的剪力Si,j小于长期抗剪力时,得到第[i,j]条柱底面的剪力与该条柱的剪切位移之间的关系为:
式中Δi,j(t)为第[i,j]条柱底面的剪切位移,Si,j为第[i,j]条柱底面的剪力,G1、G2、η1、η2为西原体模型中的参数,t是发生蠕变的时间,b是条柱行方向的宽度;
(2)当条柱底面的剪力Si,j大于等于长期抗剪力时,得到第[i,j]条柱底面的剪力与该条柱的剪切位移之间的关系为:
其中为长期抗剪力, 为长期抗剪强度,Ai,j为第[i,j]条柱的横截面积,为第[i,j]条柱底面法向力的单位向量在z轴上的分量。
3.根据权利要求2所述的一种三维蠕变边坡滑动位移的预测方法,其特征是:在步骤4中,联立三个力的平衡方程和条柱底面的剪力与该条柱的剪切位移之间的关系式,得到第[i,j]条柱底面的法向力和剪力与关键点竖向位移的关系:
(1)当条柱底面的剪力Si,j小于长期抗剪力时,
第[i,j]条柱底面的法向力和剪力与关键点竖向位移Δ0的关系为:
其中
(2)当条柱底面的剪力Si,j大于等于长期抗剪力时,
第[i,j]条柱底面的法向力和剪力与关键点竖向位移Δ0的关系为:
其中
以上各式中,Ni,j和Si,j分别为第[i,j]条柱底面的法向力和剪力,分别是[i,j]条柱底面剪力的单位向量在x、y、z方向上的分量,为[1,j0]条柱底面剪力的单位向量在x、z方向上的分量,分别是[i,j]条柱底面上法向力方向的单位向量在x、y、z方向上的分量,Wi,j是[i,j]条柱的重量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410470310.XA CN104268643B (zh) | 2014-09-16 | 2014-09-16 | 一种三维蠕变边坡滑动位移的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410470310.XA CN104268643B (zh) | 2014-09-16 | 2014-09-16 | 一种三维蠕变边坡滑动位移的预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104268643A CN104268643A (zh) | 2015-01-07 |
CN104268643B true CN104268643B (zh) | 2017-07-28 |
Family
ID=52160163
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410470310.XA Active CN104268643B (zh) | 2014-09-16 | 2014-09-16 | 一种三维蠕变边坡滑动位移的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104268643B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111639439B (zh) * | 2020-06-09 | 2022-01-25 | 四川大学 | 一种考虑长期蠕变和基质吸力的土体安全性计算方法 |
CN111611728B (zh) * | 2020-06-29 | 2023-03-28 | 中国有色金属工业昆明勘察设计研究院有限公司 | 昔格达地层岩质边坡后缘滑裂面位置的稳定性计算方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103163563A (zh) * | 2013-01-31 | 2013-06-19 | 重庆大学 | 一种三维边坡稳定性预测方法 |
CN103485353A (zh) * | 2013-09-24 | 2014-01-01 | 昆明理工大学 | 基于全局最优化的边坡稳定性分析条分法 |
-
2014
- 2014-09-16 CN CN201410470310.XA patent/CN104268643B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103163563A (zh) * | 2013-01-31 | 2013-06-19 | 重庆大学 | 一种三维边坡稳定性预测方法 |
CN103485353A (zh) * | 2013-09-24 | 2014-01-01 | 昆明理工大学 | 基于全局最优化的边坡稳定性分析条分法 |
Non-Patent Citations (4)
Title |
---|
"A three-dimensional slope stability analysis method using the upper bound theorem:Part I:thery and methods";Zuyu Chen 等;《International Journal of Rock Mechanics and Mining Sciences》;20010430;第38卷(第3期);第369-378页 * |
"Stability analysis of embankments and slopes";S.K.SARMA;《Géotechnique》;19730930(第3期);第423-433页 * |
"三维不连续变形分析理论及其在岩质边坡工程中的应用";张国新 等;《岩石力学与工程学报》;20101015;第29卷(第10期);第2116-2126页 * |
"西原体模型下圆形洞室围岩蠕变的解析解";方延强 等;《西部探矿工程》;20080615(第6期);第171-174页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104268643A (zh) | 2015-01-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
He et al. | An innovative approach for gob-side entry retaining in thick coal seam longwall mining | |
Fortsakis et al. | Anisotropic behaviour of stratified rock masses in tunnelling | |
Deng et al. | Investigation of directional hydraulic fracturing based on true tri-axial experiment and finite element modeling | |
CN104458309B (zh) | 一种用于物理模拟实验中的相似材料配比确定方法 | |
Hu et al. | Key parameters of roof cutting of gob-side entry retaining in a deep inclined thick coal seam with hard roof | |
Little et al. | Kinematics of oblique collision and ramping inferred from microstructures and strain in middle crustal rocks, central Southern Alps, New Zealand | |
CN104200100B (zh) | 一种基于滑动位移分析的三维边坡稳定性预测方法 | |
CN103135128A (zh) | 一种地震荷载作用下三维边坡稳定性预测方法 | |
CN104268380B (zh) | 一种三维蠕变边坡的长期稳定性预测方法 | |
CN104268643B (zh) | 一种三维蠕变边坡滑动位移的预测方法 | |
DING et al. | RESISTANCE TO PERMEABILITY DAMAGE AND PILLAR SAFETY OF COAL MINING UNDER GIANT THICK LOOSE SEAM—TAKE KOUZI ̄ DONG MINE AS AN EXAMPLE | |
Xu et al. | Prediction of the height of water-conductive fractured zone under continuous extraction and partial backfill mining method—A case study | |
Dong et al. | Quantitative characterization of single-phase flow through rough-walled fractures with variable apertures | |
Gao et al. | Simulation parameter selection and steady seepage analysis of binary structure slope | |
Ding et al. | Study on the damage process and numerical simulation of tunnel excavation in water-rich soft rock | |
Liu et al. | A method of backfill mining crossing the interchange bridge and application of a ground subsidence prediction model | |
Guo et al. | Research on deformation stability of soft rock slope under excavation based on FLAC3D | |
Shi et al. | The stress evolution of adjacent working faces passing through an abandoned roadway and the damage depth of the floor | |
CN116384267A (zh) | 致密储层压裂水平井的最终可采储量的确定方法及设备 | |
Wang et al. | Study of Overlying Rock Structure and Intensive Pressure Control Technology of Island Longwall Panel in Extra-Thick Coal Seams | |
Kai et al. | Computer simulation of dynamic influence of outburst gas flow on mine ventilation network | |
Qin et al. | Study on the Effect of the Undercut Area on the Movement Law of Overburden Rock Layers in the Block Caving Method | |
He et al. | Reasonable coal pillar width and surrounding rock control of gob-side entry driving in inclined short-distance coal seams | |
Yu et al. | Prediction of instability and mechanism of multi-factor comprehensive action on mine goaf | |
Yang et al. | Simulated prediction of roof water breakout for high-intensity mining under reservoirs in mining areas in Western China |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |