力学与实践, 2022, 44(1): 12-21 DOI: 10.6052/1000-0879-21-428

应用研究

含水率梯度作用下填埋场压实黏土层开裂特性试验与机理分析1)

刘瑞琪*,, 雷学文*, 万勇,,2), 刘磊

*武汉科技大学 城市建设学院,武汉 430065

中国科学院武汉岩土力学研究所 岩土力学与工程国家重点实验室,武汉 430071

EXPERIMENTS ON CRACKING CHARACTERISTICS AND MECHANISM ANALYSIS OF COMPACTED CLAY LAYERS IN LANDFILLS UNDER THE ACTION OF WATER CONTENT GRADIENT1)

LIU Ruiqi*,, LEI Xuewen*, WAN Yong,,2), LIU Lei

*School of Urban Construction, Wuhan University of Science and Technology, Wuhan 430065, China

State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan 430071, China

通讯作者: 2)万勇,副研究员。主要从事固体废弃物安全处置与填埋场灾变机理方面的研究。E-mail:ywan@whrsm.ac.cn

责任编辑: 胡漫

收稿日期: 2021-10-12   修回日期: 2021-11-17  

基金资助: 1)国家自然科学基金资助项目(51979267)
国家自然科学基金资助项目(41772342)

Received: 2021-10-12   Revised: 2021-11-17  

作者简介 About authors

摘要

为探究压实黏土防渗层干缩开裂的力学机制,开展了压实黏土层含水率梯度变化与干缩裂隙演化规律试验研究,并构建了含水率梯度作用下压实黏土层开裂分析模型,研究结果表明:蒸发过程中,压实黏土层失水率随深度指数递减,高液限黏土含水率梯度明显高于低液限黏土。压实黏土层表层拉应力随含水率梯度的增加而增加。蒸发过程中,低液限压实黏土层含水率梯度小,表层拉应力低,裂隙较少;而高液限黏土试验和理论分析结果与之相反。

关键词: 压实黏土; 含水率梯度; 干缩开裂; 应力分布函数; 解析解

Abstract

In order to investigate the mechanical mechanism of dry shrinkage cracking in compacted clay impermeable layer, the experimental study of water content gradient change and dry shrinkage crack evolution law of compacted clay layer was carried out. Meanwhile, the analytical model of compacted clay layer cracking under the action of water content gradient was constructed. The research results showed that the water loss rate of compacted clay layer decreases with depth index during evaporation, and the water content gradient of high liquid limit clay is significantly higher than that of low liquid limit clay. The tensile stress in the surface layer of compacted clay layer increases with the increase of water content gradient. During evaporation, the water content gradient of low liquid-limit compacted clay layer is small, the surface tensile stress is low, and there are fewer fissures; while the results of experimental and theoretical analysis of high-liquid-limit clay are opposite.

Keywords: compacted clay liner; water content gradient; dry shrinkage cracking; stress distribution function; analytical solution

PDF (5441KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

刘瑞琪, 雷学文, 万勇, 刘磊. 含水率梯度作用下填埋场压实黏土层开裂特性试验与机理分析1). 力学与实践, 2022, 44(1): 12-21 DOI:10.6052/1000-0879-21-428

LIU Ruiqi, LEI Xuewen, WAN Yong, LIU Lei. EXPERIMENTS ON CRACKING CHARACTERISTICS AND MECHANISM ANALYSIS OF COMPACTED CLAY LAYERS IN LANDFILLS UNDER THE ACTION OF WATER CONTENT GRADIENT1). Mechanics in Engineering, 2022, 44(1): 12-21 DOI:10.6052/1000-0879-21-428

压实黏土层因其良好的低渗透特性被广泛应用于垃圾填埋场顶部防渗结构中。我国《垃圾卫生填埋场岩土工程技术规范》要求压实黏土渗透系数小于10$^{-7}$ cm/s,但因其置于填埋场顶部,容易受到大气干湿循环影响而遭到破坏。研究表明,干湿循环作用下,压实黏土层中发生大量干缩裂隙,渗透系数增加2$\sim$3个数量级。填埋场顶部压实黏土防渗结构一旦破坏,将导致雨水下渗引发内部水位雍高,安全性下降;同时填埋气体泄露,导致周边大气环境污染事故[1-2]

已有研究认为:非饱和土内部结构复杂,在外界条件发生改变时其土骨架结构会发生变形[3]。压实黏土蒸发脱水过程中,基质吸力持续增加,当基质吸力大于土壤的抗拉强度时,干缩裂隙出现[4]。根据非饱和土力学理论,在特定条件下,土壤基质吸力与含水率存在确定的函数关系,即土水特征曲线(soil water characteristic curve, SWCC)[5],因此,目前针对含水率循环变化对压实黏土宏微观结构与岩土工程参数影响的研究包括:干湿循环作用下压实黏土层宏观裂隙变化规律[6-11]以及渗透系数的变化规律[12-13]。除了黏土脱水引发的拉伸应力大于黏土抗拉强度外,土壤收缩变形受到限制也是黏土干缩开裂的必要条件之一。收缩变形的限制来源于多方面,何俊等[14]通过胶水粘贴黏土试样和环刀侧壁的方法研究了压实黏土试样边界收缩限制对黏土干缩开裂的影响。曾浩等[15]和林朱元等[16]通过试样与数值模拟方法研究了试样底部摩擦对黏土干缩开裂的影响。

