黄土盾构隧道开挖围岩扰动特性模型试验研究

韩兴博; 陈子明; 叶飞; 梁晓明; 冯浩岚; 夏天晗

doi:10.11779/CJGE20230054

黄土盾构隧道开挖围岩扰动特性模型试验研究 English Version

韩兴博^{1, 2,},
陈子明¹,
叶飞^{1, 2, ,},
梁晓明³,
冯浩岚¹,
夏天晗¹

1.
长安大学公路学院, 陕西西安 710064
2.
长安大学公路大型结构安全教育部工程研究中心, 陕西西安 710064
3.
西南交通大学土木工程学院, 四川成都 610031

基金项目:

国家自然科学基金项目 52108360

长安大学中央高校基本科研业务费专项资金资助项目 300102213206

中国博士后科学基金项目 2020M683398

详细信息

作者简介:
韩兴博(1991—)，男，博士，副教授，主要从事黄土盾构隧道结构施工及运营期力学响应等方面的研究工作。E-mail: xingbo.han@chd.edu.cn

通讯作者:
叶飞, E-mail：xianyefei@126.com

中图分类号: TU411
计量
- 文章访问数: 400
- HTML全文浏览量: 50
- PDF下载量: 106
出版历程
- 收稿日期: 2023-01-29
- 网络出版日期: 2023-06-01
- 刊出日期: 2024-04-30

Model tests on disturbance characteristics of surrounding rock of loess shield tunnels during excavation

1.
School of Highway Engineering, Chang'an University, Xi'an 710064, China
2.
Research Center on Safety of Highway Large Structure Engineering, Ministry of Education, Chang'an University, Xi'an 710064, China
3.
School of Civil Engineering, Southwest Jiaotong University, Chengdu 610031, China

摘要

摘要: 为了探究不同含水率和埋深等工况下的黄土盾构隧道开挖对围岩的扰动作用特征，采用相似模型试验，首先通过三轴试验测定了不同含水率下原状黄土物理参数及强度指标，通过大量试配得到了不同含水率下与原状黄土强度指标参数相似的试验用土。其后，考虑盾构隧道盾尾间隙特征，开展不同含水率及埋深等组合工况下盾构开挖模拟试验，通过微型土压力传感器、百分表等监测元件以及数字图像技术等手段，分析了不同工况下地层变形、地表沉降和围岩应力的变化规律。研究发现：含水率较低时，围岩的自稳能力强；含水率越高，埋深越大，开挖对地层的扰动作用也越大；地表最大沉降值与含水率呈正相关，与埋深呈负相关；围岩剩余应力随含水率和埋深的增大而增大，但当含水率增大至一定值（26.6%）时，围岩剩余应力骤减；基于不同含水率下的围岩应力及变形发展规律，可将黄土含水率分为自稳、形变以及松动含水率，基于含水率和埋深引起的围岩松动情况，可将不同工况下的黄土盾构隧道围岩压力作用模式按照形变压力和松动压力计算，以更加适应黄土地层的特性。
- 黄土盾构隧道 /
- 相似模型试验 /
- 含水率 /
- 埋深 /
- 围岩变形 /
- 应力
Abstract: In order to explore the disturbance characteristics of the surrounding rock by excavation of loess shield tunnels under different water contents and buried depths, the similarity model tests are carried out. Firstly, the physical parameters and strength indexes of the undisturbed loess under different water contents are measured by the triaxial tests. Through a large number of tests, the test soil with similar strength indexes of the undisturbed loess with different water contents is obtained. Then, considering the characteristics of shield tail gap of shield tunnels, the simulated tests on the shield excavation under different water contents and buried depths are carried out. The variation laws of stratum deformation, surface settlement and rock stress under different working conditions are analyzed by means of the monitoring components such as micro-earth pressure sensors, dial indicator and digital image technology. It is found that the self-stabilizing capability of the surrounding rock is strong when the water content is low. The higher the water content, the greater the buried depth, the greater the disturbance of excavation to the stratum. The maximum ground settlement is positively correlated with the water content and negatively correlated with the buried depth. The residual stress of the surrounding rock increases with the increase of the water content and buried depth, but when the water content increases to a certain value (26.6 %), the residual stress of the surrounding rock decreases sharply. Based on the development law of the stress and deformation of the surrounding rock under different water contents, the water content of loess can be divided into self-stabilization, deformation and loosening water contents. Based on the loosening of the surrounding rock caused by the water content and buried depth, the pressure action mode of the surrounding rock of loess shield tunnel under different working conditions can be calculated according to the deformation pressure and loosening pressure, so as to better adapt to the characteristics of loess strata.
- loess shield tunnel /
- similarity model test /
- water content /
- buried depth /
- surrounding rock deformation /
- stress

HTML全文

0. 引言

近年来，海上风电迅速成为重点发展的能源形式之一。筒型基础凭借其成本低、工期短且对地基扰动相对较小的优势，已成为中国海上风电的重要基础形式之一。竖向承载力验算是海上风电基础稳定性校核的首要环节，筒型基础通常自重较大且长期承受海洋平台结构自重及一定程度的竖向荷载作用，因此了解其地基竖向承载特性对于实际工程设计而言十分重要。

针对筒型基础地基竖向承载力的研究较多，已有较多国内外的学者运用试验以及数值模拟等方式进行研究。Houlsby等^[1]和Byrne等^[2]在21世纪初期探讨了传统的吸力筒型基础作为海上风电基础的应用前景。Fu等^[3]、Gourvenec等^[4-5]、Chen等^[6]、刘梅梅等^[7]通过开展一系列室内模型试验及离心模型试验，得到了筒型基础的竖向承载模式、承载力p-s曲线、基础各部分载荷分担情况及筒内外土压力分布情况。Hung等^[8]、Vulpe^[9]、Mehravar等^[10]、Mana等^[11]基于大量有限元数值模拟，讨论了多种因素对地基竖向承载力的影响，得到了相关规律及半经验计算公式。然而由于筒型基础筒-土相互作用的特殊性，其地基竖向承载力的理论计算一直未能完全解决。

现有的筒型基础地基竖向承载力计算方法大多是利用极限分析上限法推导，但该方法得到计算结果往往高于地基承载力的真值，在工程应用中因偏于危险而无法推广。极限分析下限法从理论上更符合工程设计中“安全”概念，但传统方式构建容许应力场困难且具有主观性，导致下限法运用范围受限。直到1988年Sloan^[12]提出了适用于平面应变问题的有限元极限分析下限法，为极限分析下限法的应用提供了全新的思路，极大推动了其在岩土工程中的发展与应用。

由于筒型基础相互作用的特殊性，对于筒型基础的有限元极限分析下限法鲜有学者涉足，本文将筒土相互作用转化为筒土间摩擦的不等式约束条件，建立了筒型基础地基竖向承载力的下限解求解程序。利用离心模型试验结果验证了本方法的正确性，在此基础上分析了土体强度不均匀性、筒型基础长径比及筒外壁-土间摩擦对筒型基础地基竖向承载力的影响，提出了综合考虑以上因素的筒型基础地基竖向承载力计算公式。

1. 轴对称有限元极限分析下限法

本节在Sloan^[12]的有限元下限法求解格式基础上，结合轴对称问题相关理论，提出了针对筒型基础这一轴对称结构的单元平衡条件、应力间断面约束条件、土体边界条件、筒-土界面约束条件、土体屈服条件及相应的目标函数求解格式。

1.1 单元平衡约束条件

