交通荷载作用下<i>K</i><sub>0</sub>固结硬黏土动力特性研究

刘旭明; 潘坤; 杨仲轩; 蔡袁强

doi:10.11779/CJGE2020S1009

交通荷载作用下K₀固结硬黏土动力特性研究 English Version

刘旭明^1,,
潘坤²,
杨仲轩^1, ,,
蔡袁强^{1, 2}

1.
浙江大学建筑工程学院,浙江　杭州　310058
2.
浙江工业大学建筑工程学院,浙江　杭州　310014

基金项目:

国家自然科学基金委员会与英国皇家学会（NSFC-RS）牛顿高级学者基金项目 51761130078

国家杰出青年基金项目 51825803

详细信息

作者简介:
刘旭明（1995— ）,男,硕士研究生,主要从事土动力学方面的研究。Email：21712019@zju.edu.cn

通讯作者:
杨仲轩, E-mail：zxyang@zju.edu.cn

中图分类号: TU411
计量
- 文章访问数: 186
- HTML全文浏览量: 19
- PDF下载量: 92
出版历程
- 收稿日期: 2020-05-05
- 网络出版日期: 2022-12-07
- 刊出日期: 2020-10-31

Cyclic response of K₀-consolidated stiff clay under traffic loading

1.
College of Civil Engineering and Architecture, Zhejiang University, Hangzhou 310058, China
2.
College of Civil Engineering and Architecture, Zhejiang University of Technology, Hangzhou 310014, China

摘要

摘要: 利用GDS动三轴仪及空心圆柱扭剪仪对K₀ > 1的英国Gault硬黏土和Kimmeridge硬黏土进行了不排水条件下大周数（1万次）循环加载试验,研究了交通荷载作用下循环应力比及主应力轴旋转对K₀固结原状硬黏土动力特性的影响。试验结果表明：与软黏土相比,硬黏土具有较为明显的脆性特征；随着循环应力比的增大,其动力变形响应可分为“稳态”、“亚稳态”和“非稳态”等三种类型。在主应力轴旋转条件下,硬黏土轴向应变在拉伸侧不断发展,回弹模量逐渐减小,土体循环弱化程度较三轴加载条件下更为明显。
- 硬黏土 /
- 交通荷载 /
- 累积应变 /
- 变形模式 /
- 回弹模量
Abstract: To study the effects of cyclic stress ratio and principle stress axis rotation on the undrained cyclic response characteristics of K₀-consolidated stiff clay whose K₀ is larger than 1, a series of traffic loading tests with a large number of cycles are conducted on the intact Gault and Kimmeridge clays by using the dynamic triaxial test system (DYNTTS) and the hollow cylinder apparatus (HCA). The test results indicate that the stiff clay has strong brittleness compared with the soft clay. With the increasing cyclic stress ratio, the cyclic deformation responses can be categorized into stable, metastable, and unstable patterns. Under the continuous rotation of the principle stress axis, the axial strain of the stiff clay is gradually accumulated at the extension side, and its resilient modulus is rapidly degraded.
- stiff clay /
- traffic loading /
- accumulative strain /
- deformation pattern /
- resilient modulus

HTML全文

0. 引言

膨胀土具有较强的亲水性,对于处于季节冻土区的渠道工程、边坡工程等,反复冻融循环对膨胀土的劣化作用会极大影响工程稳定性^[1-2]。已有的研究表明,冻融作用对土体的性质改变相对复杂,部分特征并无统一认识^[3-4],典型的如内摩擦角,变大、变小、基本不变的结论均有。由于冻融循环作用下,不同土体响应机制的复杂性,针对膨胀土特定土体的冻融分析在工程实践中具有更好的指导意义。

目前针对膨胀土的冻融循环研究已有了较大进展,蔡正银等^[2]以高寒渠道膨胀土为研究对象,在干湿冻融耦合的试验条件下系统讨论了冻融作用下膨胀土的裂隙演化规律。Lu等^[5]、Tang等^[6]则讨论了不同状态下的膨胀土与循环冻融次数相关的劣化规律,研究表明,随着冻融循环次数的增加,膨胀土的破坏强度、弹性模量、抗剪强度等力学指标劣化加剧,且第一次冻融循环下其劣化最为明显。

