二期围堰稳定计算

小南海二期混凝土纵向围堰稳定及应力计算报告1计算背景小南海水电站二期导流拆除一期上、下游横向土石围堰进行左侧主河床截流，修建二期上、下游横向土石围堰形成二期基坑。二期导流期间，江水由导流明渠下泄，船舶从导流明渠通行。二期导流明渠左岸与坝身连接段为纵向混凝土围堰，基本体型见图1-1。三期导流拆除二期上、下游横向土石围堰，封堵导流明渠，修建三期上、下游横向土石围堰形成三期基坑。三期围堰与左岸大坝共同挡水后，第一批机组发电，船闸通航；江水由左岸24孔溢流坝下泄。混凝土纵向围堰在三期机组发电时参与挡水，其中坝轴线下游纵向围堰与永久建筑物结合参与挡水。二期混凝土纵向围堰在二期和三期导流时均要挡水，为了保证其安全性，有必要对其进行不同挡水条件下的围堰稳定计算及建基面应力分析，看其能否满足设计要求。图1-1纵向围堰剖面及剪切带分布2计算资料计算工况下的水位见表2-1，根据地质提供的资料建坝岩体向下游微倾或平缓卧于坝基下，层间剪切带发育，断层、裂隙等缓倾角结构面也有分布，坝基抗滑稳定问题突出。从单滑面滑移模式来看，坝基范围分布的剪切带主要有3S603、3S602、3S601、3S502、3S501、3S418、3S417、3S416、3S415、3S414、3S413、3S412等。其中多为剪切不充分的Ⅱ类剪切带，仅3S602、3S502、3S412属剪切充分的Ⅰ类剪切带。纵向围堰建基面下面的剪切带主要有3S603、3S602、3S601、3S502、3S501、3S418等。中坝址剪切带及结构面抗剪强度参数建议值见表2-2。地质提供中坝址岩体物理力学参数建议值表见表2-3。表2-1纵向围堰挡水时各计算水位序号项目水位（m）1P=5%（二期过流上纵）192.572P=2%（二期过流上纵）193.973三期过流发电水位197.004P=5%（三期过流下纵）193.195P=2%（三期过流下纵）194.60表2-2中坝址剪切带及结构面抗剪强度参数建议值剪切带或结构面类型抗剪断强度抗剪强度特征描述f＇c＇（MPa）fc（MPa）Ⅰ类泥夹岩屑或薄泥型0.28～0.300.06～0.080.22～0.250具泥化特征Ⅱ类岩屑夹泥型0.32～0.350.10.30岩屑、部公泥化岩块岩屑型0.35～0.40.120.35破碎岩块、岩屑层理面0.45～0.550.15～0.20.35～0.40软岩层面取低值缓倾角断层0.25～0.280.050.2～0.250粘土岩中断层具泥化0.30～0.320.08～0.10.28～0.300砂岩中断层破碎带表2－3中坝址岩体物理力学参数建议值表岩石名称天然重度（ρ）单轴湿抗压强度抗拉强度变模（E0）弹模（Ee）岩体本身抗剪断强度岩体抗剪强度混凝土/岩体抗剪断强度混凝土/岩体抗剪强度泊松比代表性岩组微、新岩体弱风化强风化微、新岩体kN/m3MPaMPa（GPa）（GPa）f＇C＇（MPa）f＇C＇（MPa）f＇C＇（MPa）fC（MPa）f＇C＇（MPa）fC（MPa）μ粘土岩、粉砂质粘土岩24.5～255～60.6～0.751.5～22～2.50.650.30.450.20.300.050.4000.600.250.400.32～0.36J3s2、J3s4－1、J3s6、J2s12泥质粉砂岩、粉砂岩25.5～2610～151～1.22.5～3.03.0～3.50.800.50.550.300.350.10.5000.650.30.4500.28～0.30J3s3、J3s4－2长石石英、砂岩、岩屑长石砂岩25.540～504～56～810～121.01.00.700.800.450.40.6500.80.80.600.22～0.24J3s1、J3s53计算原理3.1围堰建基面抗滑稳定计算围堰建基面抗滑稳定根据承载力极限状态确定，采用《混凝土重力坝设计规范》（DL5108-1999）及《水工建筑物荷载设计规范》（DL5077-1997）推荐的公式计算：a.基本组合采用下列极限状态设计表达式：),,(0KKQKGaQGS≤),(11KmKdafR式中：γ0——结构重要性系数，对于结构安全级别为Ⅰ、Ⅱ和Ⅲ级的结构或构件，可分别采用1.1、1.0、0.9；ψ——设计状况系数，对应于持久状况、短暂状况、偶然状况，分别取用1.0、0.95、0.85；作用效应函数RPS)(抗滑稳定抗力函数RRRRAcWfR式中：∑RP——坝基面上全部切向作用之和，kN；Rf——坝基面抗剪断摩擦系数；Rc——坝基面抗剪断粘聚力，kPa。G——永久作用分项系数，见表3-1；Q——可变作用分项系数，见表3-1；kG——永久作用标准值；kQ——可变作用标准值；ka——几何参数的标准值；kf——材料性能的标准值；γm——材料性能分项系数，见表3-2；γd1——基本组合结构系数，见表3-3；b.偶然组合采用下列极限状态设计表达式：),,,(0KKKQKGaAQGS≤),(12KmKdafR式中：AK——偶然作用代表值；γd2——偶然组合结构系数，见表3-3。表3-1作用分项系数序号作用类别分项系数1自重1.02水压力1）静水压力1.02）动水压力：时均压力、离心力、冲击力、脉动压力1.05、1.1、1.1、1.33扬压力1）渗透压力1.22）浮托力1.03）扬压力（有抽排）1.1（主排水孔之前）4）残余扬压力（有抽排）1.2（主排水孔之后）4淤沙压力1.25浪压力1.2表3-2材料性能分项系数序号材料性能分项系数备注1抗剪断强度1）混凝土/基岩摩擦系数f′R黏聚力c′R1.33.02）混凝土/混凝土摩擦系数f′c黏聚力c′c1.33.0常态混凝土层面3）基岩/基岩摩擦系数f′d黏聚力c′d1.43.24）软弱结构面摩擦系数f′d黏聚力c′d1.53.42混凝土强度抗压强度fc1.5表3-3结构系数序号项目组合类型结构系数备注1抗滑稳定极限状态设计式基本组合偶然组合（校核洪水）1.21.2包括建基面、层面、深层滑动面2混凝土抗压极限状态设计式基本组合偶然组合（校核洪水）1.81.83.2建基面（层面）应力围堰体型跟重力坝体型相似，故围堰的稳定和应力计算参照重力坝规范进行计算。围堰建面应力根据正常使用极限状态确定，采用《混凝土重力坝设计规范》（DL5108-1999）推荐的公式计算：正常使用极限状态作用效应的短期组合采用下列设计表达式：),,,(0KKKKSafQGS≤31/dC正常使用极限状态作用效应的长期组合采用下列设计表达式：),,,(10KKKKafQGS≤42/dC式中：C1、C2——结构的功能限值；Ss（·）、Sl（·）——作用效应的短期组合、长期组合时的效应函数；3d、4d——正常使用极限状态短期组合、长期组合时的结构系数；ρ——可变作用标准值的长期组合系数，本工程取ρ＝1。长期、短期组合下坝踵垂直应力不出现拉应力（计扬压力），计算公式为RRRRRJTMAW≥0式中：∑WR——围堰建基面法向作用之和，kN，向下为正；∑MR——全部作用对建基面形心的力矩之和，kN.m逆时针方向为正；AR——建基面的面积，m2；TR′——建基面形心轴到上游面的距离，m；JR——建基面对形心轴的惯性矩，m4。长期组合下上游坝面的垂直应力不出现拉应力（计扬压力），计算公式为JcMcTcAcWc'≥0式中：∑WC——计算截面上全部法向作用之和，kN，向下为正；∑WC——全部作用对计算截面形心的力矩之和，kN.m，逆时针方向为正；AC——计算截面的面积，m2；Jc——计算截面面积对形心轴的惯性矩，m4；Tc——计算截面形心轴到上游面的距离，m。短期组合下游坝面的垂直拉应力计算公式为：JcTcMcAcWc≤100（kPa）式中：TC——计算截面形心轴到下游面的距离。3.3大坝沿软弱结构面深层抗滑稳定计算坝基深层存在缓倾角软弱结构面时，具有单滑动面、双滑动面和多滑动面，可根据地质结构模型分析确定控制性滑动面，进行极限状态抗滑稳定分析。双斜滑动面为最常见情况，如下示意图8－1。图8－1双斜滑动面示意图采用刚体极限平衡等K法，根据《混凝土重力坝设计规范》（SL319-2005）,按抗剪断强度公式和抗剪强度公式，计算抗滑稳定安全系数。1）采用抗剪断强度公式计算。考虑ABD块的稳定，则有：coscoscossin)(csinsinsincos311'1311'1'1QUHGWAUUQHGWfK(E.0.2－1)考虑BCD块的稳定，则有：cossin)cos(Qsinsincos322'2322'2'2UGAcUUQGfK(E.0.2－2)式中：K´1、K´2—按抗剪断强度计算的抗滑稳定安全系数；基本组合：K≥3.0，特殊组合（1）：K≥2.5，特殊组合（2）：K≥2.3；W—作用于坝体上全部荷载（不包括扬压力，下同）的垂直分值,kN；H—作用于坝体上全部荷载的水平分值,kN；G1、G2—分别为岩体ABD、BCD重量的垂直作用力,kN；1'f、2'f—分别为AB、BC滑动面的抗剪断摩擦系数；1c'、2c'—分别为AB、BC滑动面的抗剪断凝聚力,kPa；A1、A2—分别为AB、BC面的面积,m2；α、β—分别为AB、BC面与水平面的夹角；U1、U2、U3—分别为AB、BC、BD面上的扬压力,kN；Q—分别为BD面上的作用力，KN；—BD面上的作用力Q与水平面的夹角。夹角值需经论证后选用,从偏于安全考虑可取00。2）采用抗剪强度公式计算。对于采取工程措施后应用抗剪断强度公式计算仍无法满足要求的坝段，可采用抗剪强度公式（E.0.3-1），公式（E.0.3-2）计算抗滑稳定安全系数。考虑ABD块的稳定，则有：coscoscossin)(sinsinsincos3131111QUHGWUUQHGWfK(E.0.3－1)考虑BCD块的稳定，则有：cossin)cos(Qsinsincos3232222UGUUQGfK(E.0.3－2)式中：K1、K2——按抗剪强度计算的抗滑稳定安全系数；基本组合：K≥1.35，特殊组合（1）：K≥1.2，特殊组合（2）：K≥1.1；1f、f2——分别为AB、BC滑动面的抗剪摩擦系数4计算结果及分析根据确定的宽顶堰过流量，按宽顶堰和台形堰两种计算原理进行上下游围堰拆除高程计算，以确定三期上下游围堰拆除的高程。上游围堰堰前水位取大坝坝前水位，上游围堰堰后水位取溢流坝下游水位；下游围堰堰前水位取大坝下游水位，下游围堰堰后水位取大中坝（中坝址）州尾水位进行初略计算。各平率洪水上下围堰过流流量见表4-1。表4-1各平率下围堰过流流量洪水平率总流量Q(m3/s)上游围堰过流量Q(m3/s)下游围堰过流量Q(m3/s)542600403940391047800480848082052600505050504.1上游围堰拆除计算按照表4-1的泄流能力对三期上游围堰进行拆除，表4-2、4-3为分别用宽顶堰和台形堰计算原理对上游围堰拆除高程的计算参数表，表4-4、4-5为计算结果表。表4-2上游围堰拆除高程水力计算参数（宽顶堰原理）洪水平率(年)堰前水位（m）堰后水位（m）流量系数m淹没系数收缩系数过水宽度B（m）5190.69190.140.340.971435.9010192.07191.460.340.971440.6320193.35192.680.340.981445.61表4-3上游围堰拆除高程水力计算参数（台形堰原理）洪水平率（年）堰前水位（m）堰后水位（m）自由出流系数pm淹没出流系数p过水宽度B（m）5190.69190.140.3510.912435.9010192.07191.460.3550.922440.6320193.35192.680.3510