采用单元-结构面离散方式将结构物离散为有限个单元体，假设单元为刚性单元。相比于平面应变中以（ ${\sigma _x}$ ， ${\sigma _y}$ ， ${\tau _{xy}}$ ）为未知变量，轴对称有限元法引入环向应力 ${\sigma _\theta }$ ，将节点应力（ ${\sigma _r}$ ， ${\sigma _z}$ ， ${\tau _{rz}}$ ， ${\sigma _\theta }$ ）作为基本未知变量，模拟应力场的典型三角形单元如图 1所示。允许沿相邻三角形单元之间的所有界面发生静态允许的应力间断，每个节点都与唯一单元相关联。对于轴对称结构，各子午面上的应力状态均相同，假定当r-z平面（即子午面）上的应力满足极限平衡条件时，地基达到极限承载状态。

图 1 土体有限元离散

Figure 1. Finite element discretization of soil mass

下载: 全尺寸图片幻灯片

对于筒型基础地基竖向承载力分析模型，土体仅受重力作用，柱坐标系下单元分析示意图如图 2所示。

图 2 柱坐标系中单元体

Figure 2. Unit in a bucket coordinate system

下载: 全尺寸图片幻灯片

柱坐标系中单元体平衡方程可表示为

$\left.\begin{array}{l}\frac{\partial \sigma }{{\partial }_{r}}+\frac{\partial {\tau }_{rz}}{{\partial }_{z}}\text{+}\frac{{\sigma }_{r}-{\sigma }_{\theta }}{r}=\text{0 }\text{，}\\ \frac{\partial {\tau }_{rz}}{{\partial }_{r}}\text{+}\frac{\partial {\sigma }_{z}}{{\partial }_{z}}\text{+}\frac{{\tau }_{rz}}{r}=\gamma \text{ }。\end{array}\right\}$

(1)

假设土单元为线性单元，单元中3个节点满足平衡方程即代表单元平衡。在此基础上，将式（1）平衡方程转化为离散形式产生如下等式约束：

${\left[ {A_{{\text{equil}}}^{\text{e}}} \right]_{2 \times 12}}{\left\{ {{\sigma ^{\text{e}}}} \right\}_{12 \times 1}} = {\left\{ {b_{{\text{equil}}}^{\text{e}}} \right\}_{2 \times 1}} 。$

(2)

式中： ${\left[{A}_{\text{equil}}^{\text{e}}\right]}_{\text{2}\times \text{12}}\text{ }$

${\text{ = }}\left[ \begin{gathered} \frac{{{z_{23}}}}{{2A}} + \frac{1}{{3\bar r}}{\text{ 0 }}\frac{{{r_{32}}}}{{2A}}{\text{ }} - \frac{1}{{3\bar r}}{\text{ }} \hfill \\ {\text{ }}0{\text{ }}\frac{{{r_{32}}}}{{2A}}{\text{ }}\frac{{{z_{23}}}}{{2A}} + \frac{1}{{3\bar r}}{\text{ }}0 \hfill \\ \end{gathered} \right.$

$\left.\begin{array}{cc}\frac{{z}_{31}}{2A}+\frac{1}{3\overline{r}}& 0\\ 0& \frac{{r}_{13}}{2A}\end{array}\text{ }\begin{array}{cc}\frac{{r}_{13}}{2A}& -\frac{1}{3\overline{r}}\\ \frac{{z}_{31}}{2A}+\frac{1}{3\overline{r}}& \text{0}\end{array}\text{ }\begin{array}{cc}\frac{{z}_{12}}{2A}+\frac{1}{3\overline{r}}& \text{0}\\ 0& \frac{{r}_{21}}{2A}\end{array}\text{ }\begin{array}{cc}\frac{{r}_{21}}{2A}& -\frac{1}{3\overline{r}}\\ \frac{{z}_{12}}{2A}+\frac{1}{3\overline{r}}& 0\end{array}\right]\text{，} \\ {\left\{ {{\sigma ^{\text{e}}}} \right\}^{\text{T}}} = \{ \sigma _{r{\text{,1}}}^{\text{e}}{\text{ }}\sigma _{z{\text{,1}}}^{\text{e}}{\text{ }}\tau _{rz{\text{,1}}}^{\text{e}}{\text{ }}\sigma _{\theta {\text{,1}}}^{\text{e}}{\text{ }}\sigma _{r{\text{,2}}}^{\text{e}}{\text{ }}\sigma _{z{\text{,2}}}^{\text{e}}{\text{ }}\tau _{rz{\text{,2}}}^{\text{e}}{\text{ }}\sigma _{\theta {\text{,2}}}^{\text{e}}\sigma _{r{\text{,3}}}^{\text{e}}{\text{ }}\sigma _{z{\text{,3}}}^{\text{e}} \\ {\tau }_{rz\text{,3}}^{\text{e}} {\sigma }_{\theta \text{,3}}^{\text{e}}{\}}_{1\times 12} \text{，}{\left\{ {b_{{\text{equil}}}^{\text{e}}} \right\}^{\text{T}}} = {\left\{ {0{\text{ }}\gamma } \right\}_{1 \times 2}}。$

式中： ${r_{ij}} = {r_i} - {r_j},{z_{ij}} = {z_i} - {z_j}$ ； $\bar r = ({r_1} + {r_2} + {r_3})/3$ 。

1.2 应力间断面约束条件

将相邻独立单元的共同面定义为应力间断面，如图 3所示。在应力允许的情况下，控制各节点的正应力和剪应力对应相等。通过一对共轭节点（1, 2）和（3, 4）定义了一个典型的应力间断线。间断线与r轴正向形成的角度ω，定义相对于r轴正向逆时针方向转动为正方向。

图 3 应力间断面

Figure 3. Stress discontinuity

下载: 全尺寸图片幻灯片

为了满足平衡条件，作用于该路径上的正应力（ $\sigma _{\text{n}}^{}$ ）和剪应力（ $\tau _{\text{n}}^{}$ ）的值应始终保持不变，即

${\sigma }_{\text{n,1}}^{\text{a}}={\sigma }_{\text{n,2}}^{\text{b}} \text{，}{\sigma }_{\text{n,3}}^{\text{a}}={\sigma }_{\text{n,4}}^{\text{b}}\text{，}{\tau }_{\text{n,1}}^{\text{a}}={\tau }_{\text{n,2}}^{\text{b}}\text{，}{\tau }_{\text{n,3}}^{\text{a}}={\tau }_{\text{n,4}}^{\text{b}} 。$

(3)

通过r-z平面某点的应力状态( $\sigma _r^{}$ ， $\sigma _z^{}$ 和 $\tau _{rz}^{}$ )，利用应力莫尔圆得到正应力和剪应力的值， $\sigma _{\text{n}}^{}$ 和 $\tau _{\text{n}}^{}$ 的最终表达式如下：

${\sigma _{\text{n}}} = {\sin ^2}\omega {\sigma _r} + {\cos ^2}\omega {\sigma _z} - \sin 2\omega {\tau _{rz}}\text{，}$

(4)

${\tau _{\text{n}}} = - \frac{1}{2}\sin 2\omega {\sigma _r} + \frac{1}{2}\sin 2\omega {\sigma _z} + \cos 2\omega {\tau _{rz}}。$

(5)

将式（3）的约束条件转化为离散形式：

${\left[ {A_{{\text{stat}}}^{{\text{dc}}}} \right]_{4 \times 16}}{\left\{ {{\sigma ^{{\text{dc}}}}} \right\}_{16 \times 1}} = {\left\{ {b_{{\text{stat}}}^{{\text{dc}}}} \right\}_{4 \times 1}} 。$

(6)

式中： ${\left[ {A_{{\text{stat}}}^{{\text{dc}}}} \right]_{{\text{4}} \times {\text{16}}}} = {\left[ {\begin{array}{*{20}{c}} T&{ - T} \\ 0&0 \end{array}{\text{ }}\begin{array}{*{20}{c}} {\text{0}}&{\text{0}} \\ T&{ - T} \end{array}} \right]_{{\text{4}} \times {\text{16}}}}$