笔者认为,除试样侧边和底部边界变形限制对黏土干缩开裂影响外,黏土内部因含水率非均匀分布(含水率梯度)导致的变形相互抑制也是黏土收缩变形受限的主要形式之一。对小尺寸黏土试样,黏土块体在脱水过程中尽管收缩明显但并不产生裂缝,这主要因为当块体尺寸较小时,黏土块体脱水相对均匀[12-13]。而对大尺寸黏土试样,在蒸发时段,试样表层开始失水并收缩,但由于试样底部还来不及失水而未发生与表层土壤等同的收缩变形,表层土壤的收缩变形将受到底部土壤的抑制(需满足变形协调条件),工程体土层间将产生层间错动力(剪应力),同时,在工程体中部将产生一个水平拉应力与此平衡,当水平拉应力大于抗拉强度时,裂隙产生。对黏土工程体裂隙单元受力分析可知(如图1所示),压实黏土在干湿循环过程中开裂必须达到两个条件,一是压实黏土层必须大于某一尺寸[17],同时压实黏土层间剪应力梯度(因含水率梯度产生)必须达到某一数值[18-19]

图1

图1   含水率梯度作用下压实黏土层开裂力学机制[17-19]

Fig. 1   Mechanical mechanism of cracking in compacted clay layers under the action of water content gradient[17-19]


为进一步探究压实黏土层干缩开裂的主要影响因素,本文选择了高液限和低液限两种黏土开展了压实黏土层干缩开裂特征与含水率梯度变化试验研究,获得了压实黏土试样干缩开裂特征参数与含水率梯度的变化规律,构建了含水率梯度作用下压实黏土层应力分布函数及解析解,运用该模型对含水率梯度作用下压实黏土干缩开裂力学机理及主要影响因素进行了分析,研究成果对填埋场顶部压实黏土防渗层安全分析与开裂失效机理揭示具有重要意义。

1 材料和方法

1.1 试验材料

本次试验选取的两种不同塑性指标土壤,均取自湖北地区,被广泛用作该地区填埋场顶部防渗材料,土壤物理性质如表1所示,根据土壤分类方法,Soil 1为粉质黏土,Soil 2为黏土。矿物成分如表2所示,由表2可知:两土壤所含矿物类型相同,其中Soil 2中黏土矿物含量(illite, montmirillonite and kaolinite)明显高于Soil 1。

表1   土壤物理性质

Table 1  Physical characteristics of soil

新窗口打开| 下载CSV


表2   土壤中矿物成分

Table 2  Mineral composition in soil

新窗口打开| 下载CSV


1.2 试验方法

1.2.1 压实黏土干缩裂隙与含水率梯度测试

压实黏土层表层裂隙与含水率梯度分布测试试验所用试样为圆柱形,试样高10 cm,直径28.6 cm,土壤干密度均为1.5 g/cm$^{3}$。每种土壤准备两个试样,一个用于表面裂隙观测(试样A),一个用于含水率梯度分布测试(试样B)。试验步骤如下。

(1) 将土壤风干,碾碎,过2 mm筛,洒水至塑限含水率,然后按预设干密度和试样尺寸称量好土壤质量,分5层压入圆形有机玻璃柱中(尺寸为:外径30 cm,高10 cm,壁厚0.7 cm),如图2(a)。

图2

图2   压实黏土层试样制作与测试分析

Fig. 2   Fabrication and test analysis of compacted clay layer specimens


(2) 将试样真空抽气饱和48 h,然后将A、B两试样放入烘箱中干燥(温度50${^\circ}$C,与垃圾填埋场顶部接近),试样底部由有机玻璃板密封,水分仅从顶部蒸发。在不同时刻,通过特制环刀(高10 cm,内径3.91 cm)从压实黏土试样A中取样1个,分层测试含水率(1 cm,3 cm,5 cm,7 cm和9 cm处),如图2(b)所示。

(3) 采用数码相机对压实黏土试样B进行定期拍照,并将图片二值化,如图2(c)所示。通过专业软件统计内部裂隙率($SCR$),边界收缩率($BSR$)和总裂隙率($TSR$),如式(1)$\sim$式(3)所示。

$SCR=\dfrac{A_{{\rm c}} }{A_{{\rm o}} }\times 100\%$
$BSR=\dfrac{A_{{\rm B}} }{A_{{\rm o}} }\times 100\%$
$TSR=BSR+SCR$

