单位文秘网 2021-08-30 09:03:09 点击: 次
摘要:受地形条件的限制,高速公路建设中出现了大量的高填方路基。建设这类填方很高路堤时,路堤的稳定就成为一项直接影响公路建设质量的因素,本文采用三维显式有限差分法程序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——试样含水量(%)
K=pd /pc ×100 (3-2)
式中:K——测试点的施工压实度(%)
pd——试样干密度(g/cm3)
pc——由击实试验得到的试样最大干密度(g/cm3)
3.2.2结论
3.2.2.1该段路基表层:表面松散,未碾压密实;含水量偏高。
3.2.2.2抽捡路段共取土样8个,土样含水量普遍稍高,具体数据见表1。
压实度大于等于94%的点为3个,所占百分率:37.5%;压实度大于等于93%的点为1个,所占百分率:12.5%;压实度大于等于90%的点为4个,所占百分率:50.0%;压实度达不到94%的设计要求。从压实度方面看,该段高填方处于未压实的状态,验证了高填方路基的稳定性计算。
4结论
4.1 FLAC-3D基于显式有限差分法的拉格朗日算法,内置了包含1个“空”模型、3个弹性模型和6个塑性模型等l0种材料模型和梁、锚杆等4种结构单元,特别适用于岩上工程结构的弹塑性力学行为、模拟施工过程等方面。
4.2 FLAC-3D抗剪强度折减法在求解过程中采用动态松弛技术,不需要求解大型联立方程组,计算简单、方便,节省了大量的内存和计算时间,同极限平衡法相比不仅能满足力的平衡方程而目也满足土体的应力和应变的关系。
4.3本文在用FLAC-3D对实际高填方路基工程进行数值模拟,采用强度折减法计算高填方的稳定后,对影响高填方路基稳定性的压实度进行检测,验证了理论计算的准确性。
参考文献:
[1] Bishop,A,W.. The use of the slip circle in the stability analysis of slopes,Geotechnique,1955,5(1):7~17.
[2] 陈祖煌. 边坡稳定分析的塑性力学上限解和下限解.中国水利电力技术发展和成就. 北京:中国电力出版社.1997.198~208.
[3] Greco,V.R.. Efficient Monte Carlo technique for locating critical slip surface. Journal of Geotechnical Engineering, ASCE,1996,122(7):517~525.
[4] 黄岩,文江泉. 昆明地区土质高边坡稳定性分析.城市道桥与防洪, 2006(1):17~20.
[5] 朱宗明,FLAC在边坡整治工程中的应用.广州建筑,2003(4):40~42.
[6] 李新平,郭运华,彭元平等.基于 FLAC-3D的改进边坡极限状态确定方法. 岩石力学与工程学报, 2005.24(2):5287~5291.
[7] 胡焕校,刘静,佘重九,龚继文.高填方公路路堑边坡稳定模拟分析及治理措施. 中南大学学报(自然科学版),2004,35(5): 856~859.
[8]谢春庆.山区机场高填方块碎石夯实地基性状及变形研究(博士学位论文),成都理工大学,2001.
[9] 陈祥军,汤劲松. 用FLAC进行马崖高边坡稳定性分析. 石家庄铁道学院学报,2002.15(3).76~79.
[10] 尹紫红,周志林,姚令侃.用软件分析高速公路边坡的稳定性. 路基工程,2006(1).45~47.
[11] 魏继红,吴继敏,孙少悦.FLAC3D在边坡稳定性分析中的应用. 勘察科学技术, 2005(2):27~30.
[12] 刘波,韩彦辉(美国). FLAC原理、实例与应用指南. 北京:人民交通出版社, 2005.
[13] 沈珠江.理论土力学. 中国水利水电出版社, 北京, 2000.
[14] Itasca Consulting Group Inc, FLAC3D (Fast Lagrangian Analysis of Continua in 3 Dimensions) version 2.1 User"sManual,USA: 2002.
作者单位:1.湖南中大建设工程检测技术有限公司
2.中南大学资源与安全工程学院
3.中南大学地学与环境工程学院
(责任编辑:单位文秘网) )地址:https://www.kgf8887.com/show-143-90332-1.html
上一篇:裂隙水流温度变化情况下的岩体非稳定温度场有限元数值分析
下一篇:空心楼板性能分析
版权声明:
本站由单位文秘网原创策划制作,欢迎订阅或转载,但请注明出处。违者必究。单位文秘网独家运营 版权所有 未经许可不得转载使用