${\left\{ {{\sigma ^{{\text{dc}}}}} \right\}^{\text{T}}}{\text{ = }}\left\{ {\sigma _{r{\text{,1}}}^{\text{a}}{\text{ }}\sigma _{z{\text{,1}}}^{\text{a}}{\text{ }}\tau _{rz{\text{,1}}}^{\text{a}}{\text{ }}\sigma _{\theta {\text{,1}}}^{\text{a}}{\text{ }}\sigma _{r{\text{,2}}}^{\text{b}}{\text{ }}\sigma _{z{\text{,2}}}^{\text{b}}{\text{ }}\tau _{rz{\text{,2}}}^{\text{b}}{\text{ }}\sigma _{\theta {\text{,2}}}^{\text{b}}} \right.$

${{\sigma }_{r\text{,3}}^{\text{a}}\text{ }{\sigma }_{z\text{,3}}^{\text{a}}\text{ }{\tau }_{rz\text{,3}}^{\text{a}}\text{ }{\sigma }_{\theta \text{,3}}^{\text{a}}\text{ }{\sigma }_{r\text{,4}}^{\text{b}}\text{ }{\sigma }_{z\text{,4}}^{\text{b}}\text{ }{\tau }_{rz\text{,4}}^{\text{b}}\text{ }{\sigma }_{\theta \text{,4}}^{\text{b}}\text{ }\}}_{\text{1}\times \text{16}} \text{，}$

${\left\{ {b_{{\text{stat}}}^{{\text{dc}}}} \right\}^{\text{T}}} = \left\{ {{\text{0 0 0 0}}} \right\}1 \times 4。$

1.3 土体边界条件

在任何土体边界l（如图 4所示）的边界条件可以用以下方式来定义：

${\sigma }_{\text{n,1}}^{\text{l}}={q}_{1}\text{，}{\sigma }_{\text{n,2}}^{\text{l}}={q}_{2}\text{，}{\tau }_{\text{n,1}}^{\text{l}}={t}_{1}\text{，}{\tau }_{\text{n,1}}^{\text{l}}={t}_{2} 。$

(7)

图 4 应力边界条件

Figure 4. Stress boundary conditions

下载: 全尺寸图片幻灯片

将任意土体边界l中两节点所对应公式（7）的4个等式约束转化为离散形式如下：

${\left[ {A_{{\text{bound}}}^{\text{l}}} \right]_{4 \times 8}}{\left\{ {{\sigma ^{\text{l}}}} \right\}_{8 \times 1}} = {\left\{ {b_{{\text{bound}}}^{\text{l}}} \right\}_{4 \times 1}}。$

(8)

式中： ${\left[{A}_{\text{bound}}^{\text{l}}\right]}_{4\times 8}={\left[\begin{array}{c}M 0\\ 0\text{ }M\end{array}\right]}_{4\times 8}$ ，

${\left[M\text{ }\right]}_{\text{2}\times \text{4}}=\left[\begin{array}{cc}{\mathrm{sin}}^{2}\alpha & {\mathrm{cos}}^{2}\alpha \\ -\frac{1}{2}\mathrm{sin}2\alpha & \frac{1}{2}\mathrm{sin}2\alpha \end{array}\text{ }\begin{array}{cc}-\mathrm{sin}2\alpha & 0\\ \mathrm{cos}2\alpha & 0\end{array}\right] \text{，}$

${\left\{{\sigma }^{\text{l}}\right\}}^{\text{T}}={\left\{{\sigma }_{r\text{,1}}^{\text{l}} {\sigma }_{z\text{,1}}^{\text{l}} {\tau }_{rz\text{,1}}^{\text{l}} {\sigma }_{\theta \text{,1}}^{\text{l}} {\sigma }_{r\text{,2}}^{\text{l}} {\sigma }_{z\text{,2}}^{\text{l}} {\tau }_{rz\text{,2}}^{\text{l}} {\sigma }_{\theta \text{,2}}^{\text{l}}\right\}}_{1\times 8} \text{，}$

${\left\{{b}_{\text{bound}}^{\text{l}}\right\}}^{T}={\left\{{q}_{1} {t}_{1} {q}_{2} {t}_{2} \right\}}_{1\times 4} 。$

式中， $\alpha$ 为边界l与r轴之间夹角，定义相对于r轴逆时针方向转动为正方向。

1.4 筒-土界面约束条件

在竖向荷载作用下，与筒壁内外侧接触的土体将产生沿筒壁向上的摩擦力，与筒顶及筒壁端部接触的土体将产生垂直端部向上的土抗力。因此，同时对筒顶-土接触面、筒端-土接触面及筒壁-土接触面的进行边界约束。

定义发生剪切破坏时沿筒顶-土界面、筒端-土界面和筒壁-土界面切向抗剪强度均满足以下条件：

$\left| {{\tau _s}} \right| \leqslant {\alpha _{\text{s}}}{s_{{\text{us}}}}\text{；}\left| {{\tau _{\text{b}}}} \right| \leqslant {\alpha _{\text{b}}}{s_{{\text{ub}}}}\text{；}\left| {{\tau _{{\text{f}}i}}} \right| \leqslant {\alpha _{{\text{f}}i}}{s_{{\text{u}}i}}。$

(9)

式中： ${\alpha _{\text{s}}}$ ， ${\alpha _{\text{b}}}$ 和 ${\alpha _{{\text{f}}i}}$ 分别代表筒顶-土界面、筒端-土界面和筒壁-土界面接触强度的折减系数； ${s_{{\text{us}}}}$ 和 ${s_{{\text{ub}}}}$ 分别为筒顶和筒端处的土体不排水抗剪强度， ${s_{{\text{u}}i}}$ 为筒壁-土界面上任一应力点所在深度处土体的不排水抗剪强度。

式（9）可离散转化为

${\left[ {{A_{{\text{int}}}}} \right]_{{\text{2}} \times {\text{4}}}}{\left\{ \sigma \right\}_{{\text{4}} \times {\text{1}}}} \leqslant {\left\{ {{b_{{\text{int}}}}} \right\}_{{\text{2}} \times {\text{1}}}}。$

(10)

式中： ${\left[ {{A_{{\text{int}}}}} \right]_{{\text{2}} \times {\text{4}}}}{\text{ = }}\left[ {\begin{array}{*{20}{c}} {\text{0}}&{\text{0}} \\ 0&0 \end{array}{\text{ }}\begin{array}{*{20}{c}} { - {\text{1}}}&0 \\ 1&0 \end{array}} \right]$ ，

${\left\{ \sigma \right\}^{\text{T}}}{\text{ = }}\left\{ {{\sigma _{r{\text{,}}i}}{\text{ }}{\sigma _{z{\text{,}}i}}{\text{ }}{\tau _{rz{\text{,}}i}}{\text{ }}{\sigma _{\theta {\text{,}}i}}} \right\}1 \times 4 \text{，}$

${\left\{ {{b_{{\text{int}}}}} \right\}^{\text{T}}}{\text{ = }}{\left\{ {{\alpha _{{\text{su}}}}{\text{ }}{\alpha _{{\text{su}}}}} \right\}_{{\text{1}} \times {\text{2}}}} 。$

图 5展示了筒顶-土界面（EP区域）、筒端-土界面（OQ区域）、筒壁-土界面（PQ区域）的剪应力约束形式，及各筒-土界面的剪应力正方向。

图 5 筒-土边界示意图

Figure 5. Boundary between bucket and soil

