《边坡地震稳定性分析.docx》由会员分享,可在线阅读,更多相关《边坡地震稳定性分析.docx(11页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、边坡地震稳定性分析基于地震作用下黏性土坡失稳滑动特点,以土体应力状态及其变化分析边坡失稳经过。通过分析地震作用下边坡不同部位土体应力状态和剪应力变化,结合实际地震边坡失稳毁坏特征,提出黏性土坡地震三段式滑动失稳机制。在分析该滑动失稳机制与有限元强度折减法之间应力关联的基础上,将两者结合应用于实际黄土地震滑坡动力稳定性分析。根据此考虑得到的动力安全系数相比拟其他方法,与极限平衡法得到的结果更为接近。关键词:黏性土坡滑动失稳机制;边坡地震稳定性分析;有限元强度折减法.引言地震边坡失稳机理和边坡地震稳定性分析是构成边坡地震稳定性研究的两大主要内容。随着室内外土体测试技术与土工数值分析方法的进步,两者
2、都得到了长足的发展并构成了各自的理论体系。然而,目前对于边坡地震失稳机理与稳定性分析的研究少有穿插和融合,将二者综合考虑用于边坡稳定性研究的理论和方法为数更少。地震边坡失稳机理往往侧重于以地震惯性力和孔隙水压力为根本诱因进行定性分析,对坡体应力分布及变化特征关注较少,导致大多数的地震边坡失稳机理很难与目前边坡地震稳定性分析法中较为成熟的有限元方法相结合。以极限平衡法和有限元法为代表的边坡地震稳定性分析方法大都专注于给出边坡动力安全系数和动力响应特征,而忽视对地震边坡失稳机理的研究。因而,探寻适宜的地震边坡失稳机理,开掘其与边坡地震稳定性分析方法间的内在联络,提出考虑地震边坡失稳机理的边坡地震稳
3、定性分析方法,具有重要的理论意义和应用价值。.反映边坡土体应力变化特征的地震边坡失稳机制边坡土体应力是地震边坡失稳机理与边坡地震稳定性分析方法的重要联络,要将两者结合应用于分析边坡地震稳定性,则需分析边坡失稳经过中土体应力变化特征。地震作用下黏性土坡的受力特点及滑动形式由于黏聚力的存在,黏性土坡不会像无黏性土坡一样沿坡面滑动,在地震作用下其滑动面多表现为上陡下缓,毁坏和滑动通常呈现出滑裂面高于坡脚、滑裂面穿过坡脚和滑裂面低于坡脚三种类型与静力作用下不同,处在地震作用下的土体,其应力状态是由地震构成的动应力和自重构成的静应力的叠加,且动应力的大小和方向时刻都在变化。地震作用下的垂直动应力和水平动
4、应力定义为:槡槡式中:为土体密度;、分别为土体中质点的竖向及水平向振动速度;为纵波波速;为剪切波速;为体积变形模量;为剪切模量。取动应力最大值和最不利于稳定的方向,则边坡一点土体应力状态如图所示。一般以为地震以横向作用为主,故取。任意的黏性土坡都能简化为图所示的平面斜坡形式。借鉴王恭先提出的滑坡三段式滑动形式,以斜坡不同位置土体应力状态及变化来描绘斜坡失稳状态,将地震作用以上述动应力的形式引入,按几何形态特征将黏性土边坡划分为坡底、坡面、坡顶三段,分别分析地震作用下黏性土坡受力特点。如图所示,对于三种类型的黏性土滑坡在失稳经过中受滑裂面构成及滑动的影响,坡顶段土体的竖直向大主应力基本保持不变,
5、水平向小主应力减小。在坡顶段土体的挤压和坡底段土体的阻碍作用下,坡面段土体大致与坡面平行的大主应力和与其垂直的小主应力均增大。对于前两种类型的滑坡,滑裂面在坡面出露,坡底段土体只会遭到微弱的挤压作用。而对于第三种类型的滑坡,滑裂面在坡底段外表出露,坡底段土体由于遭到坡面段土体和坡顶段土体的推力,致使其分别平行和垂直于滑动面的大主应力和小主应力均增大。受地震横竖向拉张作用和影响,坡顶段土体大主应力和小主应力均被削弱,坡面段土体大主应力被削弱,小主应力转变为拉应力后得到加强,坡底段土体大主应力和小主应力也都遭到削弱。地震作用下黏性土坡失稳土体应力变化特征为半定量化分析地震作用下黏性土坡失稳滑动经过
6、中不同受力段土体的应力状态变化,以图所示的一般平面土坡为分析对象,分别选取位于坡底段、坡面段及坡顶段内且埋深均为,铅垂向应力相等的、三点处的土体图,分析计算不同受力段土体在地震作用下边坡失稳前后的应力状态,揭示其剪应力变化规律。图为边坡稳定时和地震作用下开场失稳时不同受力段的应力状态。受滑裂面构成及滑动的影响,坡顶段点土体小主应力逐步减小;受坡顶段的挤压和坡底段的阻碍作用,坡面段点土体大主应力增大,小主应力则由压应力逐步减小并转变为拉应力;在坡面段蠕动土体的挤压下,坡底段点土体大主应力和小主应力都逐步增大。受地震横竖向拉张作用和影响,坡顶段点土体大主应力和小主应力均被削弱,坡面段点土体大主应力
7、被削弱,小主应力转变为拉应力后得到加强,坡底段点土体大主应力和小主应力也都遭到削弱。考虑到边坡土体应力状态的复杂性,为方便定量计算,假定土体均质且知足摩尔库伦强度准则,将位于坡面附近土体水平向应力定义为与土体埋深有关的函数:式中:为侧应力系数;为土体重度;为坡角。边坡坡顶段点剪应力计算结果为:图地震作用下边坡失稳前后各段土体应力变化为坡底段土体大主应力方向与竖直方向的夹角。根据以上剪应力计算结果,分别作出三种状态时不同受力段土体的摩尔圆图。从图中能够看出,与无地震作用下边坡失稳时相比,地震作用下不同受力部位土体的剪应力遭到的影响不同:坡顶段得到较大幅度增长,坡面段呈现小幅增长,坡底段则表现为小
8、幅降低。总体来看,地震作用下边坡失稳经过中剪应力大小及变化趋势表现为:坡顶段初始剪应力较小但增长迅速;坡面段初始剪应力较大且增长较快;坡底段初始剪应力较小且增长较慢。在地震荷载长持时作用下,边坡土体更容易发生拉毁坏,。上述分析中,坡顶段土体受地震作用影响发生因小主应力由压应力转变为拉应力的拉张型剪切毁坏,进而较好地解释了地震边坡上部往往呈现出拉毁坏特征且存在拉张型裂缝的现象。黏性土坡三段式滑动地震失稳机制边坡失稳滑动是由于其内部分土体发生剪切毁坏所致,而这种剪切毁坏的发生与否则取决于土体的剪应力大小和变化趋势。位于边坡体内不同部位的土体,其受力特点和应力水平本就不同。在边坡失稳前其剪应力水平不
9、同,在边坡失稳经过中,随着应力状态变化的差异,不同位置土体剪应力的变化趋势更是不尽一样。但就土体毁坏发展程度而言,剪应力在失稳前本来就大,且在失稳经过中增长快的土体接近毁坏的速度更快。图所示的地震作用下不同受力段土体剪应力初始水安然平静变化趋势表明,坡顶段可能同时甚至先于坡面段发生毁坏。一些震后顶部存在大量深长拉张型裂缝但整体却未发生失稳的黏性土坡的存在,也能讲明地震作用下此类边坡的失稳可能始于坡顶段的毁坏。综合考虑地震作用下黏性土坡不同受力段应力变化特征和失稳毁坏特征,将地震作用以应力的形式表征,以斜坡不同位置土体应力状态及变化来描绘斜坡地震失稳的特征,研究提出黏性土坡三段式滑动地震失稳机制
10、:地震作用下,坡顶段受地震拉张作用的影响,侧向支撑力迅速减小进而发生主动土压毁坏,构成后缘张裂缝;与此同时,坡面段因无法承受坡顶段下滑土体逐步增大的推力、本身上覆岩土产生的下滑力和地震惯性力而沿本来就已存在的强度较低的相对软弱带面发生蠕动;在坡面段和坡顶段剩余下滑力的推挤下,若坡底段土体不发生毁坏,则滑动面在坡面出露,滑体在坡面剪出,若坡底段土体产生被动土压毁坏并构成毁坏面,则滑动面空间得以完全贯穿三段土体,滑坡开场加速整体滑移。.考虑黏性土坡三段式滑动地震失稳机制的有限元强度折减法有限元强度折减法是目前较为成熟的边坡地震稳定性分析方法,若要将其与黏性土坡地震三段式滑动失稳机制相结合,则应找寻
11、该方法中与剪应力有关的定义并加以分析和利用。在有限元强度折减法中,边坡失稳是通过土体强度的折减变低导致土体毁坏来实现的。折减系数增幅的大小与土体接近毁坏速度的快慢有关。既然实际边坡失稳时土体应力变化下的剪应力初始大小和变化趋势决定了土体接近毁坏的速度,而在有限元强度折减法中土体接近毁坏的速度又取决于折减系数的增长速度,因而能够通过对不同受力段折减系数进行调整来实现黏性土坡三段式滑动地震失稳机制和有限元强度折减法的结合。在运用有限元强度折减法进行边坡地震稳定性分析时,由于不同受力段都在边坡受力到达极限平衡时停止折减,即折减时间一样,所以要知足不同受力段土体接近毁坏时速度不同的特征,三段土体的折减
12、系数初值应相等,终值应保证在使边坡受力到达极限平衡状态的前提下,知足坡顶段和坡面段较大、坡底段较小的条件。.数值模拟验证运用上述考虑黏性土坡三段式滑动地震失稳机制的有限元强度折减法进行边坡地震稳定性分析,选用定义更为明确合理的斜坡土体位移突变点法,并结合破裂面贯穿及计算收敛两种失稳判据来确定边坡动力安全系数。选用具有集合定义和场变量功能的有限元软件,利用其集合定义功能完成对边坡不同受力段的划分,将坡顶段、坡面段和坡底段折减系数、和定义为同一场变量。考虑到坡顶段与坡面段毁坏先后顺序难以确定,在保证边坡能够到达极限平衡状态的条件下,分为、及三种情况进行试算。选取年岷县漳县地震诱发的永光村西侧滑坡,
13、以验证拓展的有限元强度折减法的可靠性。该滑坡为地震诱发的典型黄土滑坡,坡高,平均坡度为,斜坡简化模型如图所示。根据已有的现场测试和室内土体力学实验结果,滑坡基岩上覆土体分为三层,各层土体物理力学参数如表所列。土体计算采用理想弹塑性本构模型及莫尔库仑屈从准则。边坡的有限元分析模型和网格划分如图所示。模型土体均采用四边形四节点平面应变单元,网格划分为个单元、个节点。模型左右边界采用黏弹性边界条件,底部为人工边界。在动力计算经过中,水平向地震荷载从模型底部由左到右输入。输入地震波采用岷县漳县级地震中距震中的岷县台记录的主震加速度时程图。其地震峰值加速度为,卓越频谱为,输入地震动作用时间为。有限元强度
14、折减法折减系数和水平向位移计算结果如图所示。折减系数的取值范围为,当折减系数取时,计算不收敛,边坡到达极限平衡状态,即该边坡以折减系数定义的动力安全系数为。该拓展方法首先应保证极限平衡状态时各段折减系数的取值范围能够到达原方法所得的折减系数;其次,由于边坡滑动时土体中同时存在残余强度和峰值强度,且前者约为后者的,取作为极限平衡状态时折减系数变化幅度的上限,即折减系数最大取值不超过。选取符合上述条件,且知足、及三种不同条件的三组折减系数表进行试算。图为拓展方法折减系数和水平向位移的计算结果,分别画出两种方法得到的折减系数位移曲线图。采用以边坡土体位移突变点法为主的综合判据法确定边坡动力安全系数,
15、并对两种方法得到的动力安全系数加以分析比拟。通过比照不同组计算结果如表所列能够看出:当时拓展方法所得动力安全系数较大;当时两种方法所得动力安全系数接近;当时拓展方法所得动力安全系数较小。该滑坡采用极限平衡法中简化法计算得到的动力安全系数为,采用瑞典条分法及法得到的动力安全系数分别为和,拓展后折减系数知足条件的有限元强度折减法得到的结果为。该结果显然与极限平衡法得出的结果更为接近,从数值模拟层面验证了黏性土坡三段式滑动地震失稳机制的合理性。这讲明考虑黏性土坡三段式滑动地震失稳机制的动力有限元强度折减法应用于黄土斜坡地震稳定性分析是合理可行的,以其为根据确定的动力安全系数是安全可靠的。.结论边坡不
16、同受力段土体剪应力遭到地震作用的影响不同:坡顶段得到较大幅度增长,坡面段呈现出小幅增长,坡底段则表现为小幅降低。总体来看,地震作用下边坡失稳经过中剪应力大小及变化趋势表现为:坡顶段初始剪应力较小但增长迅速,坡面段初始剪应力较大且增长较快,坡底段初始剪应力较小且增长较慢。结合黏性土坡滑动形式和地震失稳毁坏特征,提出了黏性土坡三段式滑动地震失稳机制。地震作用下,坡顶段受地震拉张作用的影响,侧向支撑力迅速减小进而发生主动土压毁坏,构成后缘张裂缝。与此同时,坡面段因无法承受坡顶段下滑土体逐步增大的推力、本身上覆岩土产生的下滑力和地震惯性力而沿本来就已存在的强度较低的相对软弱带面发生蠕动。在坡面段和坡顶段剩余下滑力的推挤下,若坡底段土体不发生毁坏,则滑动面在坡面出露,滑体在坡面剪出,若坡底段土体产生被动土压毁坏并构成毁坏面,则滑动面空间得以完全贯穿三段土体,滑坡开场加速整体滑移。边坡土体剪应力水平及其变化趋势是地震边坡失稳机理与边坡地震稳定性分析方法的重要联络,通过对有限元强度折减法中折减系数进行调整,能够充分体现黏性土坡三段式滑动地震失稳机制。考虑黏性土坡三段式滑动地震失稳机制的有限元强度折减法所得的动力安全系数与极限平衡法所得结果特别接近。此方法的提出,实现了地震边坡失稳机理与边坡地震稳定性分析方法的初步结合,为实际工程的黏性土坡地震稳定性分析提供了较为可靠的数值模拟根据。