目前有关冻融循环下的损伤模型,更多是基于连续介质理论展开。损伤力学中的损伤因子D部分研究需要借助其他测试手段进行确定,如卢再华等^[7]通过CT扫描技术定义损伤因子,其构建的损伤模型包含23个参数,均需测试确定。目前主流的损伤因子定义是在已有的应力应变测试成果上进行,其中假设土体微元满足weibull分布,基于应变等价原理开展损伤模型构建获得了广泛应用。杨明辉等^[8]、李向东等^[9]进行了基质吸力和外部荷载耦合状态下的非饱和土损伤模型构建。Tang等^[6]基于多种技术手段综合研究了损伤程度与冻融循环次数之间的关系,并构建了冻融循环作用下的统一损伤模型。总结目前的土的损伤模型,其损伤定义时,主要考虑了弹性模量的退化,核心均集中到应力应变曲线的模拟,但是对于土体本身的性质如压硬性等考虑不足。

传统弹性力学中的应变是微应变,但是土工试验中的土体变形量较高,因此完全按照传统应变去研究其土体性质无法充分考虑土体的内在特质。基于此,通过引入自然应变的概念、weibull统计分布规律、应变等价原理等,推导构建适用于膨胀土的损伤模型,并以损伤模型为基础进行冻融损伤分析等工作。

1. 基于自然应变的冻融损伤模型

1.1 自然应变适用性

以图1中所示的岩土体为分析模型,其满足均质各向同性的基本假定,受荷时,假定无新的损伤发生,其应力应变关系用HOOK定律进行描述,

图 1 土体分析模型

Figure 1. Soil analysis model

下载: 全尺寸图片幻灯片

$d σ = K d ε$ ,

(1)

式中, $σ$ 为岩土材料截面应力（压为正）,K为岩土体的刚度, $ε$ 为当前应力状态下的应变量。

根据弹性力学理论,岩土体的应变可为

$d ε = - d H / H$ ,

(2)

式中,H为当前应力状态下的试件高度。

式（2）计算的应变,Freed将其定义为自然应变或真实应变^[10]。具体到应用阶段通过下式计算岩土体应变量更为广泛：

$d ε = - d H / H_{0}$ ,

(3)

式中,H₀为未受荷时的试件高度。对于高刚度材料,通常变形量较原来的形状较为微小,数量级相差在2个以上,此时式（2）,（3）基本等效。式（3）由于更简洁的表达形式,在岩土工程建设中具有广泛的应用场景。但对于部分软岩、松散土体、断裂结构、多孔介质等,全部采用式（3）去分析其力学特征,则可能会阻碍对岩土本质特征的正确认识。

对于土体而言,土颗粒是高刚度材料,由于空气和水的存在,其形成的多介质土体则不完全满足式（3）的适用条件。尤其是对于亲水性矿物含量高的膨胀土而言,饱和时,会极大的削弱土体的刚度,干湿循环、冻融循环的条件下,会进一步导致其刚度的降低。因此在分析其力学特征时,式（2）具有更好的适用性。

1.2 基于应变等价原理的冻融损伤变量

岩土体的损伤产生主要有受荷和环境影响2部分。针对产生损伤岩土体,目前常用的损伤变量定义方法是基于应变等价原理实现,即可将损伤状态的材料在有效应力作用下的应变与无损状态下的应变进行等效,其无损状态的材料所受的应力称为名义应力。原理图如图1所示,表达式为

$ε = \frac{σ}{E^{'}} = \frac{σ^{'}}{E}$ ,

(4)

$σ A = σ^{'} A^{'}$ ,

(5)

式中, $σ$ 为损伤状态下的有效应力,A为损伤状态下有效承载面积,E'为损伤状态下的岩土体模量,σ'等效无损状态下的名义应力,A'等效无损状态下有效承载面积,E为等效无损状态下的岩土体模量。

损伤变量可用岩土体模量来表示,

$D = 1 - E' / E$ 。

(6)

更一般的情况下,岩土体的损伤可用岩土体变形模量的损伤来定义,目前多数文献均采用此定义^{[6, 8-9]}。

基于式（2）的自然应变概念,可进一步建立岩土体刚度和传统方法变形模量的关系,

$E = \frac{σ}{1 - \exp (- σ / K)}$ 。

(7)

