地震对边坡失稳滑动面积的影响研究.pdf
《地震对边坡失稳滑动面积的影响研究.pdf》由会员分享,可在线阅读,更多相关《地震对边坡失稳滑动面积的影响研究.pdf(9页珍藏版)》请在咨信网上搜索。
1、 第4 4卷第4期V o l.4 4 N o.4 2 0 2 3青 岛 理 工 大 学 学 报J o u r n a l o f Q i n g d a o U n i v e r s i t y o f T e c h n o l o g y 地震对边坡失稳滑动面积的影响研究李春立1,李东贤1,李 亮1,*,徐 亮1,孟凯琪1,曹吉昌2(1.青岛理工大学 土木工程学院,青岛 2 6 6 5 2 5;2.住房和城乡建设部科技与产业化发展中心,北京 1 0 0 8 3 5)摘 要:提出采用S P H方法和极限平衡法计算的边坡失稳滑动面积来量化地震对边坡失稳后果的影响。首先在G e o-S t u
2、 d i o软件建立典型边坡模型,然后通过改变折减系数得到不同的最小安全系数,在此条件下对比分析了同一最小安全系数下地震和非地震边坡的最终和最初滑动面积,并进一步考虑了边坡坡高和临界位移取值对边坡失稳后果的影响。结果表明:地震对边坡的初始滑动面积提升幅度较小,对最终滑动面积随最小安全系数的减小而快速增大;随边坡坡高的增加,边坡最终滑动面积急剧增加,而对初始滑动面积增幅不大;临界位移值增大,对非地震边坡失稳后果影响不大,而地震边坡失稳后果随着最小安全系数的减小而愈加严重。因此在考虑地震对边坡失稳后果的影响时,建议综合考虑最小安全系数、坡高和临界位移值的影响。关键词:边坡失稳滑动面积;失稳后果;S
3、 P H方法;极限平衡拟静力法;地震边坡中图分类号:T U 4 3 4 文献标志码:A 文章编号:1 6 7 3-4 6 0 2(2 0 2 3)0 4-0 0 1 9-0 9收稿日期:2 0 2 2-0 7-0 6基金项目:国家自然科学基金资助项目(5 1 7 7 8 3 1 3)作者简介:李春立(1 9 9 9-),女,山东德州人。硕士,研究方向为边坡稳定性。E-m a i l:l i 1 7 2 8 5 0 6 0 4 2 1 6 3.c o m。*通信作者:李 亮(1 9 7 7-),男,山东济南人。博士,教授,主要从事岩土工程防灾减灾方面的研究。E-m a i l:l i l i a
4、 n g q u t.e d u.c n。S t u d y o f t h e i n f l u e n c e o f e a r t h q u a k e o n s l o p e f a i l u r e s l i d i n g a r e aL I C h u n l i1,L I D o n g x i a n1,L I L i a n g1,*,X U L i a n g1,ME N G K a i q i1,C A O J i c h a n g2(1.S c h o o l o f C i v i l E n g i n e e r i n g,Q i n g d
5、 a o U n i v e r s i t y o f T e c h n o l o g y,Q i n g d a o 2 6 6 5 2 5,C h i n a;2.C e n t e r o f S c i e n c e a n d T e c h n o l o g y&I n d u s t r i a l i z a t i o n D e v e l o p m e n t,M i n i s t r y o f H o u s i n g a n d U r b a n-R u r a l D e v e l o p m e n t o f t h e P e o p l
6、 es R e p u b l i c o f C h i n a,B e i j i n g 1 0 0 8 3 5,C h i n a)A b s t r a c t:T h e s l o p e f a i l u r e s l i d i n g a r e a i s c a l c u l a t e d b y u s i n g S P H m e t h o d a n d l i m i t e q u i l i b-r i u m m e t h o d t o q u a n t i f y t h e i n f l u e n c e o f e a r t
7、h q u a k e o n s l o p e f a i l u r e c o n s e q u e n c e.A t f i r s t,G e o-S t u d i o s o f t w a r e i s u s e d t o e s t a b l i s h t h e t y p i c a l s l o p e m o d e l,a n d t h e n t h e d i f f e r e n t m i n i-m u m s a f e t y f a c t o r s a r e o b t a i n e d b y c h a n g i
8、n g t h e r e d u c t i o n f a c t o r.U n d e r t h e s e c o n d i t i o n s,t h e f i n a l a n d i n i t i a l s l i d i n g a r e a s o f s e i s m i c a n d n o n-s e i s m i c s l o p e s u n d e r t h e s a m e m i n i m u m s a f e t y f a c t o r a r e c o m p a r e d a n d a n a l y z e
9、d.F u r t h e r m o r e,t h e i n f l u e n c e o f s l o p e h e i g h t a n d c r i t-i c a l d i s p l a c e m e n t v a l u e o n s l o p e f a i l u r e c o n s e q u e n c e a r e c o n s i d e r e d.T h e r e s u l t s s h o w t h a t t h e i n i t i a l s l i d i n g a r e a s o f t h e s l
10、o p e a r e i n c r e a s e d s l i g h t l y b y t h e e a r t h q u a k e,b u t t h e f i n a l s l i d i n g a r e a s i n c r e a s e r a p i d l y w i t h t h e d e c r e a s e o f t h e m i n i m u m s a f e t y f a c t o r s;w i t h t h e i n-c r e a s e o f s l o p e h e i g h t,t h e f i n
11、a l s l i d i n g a r e a s o f t h e s l o p e i n c r e a s e s h a r p l y,b u t t h e i n i t i a l s l i d i n g a r e a s d o n o t i n c r e a s e m u c h;w h e n t h e c r i t i c a l d i s p l a c e m e n t v a l u e i n c r e a s e s,i t h a s l i t-t l e e f f e c t o n t h e i n s t a b
12、i l i t y o f t h e n o n-s e i s m i c s l o p e,w h i l e t h e i n s t a b i l i t y o f t h e s e i s m i c s l o p e b e c o m e s m o r e s e r i o u s w i t h t h e d e c r e a s e o f t h e m i n i m u m s a f e t y f a c t o r s.T h e r e f o r e,i t i s s u g-青 岛 理 工 大 学 学 报第4 4卷g e s t e
13、d t h a t w h e n c o n s i d e r i n g t h e i n f l u e n c e o f e a r t h q u a k e o n s l o p e i n s t a b i l i t y,t h e i n f l u e n c e o f t h e m i n i m u m s a f e t y f a c t o r,s l o p e h e i g h t a n d c r i t i c a l d i s p l a c e m e n t v a l u e s h o u l d b e c o m p r
14、e-h e n s i v e l y c o n s i d e r e d.K e y w o r d s:s l o p e f a i l u r e s l i d i n g a r e a;f a i l u r e c o n s e q u e n c e;S P H m e t h o d;l i m i t e q u i l i b r i u m q u a s i-s t a t i c m e t h o d;s e i s m i c s l o p e图1 第i个土条受力分析Wi-土条重力;Qi-作用于土条质心的水平向地震力;Ni-土条底部法向力;Si-土 条
15、底部 剪切力;Ei,Ei+1-土条间法向应力;li-土条基底长度;i-滑动面与水平方向的夹角;R-圆弧滑动半径;O-圆弧滑动圆心 在众多地质灾害中,边坡失稳引发的后果尤为严重,造成了巨大的人身伤亡和财产损失。边坡稳定性的研究一直以来都是岩土工程领域中的热点问题之一,国内外学者在此方面开展了很多研究,譬如,S UN G K A R等1结合B i s h o p方法和有限元法,通过确定边坡在现有条件下和加固条件下的安全系数对边坡稳定性进行分析并提出了合理的加固措施;J I A N G等2采用强度折减有限元法对金平一级水电站左坝肩边坡进行了稳定性分析,计算了边坡开挖到不同高度时潜在破坏块体的安全系数
16、。作为诱发边坡失稳的主要诱因,地震作用对边坡稳定性的影响得到了广泛关注。地震边坡稳定性分析主要有以安全系数为评价指标的拟静力法3和时程分析法4,以永久位移为评价指标的N e wm a r k滑块分析法5和基于应力-应变分析的数值模拟方法6,如限元法、有限差分法、离散元法等。譬如,L I7等基于边坡的土体参数、几何形状和地震系数,利用有限元极限分析和莫尔-库仑失效准则,提出了一种估算土壤边坡静态、拟静力安全系数的新模型,以此来高效地评价边坡稳定性;李亮等3采用极限平衡拟静力法计算地震边坡的安全系数,研究了地震水平设计加速度值分别对深浅层滑坡临界滑动面的影响。除安全系数外,边坡失稳的严重性或者危险
17、性可通过边坡失稳滑动面积来衡量。C HU等8结合极限平衡法和安全系数等于1对应滑动体最大面积的假设对边坡的失效风险进行了量化并识别了边坡失效的来源;HUA N G等9基于极限分析法提出了一种新的滑坡风险定量评估框架,通过量化边坡失稳滑动面积对边坡的多种失效模式进行了稳定性评价;L I等1 0将随机有限元法与蒙特卡罗模拟(MC S)方法相结合,量化了边坡的失稳滑动面积并对土壤边坡进行了可靠性分析和风险评估。在已有的研究工作中,边坡失稳滑动面积通常由极限平衡方法计算所得,在有些情况下,非地震边坡(岩土体强度衰减引发边坡失稳)与地震边坡(地震引发边坡失稳)可能会具有相同的最小安全系数及对应的失稳滑动
18、面积,传统方法认为两种边坡具有相同的危险程度,考虑到极限平衡方法所得失稳滑动面积仅仅为滑坡的初始滑动面积,不能合理量化边坡的失稳后果,本文利用光滑粒子流体动力学方法1 1-1 2对相同最小安全系数下,地震边坡和非地震边坡的滑动全过程进行了模拟,对比分析了地震和非地震边坡最终滑动面积的不同之处。1 地震边坡稳定性计算方法在地震边坡稳定性分析中,常常采用拟静力法对地震力进行简化,即将地震力视为地震作用系数与滑坡体重力乘积的静态荷载。基于极限平衡法并考虑圆弧滑动面,边坡滑动体划分的竖直土条受力分析如图1所示。对边坡进行静态稳定性分析,在极限平衡B i s h o p法迭代公式的基础上,增加水平地震力
19、,可得极限平衡拟静力法的计算公式:Fs=1ni=1Wis i ni+Qic o si ni=1Wi-Qis i ni t a ni+c li-c liFss i nit a nim(1)02第4期 李春立,等:地震对边坡失稳滑动面积的影响研究m=c o si+s i nit a niFs式中:Fs为潜在滑动面的抗滑稳定安全系数;i为土条底部倾角,();Wi为土条重力,k N;i为土条底部的内摩擦角,();li为土条长度,m;c为土条黏聚力,k N;Qi为水平向地震力,k N。在S L O P E/W中,水平地震力惯性效应仅由水平地震力系数kh来表示,为无量纲常数,因此水平地震力又可表示为Qi=
20、ahWig=khWi(2)式中:ah为水平地震加速度,m/s2;g为重力常数,m/s2;kh为水平地震力系数。2 光滑粒子流体动力学方法(S P H)光滑粒子流体动力学(S m o o t h e d P a r t i c l e H y d r o d y n a m i c s,S P H)是一种纯拉格朗日无网格数值模拟方法,能够模拟边坡引发、传播以及终止的整个滑动过程。与有限元等需要进行网格划分的方法不同,S P H数值模拟是将计算域离散为有限数量的粒子,可以模拟岩土领域中的大变形问题,而不受网格畸变的影响。图2 S P H数值模拟理论如图2所示,离散域中的每个粒子都被赋予状态变量(包
21、括速度、位移、应力、应变等)以及材料属性(包括重度、黏聚力、内摩擦角、弹性模量、泊松比等)。在此计算域内,粒子i与影响区域内其他粒子j相互作用,其作用通过核函数W()来量化,核函数W()的计算域在半径为h的影响区域内。其中核函数可采用三次样条函数来表示:W r,h =ad23-q2+q32 0q1162-q 3 1q20 q2 (3)式中:ad=1 57 h2为二维空间中的归一化因子;q=rh为粒子i与j之间的归一化距离;r为粒子i与j之间的距离,m;h=hi+hj2,hi和hj分别为粒子i与j到影响域边界的距离,m。根据质量守恒方程和动量守恒方程,S P H方法中粒子的运动信息可以表示为Di
22、Dt=-ivixi=NSk=1mjvi-vj Wi jxi(4)DviDt=1i ixi+fi=NS j=1mj i2i +j2j +Ci j Wi j xi+fi(5)式中:NS为粒子i影响区域内的粒子总数;为粒子中土体密度,k g/m3;m为粒子质量,k g;v为粒子速度,m/s;t为时间,s;Wi j为粒子i与j间的光滑核函数;为单个粒子的总应力张量,其表达式由土体的12青 岛 理 工 大 学 学 报第4 4卷本构模型决定;x为粒子坐标;,为不同的坐标轴;f为外力引起的加速度分量;Ci j为提高S P H计算数值稳定性的项。S P H通过一系列的时间步长来模拟边坡滑动。在计算时间t=0时
23、,每个粒子的速度、应力和应变都设为0,之后通过外力加速度分量对每个粒子施加重力,随着计算时间的增加t=0+t,所有施加了重力的粒子在引力作用下均以相同的加速度向下移动。通过设置边界附近的粒子不能超出边界,从而导致粒子间产生应力,根据D-P(D r u c k e r-P r a g e r)模型的应力-应变关系,由粒子间的应力产生应变,随后更新粒子的速度和坐标并迭代计算,直至粒子达到稳定状态或满足收敛准则,计算结束。3 基于S P H量化边坡的最终滑动面积在S P H中,为了判断边坡是否发生失稳,需确定判断粒子滑动的临界位移值,一般取坡高的1%5%。当粒子的滑动位移大于临界位移值时则判断为边坡
24、失稳,将大于临界位移值的粒子面积累加即为边坡的最终滑动面积。图3给出了S P H方法量化边坡最终滑动面积的流程:图3 S P H量化边坡最终滑动面积的计算流程图4 边坡模型1)确定计算边坡的几何尺寸、土粒间隔大小和材料属性等,采用编写的F o r t r a n程序生成边坡粒子离散区域,然后使用S P H方法模拟边坡的整个滑动过程;2)在整个边坡的滑动计算过程中,记录并保存每个粒子的滑动位移Ygj(j取值为1n,n为边坡离散粒子总数)并降序排列;3)设定临界位移值,然后将最大粒子位移值Yg 1与临界位移值进行比较,若Ygj,则判定该粒子发生滑动并将其面积记录为s(s=r2),将所有滑动粒子的面
25、积累加得到边坡最终滑动面积S。4 算例分析4.1 边坡模型与土体参数研究表明3,当边坡黏聚力较大、内摩擦角较小时,边坡将发生深层滑坡,采用极限平衡拟静力法对边坡进行分析,随水平地震力系数从0开始增大,临界滑动面加深加宽;当边坡黏聚力较小、内摩擦角较大时,边坡将发生浅层滑坡,随水平地震力系数从0开始增大,临界滑动面的位置基本不变。T a y l o rs c h a r t中显示当坡角 6 0 时,边坡将发生浅层坡趾滑坡1 3。为便于计算并对比研究地震边坡与非地震边坡的失稳滑动面积,本文根据以上结论选取滑动模式为浅层坡趾滑坡的均质边坡模型。如图4所示:坡角6 0,宽度为3 5 m,坡高H为5 1
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 地震 失稳 滑动 面积 影响 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。