动荷载下填方体压实黄土动应变与动孔压变化规律研究

代倩; 廖红建; 康孝森; 董琪

doi:10.11779/CJGE2021S1043

动荷载下填方体压实黄土动应变与动孔压变化规律研究 English Version

代倩^{1, 2,},
廖红建^1, ,,
康孝森³,
董琪⁴

1.
西安交通大学土木工程系,陕西　西安　710049
2.
西安交通大学城市学院,陕西　西安　710018
3.
长安大学地质工程系,陕西　西安　710064
4.
陕西科技控股集团有限责任公司,陕西　西安　710077

基金项目:

陕西省自然科学基础研究计划项目 2020JQ-988

国家自然科学基金项目 51879212

国家自然科学基金项目 41630639

陕西省重点研发计划项目 2019KWZ-09

中央高校基金项目 300102261303

详细信息

作者简介:
代倩（1989— ）,女,博士研究生,主要从事黄土强度变形方面的试验研究。E-mail：daiqian68@stu.xjtu.edu.cn

通讯作者:
廖红建, E-mail：hjliao@mail.xjtu.edu.cn

中图分类号: TU43
计量
- 文章访问数: 178
- HTML全文浏览量: 18
- PDF下载量: 114
出版历程
- 收稿日期: 2020-12-14
- 网络出版日期: 2022-12-05
- 刊出日期: 2021-06-30

Behaviors of dynamic strain and pore pressure of compacted loess in loess-filled foundation induced by dynamic loading

1.
Department of Civil Engineering, Xi'an Jiaotong University, Xi'an 710049, China
2.
Xi'an Jiaotong University College, Xi'an 710018, China
3.
Department of Geological Engineering, Chang'an University, Xi'an 710064, China
4.
Shaanxi Science and Technology Holding Group, Xi'an 710077, China

摘要

摘要: 某黄土填方体不同深度含水率差异明显,部分层位体积含水率达到50%,高含水率影响该填方体在长期动荷载作用下的变形与应力状态,需开展动荷载下填方体压实黄土动应变与动孔压变化研究。制备接近填方体实际干密实度的三轴试样,开展不同围压与不同循环应力比条件下饱和压实黄土动三轴试验,研究循环应力比与循环次数对饱和压实黄土塑性应变累积、应变变化率、应变幅值、孔压累积、孔压变化率与孔压幅值的影响规律,进一步基于边界面塑性理论揭示循环应力比和围压的这一影响机制。研究结果有利于黄土填方体动力分析,对该高填方体长期变形评估具有重要参考意义。
- 黄土高填方 /
- 压实黄土 /
- 动三轴试验 /
- 动孔压 /
- 动应变
Abstract: The water content of different layers is various in a typical high loess-filled foundation.The water content of some layers even comes to 50%, which may greatly affect the strain and stress state of the high loess-filled foundation under dynamic loading.Thus, it's necessary to study the dynamic strain and pore pressure of the compacted loess in the high loess-filled foundation.The compacted loess samples are made according to the actual dry density of the filled-loess foundation, then series of dynamic triaxial tests are conducted on the compacted loess samples subjected to different confining pressures and cyclic stress ratios.Based on the test results, the influences of cyclic stress ratio and cyclic number on deviatoric strain, deviatoric strain rate, deviatoric strain amplitude, pore pressure, pore pressure rate, and pore pressure amplitude are analyzed.The influences are further explored in the framework of the bounding surface plasticity theory.The conclusions are meaningful in analyzing the high filled-loess foundation under dynamic loading and in evaluating the deformation of the loess-filled foundation.
- high loess-filled foundation /
- compacted loess /
- dynamic triaxial test /
- dynamic pore pressure /
- dynamic strain

HTML全文

0. 引言

土工离心模型试验技术是一项崭新的土工物理模型技术，通过施加在模型上的离心惯性力使模型的重度变大，从而使模型的应力与原型一致，这样就可以用模型反映、表示原型。离心模型是各类物理模型中相似性最好的模型，在国内外受到广泛重视，试验技术飞速发展与进步，研究内容涉及几乎所有的岩土工程研究领域，已成为岩土工程技术研究中的最主要、最有效、最广泛应用的研究手段^[1]。

