不同水位的作用下堆石坝坝坡有限元强度折减法分析お.docVIP

不同水位的作用下堆石坝坝坡有限元强度折减法分析お 　　摘要：[HT5”K]以某加固后的堆石坝工程为例，采用有限元强度折减法，研究了在不同水位组合和坝顶通车条件下的坝坡稳定性。计算结果表明：高水位对应侧的坝坡稳定性相对较高；水位越低，上下游水位差越大，坝坡的安全系数越大，反之则越小；坝顶通车对坝坡的稳定性不利，尤其是水位组合17367～14330m和15.200～16730m时；该堆石坝在加固后的坝坡稳定性系数基本满足相关规范要求，最小安全系数出现在水位组合15.200～16730m时的上游坝坡，此时坝坡属于稳定状态，但安全储备较小。 　　关键词：[HT5”K]有限元分析；强度折减；水位变化；坝坡；稳定性 　　中图分类号：[HT5”SS]TU435 　　文献标识码： 　　[JY]文章编号：[HT5”SS]1674994.4（2015HK] 　　1引言 　　随着我国西部大开发战略的推进，在基础建设上面实现了大力发展。但在大型基础设施的建设过程中，不可避免的出现了许多大型边坡，如基坑边坡、大坝边坡、码头边坡以及路堑边坡等。边坡是一个古老而又经典的问题，其稳定性分析方法也一直是岩土领域的一个热点问题。目前，稳定性分析方法主要有极限平衡法、有限元极限平衡法、有限元法以及有限元差分法等。随着计算理论及计算机应用的逐步发展，强度折减法在工程中不断得到重视，其不需事先假设滑裂面的形状，可自动搜索滑裂面，且不需事先假设条块之间的条间力，并考虑了土体材料的弹塑性，可以反映边坡失稳破坏和塑性变形区的发展过程。 　　2工程概况 　　西南地区某拦沙坝工程由右岸混凝土面板堆石坝段、河床埋石混凝土溢流重力坝段和左岸埋石混凝土非溢流重力坝段组成，坝轴线总长27560m。其中，右岸混凝土面板堆石坝段轴线长1.2690m，河床坝段和左岸坝段总长100m。右岸混凝土面板堆石坝坝顶宽8m，坝顶高程175.30m，上游坝坡为1[DK]∶14，下游坝坡从下至上分别为1[DK]∶1、1[DK]∶2.1、1[DK]∶15。大坝填筑至坝顶高程175.3m，在长江三峡水库降低水位运行后，库区水位内外高差为4m时，下游坝坡出现了渗漏现象，之后又数次发现严重渗漏现象，经相关专家讨论研究，采用了坝顶充填灌浆对坝体渗漏进行加固处理。同时又因城区交通拥堵问题，为缓解交通压力，相关决策部门决定将该堆石坝作为左右两岸的临时交通道路使用，大坝在不同运行水位下，可能造成堆石坝灌浆加固后形成的心墙防渗体破坏、坝坡失稳等问题。鉴于此，有必要对该堆石坝的上下游坝坡进行稳定性分析。 　　3有限元模型 　　31计算模型与边界条件 　　根据王俊杰等提出的简化计算方法，将坝体的剖面简化，简化如图1所示。在Abaqus软件中，有限元强度折减计算分析采用平面应变问题来处理，单元划分采用Abaqus中的受控四边形网格（Quad-dominated）自动网格划分功能。有限元网格划分后生成的结点总数为1386，单元总数为1305，其中三结点三角形平面应变元（CPE3）有24个，四结点四边形平面应变元（CPE4）有1.281个。模型底部约束x、y两个方向的位移，模型两侧约束x方向的位移。[FL）] 　　[XC062TIF][TS（][HT6H]图[ST6HZ]1[ST6BZ]堆石坝坝体典型断面[ST][TS）] 　　[FL（2K2] 　　32材料参数 　　[JP3]筑坝材料参数详见表1，本构模型采用Mohr-Coulomb屈服准则，并采用非关联流动法则（剪胀角=0）。[JP] 　　33失稳破坏依据 　　有限元强度折减法将强度折减、极限平衡原理与弹塑性有限元计算原理相结合，把土体材料的粘聚力c和内摩擦角φ，按式（1）、式（2）所示的形式进行折减，其中Fr为所设折减系数，随后把折减以后的粘聚力cm和内摩擦角φm替代原来的粘聚力c和内摩擦角φ，再次迭代计算，通过不断调整折减系数，直到土体达到临界破坏，此状态下的折减系数即为所求的安全系数。 　　cm=c/Fr[JY]（1） 　　φm=arctan（tanφ/Fr）[JY]（2） 　　式中：c和φ是抗剪强度参数；Fr是强度折减系数；cm和φm是土体实际发挥的抗剪强度参数。 　　目前，有限元强度折减法中关于边坡失稳的的判断依据有3种：以有限元计算的不收敛为标准；将边坡特征部位的位移突变点作为标准；以计算过程中边坡出现贯通且连续的塑性区为标准。本文选用第2种评判标准，以坝顶点位移突变为失稳判据，当位移发生突变时对应的强度折减系数即为边坡的安全系数。 　　4稳定性计算分析 　　根据库区的不同运行水位，对加固后的坝体在通车和不通车作用下的上下游坝坡进行了稳定性分析，其中，通车荷载为247kN时表示坝顶通车时的作用的荷载