下载: 全尺寸图片幻灯片

1.5 土体屈服条件

根据Harr-VonKarman假设， ${\sigma _\theta }$ 为在r-z平面上的 ${\sigma _{\text{3}}}$ （最小主应力），即假定环向应力不大于大主应来求解轴对称稳定性问题^[13]。 ${\sigma _\theta }$ 的取值范围在是通过绘制两个莫尔圆确定，破坏面土体受力莫尔圆如所示。为了确保环向应力 ${\sigma _\theta }$ 的值始终接近于 ${\sigma _{\text{3}}}$ 的值，确定以下3个不等式约束：

图 6 破坏面土体受力莫尔圆

Figure 6. Mohr circle of failure surface

下载: 全尺寸图片幻灯片

$\left.\begin{array}{l}{\sigma }_{\theta \text{,}i}\ge {\sigma }_{r\text{,}i}\text{ }\text{，}\\ {\sigma }_{\theta \text{,}i}\ge {\sigma }_{z\text{,}i}\text{ }\text{，}\\ {\sigma }_{\theta \text{,}i}\le {\sigma }_{\text{3f,}i}\text{ }。\end{array}\right\}$

(11)

式中： ${\sigma _{\theta {\text{,}}i}}$ ， ${\sigma _{r{\text{,}}i}}$ 和 ${\sigma _{z{\text{,}}i}}$ 为与节点i相关的应力， ${\sigma _{{\text{3f}}}}$ 为屈服时的小主应力。这些约束可以用矩阵的形式表示为

${\left[ {A_{r\theta }^i} \right]_{{\text{3}} \times {\text{4}}}}{\left\{ {{\sigma ^i}} \right\}_{{\text{4}} \times {\text{1}}}} \leqslant {\left\{ {b_{r\theta }^i} \right\}_{{\text{3}} \times {\text{1}}}}。$

(12)

式中：

${\left[{A}_{r\theta }^{i}\right]}_{\text{3}\times \text{4}}=\left[\begin{array}{cc}1& 0\\ \begin{array}{l}\text{ }0\\ -{\scriptstyle \frac{1}{2}}\text{(}1-\mathrm{sin}\varphi \text{)}\end{array}& \begin{array}{l}\text{ 1}\\ \text{ }-{\scriptstyle \frac{1}{2}}(1-\mathrm{sin}\varphi )\end{array}\end{array}\text{ }\begin{array}{cc}\text{0}& -1\\ \begin{array}{l}0\\ \text{0}\end{array}& \begin{array}{l}-1\\ \text{ 1}\end{array}\end{array}\right]，$

${\left\{{\sigma }^{i}\right\}}^{\text{T}}={\left\{{\sigma }_{r\text{,}i}\text{ }{\sigma }_{z\text{,}i}\text{ }{\tau }_{rz\text{,}i}\text{ }{\sigma }_{\theta \text{,}i} \right\}}_{1\times 4} \text{，}$

${\left\{ {b_{r\theta }^i} \right\}^{\text{T}}} = {\left\{ {0{\text{ }}0{\text{ }}c\cos \varphi } \right\}_{1 \times 3}}。$

在下限解中用内接多边形来逼近Mohr-Coulomb准则，从而方便求解计算的线性规划问题。为保证任一点上的应力状态均满足屈服条件，需在r-z平面上保证满足以下线性化的不等式条件：

${A_k}{\sigma _r} + {B_k}{\sigma _z} + {C_k}{\tau _{rz}} \leqslant \varepsilon {\text{ (}}k = 1,2, \cdots ,p) 。$

(13)

式中： ${A_k} = \cos (2{\rm{\mathsf{π}}}k/p) + \sin \varphi \cos ({\rm{\mathsf{π}}}/p)$ ， ${B_k} = \sin \varphi \cdot$ $\mathrm{cos}({\rm{\mathsf{π}}}/p)-\mathrm{cos}(2{\rm{\mathsf{π}}}k/p)，{C}_{\text{k}}=2\mathrm{sin}(2{\rm{\mathsf{π}}}k/p)，\varepsilon =2c\mathrm{cos}\varphi \cdot$ $\cos ({\rm{\mathsf{π}}}/p)$ 。其中，c为节点i处土体的不排水抗剪强度。

由线性化屈服条件对节点i的4个应力张量施加的不等式约束（数目p）可表示为

${\left[ {A_{{\text{yield}}}^i} \right]_{p \times 4}}{\left\{ {{\sigma ^i}} \right\}_{4 \times 1}} \leqslant {\left\{ {b_{{\text{yield}}}^i} \right\}_{p \times 1}}。$

(14)

式中，

${\left[{A}_{\text{yield}}^{i}\right]}_{p\times 4}=\left[\begin{array}{c}{A}_{1} {B}_{1} {C}_{1} 0\\ {A}_{2} {B}_{2} {C}_{2} 0\\ .\text{ }.\text{ }.\\ {A}_{k} {B}_{k} {C}_{k} 0\\ .\text{ }.\text{ }.\\ {A}_{p} {B}_{p} {C}_{p} 0\end{array} \right] 。$

1.6 最优化函数

通过对筒型基础顶盖部分圆形区域、筒端部分圆环区域的垂直法向应力和筒壁区域切向应力积分。地基竖向承载力（目标函数）可通过以下表达式来定义：

$Q = {Q_{\text{s}}}{\text{ + }}{Q_{\text{f}}}{\text{ + }}{Q_{\text{b}}}{\text{ = }}2{\rm{\mathsf{π}}}\int\limits_{{L_{\text{s}}}} {{\sigma _{\text{s}}}} r{\text{d}}s{\text{ + }}2{\rm{\mathsf{π}}}\int\limits_{{L_{\text{f}}}} {{\tau _{\text{f}}}} z{\text{d}}s{\text{ + }}2{\rm{\mathsf{π}}}\int\limits_{{L_{\text{b}}}} {{\sigma _{\text{b}}}} r{\text{d}}s 。$

(15)

式中：Q为地基竖向承载力，Q_s，Q_f，Q_b分别为筒顶、筒壁及筒端承载部分。 ${\sigma _{\text{s}}}$ 和 ${\sigma _{\text{b}}}$ 分别为作用在筒顶-土和筒端-土边界面上的法向应力， ${\tau _{\text{f}}}$ 为作用在筒壁-土边界面上的切向抗剪强度。由于假设应力在每个单元中为线性变化，筒型基础各部分产生的极限承载力的分量分别为

${Q_{\text{s}}}{\text{ = }}2{\rm{\mathsf{π}}}\frac{{L_{\text{s}}^i}}{2}({\sigma _{{\text{s,1}}}}{\text{ + }}{\sigma _{{\text{s,2}}}})r_{\text{s}}^i \text{，}$

(16)

${Q_{\text{f}}}{\text{ = }}2{\rm{\mathsf{π}}}\frac{{L_{\text{f}}^i}}{2}({\tau _{{\text{f,1}}}}{\text{ + }}{\tau _{{\text{f,2}}}})z_{\text{f}}^i \text{，}$

(17)

${Q_{\text{b}}}{\text{ = }}2{\rm{\mathsf{π}}}\frac{{L_{\text{b}}^i}}{2}({\sigma _{{\text{b,1}}}}{\text{ + }}{\sigma _{{\text{b,2}}}})r_{\text{b}}^i 。$

(18)

式中： $L_{}^i$ 为节点i所在单元边界的长度； $r_{}^i$ 为边界 $L_{}^i$ 面的平均半径； $L_{\text{f}}^i$ 为筒壁节点上i所在单元边界的长度； $z_{\text{f}}^i$ 为边界 $L_{\text{f}}^i$ 面的平均长度。

