基于三滑裂面等可靠余度重力坝深层抗滑稳定计算模型.docxVIP

基于三滑裂面等可靠余度重力坝深层抗滑稳定计算模型 1 重力坝深层抗滑稳定计算方法研究 近年来，重要水库深层地震稳定分析引起了人们的关注。另一方面，由于水库底板岩石结构的薄结构，在施工和建设项目中普遍存在着客观要求。另一方面，重力水库设计中采用的双侧面深层地震稳定计算方法存在不确定因素，需要对其进行进一步的改进研究。现行规范主要问题有:(1)坝基岩体人为采用铅直分界面进行分割的合理性及其对稳定安全的影响;(2)分界面上的作用力是否满足平衡条件;(3)分界面上作用力的方向人为简化,降低了评价结论的客观合理性和科学可靠性。此外,DL5108-1999给出的分项系数双斜面稳定计算公式矢量力平衡物理意义不明确,分界滑裂面及其上下游滑体的可靠盈余程度不同。因此,探讨贴近实际工程客观条件、更科学合理和运用简单方便的重力坝深层抗滑稳定分析方法具有重大理论研究意义和工程运用经济价值。 本文根据已有模型试验和有限元分析结果的破裂面形态,把分界滑裂面倾角及其抗力方向作为未知变量,按照实际工程具体力学边界条件,建立三滑裂面同时满足分项系数极限状态平衡方程的等可靠余度重力坝深层抗滑稳定计算模型,推导出相应的计算公式,提出计算机数值求解方法步骤,并利用编制的通用计算程序与双斜面法计算结果进行对比分析,阐明通常条件下三滑裂面法计算结果是客观确定的极小值,而双斜面法具有简化性和不确定性。三滑裂面等可靠余度重力坝深层抗滑稳定计算模型最大限度地避免了人为假定因素影响,最小可靠余度计算结果揭示了相应客观边界条件确定的深层抗滑稳定内在本质。 有许多研究者对三滑裂面刚体极限平衡深层抗滑稳定进行了分析,王志良、陈国立结合模型试验成果破裂面的形态,提出三斜面深层抗滑稳定计算方法,得出简化边界条件的安全系数极值滑裂面倾角解析方法;陈际唐等对三斜面等安全系数法与双斜面法、有限元强度折减法深层抗滑稳定计算结果进行了对比分析;孙冠华等采用铅直分界面的三滑面等安全系数法,对抗力角的取值进行了分析。但后两位学者建立的分界面平衡方程式中扬压力的力学考虑均有一些错误。 分项系数法的深层抗滑稳定研究结果主要集中于有限元分析法中运用研究和双斜面等安全系数法的分项系数表达式研究,基于分项系数法三滑裂面等可靠余度重力坝深层抗滑稳定计算模型和方法的系统研究结果的相关报道尚未查阅到。 2 三段式抗滑稳定计算模型，包括三段式抗滑稳定层 2.1 承载力极限状态 分项系数法是根据可靠度原理,采用概率统计方式把荷载性质、材料性能的差异性和不确定性予以可靠度量化。因此,用分项系数法进行重力坝深层抗滑稳定性评价,计算结果暗含了一定的保证概率或可靠程度,比单一安全系数法更具有合理性。在建立重力坝深层抗滑稳定计算模型时,使用分项系数承载力极限状态设计表达式,引入余度因子λ把表达式变成等式方程形式: 则,当余度因子λ≥1.0时,表示满足可靠度要求的承载力极限状态式盈余程度;当余度因子λ1.0时,表示不满足可靠度要求的承载力极限状态式亏欠程度。 为了区别于安全系数的概念,把λ称为满足可靠度要求的余度因子,简称可靠余度。 2.2 裂面与水平线间的夹角 图中AB为坝基软弱结构面,BC为下游岩体滑裂面,BD为倾斜分界滑裂面,分别编号为1、2和3,滑裂面上的作用力及力学参数指标用角标形式表示。α、β分别为AB、BC滑裂面与水平线间的夹角,θ为BD滑裂面与铅直线的夹角,φ为BD滑裂面上的抗力R3与法线的夹角,h为AB滑裂面在上游端的深度,T为坝基处滑块的水平长度,Ni、Ti分别为滑动面上的法向反力、切向反力;Ri、Ui分别为作用在各滑动面上的反力和扬压力设计值;分别为坝基面以上所有荷载的设计值在铅直和水平方向的代数和;PS、G1分别为坝基滑动块岩体上游水平压力和自重设计值;G2、GW分别为下游滑动块可利用岩层面以下的岩体自重和上部的盖重水重设计值;H为下游滑动块滑动面以上所有水平作用力的设计值代数和。 2.3 抗力函数的选取 根据坝基块脱离体ABDE在垂直于AB面的平衡条件可得N1,代入莫尔库伦强度准则方程得T1=f1N1+c1A1。由平行于AB滑动方向的极限平衡条件,可知AB滑裂面以上所有荷载在滑动方向的分力代数和即为作用效应函数S1(·),而结构面上的切向反力T1即为抗力函数R1(·),将作用效应函数和抗力函数代入式(1)得: 同理,依据下游块脱离体BCD的平衡条件,可得沿BC滑裂面方向的作用效应函数S2(·)和抗力函数R2(·),并代入式(1)得: 由滑裂面BD沿切向平衡条件,可得沿面的作用效应函数S3(·)和抗力函数R3(·),代入式(1)得: 式中:Ai为各滑动面面积;fi、ci分别为各滑动面的抗剪断摩擦系数和黏聚力设计值。 以上3个分项系数承载力极限状态方程中共有7个未知量,当3个滑裂面具有相等可靠余度时,即λ1=