Theoretical basis and application verification of scale effects of deformation characteristics of rockfill
-
摘要: 针对筑坝料变形特性的尺度效应,以接触变形理论为基础,从粒径尺度、试件尺度及外部影响因素角度,结合试验的方法开展了研究。建立了粒径尺度与归一化压缩模量的指数相关模型,论证了试件尺度径径比大于5时压缩试验的变形可靠性,揭示了侧限条件、干密度、边缘破碎及荷载等级5个外部因素对尺度效应的影响规律和内在机理,并分别建立了各因素与变形的相关模型。基于试验回溯现场过程,建立了综合考虑可破碎材料(堆石料)与不破碎材料(砂砾石)的多尺度统一修正模型。通过对大坝堆石料EB模型计算参数的修正及变形计算,验证了多尺度统一修正模型对变形计算的适用性与可靠性,较直接采用缩尺试验结果准确度提升了约20%,为解决坝体变形预测不准确提供了新的思考方法与解决路径。Abstract: In view of the scale effects of deformation characteristics of rockfill, based on the contact deformation theory, the research is carried out from the perspective of particle sizes, specimen sizes and external factors, combined with the test method. The exponential correlation model between the particle sizes and the normalized compression modulus is established. The deformation reliability of compression tests is demonstrated when the particle size ratio is greater than 5. The influence rules and internal mechanism of five external factors on the scale effects, namely, the lateral limit condition, the dry density, the edge breakage and the load level, are revealed, and the relevant models for each factor and deformation are established respectively. Based on the field experiments, a multi-scale unified correction model is established, which takes into account the crushed material (rockfill) and the non-crushed materials (gravel). The applicability and reliability of the multi-scale unified correction model for deformation calculation are verified by modifying the parameters and deformation calculation of the EB model for dam rockfill. The accuracy of the model is about 20% higher than that of the direct scale tests. This provides a new thinking method and solution path to solve the inaccuracy of dam deformation prediction.
-
Keywords:
- scale effect /
- compression modulus /
- influencing factor /
- correction model /
- dam deformation
-
0. 引言
随着筑坝材料的最大粒径越来越大,高坝筑坝料的最大粒径已经超过800 mm[1],常规土力学试验设备无法直接用于研究原筑坝料的变形特性[2-3]。为了解决试验粒径限制的问题,通常将原级配堆石材料进行缩尺得到替代原级配材料,然后对缩尺料进行试验,进而确定原型坝的变形特性[4]。然而,一些研究成果及监测资料都表明,缩尺前后粗粒料的力学特性有较大差异,缩尺后粗粒料的力学特性并不能完全反应原级配粗粒料的力学性质[5]。因此,如何确定原级配筑坝材料力学参数,是准确预测坝体变形的一个重要因素。
对此,国内外学者进行大量的大尺寸三轴试验。Marachi等[6]对不同最大粒径的堆石料进行三轴试验得出,大粒径较小粒径压缩性更大。Hunter等[7]积累了大量的不同粒径堆石料坝体的变形监测结果,提出堆石料的压缩模量是随粒径的增大而减小的。Varadarajan等[8-9]对堆石料和砂砾石料分别进行不同缩尺条件下的对比试验,得出堆石料模量随粒径增加而降低,而砂砾石料的模量则随粒径增大而增加。汪小刚[10]根据堆石料在压缩中可破碎和砂砾石料不可破碎的特征,进行了石膏和钢材的不同尺度下的压缩试验,得出与Varadarajan相同的规律。
不同文献所反映的缩尺试验结果可能与缩尺方法、密度、侧限条件、荷载及母岩等有关,但笔者认为将筑坝材料整体上分为可破碎的堆石料和不可破碎的砂砾石料分开研究是理清目前尺度效应研究的途径,同时为避免各种因素的交叉影响,本文基于接触变形理论、室内多尺度试验,研究不同粒径尺度、不同试样尺度、不同侧限条件、不同干密度、不同破碎条件、不同荷载等级等多种组合影响下变形模量的演化规律和内在机理,以此为基础建立多尺度变形参数的统一修正模型,为解决坝体变形预测不准确提供新的思考方法与解决途径。
1. 颗粒材料尺度效应理论分析
1.1 基本理论
(1)基本假定
以理想弹性力学作为计算理论基础;接触表面的摩擦力可忽略不计。
(2)变形方程
点接触物体受力后其接触表面为椭圆;线接触物体受力后其接触表面为矩形。符合变形连续条件。
(3)物理方程
服从弹性虎克定律,接触表面上应力的变化规律与接触体的应变成线性关系,应变最大的接触面中心压应力最大。根据接触表面压应力分布规律求得表面接触压力所组成的合力应等于外加载荷。
1.2 弹性体接触变形分析
设两圆球体的半径分别为R1,R2,开始时在公切平面上的O点相互接触,如图1所示。
在两球的子午面截线上,与z1和z2相距距离r处的两点M和N两点的距离为
Z1+Z2=12(1R1+1R2)。 (1) 设w1表示球体1面上的点M由于局部变形所产生的沿Z1轴方向的位移,w2表示球2面上的点N由于局部变形所产生的沿Z2轴方向的位移,两球的中心O1、O2彼此接近的距离为δ。由于对称性,由接触产生的压力q和位移w对于接触中心O都是轴对称的,计算出最终。由此可以计算出:
q0=3P2πa2, (2) a3=34⋅R1R2R1+R2(1−μ21E1+1−μ22E2)p, (3) δ3=916⋅R1+R2R1R2(1−μ21E1+1−μ22E2)2p2。 (4) 式中 q0为接触中心的压力最大;a是接触圆半径;E1,μ1,E2,μ2分别是两个弹性体弹性常数;p为压力分布函数;δ为变形量;R1,R2分别是两个弹性体半径。
由上式得出,最大接触压应力与载荷不是线性关系,而是与载荷的立方根成正比,这时候因为随着载荷的增加,接触面积也在增大,其结果使接触面上的最大压应力的增长较载荷的增长为慢。
当弹性体与平面接触变形时,则可以将R1=R,R1=∞,E1=E2=E,μ1=μ2=0.3代入式(2)~(4)求解,得到下式:
q0=0.3883√PE2R2, (5) a=1.1093√PRE, (6) δ=1.2313√P2E2R。 (7) 当计算两等直径弹性体接触时,则可将R1=R,R1=∞,E1=E2=E,μ1=μ2=0.3代入式(2)~(4)求解。
a=1.1093√PRE, (8) δ=1.2313√2P2E2R。 (9) 1.3 粒径尺度对压缩模量影响研究
首先研究固定试样尺寸,改变材料的最大粒径情况下对压缩模量的影响,从材料的粒径改变揭示压缩模量的变化规律。推导假定为:在相同加载试样尺寸内,材料从单个半径为R的弹性体逐渐减小其半径,分别为R/2,R/4,R/6,…,R/1000,如图2所示。在此过程中密度保持不变,空间内的孔隙率保持不变,荷载等级保持不变,无侧限条件,受荷载P的作用下,逐步分析粒径逐渐缩小所引起的弹性模量的变化。
利用1.2节理论推导结果,计算图2中不同粒径下的变形量,进而计算更小的最大粒径下变形量,形成归一化压缩模量。此处定义归一化压缩模量含义为:设未进行缩尺时试件的压缩模量初始值为1,缩尺后压缩模量与初始压缩模量之比为归一化压缩模量。
为更直观表达级配组构条件下粒径尺度对压缩模量的影响,将初始最大粒径设为1000 mm,对应的压缩模量归一化为1,并化为无量刚量,则随着粒径的减小,建立如图3所示关系。
建立粒径尺度与归一化压缩模量的关系如下式所示:
Egd=0.5691+0.2116e−ξ9.9617+0.2374e−ξ212.3621, (10) 式中,Egd为粒径尺度归一化模量,ξ为固定试样尺寸条件下的归一化粒径。
1.4 试件尺度对压缩模量影响研究
对于有级配组构情况下,首先确定组构规则,根据圆形几何组构规则,空隙部分首先用最大半径进行组构,组构后每个空间为三边形,在此三边空隙内继续组构,直至R/1000。组构弹性体的粒径符合指数递减规律,并形成悬浮型颗粒构架体系,同时构建最接近实际情况下的级配曲线,悬浮级配构建如图4所示,与之对应的级配曲线如图5。
为从材料的试样尺度角度揭示加载试验的面积改变对压缩模量的影响,逐渐增大试样尺度,分别为2R,4R,6R,…,1000R,求解在不同荷载加载面积下的压缩模量,如图6所示。
根据不同试样尺度条件下的归一化压缩模量,引入径径比的概念,既试样尺寸与最大粒径之比,建立径径比与压缩模量的关系如图7所示。可得出,当径径比大于5时,归一化压缩模量变稳定,从理论角度论证了压缩试验可靠性的必要条件。
2. 变形特性的外部影响因素研究
2.1 侧限条件对压缩模量的影响
对侧向不允许变形的压缩弹性体,取一个单元进行分析其弹性阶段受力变形。该单元受到三向应力
σx ,σy ,σz 作用,在z方向的压力作用下,试样中的竖向有效应力为σz ,又因为试样的受力条件属于轴对称问题,故水平向有效应力σx =σy ,由于土体的侧向变形被限制住,故εx =εy = 0,于是:εx=σxEs−μσyEs−μσzEs=0, (11) 式中,E,Es分别为无侧线和完全侧限条件下的模量,μ为组构材料的泊松比。μ取0.3时,建立完全侧向与无侧限模量的关系为
Es=1.346E。 (12) 由此可得完全侧向与无侧限模量随粒径整体的增大倍数,在μ取0.3时增大约35%。
2.2 干密度与压缩模量的关系研究
本研究推导成立的边界条件是从原型尺度到实验室尺度均保持同样的干密度或孔隙比,而在实际室内试验中,由于粒径的缩小通常表现为筑坝原型干密度大,而室内实验试件的干密度小,且很难使压缩试件的干密度达到现场的大小,这就为实现同密度或孔隙率条件下的室内试验造成一定的障碍。为解决这一障碍,找出室内试验密度或孔隙率与压缩模量的关系,采用此规律修正后即便在密度未达到原型密度下亦可进行实验并进行回溯。
对若干个工程的材料均进行不同干密度下的压缩试验,建立如图8所示趋势图。
建立如下式所示的关系。
f(ρd)=2.331+0.702−2.3311+(ρd1.996)29.891, (13) 式中,
f(ρd) 为受密度影响下试样的压缩模量,ρd 为干密度。2.3 颗粒边缘破碎影响分析
为分析堆石料破碎与粒径间的关系,在单一粒径情况下假定大小颗粒的粒径在特定比例下,某一应力状态中发生破碎,保持该应力状态不变,不断减小最大粒径。计算基本参数为:大粒径与小粒径之比为20,在缩小大粒径过程中,保持该比例不变。
图9不同粒径下的试样破碎变形计算。由式(2)~(4)可计算小粒径上的最大应力。在粒径不断缩小的过程中,单个小颗粒承受的外部应力会逐渐减小,化为等效破碎变形量,建立破碎量与粒径变化如图10所示。
该曲线采用函数表示为
δ=1.023e−D42.848。 (14) 2.4 荷载等级对压缩模量影响研究
为了解荷载等级对压缩模量的影响关系,在有填充条件下,分别计算半径为R/n的弹性体,在有填充作用下分别受荷载2p,3p,4p,…,10p作用下的压缩模量。压缩模量的增大倍数仅与荷载大小直接相关,绘制压缩模量与荷载等级趋势关系如图11所示。并可用下式:
f(p)=2.7354−1.9981e−Pi6.4893, (15) 式中,
f(p) 为荷载影响下的压缩模量放大倍数,pi为荷载等级倍数,为当前加压等级与初始加压等级的比值。3. 堆石料多尺度统一修正模型的建立
3.1 模型的提出
多尺度统一修正模型在于将缩尺材料的压缩模量求解为筑坝原型材料的压缩模量,在此过程中综合考虑缩尺过程、荷载等级、干密度和侧限状态、破碎效应的影响,对此提出涵盖此参数的统一修正模型公式为
F(ξ,p,ρd,μ,D)=f1(ξ)⋅f2(pi)⋅f3(ρd)⋅f4(μ)⋅f5(D)。 (16) 式中
f1(ξ) 为粒径影响下的的回溯函数;f2(pi) 为荷载等级影响下的修正函数;f3(ρd) 为干密度影响下的修正函数;f4(μ) 为泊松比影响下的侧限修正函数。f5(D) 为粒径影响下的破碎修正函数。根据回溯过程,建立粒径尺度和试样尺度影响下归一化压缩模量的关系,如图12所示。具体函数关系为
f1(ξ)=A1+A2−A11+(ξ1/A3)A4A1+A2−A11+(ξ0/A3)A4⋅B1+B2e−ξ2B3+B4e−ξ2B5B1+B2e−ξ1B3+B4e−ξ1B5⋅A1+A2−A11+(ξ3/A3)A4A1+A2−A11+(ξ2/A3)A4。 (17) 式中 Ai表示固定材料最大粒径下,试样尺寸对归一化压缩模量的影响系数;Bi表示固定试样尺寸条件下材料粒径变化与归一化压缩模量的系数;ξ为转换过程的粒径。
3.2 基于外部影响因素的统一修正
(1)泊松比影响下的侧限修正
侧限条件是筑坝料变形应力状态的具体描述,侧限所引起的压缩模量改变可通过泊松比表示,泊松比对压缩模量的影响关系,可用下式表示:
f4(μ)=1−2μ21−μ, (18) 式中,μ为筑坝料的泊松比,可通过压缩试验获取。
(2)干密度影响下的修正
由2.4节所提出的干密度对变形模量的影响模型,可得到在未达到原型密度下的回溯。如下式所示:
f3(ρd)=D1+D2−D11+(ρd/D3)D4, (19) 式中,Di(i=1,2,3,4,)为待定系数,
ρd 为干密度。(3)破碎效应影响修正
筑坝材料整体上可分为可压缩破碎的堆石料和不破碎的砾石料,当采用堆石料筑坝时,往往表现出粒径越大压缩模量越低的情况,由此引入破碎效应的修正,如下式。当采用砂砾石料时不发生材料颗粒破碎,可不考虑破碎效应。
f5(D)=αeD/r, (20) 式中,D为最大粒径,
α ,r 为待定系数。(4)荷载等级影响下的修正函数
对于筑坝料,其压缩模量与试验施加的荷载应力大小有关,表现为荷载应力水平越大压缩模量越大,呈非线性正相关,具体如下式表示:
f2(pi)=C1+C2eC3ρd, (21) 式中,Ci为受荷载等级影响下的压缩模量关系的非线性表达系数。
综合以上各影响因素,选取实际试验中所涉及到的情况进行选用。
3.3 多尺度统一修正模型的应用
以河口村大坝为例,坝体为混凝土面板堆石坝,最大坝高122.5 m,面板堆石坝设计主要考虑的因素是大坝的应力变形状态,而大坝应力变形分析的主要手段是有限元数值分析。将室内缩尺三轴试验所得E-B模型参数进行的坝体沉降计算与经统一修正模型修正后的模型参数计算结果进行对比,并与工程监测值相校核。
主要考虑主堆石料计算参数的影响,主堆石区填筑材料为级配碎石,室内三轴试验所得主应力差与轴向应变如图13所示。根据该三轴试验结果计算主堆石材料的E-B模型参数如表1所示。采用该参数计算所得最大竖向沉降竣工时61.2 cm,蓄水期68.1 cm。
表 1 修正前后河口村主堆石材料E-B模型参数Table 1. E-B model parameters of main rockfill materials in Hekou Village分类 干重度/(kN·m-3) K n Rf Kur c/kPa φ0/(°) Δφ/(°) Kb m 修正前 20.7 1660 0.21 0.85 3320 0 54 10.6 380 0.14 修正后 20.7 1312 0.25 0.82 3320 0 54 10.6 581 0.14 本文的模型最要考虑变形参数的影响,E-B模型中与变形相关的过程量为Ei,根据三轴试验过程和材料属性,采用统一修正模型进行尺度及各影响因素的影响系数计算,如表2所示。
表 2 河口村主堆石材料修正计算Table 2. Correction calculation of main rockfill materials in Hekou Village项目 f1(ξ) f2(μ) f3(ρd) f4(D) f5(p) 综合影响F(·) 修正值 1.19 1.0 1.47 0.56 1.0 0.84 采用表2中修正后的E-B模型参数,计算总体沉降量对比如表3所示。
表 3 河口村坝体沉降计算对比Table 3. Comparison of settlement calculation of Hekou Village Dam项目 计算最大沉降/cm 提升准确度/% 竣工 蓄水 室内缩尺试验 61.2 68.1 监测数据 78.9 109.7 17.43~28.89 统一修正模型 75.6 99.7 4. 结论
针对筑坝料的变形特性的尺寸效应,以接触变形理论为基础,从粒径尺度、试件尺度及外部影响因素角度,结合试验的方法进行了研究,主要得到以下结论:
(1)基于弹性力学接触变形理论,建立了粒径尺度与归一化压缩模量的指数相关模型,论证了试件尺度中径径比大于5时压缩试验的变形可靠性。
(2)揭示了与尺度效应相关的5个外部因素,侧限条件、干密度、边缘破碎及荷载等级等多种组合对变形特性的影响规律和内在机理,并分别建立了各因素与变形的相关模式。
(3)基于试验回溯现场过程,建立了综合考虑缩尺过程、侧限条件、干密度、边缘破碎及荷载等级因素下的多尺度变形参数统一修正模型。
(4)通过对大坝堆石料变形E-B模型计算参数的修正及变形计算,验证了多尺度变形参数统一修正模型的适用性与可靠性,较直接采用缩尺试验结果准确度提升了约20%。
-
表 1 修正前后河口村主堆石材料E-B模型参数
Table 1 E-B model parameters of main rockfill materials in Hekou Village
分类 干重度/(kN·m-3) K n Rf Kur c/kPa φ0/(°) Δφ/(°) Kb m 修正前 20.7 1660 0.21 0.85 3320 0 54 10.6 380 0.14 修正后 20.7 1312 0.25 0.82 3320 0 54 10.6 581 0.14 表 2 河口村主堆石材料修正计算
Table 2 Correction calculation of main rockfill materials in Hekou Village
项目 f1(ξ) f2(μ) f3(ρd) f4(D) f5(p) 综合影响F(·) 修正值 1.19 1.0 1.47 0.56 1.0 0.84 表 3 河口村坝体沉降计算对比
Table 3 Comparison of settlement calculation of Hekou Village Dam
项目 计算最大沉降/cm 提升准确度/% 竣工 蓄水 室内缩尺试验 61.2 68.1 监测数据 78.9 109.7 17.43~28.89 统一修正模型 75.6 99.7 -
[1] IJHD. 2018 World Atlas & Industry Guide[R]. Wallington: International Journal of Hydropower and Dams, 2018.
[2] MA H Q, CHI F D. Technical progress on researches for the safety of high CFRDs[J]. Engineering, 2016, 2(4): 332-339. doi: 10.3969/j.issn.1009-444X.2016.04.011
[3] 秦红玉, 刘汉龙, 高玉峰, 等. 粗粒料强度和变形的大型三轴试验研究[J]. 岩土力学, 2004, 25(10): 1575-1580. doi: 10.3969/j.issn.1000-7598.2004.10.013 QIN Hong-yu, LIU Han-long, GAO Yu-feng, et al. Research on strength and deformation behavior of coarse aggregates based on large-scale triaxial tests[J]. Rock and Soil Mechanics, 2004, 25(10): 1575-1580. (in Chinese) doi: 10.3969/j.issn.1000-7598.2004.10.013
[4] 徐志华, 孙大伟, 张国栋. 堆石料应力-应变特性大型三轴试验研究[J]. 岩土力学, 2017, 38(6): 1565-1572. XU Zhi-hua, SUN Da-wei, ZHANG Guo-dong. Study on stress-strain behavior of rockfill using large-scale triaxial tests[J]. Rock and Soil Mechanics, 2017, 38(6): 1565-1572. (in Chinese)
[5] KONG X, LIU J, ZOU D, et al. Stress-dilatancy relationship of Zipingpu gravel under cyclic loading in triaxial stress states[J]. International Journal of Geomechanics, 2016, 16(4): 4016001. doi: 10.1061/(ASCE)GM.1943-5622.0000584
[6] MARACHI N D, CHAN C K, SEED H B. Evaluation of properties of rockfill materials[J]. Journal of Soil Mechanics & Foundations Div, 1972, 98(1): 95-114.
[7] HUNTER G, FELL R. Rockfill modulus and settlement of concrete face rockfill dams[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2003, 129(10): 909-917. doi: 10.1061/(ASCE)1090-0241(2003)129:10(909)
[8] VARADARAJAN A, SHARMA K G, ABBAS S M, et al. Constitutive model for rockfill materials and determination of material constants[J]. International Journal of Geomechanics, 2006, 6(4): 226-237. doi: 10.1061/(ASCE)1532-3641(2006)6:4(226)
[9] VARADARAJAN A, SHARMA K G, VENKATACHALAM K, et al. Testing and modeling two rockfill materials[J]. Journal of Geotechnical &Geoenvironmental Engineering, 2003, 129(3): 206-218.
[10] 汪小刚. 高土石坝几个问题探讨[J]. 岩土工程学报, 2018, 40(2): 203-222. WANG Xiao-gang. Discussion on some problems observed in high earth-rockfill dams[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(2): 203-222. (in Chinese)