将式（7）代入式（6）,可得自然应变下的岩土体损伤变量为

$D = 1 - K' / K$ ,

(8)

式中,K'为损伤状态下的土体刚度。式（6）,（8）的形式基本一致,但是代表的是两种不同应变概念下的损伤定义方式。对于高刚度材料,E=K,传统的应变定义是自然应变定义下的一种特殊情况。

对于环境损伤D_f（冻融循环、干湿循环）和受荷损伤D_s耦合的情况,总的损伤变量D可表示为

$D = 1 - (1 - D_{f}) (1 - D_{s})$ 。

(9)

1.3 受荷损伤模型构建

假设土体内部微元服从概率分布,可用Weibull分布定量描述微元强度的分布特征^[11],

$φ (F) = \frac{m}{F_{0}} {(\frac{F}{F_{0}})}^{m - 1} e^{- {(\frac{F}{F_{0}})}^{m}} = \frac{d D_{s}}{d F}$ ,

(10)

式中,φ(F)为土体微元强度概率分布函数,F₀,m为Weibull分布的参数,F为微元的强度。

对式（10）进行积分可得,受荷损伤的损伤变量D_s与土体微元强度F的关系,

$D_{s} = 1 - \exp [- {(F / F_{0})}^{m}]$ 。

(11)

土体微元的强度准则采用D-P准则进行计算^[12]。土体微元的强度F为

$F = α {I^{'}}_{1} + \sqrt{{J^{'}}_{2}}$ ,

(12)

