遇水湿化对堆石体颗粒破碎和压缩特性的影响研究

马刚; 王渊; 程家林; 叶晓峰; 邹宇雄; 周伟; 刘其文

doi:10.11779/CJGE202109008

遇水湿化对堆石体颗粒破碎和压缩特性的影响研究 English Version

马刚^{1, 2,},
王渊^{1, 2},
程家林^{1, 2},
叶晓峰^{1, 2, 3},
邹宇雄^{1, 2},
周伟^{1, 2, ,},
刘其文⁴

1.
武汉大学水资源与水电工程科学国家重点实验室,湖北　武汉　430072
2.
武汉大学水工岩石力学教育部重点实验室,湖北　武汉　430072
3.
中铁第四勘察设计院集团有限公司,湖北　武汉　430063
4.
贵州省水利水电勘测设计研究院,贵州　贵阳　550002

基金项目:

国家重点研发计划项目课题项目 2018YFC1508503

国家自然科学基金项目 51825905

国家自然科学基金项目 U1865204

华能集团科技项目 HNKJ18-H26

贵州省科技重大专项项目黔科合重大专项字[2017]3005-2号

详细信息

作者简介:
马刚(1985— ),男,副教授,主要从事结构工程数值模拟及仿真方面的研究工作。E-mail：magang630@whu.edu.cn

通讯作者:
周伟, E-mail：zw_mxx@whu.edu.cn

中图分类号: TV41
计量
- 文章访问数: 312
- HTML全文浏览量: 22
- PDF下载量: 170
出版历程
- 收稿日期: 2020-12-16
- 网络出版日期: 2022-12-02
- 刊出日期: 2021-08-31

Influences of wetting on crushing and compression characteristics of rockfill particles

MA Gang^{1, 2,},
WANG Yuan^{1, 2},
CHENG Jia-lin^{1, 2},
YE Xiao-feng^{1, 2, 3},
ZOU Yu-xiong^{1, 2},
ZHOU Wei^{1, 2, ,},
LIU Qi-wen⁴

1.
State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan 430072, China
2.
Key Laboratory of Rock Mechanics in Hydraulic Structural Engineering of Ministry of Education, Wuhan University, Wuhan 430072, China
3.
China Railway Siyuan Survey and Design Group Co., Ltd., Wuhan 430063, China
4.
Guizhou Survey & Design Research Institute for Water Resources and Hydropower, Guiyang 550002, China

摘要

摘要: 堆石体遇水后产生的湿化变形会影响堆石坝的安全稳定运行,一直是工程界和学术界关注的焦点问题,但目前对湿化变形细观机理的认识还不够充分。采用室内试验、电镜扫描、连续–离散耦合分析方法,进行了湿化对堆石体破碎特性和压缩特性的影响研究。干燥和饱和状态的单颗粒破碎试验和SEM电镜扫描表明,遇水饱和导致颗粒的岩石矿物成分和微观结构发生变化,其孔隙、微裂纹等缺陷增多,颗粒破碎强度和表面摩擦系数降低。率定了干燥与饱和颗粒的细观参数,开展了干燥和饱和状态的单颗粒破碎和颗粒集合体单轴压缩的数值试验。遇水湿化对堆石体有润滑、软化和劣化作用,使堆石颗粒的摩擦系数、变形模量和破碎强度均出现了不同程度的降低。湿化使颗粒破碎强度降低了22.6%,使颗粒间摩擦系数降低了11.7%,使颗粒集合体的压缩模量降低了13.1%。破碎强度降低导致的颗粒破碎是堆石体产生湿化变形的最主要原因。
- 堆石体 /
- 湿化变形 /
- 颗粒破碎 /
- 连续–离散耦合分析方法
Abstract: Wetting deformation of rockfills affects the safety and stability of rockfill dams, which has always been the focus of engineering and academic circles. However, the mechanism of wetting deformation is still not fully understood at present. The laboratory tests, scanning electron microscope and combined finite-discrete element method are used to study the effects of wetting on the crushing and compression characteristics of the rockfills. The single particle crushing tests and SEM of dry and saturated particles show that the wetting leads to the change of rock mineral composition and microstructure of particles. The pores, microcracks and other defects increase owing to the wetting, which reduce the crushing strength of particles and surface friction coefficient. The mesoscopic parameters of the dry and saturated particles are determined, and the numerical tests on the single particle crushing and the uniaxial compression of the particle assembly in the dry and saturated states are conducted. The wetting has effects of lubricating, softening and degrading on the rockfills, which make the friction coefficient, deformation modulus and crushing strength of the rockfill particles decrease to different degrees. The wetting reduces the crushing strength of the particles by 22.6%, the friction coefficient between the particles by 11.7%, and the compression modulus of the particle assembly by 13.1%. The main reason of wetting deformation of the rockfills is the particle breakage caused by the reduction of crushing strength.
- rockfill /
- wetting deformation /
- particle breakage /
- combined finite-discrete element method

HTML全文

0. 引言

土工格室等三维材料由于具有较为优良的工程性质而被较广泛的应用于双向增强复合地基中,深入分析格室体的变形情况可以得知复合地基的工作状况,从而得出桩土应力比等重要参数。在实际工程中,土工格室不仅受到来自路堤的竖向压力、桩与土的支持力,其上下表面还会与垫层材料发生摩擦,因此分析起来较为复杂,而目前的分析方法仍存在相应的局限性,故有必要进行进一步的深入探讨（图1）。

图 1 路堤-土工格室加筋体-桩土加固区整体示意图

Figure 1. Overall diagram of embankment-geocell-reinforced body-pile-soil reinforcement area

下载: 全尺寸图片幻灯片

