CN106951504A - 基于固有振动频率动态分析边坡稳定的计算方法 - Google Patents
基于固有振动频率动态分析边坡稳定的计算方法 Download PDFInfo
- Publication number
- CN106951504A CN106951504A CN201710157111.7A CN201710157111A CN106951504A CN 106951504 A CN106951504 A CN 106951504A CN 201710157111 A CN201710157111 A CN 201710157111A CN 106951504 A CN106951504 A CN 106951504A
- Authority
- CN
- China
- Prior art keywords
- slope
- eigentone
- initial
- dynamic
- computational methods
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明公开了一种基于固有振动频率动态分析边坡稳定的计算方法,涉及岩土工程监测领域,用于真实地反映区域的边坡稳定性。其中,该计算方法包括:步骤一、选择边坡重点区域作为观测点;步骤二、根据选定的观测点建立三维分析模型;步骤三、计算该观测点的初始固有振动频率;步骤四、根据计算出的初始固有振动频率,得到动态黏聚力;步骤五、根据动态黏聚力,以及三维分析模型计算出该边坡重点区域的安全系数;步骤六、根据计算出的安全系数评判该重点区域的边坡稳定性。本发明用于对评价边坡的稳定性。
Description
技术领域
本发明涉及岩土工程监测领域,尤其涉及一种基于固有振动频率动态分析边坡稳定的计算方法。
背景技术
边坡稳定性问题作为土力学三大经典问题之一,其稳定分析的方法一直是科学研究者和工程实践者研究的重要领域。
边坡稳定分析的方法很多,不同边坡都有与之对应的力学模型和分析方法,从早期土质边坡的圆弧分析法,到岩质边坡的Sarma法,从连续介质力学方法到基于渗流等的耦合分析,分析的精度和时效性都得到了大幅提高。
其实,边坡实际是一个充满复杂性的系统,以某一局部或某一时间段内测量的相对静态的模型参数来反应整体或全周期下滑坡的相对动态的模态参量,必然会有其局限性和不完善性。
国外研究表明,滑坡的失稳往往是脆性拉伸破坏的同时,也伴随着强度的实时退化。
因此,对边坡强度的实时折减,是数值分析应用实践的瓶颈所在。无论从理论研究还是从生产实践方面,对边坡稳定性模态参量的动态调整是必要的,越来越多的边坡工程提出建立基于实际工程指标的动态稳定分析的要求。
现有技术中,由于模型参数具有获取困难,静态且偏于保守的特点,其理论解在计算滑坡安全性方面过于保守,经常出现边坡安全系数小于1,实际却是稳定的现象。
发明内容
本发明的目的在于提供一种基于固有振动频率动态分析边坡稳定的计算方法,用于真实地反映区域的边坡稳定性。
为达到上述目的,本发明采用如下技术方案:
本发明提供一种基于固有振动频率动态分析边坡稳定的计算方法,该基于固有振动频率动态分析边坡稳定的计算方法包括:
步骤一、选择边坡重点区域作为观测点。
步骤二、根据选定的观测点建立三维分析模型。
步骤三、计算该观测点的初始固有振动频率。
步骤四、根据计算出的初始固有振动频率,得到动态黏聚力。
步骤五、根据动态黏聚力,以及三维分析模型计算出该边坡重点区域的安全系数。
步骤六、根据计算出的安全系数评判该重点区域的边坡稳定性。
步骤三中,根据公式(1)和公式(2)计算初始固有振动频率;
式中,M为试块重量(单位为N),L为三维分析模型的型心到原点O的距离(单位为mm),μ为黏结系数(单位为N/mm2),S为黏结宽度(单位为mm),l为黏结长度(单位为mm),θ为转角(单位为°),fo为初始固有振动频率(单位为Hz)。
根据初始黏聚力c0和初始振动频率f0,由公式(3)计算得到动态黏聚力c:
式中,c0为初始黏聚力(单位为kPa),f0为初始黏聚力对应的初始固有振动频率(Hz)。
步骤六中,根据计算出的安全系数评判该重点区域的边坡稳定性方法为:若安全系数小于或者等于1则该边坡不稳定;若安全系数大于1则表明该边坡稳定。
三维分析模型为基于GIS系统的三维分析模型。
根据上述基于固有振动频率动态分析边坡稳定的计算方法得到的安全系数,可以真实的反映该区域的边坡稳定性。另外,该方法模型参数获取较简单,从而减少了节省了计算时间。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例中计算方法流程图;
图2本发明实施例中边坡振动历史曲线和其振动速度谱示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例一
本实施例提供一种基于固有振动频率动态分析边坡稳定的计算方法,该基于固有振动频率动态分析边坡稳定的计算方法包括:
步骤一、选择边坡重点区域作为观测点。
步骤二、根据选定的观测点建立三维分析模型。
步骤三、计算该观测点的初始固有振动频率。
步骤四、根据计算出的初始固有振动频率,得到动态黏聚力。
步骤五、根据动态黏聚力,以及三维分析模型计算出该边坡重点区域的安全系数。
步骤六、根据计算出的安全系数评判该重点区域的边坡稳定性。
根据上述基于固有振动频率动态分析边坡稳定的计算方法可真实地计算出安全系数,从而真实地反映该区域的边坡稳定性。
进一步地,对于固有振动频率的获得方式,可从相应的设备上测出相应的值,或者根据公式相应的计算出相应的值。
步骤三中,根据公式(1)和公式(2)计算初始固有振动频率;
式中,M为试块重量(单位为N),L为三维分析模型的型心到原点O的距离(单位为mm),μ为黏结系数(单位为N/mm2),S为黏结宽度(单位为mm),l为黏结长度(单位为mm),θ为转角(单位为°),fo为初始固有振动频率(单位为Hz)。
根据初始黏聚力c0和初始振动频率f0,由公式(3)计算得到动态黏聚力c:
式中,c0为初始黏聚力(单位为kPa),f0为初始黏聚力对应的初始固有振动频率(Hz)。
步骤六中,根据计算出的安全系数评判该重点区域的边坡稳定性方法为:若安全系数小于或者等于1则该边坡不稳定;若安全系数大于1则表明该边坡稳定。
三维分析模型为基于GIS系统的三维分析模型。
实施例二
根据上述实施例一基于固有振动频率动态分析边坡稳定的计算方法,下面一具体的例子为例,详细的说明实施例一的计算过程。
研究区域位于日本长崎县西北向南东倾斜的小型斜坡上。研究时间为6月21日~7月1日,地区总降雨量为162.7mm。设备于雨季前后分别监测边坡的振动状况,动态分析边坡稳定性。具体的实施步骤如下:
第一步,选择边坡重点区域观测点。
第二步,建立基于GIS系统的三维分析模型。
第二步,计算该观测点的初始固有振动频率,本实施例中对于初始固有振动频率的计算,通过激光测振仪得出。参见图2,图2为激光测振仪分别在3月16日和7月16日监测的边坡表面振动情况。左侧分别为两次时间段测得的振动速度历时曲线,右侧为经过变换后的振动速度谱。如图可知,3月16日边坡初始固有振动频率为10.74Hz,而经过暴雨之后,7月16日的初始固有振动频率降低为8.789Hz。
第四步,通过初始固有振动频率,修正滑移面的粘聚力,进一步的得到动态黏聚力。计算所需参数见表1所示。由于大雨对潜在滑移面的影响最大,因此只对不易测得的潜在滑移面的黏聚力进行折减。
表1稳定性计算参数表
第五步:根据基于GIS系统的三维分析模型和动态黏聚力计算出安全系数。
在现场测得的初始固有振动频率有效修正潜在滑移面上的初始黏聚力之后,根据基于GIS系统的三维分析模型和修正后的黏聚力(也就是动态黏聚力)分别用以下4种模型计算其安全系数,如表2所示。
表2四种不同三维分析模型得到的安全系数
第六步,根据表2得到的安全系数,判断该区域的边坡稳定性。上述表2中,安全系数均大于1,则证明边坡基本稳定。
另外,通过对该边坡进行常规二维分析,得到的安全系数小于1,现场位移测量结果显示:边坡在降雨末期,部分监测点位移数据有所抬升,分别升高3~5mm不等,然而其他监测设备均未出现异常变化。故现场认定该滑坡依然处于偏稳定状态。
通过分析上述例子,基于固定振动频率调整的GIS系统的三维分析模型的安全分析则更加贴近实际。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。
Claims (5)
1.一种基于固有振动频率动态分析边坡稳定的计算方法,其特征在于,该计算方法包括:
步骤一、选择边坡重点区域作为观测点;
步骤二、根据选定的观测点建立三维分析模型;
步骤三、计算该观测点的初始固有振动频率;
步骤四、根据计算出的初始固有振动频率,得到动态黏聚力;
步骤五、根据动态黏聚力,以及三维分析模型计算出该边坡重点区域的安全系数;
步骤六、根据计算出的安全系数评判该重点区域的边坡稳定性。
2.根据权利要求1所述的基于固有振动频率动态分析边坡稳定的计算方法,其特征在于,步骤三中,根据公式(1)和公式(2)计算初始固有振动频率;
式中,M为试块重量(单位为N),L为建模的模心到原点O的距离(单位为mm),μ为黏结系数(单位为N/mm2),S为黏结宽度(单位为mm),l为黏结长度(单位为mm),θ为转角(单位为°),fo为初始固有振动频率(单位为Hz)。
3.根据权利要求2所述的基于固有振动频率动态分析边坡稳定的计算方法,其特征在于,根据初始黏聚力c0和初始振动频率f0,由公式(3)计算得到动态黏聚力c:
式中,c0为初始黏聚力(单位为kPa),f0为初始黏聚力对应的初始固有振动频率(Hz)。
4.根据权利要求1所述的基于固有振动频率动态分析边坡稳定的计算方法,其特征在于,步骤六中,根据计算出的安全系数评判该重点区域的边坡稳定性方法为:若安全系数小于或者等于1则该边坡不稳定;若安全系数大于1则表明该边坡稳定。
5.根据权利要求1所述的基于固有振动频率动态分析边坡稳定的计算方法,其特征在于,三维分析模型为基于GIS系统的三维分析模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710157111.7A CN106951504A (zh) | 2017-03-16 | 2017-03-16 | 基于固有振动频率动态分析边坡稳定的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710157111.7A CN106951504A (zh) | 2017-03-16 | 2017-03-16 | 基于固有振动频率动态分析边坡稳定的计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106951504A true CN106951504A (zh) | 2017-07-14 |
Family
ID=59473569
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710157111.7A Pending CN106951504A (zh) | 2017-03-16 | 2017-03-16 | 基于固有振动频率动态分析边坡稳定的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106951504A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109001021A (zh) * | 2018-07-13 | 2018-12-14 | 中国水利水电科学研究院 | 采用离心模型试验评价岩土边坡安全储备的方法 |
CN109060963A (zh) * | 2018-06-13 | 2018-12-21 | 北京北科安地科技发展有限公司 | 一种基于振动特征参数的崩塌预测方法、系统及装置 |
CN110595598A (zh) * | 2019-08-09 | 2019-12-20 | 华北水利水电大学 | 基于多普勒远程激光测振技术的边坡孤石稳定性监测预警方法 |
CN113392516A (zh) * | 2021-06-04 | 2021-09-14 | 北京科技大学 | 一种基于主频带重心频率的重力坝抗滑稳定性监测评价方法 |
CN114137601A (zh) * | 2021-11-26 | 2022-03-04 | 华北水利水电大学 | 一种楔体稳定性评价方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1598857A (zh) * | 2004-09-21 | 2005-03-23 | 北京科技大学 | 一种优化露天矿山边坡的方法 |
JP2014085229A (ja) * | 2012-10-24 | 2014-05-12 | West Japan Railway Co | 振動計測に基づく斜面上転石の落石危険度評価方法 |
-
2017
- 2017-03-16 CN CN201710157111.7A patent/CN106951504A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1598857A (zh) * | 2004-09-21 | 2005-03-23 | 北京科技大学 | 一种优化露天矿山边坡的方法 |
JP2014085229A (ja) * | 2012-10-24 | 2014-05-12 | West Japan Railway Co | 振動計測に基づく斜面上転石の落石危険度評価方法 |
Non-Patent Citations (2)
Title |
---|
杜岩: ""基于固有振动频率的危岩块体稳定评价模型研究"", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
杜岩等: ""基于模态参量变化的边坡动态稳定分析新方法"", 《岩土工程学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109060963A (zh) * | 2018-06-13 | 2018-12-21 | 北京北科安地科技发展有限公司 | 一种基于振动特征参数的崩塌预测方法、系统及装置 |
CN109001021A (zh) * | 2018-07-13 | 2018-12-14 | 中国水利水电科学研究院 | 采用离心模型试验评价岩土边坡安全储备的方法 |
CN110595598A (zh) * | 2019-08-09 | 2019-12-20 | 华北水利水电大学 | 基于多普勒远程激光测振技术的边坡孤石稳定性监测预警方法 |
CN113392516A (zh) * | 2021-06-04 | 2021-09-14 | 北京科技大学 | 一种基于主频带重心频率的重力坝抗滑稳定性监测评价方法 |
CN113392516B (zh) * | 2021-06-04 | 2022-01-25 | 北京科技大学 | 一种基于主频带重心频率的重力坝抗滑稳定性监测评价方法 |
CN114137601A (zh) * | 2021-11-26 | 2022-03-04 | 华北水利水电大学 | 一种楔体稳定性评价方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106951504A (zh) | 基于固有振动频率动态分析边坡稳定的计算方法 | |
CN103455702A (zh) | 一种确定区域电离层延迟的方法 | |
Niu et al. | Analysis of the dynamic characteristics of a suspension bridge based on RTK-GNSS measurement combining EEMD and a wavelet packet technique | |
Zhang et al. | A review on roller compaction quality control and assurance methods for earthwork in five application scenarios | |
CN107943823A (zh) | 数字化施工实时大数据处理平台 | |
Depeng et al. | Monitoring slope deformation using a 3-D laser image scanning system: a case study | |
Pagliaroli et al. | Numerical and experimental evaluation of site effects at ridges characterized by complex geological setting | |
Sengani et al. | Influence of rainfall intensity on the stability of unsaturated soil slope: Case Study of R523 road in Thulamela Municipality, Limpopo Province, South Africa | |
CN115808711A (zh) | 基于烈度速报的地震灾害损失评估方法及系统 | |
Rafiei et al. | Aquifer and land subsidence interaction assessment using sentinel-1 data and DInSAR technique | |
Jia et al. | An Improved Low‐Cost Continuous Compaction Detection Method for the Construction of Asphalt Pavement | |
Sui et al. | Application of 3D laser scanner for monitoring of landslide hazards | |
CN116863029A (zh) | 一种基于多源遥感数据的地震灾害应急制图方法及装置 | |
CN106872081A (zh) | 基于自振频率的应力监测预警方法 | |
Alonso-Pandavenes et al. | Early investigation of a landslide sliding surface by HVSR and VES geophysical techniques combined, a case study in Guarumales (Ecuador) | |
CN113204855B (zh) | 地层缩短速率确定方法、系统、设备及可读介质 | |
CN113092046A (zh) | 高陡边坡在地震及降雨作用下的稳定性研究系统 | |
Cui et al. | Extraction of vegetation fraction based on the dimidiate pixel model and vegetation index transform plan | |
CN109270594A (zh) | 一种台站气象数据分区空间拼接拟合方法 | |
Larroudé et al. | Three simple indicators of vulnerability to climate change on a Mediterranean beach: A modeling approach | |
Chuchro et al. | Reducing flood risk using computer system for monitoring river embankments | |
Young et al. | Initialization and Setup of the Coastal Model Test Bed: CSHORE | |
Chen et al. | GNSS for real-time monitoring of bridge dynamic responses | |
Baldi et al. | Geodetic networks for crustal movements studies | |
Tabyaoui et al. | HYDRAULIC MODELLING USING HEC-RAS AND GIS, APPLICATION TO OUED INAOUEN (TAZA, NORTHERN MOROCCO) |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170714 |
|
WD01 | Invention patent application deemed withdrawn after publication |