式中, $α$ 为材料参数, ${I^{'}}_{1}$ 为有效应力状态下的应力第一不变量,J₂'为有效应力状态下的偏应力第二不变量,

$\begin{array}{l} α = \frac{2 \sin φ}{\sqrt{3} (3 - \sin φ)}, \\ {I^{'}}_{1} = {σ^{'}}_{1} + 2 {σ^{'}}_{3}, \\ {J^{'}}_{2} = \frac{1}{3} {({σ^{'}}_{1} - {σ^{'}}_{3})}^{2}, \end{array}}$

(13)

其中,φ土体微元内摩擦角, ${σ^{'}}_{1}$ , ${σ^{'}}_{3}$ 为常规三轴条件下的竖向有效应力和有效围压。

通过自然应变的概念,有效应力作用下的最大主应力方向的应变 ${ε^{'}}_{1}$ 为

${ε^{'}}_{1} = \frac{1 - \exp (- \frac{{σ^{'}}_{1} - {σ^{'}}_{3}}{K'})}{({σ^{'}}_{1} - {σ^{'}}_{3})} ({σ^{'}}_{1} - 2 μ {σ^{'}}_{3})$ ,

(14)

式中,K'为土体的刚度,μ为泊松比。常规三轴试验时的固结过程,测得应力应变曲线主要是偏应力 $(σ_{1} - σ_{3})$ 和除去初始应变之后的应变 $ε_{1}$ ,基于应变等价原理,将有效应力进行转换,可将式（14）改写为

$σ_{1} - σ_{3} = - K_{0} (1 - D_{s}) \ln (1 - ε_{1})$ ,

(15)

式中,K₀为偏应力为0时的初始刚度。设孔隙水压力为P_W,则联立式（12）～（15）,土体基于自然应变的损伤本构模型如下：

$σ_{1} - σ_{3} = - K_{0} \ln (1 - ε_{1}) exp [- {(\frac{F}{F_{0}})}^{m}]$ ,

(16)

其中,

$\begin{matrix} F = {\frac{2 \sin φ}{\sqrt{3} (3 - \sin φ)} [(σ_{1} - σ_{3}) + 3 (σ_{3} - P_{W})] + \frac{1}{\sqrt{3}} (σ_{1} - σ_{3})} \cdot \\ \frac{- K_{0} \ln (1 - ε_{1})}{σ_{1} - σ_{3}} \end{matrix}$ 。

(17)

2. 基于遗传算法的模型参数确定

从式（16）,（17）中可以发现,对于损伤参量无法直接测量的参数主要有初始刚度K₀、土体微元内摩擦角φ、Weibull分布的参数F₀和m。由于土体结构的复杂性,为避免误差累积,降低模型本身的精度,具体参数确定时采用智能算法遗传算法实现。遗传算法本质上属于生物进化算法,通过模拟生物进化中的优胜劣汰,对于具体的问题给出最优解^[13]。对于给定的应力应变曲线,限定参数的取值范围,通过数值计算软件MATLAB编程,进行参数自动搜索,全局寻优,具体原理如图2所示。

图 2 遗传算法原理示意图

Figure 2. Schematic diagram of genetic algorithm

下载: 全尺寸图片幻灯片

3. 损伤模型验证

3.1 冻融循环试验数据

为进一步研究冻融循环后,膨胀土的损伤特性,利用Tang等^[6]。的试验数据进行进一步的分析。试验土样充分饱水后进行冻融循环,-10℃保持12 h后,10℃保持12 h为1次冻融循环。试验采用常规三轴压缩试验,试验条件为固结排水,因此可不考虑超静孔隙水压力影响,试验围压为100 kPa。文献[6]采用的循环冻融次数N_FT分别为0,1,3,7,11。试验结果如图3所示。

图 3 膨胀土冻融循环后的应力应变曲线

Figure 3. Stress-strain curves of expansive soils after freeze-thaw

下载: 全尺寸图片幻灯片

3.2 损伤模型验证

通过MATLAB编程,遗传算法可实现损伤模型参数的自动获取。通过实测的应变应变数据模拟,得到到实测值和模拟值的对比如图4所示,模拟参数如表1所示。从图4中可明显的发现,基于自然应变推导的损伤模型可以很好的模拟膨胀土在不同冻融循环下的应力应变关系,模拟精度较高。

图 4 损伤模型应力应变曲线模拟值与实测值对比

Figure 4. Comparison of simulated and measured values of stress-strain curves of damage model

下载: 全尺寸图片幻灯片

表 1 损伤模型拟合参数

Table 1. Fitting parameters for damage model

N_FT	K₀/MPa	φ/(°)	F₀	m
0	15.96	35.31	154.05	0.3174
1	14.89	35.83	153.06	0.3205
3	12.96	34.18	146.31	0.3200
7	12.42	35.30	145.52	0.3202
11	11.90	34.86	143.63	0.3209

下载: 导出CSV

| 显示表格

3.3 损伤模型参数分析

从表1可以明显发现,4个参数中初始刚度K₀、参数F₀呈现明显的劣化特征,土体微元内摩擦角φ,m相对变化不大。Weibull分布中的参数F₀和m通常代表着宏观平均强度和均布系数。由于土体结构本身的劣化,对土体形成一定的损伤,导致F₀的减小。由于冻融过程中,土颗粒的团聚效应,均布系数m有略微增加的趋势,验证了论文中本构模型对于分析冻融损伤的合理性。冻融损伤产生的一个重要原因就是孔隙的增加,结构劣化,因此膨胀土的冻融损伤同样表现在初始刚度K₀的明显降低。

选取受冻融影响明显的初始刚度K₀和参数F₀作为分析对象,绘制其与冻融循环次数的关系曲线如图5所示。从图5中可以明显看出,初始刚度K₀和参数F₀的劣化规律均满足一定的负指数关系,随着冻融次数的增加,劣化速度减慢,逐渐趋于一稳定值。

图 5 损伤模型参数随冻融次数的变化曲线

Figure 5. Variation of damage model parameters with number of freeze-thaw cycles

下载: 全尺寸图片幻灯片

4. 讨论

通过建立的自然应变损伤模型,高精度的模拟了膨胀土冻融循环后的应力应变曲线,以及分析了冻融损伤对模型参数的影响。为进一步研究冻融损伤的损伤演化过程,利用式（9）绘制损伤因子D的变化曲线如图6所示。通过图6可以看出,在应力应变曲线的后期,不同冻融循环下的损伤因子基本上趋于一致。在应力应变曲线前期,随着冻融循环次数的增加,冻融损伤速度也进一步加快。但是这种加快的趋势随循环次数的增加开始减弱。

图 6 损伤因子D变化曲线

Figure 6. Curves of damage factor D

下载: 全尺寸图片幻灯片

事实上,基于自然应变的膨胀土损伤模型优势不仅在于模拟精度上的提升,而且可以解决传统应变模式下对于土体本身性质认识的局限性。通过自然应变损伤模型,将非损伤变形和损伤变形进行一定区分,绘制图7的应力变形曲线。

图 7 不同冻融循环下非损伤和损伤变形曲线

Figure 7. Non-damage and damage deformation curves under different freeze-thaw cycles

下载: 全尺寸图片幻灯片

通过图7可以发现,非损伤变形部分呈现逐渐密实的状态,即对应加载过程中不存在土体损伤时的变形状态,变形曲线特征基本上与土体的侧限压缩试验对应。三轴压缩时,由于损伤变形贡献了主要的变形量,因此曲线特征基本上与损伤变形曲线一致。

5. 结论

通过引入自然应变的概念,基于Weibull统计理论和应变等价原理,构建了损伤模型。在模型验证的基础上进一步分析膨胀土的冻融损伤特征。通过研究得出以下结论。

（1）构建的自然应变损伤模型可高精度的模拟膨胀土不同冻融循环下的应力应变曲线特征。遗传算法在模型验证中自动化程度高,可实现参数的自动寻优。

（2）膨胀土初始刚度K₀和参数F₀的冻融劣化明显,与冻融循环次数呈负指数关系。内摩擦角呈现出一定的波动下降趋势,均布系数m有略微增加的趋势。

（3）从不同冻融循环下的损伤因子演化曲线可以看出,冻融损伤主要发生在应力应变曲线前期,且随着冻融循环次数增加,应力应变曲线前期的损伤速度有加快趋势。应力应变后期损伤过程则基本一致。

（4）绘制非损伤变形和损伤变形,非损伤变形呈现逐渐压密的状态,损伤变形贡献了更多的变形量。其变形特征与侧限压缩试验曲线和三轴压缩曲线具有较好的对应性。

图 1 试验加载波形

Figure 1. Loading waveforms in tests

下载: 全尺寸图片幻灯片

图 2 不同循环应力比下Gault黏土竖向应变的发展

Figure 2. Axial strain responses of Gault clay under different CSRs

下载: 全尺寸图片幻灯片

图 3 100,1000,10000圈后Gault黏土永久应变与循环应力比的关系

Figure 3. Permanent axial strain versus CSR for Gault clay after 100, 1000 and 10000 cycles

下载: 全尺寸图片幻灯片

图 4 不同循环应力比下Gault黏土回弹模量的发展

Figure 4. Resilient modulus responses of Gault clay under different CSRs

下载: 全尺寸图片幻灯片

图 5 不同循环应力比下Kimmeridge黏土竖向应变的发展

Figure 5. Axial strain responses of Kimmeridge clay under different CSRs

下载: 全尺寸图片幻灯片

图 6 不同循环应力比下Kimmeridge黏土回弹模量的发展

Figure 6. Resilient modulus responses of Kimmeridge clay under different CSRs

下载: 全尺寸图片幻灯片

表 1 原状硬黏土物理性质

Table 1 Physical properties of intact stiff clay

硬黏土	重度γ/(kN·m^-3)	相对密度G_s	含水率w%	液限w_L/%	塑限w_P/%	塑性指数I_P	初始孔隙比e₀
Gault	19.5	2.59	28.0~30.0	71.0~77.0	28.0~31.0	40~46	0.700~0.750
Kimmeridge	21.0	2.50	16.0~21.0	46.0~55.0	18.0~22.0	28~33	0.430~0.520

下载: 导出CSV

表 2 原状硬黏土原位应力状态预估值

Table 2 Estimated values of in-situ stresses for intact stiff clay

硬黏土	K₀	${σ^{'}}_{v0}$ /kPa	${σ^{'}}_{h0}$ /kPa	${p^{'}}_{0}$ /kPa	q₀/kPa
Gault	1.8	39	70	60	-31
Kimmeridge	1.8	133	239	204	-106

下载: 导出CSV

表 3 不排水循环加载试验方案

Table 3 Summary of undrained cyclic test conditions

试验分组	试验编号	试验类型	$σ_{z}^{ampl}$ /kPa	$τ_{z θ}^{ampl}$ /kPa	CSR	η	N
1	GT-1	三轴	10.5	0	0.088	0	10000
1	GT-2	三轴	19.8	0	0.165	0	10000
1	GT-3	三轴	24.6	0	0.205	0	10000
1	GT-4	三轴	29.4	0	0.245	0	10000
1	GT-5	三轴	34.3	0	0.286	0	10000
1	GT-6	三轴	36.0	0	0.300	0	10000
1	GT-7	三轴	48.0	0	0.400	0	10000
1	GT-8	三轴	54.0	0	0.450	0	10000
1	GT-9	三轴	66.0	0	0.550	0	10000
1	GT-10	三轴	78.0	0	0.650	0	10000
1	GT-11	三轴	96.0	0	0.800	0	2000
2	KT-1	三轴	30.0	0	0.074	0	10000
2	KT-2	三轴	105.0	0	0.257	0	10000
2	KT-3	三轴	150.0	0	0.368	0	10000
3	KH-1	扭剪	30.0	10	0.074	1/3	10000
3	KH-2	扭剪	96.0	32	0.235	1/3	10000
3	KH-3	扭剪	183.0	61	0.449	1/3	308

下载: 导出CSV

参考文献(17)

[1]	刘雪珠, 陈国兴. 轨道交通振动下南京片状细砂的有效应力路径及破坏模式[J]. 岩土力学, 2010, 31(3): 719-726. doi: 10.3969/j.issn.1000-7598.2010.03.009 LIU Xue-zhu, CHEN Guo-xing. Effective stress path and failure pattern of Nanjing fines and with flake-shape structure under vibrating load of rail transit[J]. Rock and Soil Mechanics, 2010, 31(3): 719-726. (in Chinese) doi: 10.3969/j.issn.1000-7598.2010.03.009
[2]	王瑞, 王雷, 胡志平, 等. 交通荷载引起的静偏应力对压实黄土动力特性的影响[J]. 铁道学报, 2019, 41(7): 110-117. doi: 10.3969/j.issn.1001-8360.2019.07.014 WANG Rui, WANG Lei, HU Zhi-ping, et al. Study on dynamic characteristic of compacted loess subjected to static deviatoric stress induced by traffic loading[J]. Journal of the China Railway Society, 2019, 41(7): 110-117. (in Chinese) doi: 10.3969/j.issn.1001-8360.2019.07.014
[3]	陈成, 周正明, 张先伟, 等. 循环荷载作用下泥炭质土动力累积特性试验研究[J]. 岩石力学与工程学报, 2017, 36(5): 1247-1255. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201705020.htm CHEN Cheng, ZHOU Zheng-ming, ZHANG Xian-wei, et al. Experimental study on accumulative behavior of peaty soil under cyclic loading[J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(5): 1247-1255. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201705020.htm
[4]	黄博, 丁浩, 陈云敏. 高速列车荷载作用的动三轴试验模拟[J]. 岩土工程学报, 2011, 33(2): 195-202. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201102007.htm HUANG Bo, DING Hao, CHEN Yun-min. Simulation of high-speed train load by dynamic triaxial tests[J]. Chinese Journal of Geotechnical Engineering, 2011, 33(2): 195-202. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201102007.htm
[5]	蔡袁强, 陈静, 王军. 循环荷载下各向异性软黏土应变-软化模型[J]. 浙江大学学报(工学版), 2008(6): 160-166. https://www.cnki.com.cn/Article/CJFDTOTAL-ZDZC200806032.htm CAI Yuan-qiang, CHEN Jing, WANG Jun. Strain-degradation model for anisotropic soft clay under cyclic loading[J]. Journal of Zhejiang University (Engineering Science), 2008(6): 160-166. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZDZC200806032.htm
[6]	王常晶, 陈云敏. 交通荷载引起的静偏应力对饱和软黏土不排水循环性状影响的试验研究[J]. 岩土工程学报, 2007, 29(11): 1742-1747. doi: 10.3321/j.issn:1000-4548.2007.11.026 WANG Chang-jin, CHEN Yun-min. Study on effect of traffic loading induced static deviator stress on undrained cyclic properties of saturated soft clay[J]. Chinese Journal of Geotechnical Engineering, 2007, 29(11): 1742-1747. (in Chinese) doi: 10.3321/j.issn:1000-4548.2007.11.026
[7]	黄茂松, 李进军, 李兴照. 饱和软黏土的不排水循环累积变形特性[J]. 岩土工程学报, 2006, 28(7): 891-895. doi: 10.3321/j.issn:1000-4548.2006.07.016 HUANG Mao-song, LI Jin-jun, LI Xing-zhao. Cumulative deformation behavior of soft clay in cyclic undrained tests[J]. Chinese Journal of Geotechnical Engineering, 2006, 28(7): 891-895. (in Chinese) doi: 10.3321/j.issn:1000-4548.2006.07.016
[8]	PUPPALA A J, SARIDE S, CHOMTID S. Experimental and modeling studies of permanent strains of subgrade soils[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2009, 135(10): 1379-1389. doi: 10.1061/(ASCE)GT.1943-5606.0000163
[9]	LEKARP F, ISACSSON U, DAWSON A. State of the art I: Resilient response of unbound aggregates[J]. Journal of Transportation Engineering, 2000, 126(1): 66-75. doi: 10.1061/(ASCE)0733-947X(2000)126:1(66)
[10]	肖军华, 许世芹, 韦凯, 等. 主应力轴旋转对地铁荷载作用下软黏土累积变形的影响[J]. 岩土力学, 2013, 34(10): 2938-2945. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201310028.htm XIAO Jun-hua, XU Shi-qin, WEI Kai, et al. Influences of rotation of principal stress axis on accumulative deformation of soft clay under subway cyclic loading[J]. Rock and Soil Mechanics, 2013, 34(10): 2938-2945. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201310028.htm
[11]	钱建固, 王永刚, 张甲峰, 等. 交通动载下饱和软黏土累计变形的不排水循环扭剪试验[J]. 岩土工程学报, 2013, 35(10): 1790-1798. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201310005.htm QIAN Jian-gu, WANG Yong-gang, ZHANG Jia-feng, et al. Undrained cyclic torsion shear tests on permanent deformation responses of soft saturated clay to traffic loadings[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(10): 1790-1798. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201310005.htm
[12]	邓鹏, 郭林, 蔡袁强, 等. 循环荷载下考虑主应力旋转的软土力学响应研究[J]. 岩土力学, 2015(增刊2): 156-164. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2015S2020.htm DENG Peng, GUO Lin, CAI Yuan-qiang, et al. Research on mechanical response of soft clay under cyclic loading involving principle stress rotation[J]. Rock and Soil Mechanics, 2015(S2): 156-164. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2015S2020.htm
[13]	熊焕, 郭林, 蔡袁强. 交通荷载应力路径下砂土地基变形特性研究[J]. 岩土工程学报, 2016, 38(4): 662-669. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201604012.htm XIONG Huan, GUO Lin, CAI Yuan-qiang. Deformation behaviors of sandy subgrade soil under traffic load-induced stress path[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(4): 662-669. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201604012.htm
[14]	BROSSE A M. Study of the Anisotropy of Three British Mudrocks Using a Hollow Cylinder Apparatus[D]. London: Imperial College, 2012.
[15]	WILKINSON S. The microstructure of UK mudrocks[D]. London: Imperial College, 2011.
[16]	HOSSEINI KAMAL R. Experimental Study of the Geotechnical Properties of UK Mudrocks[D]. London: Imperial College, 2012.
[17]	CAI Y Q, GUO L, JARDINE R J, et al. Stress-strain response of soft clay to traffic loading[J]. Géotechnique, 2017, 67(5): 446-451.

施引文献(10)

期刊类型引用(5)

1.	谢柏涵，孔元元，张学飞，孙东彦，王智慧. 冻融循环作用下盐渍土损伤模型的建立及验证. 人民长江. 2025(02): 183-189+200 . 百度学术
2.	孙超，谷鸿雁，郭浩天. 冻融循环次数对粉质黏土强度劣化影响的试验研究. 吉林建筑大学学报. 2024(05): 66-71 . 百度学术
3.	吕思清，朱杰兵，汪斌，祝永锁. 冻融荷载耦合作用下含开口裂隙砂岩宏细观损伤模型研究. 岩石力学与工程学报. 2023(05): 1124-1135 . 百度学术
4.	李新宇，凌贤长，曲娜. 考虑温度效应的冻结膨胀土统计损伤模型. 吉林大学学报(工学版). 2023(08): 2339-2349 . 百度学术
5.	张雁，薛晓雨，杨晓蕴. 低温状态下纤维加筋膨胀土数值模拟试验. 中国科技论文. 2022(12): 1351-1357 . 百度学术