Centrifugal model tests on deterioration process of canal under cyclic action of coupling wetting-drying and freeze-thaw
-
摘要: 季冻土地区复杂的水-热耦合作用是输水渠道等水工结构物发生劣化失稳的诱因。研制了一套季冻区渠道湿干冻融离心模型试验设备,实现了在离心场下进行湿润—干燥—冻结—融化过程的水-热耦合循环作用模拟。以北疆输水渠道为原型,开展了渠道在湿干冻融耦合循环作用下渠道劣化过程的离心模型试验。试验结果表明：湿干冻融耦合循环作用诱发模型渠道中的裂隙不断拓展,并最终产生一条贯穿渠顶表面的深层横向张拉裂隙,导致渠道有自渠顶张拉裂隙发生劣化失稳的趋势,这一过程中,渠水入渗深度不断增大,入渗量逐渐减小,土体的冻结速率和融化速率呈上升趋势。试验观测出渠坡处产生的累积法向变形最大,渠顶次之,渠底最小。所采用的离心模拟试验设备及测试方法对于寒冷地区渠道灾变问题的研究是适用的。Abstract: The complex action of coupling moisture and heat in seasonally frozen soil regions is the inducement of deterioration and instability of hydraulic structures such as canals. A set of centrifugal model test equipment for canal under the cyclic action of coupling wetting-drying and freeze-thaw in seasonally frozen soil regions is developed, which realizes the simulation of cyclic action of coupling moisture and heat during the process of wetting-drying and freeze-thaw in the centrifugal field. The centrifugal model tests are conducted to investigate the deterioration process of the canal in north Xinjiang under the cyclic action of coupling wetting-drying and freeze-thaw. The tests results show that the cyclic action of coupling wetting-drying and freeze-thaw makes the cracks of model canal expand, and eventually produces a deep transverse tension crack across the top surface of the canal, which leads to that the model canal has a tendency of instability from the tension cracks at the canal top. In this process, the infiltration depth of canal water increases and the infiltration volume decreases. The freezing rate and thawing rate show an upward trend. In addition, the final normal deformation at the canal slope is the largest, followed by that at the top of the canal and that at the bottom of the canal. The equipment and method adopted in the centrifugal model tests are applicable to the disaster researches on canals in cold regions.
-
Keywords:
- wetting-drying and freeze-thaw /
- canal /
- crack /
- centrifugal model test
-
0. 引言
盾构法隧道是目前软土地层中最主要的地铁隧道结构之一。地铁隧道在逐步网络状的同时其结构会不可避免受到城市建设施工影响，如邻近大面积突发堆载、新建线路近距离穿越、紧邻深基坑开挖等。盾构隧道作为复合细长结构在外部施工作用下将诱发其纵向不均匀位移，当引起的变形超过其自身结构的承载极限时，将会导致地铁隧道发生环间变形、管片渗水和破裂甚至列车脱轨等严重安全事故。因此研究地铁盾构隧道在外荷载作用下的纵向变形对地铁隧道的保护具有重要的意义。
盾构隧道是由一系列预制钢筋混凝土管片通过环间螺栓连接而成，其最明显的特征是存在大量薄弱的环间接头。这决定了在遭受外部荷载时，其纵向变形特性和承载能力与连续细长结构物（如桩基和连续管道）有较大的差异。在外部施工作用时，拼装为整体的管环段由于刚度较大，不易发生剪切和弯曲变形，往往以刚体移动为主，而环间接头由于抵抗变形能力较弱，则容易发生转动和错台[1]。这使得盾构隧道在承受外荷载时其纵向位移分布具有明显的不连续性，环间同时存在张开和错台变形，如图 1所示。
![]()
图 1 盾构隧道在外荷载下的纵向变形Figure 1. Longitudinal deformations of a shield tunnel subjected to.external loads目前分析和预测新建隧道穿越[2-5]，临近基坑开挖[6-7]和地表堆载等外部施工作用下盾构隧道的纵向变形，主要采用的是等效连续地基梁模型，如图 2（a）所示。等效连续地基梁模型是将盾构隧道视为地基上具有等效刚度的连续长梁。目前，最为常用的等效连续梁模型为Euler-bernoulli(EB)梁和Timoshenko(TM)梁模型。等效连续梁模型因其概念简单、可靠有效而在临近施工引起盾构隧道的纵向变形评估中备受青睐。Liang等[2-3]先后将既有盾构隧道简化为EB梁和TM梁模型，推导得到了新建隧道上穿时既有盾构隧道纵向解答；Zhang等[4]、Yu等[5]分别将既有盾构隧道视为EB和TM梁来预测新建隧道下穿时既有隧道的纵向行为；陈仁朋等[6]、Cheng等[7]将TM梁模型用于分析临近基坑开挖引起的既有盾构隧道纵向变形。
等效连续梁模型在分析盾构隧道纵向变形特点时，往往会忽略隧道环段和环间接头的刚度差异，且难以真实反映盾构隧道的环间张开与错台。为此部分学者提出了协同变形模型[8]，如图 2（b）所示。在该模型中，盾构隧道的管环段被视为刚体，环间接头通过剪切和拉伸弹簧来模拟，Liu等[9]进一步在接头增加了压缩弹簧以考虑轴向作用。协同变形模型能够有效地预测临近施工引起的隧道纵向变形。Liu等[9-10]应用协同变形模型分别提出新建隧道上穿和下穿作用下既有盾构隧道的纵向变形预测理论方法；魏纲等[11]、Huang等[12]分别基于协同变形模型推导得到了临近基坑开挖卸载作用下既有盾构隧道的纵向变形解答，随后张治国等[13]进一步采用该模型用于研究新建盾构下穿下考虑衬砌渗透性的既有隧道纵向变形。
相较于以上两种盾构隧道模型，Koizumi等[14]提出的纵向梁-弹簧模型能更真实地反映盾构隧道的环间变形，该模型将盾构隧道视为由一系列短梁和具有不同功能的接头弹簧连接而成，如图 2（c）所示。但由于目前尚未有简便、有效的解决方案使得该模型难以普遍采用。Huang等[15]首先基于状态空间法给出了Winkler地基上纵向梁-弹簧模型的简便解，但其方法并不能同时描述隧道环段和接头的变形，得到的隧道位移曲线是连续的。
针对前人研究的不足，采用纵向梁-弹簧模型构建盾构隧道纵向变形模型，引入Timoshenko梁表征隧道环段的变形，结合剪切弹簧和转动弹簧模拟环间接头受力特点。通过构建位于弹性地基上纵向梁-弹簧模型的有限差分方程，以解决环间接头-管环非连续变形求解问题，从而推导既有隧道在外荷载下的纵向变形公式。在此基础上，分别建立新建隧道上穿和下穿施工引起的外荷载作用下既有盾构隧道的纵向变形计算方法，并与新建隧道上、下穿越施工实测案例及既有理论方法对比以验证所提方法的合理性。
1. 简化的纵向梁-弹簧盾构隧道模型
为同时考虑盾构隧道环间张开和错台变形以及环间接头刚度的弱化，引入纵向梁-弹簧模型，如图 3所示。本文对该模型基本假设如下：①盾构隧道由管环段通过环间接头连接而成；②隧道环段变形采用Timoshenko梁模拟；③环间接头采用线性转动和剪切弹簧模拟。
基于单变量的Timoshenko梁理论[16]，可得管环段的位移、内力和几何方程：
w=wb+ws,w=wb−EIκGAd2wbdy2,ws=−EIκGAd2wbdy2； (1a) θ=dwbdy,γ=−EIκGAd3wbdy3； (1b) M=−EId2wbdy2,Q=−EId3wbdy3,EId4wbdy4=q。 (1c) 式中：θ为环段截面转角；w为环段竖向位移；ws和wb分别为环段的剪切和弯曲位移；M和Q分别为环段弯矩和剪力；EI为环段抗弯刚度；GA为环段抗剪刚度；κ为Timoshenko梁剪切系数，圆形隧道取0.53[17]；q为作用于环段上的外部荷载。
对于隧道接头处弯矩和剪力有以下表达式[1]：
M−j=M+j=Mj=kθΔθ=kθ(θ−−θ+)， (2a) Q−j=Q+j=Qj=ksδ=ks(w−b−w+b)。 (2b) 式中：M−j，M+j和Mj分别为隧道接头左侧隧道截面、右侧隧道截面以及接头处的弯矩；Q−j，Q+j和Qj分别为隧道接头左侧隧道截面、右侧隧道截面以及接头处的剪力；kθ，ks分别为隧道环间接头转动刚度和剪切刚度；Δθ，δ分别为接头段两侧隧道截面的相对转角和剪切错台；θ−，θ+分别为接头段左侧和右侧隧道截面转角；w−b，w+b分别为接头段左侧和右侧隧道截面的弯曲位移。
2. 外荷载下盾构隧道纵向变形
2.1 盾构隧道纵向变形推导
图 4为在外荷载下盾构隧道纵向响应计算模型。将盾构隧道视为Winkler地基上的简化纵向梁-弹簧模型。根据管环段受力平衡及弯矩平衡，可得到隧道环段在纵向外荷载作用下弯曲位移wb(y)的微分控制方程为
d4wbdy4−kDtκGAd2wbdy2+kDtEIwb=DtEIq(y)。 (3) ![]()
图 4 在外荷载下盾构隧道纵向响应计算模型Figure 4. Model for longitudinal responses of shield tunnel to external loads式中：Dt为既有隧道直径；k为地基反力系数。
由图 4可见，盾构隧道包含管环段和大量的环间接头，其受力变形具有不连续的特点，因此，可通过对其离散化处理以得到其数值解答。图 5为盾构隧道纵向离散化示意图。假定盾构隧道含有n2个环段，每个环段离散为n1个长度为l的节点。为考虑边界条件影响，在隧道两端各增加2个虚拟节点。因此，盾构隧道共可分为n1n2+4个节点。
根据有限差分原理，式（3）在第i个节点上可改写为以下差分表达式：
wb(i+2)−4wb(i+1)+6wb(i)−4wb(i−1)+wb(i−2)l4−kDtκGA.wb(i+1)−2wb(i)+wb(i−1)l2+kDtEIwbi=DtEIqi。 (4) 由于挠度wb(y)仅在隧道环段上连续可导，而在环间接头处既不连续也不可导，因此需在接头处附近增设虚拟节点，如图 6所示，虚拟节点表示若环间接头不存在时隧道变形后的位置。
根据式（1）和（2），第j个接头处的弯矩Mj表示为
Mj=kθ(θ−j−θ+j)=kθ2l(w∗b(j+1)+w∗b(j−1)−wb(j−1)−wb(j+1))。 (5) 而第j个接头段左右侧环段截面的弯矩M−j和M+j可分别表示为
M−j=−EId2wbdy2|i=j−=−EIwb(j−1)−2w−b(j)+w∗b(j+1)l2， (6a) M+j=−EId2wbdy2|i=j+=−EIw∗b(j−1)−2w+b(j)+wb(j+1)l2。 (6b) 结合式（2a），（5），（6a），（6b）可求得虚拟节点wb(j-1)*和wb(j+1)*的表达式分别为
w∗b(j−1)=kθlkθl+EIwb(j−1)−EIkθl+EIwb(j+1)+kθl+2EIkθl+EIw+b(j)−kθlkθl+EIw−b(j), (7a) w∗b(j+1)=kθlkθl+EIwb(j+1)−EIkθl+EIwb(j−1)+kθl+2EIkθl+EIw−b(j)−kθlkθl+EIw+b(j)。 (7b) 同理，第j个接头段左右侧环段截面的剪力Q−j和Qj+可分别表示为
Q−j=−EIw∗b(j+2)−2w∗b(j+1)+2wb(j−1)−wb(j−2)2l3， (8a) Q+j=−EIwb(j+2)−2wb(j+1)+2w∗b(j−1)−w∗b(j−2)2l3。 (8b) 结合式（2b），（7a），（7b），（8a），（8b）可求得虚拟位移节点wb(j-2)*和wb(j+2)*的表达式分别为
w∗b(j−2)=wb(j+2)−4EI+2kθlEI+kθlwb(j+1)+2kθlEI+kθlwb(j−1)+(4EI+2kθlEI+kθl−2ksl3EI)w+b(j)+(2ksl3EI−2kθlEI+kθl)w−b(j), (9a) w∗b(j+2)=wb(j−2)−4EI+2kθlEI+kθlwb(j−1)+2kθlEI+kθlwb(j+1)+(4EI+2kθlEI+kθl−2ksl3EI)w−b(j)+(2ksl3EI−2kθlEI+kθl)w+b(j)。 (9b) 结合式（7a），（7b），（9a），（9b）可以分别得到式（3）在第j个接头左右相邻侧环段节点j-和j+处的有限差分表达式：
kDtEIw−b(j)+1l4[2wb(j−2)−4EI+6kθlEI+kθlwb(j−1)−2kθlwb(j+1)EI+kθl+(4kθl+2EIkθl+EI−2ksl3EI)w−b(j)+(2ksl3EI+2kθlEI+kθl)w+b(j)]−kDtκGA(kθlkθl+EIwb(j+1)+kθlkθl+EIwb(j−1)−kθlkθl+EIw−bj−kθlkθl+EIw+bj)=DtEIq−j, (10a) kDtEIw+b(j)+1l4[2wb(j+2)−4EI+6kθlEI+kθlwb(j+1)−2kθlwb(j−1)EI+kθl+(4kθl+2EIkθl+EI−2ksl3EI)w+b(j)+(2ksl3EI+2kθlEI+kθl)w−b(j)]−kDtκGA(kθlkθl+EIwb(j+1)+kθlkθl+EIwb(j−1)−kθlkθl+EIw−bj−kθlkθl+EIw+bj)=DtEIq+j。 (10b) 同理，结合式（7a），（7b），（9a），（9b）可分别得到式（3）在第j个接头左右相邻节点j-1和j+1处的有限差分表达式：
[kθlkθl+EIwb(j+1)+(6−EIkθl+EI)wb(j−1)−kθlkθl+EIw+b(j)+(kθl+2EIkθl+EI−4)w−b(j)−4wb(j−2)+wb(j−3)]1l4−kDtκGA(w−b(j)−2wb(j−1)+wb(j−2))+kDtEIwb(j−1)=DtEIqj−1, (11a) [kθlkθl+EIwb(j−1)+(6−EIkθl+EI)wb(j+1)−kθlkθl+EIw−b(j)+(kθl+2EIkθl+EI−4)w+b(j)−4wb(j+2)+wb(j+3)]1l4−kDtκGA(w+b(j)−2wb(j+1)+wb(j+2))+kDtEIwb(j+1)=DtEIqj+1。 (11b) 整个盾构隧道两端可视为自由边界（即两端Q和M为0），则得到隧道两端4个虚拟位移节点的表达式：
wb(−2)=4wb(1)−4wb(2)+wb(3),wb(−1)=2wb1−wb2,wb(n+1)=2wb(n)−wb(n−1),wb(n+2)=4wb(n)−4wb(n−1)+wb(n−2)。} (12) 式中，n=n1n2。
结合式（12），整理式（4），（10），（11）后可写为以盾构隧道纵向弯曲位移{\boldsymbol{w}_{\text{b}}}为未知数的矩阵-向量形式：
\boldsymbol{K} \boldsymbol{w}_{\mathrm{b}}=\boldsymbol{Q}\;\;, (13) 式中：K为盾构隧道-接头刚度矩阵，Q为综合外部荷载向量，{\boldsymbol{w}_{\text{b}}}为隧道纵向弯曲位移向量。
根据式（13）可获得整个盾构隧道的纵向弯曲位移{\boldsymbol{w}_{\text{b}}}，之后由式（1a）可求得隧道纵向位移w，由式（1c）可求得隧道弯矩M和剪力Q。
2.2 隧道环间接头刚度和基床系数的确定
环间接头转动刚度kθ根据Shiba等[18]提出的纵向等效连续模型确定：
k_\mathtt{θ}=\eta \frac{E_{\mathrm{c}} I_{\mathrm{c}}}{l_{\mathrm{s}}} \frac{\cos ^3 \psi}{\cos \psi+\left(\psi+\frac{\pi}{2}\right) \sin \psi} \;。 (14a) 式中： \eta 为转动刚度系数，当不考虑环间接头的加固时，取值为1； \psi 为纵向等效连续模型中的中性轴位置，其计算公式为
\psi + \cot \psi = {\text{π }}\left( {0.5 + \frac{{{n_{\text{b}}}{k_{\text{b}}}{l_{\text{s}}}}}{{{E_{\text{c}}}{A_{\text{c}}}}}} \right) \;\;。 (14b) 式中：nb为接头处螺栓个数；kb为接头螺栓的平均线刚度，kb=EbAb/lb，Eb为螺栓的弹性模量，Ab为螺栓的横截面积，lb为螺栓长度；ls为环宽；Ec为管片的弹性模量；Ac为隧道管片横截面积；Ic隧道截面惯性矩。
根据式（14a）可进一步得到盾构隧道环间张开量 \mathit{\Delta} 的计算公式为
\mathit{\Delta}=\tan \Delta \theta(r+r \sin \psi) \approx \frac{M_j}{k_\theta}(r+r \sin \psi)\;\;, (15) 式中，r为纵向拼接螺栓中心到盾构隧道轴线的距离。
Huang等[15]基于Wu等[17]提出的等效剪切刚度得到了隧道环间接头剪切刚度ks的表达式为
{k_{\text{s}}} = \frac{{{n_{\text{b}}}{\kappa _{\text{b}}}{G_{\text{b}}}{A_{\text{b}}}\kappa GA}}{{\left( {\kappa GA - {n_{\text{b}}}{\kappa _{\text{b}}}{G_{\text{b}}}{A_{\text{b}}}} \right){l_{\text{b}}}}}\;\;。 (16) 式中：{\kappa _{\text{b}}}为螺栓的剪切系数，取0.9[16]；Gb为螺栓的剪切模量，Gb=Eb/2(1+{\nu _{\text{b}}})，{\nu _{\text{b}}}为螺栓泊松比。
根据式（2b），可以得到隧道相邻环间错台量为
\delta = w_j^ - - w_j^ + \;\;。 (17) 为考虑盾构隧道的任意埋深，地基反力系数k采用Liang[19]提出的修正值：
k=\frac{1.3 E_{\mathrm{s}}}{\omega D_{\mathrm{t}}\left(1-v_{\mathrm{s}}^2\right)} \sqrt[12]{\frac{E_{\mathrm{s}} D_{\mathrm{t}}^4}{E I}} \quad\left(\frac{H_{\mathrm{t}}}{D_{\mathrm{t}}} \geqslant 0.5\right) \;。 (18) 式中：ω为隧道埋深影响系数，ω=1+1/(1.7Ht/Dt)，Ht为隧道埋深；Es为地基弹性模量，Es=（2.5~3.5）Es0.1-0.2[19]，Es0.1-0.2为土体压缩模量；\nu _{\text{s}}^{}为土体泊松比。
3. 工程案例验证
3.1 新建隧道上穿和下穿引起的外荷载
应用上述外荷载下盾构隧道纵向变形理论于新建隧道上穿和下穿既有盾构隧道工况中，以验证所提方法的正确性与适用性。
新建隧道上穿是由于其开挖的土体重量远大于安装的管片重量，从而产生卸荷效应，进而引发下方既有隧道发生隆起变形[2-3]。可通过Mindlin应力解求得其作用于下卧盾构隧道轴线处的纵向附加荷载[2-3]。而新建隧道下穿是由于其施工造成的地层损失会引起上部地层沉降，从而引发上方既有隧道发生沉降变形。前人主要通过Peck高斯分布公式和Loganathan位移解来预测作用于既有盾构隧道上的纵向分布位移[5, 10, 13]。相关计算模型分别如图 7所示。
![]()
图 7 新建隧道上穿和下穿既有盾构隧道计算模型Figure 7. Model for existing shield tunnel due to over-crossing tunneling and under-crossing tunneling基于Mindlin应力解[2]可得新建隧道上穿引起既有隧道轴线处的附加荷载为q1(y)为
\begin{array}{l} q_1(y)=\int\limits_{-R / 2}^{R / 2} \int\limits_{-L_1}^{L_2} \frac{p \mathrm{d} \varepsilon \mathrm{d} \sigma}{8 \pi\left(1-v_{\mathrm{s}}\right)}\left[-\frac{\left(1-2 v_{\mathrm{s}}\right)\left(H_{\mathrm{t}}-H_1\right)}{R_1^3}+\right. \\ \frac{\left(1-2 v_{\mathrm{s}}\right)\left(H_{\mathrm{t}}-H_1\right)}{R_2^3}-\frac{3\left(H_{\mathrm{t}}-H_1\right)^3}{R_1^5}-\frac{30 H_1 H_{\mathrm{t}}\left(H_{\mathrm{t}}+H_1\right)^3}{R_2^7}- \\ \left.\frac{3\left(3-4 v_{\mathrm{s}}\right) H_{\mathrm{t}}\left(H_{\mathrm{t}}+H_1\right)^2-3 H_{\mathrm{t}}\left(H_{\mathrm{t}}+H_1\right)\left(5 H_{\mathrm{t}}-H_1\right)}{R_2^5}\right] 。 \end{array} (19) 式中，p为新建隧道开挖土体重量与管片及泥浆重量之差，
p=\frac{\gamma_{\mathrm{s}} \pi R_{\mathrm{s}}^2-\gamma_{\mathrm{t}} \pi\left(R_{\mathrm{o}}^2-R_{\mathrm{i}}^2\right)-\gamma_{\mathrm{n}} \pi\left(R_{\mathrm{s}}^2-R_{\mathrm{o}}^2\right)}{2 R_{\mathrm{s}}} \;。 (20) 式中：{\gamma _{\text{s}}}，{\gamma _{\text{t}}}，{\gamma _{\text{n}}}分别为开挖土体、管片和泥浆的重度；Rs，Ri，Ro分别为新建隧道开挖半径、管环内外半径；L1，L2分别为新建隧道前后端与轴线交点的距离；H1为新建隧道底部距离地面的垂直距离；{H_{\text{t}}}为既有隧道埋深；R1，R2的表达式为
{R_1} = \sqrt {{{(y\sin {\alpha _1} - \sigma )}^2} + {{(y\cos {\alpha _1} - \varepsilon )}^2} + {{({H_{\text{t}}} - {H_1})}^2}} \text{，} (21a) {R_2} = \sqrt {{{(y\sin {\alpha _1} - \sigma )}^2} + {{(y\cos {\alpha _1} - \varepsilon )}^2} + {{({H_{\text{t}}} + {H_1})}^2}} 。 (21b) 式中，{\alpha _1}为新旧隧道轴线夹角。
基于修正Loganathan位移解[10]可得新建隧道下穿开挖引起既有盾构隧道轴线处的附加荷载q2(y)为
\begin{array}{c} q_2(y)=k V_1 R_{\mathrm{x}}^2\left\{-\frac{2 H_{\mathrm{t}}\left[\left(y \cos \alpha_2\right)^2-\left(H_{\mathrm{t}}+H_2\right)^2\right]}{\left[\left(y \cos \alpha_2\right)^2+\left(H_{\mathrm{t}}+H_2\right)^2\right]^2}+\right. \\ \left.\frac{H_2-H_{\mathrm{t}}}{\left(y \cos \alpha_2\right)^2+\left(H_{\mathrm{t}}-H_2\right)^2}+\frac{\left(3-4 v_{\mathrm{s}}\right)\left(H_{\mathrm{t}}+H_2\right)}{\left(y \cos \alpha_2\right)^2+\left(H_{\mathrm{t}}+H_2\right)^2}\right\} . \\ \quad \exp \left\{-\frac{1.38\left(y \cos \alpha_2\right)^2}{\left(H_2+R_{\mathrm{x}}\right)^2}-\frac{0.69 H_{\mathrm{t}}^2}{H_2^2}\right\}\;\;。 \end{array} (22) 式中：H2为新建隧道轴线距离地表的距离；Vl为地层损失率；Rx为新建隧道半径；α2为新旧隧道轴线夹角。
将式（19），（22）分别代入式（13）即可求得新建隧道上穿和下穿引起的既有盾构隧道纵向变形和内力。
3.2 上穿案例：杭州地铁7号线上穿既有地铁1号
杭州新建地铁7号线右线隧道采用直径6.34 m的土压平衡盾构机施工，以净距1.9 m和斜角80°上穿既有地铁1号线双线隧道[20]，如图 8所示。新建隧道与既有盾构道外径均为6.2 m，内径均为5.5 m，环宽均为1.2 m，其余结构及材料参数见表 1，结合表 1和式（14a），（14b），（16）分别计算得到既有隧道接头转动刚度kθ和剪切刚度ks为6.5×107 kN·m/rad和2.2×106 kN/m。新建隧道顶部距离地表 8.2 m，既有隧道轴线距离地表 19.4 m，既有隧道左右线水平间距15.6 m。新建隧道主要穿越3-3粉砂以及3-5粉砂夹砂质粉土层，既有地铁1号线隧道主要位于3-7砂质粉土和7-2粉质黏土层，主要土层物理参数见表 2，根据表 2计算得到土层加权压缩模量Es0.1~0.2为12.04 MPa，加权泊松比 {\nu _{\text{s}}} 为0.3。既有隧道所在粉质黏土层弹性模量Es取2.5倍压缩模量值，即30 MPa。隧道掘进长度L1和L2对于左线和右线隧道分别为30，99，46，53 m。
![]()
图 8 上穿工程剖面图Figure 8. Sectional view of existing shield tunnel due to over-crossing tunneling表 1 既有盾构隧道结构参数Table 1. Structural parameters of existing shield tunnel管片弹性模量Ec/MPa 管片泊松比 {\nu _{\text{c}}} 环间螺栓数量nb 螺栓直径Db/mm 螺栓长度lb/mm 螺栓弹性模量Eb/MPa 3.45×104 0.2 17 30 400 2.06×105 表 2 主要地层物理参数Table 2. Physical parameters of main soil layers土层 重度 {\gamma _{\text{s}}} /(kN·m-3) 含水率/% 压缩模量Es0.1-0.2/MPa 泊松比 {\nu _{\text{s}}} 3-3粉砂 19.6 24.8 14.0 0.28 3-5粉砂夹粉质粉土 19.7 23.6 15.0 0.27 3-7砂质粉土 19.2 29.5 7.0 0.34 7-2粉质黏土 19.7 26.5 12.0 0.32 图 9给出了1号线双线隧道实测值与计算值对比。其中文献[2]采用EB连续梁模拟盾构隧道，结合Winkler地基建立了新建隧道上穿既有盾构隧道纵向变形解答；而文献[3]进一步考虑上穿施工中既有盾构隧道的剪切变形，建立了基于TM连续梁的盾构隧道纵向变形的计算方法。从图中可以看出本文方法和TM连续梁方法计算得到的隧道隆起值基本一致，而EB连续梁方法由于无法考虑隧道的剪切变形而略微低估了隧道隆起值。总体而言，3种方法预测得到的隧道位移与实测值较为一致。
![]()
图 9 1号线双线隧道纵向位移实测值与计算值对比Figure 9. Comparison between measured and computed displacements of tunnel of Metro Line 1进一步观察图 9中可以发现，本文方法计算得到的隧道位移曲线既不连续也不光滑，而是由一系列具有一定倾斜角度的短直线连接而成，且在接头处有明显台阶。隧道环段的位移曲线几乎为直线，这表明环段发生刚体运动，其自身弯曲和剪切变形微乎其微。环间接头部位的台阶代表环段之间的错台，前后环段存在一定的夹角，表明前后环段关于接头同时发生了转动变形。在实际工程中，隧道环段刚度较大，在上穿卸荷作用下，环段主要发生刚体位移，其自身的弯曲变形和剪切变形极小，可以忽略不计。而环间接头抵抗外部变形能力较差，上穿卸载往往会诱发其转动和错台[21]。因此，在外部荷载作用下，实际的盾构隧道纵向变形以环间接头的错台和转动为主，而环段以发生刚体位移为主。
所提计算模型可以同时捕捉隧道环段和环间接头的变形，预测结果更加符合实际情况。由于现场隧道竖向位移监测点无法密集布置，因此导致该变形现象无法监测到。相反，TM和EB连续梁方法计算得到的隧道位移曲线为光滑且连续，这是因为等效连续梁模型将盾构隧道视为一均质且连续的长梁，赋予隧道等效刚度，并不能考虑隧道环段和接头处的刚度差异。
图 10给出了基于上述3种方法计算得到的1号线右线隧道环间变形对比。从图 10（a）中可以看出，3种方法得到的隧道环间张开量的分布趋势均一致，本文方法得到的隧道环间张开量小于EB梁方法，而略高于TM连续梁方法，这是由于EB梁高估了隧道接头处的弯矩。而在图 10（b）中，本文方法得到的隧道环间错台量略小于TM连续梁方法，而EB梁方法由于剪切刚度无穷大的假设，导致错台量为0。
![]()
图 10 1号线右线隧道环间变形计算值对比Figure 10. Comparison of calculated deformations of rings of right line tunnel of Metro Line 13.3 下穿案例：上海地铁11号线下穿既有地铁4号
上海新建地铁11号线上行线依次下穿既有地铁4号线内圈和外圈隧道[4]，如图 11所示。11号线隧道与既有盾构隧道轴线水平夹角为75°，采用直径6.34 m的土压平衡盾构机施工，新旧隧道尺寸与案例3.2中一致，既有盾构隧道接头转动刚度和剪切刚度分别为6.5×107 kN·m/rad和2.2×106 kN/m。新建隧道与既有地铁4号线轴线埋深分别为25.1，17.1 m，新旧隧道竖向净距仅1.82 m。11号线隧道主要穿越⑤1黏土层和⑤3粉质黏土层，既有地铁4号线主要位于④淤泥质黏土层和⑤1a砂质粉土层，主要土层物理参数见表 3，根据文献[22]，土体加权泊松比 {\nu _{\text{s}}} 取0.26；土体弹性模量Es取20.5 MPa；土体损失率Vl取0.25%。本例中仅研究11号线隧道穿越对4号线内圈隧道的影响。
![]()
图 11 下穿工程剖面图Figure 11. Sectional view of existing shield tunnel due to under-crossing tunneling表 3 主要土层物理参数Table 3. Physical parameters of main soil layers土层 厚度/m 重度 {\gamma _{\text{s}}} /(kN·m-3) 压缩模量Es0.1-0.2/MPa 泊松比 {\nu _{\text{s}}} ④淤泥质黏土 8.4 16.7 2.09 0.33 ⑤1a砂质粉土 4.5 18.2 8.21 0.24 ⑤1黏土 4.3 17.8 3.36 0.26 ⑤3粉质黏土 5.8 18.1 4.66 0.29 图 12给出了4号线内圈隧道竖向位移实测值与计算值对比。其中Zhang等[4]采用EB连续梁模拟盾构隧道，采用Sagaseta位移公式[23]计算下穿隧道引起的附加荷载，建立了既有盾构隧道纵向变形简便解；梁荣柱等[22]将既有盾构隧道简化为TM连续梁以考虑盾构隧道的剪切变形，并结合Loganathan位移解推导了新建隧道下穿既有盾构隧道纵向变形公式；Liu等[10]将盾构隧道视为Pasternak地基上的协同变形模型，采用Loganathan位移解计算下穿隧道引起的附加荷载，建立了既有盾构隧道纵向变形解答。从图 12中可以看出本文方法和Liu等[10]的方法计算得到的隧道沉降曲线基本一致，且均较接近实测结果，这是因为所提方法和Liu等[10]的方法均可考虑盾构隧道环间接头刚度的弱化。梁荣柱等[22]的方法中由于地基反力系数的取值未考虑隧道埋深的影响，高估了既有隧道沉降量。相较于以上3种方法，Zhang等[4]的方法则预测出了更宽的隧道沉降槽以及更大的沉降值，这是因为Zhang等[4]的方法中地层沉降是采用Sagaseta位移公式计算的，而Sagaseta位移公式在计算地层沉降槽曲线时往往宽于实测结果[22]。
![]()
图 12 4号线内圈隧道竖向位移实测值与计算值对比Figure 12. Comparison between measured and calculated vertical displacements of inner tunnel of Metro Line 4进一步观察图 12发现，所提计算得到的下穿作用下隧道位移曲线既不连续也不光滑，由一系列倾斜的短直线连接而成，且在接头处有明显错台发生。而其他3种计算方法得到的隧道沉降曲线均为连续的，这是因为本文方法不仅能够考虑盾构隧道环段与接头刚度的差异，而且可以同时得到盾构隧道环段位移和环间接头的变形。环间接头处较环段更易发生错台和转动变形，而环段则发生刚体运动。以上结果表明，本文方法能够给出盾构隧道在外荷载下更符合实际纵向变形的位移曲线。
图 13给出了基于上述4种方法得到的4号线内圈隧道环间错台量对比。从图 13中可以看出，本文方法计算的隧道环间错台量小于Liu等[10]和梁荣柱等[22]的方法，且错台分布曲线较Liu等[10]的方法更加光滑。这是因为梁荣柱等[22]的方法预测出更大的隧道沉降，导致得到的错台量也大于本文方法。而在Liu等[10]的方法中假设隧道环间的错台量与其相邻环段的相对位移差为线性关系[8]，这导致尽管其方法和本文方法得到的位移曲线基本一致，但其方法却得到了更大的错台分布。Zhang等[4]的方法中使用的EB梁由于无法考虑隧道剪切变形，因此计算得到的错台量为0。
![]()
图 13 4号线内圈隧道环间错台量计算值对比Figure 13. Comparison of calculated dislocations of rings of inner tunnel of Metro Line 4图 14进一步给出了梁荣柱等[22]的方法、Liu等[10]的方法和本文方法计算得到的4号线内圈隧道剪力对比。从图 14中可以看出，3种方法得到的隧道剪力的分布趋势一致，但本文方法得到的剪力值小于Liu等[10]的方法和梁荣柱等[22]的方法，这是因为Liu等[10]的方法中剪力计算的假设条件与错台是一致的，而梁荣柱等[22]的方法高估了隧道沉降，进而也导致高估了隧道剪力。
![]()
图 14 4号线内圈隧道剪力计算值对比Figure 14. Comparison of calculated shearing forces of inner tunnel of Metro Line 44. 结论
（1）本文首先引入纵向梁-弹簧盾构隧道模型，其中隧道环段采用Timoshenko短梁模拟，环间接头采用转动和剪切弹簧模拟；然后采用有限差分法推导出了弹性地基上盾构隧道在外荷载作用下的纵向位移、内力和环间变形计算公式，并分别建立了新建隧道上穿和下穿引起既有盾构隧道纵向变形解答。
（2）上穿案例中，本文方法和Timoshenko连续梁方法预测得到的既有盾构隧道纵向位移与实测结果较为接近，而Euler-Bernoulli连续梁方法则低估了隧道隆起值；本文方法得到的隧道环间张开量略高于Timoshenko连续梁方法结果，而低于Euler-Bernoulli梁结果；Euler-Bernoulli连续梁得到的隧道环间错台量为0，而Timoshenko连续梁方法得到的隧道环间错台量略高于所提方法。
（3）下穿案例中，本文方法和协同变形模型得到的位移曲线在趋势和数值都较为一致，均接近实测结果；但本文方法得到的隧道剪力和环间错台量都小于协同变形模型，且分布曲线更加光滑；而采用Sagaseta位移公式计算地层沉降的Euler-Bernoulli连续梁方法则得到了更宽更大的隧道沉降槽。
（4）不同于既有盾构隧道模型计算得到的连续位移曲线，本文方法得到的隧道位移曲线呈现出既不光滑也不连续的特征，由一系列倾斜的短直线组成，在接头处有明显的跳跃。本文方法可以同时反映盾构隧道环段和环间接头在外荷载下的纵向变形，其中环段变形表现为刚体运动，而环间接头变形表现为转动和错台。
值得注意的是，本文提出的计算模型仅适用于外荷载作用下盾构隧道纵向变形较小的工况，当盾构隧道承受较大外荷载时，会导致接头处产生较大的弯矩和剪力，此时环间螺栓接头刚度会出现变化甚至破坏，环间变形与内力间将呈现非线性变化关系，此时本文方法不再适用。
谨对中国工程物理研究院总体工程研究所,洪建忠副主任、赵宝忠高级技师、张鸣雷工程师、汤子峰高级工程师等为设备研发的投入表示衷心感谢。 -
![]()
![]()
图 11 试验#1中孔隙水压力与时间的变化曲线
Figure 11. Variation of pore water pressure with time in test No. 1
![]()
图 12 试验#2中孔隙水压力与时间的变化曲线
Figure 12. Variation of pore water pressure with time in test No. 2
![]()
图 13 试验#1中渠基土温度与时间的变化曲线
Figure 13. Variation of soil temperature with time in test No. 1
![]()
图 14 试验#2中渠基土温度与时间的变化曲线
Figure 14. Variation of soil temperature with time in test No. 2
![]()
图 15 试验#1模型渠道法向位移与时间的变化曲线
Figure 15. Variation of normal displacement of model canal with time in test No. 1
![]()
图 16 试验#2模型渠道法向位移与时间的变化曲线
Figure 16. Variation of normal displacement of model canal with time in test No. 2
![]()
图 17 试验#1和试验#2中膨胀土渠道劣化过程
Figure 17. Deterioration process of expansive soil canals in test No. 1 and test No. 2
表 1 土工离心机技术指标
Table 1 Technical indexes of geotechnical centrifuge
技术特征 指标 离心机最大容量 400 g·t 有效半径 5.5 m 转速幅值 5~200g 有效负荷(模型箱+模型) 100g时4000 kg;200g时2000 kg 加速度稳定度 ±0.5%F·S 
下载: 导出CSV
表 2 土样基本性质
Table 2 Basic properties of soils
参数 弱膨胀土 中等膨胀土 液限/% 47.2 52.6 塑限/% 16.8 18.4 塑性指数 30.4 34.2 最优含水率/% 15.6 18.4 最大干密度/(g·cm-3) 1.80 1.70 自由膨胀率/% 52.0 71.0 注： 弱膨胀土和中等膨胀土分别为弱胀缩性的青色膨胀土和中等胀缩性的黄色膨胀土,下文以此简称。
下载: 导出CSV
表 3 试验相似率（原型/模型）
Table 3 Similarity ratios of tests (prototype/model)
物理量 相似率 物理量 相似率 加速度 1/N 位移 1 颗粒尺寸 1 温度 1 密度 1 热扩散系数 1 应力 1 导热系数 1 渗透系数 N 时间（固结） N2 渗流时间 N2 时间（热交换） N2
下载: 导出CSV
表 5 冻结过程和融化过程所需完成时间
Table 5 Time required for freeze-thaw processes
试验编号 完成时间/min 第一次循环 第二次循环 第三次循环 冻结 融化 冻结 融化 冻结 融化 #1 189 90 330 77 180 72 #2 366 104 266 93 222 77
下载: 导出CSV
表 6 湿润过程中的渠水入渗量
Table 6 Infiltration capacities of canal water during wetting process
试验编号 渠水入渗量/L 第一次湿润 第二次湿润 第三次湿润 #1 6.9 8.4 4.1 #2 8.1 7.3 5.4
下载: 导出CSV
-
[1] 邓铭江. 新疆水资源及可持续利用[M]. 北京: 中国水利水电出版社, 2005. DENG Ming-jiang. Water Resources and Sustainable Utilization in Xinjiang Uygur Autonomous Region[M]. Beijing: China Water & Power Press, 2005. (in Chinese)
[2] 蔡正银, 黄英豪. 咸寒区渠道冻害评估与处治技术[M]. 北京: 科学出版社, 2015. CAI Zheng-yin, HUANG Ying-hao. Evaluation and Treatment Technology of Frost Damage in Canals in Saline and Cold Regions[M]. Beijing: Science Press, 2015. (in Chinese)
[3] RAO S M, REDDY B V V, MUTTHARAM M. The impact of cyclic wet and dry on the swelling behaviour of stabilized expansive soils[J]. Engineering Geology, 2001, 60(1): 223-233.
[4] 许雷, 刘斯宏, 鲁洋, 等. 冻融循环下膨胀土物理力学特性研究[J]. 岩土力学, 2016, 37(增刊2): 167-174. doi: 10.16285/j.rsm.2016.S2.020 XU Lei, LIU Si-hong, LU Yang, et al. Physico-mechanical properties of expansive soil under freeze-thaw cycles[J]. Rock and Soil Mechanics, 2016, 37(S2): 167-174. (in Chinese) doi: 10.16285/j.rsm.2016.S2.020
[5] 詹良通, 吴宏伟, 包承纲, 等. 降雨入渗条件下非饱和膨胀土边坡原位监测[J]. 岩土力学, 2003, 24(2): 151-158. doi: 10.3969/j.issn.1000-7598.2003.02.034 ZHAN Liang-tong, WU Hong-wei, BAO Cheng-gang, et al. Artificial rainfall infiltration tests on a well-instrumented unsaturated expansive soil slope[J]. Rock and Soil Mechanics, 2003, 24(2): 151-158. (in Chinese) doi: 10.3969/j.issn.1000-7598.2003.02.034
[6] WANG R, ZHANG G, ZHANG J M. Centrifuge modelling of clay slope with montmorillonite weak layer under rainfall conditions[J]. Applied Clay Science, 2010, 50(3): 386-394. doi: 10.1016/j.clay.2010.09.002
[7] 程永辉, 程展林, 张元斌. 降雨条件下膨胀土边坡失稳机理的离心模型试验研究[J]. 岩土工程学报, 2011, 33(增刊1): 416-421. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2011S1079.htm CHENG Yong-hui, CHENG Zhan-lin, ZHANG Yuan-bin. Centrifugal model tests on expansive soil slope under rainfall[J]. Chinese Journal of Geotechnical Engineering, 2011, 33(S1): 416-421. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2011S1079.htm
[8] 詹良通, 刘小川, 泰培, 等. 降雨诱发粉土边坡失稳的离心模型试验及雨强-历时警戒曲线的验证[J]. 岩土工程学报, 2014, 36(10): 1784-1790. doi: 10.11779/CJGE201410004 ZHAN Liang-tong, LIU Xiao-chuan, TAI Pei, et al. Centrifuge modelling of rainfall-induced slope failure in silty soils and validation of intensity-duration curves[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(10): 1784-1790. (in Chinese) doi: 10.11779/CJGE201410004
[9] ZHANG G, YUN H, WANG L. Behaviour and mechanism of failure process of soil slopes[J]. Environmental Earth Sciences, 2015, 73(4): 1-13.
[10] 杨春宝, 朱斌, 孔令刚, 等. 水位变化诱发粉土边坡失稳离心模型试验[J]. 岩土工程学报, 2013, 35(7): 1261-1271. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201307013.htm YANG Chun-bao, ZHU Bin, KONG Ling-gang, et al. Centrifugal model tests on failure of silty slopes induced by change of water level[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(7): 1261-1271. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201307013.htm
[11] 陈湘生, 濮家骝, 罗小刚, 等. 土壤冻胀离心模拟试验[J]. 煤炭学报, 1999, 24(6): 615-619. doi: 10.3321/j.issn:0253-9993.1999.06.013 CHEN Xiang-sheng, PU Jia-liu, LUO Xiao-gang, et al. Centrifuge modelling tests of soil freezing heave[J]. Journal of China Coal society, 1999, 24(6): 615-619. (in Chinese) doi: 10.3321/j.issn:0253-9993.1999.06.013
[12] YANG D, GOODINGS D J. Climatic soil freezing modeled in centrifuge[J]. Journal of Geotechnical and Geoenviron- mental Engineering, 1998, 124(12): 1186-1194. doi: 10.1061/(ASCE)1090-0241(1998)124:12(1186)
[13] 黄英豪, 蔡正银, 张晨, 等. 渠道冻胀离心模型试验设备的研制[J]. 岩土工程学报, 2015, 37(4): 615-621. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201504007.htm HUANG Ying-hao, CAI Zheng-yin, ZHANG Chen, et al. Development of centrifugal model test facility for frost-heave of channels[J]. Chinese Journal of Geotechnical Engineering, 2015, 37(4): 615-621. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201504007.htm
[14] ZHOU J, TANG Y Q. Centrifuge experimental study of thaw settlement characteristics of mucky clay after artifical ground freezing[J]. Engineering Geology, 2015, 190: 98-108. doi: 10.1016/j.enggeo.2015.03.002
[15] 张晨, 蔡正银, 黄英豪, 等. 输水渠道冻胀离心模拟试验[J]. 岩土工程学报, 2016, 38(1): 109-117. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201601013.htm ZHANG Chen, CAI Zheng-yin, HUANG Ying-hao, et al. Centrifuge modelling of frost-heave of canals[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(1): 109-116. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201601013.htm
[16] 张晨. 渠道冻害过程模拟技术研究[D]. 南京: 南京水利科学研究院, 2016. ZHANG Chen. Simulation Technology of Frost Damage on Water Delivery Canal[D]. Nanjing: Nanjing Hydraulic Research Institute, 2016. (in Chinese)
[17] 蔡正银, 朱洵, 黄英豪, 等. 湿干冻融耦合循环作用下膨胀土裂隙演化规律[J]. 岩土工程学报, 2019, 41(8): 1381-1389. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201908002.htm CAI Zheng-yin, ZHU Xun, HUANG Ying-hao, et al. Evolution rules of fissures in expansive soils under cyclic action of coupling wet-dry and freeze-thaw[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(8): 1381-1389. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201908002.htm
[18] 唐少容, 王红雨, 潘鑫, 等. U形混凝土衬砌结构冻胀性能离心模型试验研究[J]. 农业工程学报, 2019, 35(1): 157-163. https://www.cnki.com.cn/Article/CJFDTOTAL-NYGU201901020.htm TANG Shao-rong, WANG Hong-yu, PAN Xin, et al. Frost heave performance of U-shaped canal concrete lining based on centrifuge model test[J]. Transactions of the Chinese Society of Agricultural Engineering, 2019, 35(1): 157-163. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-NYGU201901020.htm
[19] LIN L C, BENSON C H. Effect of wet-dry cycling on swelling and hydraulic conductivity of GCLs[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2000, 126(1): 40-49.
[20] 王国利, 陈生水, 徐光明. 干湿循环下膨胀土边坡稳定性的离心模型试验[J]. 水利水运工程学报, 2005(4): 6-10. https://www.cnki.com.cn/Article/CJFDTOTAL-SLSY200504001.htm WANG Guo-li, CHEN Sheng-shui, XU Guang-ming. Centrifuge model test on stability of expansive soil slope under alternation between dry and wet[J]. Hydro-Science and Engineering, 2005(4): 6-10. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLSY200504001.htm
[21] CHEN T L, ZHOU C, WANG G L, et al. Centrifuge model test on unsaturated expansive soil slopes with cyclic wet-dry and inundation at the slope toe[J]. International Journal of Civil Engineering, 2018, 18(10): 1341-1360.
[22] 张晨, 蔡正银, 徐光明, 等. 冻土离心模型试验相似准则分析[J]. 岩土力学, 2018, 39(4): 1236-1244. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201804013.htm ZHANG Chen, CAI Zheng-yin, XU Guang-ming, et al. Dimensional analysis of centrifugal modeling of frozen soil[J]. Rock and Soil Mechanics, 2018, 39(4): 1236-1244. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201804013.htm
[23] KRISHNAIAH S, SINGH D N. Determination of thermal properties of soils in a geotechnical centrifuge[J]. Journal of Testing and Evaluation, 2006, 34(4): 319-326.
[24] 曾召田, 吕海波, 赵艳林, 等. 膨胀土干湿循环效应及其对边坡稳定性的影响[J]. 工程地质学报, 2012, 20(6): 934-939. https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ201206006.htm ZENG Zhao-tian, LÜ Hai-bo, ZHAO Yan-lin, et al. Wetting-drying effect of expansive soils and its influence on slope stability[J]. Journal of Engineering Geology, 2012, 20(6): 934-939. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ201206006.htm
[25] 蔡正银, 陈皓, 黄英豪, 等. 考虑干湿循环作用的膨胀土渠道边坡破坏机理研究[J]. 岩土工程学报, 2019, 41(11): 1977-1982. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201911002.htm CAI Zheng-yin, CHEN Hao, HUANG Ying-hao, et al. Failure mechanism of canal slopes of expansive soils considering action of wetting-drying cycles[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(11): 1977-1982. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201911002.htm
[26] 齐吉琳, 马巍. 冻融作用对超固结土强度的影响[J]. 岩土工程学报, 2006, 28(12): 2082-2086. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200612006.htm QI Ji-lin, MA Wei. Influence of freezing-thawing on strength of overconsolidated soils[J]. Chinese Journal of Geotechnical Engineering, 2006, 28(12): 2082-2086. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200612006.htm
[27] 胡再强, 梁志超, 吴传意, 等. 冻融循环作用下石灰改性黄土的力学特性试验研究[J]. 土木工程学报, 2019, 52(增刊1): 211-217. https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC2019S1027.htm HU Zai-qiang, LIANG Zhi-chao, WU Chuan-yi, et al. Experimental study on mechanical properties of lime modified loess under freeze-thaw cycle[J]. China Civil Engineering Journal, 2019, 52(S1): 211-217. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC2019S1027.htm
[28] 叶万军, 杨更社, 彭建兵, 等. 冻融循环导致洛川黄土边坡剥落病害产生机制的试验研究[J]. 岩石力学与工程学报, 2012, 31(1): 199-205. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201201025.htm YE Wan-jun, YANG Geng-she, PENG Jian-bing, et al. Test research on mechanism of freezing and thawing cycle resulting in loess slope spalling hazards in Luochuan[J]. Chinese Journal of Geotechnical Engineering, 2012, 31(1): 199-205. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201201025.htm
计量
- 文章访问数: 361
- HTML全文浏览量: 25
- PDF下载量: 285