式中$A_{{\rm c}} $, $A_{{\rm B}} $和$A_{{\rm o}}$分别为内部裂隙面积、边界收缩面积和压实黏土层试样总面积($A_{{\rm o}}={\pi }/{4}\times 28.6^{2}$ cm$^{2})$。

1.2.2 压实黏土收缩特征曲线测试

压实黏土收缩特征曲线(soil shrinkage characteristic curve, SSCC)对评价压实黏土干缩开裂与含水率的直接关系至关重要。本文采用环刀试样(直径为6.18 cm,高2 cm)作为压实黏土的收缩特征曲线的初始试样,将黏土按预设干密度压实成型后(每种类型黏土环刀试样3个),真空饱和48 h,然后放入烘箱中(温度50${^\circ}$C),定期称重并用游标卡尺测量试样直径,用3个试样的平均值作为最终的试样结果。压实黏土自由收缩率($FRS$)通过式(4)计算

$FSR=\dfrac{A_{{\rm o}} -A_{\omega } }{A_{{\rm o}} }\times 100\%$

式中,$A_{\omega } $是含水率为$\omega$的环刀试验面积。

2 结果与分析

2.1 表层收缩变形观测结果

图3为压实黏土表层裂隙照片的二值化处理结果。图中黑色部分为裂隙区域,白色部分为土壤区域,其中黑色外圈为压实黏土试样边界收缩面积,圈内黑色部分为压实黏土干缩开裂面积。由图3可知:干燥过程中,压实黏土试样表层出现裂缝,同时试样整体也发生收缩。干燥过程中,Soil 1仅有少量不规则裂隙出现在土壤表层,但试样边界大幅度地收缩。而土壤Soil 2表面出现龟裂现象,试样边界收缩较Soil 1小。

图3

图3   不同干燥时间压实黏土表层特性特征

Fig. 3   Characteristics of compacted clay surface layer with different drying time


图4为压实黏土表层裂隙特征参数与时间的变化关系。由图4可知:在半对数坐标中,压实黏土表层裂隙率(边界收缩裂隙率$BCR$,内部裂隙率$ICR$,总裂隙率$TCR$)均与干燥时间呈分段线性关系。当边界开始收缩时,三种压实黏土边界收缩率随着时间对数线性增加;而当内部裂隙出现后,内部裂隙率随时间对数先线性增加,某时刻达到峰值后线性递减;内部裂隙率达到峰值后,总裂隙的变化不再明显,其变化多是内部裂隙向边界收缩的转化。土壤类型不同,边界收缩和内部裂隙出现的先后和时间点不同,内部裂隙达到峰值的时间也有差异。Soil 1试样一开始仅发生边界收缩,约12 h后,表面开始出现裂隙,约25 h后,内部裂隙率达到最大值,随后内部裂隙率持续降低,而边界收缩裂隙率一直保持增加,整个干燥过程中,边界收缩率均明显大于内部裂隙率。对Soil 2,裂隙率的变化规律与Soil 1完全相反。

图4

图4   压实黏土表面裂隙率与时间的变化关系

Fig. 4   Surface crack ratio of fractality of compacted clay with time


2.2 含水率分布测试结果

图5为压实黏土不同层面含水率与干燥时间的变化关系。采用土层的中心深度作为土层的标识。由图5可知:压实黏土不同层面含水率均随干燥时间的增加逐渐递减,相同时刻,压实黏土含水率随土层深度增加逐渐增加。对Soil 1,在干燥前期,土壤含水率迅速降低,不同层面土壤失水速率差别不大,约30 h后出现拐点,此后土壤含水率缓慢减少,干燥时间达到100 h时,不同含水率均降至塑限含水率以下。Soil 2最表层土壤含水率变化趋势与Soil 1相似,但下四层失水速率明显低于表层失水速率,随着干燥时间增加,下四层脱水速率有降低趋势,但无明显拐点,干燥100 h后,Soil 2下四层含水率依然高于缩限含水率。

图5

图5   压实黏土不同层面含水率与干燥时间的变化关系

Fig. 5   Variation of moisture content of different layers of compacted clay in relation to drying time


指数函数常被用于表征表层蒸发作用下岩土层不同深度失水率$\Delta \theta_{z} $与深度$z$变化关系,其关系为

$\Delta \theta_{z} =\Delta \theta_{0} {\rm e}^{-kz}$

式中,$\Delta \theta_{0}$为表土层的失水率,%;$k$为拟和参数,反映压实黏土层中含水率梯度大小。

图6为不同干燥时间压实黏土失水率与测试点深度的指数拟合结果,拟合参数如表3表4所示。由图6表3表4可知:各拟合曲线相关系数$R$均大于90%,指数函数可很好地反映表面蒸发条件下压实黏土失水率随深度的梯度变化。同时,相同蒸发时间条件下,Soil 2和Soil 1的表层失水率$Wt$(0)相差不大,但是Soil 2的拟合参数$k$明显大于Soil 1,说明蒸发过程中,Soil 2中失水率梯度明显大于Soil 1。