Biot动力固结方程是土力学的基本方程，包括动力平衡方程和连续方程。从连续方程出发，得到固结与孔隙水压力的扩散问题的时间比尺关系t_p/t_m=n^2[2]，这个关系的前提假设是，模型的渗透系数k_m是原型渗透系数k_p的n倍，模型与原型的固结系数相等，这一假设存在争议，未得到试验完全证实。

Pokrovsky & Fyodorov认为离心模型的水力坡降i_m与原型的水力坡降i_p相等，模型的渗透系数k_m是原型渗透系数k_p的n倍，得到了t_p/t_m=n²的时间比尺关系；而Cargill等^[3]则认为离心模型的水压力坡降i_m是原型的水压力坡降i_p的n倍，模型的渗透系数k_m与原型渗透系数k_p相等，令人奇怪的是，他们利用相互矛盾、相互对立的假定得到了相同的结论。

其实，Cargill等把达西定律中的“水力梯度（h/L）”用“水压力梯度（p/L）”来代替是错误的，因为水力梯度（无量纲）与水压力梯度（有量纲）并不相同。因而，模型渗透系数k_m与原型渗透系数k_p的关系，成了问题的关键。即土的渗透系数（达西渗透系数）在高重力场下变不变？是否随加速度而增大？

从不同的研究目的出发，世界各国的离心模型试验专家进行了不懈的努力。Goodings^[4]通过试验研究了灰坝的浸润线，Croce等^[5]研究了土样的固结问题，Singh等^[6]和Khalifa等^[7]都研究过土的渗透系数随加速度的变化规律。令人不解的是，Singh等的结果表明，土的渗透系数随加速度而增大，而Khalifa等的结果却表明，在50g以下，土的渗透系数不随加速度变化。Singh等的试验装置不能防止水的蒸发，结果令人怀疑；Khalifa等的试验只进行到50g，也难以说明问题，不能令人信服。Sharma等^[8]分析了离心机半径对变水头试验测得的渗透系数的影响。隋海宾等^[9]总结了国外离心模型试验中土的渗透相似性研究成果，但没有自己的试验。Anderson等^[10]研制了在离心机上进行土的渗透系数试验的装置，可以大大缩短试验时间。

国内外离心模型试验发展很快，离心模型试验中固结与孔隙水压力扩散问题的时间比尺问题一直未能得到有效的试验验证。固结与孔隙水压力扩散问题的时间比尺关系t_p/t_m=n²正确与否，归根结底就是高重力场条件下土的渗透特性与普通重力场条件下是否一致？如果不一致，土的渗透特性随离心加速度是怎样变化的？本文通过离心模型渗透试验，研究了土的渗透系数随离心加速度的变化规律，验证离心模型渗透相似理论，澄清固结与孔隙水压力扩散问题的时间比尺关系的一些模糊不清问题，对离心模型试验技术的发展具有十分重要的理论意义和应用价值。

1. 渗透相似理论

水在土体孔隙中的渗透采用达西定律描述，

$v = k\frac{{\Delta h}}{L} = ki \text{，}$

(1)

式中：v为渗透速度； $k$ 为渗透系数，其物理意义是当水力梯度等于1时的渗透速度；L为渗径长度；Δh为试样两端的水位差；i为水力梯度，它是沿渗透方向单位距离的水头损失，无量纲。

定义任一模型变量x_m与原型变量x_p的相似常数：

${\eta _x} = {x_p}/x{}_m 。$

(2)

将式（2）代入式（1），可以得到

${\eta _v}v = {\eta _k}{\eta _i}ki 。$

(3)

达西定律可以描述原型和模型，因此它对原型与模型来说都应当是正确的。只有当所有的η项系数都相等时，式（3）才与原方程式（1）相同，也就是，

