《武汉工程大学学报》 2013年09期
18-23
出版日期：2013-10-10
ISSN:1674-2869
CN:42-1779/TQ
含水率对花岗岩残积土邓肯\|张模型参数影响
0引言 花岗岩残积土是花岗岩经风化后残留在原地的碎屑物，具有很强的结构性，在原位状态或土的微结构未受损伤时呈现高强度、中低压缩性［1］.但花岗岩残积土经过浸水饱和后，承载力降低，压缩性增大，在具有临空面的浸水条件下，花岗岩残积土会因软化崩解而坍塌［2］.事实上，含水率对所有岩土的力学性质都有很大影响，因此也是研究热点之一，如李险峰［3］通过三轴排水剪切试验，研究了含水率对滑坡碎屑土力学特性的影响.土的本构模型是研究岩土力学性质的核心问题，显然，研究含水率对花岗岩残积土本构模型参数的影响具有理论与工程实际意义.近几十年来，国内外学者提出描述土体应力\|应变关系的本构模型多达数百个，这些模型基本上可分为弹性、弹塑性及其他各种新型模型.邓肯\|张模型因其概念清晰，模型参数明确易得，自提出后40多年来得到了不断完善，被广泛应用于工程实践中［4］. 本研究进行了一系列不同含水率条件下花岗岩残积土的固结排水三轴试验，分析不同围压下花岗岩残积土的应力\|应变曲线特点，通过计算得到不同含水率花岗岩残积土的邓肯\|张模型参数，分析并归纳含水率对花岗岩残积土邓肯\|张模型参数影响的基本规律.1 模型简介 1963年，康纳(Kondner)［5］在大量三轴试验基础上，提出围压为常数的常规三轴固结试验的加工硬化型应力\|应变关系可近似地用双曲线来表示，即：σ1-σ3=ε1a+bε1(1)式(1)中，a、b为试验常数，ε1为轴向应变，σ3为围压，σ1为轴向主应力. 可将式(1)改写为：ε1σ1-σ3=a+bε1(2)式(2)是一直线方程，a、b分别为该直线的截距、斜率，可通过试验数据拟合直线获得. 1970年，邓肯(Duncan)和张(Chang)等［6］根据式(1)提出了包括切线模量Et和切线泊松比μt的邓肯\|张本构方程式，简称邓肯\|张模型E\|μ模型.其中切线模量Et表达式为：Et=Ei1-Rf(σ1-σ3)(1-sinφ)2ccosφ+2σ3sinφ2(3)Ei=KPa(σ3Pa)n (4)式(3)(4)中：Ei为初始切线模量(kPa)，c为土的粘聚力(kPa)，φ为土的内摩擦角(°)，Rf为破坏应力比，K、n为试验常数，Pa为大气压力(取101.4 kPa). Ei、Rf与a、b的关系为：Ei=1a(5)Rf=b(σ1-σ3)f(6)式(5)(6)中，(σ1-σ3)f为土样破坏时的偏应力. 将式(4)进行坐标变换，可变为lgEi~lgσ3直线，则该直线斜率为n，当lgσ3=1时的lgEi值即为lgK.2试验研究2.1试验土样 试验土样取自福州晋安区鹤林新城二区某地块，取土深度为12.3~16.8 m，土样呈褐黄、肉红色，湿\|稍湿，硬塑状态，土样的物理性质指标如下：天然密度ρ=1.84 g/cm3，天然含水率ω=30.1%，土粒相对密度ds=2.72，天然孔隙比e=0.976，液限ωL=41.1%，塑限ωP=27.3%，塑性指数Ip=13.8.经筛分，土样各粒组的分布情况见表1.土样击实试验的最优含水率ωop=17.4%.表1花岗岩残积土粒组分布情况Table 1Granularity distribution of the granite residual soil粒组/mm10~55~22~11~0.50.5~0.250.25~0.1<0.1w/%1.417.47.39.88.210.645.3第9期吴能森，等：含水率对花岗岩残积土邓肯\|张模型参数影响武汉工程大学学报第35卷2.2试验方案 试验仪器为SLB\|1型应力\|应变控制式三轴剪切渗透试验仪，采用固结排水三轴压缩试验(CD试验) ［7］，花岗岩残积土经过风干、碾碎、过孔径为2 mm筛后，测定风干后土样的含水率，然后按照设计的含水率制作试样，试样按击实法制备扰动土制样方法制作，试样直径39.1 mm，高80 mm.试样含水率以最优含水率为基准，上下相差4%左右，经取整分别为13%，17%，21%.试验固结围压σ3分别取100、200、300 kPa，当体积读数变化微小或长时间保持不变时认为试样完成固结.固结完成后，在围压σ3保持不变的情况下，采用应变控制方式，以0.012%/min剪切速率增加轴向压力进行排水剪切，直至试样轴向应变值达到15%时停止试验，试验过程中排水阀始终保持开启状态.自动数据采集系统能够在试验过程中采集试验数据并绘制相关曲线.2.3试验结果及分析2.3.1试样变形特征及应力\|应变特性各含水率试样在三种不同围压下，均在轴向应变达到15%时呈比较明显的鼓状变形，但尚未发生剪切破坏.各试样的变形情况见图1.图2为各含水率试样的主应力差(σ1-σ3)与轴向应变ε1关系曲线，可见，不同含水率花岗岩残积土试样在固结排水条件下得到的应力\|应变曲线形状相似，未出现明显的峰值，均呈非线性的加工硬化型曲线，与邓肯\|张模型的条件式(1)基本吻合.2.3.2邓肯\|张模型参数及分析根据三种围压σ3作用下试样破坏时的偏应力(σ1-σ3)f 试验数据，可绘制出不同含水率下三组试样的强度包线，见图3，得到不同含水率的花岗岩残积土抗剪强度指标c、φ 值.图1试样破坏时变形特征Fig.1Deformation of destroyed test samples图2(σ1-σ3)\| ε1关系曲线Fig.2The curves of (σ1-σ3)\| ε1图3三轴剪切摩尔强度包线Fig.3The Mohr stress circles of triaxial compression test将(σ1-σ3)\| ε1试验数据按式(2)处理并绘制直线，可得试验常数a、b，进而根据式(5)、(6)得到初始切线模量Ei和破坏应力比Rf.再利用式(4)的线性化变换和绘图，得到试验常数K、n. 以上6个邓肯\|张模型参数及主要试验数据见表2.依表2绘制得到各参数随含水率的变化情况，见图4~图8.表2不同含水率花岗岩残积土邓肯\|张模型参数表Table 2Parameters of Duncan\|Chang model under different water contentω/%σ3/kPa(σ1-σ3)f/kPaEi/kPac/kPaφ/(°)KnRf100207.299 090.9113200325.3310 309.2824.8824.187.090.220.908300455.8111 363.64100176.587 812.5017200320.349 708.7418.5123.977.630.310.919300452.5110 869.57100153.526 024.1021200294.229 259.268.6923.566.070.450.952300421.4410 416.67图4粘聚力c与内摩擦角φ随含水率变化Fig.4Cohesion c and internal friction angle φ change with the water content图5初始切线模量Ei随含水率变化Fig.5Initial tangent modulus Ei changes with the water content图6参数K随含水率变化Fig.6Parameter K changes with the water content图7参数n随含水率变化Fig.7Parameter n changes with the water content图8参数Rf随含水率变化Fig.8Parameter Rf changes with the water content由图4可见，含水率对花岗岩残积土抗剪强度指标c、φ值的影响差异显著，虽然c、φ值都随着含水率的增大而减小，但两者的减小幅度相差很大，含水率对c值影响很大而对φ值影响相对微小，该结果与文献［1］的统计分析结论一致.由表1可得：当含水率ω=13%~17%，ω值每增加1%，c值下降6.40%，而φ值仅下降0.21%；当含水率ω=17%~21%，ω值每增加1%，c值下降13.26%，而φ值仅下降0.42%.在同等条件下，c值降幅达φ值30倍以上.同时，含水率愈高对c、φ值的影响愈大：以最优含水率ωop为界，含水率大于ωop时的c、φ值降幅大约是含水率小于ωop时的2倍. 由图5及表2可知，花岗岩残积土初始变形模量Ei 随着含水率和围压σ3的变化呈现一定规律：在围压一定时，Ei 随着含水率的增大而减小，其中围压较小(σ3=100 kPa)时，Ei 随含水率增大而减小的幅度较大，而且以最优含水率ωop为界，含水率大于ωop时的Ei值减幅大约是含水率小于ωop时的1.4倍；而围压较大(σ3=200，300 kPa)时，Ei 随含水率增大而减小的幅度较小，而且基本呈线性关系.在含水率不变时，Ei值随着围压σ3的增加而增大，特别是围压较小时，围压的增加对Ei值的增大影响更显著. 同样地，其余参数也随着含水率的变化而变化.由图6~图8及表2可知：参数K随着含水率的增大而减小，而且以最优含水率ωop为界，含水率大于ωop时的K值减幅大约是含水率小于ωop时的1.2倍；参数n、Rf均随着含水率的增大而增大：对参数n，含水率大于ωop时的n值增幅大约是含水率小于ωop时的1.56倍；对参数Rf，含水率大于ωop时的Rf值增幅大约是含水率小于ωop时的3倍. 说明一点：参数Rf值一般在0.75~1.00之间［8］，表2的Rf值均在0.90以上，表明实验取轴向应变ε1=15%对应的偏应力作为破坏偏应力(σ1-σ3)f，与极限偏应力(σ1-σ3)u接近，而此实试样除了继续发生鼓胀变形外，并未出现破坏剪切面，体现了土样的应变硬化特性和应变控制破坏的特点.3结语 a. 在不同围压和不同含水率条件下，花岗岩残积土击实试样在三向应力作用下均呈鼓胀变形特征，其固结排水应力应变曲线均为加工硬化型，与邓肯\|张模型的条件式基本吻合. b. 花岗岩残积土的邓肯\|张模型参数随含水率变化的基本规律是：抗剪强度指标c、φ值，初始变形模量Ei 及参数K均随着含水率的增大而减小；参数n、Rf随含水率的增大而增大. c. 除了内摩擦角φ之外，含水率对花岗岩残积土的邓肯\|张模型参数整体上影响较大，其中粘聚力c受影响显著，参数n和低围压时的初始变形模量Ei 受影响较显著.各参数与含水率关系曲线均在最优含水率ωop处出现拐点，而且当含水率大于ωop时，其对各参数的影响不同程度地加大. d. 试样制作前筛除了2 mm以上粗颗粒(占18.8%)，必会影响试验结果，但由于土中水主要与细粒土发生作用，因此预计该影响主要表现在参数的大小上，不会改变其应力应变特性及参数的变化规律.致谢 本研究得到福建省科技厅、福建农林大学等部门资助，在此表示感谢.同时也感谢福建省现代工程勘察院在钻探取样方面的支持与帮助，感谢福建农林大学交通与土木工程学院给予本研究的帮助、关心与支持！