图6

图6   不同干燥时间压实黏土失水率与测试点深度的指数拟合结果

Fig. 6   Exponential fitting results of water loss rate of compacted clay versus depth of test points at different drying times


表3   Soil 1失水率与深度的指数拟合参数

Table 3  Exponential fitting parameters of water loss rate versus depth of soil 1

新窗口打开| 下载CSV


表4   Soil 2失水率与深度的指数拟合参数

Table 4  Exponential fitting parameters of water loss rate versus depth of soil 2

新窗口打开| 下载CSV


2.3 土壤收缩曲线测试结果

图7为压实黏土收缩率与失水率的变化关系(点为实测数据)。已有研究表明:土体失水收缩存在四个阶段,即结构性收缩阶段,正常收缩阶段,残余收缩阶段和零收缩阶段。相同失水率下Soil 2的收缩率明显大于Soil 1,最终Soil 1和Soil 2的收缩率分别为8.48%和14.67%。土体失水收缩应变$\varepsilon_{\Delta \theta } $与失水率$\Delta \theta $之间的非线性变化关系可用$m$次多项式表示

$\varepsilon_{\Delta \theta } =\sum\limits_{i=1}^m {a_{i} \left( {\Delta \theta } \right)^{i}}$

图7

图7   压实黏土收缩曲线

Fig. 7   Shrinkage curve of compacted clay


采用式(6) ($m=3$)对压实黏土收缩率与失水率进行拟合,拟合结果如图7中实线,拟合参数如表5所示。

表5   压实黏土收缩曲线拟合参数

Table 5  Fitting parameters for the shrinkage curve of compacted clay

新窗口打开| 下载CSV


3 含水率梯度作用下压实黏土层应力分布函数及解析解

3.1 压实黏土层应力分布函数

蒸发过程中,土体内一点处的应变由应力变化导致的应变和失水引起的收缩应变两部分组成。应力变化导致的应变可由广义胡克定律确定,而失水收缩应变可由土体的SSCC确定。假设土体失水收缩具有各项同性特征,且土体失水$\Delta \theta $,引发的应变为$\varepsilon_{\Delta \theta }$,则平面应力状态下,土体的总应变与应力和失水量的关系可用式(7)表示

$\left. {{\begin{array}{l} {\varepsilon_{x} =\dfrac{1}{E}\left( {\sigma_{x} -\mu \sigma_{z} } \right)+\varepsilon_{{\Delta }\theta } } \\ {\varepsilon_{z} =\dfrac{1}{E}\left( {\sigma_{z} -\mu \sigma_{x} } \right)+\varepsilon_{{\Delta }\theta } } \\ {\gamma_{xz} =\dfrac{2\left( {1+\mu } \right)}{E}\tau_{xz} } \hfill \\ \end{array} }} \right\}$

式中,$E$和$\mu $分别为土壤的弹性模量和泊松比。

引入应力函数$\varphi $,使

$\left.\begin{array}{l} {\sigma_{x} =\dfrac{\partial^{2}\varphi }{\partial z^{2}}} \\ {\sigma_{z} =\dfrac{\partial^{2}\varphi }{\partial x^{2}}} \\ {\tau_{xz} =-\dfrac{\partial^{2}\varphi }{\partial x\partial z}}\\ \end{array} \right\}$

将式(8)代入式(7),然后再将所得结果代入应力平衡方程,得到本问题应力函数$\varphi$所必须满足的变形协调方程

$\nabla^{2}\nabla^{2}\varphi =-E\nabla^{2}\varepsilon_{\Delta \theta }$

为了求解方程(9),采用无穷级数方法对本问题进行求解,令

$\varphi =\sum\limits_{n=1}^\infty {f_{n} \left( z \right)} \sin a_{n} x$

式中,$f_{n} \left( z \right)$为待定函数,$a_{n} ={n\pi }/{L}$。

先将$E\nabla^{2}\varepsilon_{\Delta \theta }$展成周期为$L$的傅里叶级数,即

$E\nabla^{2}\varepsilon_{{\Delta }\theta } =\sum\limits_{n=1}^\infty {F_{n} \left( z \right)} \sin a_{n} x$
$F_{n} \left( z \right)=-\dfrac{2}{L}\int_0^L {E\nabla^{2}\varepsilon_{{\Delta }\theta } \sin a_{n} x{\rm d}x}$

然后将式(10)和式(11)代入式(9),得到偏微分方程组

$f_{n}^{4} \left( z \right)+a_{n}^{4}f_{n} \left( z \right)-2a_{n}^{4}f_{n}^{2} \left( z \right)=F_{n} \left( z \right)$

微分方程组(13)的通解为

$f_{n} \left( z \right)=A_{n}{\rm sh}(a_{n} z)+B_{n} {\rm ch}(a_{n} z)+\\ z\left( {C_{n}{\rm sh}(a_{n} z)+D_{n} {\rm ch}(a_{n} z)} \right)+H_{n} \left( z \right)$

