CN113239434A - 一种确定锚杆衰减系数的新方法 - Google Patents

一种确定锚杆衰减系数的新方法 Download PDF

Info

Publication number
CN113239434A
CN113239434A CN202110508941.6A CN202110508941A CN113239434A CN 113239434 A CN113239434 A CN 113239434A CN 202110508941 A CN202110508941 A CN 202110508941A CN 113239434 A CN113239434 A CN 113239434A
Authority
CN
China
Prior art keywords
attenuation coefficient
determining
slope
failure
failure probability
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
Application number
CN202110508941.6A
Other languages
English (en)
Other versions
CN113239434B (zh
Inventor
李东贤
李亮
王上上
林后来
刘志良
裴华富
迟世春
贾宇峰
袁铭徽
杨正权
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Qingdao University of Technology
Original Assignee
Qingdao University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Qingdao University of Technology filed Critical Qingdao University of Technology
Priority to CN202110508941.6A priority Critical patent/CN113239434B/zh
Publication of CN113239434A publication Critical patent/CN113239434A/zh
Application granted granted Critical
Publication of CN113239434B publication Critical patent/CN113239434B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/08Probabilistic or stochastic CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/23Dune restoration or creation; Cliff stabilisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Architecture (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Piles And Underground Anchors (AREA)

Abstract

一种确定锚杆衰减系数的新方法,涉及安全风险评估技术领域,包括如下步骤:步骤1、建立加固边坡数学模型,确定一系列锚杆灌浆拔出面上单位表面摩擦力的衰减系数值;步骤2、在某一衰减系数下,计算该衰减系数对应的加固边坡失效概率Pf;步骤3、反复重复步骤2,获取一系列衰减系数与加固边坡失效概率Pf的关系曲线;步骤4、根据期望失效概率确定锚杆的衰减系数。本发明提供了一种确定锚杆衰减系数的新方法,该方法可以对加固后的边坡进行稳定程度评估,并确定锚杆衰减系数。

Description

一种确定锚杆衰减系数的新方法
技术领域
本发明涉及安全风险评估技术领域,具体涉及一种确定锚杆衰减系数的新方法。
背景技术
在对加筋土体研究过程中,锚杆防护因其施工简单,成本较低,支护效果好等优点,而受到边坡工程界的青睐。
锚杆灌浆拔出面上的单位表面摩擦力与土层情况、及锚杆的安装过程有关,是一个与特定场地相关的参数,另外,受到边坡内渗流水的化学与力学作用,锚杆灌浆拔出面上的单位表面摩擦力会降低,导致锚杆局部失效。
如上所述,由于锚杆荷载传递机理的复杂性及支护参数的多样性,锚杆的不确定性对边坡失效概率的影响研究具有重要的意义,当前对锚杆不确定研究较少,因此加固边坡仍然面临巨大挑战,亟需一种可以对加固后的边坡进行稳定程度评估并确定锚杆衰减系数的新方法。
发明内容
为了解决这些问题,本发明提出了一种确定锚杆衰减系数的新方法。
为达到上述目的,本发明的技术方案为:
一种确定锚杆衰减系数的新方法,包括如下步骤:
步骤1、建立加固边坡数学模型,确定一系列锚杆灌浆拔出面上单位表面摩擦力的衰减系数值;
步骤2、在某一衰减系数下,计算该衰减系数对应的加固边坡失效概率Pf
步骤3、反复重复步骤2,获取一系列衰减系数与加固边坡失效概率Pf的关系曲线;
步骤4、根据期望失效概率确定锚杆的衰减系数。
优选的,所述的步骤1中,通过选定锚杆失效根数,在规定范围0≤Ki≤1内,确定M个锚杆灌浆拔出面上单位表面摩擦力的衰减系数值K1,K2,......KM
优选的,所述的步骤2中,将边坡土体参数黏聚力c和内摩擦角φ视为随机变量来模拟土层参数的不确定性,利用蒙特卡罗抽样方法产生一系列抽样样本,利用极限平衡方法的摩根斯坦法求出对应的安全系数,进而利用边坡安全系数小于某一规定值确定失效样本,汇总计算该衰减系数对应的加固边坡失效概率Pf
优选的,所述的步骤2中,在衰减系数值Ki下,根据锚杆失效根数,在所有施加锚杆中随机选择失效锚杆,基于失效锚杆灌浆拔出面上单位表面摩擦力的初始值F0和Ki确定失效锚杆灌浆拔出面上单位表面摩擦力的计算值F=F0×Ki,将边坡土体参数黏聚力c、内摩擦角φ、视为随机变量,利用Matlab软件产生N个样本Z1,Z2,......ZN,基于极限平衡方法中的摩根斯坦法计算出Z1,Z2,......ZN相应的安全系数Fs1,Fs2,……,FsN,按<<建筑边坡工程技术规范GB50330-20313>>确定最小安全系数P值,若Fsi<P,则称Zi为失效样本,如此重复进行判断,共获取m个失效样本,计算加固边坡失效概率Pfi=m/N。
优选的,所述的步骤3中,重复步骤2计算出K1,K2,......KM对应的加固边坡失效概率,分别记为Pf1,Pf2,Pf3,……,PfM,利用(Ki,Pfi),i=1、2...M,绘制衰减系数与边坡失效概率关系曲线。
优选的,所述的步骤4中,通过改变锚杆失效根数,重复步骤2、3,得到不同锚杆失效根数时,衰减系数与边坡失效概率关系曲线,根据期望失效概率,锚杆失效根数,查询曲线确定锚杆衰减系数。
本发明一种确定锚杆衰减系数的新方法的有益效果为:本发明通过确定不同的锚杆灌浆拔出面上单位表面摩擦力的衰减系数Ki设计值,通过计算获取衰减系数Ki与加固边坡失效概率关系曲线,最终可以根据期望失效概率有效确定锚杆的衰减系数。
附图说明
图1、本发明流程图;
图2、深圳假日广场深基坑边坡支护图;
图3、深圳假日广场地质柱状图;
图4、不同锚杆个数局部失效模型下衰减系数与失效概率变化曲线;
具体实施方式
以下所述,是以阶梯递进的方式对本发明的实施方式详细说明,该说明仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
本发明的描述中,需要说明的是,术语“上”“下”“左”“右”“顶”“底”“内”“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以及特定的方位构造和操作,因此不能理解为对本发明的限制。
一种确定锚杆衰减系数的新方法,如图1所示,包括如下步骤:
步骤1、建立加固边坡数学模型,确定一系列锚杆灌浆拔出面上单位表面摩擦力的衰减系数值;
步骤2、在某一衰减系数下,计算该衰减系数对应的加固边坡失效概率Pf
步骤3、反复重复步骤2,获取一系列衰减系数与加固边坡失效概率Pf的关系曲线;
步骤4、根据期望失效概率确定锚杆的衰减系数;
如图1所示,所述的步骤1中,通过选定锚杆失效根数,在规定范围0≤Ki≤1内,确定M个锚杆灌浆拔出面上单位表面摩擦力的衰减系数值K1,K2,......KM
如图1所示,所述的步骤2中,将边坡土体参数黏聚力c和内摩擦角φ视为随机变量来模拟土层参数的不确定性,利用蒙特卡罗抽样方法产生一系列抽样样本,利用极限平衡方法的摩根斯坦法求出对应的安全系数,进而利用边坡安全系数小于某一规定值确定失效样本,汇总计算该衰减系数对应的加固边坡失效概率Pf
如图1所示,所述的步骤2中,在衰减系数值Ki下,根据锚杆失效根数,在所有施加锚杆中随机选择失效锚杆,基于失效锚杆灌浆拔出面上单位表面摩擦力的初始值F0和Ki确定失效锚杆灌浆拔出面上单位表面摩擦力的计算值F=F0×Ki,将边坡土体参数黏聚力c、内摩擦角φ、视为随机变量,利用Matlab软件产生N个样本Z1,Z2,......ZN,基于极限平衡方法中的摩根斯坦法计算出Z1,Z2,......ZN相应的安全系数Fs1,Fs2,……,FsN,按<<建筑边坡工程技术规范GB50330-20313>>确定最小安全系数P值,若Fsi<P,则称Zi为失效样本,如此重复进行判断,共获取m个失效样本,计算加固边坡失效概率Pfi=m/N;
如图1所示,所述的步骤3中,重复步骤2计算出K1,K2,......KM对应的加固边坡失效概率,分别记为Pf1,Pf2,Pf3,……,PfM,利用(Ki,Pfi),i=1、2...M,绘制衰减系数与边坡失效概率关系曲线;
如图1所示,所述的步骤4中,通过改变锚杆失效根数,重复步骤2、3,得到不同锚杆失效根数时,衰减系数与边坡失效概率关系曲线,根据期望失效概率,锚杆失效根数,查询曲线确定锚杆衰减系数。
下面结合图2进行实例说明:
该边坡坡高14.35m,坡面水平投影长4m,坡角75.5°,边坡坡率为1:0.28,场地地层由人工素填土层、坡洪积层、残积层和燕山期基岩组成,第一层素填土粘聚力c为10kPa、内摩擦角φ为8°、容重γ为1.96kN/m3;第二层粉质黏土粘聚力c为24.9kPa、内摩擦角φ为20°、容重γ为1.96kN/m3,第三层花岗岩粘聚力c为160kPa、内摩擦角φ为35°、容重γ为2.60kN/m3
采用本发明提出的新方法进行计算的步骤如下:
利用Geo-studio软件中Slope/w模块建立边坡数值模型,本工程共施加4根锚杆,锚杆灌浆拔出面上单位表面摩擦力的初始值F0=70kPa,设定衰减系数K1=0,K2=0.1,K3=0.3,K4=0.5,K5=0.7,K6=0.9,K7=1.0,首先选定锚杆失效根数为1,从四根锚杆中随机选择一根作为失效锚杆,利用F0与衰减系数值K1计算失效锚杆灌浆拔出面上单位表面摩擦力的计算值F=70×0=0,将三层边坡土体参数黏聚力c、内摩擦角φ视为服从对数正态分布的随机变量,变异系数均为0.3,重度保持不变,基于随机方法利用Matlab软件产生N=2000个土体参数样本Z1,Z2,……,Z2000,利用极限平衡方法中的摩根斯坦法计算Z1,Z2,……,Z2000对应2000个安全系数Fsi,i=1,2,……,2000;按<<建筑边坡工程技术规范GB50330-20313>>确定P=1.0,若安全系数Fsi<1,则认为Zi为失效样本,共统计失效样本个数为m,则衰减系数K1时边坡的失效概率Pf为m与2000之比;重复计算,可得到不同衰减系数下,加固边坡的失效概率值;通过改变锚杆失效根数为2根、3根,重复上述计算过程,可以得到一根、两根和三根锚杆失效下不同衰减系数对应的失效概率曲线;
如图4所示:
一根锚杆局部失效边坡失效概率分别是:
Pf1=21.4%,Pf2=21.1%,Pf3=17.9%,Pf4=17.8%,Pf5=13.4%,Pf6=13%,Pf7=12.1%;
两根锚杆局部失效边坡失效概率分别是:
Pf1=34.2%,Pf2=31.1%,Pf3=25%,Pf4=23.2%,Pf5=17.2%,Pf6=12.4%,Pf7=12.1%;
三根锚杆局部失效边坡失效概率分别是:
Pf1=56.7%,Pf2=50.8%,Pf3=41.2%,Pf4=31.8%,Pf5=24.6%,Pf6=12.8%,Pf7=12.1%;利用
(Ki,Pfi),i=1,2,3,......7绘制衰减系数与边坡失效概率关系曲线,根据期望失效概率,查询曲线,确定锚杆的衰减系数;譬如当期望失效概率为20%时,一根锚杆失效的衰减系数为0.17,两根锚杆失效的衰减系数为0.58,三根锚杆失效的衰减系数为0.78。
通过以上实施例与现有技术进行对比发现:已有研究成果主要集中于锚固体强度特性和锚杆时变性对锚固系统的破坏概率的影响,而对锚杆不同衰减系数对边坡失效概率的影响研究涉及较少,通过计算获得各个衰减系数下的边坡失效概率,充分验证了衰减系数对边坡失效概率的影响。通过工程实例验证了随着衰减系数的增大,边坡的失效概率逐渐减小,最终可以根据期望失效概率来查询曲线确定衰减系数,通过实例对比分析验证了本发明的有效性。

Claims (6)

1.一种确定锚杆衰减系数的新方法,其特征为:包括如下步骤:
步骤1、建立加固边坡数学模型,确定一系列锚杆灌浆拔出面上单位表面摩擦力的衰减系数值;
步骤2、在某一衰减系数下,计算该衰减系数对应的加固边坡失效概率Pf
步骤3、反复重复步骤2,获取一系列衰减系数与加固边坡失效概率Pf的关系曲线;
步骤4、根据期望失效概率确定锚杆的衰减系数。
2.如权利要求1所述的一种确定锚杆衰减系数的新方法,其特征为:所述的步骤1中,通过选定锚杆失效根数,在规定范围0≤Ki≤1内,确定M个锚杆灌浆拔出面上单位表面摩擦力的衰减系数值K1,K2,......KM
3.如权利要求2所述的一种确定锚杆衰减系数的新方法,其特征为:所述的步骤2中,将边坡土体参数黏聚力c和内摩擦角φ视为随机变量来模拟土层参数的不确定性,利用蒙特卡罗抽样方法产生一系列抽样样本,利用极限平衡方法的摩根斯坦法求出对应的安全系数,进而利用边坡安全系数小于某一规定值确定失效样本,汇总计算该衰减系数对应的加固边坡失效概率Pf
4.如权利要求3所述的一种确定锚杆衰减系数的新方法,其特征为:所述的步骤2中,在衰减系数值Ki下,根据锚杆失效根数,在所有施加锚杆中随机选择失效锚杆,基于失效锚杆灌浆拔出面上单位表面摩擦力的初始值F0和Ki确定失效锚杆灌浆拔出面上单位表面摩擦力的计算值F=F0×Ki,将边坡土体参数黏聚力c、内摩擦角
Figure FDA0003059529780000011
视为随机变量,利用Matlab软件产生N个样本Z1,Z2,......ZN,基于极限平衡方法中的摩根斯坦法计算出Z1,Z2,......ZN相应的安全系数Fs1,Fs2,……,FsN,按<<建筑边坡工程技术规范GB50330-20313>>确定最小安全系数P值,若Fsi<P,则称Zi为失效样本,如此重复进行判断,共获取m个失效样本,计算加固边坡失效概率Pfi=m/N。
5.如权利要求4所述的一种确定锚杆衰减系数的新方法,其特征为:所述的步骤3中,重复步骤2计算出K1,K2,......KM对应的加固边坡失效概率,分别记为Pf1,Pf2,Pf3,……,PfM,利用(Ki,Pfi),i=1、2...M,绘制衰减系数与边坡失效概率关系曲线。
6.如权利要求5所述的一种确定锚杆衰减系数的新方法,其特征为:所述的步骤4中,通过改变锚杆失效根数,重复步骤2、3,得到不同锚杆失效根数时,衰减系数与边坡失效概率关系曲线,根据期望失效概率,锚杆失效根数,查询曲线确定锚杆衰减系数。
CN202110508941.6A 2021-05-11 2021-05-11 一种确定锚杆衰减系数的方法 Active CN113239434B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110508941.6A CN113239434B (zh) 2021-05-11 2021-05-11 一种确定锚杆衰减系数的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110508941.6A CN113239434B (zh) 2021-05-11 2021-05-11 一种确定锚杆衰减系数的方法

Publications (2)

Publication Number Publication Date
CN113239434A true CN113239434A (zh) 2021-08-10
CN113239434B CN113239434B (zh) 2022-07-12

Family

ID=77133277

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110508941.6A Active CN113239434B (zh) 2021-05-11 2021-05-11 一种确定锚杆衰减系数的方法

Country Status (1)

Country Link
CN (1) CN113239434B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104794365A (zh) * 2015-05-06 2015-07-22 南华大学 一种基于数学模型预测锚杆极限承载力的计算方法
CN106226100A (zh) * 2016-07-15 2016-12-14 中南大学 一种边坡预应力锚索使用性能测试结构及测试方法
CN106638718A (zh) * 2016-10-19 2017-05-10 湖南工程学院 全长粘结性土锚脱粘长度及动态锚固力测定方法
CN110532617A (zh) * 2019-07-29 2019-12-03 中国地质大学(武汉) 一种锚杆及锚索的安全性评价方法
CN111985041A (zh) * 2020-09-17 2020-11-24 青岛理工大学 一种加固边坡挡墙高度确定方法及加固边坡挡墙
CN112613143A (zh) * 2020-12-26 2021-04-06 北京工业大学 一种基于时变失效的螺栓可靠性分析方法
CN112685817A (zh) * 2020-12-24 2021-04-20 青岛理工大学 一种对锚固边坡风险进行量化分析的方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104794365A (zh) * 2015-05-06 2015-07-22 南华大学 一种基于数学模型预测锚杆极限承载力的计算方法
CN106226100A (zh) * 2016-07-15 2016-12-14 中南大学 一种边坡预应力锚索使用性能测试结构及测试方法
CN106638718A (zh) * 2016-10-19 2017-05-10 湖南工程学院 全长粘结性土锚脱粘长度及动态锚固力测定方法
CN110532617A (zh) * 2019-07-29 2019-12-03 中国地质大学(武汉) 一种锚杆及锚索的安全性评价方法
CN111985041A (zh) * 2020-09-17 2020-11-24 青岛理工大学 一种加固边坡挡墙高度确定方法及加固边坡挡墙
CN112685817A (zh) * 2020-12-24 2021-04-20 青岛理工大学 一种对锚固边坡风险进行量化分析的方法
CN112613143A (zh) * 2020-12-26 2021-04-06 北京工业大学 一种基于时变失效的螺栓可靠性分析方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
许明 等: ""锚索失效对边坡稳定性影响的数值研究"", 《地下空间与工程学报》 *

Also Published As

Publication number Publication date
CN113239434B (zh) 2022-07-12

Similar Documents

Publication Publication Date Title
Poulos et al. Piled raft foundations for tall buildings
Nadim Tools and strategies for dealing with uncertainty in geotechnics
Evans et al. Three-dimensional simulations of plate anchor pullout in granular materials
CN105547526A (zh) 一种断层构造应力的监测装置和方法
Amani et al. Design of monopiles for offshore and nearshore wind turbines in seismically liquefiable soils: Methodology and validation
Demir et al. Uplift response of multi-plate helical anchors in cohesive soil
Li et al. The loading behavior of innovative monopile foundations for offshore wind turbine based on centrifuge experiments
Vicent et al. Evaluation of horizontal and vertical bearing capacities of offshore bucket work platforms in sand
Shu et al. Bearing capacity and reliability analysis of spudcan foundations embedded at various depths based on the non-stationary random finite element method
CN113239434B (zh) 一种确定锚杆衰减系数的方法
Byrne et al. Seismic liquefaction: centrifuge and numerical modeling
Niroumand et al. Analytical and numerical study of horizontal anchor plates in cohesionless soils
CN115358111A (zh) 一种基于dfm模型的蓄水期库岸裂隙岩体边坡稳定性模拟分析方法、设备及介质
Győri et al. Methods and uncertainties in liquefaction hazard assessment for nuclear power plants
CN115203793A (zh) 一种竖向受荷桩极限承载力预测方法及系统
CN112926118A (zh) 一种变电站工程深基础土-结构协同分析方法
Bhattacharya et al. Seismic design and analysis of offshore wind turbines
Juárez-Luna et al. Application of fracture mechanics to cracking problems in soils
Wu Seismic design of dams
Byrne et al. Numerical modeling of dynamic centrifuge tests
Hu et al. Time-variant seismic reliability analysis of slopes with combined uncertainties based on probability density evolution method
Peterson Development of a Performance-Based Procedure for Assessment of Liquefaction-Induced Free-Field Settlements
EP3445915A1 (en) Method and kit for mitigating the risk of liquefaction of ground to be consolidated
Strassburg Influence of Relative Compaction on Passive Resistance of Abutments with Mechanically Stabilized Earth (MSE) Wingwalls
Al-Kaabi et al. Numerical Modeling Load Displacement Behavior of Screw Piles under Seismic Loading in Soft Soil

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