将方程（1），（3）和（7）进行整合，得到全局等式约束矩阵形式如下：

$\left[ {{A_1}} \right] = \sum\limits_{{\text{e = }}1}^E {\left[ {A_{_{{\text{equil}}}}^{\text{e}}} \right]} {\text{ + }}\sum\limits_{{\text{dc}}}^{{D_{\text{c}}}} {\left[ {A_{_{{\text{stat}}}}^{{\text{dc}}}} \right]} {\text{ + }}\sum\limits_{{\text{l = }}1}^L {\left[ {A_{{\text{bound}}}^{\text{l}}} \right]} 。$

(19)

将方程（9），（10），（14）进行整合，得到全局不等式约束矩阵成如下形式：

$\left[ {{A_2}} \right] = \sum\limits_{i{\text{ = }}1}^N {\left[ {A_{_{{\text{yield}}}}^i} \right]} {\text{ + }}\sum\limits_{i{\text{ = }}1}^N {\left[ {A_{{\text{s }{\rm{\mathsf{θ}}}}}^i} \right]} {\text{ + }}\sum\limits_{{\text{ns = }}1}^{{N_{\text{s}}}} {\left[ {A_{{\text{ints}}}^{{\text{ns}}}} \right]} {\text{ + }}\sum\limits_{{\text{nf = }}1}^{{N_{\text{f}}}} {\left[ {A_{{\text{intf}}}^{{\text{nf}}}} \right]} +$

$\sum\limits_{{\text{nb = }}1}^{{N_{\text{b}}}} {\left[ {A_{{\text{intb}}}^{{\text{nb}}}} \right]} 。$

(20)

与等式及不等式约束条件对应的结果矩阵 ${b_1}$ 和 ${b_2}$ 形式如下：

$\left[ {{b_1}} \right] = \sum\limits_{{\text{e = }}1}^E {\left\{ {b_{_{{\text{equil}}}}^{\text{e}}} \right\}} {\text{ + }}\sum\limits_{{\text{dc}}}^{{D_{\text{c}}}} {\left\{ {b_{_{{\text{stat}}}}^{{\text{dc}}}} \right\}} {\text{ + }}\sum\limits_{{\text{l = }}1}^L {\left\{ {b_{{\text{bound}}}^{\text{l}}} \right\}} \text{，}$

(21)

$\left[ {{b_2}} \right] = \sum\limits_{i{\text{ = }}1}^N {\left\{ {b_{_{{\text{yield}}}}^i} \right\}} {\text{ + }}\sum\limits_{i{\text{ = }}1}^N {\left\{ {b_{{\text{s }{\rm{\mathsf{θ}}}}}^i} \right\}} {\text{ + }}\sum\limits_{{\text{ns = }}1}^{{N_{\text{s}}}} {\left\{ {b_{{\text{ints}}}^{{\text{ns}}}} \right\}{\text{ + }}\sum\limits_{{\text{nf = }}1}^{{N_{\text{f}}}} {\left\{ {b_{{\text{intf}}}^{{\text{nf}}}} \right\}} } \sum\limits_{{\text{nb = }}1}^{{N_{\text{b}}}} {\left\{ {b_{{\text{intb}}}^{{\text{nb}}}} \right\}} 。$

(22)

调用线性最优化函数对目标函数求解的方法，陈广思^[14]和Vishwas^[15]等已给出介绍，不在此赘述。

2. 离心试验验证

2.1 试验准备

为验证筒型基础地基承载力下限解的准确性，开展了饱和黏土地基离心模型试验。试验在天津大学TLJ-100A型土工离心机完成。模型箱内部尺寸为880 mm×595 mm×400 mm。将高岭土粉按照2倍液限的含水率进行真空搅拌，为模拟天然地基的自然固结过程，将土体在100g离心条件下固结10×12 h，固结完成土体厚度约为280 mm。

模型筒外径100 mm，裙板长度95 mm，壁厚10 mm，顶盖厚5 mm，筒顶处预留连接螺杆，与量程为8 kN的荷载传感器相连，模型示意图如图 7所示。

图 7 筒型基础模型示意图

Figure 7. Model for bucket foundation

下载: 全尺寸图片幻灯片

2.2 试验方案

试验前取样测得软黏土的有效重度 $\gamma '$ =4.83 kN/m³，采用十字板对模型槽内土体的不排水抗剪强度进行了测量，结果如图 8所示，模型土泥面处土体强度s_um=6 kPa且土体不排水抗剪强度以k=0.426 kPa/m随深度线性增加。

图 8 试验土体的不排水强度

Figure 8. Undrained shear strengths of soil

下载: 全尺寸图片幻灯片

试验布置方式如图 9所示，试验采用分级加载，具体加载方式见图 10。

图 9 筒型基础及传感器整体布置

Figure 9. Overall layout of bucket foundation and sensor

下载: 全尺寸图片幻灯片

图 10 荷载时间曲线

Figure 10. Load-time curves

下载: 全尺寸图片幻灯片

2.3 试验结果对比

离心机模型试验得到的p-s曲线如图 11所示，在整个沉降过程中未出现明显的拐点。利用初始直线段与末尾直线段的切线交点对应的q_u值作为地基极限承载力，如图中黑色虚线所示，此时对应的竖向承载力q_u约为93 kPa。

图 11 离心试验p-s曲线

Figure 11. Load-displacement curves

下载: 全尺寸图片幻灯片

根据试验参数建立有限元极限分析下限解分析模型，将下限解结果和试验结果进行对比见表 1。下限解结果与模型试验的结果吻合良好，下限解结果小于API规范结果，且计算结果误差均在可接受范围，验证了极限分析下限解的正确性。

表 1 结果对比

Table 1. Comparison of results

试验结果/kPa	API/kPa	下限解/kPa
93	104.84	92.02

下载: 导出CSV

| 显示表格

3. 筒型基础承载特性分析

3.1 下限分析模型的建立

以海上风电筒型基础作为主要研究对象，基础外径D=10 m，侧壁壁厚t_s=0.01 m。研究不同筒裙长度d筒型基础的承载特性，主要针对长径比d/D=0.25，0.5，0.75，1，1.5和2的基础进行分析。为降低边界效应，对土体模型的尺寸进行灵敏度分析，确定土体直径取D_soil=2.5D，土体高度d_soil=3D+d，网格加密形式及模型示意图见图 12。

图 12 筒型基础下限分析网格示意图

Figure 12. Meshes for lower bound analysis of bucke foundation

下载: 全尺寸图片幻灯片

利用基础横截面积S和筒端处土体不排水土强度s_u0对地基竖向承载力V_ult进行了无量纲化处理得到N_cV，即

$N_{\rm cV}=V_{\rm ult}/S·s_{\rm u0} 。$

(23)

3.2 土体强度变化对竖向承载力的影响

大量现场原位测试结果表明，对于大部分的正常固结土及超固结土而言，其不排水抗剪强度随埋深的增加而呈近似的线性增长趋势，假设黏土的不排水剪切强度随深度呈线性变化：s_uz=s_um+kz。其中k为土体强度随深度的增长系数，s_um为泥面处土体的不排水抗剪强度，s_uz为地基深度为z处土体的不排水抗剪强度。

为得出均质土不排水抗剪强度s_u对筒型基础地基竖向极限承载力的影响，针对s_u=10~40 kPa的均质黏土地基进行计算，分析模型中筒-土界面设置为完全粗糙接触，计算结果如图 13所示。

图 13 土体不排水强度对筒型基础竖向承载力的影响

Figure 13. Effects of s_u on vertical bearing capacity of bucket foundation

下载: 全尺寸图片幻灯片