式中,$H_{n} (z)$ 为方程(13)的特解,可由$F_{n}(z)$确定;$A_{n} $,$B_{n} $,$C_{n} $,$D_{n}$为待定参数组,可由边界条件确定。

$\left.\begin{array}{l}\sigma_{x} =\sum\limits_{n=1}^\infty \Big[ A_{n} a_{n}^{2} {\rm ch}(a_{n} z)+B_{n} a_{n}^{2}{\rm sh}(a_{n} z)+\\ C_{n} \left( {2a_{n}{\rm sh}(a_{n} z)+a_{n}^{2} {\rm ch}(a_{n} z)} \right)+\\ D_{n} \left( {2a_{n}{\rm ch}(a_{n} z)+a_{n}^{2}{\rm sh}(a_{n} z)} \right)+ \\ {\rm d}^{2}H_{n} \left( z \right)/{\rm d}z^{2}\Big]\sin (a_{n}x)\\ \sigma_{z} =-\sum\limits_{n=1}^\infty \Big[ A_{n} {\rm ch}(a_{n} z)+B_{n}{\rm sh}(a_{n} z)+\\ C_{n} z{\rm ch}(a_{n} z)+ D_{n} z{\rm sh}(a_{n} z)+\\ H_{n} \left( z \right)\Big] a_{n}^{2} \sin (a_{n} x)\\ \tau_{xz} =\sum\limits_{n=1}^\infty \Big[ A_{n} a_{n}{\rm sh}(a_{n} z)+B_{n} a_{n} {\rm ch}(a_{n} z)+\\ C_{n} \left( {\rm ch}(a_{n} z)+a_{n} z{\rm sh}(a_{n} z) \right)+\\ D_{n} \left( {\rm sh}(a_{n} z)+a_{n} z{\rm ch}(a_{n} z) \right)+\\ {\rm d}^{2}H_{n} \left( z \right)/{\rm d}z^{2}\Big] a_{n} \sin (a_{n} x) \end{array}\right\}$

将式(14)代入式(10)再代入式(8),得到裂隙单元体应力分布函数解析解,如式(15)所示。

3.2 应力分布函数参数求解

3.2.1 特解$H_{n} (z)$求解

根据2.3节中试验结果,假设土体失水收缩应变$\varepsilon_{\Delta \theta }$与失水率$\Delta \theta$之间变化关系可用式(6)表示;不同深度处失水率$\Delta \theta$可用式(5)表示。将式(5)代入式(6),得到不同深度处的干缩应变与表层失水率变化关系

$\varepsilon_{\Delta \theta } =\sum\limits_{i=1}^m {a_{i} \left( {\Delta \theta_{0} } \right)^{i}{\rm e}^{-{i}kz}}$

将式(16)代入式(12)得

$F_{n} (z)=\sum\limits_{i=1}^m {\dfrac{2Ea_{i} \left( {\Delta \theta_{0} } \right)^{i}(ik)^{2}}{La_{n} }} \left[ {(-1)^{n}-1} \right]$

由于$H_{n} (z)$为方程组(13)的特解,因此,我们令$A_n=B_n=C_n=D_n=0$,然后将式(17)和式(14)代入式(13),得到方程组式(13)的特解为

$H_{n} (z)=\sum\limits_{i=1}^m {H_{{ni}} }{\rm e}^{-{i}kz}$

式中

$H_{{ni}} =\dfrac{2Ea_{i} \left( {\varDelta \theta_{0} }\right)^{i}({i}k)^{2}}{La_{n} \left[ {({i}k)^{2}-a_{n}^{2} }\right]^{2}}\left[ {(-1)^{n}-1} \right]$

3.2.2 参数组$A_{n},B_{n},C_{n},D_{n} $求解

忽略土体的自重和底面摩擦,则在$z=0$,$h$处(裂隙单元上表面和底面处),$\left. {\sigma_{z} } \right|_{z=0} =\left. {\sigma_{z} } \right|_{z=h}=\left. {\tau_{xz} } \right|_{z=0} =\left. {\tau_{xz} }\right|_{z=h} =0$,将此条件分别代入式(15)中第二项和第三项,得

$\beta_{11} A_{n} +\beta_{12} B_{n} +\beta_{13} C_{n} +\beta_{14} D_{n} =-\sum\limits_{i=1}^3 {H_{{ni}} }\quad$

式中,$\beta_{11} =1$,$\beta_{12} =\beta_{13} =\beta_{14} =0$。

$\beta_{21} A_{n} +\beta_{22} B_{n} +\beta_{23} C_{n} +\beta_{24} D_{n} =\sum\limits_{i=1}^3 {H_{{ni}} }$

式中,$\beta_{21} =0$,$\beta_{22} =a_{n}$,$\beta_{23} =1$,$\beta_{24} =0$。