通过室内试验及数值试验研究土工格室加筋体性能是较为常见的方法：周亚梅等^[1]通过定速度压缩试验分析了不同格室形状及侧限形式对单个土工格室承载变形特性影响；郑超毅等^[2]通过多组模型试验证明了土工格室在提高地基承载力,减少地基沉降方面的作用；高昂等^[3]、邓鹏等^[4]孙州等^[5]分别通过循环加载、大比尺直剪试验及对比实验,研究了不同加筋形式、格室高度、埋深及焊距等参数在其提高地基承载性能方面的影响。侯娟等^[6]通过建立单个高强土工格室加筋和未加筋地基的有限元模型对比分析了土工格室对土体有格室侧壁的摩擦力及环箍约束作用；汪海年等^[7]利用离散连续耦合算法分析了格室体高度及焊炬对加筋效果的影响。

理论研究方面,弹性地基梁板理论是分析土工格室变形的常用方法,其中传统地基梁理论较为简便,但其忽略格室体水平摩阻效应会夸大竖向变形,因而部分学者在其基础上进行了一定的改进。张福海等^[8]、张玲等^[9]在分析格室变形时分别利用修正的双参数地基模型与假定的摩阻力分布模式来考虑摩阻力的影响,但由于其只计算摩阻力对剪力的影响而忽略了对弯矩的影响,使得计算结果与传统方法十分接近,陈仁朋等^[10]通过进一步研究分析,得出了桩土应力比、沉降受筋材抗拉模量变化的影响不大,这是由于格室体自身具有较大的刚度,在路堤荷载下只发生有限的挠曲变形,其拉应变对竖向附加应力的贡献不大。于是,为弥补上述不足,张玲等^[11-13]、赵明华等^[14-16]、马缤辉等^[17]在地基梁界面设置了水平弹簧,进一步优化了变形计算结果,其中张玲等^[13]注意到桩与土在变形刚度方面具有较大差异,采用具有不同刚度的弹簧体系表征桩与土的这种明显区别,为得到格室体变形曲线,假定桩土交界处土工格室变形协调,在二维情况下实现了桩土应力比的解析解答。然而一般情况下,当采用地基梁方法计算此类问题时,由于难以考虑桩土刚度差异,也通常忽略土拱效应影响,会出现计算结果较为保守的情况。

于是,又有部分学者采用薄板理论来计算格室体的变形,从而避免上述地基梁方法存在的不足。饶为国等^[18-20]、谭慧明等^[21]及张军等^[22]、郑俊杰等^[23]分别基于矩形薄板理论分析了土工格室的变形情况,其中郑俊杰等提出了薄板复合弹性模量的概念,但鉴于此方法对于布桩方式的适用性存在一定局限,因而不能用来分析梅花形布桩方式时的情况。赵明华等^[24]基于对上述学者的研究成果,采用小挠度弹性圆薄板模型模拟土工格室加筋体,本文在此基础上引入土拱效应的影响并构造非线性代数方程组,建立了“路堤–格室垫层–桩土加固区”共同作用模型,通过迭代求解方法得到土工格室变形情况并进一步得出计算桩土应力比、沉降及桩土差异沉降的新方法,通过实例验证和参数分析与前人方法进行了结果对比,证明了本文方法的合理性。

1. 土拱模型建立

1.1 内外土柱界面摩阻力分布

图2为路堤“土柱模型”,参考文献[25],假定其侧摩阻力发挥系数在等沉面与路堤底部之间由0变化到1,因此,土柱界面摩阻力分布为

图 2 路堤土拱分析模型

Figure 2. Analysis model for soil arch

下载: 全尺寸图片幻灯片

$τ_{e} (z) = \frac{δ}{Δ s} τ_{eu} (z)$ 。

(1)

式中 τ_e为z截面处内外土柱交界处的摩阻力；δ为z截面处内、外土柱的相对位移；Δs为路堤底部桩土差异沉降；τ_eu为z截面处内、外土柱之间的极限摩阻力。

$τ_{eu} (z) = f K_{e} (z) p_{es} = f K_{e} (z) (γ_{e} z - σ_{es})$ ,

(2)

式中,f为内外土柱之间的摩擦系数,f=tanϕ_e,K_e为z截面处内、外土柱之间土压力系数,γ_e为填土体重度,p_es为外土柱截面总应力,σ_es为z截面处因荷载转移外土柱所减少的应力,自重应力扣除σ_es即为外土柱实际应力。

在z= $H - H_{e}$ 处,K_e=K₀,而在z=H处,K_e=K_p,将土柱间相对位移δ与桩土差异沉降Δs的关系表示为

$K_{e} (z) = (K_{p} - K_{0}) \frac{δ (z)}{Δ s} + K_{0}$ ,

(3)

由土力学相关知识可知K_p=tan²(45+φ_e/2),K₀=1-sinφ_e,ϕ_e为路堤土内摩擦角。

联立上述3式,有

$τ_{e} (z) = \frac{δ (z)}{Δ s} f [(K_{p} - K_{0}) \frac{δ (z)}{Δ s} + K_{0}] (γ_{e} z - σ_{es} (z))$ ,

(4)

式中,f为内外土柱之间的摩擦系数,f=tanϕ_e。

1.2 受力机制分析

如图3所示,分析dz厚度的内土柱在z方向受力,有

$σ_{ep} A_{p} + γ_{e} z + γ_{e} d z + τ_{e} U_{p} d z = (σ_{ep} + d σ_{ep}) A_{p} + γ_{e} (z + d z),$

(5)

图 3 格室体界面位移与摩阻力的关系

Figure 3. Relationship between displacement at interface of geocell and friction force

下载: 全尺寸图片幻灯片

式中 σ_ep为z截面处内土柱附加应力,即因荷载转移内土柱所增加的应力；A_p为桩身（桩帽）横截面面积,A_p=πd_p²/4,其中,d_p为桩体（桩帽）直径；U_p为横截面周长,U_p=πd_p。

上式进一步化简为

$\frac{d σ_{ep}}{d z} = \frac{4}{d_{p}} τ_{e}$ 。

(6)

将式（4）代入式（6）,有