${\eta _v} = {\eta _k}{\eta _i} 。$

(4)

离心模型试验中，模型的渗透速度是原型的n倍，即η_v=1/n，而水力梯度无量纲，则η_i=1，得到

${\eta _k} = 1/n 。$

(5)

Muskat^[11]、Lambe等^[12]分别给出了渗透系数与流体重度的关系式：

$k = K\frac{\gamma }{\mu } = K\frac{{\rho g}}{\mu } 。$

(6)

式中：K为土体的内在渗透系数，是颗粒形状、直径和填料的函数；μ为流体的动力黏度；ρ为流体的密度；g为重力加速度。

由式（6）可以明显看出，渗透系数随土体的内在渗透系数、流体的动力黏度、流体的密度、重力加速度变化，模型和原型土体的内在渗透系数、流体的动力黏度、流体的密度相同，重力加速度为n倍，即，η_K=η_μ=η_ρ=1，η_g=1/n，因此，可得η_k=1/n。

2. 渗透离心模型试验

2.1 试验设备

试验在南京水利科学研究院NHRI400 g·t土工离心机上进行。该机最大半径5.5 m，吊篮平台的尺寸为1100 mm×1100 mm，最大加速度200g，最大负荷2 t，容量400 g·t。为在离心机中进行渗透试验，研制出离心机渗透仪，渗透仪装置由渗透容器、变水头管、水头测量装置等组成（图 1）。渗透容器由金属圆筒（内径100 mm、高100 mm）、透水石、上盖和下盖组成，变水头管由连接段（内径40 mm）和测量段（内径20 mm）组成，水头测量装置由激光位移传感器和浮标组成。

图 1 离心机渗透仪

Figure 1. Permeameter of centrifuge model tests

下载: 全尺寸图片幻灯片

2.2 试验方法

试验土样为低液限黏土（CL）、低液限粉土（ML）、粉土质砂（SM）3种，其物理性质指标见表 1。

表 1 土样的物理性质指标

Table 1. Physical properties of soils

土名		低液限黏土	低液限粉土	粉土质砂
含水率/%		29.0	26.2	20.5
密度/(g·m^-3)		1.94	2.08	2.10
试样高度/mm		40	36	69
液限/%		42.2	24.8
塑性指数		20.7	6.6
颗粒组成 /%	砂粒	12.8	39.0	55.3
	粉粒	44.7	41.6	38.5
	黏粒	42.5	19.4	6.2

下载: 导出CSV

| 显示表格

将土样制备成饱和泥浆，放入试样筒内，采用预压方法进行固结成型，固结完成后土样的含水率、密度、试样高度见表 1。将安装好的渗透仪装置放置在模型箱内，以控制下水位，从而控制渗透水头。将渗透仪装置和模型箱放入离心机吊篮，接好激光位移传感器，启动离心机到设定加速度，通过位移传感器测试上水位下降过程。每种土样的试验加速度分别1g（常规）、10g、25g、50g、75g、100g、150g，每种工况进行2组平行试验。

2.3 渗透水头修正

离心惯性加速度场是通过离心机的旋转而产生的，离心机旋转产生的离心惯性加速场是沿旋转中轴形成的一个个圆形的柱面，在每个等半径的圆柱面上，其离心惯性加速场的惯性势是相等的，离心加速度也相同，离心加速度的大小与半径成正比。这样，在模型的不同高度上，所受到的离心惯性力是不同的。那么离心模型中，上、下水位差并不是真正的渗透水头，而必须进行修正。

如图 2所示，任意旋转半径R处的水柱dR，对试样中心产生的渗透压力增量dp可表示为，

$\text{d}p={\rho }_{\text{w}}\frac{{a}_{0}}{{R}_{0}}R\text{d}R\text{ }。$

(7)

图 2 渗透水头修正示意图

Figure 2. Schematic diagram of seepage water head correction

下载: 全尺寸图片幻灯片