由图 13可知，当d/D＜1时，N_cV增长曲线随土体强度增长近似水平，故认为在此范围内筒型基础承载力与s_u无关；而当d/D＞1时，N_cV随强度增长有较为轻微下降趋势。总体来看，对于长径比较大的筒型基础N_cV虽随s_u增长而略有降低，但结合现有筒型基础工程实际情况，此部分降低可忽略。

为得出非均质土体沿深度方向强度不均匀性对筒型基础竖向承载力的影响，针对s_um=4 kPa及s_um=40 kPa且k均为0~2的地基进行了下限解求解，计算结果如所示。引入不同土体强度线性增长梯度下筒型基础地基竖向极限承载力与均质土体（k=0）条件下承载力的比值 ${\eta _k}$ 来评估土体k增长对筒型基础地基极限承载力的增益。 ${\eta _k}$ 定义如下：

${\eta _k} = \frac{{{V_{{\text{ult}}}}(k)}}{{{V_{{\text{ult}}}}(k = 0)}}。$

(24)

由图 14可知，当s_um=4时，随着土体强度增长梯度k由0增长至2， ${\eta _k}$ 增长约120%~410%；当s_um=40时，随着土体强度增长梯度k由0增长至2， ${\eta _k}$ 增长约2%~35%，不同长径比筒型基础V_ult皆随k增长而不同程度的线性增加。单独考虑土体参量s_um或k时，无法通过N_cV来表征竖向承载力特性。为使承载力系数具有普适性，通过无量纲参数κ来定量表征土体强度沿深度方向不均匀性，定义κ=kD/s_um，计算时k和D固定不变，通过改变s_um来改变κ。针对κ=0.5，2，3，5，10，20的非均质饱和黏土地基进行讨论，计算结果如图 15所示。

图 14 k对筒型基础竖向承载力影响

Figure 14. Effects of k on vertical bearing capacity of bucket foundation

下载: 全尺寸图片幻灯片

图 15 κ对筒型基础竖向承载力的影响

Figure 15. Effects of κ on vertical bearing capacity of bucket foundation

下载: 全尺寸图片幻灯片

如图 15所示，N_cV随κ的增加而减小。因为κ越高意味沿裙板上平均强度越低，筒壁侧摩阻发挥较低，所以其承载力增长率比基础端部土体强度的增长率低，从而使N_cV降低。

3.3 长径比对筒型基础竖向极限承载力的影响

对均质土中长径比为0.5~2的筒型基础地基承载力进行计算，得到其地基承载力随长径比增长的变化，结果如图 16所示。

图 16 不同长径比的筒型基础在均质土中的承载特性

Figure 16. Vertical bearing capacities of bucket foundation with different length diameter ratios in homogeneous soils

下载: 全尺寸图片幻灯片

由图 16可知，均质土中筒型基础承载力系数N_cV增率随d/D增大呈现非线性增长趋势，且规律随土体强度变化并未发生改变。将均质土中筒型基础地基竖向承载力系数N_cV随土体不排水抗剪强度s_u及长径比d/D的变化规律拟合为

${N_{{\text{cV}}}}{\text{(}}\kappa {\text{ = 0) = 5}}{\text{.91 + 8}}{\text{.41}}\left( {\frac{d}{D}} \right) - {\text{1}}{\text{.12}}{\left( {\frac{d}{D}} \right)^2} 。$

(25)

对于强度沿深度线性增加黏土地基，选取κ=0.5，2，3，5，10，20对不同长径比筒型基础地基竖向承载力进行分析，得到N_cV随长径比增长的变化规律，结果如图 17所示。

图 17 不同长径比筒型基础在非均质土中的承载特性

Figure 17. Vertical bearing capacities of bucket foundation with different length diameter ratios in nonhomogenous soil

下载: 全尺寸图片幻灯片

由图 17可知，对于强度沿深度线性增加的土体，N_cV随d/D的增长而增大。将筒型基础地基竖向承载力系数N_cV随土体强度沿深度不均匀性κ及长径比d/D的变化关系可采用公式（26）进行拟合，由图 17可知，拟合结果良好。

${N_{{\text{cV}}}}{\text{(}}\kappa {\text{) = 3}}{\text{.22 + 2}}{\text{.02}}\left( {\frac{d}{D}} \right){\text{ + 6}}{\text{.05}}{{\text{e}}^{ - 0.45\kappa }} 。$

(26)

为探究筒型基础不同部分对地基承载力的贡献率随长径比的变化，分别计算不同土体强度（s_um=10 kPa，s_um=40 kPa）及不同土体强度线性增长梯度k（k=0，2）条件下不同长径比的筒型基础地基竖向承载力，提取极限状态下筒壁内、外侧摩阻力，筒顶抗力及筒端处抗力。考虑到荷载情况仅为竖向荷载，筒内壁与土体接触良好，认为土体抗力为筒顶抗力和筒内壁摩阻之和，分析各部分在地基竖向承载力中占比，结果绘制于图 18。

图 18 筒基各部分承担荷载比例随长径比的变化

Figure 18. Variation of loading ratio in different parts of bucket foundations with d/D

下载: 全尺寸图片幻灯片

由图 18可知，在不同土体强度条件下，随着基础长径比的增加，筒端抗力贡献率稳定在10%以下，土体抗力的贡献率逐渐降低，筒壁外侧摩阻逐渐增大。分析原因认为，随d/D的增加，筒土接触面积逐步增大，外侧摩阻增长显著，筒端处虽应力较高，但由于筒壁较薄，故对承载力贡献率不高。在相同长径比条件下，随κ增大，土体抗力的贡献率逐步增加。由于κ越高意味着筒端下方参与承载的土体强度越高，筒端以上和筒壁接触的土体强度相对较弱，无法有效地为筒壁-土接触面提供剪切力，故而筒壁外侧摩阻的增长相对较低。

3.4 ${\alpha _{{\text{f}}i}}$ 沿筒壁分布规律

为探究土体达到极限承载状态时筒外壁-土界面接触强度折减程度沿深度方向的实际分布规律，以均质土（s_um=10 kPa，k=0，即κ=0）和强度随深度线性增长黏土（s_um=10 kPa，k=2，即κ=2）地基为例，计算筒外壁-土间完全粗糙接触条件下的地基竖向承载力。提取极限竖向荷载状态下筒外壁-土界面切向抗剪强度 ${\tau _{{\text{f}}i}}$ 及应力点处不排水抗剪强度s_ui，参考式（9）进行归一化处理，即 ${\alpha _{{\text{f}}i}}$ = ${\tau _{{\text{f}}i}}$ /s_ui，计算结果如图 19所示。

图 19

${\alpha _{{\text{f}}i}}$ 沿深度分布规律

Figure 19. Distribution laws of

${\alpha _{{\text{f}}i}}$ along depth

下载: 全尺寸图片幻灯片

图 19中纵坐标w/D为应力点距筒顶的距离w与筒径D之比。分析可知，极限状态时，尽管筒外壁-土间接触关系设定为完全粗糙， ${\alpha _{{\text{f}}i}}$ 沿深度分布也并非定值。

均质土中，筒壁上部w/D=0~0.5范围内，d/D＜1的筒型基础 ${\alpha _{{\text{f}}i}}$ 随深度增加而增大，d/D＞1的筒型基础筒壁侧摩阻力不变且 ${\alpha _{{\text{f}}i}}$ 近似等于1。而当w/D＞0.5，随着应力点靠近筒端 ${\alpha _{{\text{f}}i}}$ 逐渐减小。

