季节冻土区边坡冻融稳定性剖析.docVIP

季节冻土区边坡冻融稳定性剖析 　　摘要：随着西部大开发和振兴东北战略的实施，季节冻土区的基础设施建设将得到很大的发展，然而，季节冻土区边坡的冻融灾害是制约季节冻土去工程建设的主要问题之一。基于有限单元法，建立了季节冻土区边坡冻融稳定性的数值模拟途径，以东北某一路堑边坡为背景，对初始含水量和冻结负温对融化过程中季节冻土区边坡稳定性的影响进行定量分析。研究表明，土的初始含水量越大、冻结温度越低，边坡冻融稳定性越差。分析成果可进一步充实季节冻土区边坡冻融稳定性的研究内容，对工程实践提供一定指导作用。 　　关键词：季节冻土区；边坡稳定性；数值模拟 　　中图分类号：TU475文献标志码：A文章编号： 　2015）002018903 　　我国是世界第三冻土大国，季节冻土面积约占国土面积的535%，多年冻土面积约占国土面积的215%[1]。由于经济与国防建设的需要，我国在季节冻土区的开展了一大批基础设施建设，主要包括东北地区的中俄输油管道工程、哈尔滨至大连的高速铁路工程、301国道等。目前，我国现有三分之二以上干线交通位于季节冻土区。随着西部大开发和振兴东北战略的实施，季节冻土区的基础设施建设将得到更大的发展，然而，季节冻土区边坡的冻融灾害是制约季节冻土去工程建设的主要问题之一。目前，我国的科研工作者对常温地区的边坡稳定性问题进行了深入的研究，取得了丰硕的成果，但是针对季节冻土区边坡的研究很少且资料分散，致使理论研究严重滞后于工程实践[25]。因此，季节冻土区边坡冻融稳定性研究迫在眉睫。 　　鉴于上述，本文根据季节冻土区边坡春融期易发生滑坡的现象，基于有限单元法建立了季节冻土区边坡冻融稳定性的数值模拟途径，以哈同公路某一路堑边坡为背景，对初始含水量和冻结负温对融化过程中季节冻土区边坡稳定性的影响进行定量分析，并对计算结果进行了讨论。本文的研究对有效解决寒区筑路技术的棘手难题、避免或减轻道路冻融病害无疑具有极其重要的现实意义。 　　1季节冻土区边坡的数值分析 　　1.1条件简化 　　（1） 土体为各向同性均质体。由于中国东北季节冻土区自然地形无大的起伏，因而自然边坡和人工边坡的规模均不算大，加之盆地内部各地第四系沉积变化较小，故此同一边坡中土体一般为各向同性的均质体，所以这种简化假定基本合理。 　　（2） 土体中无盐分影响。由于中国东北盐分含量较高的季节冻土均分布于油田及其附近地区，而这些地区基本无边坡，存在边坡地区的土体中盐分含量一般很低，所以这种简化假定基本合理。 　　（3）不考虑地下水动水压力与静水压力的影响，并且无边界水份补给与排泄。 　　（4） 按照平面应变问题处理。由于一般边坡的横剖面较纵剖面尺寸小得多，所以这种简化假定基本合理。 　　1.2本构关系的选取 　　边坡冻融稳定与失稳破坏的关键在于土的抗剪强度，而土的初始含水量、冻结温度直接强烈影响冻土的抗剪强度，当土的最大剪应力达到破坏极限时，边坡便失稳破坏。鉴于此，采用以MohrCoulomb屈服准则为破坏准则的弹-塑性模型刻画土的本构行为。 　　1.3边坡稳定性判据 　　边坡稳定性安全系数定义如下： 　　Fs=[SX（]τf[]τM[SX）][JY]（1） 　　式中：Fs为安全系数；τf土的实际抗剪强度；τM为土临界状态下的最大抗剪强度。 　　引入MohrCoulomb屈服条件，则安全系数可改写如下： 　　Fs=[SX（]c+σntanφ[]cc+σntanφc[SX）][JY]（2） 　　式中：c为土的实际凝聚力；φ为土的实际内摩擦角；cc为土临界状态下凝聚力；φc为土临界状态下内摩擦角；σn为实际法向应力。 　　以坡体中塑性区贯通及特征结点位移突变的综合考虑，作为边坡冻融稳定性的评价判据。 　　1.4计算模型与网格划分 　　基于中国东北季节冻土区某一公路路堑边坡算例，阐述有限单元数值模拟方法在边坡冻融稳定性分析中的具体应用。算例边坡为一匀质的粉质黏土边坡，计算剖面见图1，坡高H=20 m、坡度β=45°。计算域分为冻融区与非冻土区，由8结点等参四边形单元剖分计算域。采用位移约束边界条件，即计算域左、右边界允许竖向位移、限制水平位移，下部为固定边界。基于虚功原理形成单刚矩阵，通过直接刚度迭加法组装总刚矩阵，由威尔逊逐步积分法求解。利用虎克定律求解单元应力。 　　图1粉质黏土边坡剖面 　　1.5土体物理力学参数 　　冻土与非冻土的区别在于冻土的物理力学参数是随着温度及含水量变化的，具体参数见表1。 　　2影响因素分析 　　实践证明，在季节冻土区，土的初始含水量和冻结温度为影响边坡冻融稳定性的两个主要因素。以下基于上述中国东北季节冻土区某一公路路堑边坡算例，通过有限单元数值模拟方法，进行季节冻土区边坡冻