式中dp为试样中心的渗透压力增量； ${\rho _{\text{w}}}$ 为水的密度；a₀为试样中心的加速度；R₀为试样中心的旋转半径；R为任意水柱的旋转半径；dR为任意水柱的高度。

对式（7）积分，可得土样中心产生的渗透压力为

$p = \int_{{R_1}}^{{R_2}} {{\rho _{\text{w}}}\frac{{{a_0}}}{{{R_0}}}R{\text{d}}R} = \frac{1}{2}{\rho _{\text{w}}}\frac{{{a_0}}}{{{R_0}}}(R_2^2 - R_1^2) 。$

(8)

式中，p为试样中心的渗透压力；R₂为下水位的旋转半径；R₁为上水位的旋转半径；其他符号意义同前。

因此，渗透水头为

$\Delta h = \frac{p}{{{\rho _{\text{w}}}{a_0}}} = \frac{{R_2^2 - R_1^2}}{{2{R_0}}} \text{，}$

(9)

式中，Δh为试样的渗透水头，其他符号意义同前。

变水头渗透系数按式（10）计算

${k_T} = \frac{{BL}}{{A({t_2} - {t_1})}}\ln \frac{{\Delta {h_1}}}{{\Delta {h_2}}} 。$

(10)

式中，k_T为水温为T℃时试样的渗透系数；B为测量段变水头管的断面积；L为渗径，即试样高度；A为试样的断面积；t₁，t₂分别为测读水头的起始和终止时间；Δh₁，Δh₂分别为起始和终止渗透水头。

标准温度（20℃）下的渗透系数按式（11）计算

${k_{20}} = {k_T}\frac{{{\eta _T}}}{{{\eta _{20}}}} 。$

(11)

式中， ${k_{20}}$ 为标准温度时试样的渗透系数； ${\eta _T}$ 为T℃时水的动力黏滞系数； ${\eta _{20}}$ 为20℃时水的动力黏滞系数。

3. 试验结果分析

3.1 渗透系数的变化规律

表 2列出了不同加速度下3种土样的渗透系数试验结果，图 3给出了3种土样的渗透系数随加速度的变化，从图中可以看出，3种土样的渗透系数随加速度的增加而线性增大，表明渗透系数与加速度呈正线性关系。把图 3的纵坐标改为ng时渗透系数与1g时渗透系数之比k_m/k_p，如图 4所示，可以看出，k_m/k_p与加速度具有良好的线性变化关系，可表示为

${k_{\text{m}}}/{k_{\text{p}}} = \chi n \text{，}$

(12)

表 2 渗透系数试验结果

Table 2. Test results of permeability coefficient (cm/s)

土名	加速度/g
土名	1		10	25	50	75	100	150
低液限黏土	1.21×10^-8	1.24×10^-8	1.28×10^-7	3.58×10^-7	6.16×10^-7	8.50×10^-7	1.40×10^-6	1.79×10^-6
低液限黏土	1.29×10^-8	1.25×10^-8	1.14×10^-7	2.97×10^-7	5.26×10^-7	7.57×10^-7	1.26×10^-6	1.73×10^-6
低液限粉土	1.21×10^-7	1.37×10^-7	1.21×10^-6	2.92×10^-6	6.06×10^-6	8.48×10^-6	1.07×10^-5	1.73×10^-5
低液限粉土	1.26×10^-7	1.18×10^-7	1.17×10^-6	2.67×10^-6	5.77×10^-6	7.91×10^-6	9.97×10^-6	1.52×10^-5
粉土质砂	1.32×10^-6	1.18×10^-6	1.29×10^-5	2.99×10^-5	5.99×10^-5	9.44×10^-5	1.22×10^-4	1.85×10^-4
粉土质砂	1.50×10^-6	1.09×10^-6	1.13×10^-5	3.10×10^-5	6.16×10^-5	9.81×10^-5	1.26×10^-4	1.79×10^-4

下载: 导出CSV

| 显示表格

图 3 渗透系数随加速度的变化关系

Figure 3. Variation of permeability coefficient with acceleration

下载: 全尺寸图片幻灯片