$\beta_{31} A_{n} +\beta_{32} B_{n} +\beta_{33} C_{n} +\beta_{34} D_{n} =\\ -\sum\limits_{i=1}^3 {H_{{ni}} }{\rm e}^{-{i}kh}$

式中,$\beta_{31} ={\rm ch}(a_{n} h)$,$\beta_{32} ={\rm sh}(a_{n} h)$,$\beta_{33} =h{\rm ch}(a_{n} h)$,$\beta_{34} =h{\rm sh}(a_{n} h)$

$\beta_{41} A_{n} +\beta_{42} B_{n} +\beta_{43} C_{n} +\beta_{44} D_{n} =\\ \sum\limits_{i=1}^3 {({i}k)^{2}} H_{{ni}}{\rm e}^{-{i}kh}$

式中,$\beta_{41} =a_{n}^{2} {\rm ch}(a_{n} h)$,$\beta_{42} ={\rm a}_{n}^{2}{\rm sh}(a_{n} h)$,$\beta_{43}=2a_{n} {\rm sh}(a_{n} h)+a_{n}^{2}h{\rm ch}(a_{n} h)$,$\beta_{44} =2a_{n} {\rm ch}(a_{n} h)+{\rm a}_{n}^{2} h{\rm sh}(a_{n} h)$。

联立方程(19a)$\sim$(19d),求解得到参数$A_{n}$, $B_{n}$, $C_{n}$, $D_{n}$, 如式(20)所示

$\left.\begin{array}{l} {A_{n} =-\displaystyle\sum\limits_{i=1}^m {H_{{ni}} } } \\ {B_{n} =\dfrac{\beta_{44} C_{n1} -\beta_{34} C_{n2} }{\varDelta_{n} }} \\ {{\rm C}_{n} =\displaystyle\sum\limits_{i=1}^m {H_{{ni}} } -a_{n} \dfrac{\beta_{44} C_{n1} -\beta_{34} C_{n2} }{\varDelta_{n} }} \\ {D_{n} =\dfrac{\left( {\beta_{42} -a_{n} \beta_{43} } \right)C_{n1} -\left( {\beta_{32} -a_{n} \beta_{33} } \right)C_{n2} }{\varDelta_{n} }} \\\end{array}\right\}\quad$

式中

$\begin{align} \varDelta_{n} =\left[ {\left( {a_{n} h} \right)^{2}+1} \right]\left( {\rm sh}(a_{n} h) \right)^{2}-\left( {a_{n} h{\rm ch}(a_{n} h)} \right)^{2} \\ C_{{\rm n1}} =\sum\limits_{i=1}^m {H_{{ni}} } \left( {\beta_{33} +1-{\rm e}^{-{i}kh}} \right) \\ C_{n2} =\sum\limits_{i=1}^m {H_{{ni}} } \left[ {\beta_{43} +1+({i}k)^{2}{\rm e}^{-{i}kh}} \right] \end{align}$

在$x=\pm L$ 处,有$\displaystyle\int_0^h \sigma_{x}\bigg|_{x=\pm L} {\rm d}z=0$, $\displaystyle\int_0^h \tau_{xz}\bigg|_{x=\pm L}{\rm d}z=0$,由于$x=\pm L$时,$\cos a_{n} x=0$,因此条件自然满足。

4 含水率梯度作用下压实黏土层开裂机理分析

4.1 含水率梯度作用下压实黏土层应力分析

图8为不同含水率梯度条件下压实黏土水平应力随深度的变化关系。压实黏土层参数选取如下:$h=0.1$ m,$E=6000$ kPa,$a=0.002 8$,$\Delta\theta_{0}=-20$,$k=1$,5,10和15。由图8可知:压实黏土层水平应力随深度的变化趋势分三个区间:表层拉伸应力区间,中层压缩应力区间和底层拉伸应力区间。最表层出现最大的拉伸应力,随着深度的增加,拉伸应力逐渐减小至零;随后进入中层压缩应力区间,压缩应力先增后又减至零;此后又出现拉伸应力,且随深度增加逐渐增加。三应力区间的相互作用表现在:中层区间对表层区间失水收缩的抑制,致使表层区间出现拉伸应力,同时中层区间出现压缩应力,而底层区间对中层区间收缩的抑制致使本区间又出现拉伸应力(由应力改变而引起的应变变化如图9所示),因此表面蒸发过程中,压实黏土顶部和底部均匀出现干缩裂隙,与文献[15,16]中的现象一致。表层失水率相同条件下,含水率梯度越大($k$值越大),表层的拉伸应力越大。

图8

图8   不同含水率梯度作用下压实黏土水平应力随深度的变化关系

Fig. 8   Variation of horizontal stress with depth in compacted clay under the action of different water content gradients


图9

图9   不同应力层内压实黏土的应变随深度变化($\Delta \theta_0=-20$,$k=10$)

Fig. 9   Variation of strain with depth in compacted clay within different stress layers ($\Delta \theta_0=-20$,$k=10$)


