高填方路基稳定性FLAC-3D分析.docVIP

高填方路基稳定性FLAC-3D分析 　　摘要：受地形条件的限制，高速公路建设中出现了大量的高填方路基。建设这类填方很高路堤时，路堤的稳定就成为一项直接影响公路建设质量的因素，本文采用三维显式有限差分法程序FLAC-3D模拟了长沙某段高填方路基的三维力学行为。通过现场取样试验得到该段高填方路基土体的材料属性，采用强度折减法对该段路基进行稳定性分析，并计算影响高填方路基稳定性的路基压实度等指标，作为路基稳定性计算的验证。 　　关键词：高填方路基； 稳定性； 压实度； 强度折减 　　中图分类号：U416.1 　　文献标识码：B 　　文章编号：1008-0422（2007）06-0072-04 　　收稿日期：2007-04-09 　　作者简介：秦敏（1976-），男（汉族），湖南岳阳人，工程师。 　　 　　1引言 　　 　　高填方路基稳定性分析采用和边坡稳定性分析相似的方法，主要有极限平衡法[1，2，3]、现代岩土数值方法[4~11]和概率法等。极限平衡法已经发展得相对成熟， 目前大多数规范均规定采用极限平衡法进行边坡稳定分析，但是极限平衡分析法不能解决分析高填方路基应力和应变的问题。后来出现有限元分析高填方路基的稳定性，有限元在解决小变形方面有其优越性，但通常的高填方路基破坏为大变形问题，有限元在解决大变形方面不十分方便，而基于拉格朗日差分原理的FLAC Fast Lagrangian Analysis of Continua）在计算过程允许材料发生屈服及流变，适合于解决岩土工程中经常遇到的大应变，因此本文采用FLAC-3D结合长沙某段高填方路基来分析高填方路基的大变形问题。 　　 　　2FLAC-3D基本原理 　　 　　2.1 FLAC-3D求解时采用了下面的计算方法 　　2.1.1模型的混合离散化:模型经过网格划分，物理网格映射成数学网格，连续介质被离散为若干互相连接的六面体单元，作用力均被集中在结点上，而在计算过程中，程序内部又将每个六面体分为以六面体角点作为角点的两个四面体，变量均在四面体上进行计算，六面体单元的应力――应变取值为其内四面体的体积加权平均值； 　　2.1.2显式有限差分方法:变量关于空间和时间的一阶导数均用有限差分来近似，运动方程和动力方程均采用显式方法求解； 　　2.1.3动态松弛方法:应用质点运动方程求解，通过增加非粘性阻尼力使系统运动衰减至平衡状态。 　　 　　2.2FLAC的计算方法 　　2.2.1连续介质几何方程 　　 　　 　　可以联立(2－1)、(2－2)、(2－3)三式即可求得各单元的应力速度张量、应变速度张量和速度矢量，求解过程见图1。 　　 　　3工程实例 　　 　　3.1 理论分析 　　3.1.1 强度折减法 　　莫尔库伦塑性模型适用于松散状和粘结状的散粒材料，对一般的岩土体力学行为都有很好的模拟[12]。根据文献[4,5]，强度折减法计算所得的安全系数同简化 Bishop法相差不大，因此本文采用莫尔库伦模型作为高填方路基的模拟，通过强度折减法求解安全系数。 　　 　　 　　3.1.2 某高填方路基稳定性分析 　　本文结合长沙某段高填方路基，该路基采用粉质粘土作为填充材料，最高填方达到10m，上部6m坡率为1:2，下部4m坡率为1:1.5，下底宽度为62m，对该段路基进行适当的简化，得到如图2所示的FLAC-3D模型。 　　图3~图6为内聚力和内摩擦角折减后得到的模型的计算参数。 　　图3为利用FLAC3D计算出的最大不平衡力随时间步长关系曲线，可以此次计算是收敛的，通过近10000步的计算，最终的最大不平衡力为20.49N，节点近似达到平衡。 　　图4是对高填方路基边坡坡率变换的点(9 0 5)点X方向位移进行监视得到的该点X方向位移随时步的变化的关系，其最大的位移为-1.37cm(负号表示向左边偏移)，该处位移较大。 　　图5中显示高填方路基的最大位移绝对值均超过1.57cm，图6中显示XZ方向的剪力最大值为6.7kpa，图6左侧窗口中显示了采用强度折减法计算得到的高填方路基稳定性安全系数fos为1.09，可以判定该高填方路基处于不稳定状态。 　　 　　3.2 压实度检测 　　3.2.1弯沉试验 　　首先要求施工单位对路基进行整平，整平后再对其全线进行弯沉试验。弯沉试验检测频率为每车道20m检测一次。 　　3.2.2钻芯试验 　　用φ110标准薄壁取土器取样，再用环刀法或蜡封法测试湿密度；用酒精燃烧法测试土样含水量。 　　压实度计算： 　　pd=p/(1+0.01w)(3－1) 　　式中：pd――试样干密度(g/cm3) 　　Ρ――试样湿密度(g/cm3) 　　w――试样含水量(%)