图 4 k_m/k_p随加速度的变化关系

Figure 4. Relationship between k_m/k_p and acceleration

下载: 全尺寸图片幻灯片

式中， $\chi$ 为比例系数。试验数据拟合得出，3种土样的 $\chi$ 分别0.984，0.937，0.992。除去试验仪器以及固结等因素的影响，比例系数 $\chi$ ≅1。因此，试验证明了离心模型试验渗透系数比尺为η_k=1/n。Singh等^[6]利用小离心机变水头试验，也得出了类似的试验结果。

3.2 离心模型中达西定律的适用性

进行渗透试验时，渗透水头随时间的变化直接反映了试样的渗透稳定情况，图 5给出了不同加速度条件下3种土样的渗透水头过程线。试验结果表明，3种土样的渗透水头随时间的延长而逐渐下降，且下降速率基本一致，说明试样的渗透稳定。下降速率与离心机加速度和土性有关，加速度越大，下降速率越大，土样黏粒含量越大，下降速率越慢。

图 5 不同加速度下渗透水头随时间的变化

Figure 5. Variation of seepage water head with time under different accelerations

下载: 全尺寸图片幻灯片

图 6给出了3种土样的渗透速度随加速度与水力梯度乘积的关系曲线，从图中可以看出，3种土样的渗透速度随加速度与水力梯度乘积的增加而线性增大，且线性关系相当好，增长速率与试样的黏粒含量有关，试样黏粒含量越大，增长速率越大，增长速率即为1g条件下试样的渗透系数。表 3列出了3种土样渗透试验最大加速度与水力梯度乘积，表明在这个范围内，离心模型试验渗透符合达西定律，试样黏粒含量越大，加速度与水力梯度乘积也越大。

图 6 渗透速度随加速度与水力梯度乘积的变化

Figure 6. Variation of seepage velocity with product of acceleration and hydraulic gradient

下载: 全尺寸图片幻灯片

表 3 3种土样渗透试验最大加速度与水力梯度乘积

Table 3. Products of maximum acceleration and hydraulic gradients of three soil samples in permeability tests

土名	d₁₀/mm	k/(cm⋅s^-1)	v/(cm⋅s^-1)	试验最大加速度×水力梯度
低液限黏土	0.0014	1.22×10^-8	4.81×10^-5	4033
低液限粉土	0.0021	1.15×10^-7	3.50×10^-4	3037
粉土质砂	0.0112	1.24×10^-6	8.14×10^-4	661

下载: 导出CSV

| 显示表格

4. 结语

离心模型渗透试验结果表明，渗透系数随加速度的增加而线性增大，离心模型试验渗透符合达西定律，模型的渗透系数与原型渗透系数之比k_m/k_p随加速度成比例线性增大，比例系数≈1，从而验证了离心模型试验渗透系数比尺为η_k=1/n。虽然国外也得出了类似的试验结果，但我们试验的渗透水头比他们的大，土样类型也多。研究成果对离心模型试验具有重大理论和应用价值。

图 1 黄土取样点

Figure 1. Sampling location of loess

下载: 全尺寸图片幻灯片

图 2 基本指标的确定

Figure 2. Determination of basic indices

下载: 全尺寸图片幻灯片

图 3 压实黄土塑性偏应变与循环次数关系

Figure 3. Relationship between plastic deviatoric strain and cyclic number of compacted loess

下载: 全尺寸图片幻灯片

图 4 塑性偏应变变化率与循环次数关系

Figure 4. Relationship between plastic deviatoric strain rate and cyclic number

下载: 全尺寸图片幻灯片

图 5 偏应变波动幅值与循环次数关系

Figure 5. Relationship between amplitude of plastic deviatoric strain and cyclic number

下载: 全尺寸图片幻灯片

图 6 不同循环应力比下压实黄土孔隙水压力与循环次数

Figure 6. Pore pressures and cyclic number of compacted loess subjected to different CSRs

下载: 全尺寸图片幻灯片

