冻融循环作用下边坡地震动稳定性研究震灾防御技术.PDF

第 8 卷 第 2 期 震灾防御技术 Vol. 8, No. 2 2013 年 6 月 Technology for Earthquake Disaster Prevention Jun., 2013 王文丽，王兰民，郑龙，2013. 冻融循环作用下边坡地震动稳定性研究. 震灾防御技术，8 （2 ）：156—163. 冻融循环作用下边坡地震动稳定性研究 1 王文丽1, 2) 王兰民1, 2) 郑 龙1) 1）中国地震局兰州地震研究所，兰州 730000 2 ）兰州大学土木工程与力学学院，兰州 73000 摘要 季节冻土区和多年冻土区边坡在经历冻融循环后其力学性质产生较大变化，边坡稳定性随之改 变，在地震动荷载作用下存在着巨大的安全隐患。本文基于 FLAC 数值软件，将考虑冻土相变的显热 熔算法引入热学计算模块，并结合静态边界动力计算模块，建立了适用于冻融循环条件下的边坡地震 动稳定性分析平台。研究了冻融循环对边坡土体地震动稳定性的影响。结果表明，在发生冻融循环的 季节冻土区或多年冻土已经退化的区域，地震动后边坡的塑性破坏形式主要取决于发生冻融循环的活 动层的范围。 关键词：冻融循环 边坡 地震动稳定性 引言 随着我国在寒冷地区大规模工程建设的深入开展，由工程扰动诱发的冻土退化与冻融循 环所引起的工程稳定性，以及构筑物病害等工程问题日益突出（王兰民等，2007 ）。尤其在季 节性冻土区，新修边坡会经历剧烈的冻融循环，其力学性质会产生相应变化。在此前提下， 边坡地震动稳定性也会随之改变（孟庆生，2011 ）。冻融循环正是通过改变土体强度参数来影 响边坡稳定性的。因此需要准确掌握岩土强度在冻融前后的变化规律，来研究冻融作用对地 震动稳定性的影响。 根据前人的研究，对于不同物理性质的土体，冻融作用对土体强度的影响结果存在较大 的差异（Swan 等，1998；Aoyama 等，1985；Ogata 等，1985；齐吉琳等，2003 ）。不同地区 的土体，其物理性质存在较大差异，在不同试验条件下会产生差异性较大的研究结果。因此 在对具体的边坡稳定性计算分析时，需要根据实际研究区域的土体性质来获取相关强度参数， 并借鉴已有的关于冻融循环后土体强度变化的研究成果，来研究冻融后土体地震动稳定性的 变化规律。 目前在岩土工程问题的计算中，广为使用的是有限元方法（李晶岩等，2011 ）。在融土 地区，土体的力学性质不会随地温的变化而发生较大改变，因此进行岩土工程问题的有限 1 基金项目 中国地震局兰州地震研究所地震科技发展青年基金“2011Q05 ”资助 [收稿日期] 2012-12-28 [作者简介] 王文丽，女，生于 1981 年。助理研究员，兰州大学在读博士研究生。主要从事地震动工程安全稳定性方面的 研究。E-mail ：tywangwenli_2005@163.com 2 期 王文丽等：冻融循环作用下边坡地震动稳定性研究 157 元计算时不需要考虑温度场影响。在季节冻土区或多年冻土区，浅层土体随着季节的交替 变化会经历反复的冻融作用，其力学性质也会随之发生较大改变（Leroueil 等，1991； Eigenbrod 等，1996；Ono 等，1997）。因此，在计算寒冷地区的相关岩土工程问题时，需 要考虑温度场变化的影响，以此来确定土体内发生冻融循环的范围（即：活动层）。地震动 作用下的边坡稳定问题一直是岩土工程难题之一，边坡在动力作用下的稳定性研究仍处于 探索阶段（赵涛等，2012 ）。在前人的边坡地震动稳定性的研究中，大多并未考虑活动层范 围的影响（Chuvilin 等，1988），这将在很大程度上影响到数值计算的准确性，增加边坡工 程的潜在危险。 本文基于 FLAC 软件开发，考虑到经历冻融循环边坡的地震动稳定性分析平台，将冻土 相变的显热熔法引入热学计算模块，并结合动力学计