非线性抗剪强度与线性抗剪强度的堆石坝抗震稳定性分析.docxVIP

非线性抗剪强度与线性抗剪强度的堆石坝抗震稳定性分析 1 抗剪强度的计算方法 目前，土壤侵蚀稳定分析广泛应用于瑞典圆弧法和简化法中的离散法。根据试验,堆石材料的抗剪强度更符合非线性表达式,而以往多采用莫尔-库仑准则进行边坡稳定分析的方法。故拟静力法分析堆石坝下游坝坡的稳定性,分别采用瑞典法和毕肖普法,对每一种方法又分别计算了采用线性抗剪强度与采用非线性抗剪强度的坝坡稳定性。 1.1 竖向地震堆石地震我国[12 瑞典圆弧法 Fs=∑[(W±Q′)cosβ-Qsinβ]tgφ∑[Qcosβ-Qh/2R+(W±Q′)sinβ](1) 式中:W为土条重量;Q、Q′分别为水平、竖向地震惯性力;R为滑弧半径;h为土条高度;β为土条底滑弧面与水平面的夹角;φ为堆石内摩擦角,堆石的凝聚力为零。 简化的毕肖普法 Fs=∑(W±Q′)tgφsecβ1+tgβtgφ/Fs∑[Qcosβ-Qh/2R+(W±Q′)sinβ](2) 1.2 土条底部法向应力 邓肯非线性抗剪强度公式的内摩擦角 φ=φ0-Δφlg(σ3Ρa)(3) 式中:φ0为堆石在法向应力为100kPa时的内摩擦角;Δφ为试验常数;σ3为小主应力;Pa为大气压力。 σ3=σn(sec2φe-tgφesecφe)(4) 式中:σn为土条底部的法向应力;tgφe=tgφ/Fs。 用邓肯非线性抗剪强度,瑞典法和简化毕肖普法的表达式是一样的。 Fs=∑Νtgφ∑[Qcosβ-Qh/2R+(W±Q′)sinβ](5) 用瑞典法求解Fs,式中土条底部法向力N已知, Ν=(W±Q′)cosβ-Qsinβ 由于按(3)式求φ需σ3,按(4)式求σ3则需要φe,φe中又包含Fs,故需迭代计算。 用简化的毕肖普法按(5)式求解安全系数,其中的N也未知 Ν=W±Q′cosβ+tgφsinβ/Fs 故也需迭代计算。 1.3 明确n的法定地位 堆石材料的非线性强度准则除了邓肯公式外,还有德迈洛(de Mello)公式应用也较普遍。 τf=Aσmn(6) 式中:A、m分别为堆石的抗剪强度系数与指数。式(6)对应堆石的摩擦角为 φ=arctg[Aσ(m-1)n](7) 按德迈洛抗剪强度计算坝坡稳定安全系数的公式为 瑞典圆弧法 Fs=∑A[(W±Q′)cosβ-Qsinβ]mt1-mcosm-1β∑[Qcosβ-Qh/2R+(W±Q′)sinβ](8) 式中,t为土条宽度。 简化的毕肖普法 Fs=∑Aσmntsecβ∑[Qcosβ-Qh/2R+(W±Q′)sinβ](9) 由于简化的毕肖普法考虑了水平条间力,故σn无法像瑞典法那样直接表达出来 σn=Νtcosβ 式中:N为土条底部滑弧面上的法向力。列土条的竖向平衡方程式,并用牛顿法迭代求N Νk+1=Νk-f(Νk)f′(Νk) 式中f(Ν)=(W±Q′)-Νcosβ-AFsΝmt1-mcosmβtgβ。迭代Fs和N即可求出安全系数。 2 单元应力、动剪应力及一阶段i 运用有限元法计算出坝坡的静应力和每一瞬时的动应力,可利用其成果分析坝坡的稳定性。作用于滑动面的法向应力σN和切向应力τN为 σΝ=σy+σx2+σy-σx2cos2β-τxysin2β τΝ=σy-σx2sin2β+τzycos2β 式中:σx=(σxs+σxd),σxs为单元的静水平应力,σxd为单元的动水平应力;σy=(σys+σyd);τxy=(τxys+τxyd),σys、τxys分别为单元的静竖向应力、静剪应力,σyd、τxyd分别为单元的动竖向应力、动剪应力。 安全系数Fs为 Fs=∑σΝitgφili∑τΝili(10) 式中:φi、li分别为滑弧面上第i单元的内摩擦角和弧长。φi按线性的莫尔-库仑准则为常数,也可按邓肯非线性抗剪强度公式计算。 德迈洛公式计算抗剪强度的安全系数为 Fs=∑AσmΝili∑τΝili(11) 3 水平试算线v-v转换 先假设弧面滑出点B,如图1。求通过B点的最危险滑弧及其相应的安全系数。寻求的方法是:初选水平试算线H-H,作为第一根试算线,将H-H分成n等份即有n+1个分点,以每个分点为圆心,以分点至滑出点B的距离为半径,作画弧计算。比较n+1个分点的安全系数Fs值,即可找到H-H线上最小的抗滑安全系数Fs1及圆心坐标PH。然后过PH作一根垂直试算线V-V,用同样的方法确定Fs2及PV。为保证计算精度,V-V线段的长度是H-H线段的0.8倍,以下类推。过PV点再做水平试算线H′-H′,求Fs3及PH′点。以此循环,直到满足精度要求,即 (Fsi-Fsi-1)/Fsi≤ε(12) 式中:Fsi、Fsi-1分别为第i、i-1次试算线上各分点中最小的安全系数;一般可取ε=0.001。 最危险滑出点寻求采用0.618优选法,如图2。设函数f(x)在[A,B]内有极小