对于强度沿深度线性增长的黏土地基， ${\alpha _{{\text{f}}i}}$ 随深度增长而降低，w/D=0~0.5范围内降低速率较小，而距筒顶0.5D后 ${\alpha _{{\text{f}}i}}$ 开始加速降低。对于同一深度的筒壁上各应力点，强度线性增长黏土中 ${\alpha _{{\text{f}}i}}$ 均小于均质土，因此，对于黏土地基，筒壁-土界面处土体不排水抗剪强度越高， ${\alpha _{{\text{f}}i}}$ 越小。

3.5 ${\alpha _{{\text{f}}i}}$ 对筒型基础地基竖向承载力的影响

以均质土（κ=0）和强度沿深度线性增长的黏土地基（κ=20）为例，计算不同 ${\alpha _{{\text{f}}i}}$ 的筒型基础地基竖向承载力，结果绘制于图 20。

图 20 α_fi对筒型基础竖向承载力的影响

Figure 20. Effects of α_fi on vertical bearing capacity of bucket foundation

下载: 全尺寸图片幻灯片

由图 20可知，对于均质黏土和强度随深度线性增长的黏土地基，N_cV均随 ${\alpha _{{\text{f}}i}}$ 的增大而增大，并且长径比大的筒型基础N_cV增率更为显著。其原因在于长径比大的筒型基础筒外壁-土接触面积更大，筒外壁侧摩阻对承载力的贡献也更大。

综合考虑土体不均匀性κ、筒型基础长径比d/D和筒外壁-土界面接触强度折减系数 ${\alpha _{{\text{f}}i}}$ 的影响，均质黏土及强度沿深度线性增长的黏土地基中筒型基础地基竖向承载力下限解可拟合为

${N}_{\text{cV}}{}_{（\kappa =0）}\text{=}\left(0.72+0.35{\alpha }_{\text{f}}-0.07\frac{d}{D}\right)\cdot$

$\left[ {5.91 + 8.41\frac{d}{D} - 1.12{{\frac{d}{D}}^2}} \right] \text{，}_{ }$

(27)

${N}_{\text{cV}}（\kappa ）=\left(0.79+0.29{\alpha }_{\text{f}}-0.07\frac{d}{D}\right)\cdot$

$\left[ {3.22 + 2.02\frac{d}{D} + 6.05{{\text{e}}^{ - 0.45\kappa }}} \right] \text{，}$

(28)

$V_{{\rm ult}（κ=0）}=N_{{\rm cV}（κ=0）}·As_{\rm u} \text{，}$

(29)

$V_{{\rm ult}（κ）}=N_{{\rm cV}（κ）}·As_{\rm u0} 。$

(30)

4. 结论

本文针对黏性土地基上筒型基础地基竖向承载力问题，采用轴对称有限元极限分析下限法进行分析，主要得到以下4点结论。

（1）将筒型基础简化为轴对称模型并考虑筒-土间摩擦，求得了筒型基础地基竖向承载力下限解，并对计算结果进行了离心模型试验验证。

（2）在饱和黏土地基中，筒型基础地基竖向承载力系数随土体强度不均匀系数的增大而减小，随筒型基础长径比增长而增长，并给出了拟合公式。

（3）土体抗力和侧摩阻力是筒型基础主要承力部位，随长径比的增加或土体强度不均匀性的降低，土体抗力对承载力贡献率逐渐下降，侧壁摩阻贡献逐渐增长。

（4）对于长径比小于1的筒型基础，筒外壁-土间侧摩阻力随深度呈现先增大后减小趋势；对于长径比大于1的筒型基础筒顶下0.5D范围内外壁侧摩阻发挥较高，而后随深度增加而逐渐降低至近筒端处趋于0。

图 1 原状土参数测定过程

Figure 1. Determination process of parameters of undisturbed soil

下载: 全尺寸图片幻灯片

图 2 相似土样制备

Figure 2. Preparation of similar soil samples

下载: 全尺寸图片幻灯片

图 3 不同配比土样的强度参数

Figure 3. Strength parameters of soil samples with different proportions

下载: 全尺寸图片幻灯片

图 4 土压力传感器布置

Figure 4. Arrangement of earth pressure sensor

下载: 全尺寸图片幻灯片

图 5 地表、地层位移监测装置

Figure 5. Monitoring devices for surface and stratum displacements

下载: 全尺寸图片幻灯片

图 6 不同含水率观察面情况（3D）

Figure 6. Observation surfaces under different water contents (3D)

下载: 全尺寸图片幻灯片

图 7 不同埋深观察面情况（26.6%）

Figure 7. Situations of observation surface under different buried depths (26.6%)

下载: 全尺寸图片幻灯片

图 8 不同含水率地层变形（3D）

Figure 8. Deformations of strata with different water contents (3D)

下载: 全尺寸图片幻灯片

图 9 不同埋深地层变形

Figure 9. Deformations of strata at different depths

下载: 全尺寸图片幻灯片

图 10 不同含水率地表沉降（3D）

Figure 10. Surface subsidences under different water contents (3D)

下载: 全尺寸图片幻灯片

图 11 不同埋深地表沉降（22.3%）

Figure 11. Surface subsidences at different buried depths (22.3%)

下载: 全尺寸图片幻灯片

图 12 不同埋深地表沉降（26.6%）

Figure 12. Surface subsidences at different buried depths (26.6%)

下载: 全尺寸图片幻灯片

图 13 不同含水率应力释放

Figure 13. Stress release under different water contents

下载: 全尺寸图片幻灯片

图 14 不同埋深应力释放（22.3%）

Figure 14. Stress release at different depths (22.3 %)

下载: 全尺寸图片幻灯片

图 15 不同埋深应力释放（26.6%）

Figure 15. Stress release at different depths (26.6 %)

下载: 全尺寸图片幻灯片

图 16 不同含水率围岩应力与释放率（3D）

Figure 16. Stresses and release rates of surrounding rock under different water contents (3D)

下载: 全尺寸图片幻灯片

图 17 不同埋深围岩应力与释放率

Figure 17. Stresses and release rates of surrounding rock under different buried depths

下载: 全尺寸图片幻灯片

表 1 原状土参数

Table 1 Parameters of undisturbed soil

含水率 w/%	天然密度 ρ₀/(g·cm^-3)	干密度 ρ_d/(g·cm^-3)	土粒相对质量密度G_s	黏聚力 c/kPa	内摩擦角 φ/(°)
5.9	1.44	1.36	2.71	53.86	26.73
7.8	1.47	1.36	2.71	50.99	25.88
11.8	1.52	1.36	2.71	45.13	23.89
15.1	1.56	1.36	2.71	38.77	22.46
16.8	1.59	1.36	2.71	32.00	20.94
20.2	1.63	1.36	2.71	30.12	19.76
25.4	1.70	1.36	2.71	29.01	18.35

下载: 导出CSV

表 2 试验工况

Table 2 Test conditions

组号	含水率/%	埋深
1	8.2	3 D
2	16.7	3 D
3	22.3	1 D
4	22.3	2 D
5	22.3	3 D
6	26.6	1 D
7	26.6	2 D
8	26.6	3 D

下载: 导出CSV

参考文献(22)

