CN110569609B - 一种边坡失稳后粒子位移临界值的确定方法 - Google Patents
一种边坡失稳后粒子位移临界值的确定方法 Download PDFInfo
- Publication number
- CN110569609B CN110569609B CN201910861265.3A CN201910861265A CN110569609B CN 110569609 B CN110569609 B CN 110569609B CN 201910861265 A CN201910861265 A CN 201910861265A CN 110569609 B CN110569609 B CN 110569609B
- Authority
- CN
- China
- Prior art keywords
- slope
- particle
- group
- equal
- displacement
- 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
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
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
本发明属于边坡稳定和风险分析技术领域,尤其涉及一种边坡失稳后粒子位移临界值的确定方法,包括设计n组粘性土边坡的坡高和坡角;对于第i组边坡,用有限差分强度折减法分析,得到安全系数Fsi和塑性区贯通带;画出塑性区贯通带的中间线,计算出有限差分法下的滑动土方量Ai;用光滑粒子流体动力学方法分析第i组边坡,将得到的安全系数Fsi作为土体强度的折减系数,计算得到第i组粘性土边坡失稳后的粒子位移文件和边坡位移云图;计算第i组粘性土边坡失稳后的粒子位移临界值d*i;令i=i+1,且1≤i≤n,重复步骤2~步骤5,得到所有坡高、坡角相对应的粒子位移临界值将得到的粒子位移临界值与其对应的坡高h、坡角α制作成实用表格,供边坡风险分析用。本发明使得临界值的获得快捷、有效。
Description
技术领域
本发明属于边坡稳定和风险分析技术领域,尤其涉及一种边坡失稳后粒子位移临界值的确定方法。
背景技术
边坡稳定作为土力学三大经典问题之一,确定性因素导致的应力、变形和稳定性分析方法已渐趋于完善和成熟,而分析过程中所包含的不确定性因素致使的风险暴露的愈加明显。所以边坡稳定的可靠度和风险分析成为科研工作者和工程师越发关注的课题。而进行风险分析时边坡失稳后粒子滑动的位移临界值亟需一种方便快捷、行之有效的确定方法。
发明内容
根据以上现有技术的不足,本发明提供了一种边坡失稳后粒子位移临界值的确定方法,其能快捷、有效地确定失稳后粒子的位移临界值,整理数据得到实用表格,并进一步用于边坡稳定的风险分析。
本发明解决的技术问题采用的技术方案为:一种边坡失稳后粒子位移临界值的确定方法,目的在于确定光滑粒子流体动力学方法(SPH)分析粘性土边坡稳定问题时土体粒子的位移临界值,以便应用于风险分析。具体包括步骤:
步骤1:边坡失稳后的滑动土方量只与坡高h和坡角α相关,故按照均分思想合理设计n组粘性土边坡的坡高h1,h2,...,hn和坡角a1,a2,...,an,n为正整数;
步骤2、对于第i组粘性土边坡,1≤i≤n,且i为正整数,使用有限差分强度折减法进行分析,得到其安全系数Fsi和塑性区贯通带;
步骤3、画出塑性区贯通带的中间线,并基于此线计算出第i组粘性土边坡有限差分法下的滑动土方量Ai;
步骤4、用光滑粒子流体动力学方法分析第i组粘性土边坡,将有限差分法得到的安全系数Fsi作为土体强度的折减系数,计算得到第i组粘性土边坡失稳后的粒子位移文件和边坡位移云图;
步骤5、用MATLAB软件计算第i组粘性土边坡失稳后的粒子位移临界值d*i;
用光滑粒子流体动力学方法对边坡进行稳定性分析时,每组边坡的土体将由N个粒子组成,在计算得到的粒子位移文件中,各个粒子的顺序按位移由大到小排列,依次将每个粒子的位移记为d1,d2,…,dj,…,dN;
记参数A*i为光滑粒子流体动力学方法计算得到的第i组边坡的滑动土方量,参数d*i为第i组边坡的粒子位移临界值;
5.1、将A*i的初值赋为0,d*i的初值赋为该组边坡最大的粒子位移,即d*i=d1,令j=1;
5.2、判断d*i≤dj是否成立,其中1≤j≤N,j为正整数;
若成立,则进行步骤5.3;
若不成立,则进行步骤5.4;
5.3、令A*i=A*i+0.04,其中0.04m2为一个粒子的土方量,令j=j+1,跳至步骤5.2继续进行;
5.4、判断|A*i-Ai|/Ai<1%是否成立,即判断两种方法所计算的土方量误差是否在1%以内;
若成立,即输出粒子位移临界值d*i;
若不成立,则令A*i=0,d*i=d*i-0.001,令j=1,因为粒子位移最小精确到千分位,跳至步骤5.2继续进行。
本发明具有以下有益效果:本发明通过对比有限差分强度折减法和光滑粒子流体动力学方法(SPH)计算得到的滑坡土方量,从而确定失稳后的位移临界值,整理数据得到实用表格,进一步用于边坡稳定的风险分析,本发明所述方法使得临界值的获得快捷、有效。
附图说明
图1是本发明所提供实施例的流程框图;
图2是本发明所提供实施例中步骤5的流程图;
图3是本发明所提供实施例中具体算例模型示意图;
图4是本发明方法的验证图;
具体实施方式
下面结合附图对本发明做进一步描述。
实施例一:
如图1~图4所示,本发明所述的一种边坡失稳后粒子位移临界值的确定方法,本发明的目的在于确定光滑粒子流体动力学方法(SPH)分析粘性土边坡稳定问题时土体粒子的位移临界值,以便应用于风险分析。
首先,用有限差分强度折减法计算得到边坡的安全系数以及塑性发展区贯通带,并画出贯通带的中间线,进而计算出边坡失稳后的滑动土方量;
然后,用光滑粒子流体动力学方法(SPH)计算得到边坡失稳后土体粒子的滑动位移和位移云图,并按粒子位移的大小顺序逐个累加粒子的土方量,并与有限差分强度折减法的土方量作对比,得到该工况下的粒子位移临界值,最后,重复以上步骤,计算不同工况下的粒子位移临界值,并整理成实用表格。
具体包括如下步骤:
步骤1:在坡下土层厚度一定的情况下,土体的材料参数变化只会改变拟研究边坡的安全系数,而不会改变滑动破坏面的形状和位置,边坡失稳后的滑动土方量只与坡高h和坡角α相关,故按照均分思想合理设计n组粘性土边坡的坡高h1,h2,...,hn和坡角a1,a2,...,an,n为正整数;
步骤2、对于第i组粘性土边坡,1≤i≤n,且i为正整数,使用有限差分强度折减法进行分析,得到其安全系数Fsi和塑性区贯通带;
步骤3、画出塑性区贯通带的中间线,并基于此线计算出第i组粘性土边坡有限差分法下的滑动土方量Ai;
步骤4、用光滑粒子流体动力学方法(SPH)分析第i组粘性土边坡,将有限差分法得到的安全系数Fsi作为土体强度的折减系数,计算得到第i组粘性土边坡失稳后的粒子位移文件和边坡位移云图;
步骤5、用MATLAB软件计算第i组粘性土边坡失稳后的粒子位移临界值d*i;
用光滑粒子流体动力学方法(SPH)对边坡进行稳定性分析时,每组边坡的土体将由N个粒子组成,在计算得到的粒子位移文件中,各个粒子的顺序按位移由大到小排列,依次将每个粒子的位移记为d1,d2,…,dj,…,dN;
记参数A*i为光滑粒子流体动力学方法(SPH)计算得到的第i组边坡的滑动土方量,参数d*i为第i组边坡的粒子位移临界值;
步骤5.1、将A*i的初值赋为0,d*i的初值赋为该组边坡最大的粒子位移,即d*i=d1,令j=1,因为粒子位移文件中,各个粒子的顺序按位移由大到小排列,依次将每个粒子的位移记为d1,d2,…,dN,最大的粒子位移即为d1;
步骤5.2、判断d*i≤dj是否成立,其中1≤j≤N,j为正整数;
若成立,则进行步骤5.3;
若不成立,则进行步骤5.4;
步骤5.3、令A*i=A*i+0.04,其中0.04m2为一个粒子的土方量,令j=j+1,跳至步骤5.2继续进行;
步骤5.4、判断|A*i-Ai|/Ai<1%是否成立,即判断两种方法所计算的土方量误差是否在1%以内;
若成立,即输出粒子位移临界值d*i;
若不成立,则令A*i=0,d*i=d*i-0.001,j=1,因为粒子位移最小精确到千分位,跳至步骤5.2继续进行。
根据表1,以一实际算例说明本发明的可行性。
首先,由极限平衡法易知该边坡安全系数为1.204,是一稳定边坡,并可以得到其滑动破坏面。
然后查表1可知,该边坡失稳后其粒子位移临界值为0.7363m,图4为本发明用于该算例的验证图,可以看出大于位移临界值的粒子组成的滑动土体与极限平衡滑动面贴合得十分一致,证明了本发明是一种行之有效的方法。
表1、边坡失稳后粒子位移临界值的实用表(单位:m)
以上所述为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书以及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (1)
1.一种边坡失稳后粒子位移临界值的确定方法,其特征在于,包括如下步骤:
步骤1:边坡失稳后的滑动土方量只与坡高h和坡角α相关,故按照均分思想合理设计n组粘性土边坡的坡高h1,h2,...,hn和坡角a1,a2,...,an,n为正整数;
步骤2、对于第i组粘性土边坡,1≤i≤n,且i为正整数,使用有限差分强度折减法进行分析,得到其安全系数Fsi和塑性区贯通带;
步骤3、画出塑性区贯通带的中间线,并基于此线计算出第i组粘性土边坡有限差分法下的滑动土方量Ai;
步骤4、用光滑粒子流体动力学方法分析第i组粘性土边坡,将有限差分法得到的安全系数Fsi作为土体强度的折减系数,计算得到第i组粘性土边坡失稳后的粒子位移文件和边坡位移云图;
步骤5、用MATLAB软件计算第i组粘性土边坡失稳后的粒子位移临界值d*i;
用光滑粒子流体动力学方法对边坡进行稳定性分析时,每组边坡的土体将由N个粒子组成,N为正整数,在计算得到的粒子位移文件中,各个粒子的顺序按位移由大到小排列,依次将每个粒子的位移记为d1,d2,…,dj,…,dN;
记参数A*i为光滑粒子流体动力学方法计算得到的第i组边坡的滑动土方量,参数d*i为第i组边坡的粒子位移临界值;
步骤如下:
步骤5.1、将A*i的初值赋为0,d*i的初值赋为该组边坡最大的粒子位移,即d*i=d1,令j=1;
步骤5.2、判断d*i≤dj是否成立,其中1≤j≤N,j为正整数;
若成立,则进行步骤5.3;
若不成立,则进行步骤5.4;
步骤5.3、令A*i=A*i+0.04,其中0.04m2为一个粒子的土方量,令j=j+1,跳至步骤5.2继续进行;
步骤5.4、判断|A*i-Ai|/Ai<1%是否成立,即判断两种方法所计算的土方量误差是否在1%以内;
若成立,即输出粒子位移临界值d*i;
若不成立,则令A*i=0,d*i=d*i-0.001,j=1,跳至步骤5.2继续进行;
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910861265.3A CN110569609B (zh) | 2019-09-12 | 2019-09-12 | 一种边坡失稳后粒子位移临界值的确定方法 |
PCT/CN2019/130226 WO2021047121A1 (zh) | 2019-09-12 | 2019-12-31 | 一种边坡失稳后粒子位移临界值的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910861265.3A CN110569609B (zh) | 2019-09-12 | 2019-09-12 | 一种边坡失稳后粒子位移临界值的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110569609A CN110569609A (zh) | 2019-12-13 |
CN110569609B true CN110569609B (zh) | 2020-04-17 |
Family
ID=68779429
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910861265.3A Active CN110569609B (zh) | 2019-09-12 | 2019-09-12 | 一种边坡失稳后粒子位移临界值的确定方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN110569609B (zh) |
WO (1) | WO2021047121A1 (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110569609B (zh) * | 2019-09-12 | 2020-04-17 | 青岛理工大学 | 一种边坡失稳后粒子位移临界值的确定方法 |
CN111310392B (zh) * | 2020-03-04 | 2021-04-27 | 青岛理工大学 | 一种用于评估基坑开挖边坡失稳面积放大效应的方法 |
CN112016224B (zh) * | 2020-07-28 | 2022-11-18 | 西南大学 | 基于sph的土坡滑面分析判定方法、系统、终端及介质 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101547070B1 (ko) * | 2014-12-05 | 2015-08-26 | 연세대학교 산학협력단 | 선행 강우 효과를 고려한 비탈면 해석 방법 및 시스템 |
JP2017116971A (ja) * | 2015-12-21 | 2017-06-29 | 公益財団法人鉄道総合技術研究所 | 降雨による斜面崩壊と地震の複合災害に対する運転規制の基準値変更システム |
CN106227904B (zh) * | 2016-02-06 | 2019-05-24 | 湖北工业大学 | 一种基于边坡力分布特征的防护措施设计方法 |
CN108334719B (zh) * | 2018-03-26 | 2021-08-17 | 四川理工学院 | 一种基于sph法的土质边坡稳定性及滑坡运动过程分析方法 |
CN109359361B (zh) * | 2018-09-30 | 2019-06-28 | 青岛理工大学 | 边坡失稳后果量化分析方法 |
CN109457739B (zh) * | 2018-11-08 | 2019-06-18 | 青岛理工大学 | 一种基于下游构筑物受损程度的边坡安全度评价方法 |
CN109977554B (zh) * | 2019-03-28 | 2020-06-30 | 青岛理工大学 | 一种边坡滑动面积的评估方法 |
CN110210175A (zh) * | 2019-06-21 | 2019-09-06 | 贵州正业工程技术投资有限公司 | 一种复合地基填方边坡有限元稳定性计算方法 |
CN110569609B (zh) * | 2019-09-12 | 2020-04-17 | 青岛理工大学 | 一种边坡失稳后粒子位移临界值的确定方法 |
-
2019
- 2019-09-12 CN CN201910861265.3A patent/CN110569609B/zh active Active
- 2019-12-31 WO PCT/CN2019/130226 patent/WO2021047121A1/zh active Application Filing
Also Published As
Publication number | Publication date |
---|---|
CN110569609A (zh) | 2019-12-13 |
WO2021047121A1 (zh) | 2021-03-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110569609B (zh) | 一种边坡失稳后粒子位移临界值的确定方法 | |
CN104899380B (zh) | 一种基于蒙特卡洛模拟的边坡稳定可靠度敏感性分析方法 | |
CN109165425B (zh) | 一种齿轮接触疲劳可靠性分析方法 | |
Alizadeh et al. | A comparison of two and three-dimensional analyses of fatigue crack closure | |
CN111310392A (zh) | 一种用于评估基坑开挖边坡失稳面积放大效应的方法 | |
CN109359361B (zh) | 边坡失稳后果量化分析方法 | |
CN111931369B (zh) | 降雨型滑坡稳定性分析及运动距离测算方法、设备及介质 | |
CN111914457B (zh) | 一种输电线路塔基边坡稳定性判别方法、装置和存储介质 | |
CN105603955B (zh) | 一种基于实测数据的剩余沉降计算方法及应用 | |
Erami et al. | Numerical solution of bed load transport equations using discrete least squares meshless (DLSM) method | |
CN113338904B (zh) | 一种海上油井产能评价方法 | |
CN112464584A (zh) | 自由表面流的水位和流量推求方法 | |
CN111259473B (zh) | 一种基于混合极限平衡法的涉水岸坡安全系数计算方法 | |
CN111104746B (zh) | 一种基于小波分析的河漫滩高程确定方法 | |
Ndengna et al. | A Novel Sediment Transport Model (STM) Accounting Phase Lag Effect. A Resonance Condition | |
CN111931272B (zh) | 任意尺度边坡等精度安全系数计算方法及网格划分方法 | |
CN113591407A (zh) | 一种基于变维转化算法的流固区域的流体动力学模型构建及分析方法、系统和可读存储介质 | |
CN113360983A (zh) | 一种边坡可靠度分析与风险评估方法 | |
CN108760119B (zh) | 一种考虑粒径尺寸效应的岩堆主动土压力作用点修正方法 | |
CN104089859A (zh) | 一种支撑剂有效沉降粒径的测量方法 | |
CN111651906B (zh) | 一种高效的大变形尾矿固结渗透性能求解方法 | |
CN104021280B (zh) | 一种适用于尾粉砂发生管涌的临界水力梯度的计算方法 | |
CN107090798A (zh) | 一种水库超饱和输沙状态下恢复饱和系数计算方法 | |
CN113947043A (zh) | 基于cfd模拟的海底管道的振动频率监测方法和系统 | |
CN112113890B (zh) | 测定高含水量沉积物从悬浮到固结状态渗透系数的方法 |
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 |