4.2 高低液限黏土干缩开裂特性对比分析

图10图11分别为两种试验土壤的应力和应变随深度的变化关系。由图10可知,蒸发过程中,Soil 1 (低液限黏土)的表层拉伸应力较低,未能达到土壤的抗拉强度,因此蒸发过程中,Soil 1表层仅出现少量的裂隙(如图3所示)。从图11可知,蒸发过程中,Soil 1收缩变形均随深度的增加而减小,但顶部和底部收缩应变相差不大,即压实黏土试样主要发生整体的收缩变形(如图3所示,蒸发过程中Soil 1表面以边界收缩为主)。而Soil 2 (高液限黏土)表层拉伸应力大,因此蒸发过程中,Soil 2表层出现大量拉伸裂隙(如图3所示);同时,Soil 2表层收缩明显大于底部,试样整体收缩不明显(如图3所示,蒸发过程中Soil 2表面的收缩以内部裂隙为主,边界收缩为辅)。

图10

图10   不同蒸发时刻两种试验土壤的应力随深度的变化关系

Fig. 10   Variation of stress with depth for two test soils at different evaporation moments


图11

图11   不同蒸发时刻两种试验土壤的应变随深度的变化关系

Fig. 11   Variation of strain with depth for two test soils at different evaporation moments


5 主要结论

(1) 开展了表面蒸发过程中高低液限两种压实黏土试样表层裂隙演化规律试验研究,结果表明:表层蒸发过程中,压实黏土含水率从表层至底层依次持续减少;高低液压实黏土表层失水速率相差不大,但高液限黏土底层失水速率明显低于低液限黏土;含水率随深度的变化符合指数函数分布。干缩开裂和边界收缩是压实黏土试样的两种收缩变形,边界收缩率随蒸发时间持续增加,而裂隙率随蒸发时间先增后减。低液限黏土主要以边界收缩为主,而高液限黏土主要以开裂变形为主。

(2) 基于脱水过程中压实黏土含水率梯度变化特征方程,压实黏土应力-应变本构关系和土层间变形协调方程,建立含水率梯度作用下压实黏土层应力分布模型,并推导了该模型解析解。研究结果表明,表层蒸发过程中,压实黏土试样水平应力随深度的变化趋势分三个区间:表层拉伸应力区间,中层压缩应力区间和底层拉伸应力区间。三应力区间的相互作用表现在:中层区间对表层区间脱水收缩的抑制,致使表层区间出现拉伸应力,同时中层区间出现压缩应力,而底层区间对中层区间收缩的抑制致使本区间又出现拉伸应力。蒸发过程中,低液限压实黏土层脱水均匀,层间相互作用不能产生较大拉伸应力致使压实黏土产生拉伸裂缝,因此主要以蒸发过程中表面干缩裂隙较少,试样整体收缩变形为主。

参考文献

薛强, 赵颖, 刘磊 .

垃圾填埋场灾变过程的温度-渗流-应力-化学耦合效应研究

岩石力学与工程学报, 2011, 30(10): 1970-1988

[本文引用: 1]

Xue Qiang, Zhao Ying, Liu Lei, et al.

Study of thermo-hydro-mechanical-chemical coupling effect of catastrophe process of landfill

Chinese Journal of Rock Mechanics and Engineering, 2011, 30(10): 1970-1988 (in Chinese)

[本文引用: 1]

陈云敏, 施建勇, 朱伟 .

环境岩土工程研究综述

土木工程学报, 2012, 45(4): 165-180

[本文引用: 1]

Chen Yunmin, Shi Jianyong, Zhu Wei, et al.

A review of geoenvironmental engineering

China Civil Engineering Journal, 2012, 45(4): 165-180 (in Chinese)

[本文引用: 1]

刘艳, 赵成刚, 蔡国庆 .

非饱和土力学理论的研究进展

力学与实践, 2015, 37(4): 457-465

[本文引用: 1]

Liu Yan, Zhao Chenggang, Cai Guoqing, et al.

Research progress of unsaturated soil mechanics

Mechanics in Engineering, 2015, 37(4): 457-465 (in Chinese)

[本文引用: 1]

郑少河, 金剑亮, 姚海林 .

地表蒸发条件下的膨胀土初始开裂分析

岩土力学, 2006, 27(12): 2229-2233

[本文引用: 1]

Zheng Shaohe, Jin Jianliang, Yao Hailin, et al.

Analysis of initial cracking behavior of expansive soil due to evaporation

Rock and Soil Mechanics, 2006, 27(12): 2229-2233 (in Chinese)

[本文引用: 1]

Lu N, Likos WJ.

Unsaturated Soil Mechanics

New York: John Wiley & Sons, 2004

[本文引用: 1]

吕海波, 曾召田, 赵艳林 .

膨胀土强度干湿循环试验研究

岩土力学, 2009, 30(12): 3797-3802

[本文引用: 1]

Haibo, Zeng Zhaotian, Zhao Yanlin, et al.