[1]	王明年, 郭军, 罗禄森, 等. 高速铁路大断面黄土隧道深浅埋分界深度研究[J]. 岩土力学, 2010, 31(4): 1157-1162. doi: 10.3969/j.issn.1000-7598.2010.04.027 WANG Mingnian, GUO Jun, LUO Lusen, et al. Study of critical buried depth of large cross-section loess tunnel for high speed railway[J]. Rock and Soil Mechanics, 2010, 31(4): 1157-1162. (in Chinese) doi: 10.3969/j.issn.1000-7598.2010.04.027
[2]	朱才辉, 贺豪楠, 兰开江, 等. 黄土地层含水率增大对大跨度隧道围岩压力影响[J]. 岩土工程学报, 2021, 43(增刊1): 93-98. doi: 10.11779/CJGE2021S1017 ZHU Caihui, HE Haonan, LAN Kaijiang, et al. Influences of increase of moisture content on surrounding soil pressure of large-span tunnels in loess[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(S1): 93-98. (in Chinese) doi: 10.11779/CJGE2021S1017
[3]	陈福江. 黄土隧道围岩含水量变化对隧道形态影响的研究[D]. 成都: 西南交通大学, 2008. CHEN Fujiang. Study on the Tunnel's Modality Impact of Changes of Surrounding Rock Water Content at Loess Tunnel[D]. Chengdu: Southwest Jiaotong University, 2008. (in Chinese)
[4]	SHI W, QIU J L, ZHANG C P, et al. Immersion mode and spatiotemporal distribution characteristic of water migration in loess tunnel[J]. Arabian Journal of Geosciences, 2022, 15(7): 654. doi: 10.1007/s12517-022-09565-7
[5]	苗贺朝, 邵生俊, 李骏, 等. 隧道穿越非饱和黄土场地浸水试验初探[J]. 地下空间与工程学报, 2019, 15(4): 1066-1072. https://www.cnki.com.cn/Article/CJFDTOTAL-BASE201904013.htm MIAO Hechao, SHAO Shengjun, LI Jun, et al. Primary investigation of infiltration field tests on unsaturated loess tunnel[J]. Chinese Journal of Underground Space and Engineering, 2019, 15(4): 1066-1072. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BASE201904013.htm
[6]	邵生俊, 李骏, 李国良, 等. 大厚度湿陷性黄土隧道现场浸水试验研究[J]. 岩土工程学报, 2018, 40(8): 1395-1404. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201808005.htm SHAO Shengjun, LI Jun, LI Guoliang, et al. Field immersion tests on tunnel in large-thickness collapsible loess[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(8): 1395-1404. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201808005.htm
[7]	邵生俊, 李骏, 邵将, 等. 大厚度湿陷性黄土地层的现场砂井浸水试验研究[J]. 岩土工程学报, 2016, 38(9): 1549-1558. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201609002.htm SHAO Shengjun, LI Jun, SHAO Jiang, et al. In-situ sand well immersion tests on self-weight collapsible loess site with large depth[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(9): 1549-1558. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201609002.htm
[8]	苏忍, 张恒睿, 张稳军, 等. 兰州地铁大厚度湿陷性黄土地层的现场浸水试验研究[J]. 土木工程学报, 2020, 53(增刊1): 186-193. https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC2020S1030.htm SU Ren, ZHANG Hengrui, ZHANG Wenjun, et al. Immersion tests on self-weight collapsible loess site with large depth of Lanzhou metro line[J]. China Civil Engineering Journal, 2020, 53(S1): 186-193. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC2020S1030.htm
[9]	李骏, 邵生俊, 李国良, 等. 黄土隧道的湿陷变形规律及其对衬砌结构的作用[J]. 岩石力学与工程学报, 2018, 37(1): 251-260. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201801025.htm LI Jun, SHAO Shengjun, LI Guoliang, et al. Collapse deformation of loess tunnel and its effect[J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(1): 251-260. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201801025.htm
[10]	LI J, SHAO S J, SHAO S. Collapsible characteristics of loess tunnel site and their effects on tunnel structure[J]. Tunnelling and Underground Space Technology, 2019, 83: 509-519. doi: 10.1016/j.tust.2018.08.035
[11]	翁效林, 王俊, 王立新, 等. 黄土地层浸水湿陷对地铁隧道影响试验研究[J]. 岩土工程学报, 2016, 38(8): 1374-1380. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201608003.htm WENG Xiaolin, WANG Jun, WANG Lixin, et al. Experimental research on influence of loess collapsibility on subway tunnels[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(8): 1374-1380. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201608003.htm
[12]	WENG X L, ZHOU R M, RAO W, et al. Research on subway shield tunnel induced by local water immersion of collapsible loess[J]. Natural Hazards, 2021, 108(1): 1197-1219. doi: 10.1007/s11069-021-04727-4
[13]	WENG X L, SUN Y F, ZHANG Y W, et al. Physical modeling of wetting-induced collapse of shield tunneling in loess strata[J]. Tunnelling and Underground Space Technology, 2019, 90: 208-219. doi: 10.1016/j.tust.2019.05.004
[14]	张玉伟, 宋战平, 翁效林, 等. 大厚度黄土地层浸水湿陷对地铁隧道影响的模型试验研究[J]. 岩石力学与工程学报, 2019, 38(5): 1030-1040. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201905017.htm ZHANG Yuwei, SONG Zhanping, WENG Xiaolin, et al. Model test study on influence of the collapsibility of large thickness loess stratum on subway tunnels[J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(5): 1030-1040. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201905017.htm
[15]	王二磊, 梁庆国, 王飞, 等. 不同浸水方式对黄土地铁隧道变形影响研究[J]. 铁道科学与工程学报, 2018, 15(1): 156-162. doi: 10.3969/j.issn.1672-7029.2018.01.022 WANG Erlei, LIANG Qingguo, WANG Fei, et al. Study on the influence of different flooding methods on the deformation of loess metro tunnel[J]. Journal of Railway Science and Engineering, 2018, 15(1): 156-162. (in Chinese) doi: 10.3969/j.issn.1672-7029.2018.01.022
[16]	饶伟. 湿陷性黄土地层局部浸水对地铁盾构隧道的影响研究[D]. 西安: 长安大学, 2019. RAO Wei. Research on Effect of Local Immersion of Collapsible Loess Stratum on Metro Shield Tunnel[D]. Xi'an: Changan University, 2019. (in Chinese)
[17]	董岳林. 深厚黄土地层湿陷对地铁盾构隧道结构变形的影响研究[D]. 西安: 长安大学, 2018. DONG Yuelin. Effect of Tunnel Deformation on Metro Shield Tunnel Structure in Deep Loess Areas[D]. Xi'an: Changan University, 2018. (in Chinese)
[18]	黄文鹏. 黄土地层浸水湿陷对盾构隧道结构影响机制及处治技术研究[D]. 西安: 长安大学, 2021. HUANG Wenpeng. Research on Mechanism and Treatment Technology of Shield Tunnel Structures Affected by Submersion and Collapse of Loess Layer[D]. Xi'an: Changan University, 2021. (in Chinese)
[19]	王利明. 黄土湿陷对地铁结构稳定性影响机制及剩余湿陷量控制标准研究[D]. 西安: 长安大学, 2020. WANG Liming. Study on the Influence Mechanism of Loess Collapsing on Metro Structure Stability and the Control Standard of Remaining Collapsing Amount[D]. Xi'an: Changan University, 2020. (in Chinese)
[20]	QIU J L, LU Y Q, LAI J X, et al. Experimental study on the effect of water gushing on loess metro tunnel[J]. Environmental Earth Sciences, 2020, 79(11): 1-19.
[21]	QIU J L, LU Y Q, LAI J X, et al. Failure behavior investigation of loess metro tunnel under local-high-pressure water environment[J]. Engineering Failure Analysis, 2020, 115: 104631. doi: 10.1016/j.engfailanal.2020.104631
[22]	扈世民, 张顶立, 王梦恕. 大断面黄土隧道开挖引起的围岩力学响应[J]. 中国铁道科学, 2011, 32(5): 50-55. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGTK201105010.htm HU Shimin, ZHANG Dingli, WANG Mengshu. The mechanical response of the surrounding rock caused by the excavation of loess tunnel with large cross section[J]. China Railway Science, 2011, 32(5): 50-55. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZGTK201105010.htm