面板堆石坝堆料的本构模型研究.docxVIP

面板堆石坝堆料的本构模型研究 1 基于堆石材料特性的实际堆石坝地震资料 在设计和建造环保石坝时，堆石体的变形控制是影响水库和板的工作性能的重要因素。就目前而言,有限元计算分析是进行面板堆石坝应力变形分析的主要手段,而材料的本构关系又是有限元计算分析的核心问题之一。在目前国内的面板堆石坝数值计算分析中,最常用的堆石本构模型是Duncan-Chang模型、双屈服面弹塑性模型、非线性解耦K-G模型等。但这些模型的参数都需要试验测得,增加了工作难度、费用。 由于堆石材料的变形特性与其岩性、级配、空隙率和密度有着极为密切的相关关系,在岩性和级配确定的情况下,能否不通过三轴试验而根据密度或空隙率直接得到堆石料的本构模型的参数呢?本文就是基于这个思想,通过对国内典型面板堆石坝邓肯张E-B模型的参数进行统计分析,得出了邓肯张模型各参数与空隙率、密度的关系,从而不需试验即可充分了解模型参数的变化趋势与数值范围。由于邓肯模型的广泛应用,已经积累了相当多的试验资料,为本文的研究提供了基础。 2 坝体堆弹性模量模型 1963年,康纳(Kondner)根据大量土的三轴试验的应力应变关系曲线,提出可以用双曲线拟合一般土的三轴试验曲线,即: σ1-σ3=εaa+bεa(1)σ1?σ3=εaa+bεa(1) 式中:a、b为试验常数。 邓肯等人根据这一关系提出了一种目前被广泛应用的增量弹性模型,即邓肯张模型( E-B 模型) ,其切线弹模和体积模量分别按下式计算: E=ΚΡa(σ3Ρa)n[1-Rf(1-sinφ)(σ1-σ3)2Ccosφ+2σ3sinφ](2)Bt=ΚbΡa(σ3/Ρa)m(3)E=KPa(σ3Pa)n[1?Rf(1?sinφ)(σ1?σ3)2Ccosφ+2σ3sinφ](2)Bt=KbPa(σ3/Pa)m(3) 同时模型还考虑内摩擦角φ随围压σ3的变化: φ=φ0-Δφlg(σ3/Ρa)(4) 对于卸载情况,该模型采用回弹模量进行计算: Eur=ΚurΡa(σ3Ρa)nur(5) 式中:C、φ为坝体堆石料强度指标;K、n为弹性模量数和弹性模量指数;Rf为破坏比;Kb、m为体积模量数和体积模量指数;Pa为大气压力;σ1和σ3为最大和最小主应力。 由于参数较多,首先进行参数敏感性分析,得出对数值计算结果影响较大的参数,然后对其进行统计分析,进而得出参数与空隙率、密度的关系式。 3 参数集的确定 参数敏感性分析的思路如下:设有一结构系统,其系统特性P主要由n个因素α={α1,α2,…,αn}所决定,系统模型为P=f(α1,α2,…,αn)。在某一基准状态α*={α*1,α*2,…,α*n}下,系统特性为P*。各因素在各自的可能范围内变动时,引起系统特性P偏离基准状态P*的趋势为ΔP,其相对于各自变因素的变化便被表征为它的敏感程度,即敏感性。在基准参数状态所决定的基准参数集α*={α*1,α*2,…,α*n}下,分析参数αk对特性P的影响时,可固定其他参数为基准参数,而令αk在其可能的范围内变动。则将系统P的相对偏差δP=|ΔP|/P和参数αk的相对偏差δαk=|Δαk|/αk的关系定义为参数αk的敏感度Sk(αk)。 Sk(αk)=(|ΔΡ|/Ρ)/(|Δαk|/αk) 堆石坝中主堆石区的堆石料占坝体体积的比例较大,是面板堆石坝主要的承载体,它的变形不仅是考核堆石坝性能的重要指标,而且还决定了混凝土面板的变形性态,因此对面板堆石坝的主堆石材料参数予以考察,进行灵敏度性分析。 在面板堆石坝非线性计算中,对在岩土工程测试上技术比较成熟可以较为准确量测的材料参数数据如容重γ,一般不必进行灵敏度分析;堆石材料内凝聚力参数c一般取为0;卸荷试验确定的参数nur、kur在数值计算中分别取加荷情况下n的1倍和k的1.2～3.0倍的值,因此对这3个参数也不进行灵敏度分析;对用于参数Δφ,一般它的试验值较小(5°～14°),因此对Δφ也不加考虑。 在对某面板堆石坝静力计算分析的基础上,利用变化的参数变量求大坝相应的力学响应。从而求出在各力学响应下的参数灵敏度,结果如表1所示。 通过对参数的敏感性分析得出,对位移的变化灵敏度较大的参数为Kb,φ0,K;对应力和应力水平影响较大的参数为Kb,φ0,K,Rf。 4 材料参数与耦合参数关系 由于面板坝设计中坝体各分区岩体的风化程度、级配和填筑干密度均有着一个相对统一的基本要求,因此对于堆石岩性相同的不同工程,其材料参数有着一定的可类比性。由于级配与模型参数关系的资料缺乏,因此,本文重点就Kb,φ0,K,Rf四个参数进行与密度、空隙率之间变化关系的统计分析。 就面板堆石坝的筑坝材料而言,我国应用较多的是沉积岩中的灰岩、砂岩以及火成岩中的花岗岩、凝灰岩、安山岩等5种,此外,随着面板堆石坝技术的不断发展,采用砂砾岩和软岩的工