基于ABAQUS的边坡临界滑动面的确定及其应用.docVIP

基于ABAQUS的边坡临界滑动面的确定及其应用 　　摘要： 利用大型有限元软件ABAQUS对湖南省某高速公路边坡滑坡段进行了数值模拟分析，确定边坡的最危险潜在滑动面和边坡稳定性安全系数，对出现滑坡的原因进行了总结。同时进一步对治理后的边坡进行了模拟分析，得出治理后的边坡稳定性安全系数有较大提高，能保证更改设计后的边坡稳定。 　　Abstract: An analysis was made on the subway slope in Hunan province by using the finite element software ABAQUS， and using the way to get the position and average stability safety factor of the most dangerous potential slip curve. Numerical simulation and analysis was made to the slope after management and the stability safety factor was greatly improved. It is validated that it can ensure the slope stability. 　　关键词： 数值模拟；临界滑动面；边坡稳定性安全系数 　　Key words: numerical simulation；critical slip surface；stability safety factor of slope 　　 中图分类号：P5 文献标识码：A文章编号：1006-4311（2011）18-0279-02 　　0引言 　　临界滑动面的位置与相应的最小安全系数的确定是分析边坡稳定的关键问题，其求解主要依赖于可靠的边坡稳定分析方法和高效的优化搜索方法。现有的边坡稳定分析方法主要分为极限平衡法和有限单元法，与传统的极限平衡法相比，有限单元法全面满足静力许可、应变相容和应力―应变之间的本构关系，可以作为一种理论体系跟为严格的方法用于土体稳定分析。同时因为是采用数值分析方法，可以不受边坡几何形状的不规则和材料的不均匀性的限制，因此，是一种比较理想的分析边坡应力、变形、稳定性的手段。 　　本文利用大型有限元软件ABAQUS对某高速公路边坡滑坡段进行了数值模拟分析，确定了边坡的最危险潜在滑动面和边坡稳定性安全系数，对快速、准确评价复杂边坡稳定性具有很强工程实用价值。 　　1工程实例 　　此边坡属黄土地区路堑边坡，设计采用锚喷支护，锚杆按梅花形（1.5m×2.0m）布置，上部锚杆设计长15m，上部锚杆设计长10m。实际施工中分三次开挖，第一次开挖至标高112.8m，第二次开挖至标高109.5m，第三次开挖至坡底。第二次开挖后不久便出现两处滑坡地段，滑坡体长达150m，体积达万余立方，在坡顶出现贯通的张裂缝，且在第二级坡和第一级坡间的小平台上产生局部塌方，部分小平台已滑塌。坡顶裂缝离坡底线的水平距离11.8m。滑坡治理费达20万元，补救措施工程量大且实施困难，造成了人力和财力极大损失。为避免类似的边坡灾害发生，本文运用大型有限元软件ABAQUS建立边坡稳定性的预测模型，总结出滑坡的原因，同时对治理后的边坡也进行了分析和评价。 　　滑坡段的物理力学参数依据勘察设计院提供的数据取值。具体计算物理力学参数取值见表1： 　　以滑坡段作为分析对象，坡高13.6m，共分为3级。坡纵断面长度为21.3m，其中平台宽2-3m，坡形图见图1。 　　2有限元分析计算 　　2.1 模型建立现采用大型有限元软件ABAQUS对该边坡滑坡段建立模型进行有限元分析，取边坡坡面的纵断面方向作为X方向，边坡高度方向作为Y方向。模型X方向宜取4倍坡纵断面长度，Y方向取为2倍坡高作为分析的边界。采用平面应变有限元对边坡进行分析，土体材料特性按均质弹塑性考虑，采用Drucker-Prager屈服准则，土体单元采用8节点平面应变缩减单元（CPE8R），单元数为2768个。计算模型边界条件除上部为垂直荷载边界外，其余各侧面和底面为法向约束边界，计算时按自重应力场考虑。 　　 　　2.2 计算结果及讨论经数值模拟计算，边坡土体的塑性应变等值线如图2所示，从图中可以看出边坡土体中产生两条明显的塑性带。其中一条坡脚与坡顶两处塑性应变分布区域恰好贯通，因此可以认为，该边坡的最小安全系数为1.03，这与采用简化毕肖普（Bishop）极限平衡分析法得出的安全系数小0.04，误差在可接受范围内，表明此法是合理、可行的；同时从ABAQUS可视化截图上可以清楚地看出临界滑动面的位置及塑性值大小，坡顶塑性带坡顶塑