Experimental studies of strength of expansive soil in drying and wetting cycle

Rock and Soil Mechanics, 2009, 30(12): 3797-3802 (in Chinese)

[本文引用: 1]

龚壁卫, 程展林, 胡波 .

膨胀土裂隙的工程特性研究

岩土力学, 2014, 35(7): 1825-1831

Gong Biwei, Cheng Zhanlin, Hu Bo, et al.

Research on engineering properties of fissures in expansive soil

Rock and Soil Mechanics, 2014, 35(7): 1825-1831 (in Chinese)

Tang C, Shi B, Liu C, et al.

Influencing factors of geometrical structure of surface shrinkage cracks in clayey soils

Engineering Geology, 2008, 101(3-4): 204-217

DOI      URL    

Xue Q, Zhao Y, Li ZZ, et al.

Numerical simulation on the cracking and failure law of compacted clay lining in landfill closure cover system

International Journal for Numerical and Analytical Methods in Geomechanics 2014, 38(15): 1556-1584

DOI      URL    

Li JH, Li L, Chen R, et al.

Cracking and vertical preferential flow through landfill clay liners

Engineering Geology, 2016, 206: 33-41

DOI      URL    

Wan Y, Wu CX, Xue Q, et al.

Effects of plastic contamination on water evaporation and desiccation cracking in soil

Science of the Total Environment, 2019, 654(1): 576-582

DOI      URL     [本文引用: 1]

万勇, 薛强, 赵立业 .

干湿循环对填埋场压实黏土盖层渗透系数影响研究

岩土力学, 2015, 36(3): 679-686, 693

[本文引用: 2]

Wan Yong, Xue Qiang, Zhao Liye, et al.

Effects of wetting-drying cycles on permeability of compacted clay cover at landfill site

Rock and Soil Mechanics, 2015, 36(3): 679-686, 693 (in Chinese)

[本文引用: 2]

赵立业, 薛强, 万勇 .

干湿循环作用下高低液限黏土防渗性能对比研究

岩土力学, 2016, 37(2): 446-452, 464

[本文引用: 2]

Zhao Liye, Xue Qiang, Wan Yong, et al.

A comparative study of anti-seepage performance of clays with high and low liquid limits under drying-wetting cycles

Rock and Soil Mechanics, 2016, 37(2): 446-452, 464 (in Chinese)

[本文引用: 2]

何俊, 万娟, 王宇.

压实黏土干燥裂隙及渗透性能研究

工程地质学报, 2012, 20(3): 397-402

[本文引用: 1]

He Jun, Wan Juan, Wang Yu.

Desiccation cracks on hydraulic performance of compacted clay liner via laboratory wet-dry cycling tests

Chinese Journal of Engineering Geology, 2012, 20(3): 397-402 (in Chinese)

[本文引用: 1]

曾浩, 唐朝生, 林銮 .

土体干缩裂隙发育方向及演化特征的层间摩擦效应研究

岩土工程学报, 2019, 46(6): 1172-1181

[本文引用: 2]

Zeng Hao, Tang Chaosheng, Lin Luan, et al.

Interfacial friction dependence of propagation direction and evolution characteristics of soil desiccation cracks

Chinese Journal of Geotechnical Engineering, 2019, 46(6): 1172-1181 (in Chinese)

[本文引用: 2]

林朱元, 唐朝生, 曾浩 .

土体干缩开裂过程的边界效应试验与离散元模拟

岩土工程学报, 2020, 42(2): 372-380

[本文引用: 2]

Lin Zhuyuan, Tang Chaosheng, Zeng Hao, et al.

Laboratory characterization and discrete element modeling of desiccation cracking behavior of soils under different boundary conditions

Chinese Journal of Geotechnical Engineering, 2020, 42(2): 372-380 (in Chinese)

[本文引用: 2]

孙即超, 王光谦, 孙其诚.

临界状态下非饱和土裂缝开裂间距

科学通报, 2009, 54(9): 1311-1314

[本文引用: 3]

Sun Jichao, Wang Guangqian, Sun Qicheng.

Crack spacing of unsaturated soil under critical condition

Chinese Science Bulletin, 2009, 54(9): 1311-1314 (in Chinese)

[本文引用: 3]

吴珺华, 袁俊平, 卢廷浩.

基于变湿应力概念的膨胀土初始开裂分析

岩土力学, 2011, 32(6): 1631-1636

[本文引用: 1]

Wu Junhua, Yuan Junping, Lu Tinghao.

Analysis of initial cracking behavior of expansive soil due to moisture change stress

Rock and Soil Mechanics, 2011, 32(6): 1631-1636 (in Chinese)

[本文引用: 1]

Wan Y, Xue Q, Liu L.

Relationship between the shrinkage crack characteristics of compacted clay liner and moisture ratio or moisture ratio gradient

Soil and Foundation, 2018, 58(6): 1435-1445

DOI      URL     [本文引用: 3]

/