$\frac{d σ_{e}_{p}}{d z} = [\frac{4 f (K_{p} - K_{0})}{d_{p} Δ s^{2}} δ^{2} + \frac{4 f K_{0}}{d_{p} Δ s} δ] (γ_{e} z - σ_{e}_{s})$ 。

(7)

同时根据z截面处总附加应力为0,有

$(1 - m) σ_{es} - m σ_{ep} = 0$ ,

(8)

式中,m为置换率,m=A_p/A_e。

联立式（7）,（8）得

$\frac{d σ_{e}_{p}}{d z} = [\frac{4 f (K_{p} - K_{0})}{d_{p} Δ s^{2}} δ^{2} + \frac{4 f K_{0}}{d_{p} Δ s} δ] (γ_{e} z - \frac{m}{1 - m} σ_{e}_{p})$ 。

(9)

截面的差异压缩量δ：

$δ = \int_{0}^{z} \frac{σ_{ep}}{E_{e}} d z + \int_{0}^{z} \frac{σ_{es}}{E_{e}^{'}} d z$ ,

(10)

式中,Ε_e和Ε_e'分别为路堤填土的压缩模量与回弹模量,此处令二者相等来使计算更加简便。式中第一部分为内、外土柱横截面竖向应力增加产生压缩变形,式中第二部分为内、外土柱横截面竖向应力减小而产生的回弹变形。

对式（10）进一步微分,得

$\frac{d δ}{d z} = \frac{σ_{ep}}{E_{e}} + \frac{σ_{es}}{E_{e}}$ 。

(11)

与式（8）联立得

$\frac{d z}{d δ} = \frac{(1 - m) E_{e}}{σ_{ep}}$ 。

(12)

结合式（9）与式（12）并进行整理,而后采用分离变量法可最终得到

$\begin{array}{l} - \frac{{(1 - m)}^{2}}{m^{2}} γ_{e} z \ln (\frac{1 - m}{m} γ_{e} z - σ_{ep}) - \frac{1 - m}{m} σ_{ep} \\ = \frac{4 (1 - m) f (K_{p} - K_{0}) E_{e}}{3 d_{p} Δ s^{2}} δ^{3} + \frac{2 (1 - m) f K_{0} E_{e}}{d_{p} Δ s} δ^{2} + C_{e} (z), \end{array}$

(13)

式中,C_e(z)是关于z的待定函数。此式表示内土柱附加应力σ_ep与土柱截面差异压缩量δ之间的关系。

1.3 差异沉降与桩土应力比的关系

根据实际情况,当δ=0时

$σ_{ep} = 0, σ_{es} = 0$ 。

(14)

将这一条件代入式（13）中,可得

$C_{e} (z) = - \frac{{(1 - m)}^{2}}{m^{2}} γ_{e} z \ln (\frac{1 - m}{m} γ_{e} z)$ 。

(15)

将式（15）代入式（13）,并考虑路堤底面处,即z=H时,差异压缩量δ=Δs,又有

$p_{ep} = γ_{e} H + σ_{ep}$ ,

(16)

式中,p_ep是桩顶应力,H为路堤高度。

整理可得

$\begin{array}{l} Δ s = \frac{3 d_{p}}{4 f K_{p} E_{e} + 2 f K_{0} E_{e}} \cdot \\ [\frac{1 - m}{m^{2}} γ_{e} H \ln (\frac{(1 - m) γ_{e} H}{γ_{e} H - m p_{ep}}) - \frac{1}{m} (p_{ep} - γ_{e} H)] \end{array}$ 。

(17)

又因为路堤底面桩土应力比n_e与p_ep的关系为

$p_{ep} = \frac{n_{e} γ_{e} H}{n m - m + 1}$ 。

(18)

将式（18）代入式（17）,可得

$\begin{array}{l} Δ s = \frac{3 d_{p}}{4 f K_{p} E_{e} + 2 f K_{0} E_{e}} \cdot \\ [\frac{1 - m}{m^{2}} γ_{e} H \ln (n_{e} m - m + 1) - \frac{γ_{e} H}{m} (\frac{n_{e} - n_{e} m + m - 1}{n_{e} m - m + 1})], \end{array}$

(19)

由此表示出桩土应力比n_e与差异沉降Δs的关系。

2. 考虑土拱效应的土工格室加筋体变形分析

参考文献[24]已有研究,未考虑土拱效应时,桩顶及桩间土控制方程分别建立如下所示：

$D \nabla^{4} w_{p} - \frac{(k_{xp,u} + k_{xp,d}) h^{2}}{4} \nabla^{2} w_{p} + p_{ep} = q_{p}$ ,

(20)

$D \nabla^{4} w_{s} - \frac{(k_{xs,u} + k_{xs,d}) h^{2}}{4} \nabla^{2} w_{s} + p_{es} = q_{s}$ ,

(21)

式中,D为薄板的弯曲刚度

$D = \frac{E h^{3}}{12 (1 - ν^{2})}$ ,

(22)

式中,E为板的弹性模量, $ν$ 为板的泊松比, $\nabla^{2}$ 为拉普拉斯算子

$\nabla^{2} = \frac{d^{2}}{d ρ^{2}} + \frac{1}{ρ} \frac{d}{d ρ}$ ,

(23)

w_p和w_s分别表示格室体在桩顶和桩间土部分的挠曲函数,由Winkler假定,p_es=k_s w_s,k_s为桩间土基床系数。k_xp,u与k_xp,d分别为桩顶范围内格室体上、下界面水平摩阻系数,k_xs,u与k_xs,d分别为桩间土范围内格室体上、下界面水平摩阻系数,h为格室体厚度。q_p和q_s分别表示桩顶范围及桩间土范围内所受上部荷载。

桩顶范围内格室体挠度求解为

$w_{p} (ρ) = C_{1} I_{0} (\sqrt{λ_{1}} ρ) + C_{2} - \frac{(q_{p} - p_{p}) ρ^{2}}{(k_{xp,u} + k_{xp,d}) h^{2}}$ ,

(24)

$λ_{1} = \frac{(k_{xp,u} + k_{xp,d}) h^{2}}{4 D}$ ,

(25)

式中,I_N为第一类N阶虚变量Bessel函数,在此解中N=0；C₁,C₂为待定系数,ρ表示距离圆板中心的水平距离。

$Δ = \frac{{(k_{xs,u} + k_{xs,d})}^{2} h^{4}}{16 D^{2}} - \frac{4 k_{s}}{D}$ 。

(26)

当Δ<0时,桩间土范围内的格室体的挠度求解为

$w_{s} (ρ) = C_{3} α_{0} (ρ) + C_{4} β_{0} (ρ) + C_{5} χ_{0} (ρ) + C_{6} κ_{0} (ρ) + \frac{q_{s}}{k_{s}}$ ,

(27)

式中,α₀,β₀,χ₀,κ₀为线性独立函数,C₃,C₄,C₅,C₆为待定系数。

当Δ>0时,桩间土范围内的格室体的挠度求解为

$\begin{array}{l} w_{s} (ρ) = C_{3} I_{0} (\sqrt{- λ_{21}} ρ) + C_{4} K_{0} (\sqrt{- λ_{21}} ρ) + \\ C_{5} I_{0} (\sqrt{- λ_{22}} ρ) + C_{6} K_{0} (\sqrt{- λ_{22}} ρ) + \frac{q_{s}}{k_{s}} \end{array}$ ,

(28)

${\begin{cases} λ_{21} \\ λ_{22} \end{cases} = - \frac{1}{2} [\frac{(k_{xs,u} + k_{xs,d}) h^{2}}{4 D} \pm \sqrt{Δ}]$ ,

(29)

式中,K_N为第二类N阶虚变量Bessel函数,在此解中N=0。

由上述各式可进一步得出桩顶及桩间土范围内格室体的转角、弯矩及剪力的表达式。

薄板在ρ=d_p/2处的连续条件为

${\begin{cases} {w_{p} |}_{ρ = \frac{d_{p}}{2}} = {w_{s} |}_{ρ = \frac{d_{p}}{2}}, \\ {θ_{ρ,p} |}_{ρ = \frac{d_{p}}{2}} = {θ_{ρ,s} |}_{ρ = \frac{d_{p}}{2}}, \\ {M_{ρ,p} |}_{ρ = \frac{d_{p}}{2}} = {M_{ρ,s} |}_{ρ = \frac{d_{p}}{2}}, \\ {F_{Sρ,p} |}_{ρ = \frac{d_{p}}{2}} {- F_{Sf,p} |}_{ρ = \frac{d_{p}}{2}} = {F_{Sρ,s} |}_{ρ = \frac{d_{p}}{2}} {- F_{Sf,s} |}_{ρ = \frac{d_{p}}{2}} 。 \end{cases}$

(30)

桩顶范围内,θ_ρ,p,M_ρ,p及F_Sρ,p分别表示土工格室的径向转角、径向弯矩及径向剪力；桩间土范围内则用θ_ρ,s,M_ρ,s及F_Sρ,s分别表示。F_Sf,s为桩间土范围内界面摩阻引起的格室体径向截面剪力。

在ρ=d_e/2处有边界条件：

${\begin{cases} {θ_{ρ,s} |}_{ρ = \frac{d_{e}}{2}} = 0, \\ {F_{Sρ,s} |}_{ρ = \frac{d_{e}}{2}} = 0 。 \end{cases}$

(31)

假设q_p已知,还有

${w_{p} |}_{ρ = 0} = \frac{p_{ep}}{k_{p}}$ ,

(32)

式中,k_p为桩体变形刚度系数。

联立式（30）～（32）,即可求出参数C₁～C₆及p_ep。实际上,桩土刚度差异将导致在路堤荷载作用下桩土上部的路堤土发生差异沉降,因而产生土拱效应,土拱效应会进一步改变桩与桩间土上的荷载,桩顶会形成应力集中,即p_ep>p_es,前文所导得的土拱效应公式（17）可表示p_ep与Δs的关系。

文献[26]指出,格室体界面与填料之间的摩阻与剪切位移会在受到上部力的作用时显示为近似双曲线关系,通常将其简化为理想弹塑性模型,据此文献[27]选用如图3所示进行描述。

${\begin{cases} τ = \frac{σ_{n} \tan φ_{g}}{u_{u}} u (u < u_{u}), \\ τ = τ_{\max} = σ_{n} \tan φ_{g} (u \geq u_{u}) 。 \end{cases}$

(33)

式中 τ为摩阻力； $σ_{n}$ 摩擦面正应力；ϕ_g界面摩擦角；u_u为极限相对位移；τ_max为界面极限摩阻力。

又有摩阻力与位移间的近似线性关系,得到

${\begin{cases} k_{xp,u} = \frac{\tan φ_{g}}{u_{u}} q_{p}, k_{xp,d} = \frac{\tan φ_{g}}{u_{u}} p_{ep}, \\ k_{xs,u} = \frac{\tan φ_{g}}{u_{u}} q_{s}, k_{xs,d} = \frac{\tan φ_{g}}{u_{u}} p_{es} 。 \end{cases}$

(34)

根据上述分析,将式（30）～（33）,（35）联立,从而构成非线性代数方程组,用来分析土拱效应影响下,土工格室加筋垫层挠曲变形情况。由于直接求解较为困难,为简化求解过程,本文采用迭代求解方法,步骤如下：

（1）首先初始计算时暂且忽略土拱效应影响,认为q_p,1与q_s,1相等并均等于q,根据式（34）进行计算,获得k_xp,d,k_xp,u,k_xs,u,k_xs,d的具体结果,同时由于p_ep与p_es未知,故也近似为q进行计算并代入式（34）计算出k_xp,d与k_xs,d的值,联立方程组（30）～（32）得Δs₁,p_p,1,p_s,1等参数。

（2）将步骤（1）中求得的Δs₁代入式（17）中进行计算,反算出q_p,2和q_s,2,再利用步骤（1）中的方式将这两个值代入（34）并联立（30）～（32）求得各个参数。

（3）当到达某次运算如第j次运算时,类似步骤（1）、（2）的方式,q_p,j与q_s,j利用Δs_j-1得到,此时式（34）具体表现为式（35）,以此计算界面摩阻系数,然后联立方程组（30）～（32）得到格室体挠度曲线及Δs_j等相关参数。

${\begin{cases} k_{xp,u} = \frac{\tan φ_{g}}{u_{u}} q_{p, j}, k_{xp,d} = \frac{\tan φ_{g}}{u_{u}} p_{ep, j - 1}, \\ k_{xs,u} = \frac{\tan φ_{g}}{u_{u}} q_{s, j}, k_{xs,d} = \frac{\tan φ_{g}}{u_{u}} p_{es, j - 1} 。 \end{cases}$

(35)

（4）按照步骤（1）,（2）,（3）的方式对每次数据进行反复迭代计算,直到Δs_j与Δs_j-1差值满足初始设置的误差要求时结束运算,输出所需结果。

3. 算例验证

本文选取京珠高速公路临长段k111+620～750段中k111+720中心格室体进行变形计算,此段公路为填方区,且下卧较深厚软土层。为处理大面积软土区域,使其承载性能达到预期效果,对软土深厚>5 m的部分路段,采用“沉管碎石桩+土工格室”双向增强复合地基进行处治。本算例选取路段采用梅花形布桩方式,碎石采用未风化干净砾石,砾石粒径20～40 mm,含泥量<5%,自然级配。图4为路基典型设计断面,表1为软基主要物理力学指标,表2为其他路段相关数据。

图 4 路基典型横断面设计图

Figure 4. Typical cross-sectional design of roadbed

下载: 全尺寸图片幻灯片

表 1 软土主要物理力学性质指标

Table 1. Main physical and mechanical properties of soft soil

含水率/%	孔隙比e	液限/%	塑性指数I_P	黏聚力c/kPa	内摩擦角φ/(°)	压缩系数/MPa^-1	密度/(g·cm^-3)
33~50	0.8~1.2	40~45	10~25	3~25	10~28	0.3~1.0	1.7~2.0

下载: 导出CSV

| 显示表格

表 2 路段相关数据

Table 2. Relevant statistics of road

天然地基承载力/kPa	碎石桩设计直径d_p/cm	桩间距s_a/m	桩长L_p/m	格栅屈服强度σ/MPa	格栅厚度h_g/cm	垫层厚度h/cm	路堤土填土重度γ/(kN·m^-3)	内摩擦角/(°)	路堤高度H/m	总沉降S/cm	桩土应力比n
55	38.5	1.5	10	22.5	10	50	20	30	6	32.5	4~6

下载: 导出CSV

| 显示表格

结合现有研究,桩体变形刚度和基床系数通常由静载荷试验、理论计算或参考规范获取,而界面摩阻系数、格室垫层变形刚度及弹性模量等可参考前人研究^[27-28]得出,如表3所示。

表 3 相关计算参数取值

Table 3. Values of relevant parameters

格室厚度h/m	格室复合弹性模量E/MPa	格室复合泊松比ν	界面摩擦角ϕ_g	界面屈服位移u_u/mm	桩土变形刚度k_p/(kN·m^-1)	基床系数k_s/(kPa·m^-1)
0.1	55	0.35	40°	2	200	300

下载: 导出CSV

| 显示表格

根据上述参数,采用本文方法,计算所选取路段下格室体变形,进一步得到了沉降与桩土应力比,如表4所示。由表可知,格室上下界面摩阻力及土拱效应可明显增大桩土应力比,而采用本文方法所得的沉降计算值与实测值均较为接近,证明本文方法具有合理性。通过对比表中网上、网下桩土应力比的计算结果可知土工格室加筋体可明显起到调节荷载分配的积极作用。

表 4 计算结果与实测结果对比

Table 4. Comparison between calculated and measured results

	方法对比	网上桩土应力比n_e	网下桩土应力比n	沉降S/cm
	现场实测	—	4~6	32.50
	忽略k_x与土拱	—	5.00	33.51
本文方法	考虑k_x忽略土拱	—	5.49	31.99
	考虑k_x与土拱	4.71	5.81	31.20

下载: 导出CSV

| 显示表格

4. 桩土刚度比k_p/k_s、路堤填土压缩模量E_e及路堤填土内摩擦角ϕ_e对桩土应力比n_e与n的影响

限于文章篇幅,本文将桩体刚度系数K_p取为固定值500 kN/m,通过分析部分重要计算参数对网上网下桩土应力比n_e与n和对格室体变形的影响情况进而得出土拱、垫层等对路堤荷载在桩与土部分的分配情况的影响。

由图5可知,k_p/k_s的值在很大程度上决定n_e与n的计算结果,格室垫层对荷载的调节作用在桩土刚度比较小时不明显,而随着k_p/k_s的值的增大几乎呈线性增加,且这个作用效果与垫层模量正相关；图6表示E_e影响荷载分配的情况,n_e与n随着E_e的增加呈现近似的非线性增加,且速率变缓,此外,格室垫层模量E对n与n_e的影响分别在E_e较小与较大时表现明显；图7描述的是ϕ_e影响n_e与n的情况,其中,n_e随着ϕ_e的增大而增长且趋势基本不变,n随着ϕ_e的增大而增长且趋势渐缓,此外,在ϕ_e较小时,格室垫层模量对n的影响较大。

图 5 不同E条件下k_p/k_s对n_e与n的影响

Figure 5. Effects of k_p/k_s on n_e and n under different values of E

下载: 全尺寸图片幻灯片

图 6 不同E条件下E_e对n_e与n的影响

Figure 6. Effects of E_e on n_e and n under different values of E

下载: 全尺寸图片幻灯片

图 7 不同E条件下φ_e对n_e与n的影响

Figure 7. Effects of φ_eon n_e and n under different values of E

下载: 全尺寸图片幻灯片

将上述分析结果与未考虑土工格室垫层作用时所计算的桩土应力比值（约为12左右）相比较,可以明显看出,土工格室垫层在降低土拱效应及调节荷载向桩顶集中方面有显著效果,且此效果随着填土性质的下降而愈发明显。

5. 结论

（1）本文在假定内、外土柱界面侧土压力系数与摩阻力发挥程度均与界面相对位移相关的基础上,对传统土柱模型进行了改进,得出本文土拱效应分析模型,进而运用数学方法建立路堤底部桩土应力比与桩土差异沉降的函数关系,这一改进土柱模型能够在土柱间摩阻力与差异变形之间建立联系,并能够得出相较其他方法而言更为简单清晰的s–n_e关系式。

（2）土工格室垫层在实际工程中为三维立体结构,会受到水平方向的界面摩阻力及竖直方向的路堤荷载与地基反力,本文基于此受力特点视格室垫层为薄板模型并建立挠曲控制方程获得挠曲函数,结合土拱效应的影响,建立了考虑变形协调的路堤–土工格室–桩土加固区共同作用模型,并给出了迭代求解步骤。

（3）由参数分析可知,格室垫层发挥作用会降低路堤土拱效应,并调节荷载向桩顶分配,这一效果随着土性质的变差而愈发显著。格室体上下界面的摩阻效应有利于降低其变形,具体而言,随着摩阻系数增大到一定值后,界面摩阻力会对结构产生附加的弯矩和剪力,这些作用将主导荷载重新分配。

（4）通常格室垫层上下表面与路堤填土的摩阻力存在差异,这将使格室体中面产生位移,为简化分析本文未考虑这一情况,也忽略了其由于挠曲所造成拉伸,这样的简化方法虽然不会明显影响最终桩土应力比及沉降的结果,但在加筋体内力分析,尤其是在求解最大拉应力时,难免产生一定误差,因此尚需深入研究寻求更为合理的方法。

图 1 大坝沉降与水位变化关系图

Figure 1. Relationship between dam subsidence and change of water level

下载: 全尺寸图片幻灯片

图 2 单颗粒压缩试验

Figure 2. Single particle compression test

下载: 全尺寸图片幻灯片

图 3 干燥与饱和颗粒的典型荷载–位移曲线

Figure 3. Typical load-displacement curves for dry and saturated particles

下载: 全尺寸图片幻灯片

图 4 干燥与饱和试样的SEM电镜扫描图

Figure 4. SEM images of dried and saturated samples

下载: 全尺寸图片幻灯片

图 5 界面单元本构模型示意图

Figure 5. Constitutive relations for interface elements

下载: 全尺寸图片幻灯片

图 6 颗粒形状分布图

Figure 6. Distribution of particle shape

下载: 全尺寸图片幻灯片

图 7 部分颗粒和单颗粒压缩试验加载示意图

Figure 7. Partial particle samples and diagram of loading for single particle crushing tests

下载: 全尺寸图片幻灯片

图 8 岩石单轴压缩室内及数值试验应力–应变曲线

Figure 8. Stress-strain curves of rock in uniaxial compression physical and numerical tests

下载: 全尺寸图片幻灯片

图 9 典型荷载–位移曲线和颗粒渐进破坏过程

Figure 9. Typical load-displacement curve and progressive failure process of particles

下载: 全尺寸图片幻灯片

图 10 数值与物理试验结果及Weibull拟合曲线

Figure 10. Numerical and physical test results and Weibull fitting curves

下载: 全尺寸图片幻灯片

图 11 颗粒集合体数值试样和单轴压缩试验加载示意图

Figure 11. Schematic diagram of numerical samples and uniaxial compression test loading of particle assembly

下载: 全尺寸图片幻灯片

图 12 湿化变形特征曲线

Figure 12. Characteristic curves of wetting deformation

下载: 全尺寸图片幻灯片

图 13 单轴压缩过程中的新增颗粒数量

Figure 13. Number of new particles during uniaxial compression

下载: 全尺寸图片幻灯片

图 14 试样级配曲线( $σ_{v}$ =11.25 MPa)

Figure 14. Gradation curves of samples ( $σ_{v}$ =11.25 MPa)

下载: 全尺寸图片幻灯片

表 1 颗粒形状参数表

Table 1 Parameters of particle shape

形状参数	公式	平均值	标准差
伸长率	EI=I/L	0.80	0.115
扁平率	FI=S/I	0.79	0.099
Domokos因子	$\begin{array}{l} S_{f} = (1 / S + 1 / I + 1 / L) \times \\ \sqrt{S^{2} + I^{2} + L^{2}} / \sqrt{3} \end{array}$	3.21	0.088
凸度	$C_{X} = \frac{V}{V_{ch}}$	0.88	0.031
球度	$ψ_{3D} = \frac{\sqrt[3]{36 π V^{2}}}{A}$	0.81	0.031
注：S,I,L分别为颗粒的短轴、中轴、长轴,V和A为颗粒的实际体积和表面积, $V_{ch}$ 为颗粒外围凸包体的体积。

下载: 导出CSV

表 2 细观参数取值

Table 2 Values of mesoscopic parameters

细观参数		干燥试样	饱和试样
实体单元	密度 $ρ$ /(kg·m^-3)	2240	2400
	弹性模量 $E$ /GPa	13.2	11.88
	泊松比 $v$	0.2	0.2
界面单元	抗压强度 $f_{c}^{e}$ /MPa	110	96
	压拉强度比 $f_{c} / f_{t}$	18.96	21.30
	Weibull模数 $m$	10	10
	内摩擦角 $φ$ /(°)	40	37
	Ⅰ型断裂能 $G_{Ι C}$ /(N·m^-1)	85.5	72.0
	Ⅱ型断裂能 $G_{Ⅱ C}$ /(N·m^-1)	427.5	360.0
接触模型	摩擦系数 $μ_{f}$	0.57700	0.50925
	法向刚度 $k_{n}$ /(N·m^-3)	9.9×10¹³	8.91×10¹³
	切向刚度 $k_{s}$ /(N·m^-3)	4.125×10¹³	3.7125×10¹³

下载: 导出CSV

参考文献(33)

[1]	NAYLOR D J, MARANHA J R, NEVES E M D, et al. A back-analysis of Beliche Dam[J]. Géotechnique, 1997, 47(2): 221-233. doi: 10.1680/geot.1997.47.2.221
[2]	胡超, 周伟, 常晓林, 等. 基于内聚力模型的高心墙堆石坝坝顶裂缝模拟及其成因分析[J]. 中南大学学报(自然科学版), 2014, 45(7): 2303-2310. https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD201407023.htm HU Chao, ZHOU Wei, CHANG Xiao-lin, et al. Modeling of dam crest cracks of high core rockfill dam based on cohesive zone model and crack formation analysis[J]. Journal of Central South University (Science and Technology), 2014, 45(7): 2303-2310. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD201407023.htm
[3]	张丙印, 师瑞锋. 流变变形对高面板堆石坝面板脱空的影响分析[J]. 岩土力学, 2004, 25(8): 1179-1184. doi: 10.3969/j.issn.1000-7598.2004.08.001 ZHANG Bing-yin, SHI Rui-feng. Influence of creeping on separation between concrete slab and cushion layer in high concrete face rock-fill dam[J]. Rock and Soil Mechanics, 2004, 25(8): 1179-1184. (in Chinese) doi: 10.3969/j.issn.1000-7598.2004.08.001
[4]	BAUD P, ZHU W, WONG T F. Failure Mode and Weakening Effect of Water on Sandstone[M]. New York: John Wiley & Sons, Ltd, 2000.
[5]	HAM T G, NAKATA Y, ORENSE R, et al. Influence of water on the compression behavior of decomposed granite soil[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2010, 136(5): 697-705. doi: 10.1061/(ASCE)GT.1943-5606.0000274
[6]	魏松, 朱俊高. 粗粒料三轴湿化颗粒破碎试验研究[J]. 岩石力学与工程学报, 2006, 25(6): 1252-1258. doi: 10.3321/j.issn:1000-6915.2006.06.025 WEI Song, ZHU Jun-gao. Study on wetting breakage of coarse-grained materials in triaxial test[J]. Chinese Journal of Rock Mechanics and Engineering, 2006, 25(6): 1252-1258. (in Chinese) doi: 10.3321/j.issn:1000-6915.2006.06.025
[7]	丁艳辉, 袁会娜, 张丙印. 堆石料非饱和湿化变形特性试验研究[J]. 工程力学, 2013, 30(9): 139-143. https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX201309022.htm DING Yan-hui, YUAN Hui-na, ZHANG Bing-yin. Unsaturated wetting deformation characteristics of rockfill materials[J]. Engineering Mechanics, 2013, 30(9): 139-143. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX201309022.htm
[8]	丁艳辉, 张丙印, 钱晓翔, 等. 堆石料湿化变形特性试验研究[J]. 岩土力学, 2019, 40(8): 2975-2981. doi: 10.16285/j.rsm.2018.0891 DING Yan-hui, ZHANG Bing-yin, QIAN Xiao-xiang, et al. Experimental study of the characteristics of wetting deformation of rockfill materials[J]. Rock and Soil Mechanics, 2019, 40(8): 2975-2981. (in Chinese) doi: 10.16285/j.rsm.2018.0891
[9]	贾宇峰, 姚世恩, 迟世春. 等应力比路径下粗粒土湿化试验研究[J]. 岩土工程学报, 2019, 41(4): 648-654. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201904009.htm JIA Yu-feng, YAO Shi-en, CHI Shi-chun. Wetting of coarse-grained soil under equal stress ratio path[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(4): 648-654. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201904009.htm
[10]	JIA Y F, XU B, DESAI C S, et al. Rockfill particle breakage generated by wetting deformation under the complex stress path[J]. International Journal of Geomechanics, 2020, 20(10): 1-11.
[11]	ULUSAY R, KARAKUL H. Assessment of basic friction angles of various rock types from Turkey under dry, wet and submerged conditions and some considerations on tilt testing[J]. Bulletin of Engineering Geology and the Environment, 2016, 75(4): 1683-1699. doi: 10.1007/s10064-015-0828-4
[12]	沈广军, 殷宗泽. 粗粒料浸水变形分析方法的改进[J]. 岩石力学与工程学报, 2009, 28(12): 2437-2444. doi: 10.3321/j.issn:1000-6915.2009.12.008 SHEN Guang-jun, YIN Zong-ze. Improvement of wetting deformation analysis method of coarse-grained materials[J]. Chinese Journal of Rock Mechanics and Engineering, 2009, 28(12): 2437-2444. (in Chinese) doi: 10.3321/j.issn:1000-6915.2009.12.008
[13]	程展林, 左永振, 丁红顺, 等. 堆石料湿化特性试验研究[J]. 岩土工程学报, 2010, 32(2): 243-247. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201002016.htm CENG Zhan-lin, ZUO Yong-zhen, Ding Hong-shun, et al. Wetting characteristics of coarse-grained materials[J]. Chinese Journal of Geotechnical Engineering, 2010, 32(2): 243-247. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201002016.htm
[14]	曹光栩, 宋二祥, 徐明. 碎石料干湿循环变形试验及计算方法[J]. 哈尔滨工业大学学报, 2011, 43(10): 98-104. https://www.cnki.com.cn/Article/CJFDTOTAL-HEBX201110022.htm CAO Guang-xu, SONG Er-xiang, XU Ming. Study on experiment and calculation method of dry-wet cycle characteristics of rockfills[J]. Journal of Harbin Institute of Technology, 2011, 43(10): 98-104. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-HEBX201110022.htm
[15]	迟世春, 周雄雄. 堆石料的湿化变形模型[J]. 岩土工程学报, 2017, 39(1): 48-55. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201701003.htm CHI Shi-chun, ZHOU Xiong-xiong. Slaking deformation model for rockfill materials[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(1): 48-55. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201701003.htm
[16]	ZHOU X X, CHI S C, JIA Y F, et al. A new wetting deformation simulation method based on changes in mechanical properties[J]. Computers and Geotechnics, 2020, 117: 1-21.
[17]	黄绪武, 周伟, 马刚, 等. 考虑摩擦系数和颗粒强度劣化效应的堆石体湿化细观数值模拟[J]. 中国农村水利水电, 2017(9): 125-131. https://www.cnki.com.cn/Article/CJFDTOTAL-ZNSD201709029.htm HUANG Xu-wu, ZHOU Wei, MA Gang, et al. Numerical simulation of rockfill wetting considering deterioration of friction coefficient and bond strength[J]. China Rural Water and Hydropower, 2017(9): 125-131. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZNSD201709029.htm
[18]	杨贵, 刘汉龙, 朱俊高. 粗粒料湿化变形数值模拟研究[J]. 防灾减灾工程学报, 2012, 32(5): 535-538. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXK201205004.htm YANG Gui, LIU Han-long, ZHU Jun-gao. Numerical simulation research on coarse material wetting deformation[J]. Journal of Disaster Prevention and Mitigation Engineering, 2012, 32(5): 535-538. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DZXK201205004.htm
[19]	王蕴嘉, 周梦佳, 宋二祥. 考虑颗粒破碎的堆石料湿化变形特性离散元模拟研究[J]. 工程力学, 2018, 35(增刊1): 217-222. https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX2018S1038.htm WANG Yun-jia, ZHOU Meng-jia, SONG Er-xiang. DEM simulation of wetting deformation characteristics of rockfill considering particle breakage[J]. Engineering Mechanics, 2018, 35(S1): 217-222. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX2018S1038.htm
[20]	SHAO X Q, CHI S C, TAO Y, et al. DEM simulation of the size effect on the wetting deformation of rockfill materials based on single-particle crushing tests[J]. Computers and Geotechnics, 2020, 123(3): 1-12.
[21]	FENG R, HE Y L, CAO X X, et al. Different deformation patterns in high core wall rockfill dams: a case study of the Maoergai and Qiaoqi Dams[J]. Advances in Civil Engineering, 2019(2): 1-17.
[22]	MA G, ZHOU W, CHANG X L, et al. A hybrid approach for modeling of breakable granular materials using combined finite-discrete element method[J]. Granular Matter, 2016, 18(1): 1-17.
[23]	邓璇璇, 马刚, 周伟, 等. 局部约束模式对单颗粒破碎强度的影响[J]. 浙江大学学报(工学版), 2018, 52(7): 1329-1337. https://www.cnki.com.cn/Article/CJFDTOTAL-ZDZC201807013.htm DENG Xuan-xuan, MA Gang, ZHOU Wei, et al. Effect of local constraints patterns on fragmentation of single grain[J]. Journal of Zhejiang University (Engineering Science), 2018, 52(7): 1329-1337. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZDZC201807013.htm
[24]	陈兴, 马刚, 周伟, 等. 无序性对脆性材料冲击破碎的影响[J]. 物理学报, 2018, 67(14): 219-228. https://www.cnki.com.cn/Article/CJFDTOTAL-WLXB201814024.htm CHEN Xing, MA Gang, ZHOU Wei, et al. The influence of disorder on the impact crushing of brittle materials[J]. Acta Physica Sinica, 2018, 67(14): 219-228. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-WLXB201814024.htm
[25]	MCDOWELL G R. Statistics of soil particle strength[J]. Géotechnique, 2001, 51(10): 897-900.
[26]	WEI D, ZHAO B, DIAS-DA-COSTA D, et al. An FDEM study of particle breakage under rotational point loading[J]. Engineering Fracture Mechanics, 2019, 212: 221-237.
[27]	OVALLE C, FROSSARD E, DANO C, et al. The effect of size on the strength of coarse rock aggregates and large rockfill samples through experimental data[J]. Acta Mechanica, 2014, 225(8): 2199-2216.
[28]	周海娟, 马刚, 袁葳, 等. 堆石颗粒压缩破碎强度的尺寸效应[J]. 岩土力学, 2017, 38(8): 2425-2433. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201708035.htm ZHOU Hai-juan, MA Gang, YUAN Wei, et al. Size effect on the crushing strengths of rock particles[J]. Rock and Soil Mechanics, 2017, 38(8): 2425-2433. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201708035.htm
[29]	MCDOWELL G R, AMON A. The application of Weibull statistics to the fracture of soil particles[J]. Soils and Foundations, 2000, 40(5): 133-141.
[30]	孙壮壮, 马刚, 周伟, 等. 颗粒形状对堆石颗粒破碎强度尺寸效应的影响[J]. 岩土力学, 2021, 42(2): 1-10. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202102015.htm SUN Zhuang-zhuang, MA Gang, ZHOU Wei, et al. Effect of particle shape on size effect of crushing strength of rockfill particles[J]. Rock and Soil Mechanics, 2021, 42(2): 1-10. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202102015.htm
[31]	MA G, ZHOU W, REGUEIRO R A, et al. Modeling the fragmentation of rock grains using computed tomography and combined FDEM[J]. Powder Technology, 2017, 308: 388-397.
[32]	MCDOWELL G R, DE BONO J P. On the micro mechanics of one-dimensional normal compression[J]. Géotechnique, 2013, 63(11): 895-908.
[33]	MA G, CHEN Y, YAO F H, et al. Evolution of particle size and shape towards a steady state: insights from FDEM simulations of crushable granular materials[J]. Computers and Geotechnics, 2019, 112: 147-158.