滑坡预报与有限元实例分析对比研究 (1).pdf
《滑坡预报与有限元实例分析对比研究 (1).pdf》由会员分享,可在线阅读,更多相关《滑坡预报与有限元实例分析对比研究 (1).pdf(7页珍藏版)》请在咨信网上搜索。
1、收稿日期:2023-03-10;修订日期:2023-03-29作者简介:曹 伟(1981),男,硕士研究生,工程师,主要从事岩土工程设计及地质灾害研究。第 41 卷 第 4 期2023 年 8 月江 西 科 学JIANGXI SCIENCEVol.41 No.4Aug.2023 doi:10.13990/j.issn1001-3679.2023.04.032滑坡预报与有限元实例分析对比研究曹 伟,陈仲如(国家林业和草原局西南调查规划院,650031,昆明)摘要:为对滑坡的分析预报工作总结一套行之有效的方法,以三峡库区内滑坡为例展开研究。采用基于有限元原理的 ABAQUS 软件及摩尔库伦本构模型
2、,建立三维地质模型并进行各工况下数值模拟计算。运用MATLAB 自编程序滑坡预测集成系统对滑坡既有监测数据进行分析与预测。最后,将数值模拟计算结果、实际监测数据及预测结果进行对比分析,结果显示该分析预报方法具有较好的可靠性和实用性。对滑坡稳定性及流固耦合分析具有积极作用,对滑坡预测及其他领域的预测预报研究也具有广泛的参考意义。关键词:涉水滑坡;监测预报;有限元模拟;稳定性分析中图分类号:P642 文献标识码:A 文章编号:1001-3679(2023)04-794-07Comparative Study of Landslide Forecasting and Finite Element C
3、ase AnalysisCAO Wei,CHEN Zhongru(Southwest Survey and Planning Institute of National Forestry and Grassland Administration,650031,Kunming,PRC)Abstract:In order to summarize a set of effective methods for the analysis and forecasting of land-slides,landslide in the Three Gorges Reservoir area is take
4、n as an example.Based on the finite element theory of ABAQUS software and the Moll Coulomb constitutive model are used.The three-dimensional geological model is established,and the numerical simulation is carried out under vari-ous working conditions.Then,the landslide prediction integrated system p
5、rogrammed by MATLAB is used to analyze and forecast landslide.Finally,the numerical simulation results,actual moni-toring data and prediction results are compared and analyzed.The analysis and prediction method is proved to be reliable and practical by these conclusions.This paper plays an active ro
6、le in landslide stability and fluid-structure coupling analysis.At the same time,this paper also has extensive refer-ence significance for landslide prediction and prediction research in other fields.Key words:fording landslide;monitoring and prediction;finite element simulation;stability analysis0
7、引言目前,如何对滑坡进行准确的预测预报是世界性难题。根据我国目前大多数已运营水库滑坡的发育情况来看,水库高水位运营期间,特别是蓄水期及降落期间极大概率发生大量的滑坡。三峡库区现处高水位运行期,极可能发生大量的滑坡,这些滑坡一旦下滑入江,除了对崩滑体上的居民造成危害之外,入江涌浪将产生派生灾难。因此,对复活型崩塌滑坡预测预报开展深入研究,是十分必要的,且具有重要的工程意义。1988 年,陈明东、王兰生1两位教授合作,首次在滑坡位移-时间曲线拟合中将灰色系统理论的 GM(1,1)模型纳入,并将拟合结果深入后推,进而提出了一个利用滤波灰色系统分析法模型来实现 对 滑 坡 的 中 期 预 报 的 方
8、法。廖 小 平2(1994)以弹塑性力学为理论基础提出了用于滑坡预报的功率模型。王尚庆3提出了一套完整的预报模型,包括安全系数、可靠性、变形速率和综合信息,且在三峡库区的实践中,更加系统地提出了一种基于滑坡变形破坏迹象、诱发因素等的宏观和综合预报方法。过去几十年,虽然我国在滑坡监测预报方面取得了显著进展,但大多数滑坡仍然发生在未被发现的地区,这一现象令人担忧4-5。本文以三峡库区内典型滑坡为例,筛选出适合该滑坡的预报模型,建立变形破坏的数学模型;综合分析数学模型和数值模拟试验的结果,并与实际变形破坏状况及滑坡现场监测资料验证分析,建立滑坡的预报模型和预报判据。1 预测模型及判据1.1 预测模型
9、简介回归分析6就是通过对观察对象数据的统计分析和处理,研究认识和确定事物之间内在关系及形式的一种有效方法。按照自变量和因变量之间的关系类型,回归分析可分为线性回归分析和非线性回归分析。通过非线性回归技术,构建出滑坡位移-时间关系数学模型,从而更好地探索和预测滑坡未来的变形特征。1.2 预测判据滑坡时间预报和临界判据7主要系指滑坡出现快速变形滑动或滑坡开始时的征兆、状态及主要诱发滑坡因素,统计和研究分析提出滑坡宏观破裂长度值和宽度值以及滑坡边界长度之间的定量关系值作为判定滑坡临界变形滑动的依据。具体归纳概括出有下列几种时间预报的临界判据。1)滑坡的总体稳定性系数等于 1,仅局部小于 1;2)在滑
10、坡前缘区域分布有带状的地下泉水涌出区;3)滑坡体裂隙带处于扩张变形状态,后缘裂隙带宽度应超过滑坡体长度的 1/2;4)滑坡蠕滑移动的速率持续增大,且变形显著,或滑动平均速率明显超过 10 cm/月,且受季节变化影响显著;5)雨季或单点暴雨,滑动面处饱和状态,且地表水冲蚀处出现明显崩塌、滑塌,或局部破裂变形加剧滑动、塌陷变形;6)地震动峰值加速度为 0.15 g,或滑坡场地地震烈度达到;7)滑坡周界轨迹连续不间断,轨迹分维值在1.401.60 区间。从以上滑坡的评判标准8和确定原则来看,要建立一种恒定的、又普遍适用的评判标准是不大可能。滑坡变形是其自身稳定状态最为直观的表征,加上位移监测信息相对
11、容易获得,故位移变化信息也可用来监测。以与位移变形相关的数据信息作为判据包括以下几个因素:位移速率(mm/d)、最大位移量(mm)、位移加速度(mm/d2)、变形曲线切线角和位移矢量角变化。本文通过采用 MATLAB 自编应用程序作为位移数据信息的整理、预处理和预测分析工具,自动计算生成预测位移、位移速率、位移加速度和变形曲线切线角等位移判据信息的数据、图像文件。软件主界面如图 1 所示。清除图导入监测表格线性图柱形图导出表格和图片批量导出更新网格查看监测表格开始预测生成Word文档批量生成加密曲线回归分析模型时间位移图位移/mm02004006008001 000 1 200 1 400 1
12、 600 1 800 2 000050100150200250观测期数图 1 预测分析集成系统主界面2 流固耦合求解流固耦合的方法一般主要有如下 2种9。1)解析法:由于岩土体材料在天然状态下很难运用解析法对流固耦合问题求得解析解。因此,流固耦合问题通常不采用该方法求解。2)数值分析:针对流固耦合问题的有限元数值分析方法一般主流方法包括有限差分法和有限元法。总的来说,从数学理论这一方向来看,有限差分法的近似化程度要比有限元法的高,但若从597第 4 期 曹 伟等:滑坡预报与有限元实例分析对比研究计算实际和应用层面来看,有限元法往往要更简便、灵活。因此,本文选取有限元法计算模拟。3 工程实例3.
13、1 工程概述滑坡位于秭归县台子湾村,地处长江右岸,距离三峡大坝坝址约 53 km,地理坐标为:经度1103331.8,纬度 310031.1。滑坡位于百福坪背斜北侧,地层为三迭系中统巴东组泥岩、灰岩及砂岩,地层产状 2532,自然斜坡走向北西西,位于长江凹岸,滑坡全貌如图 2 所示。图 2 滑坡全貌图滑坡体前缘高程 6065 m,宽 570 m,后缘高程 560584 m,宽 440 m,纵向长 1 250 m,平面上总体呈长条形,面积约为 0.7 km2,估算方量为2 300104 m3,滑坡平面如图 3 所示。1ZG286ZG285ZG284ZG283ZG282ZG281长0100200
14、mN剖面及编号图例滑坡边界GPS监测点及编号江1图 3 滑坡平面图3.2 滑坡物质组成及结构特征1)滑体:据钻孔揭露,滑坡体物质平均厚度3035 m,滑坡中部最厚达71 m,滑体组成物质以残积和崩坡积物为主,滑体土层以块石土、碎石土为主,块石、碎石母岩成分以泥岩、灰岩及砂岩为主,滑坡地层如图 4 所示。125175225275325375425475高程/m174路公路长江175 m145 mZG286ZG285QdelZG284T2b3225路高程/m间距/m50100 150 200750700650600550500450400350300250900850800950 1000 105
15、0 1100125175225275325375425475图例Qdel第四系滑坡堆积物碎石土T2b三叠系中统巴东组泥岩、砂质、灰岩滑面岩层产状GPS监测点编号3225ZG285图 4 滑坡剖面图2)滑床:高程在 440 m 以下地层为紫红色粉砂岩及泥岩,高程 440 m 以上地层为灰黄色泥灰岩。滑床总体形态为一斜坡地形,总体坡度为 2022,中间凹陷,侧缘突起,呈凹形。3)滑带:地表未见明显的滑带露头,上覆崩坡积物与下伏基岩接触带为滑带。4 库水作用下滑坡的稳定性分析4.1 计算模型滑坡计算模型采用有限元划分网格,共有126 000 个节点,116 127 个单元。模型单元类型采用减缩积分单
16、元 C3D8RP,该类型单元可进行渗流-应力耦合,可以将所有节点孔隙压力作为一个自由度,保证了模型计算精度及速度。模型计算网格如图 5 所示。图 5 模型计算网格模型中考虑滑坡内的 3 种典型土体:滑体、滑带和滑床,各种材料均假定为各向同性体,综合对滑坡试验数据的统计、类比与反算分析,其抗剪强度参数取值及滑坡各相关参数取值见表 1。由于滑坡中的岩土体物理性质及地质条件较为复杂,其各项强度指标具有高度的非线性分布特征,因此无法采用线性本构模型对其进行模拟,而应采用非线性模型对其进行有限元模拟分析。经多次对比分析论证后,滑坡采用 Mohr-coulumb 模型来研究分析。697江 西 科 学202
17、3 年第 41 卷表 1 滑坡典型土层相关参数典型土层渗透系数K/cms-1弹性模量E/Mpa泊松量 孔隙比 eC/KpaC/Kpa/r/KNm-3r/KNm-3滑体6.010-312.218.5 0.250.540.6614.826.517.518.416.218.912.516.520.122.522.523.4滑带2.210-62.713.40.200.450.6716.417.813.515.415.816.912.113.816.419.621.123.5滑床3.210-67.21040.180.11 1501 25060070029.535.126.130.323.225.326.
18、428.4 计算模型应力分析的边界条件为:底部为三向约束,四周均采用法向约束;渗流分析的边界条件为:底部为不透水边界,四周和坡面均为自由透水面。4.2 计算工况滑坡位于三峡库区内,对滑坡稳定性产生影响的外部因素主要是库水位变化,因此本章主要对以下 3 种工况进行数值模拟,以此分析滑坡在库水位变化影响下的稳定性。工况 1:自重+地表荷载+初始水位(135 m 静水位);工况 2:自重+地表荷载+145 m 水位;工况 3:自重+地表荷载+175 m 水位。4.3 模型计算结果及分析以滑坡及滑坡的纵断面 1-1断面为研究对象,分别从应力场、孔隙水压力、位移及塑性区等方面论述 3 种工况下滑坡在库水
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 滑坡预报与有限元实例分析对比研究 1 滑坡 预报 有限元 实例 分析 对比 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。