图 7 孔压变化率与循环次数关系

Figure 7. Relationship between rate of pore pressure and cyclic number

下载: 全尺寸图片幻灯片

图 8 孔压波动幅值与循环次数关系

Figure 8. Relationship between amplitude of pore pressure and cyclic number

下载: 全尺寸图片幻灯片

图 9 有效应力路径与边界面和临界状态线的关系

Figure 9. Relationship among effective stress paths, bounding surface and critical state line

下载: 全尺寸图片幻灯片

图 10 半周期内塑性模量变化

Figure 10. Change of plastic modulus in half cyclic number

下载: 全尺寸图片幻灯片

图 11 动杨氏模量与循环加载次数关系

Figure 11. Relationship between dynamic Young's modulus and cyclic number

下载: 全尺寸图片幻灯片

图 12 循环应力比对孔压和偏应变的影响

Figure 12. Influences of cyclic stress ratio on pore pressure and plastic deviatoric strain

下载: 全尺寸图片幻灯片

参考文献(19)

[1]	郑建国, 曹杰, 张继文, 等. 基于离心模型试验的黄土高填方沉降影响因素分析[J]. 岩石力学与工程学报, 2019, 38(3): 560-571. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201903012.htm ZHENG Jian-guo, CAO Jie, ZHANG Ji-wen, et al. Analysis of influencing factors of high loessfilled foundations based on centrifugal model tests[J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(3): 560-571. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201903012.htm
[2]	葛苗苗, 李宁, 张炜, 等. 黄土高填方沉降规律分析及工后沉降反演预测[J]. 岩石力学与工程学报, 2017, 36(3): 745-753. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201703024.htm GE Miao-miao, LI Ning, ZHANG Wei, et al. Settlement behavior and inverse prediction of postconstruction settlement of high filled loess embankment[J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(3): 745-753. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201703024.htm
[3]	段旭, 董琪, 门玉明, 等. 黄土沟壑高填方工后地下水与土体含水率变化研究[J]. 岩土工程学报, 2018, 40(9): 1753-1758. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201809027.htm DUAN Xu, DONG Qi, MEN Yu-ming, et al. Change of groundwater and water content of loess high fill in gullyregions[J]. ChineseJournalofGeotechnical Engineering, 2018, 40(9): 1753-1758. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201809027.htm
[4]	王兰民, 刘红玫, 李兰, 等. 饱和黄土液化机理与特性的试验研究[J]. 岩土工程学报, 2000, 22(1): 92-97. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200001015.htm WANG Lan-min, LIU Hong-mei, LI Lan, et al. Laboratory study on the mechanism and behaviors of saturated loess liquefaction[J]. ChineseJournalofGeotechnical Engineering, 2000, 22(1): 89-94. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200001015.htm
[5]	骆亚生, 李瑞, 田堪良. 非饱和黄土动力特性试验方法研究[J]. 地下空间与工程学报, 2007(6): 1041-1046. doi: 10.3969/j.issn.1673-0836.2007.06.013 LUO Ya-sheng, LI Rui, TIAN Kan-liang. Research on method of dynamic characteristics tests for unsaturated loess[J]. Chinese Journal of Underground Space and Engineering, 2007(6): 1041-1046. (in Chinese) doi: 10.3969/j.issn.1673-0836.2007.06.013
[6]	孙海妹, 王兰民, 王平, 等. 振动荷载作用下黄土应变增长特征试验研究[J]. 地下空间与工程学报, 2011, 7(6): 1103-1107. doi: 10.3969/j.issn.1673-0836.2011.06.010 SUN Hai-mei, WANG Lan-min, WANG Ping, etal. Experimentalstudyonthedeformation characteristics of loess under cyclic load[J]. Chinese Journal of Underground Space and Engineering, 2011, 7(6): 1103-1107. (in Chinese) doi: 10.3969/j.issn.1673-0836.2011.06.010
[7]	WANG S, ZHONG Z L, LIU X R, et al. Influences of principal stress rotation on the deformation of saturated loess under traffic loading[J]. KSCE Journal of Civil Engineering, 2019, 23(5): 2036-2048.
[8]	HU Z P, WANG R, REN X, et al. Permanent deformation behavior of compacted loess under long-term traffic loading[J]. Journal of Materials in Civil Engineering, 2019, 31(8): 04019157.
[9]	廖红建, 肖正华, 刘健. 动载下饱和重塑黄土的骨干曲线变化研究[J]. 岩土力学, 2011, 32(2): 375-379. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201102013.htm LIAO Hong-jian, XIAO Zheng-hua, LIU Jian. Study of variation of backbone curve of saturated remolded loess under dynamic loading[J]. Rock and Soil Mechanics, 2011, 32(2): 375-379(in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201102013.htm
[10]	XIAO Z H, HAN B, LIAO H J, et al. Experiment study on dynamic strengthofloessunderrepeatedload[J]. International Journal of Modern Physics B, Singapore, 2008, 22(31/32): 5825-5830.
[11]	王谦, 马金莲, 马海萍, 等. 饱和黄土动剪切模量和阻尼比的试验研究[J]. 岩石力学与工程学报, 2019, 38(9): 1919-1927. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201909020.htm WANG Qian, MA Jin-lian, MA Haiping, et al. Dynamic shear modulus and damping ratio of saturated loess[J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(9): 1919-1927. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201909020.htm
[12]	ZHANG Z L, WANG T, WU S R, et al. Dynamics stressstrain behavior of Tianshui soils[J]. Landslides, 2017, 14(1): 323-335.
[13]	SONG B H, TSINARIS A, ANASTASIADIS A, et al. Small-strain stiffness and damping of Lanzhou loess[J]. Soil Dynamics&Earthquake Engineering, 2017, 95: 96-105.
[14]	杨利国, 骆亚生, 李焱, 等. 初始应力条件对压实黄土动强度影响的研究[J]. 岩土力学, 2010, 31(1): 87-91. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201001017.htm YANG Li-guo, LUO Ya-sheng, LI Yan, et al. Research on effect of initial stress conditions on dynamic strength of compacted loess[J]. Rock and Soil Mechanics, 2010, 31(1): 87-92. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201001017.htm
[15]	王军海, 刘亚明. 基于动三轴试验的压实黄土动强度特性研究[J]. 地震工程学报, 2016, 38(3): 439-444. https://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ201603015.htm WANG Jun-hai, LIU Ya-ming. Dynamicstrength properties of compacted loessbased on dynamic triaxialtest[J]. China Earthquake Engineering Journal, 2016, 38(3): 439-444. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ201603015.htm
[16]	刘保健, 张晓荣, 程海涛. 应变控制下压实黄土的动三轴试验研究[J]. 岩土力学, 2007(6): 1073-1076. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200706001.htm LIU Bao-jian, ZHANG Xiao-rong, CHENG Hai-tao. Study on compacted loess under strain control at dynamic triaxialtest[J]. Rock and Soil Mechanics, 2007, 28(6): 1073-1076. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200706001.htm
[17]	段旭, 门玉明, 董琪, 等. 黄土沟壑高填方分层沉降监测与深层浸水试验[J]. 水土保持通报, 2017, 37(3): 222-226. https://www.cnki.com.cn/Article/CJFDTOTAL-STTB201703037.htm DUAN Xu, MENG Yu-ming, DONG Qi, et al. Layered settlement monitoring and soaking test of loess high fill in gully region, 2017, 37(3): 222-226. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-STTB201703037.htm
[18]	KANG X S, LIAO H J. Bounding surface plasticity model for jointed soft rocks considering overconsolidation and structural decay[J]. Computers and Geotechnics, 2019, 108: 295-307.
[19]	KANG X S, ZHANG J H, NING C M, LIAO H J. A bounding surface plasticity model for normally consolidated volcanic cohesive soils subjected to cyclic loading[C]//7th International Conference on Earthquake Geotechnical Engineering, 2019, Rome.