节理巷道围岩卸荷变形试验研究

侯公羽; 邵耀华; 刘春雷; 王一哲; 陈钦煌; 胡志宇; 赵庆儒

doi:10.11779/CJGE20230028

节理巷道围岩卸荷变形试验研究 English Version

中国矿业大学(北京)力学与土木工程学院, 北京 100083

基金项目:

国家自然科学基金委员会与2020年度高速铁路基础研究联合基金资助重点项目 U2034205

详细信息

作者简介:
侯公羽(1965—)，男，教授，博士生导师，主要从事煤矿开采、岩土工程和岩石力学等方面的教学与研究工作。E-mail: hgyyht@126.com

通讯作者:
邵耀华, E-mail: syh9376@126.com

中图分类号: TU432
计量
- 文章访问数: 337
- HTML全文浏览量: 50
- PDF下载量: 104
出版历程
- 收稿日期: 2023-01-08
- 网络出版日期: 2024-03-14
- 刊出日期: 2024-02-29

Unloading deformations of surrounding rock in jointed roadways

School of Mechanics and Civil Engineering, China University of Mining and Technology, Beijing 100083, China

摘要

摘要: 自然界中的岩体在长期的地质构造运动作用下，内部往往发育不同尺度的节理，而巷道围岩的稳定性与节理的发育密切相关。为研究节理对巷道围岩开挖卸荷变形的影响，采用课题组的开挖卸荷模型试验系统，对不同节理尺寸水泥砂浆围岩试件开展了系列试验。研究结果表明：①巷道围岩开挖卸荷变形主要发生在卸荷阶段，维持阶段变形继续存在，但节理对其卸荷变形无明显影响。此阶段，卸荷变形程度显著降低并逐渐趋于稳定，其应变增量大约为100 με；②与无节理试件相比，含节理试件S-1和S-2的应变平均降低幅度超过70%，远高于S-3，即节理尺寸越大，卸荷变形程度越大，卸荷作用越明显；③试件在不同方向的应变差异较大，切向应变均大于径向应变，含节理试件尤为突出。
- 小型模型试验 /
- 节理 /
- 开挖卸荷 /
- 围岩变形
Abstract: The rock mass in nature often develops joints of different scales under the action of long-term geological tectonic movement, and the stability of surrounding rock of roadways has a great relationship with the development of joints. For the sake of studying the effects of joints on the excavation unloading deformation of surrounding rock of roadways, a series of tests are carried out on the cement mortar surrounding rock specimens with different joint sizes through the excavation unloading model test system developed by the authors. The results show that: (1) The unloading deformations of surrounding rock mainly occur at the unloading stage, and the deformations continue to exist at the maintenance stage, but the joints have no obvious influences on the unloading deformations. At this stage, the unloading deformation degree decreases significantly and gradually becomes stable, with the strain increment of about 100με. (2) Compared with the specimens without joints, the average strain reduction of S-1 and S-2 specimens with joints is more than 70%, which is much higher than S-3. That is, the larger the joint size, the greater the degree of unloading deformations, and the more obvious the unloading effects. (3) The specimens have great differences in strain in different directions, and the tangential strain is greater than the radial strain, especially for the specimens with joints.
- small-scale model test /
- joint /
- excavation unloading /
- surrounding rock deformation

HTML全文

0. 引言

近年来，中国地下综合管廊快速发展，地下管廊抗震分析受到了人们的重视，学者们通过振动台模型试验或数值模拟等方法对地下管廊的抗震性能进行了大量研究。李杰等^[1]根据振动台试验中地基土的实际边界条件，提出一种变刚度方法模拟层状剪切土箱，并验证了模拟方法的可靠性。汤爱平等^[2]通过振动台试验，讨论了管廊体系及其周围土体的地震反应，结果表明管体与其内部管线的动力响应存在显著不同。郭恩栋等^[3]通过数值模拟方法研究了地下管廊体系的横断面地震响应，提出与管廊内部管线相比，应重点关注管廊自身的抗震设计。由于地下管廊多被设计为浅埋结构，需穿越不同类别场地，因此一些学者研究了面波作用下或不同类型场地中管廊的地震响应。施有志等^[4]讨论了Rayleigh波作用下，地下管廊的动力特征。蒋录珍等^[5]研究了饱和土体内地下管廊的动力响应，讨论了不同地震动峰值、地震波入射角度及土体孔隙率等因素的影响。李金奎等^[6]通过等价线性本构模拟地基土的动力非线性特性，给出了地震作用下，管廊自身的薄弱位置。阮滨等^[7]研究了可液化场地中地下管廊结构的变形响应，结果表明管廊接头张开量在液化区显著增大。此外，由于装配式管廊的预应力承插式接头与盾构隧道结构多采用的拼装式螺栓接头在力学性能上有很大不同，谷音等^[8]、梁建文等^[9]研究了管廊预制接头的建模及接头预应力的实现方法，讨论了预制接头对地下管廊自身抗震性能的影响。

值得注意的是，与地铁隧道等地下交通结构不同，地下管廊相互交叉，形成大量十字或T型交叉节点；由于交叉管廊在两个主轴方向上的抗侧移刚度相差较大，交叉节点的地震内力非常复杂，成为地下管廊抗震性能的薄弱点。然而，目前中国还没有专门的地下管廊抗震设计规范，而《地下结构抗震设计标准》等现行标准^[10-11](以下简称《标准》)也仅给出了直线型地下(隧道)结构纵向抗震设计方法。

关于交叉管廊的抗震分析，近期有一些初步研究。赵丹阳^[12]计算了某十字交叉管廊的动力响应，讨论了单向及双向地震作用下交叉节点的破坏模式。梁建文等^[13-14]基于拟静力方法和动力时程方法，计算了某T型交叉管廊的地震反应，给出了地震作用下，T型交叉节点的内力放大系数及影响范围。由于交叉管廊体系中两个主轴结构的计算长度一般较大，由地震作用引起的两个主轴结构的拉压及挠曲等纵向地震反应对交叉节点的影响不可忽略，然而上述研究^[12-14]均是横向抗震分析，研究交叉管廊交叉节点处顶板和底板之间相对变形，受限于地震动输入机制，这些研究均没有涉及交叉管廊的纵向抗震设计。

为此，本文提出地下交叉管廊纵向抗震设计的一个反应位移法，并以预制装配式T型交叉管廊为例，进行了算例分析。研究可为交叉管廊的纵向抗震设计提供一定参考。

1. T型交叉管廊纵向反应位移法

现行《标准》中规定，地震作用下，地层沿结构轴线方向的纵向位移u_x及与结构轴线垂直方向的横向位移u_y可采用正弦规律分布，如图 1所示，场地的位移分布函数u(x', z')为

图 1 《标准》中纵向反应位移法计算模型

Figure 1. Computational model for longitudinal response deformation method in Chinese standard

下载: 全尺寸图片幻灯片

$u(x',z') = {u_{\max }}(z')\sin \left( {\frac{{2{\text{π }}}}{L}x'} \right) 。$

(1)

式中：φ为地震波传播方向与地下结构轴向的夹角；u_max为地震作用下场地深度z(z=z')处的水平位移峰值；L为地层变形波长。

可以看到，纵向反应位移法基本假定下，当给定φ和u_max两个参数后，观测场地同一深度z处沿纵向将呈现出正弦模式下的单一地层变形，仅一次计算即可获得结构的地震响应峰值。值得注意的是，该方法仅适用于单线地下结构，当研究对象为交叉管廊时，上述单一地层变形思想将不再适用。

鉴于此，本文基于反应位移法基本原理，提出一种适用于交叉管廊纵向抗震设计的反应位移法。以某T型管廊为例，假定某一地层变形状态下，管廊所在场地的地层变形形式如所示。图中， $u_x^{\text{L}}$ ， $u_x^{\text{T}}$ 表示管廊结构所在位置地基土的纵向(x向)变形， $u_y^{\text{L}}$ ， $u_y^{\text{T}}$ 表示结构所在地基土的横向(y向)变形，上标L和T分别表示x向结构和y向结构。o-x'-y'-z'为入射波坐标系，o-x-y-z为场地坐标系。

图 2 T型交叉管廊地层位移分解模型

Figure 2. Decomposition of stratum deformation for T-type utility tunnels

下载: 全尺寸图片幻灯片

定义o-x'-y'-z'坐标系中任意一点( ${x'_i}$ , ${y'_j}$ , ${z'_k}$ )，经坐标变换后，其与场地坐标系o-x-y-z的关系为

$\left\{ \begin{gathered} {{x'}_i} \hfill \\ {{y'}_j} \hfill \\ {{z'}_k} \hfill \\ \end{gathered}\right\} = \left[ {\begin{array}{*{20}{c}} {\cos \varphi }&{\sin \varphi }&0 \\ { - \sin \varphi }&{\cos \varphi }&0 \\ 0&0&1 \end{array}} \right]\left\{ \begin{gathered} {x_i} \hfill \\ {y_j} \hfill \\ {z_k} \hfill \\ \end{gathered}\right\} 。$

(2)

点( ${x'_i}$ , ${y'_j}$ , ${z'_k}$ )的位移u( ${x'_i}$ , ${y'_j}$ , ${z'_k}$ )与其在o-x-y-z坐标系内的分量u_x，u_y和u_z的关系可为

$\left\{ {\begin{array}{*{20}{c}} {{\boldsymbol u_x}({x_i},{y_j},{z_k})} \\ {{\boldsymbol u_y}({x_i},{y_j},{z_k})} \\ {{\boldsymbol u_z}({x_i},{y_j},{z_k})} \end{array}} \right\} = \left\{ {\begin{array}{*{20}{c}} { - \sin \varphi } \\ {\cos \varphi } \\ 1 \end{array}} \right\}\boldsymbol u({x'_i},{y'_j},{z'_k}) \text{，}$

(3)

式中， ${z'_k}$ = ${z_k}$ 。

将式(2)，(3)代入式(1)，可得u(x', y', z')在场地坐标系o-x-y-z下的分布函数，其沿x向和y向的位移函数形式如式分别为

$u_x(x, y, z)=-u_{\max }(z) \sin \varphi \sin \left[\frac{2 \pi}{L}(y \sin \varphi+x \cos \varphi)\right]，$

(4)

$u_y(x, y, z)=u_{\max }(z) \cos \varphi \sin \left[\frac{2 \pi}{L}(y \sin \varphi+x \cos \varphi)\right] 。$

(5)

值得注意的是，式(4)，(5)仅为图 2示意单一状态下场地的x向及y向地层变形分布，不足以反映整个地层变形周期内T型交叉管廊结构的全部地震响应。

与单线地下结构相比，交叉型地下结构纵向抗震设计的一个显著区别在于，一次分析计算无法反映整个地层变形周期内结构的变形历程。由于交叉节点的存在，交叉型地下结构在处于同一周期内不同状态下的地层变形分布时，结构峰值响应可能会有很大区别。因此，明确交叉型地下结构两个主轴方向各自的地层变形计算方法，设计可反映其周期内变形历程的计算模型十分必要。

为体现出一个地层变形周期内结构的变形历程，本文引入计算步长i，通过改变位移函数相位来实现场地地层变形的周期性运动，分n步逐次令场地运动一个地层变形周期，这里约定图 2示意为整个地层变形周期的初始状态(i=0)，那么整个周期内所有状态下场地的地层变形分布为

${u}_{x}(x,y,z)=-{u}_{\mathrm{max}}\text{(}z\text{)sin}\varphi \mathrm{sin}\left[\frac{2\text{π}}{L}(y\mathrm{sin}\varphi +x\mathrm{cos}\varphi -\frac{i}{n}L)\right]\text{，}$

(6)

${u}_{y}(x,y,z)={u}_{\mathrm{max}}\text{(}z\text{)cos}\varphi \mathrm{sin}\left[\frac{2\text{π}}{L}\left(y\mathrm{sin}\varphi +x\mathrm{cos}\varphi -\frac{i}{n}L\right)\right]。$

(7)

式中，i=0, 1, 2, ..., n，n为计算步长总数。

可以看到，i=0时，由式(6)，(7)计算得到的地层变形分布即为图 2示意的初始地层变形模式；随着计算步长i的增加，两个主轴方向管廊轴线处的地层变形在一个地层变形周期内逐渐偏移，直至i=n时完成一整个周期的运动历程。图 3展示了由式(6)，(7)计算得到的场地周期性变形。图 3中，φ取30°，T型交叉管廊x向及y向尺寸分别取为2L和L，以i/n=0，1/4，1/2，3/4四个周期状态为例。

图 3 T管廊所在位置地层变形的周期性分布

Figure 3. Periodic distribution of stratum deformation at location of T-type utility tunnels

下载: 全尺寸图片幻灯片

可以看到，本文方法可以反映场地在一个地层变形周期内的运动历程。因此，仅需确定适宜的计算步长总数n的取值，便可保证获取足够丰富的场地变形细节，关于n的取值问题将在下文详细讨论。

结合上文计算方法，本文于ABAQUS有限元平台编制了用于交叉管廊抗震设计的Python程序，实现了从模型建立→地基弹簧及接头弹簧添加→周期性地层变形赋予等全流程作业，从建模到输出结果自动运行，无需外部干预。

2. T型交叉管廊纵向抗震分析

2.1 算例

以某拟建的T型交叉综合管廊为例，该T型管廊两个主轴方向横截面均为双舱形式，总宽度11000mm，高5750 mm，底板厚600 mm，顶板和侧板厚500 mm，中隔板厚300 mm。管廊混凝土强度等级为C60，密度2650 kg/m³，弹性模量35.5 GPa。管廊纵向接头通过承插口连接，预制管廊段长1.5 m，腋角处共配置6条钢绞线，有效截面积140 mm²，弹性模量195000 MPa，屈服强度1628 MPa。管廊处于弹性基岩上的软土场地，计算埋深3 m，场地总厚度为80 m，共包含4类土，具体土层参数信息及土体非线性特性参数见文献[15]。

2.2 计算模型

采用三维线性B31梁单元模拟管廊结构，单元尺寸0.75 m。结合《标准》中相关规定，模型纵向(x向)计算长度取2个地层变形波长(2L=1066 m)，横向(y向)计算长度取一个地层变形长(L=533 m)，模型中T型节点为整体现浇式，其余部分为预制装配式。

于管节间承插式接头处将梁单元断开，采用非线性弹簧单元模拟接头力学行为，其中3个方向平动弹簧用于模拟预制接头的轴向及横向行为，两个方向转动弹簧用于模拟接头弯曲行为，忽略接头扭转行为。地基土采用沿地下结构轴向均匀布置的纵向弹簧和横向弹簧模拟，如图 4所示。

图 4 计算模型示意

Figure 4. Computational model

下载: 全尺寸图片幻灯片

考虑到承插式接头构造，假定由管廊腋角6根钢绞线承担全部拉力，止水橡胶圈及接头混凝土协同承担压力和剪力，接头弹簧刚度系数由K_i=E_iA/l计算，其中K_i为弹簧不同阶段的刚度系数，E_i为接头处混凝土或止水橡胶或钢绞线的弹性模量，A为钢绞线截面积或混凝土承压面积，l为钢绞线有效长度。接头转动时，假定断面受拉区由钢绞线承担拉力，受压区由橡胶混凝土协同承担压力，结构始终处于弹性状态，采用静力有限元法^[16]计算两个方向转动弹簧刚度。接头弹簧参数具体计算方法见文献[9]。

2.3 场地地震动反应及模型地震动参数

入射波选用El Centro波，其加速度峰值调幅为0.20g，本文采用一维地震分析程序EERA^[17]计算场地地震反应，计算得到的模型地震动参数如表 1所示。

表 1 计算模型地震动参数

Table 1. Parameters of model ground motion

u_max/m	G_d/MPa	k_t/(kN·m^-1)	k_l/(kN·m^-1)
0.0934	7.3723	190758.52	370458.58

下载: 导出CSV

| 显示表格

采用本文方法于地基弹簧远端逐次输入地层位移，直至完成一个周期的地层变形，计算地震波传播方向与x向管廊轴线夹角φ=0°，30°，45°，60°四种情况下的结构响应。

2.4 计算步长总数的确定

本文通过改变位移函数的方式分步逐次令场地运动一个地层变形周期，为了确定适宜的计算步长总数n，保证获取足够丰富的结构变形细节且兼顾计算机时，计算了步长总数n=4，8，16，32四类情况(以入射角φ=0°为例)的T型交叉管廊地震反应。以T型交叉节点处的接头变形为观测对象，提取了图 5所示3个监测点的接头张开量，其中每个监测点接头张开量取为接头相邻节点的位移差Δu与节点转动弧长l(l_A+l_B)之和，具体计算方法如图 6所示。图 6中，u_x为节点位移，上标A和B为接头处梁单元节点；θ_A和θ_B分别为节点A和节点B相对初始位置的转角；H_W为管廊顶、底板宽度。

图 5 交叉节点处监测点示意

Figure 5. Selection of monitoring points at T-type cross

下载: 全尺寸图片幻灯片

图 6 接头挠曲变形示意及张开量计算方法

Figure 6. Schematic diagram of joint deflection and method for joint opening

下载: 全尺寸图片幻灯片

图 7给出了各监测点张开量在一个地层变形周期内的变化如所示。可以看到，3个监测点的张开量在不同周期状态(i/n)下存在显著差异。随着计算步长总数n的增加，各监测点的周期内变形响应逐渐趋同。当步长数n≥8后，各张开量曲线已基本重合，足以捕捉变形响应峰值，且满足计算精度要求。因此，为兼顾计算机时，本文将计算步长总数取n为8。

图 7 不同工况下各监测点的周期内变形

Figure 7. Periodic deformations of monitoring points

下载: 全尺寸图片幻灯片

3. 计算结果

一个地层变形周期内T型交叉管廊的接头变形及结构内力响应沿管廊轴线的分布(以x向管廊为例，限于篇幅，仅展示了φ=0°时，i/n=0，1/4，1/2，3/4等4个状态下的结构地震响应分布)，如图 8所示。

图 8 T型交叉管廊内力和变形沿管廊轴线的分布

Figure 8. Distribution of internal forces and deformations along axis of T-type cross utility tunnels

下载: 全尺寸图片幻灯片

可以看到，T型节点处因存在截面突变，整个地层变形周期内均表现出显著的变形及内力集中，T型交叉节点的影响范围约为20 m。地震波传播方向与x向管廊轴线夹角φ=0°时，x向管廊虽仅承担弯曲变形，但图 8中，T型交叉节点处在初始状态(i/n=0)及1/2周期状态(i/n=1/2)下仍表现出较大的轴力，其它地震响应也均出现不同程度的放大现象，体现出交叉节点的影响。随着地层变形的周期内推移，T型节点处的地震响应不断变化且不同周期状态i/n下的幅值差距较大，表明周期内不同的地层变形分布形式对交叉型地下结构的地震响应影响显著，因此确定交叉节点处结构响应达到最不利状态时的地层变形十分必要。

为获取地震作用下，由本文方法计算得到的T型交叉管廊的响应峰值及其对应的结构最不利模式，根据图 8所示T型交叉节点的变形及受力特征，提取了不同地震波入射角度φ下，T型交叉节点处的接头的变形响应以及结构的内力响应，各类响应峰值在一个地层变形周期内的变化如图 9所示(以φ=0°为例，监测点选取与图 5相同)。

图 9 T型交叉节点处内力和变形的周期内变化

Figure 9. Periodic variation of internal forces and deformations

下载: 全尺寸图片幻灯片

可以看到，在一个地层位移周期内，不同阶段、不同变形模式下T型交叉节点处的变形及内力响应存在显著不同。地震波入射角度φ=0°时，管廊的内力响应峰值均出现在初始周期状态(i/n=0)，而变形响应则在初始周期状态(i/n=0)及1/2周期状态(i/n=1/2)均表现出较大值，监测点3峰值张开量达到4.1 mm，已超出《城市综合管廊工程技术规范：GB50838— 2015》^[18]规定的限值2 mm。

不同地震波入射角度φ下，各监测点的地震响应峰值及其对应周期状态(i/n)如表 2所示。出于简便，这里定义轴力最不利状态对应步长i_f，剪力最不利状态对应步长i_s，弯矩最不利状态对应步长i_m；接头变形的不利状态对应步长i_d。

表 2 T型交叉节点最不利状态下的峰值响应

Table 2. Peak responses of T-type cross under most unfavorable state

φ/(°)	峰值轴力/MN	i_f/n	峰值剪力/MN	i_s/n	峰值弯矩/(MN·m)	i_m/n	峰值张开量/mm	i_d/n
0	-0.72(监测点2)	0	1.05(监测点3)	0	3.96(监测点3)	0	4.13(监测点3，超限)	0
30	-5.59(监测点2)	7/8	1.33(监测点1)	3/8	3.28(监测点1)	3/8	3.66(监测点1，超限)	7/8
45	-5.94(监测点2)	3/4	1.20(监测点2)	1/4	2.64(监测点2)	1/4	3.31(监测点2，超限)	1/4
60	-5.61(监测点1)	1/2	1.33(监测点2)	0	3.27(监测点2)	0	3.67(监测点2，超限)	1/2

下载: 导出CSV

| 显示表格

现有研究^[19]表明，采用纵向反应位移法进行线长形地下结构的地震反应计算时，随着入射角度φ的增大，结构的弯矩、剪力逐渐减小，轴力及张开量先增大后减小。但从表 2中可以看出，T型交叉管廊的地震响应与上述规律存在显著区别。可以看到，地震波入射角度φ=0°时，x向管廊主要承担弯曲变形，而y向管廊表现为整体平动，当两者变形耦合在一起后，主要体现为整体平动的y向管廊出现了较大的接头张开，这是由于两个主轴方向管廊结构的抗侧向刚度及结构周围地基弹簧的约束能力相差较大，x向管廊在周期内弯曲过程中，y向交叉节点出现较大接头转动导致。

随着入射角度的增大，管廊的接头变形峰值和截面弯矩峰值有所减小，而轴力和剪力峰值有所增大。其中，截面轴力峰值的变化最为显著，φ=45°时监测点2处轴力达到峰值5940 kN，表明双向地层变形的耦联运动对T型交叉节点处轴力有不利影响。

4. 敏感性分析

4.1 正交试验方案设计

为了进一步研究地震作用下，T型交叉管廊的变形及受力机理，明确不同地层条件、地震动参数等因素对本文方法计算结果及T型管廊最不利模式的影响，获取一定普适性结论，拟通过正交试验对本文方法进行敏感性分析，考虑了包括场地剪切波速v_s、入射角度φ、地震峰值加速度a_max以及管廊计算埋深h在内的4个因素的交互影响，每个因素设计为3个水平。

由于城市综合管廊处地层的等效剪切波速v_s通常小于500 m/s，且多为浅埋结构，故选取剪切波速v_s为200~400 m/s的场地进行分析，设计结构埋深h为2~4 m。同时，仅讨论入射角φ > 0°的情况，地震动峰值调幅为0.10g~0.20g。综上，设计了L9(3⁴)正交表，如表 3所示。假定场地覆土深度H为80 m，质量密度ρ为2000 kg/m³，场地地层变形波长L为400 m，暂不考虑地基土的非线性影响，各工况模型参数计算方法与上文相同。

表 3 正交试验方案

Table 3. Scheme of orthogonal experiments

工况	v_s/(m·s^-1)	φ/(°)	a_max/g	h/m
1	200	30	0.10	2
2	200	45	0.15	3
3	200	60	0.20	4
4	300	30	0.15	4
5	300	45	0.20	2
6	300	60	0.10	3
7	400	30	0.20	3
8	400	45	0.10	4
9	400	60	0.15	2

下载: 导出CSV

| 显示表格

4.2 正交试验计算结果

计算了各正交试验工况，得到的各工况下T型交叉管廊的内力和变形峰值响应及其对应周期状态(i/n)如表 4所示。

表 4 正交模型结果汇总

Table 4. Results of orthogonal model

工况	峰值轴力/MN	i_f/n	峰值剪力/MN	i_s/n	峰值弯矩/(MN·m)	i_m/n	峰值张开量/mm	i_d/n
1	-6.85(监测点2)	7/8	2.20(监测点1)	3/8	3.76(监测点1)	3/8	2.22(监测点3)	7/8
2	-10.25(监测点2)	3/4	2.42(监测点2)	1/4	3.52(监测点2)	1/4	2.10(监测点2)	1/4
3	-13.69(监测点1)	1/2	4.40(监测点2)	0	7.53(监测点2)	0	4.43(监测点3)	1/2
4	-11.72(监测点2)	7/8	3.92(监测点1)	3/8	6.02(监测点1)	3/8	2.41(监测点1)	7/8
5	-15.48(监测点2)	3/4	3.69(监测点2)	1/4	4.78(监测点2)	1/4	1.85(监测点2)	1/4
6	-7.85(监测点1)	1/2	2.63(监测点2)	0	4.03(监测点2)	0	1.61(监测点3)	1/2
7	-18.44(监测点2)	7/8	6.20(监测点1)	3/8	8.87(监测点1)	3/8	2.55(监测点3)	7/8
8	-9.05(监测点2)	3/4	2.14(监测点2)	1/4	2.58(监测点2)	1/4	0.70(监测点2)	1/4
9	-13.87(监测点2)	1/2	4.66(监测点2)	0	3.88(监测点2)	0	6.68(监测点2)	1/2

下载: 导出CSV

| 显示表格

可以看到，四类因素交互作用下，入射角φ相同时，T型交叉节点处内力和变形峰值的最不利模式一致，即同一地震波入射角度φ下，改变其它因素不会影响内力和变形出现峰值时所对应的周期状态(i/n)，表明线弹性假定下，地震波入射角度φ对T型管廊的最不利模式起控制作用。利用本文方法进行T型交叉管廊纵向抗震计算时，需着重验算表 5中最不利变形下的地震响应。

表 5 交叉节点处截面轴力峰值极差分析

Table 5. Range analysis of axial force at cross section

水平	因子
水平	v_s/(m·s^-1)	φ/(°)	a_max/g	h/m
1	10.263	12.337	7.917	12.067
2	11.683	11.593	11.947	12.180
3	13.787	11.803	15.870	11.487
极差	3.523	0.743	7.953	0.693
影响度排序	2	3	1	4

下载: 导出CSV

| 显示表格

根据正交试验结果，计算得到的每个因素在不同水平情况下各个指标的平均值及极差如表 5~8所示。可以看到，加速度峰值对T型交叉管廊的内力和变形均影响显著；入射角度是影响T型交叉管廊挠曲及接头变形的主要因素，但对交叉节点处轴力峰值的影响较小；埋深及剪切波速对交叉节点处内力和变形的影响不明显。因此，对交叉型管廊进行纵向地震反应计算时，需着重验算不同入射角φ下交叉管廊的内力和变形。

表 6 交叉节点处截面剪力峰值极差分析

Table 6. Range analysis of shear force at cross section

水平	因子
水平	v_s/(m·s^-1)	φ/(°)	a_max/g	h/m
1	3.007	4.107	2.323	3.517
2	3.413	2.750	3.667	3.750
3	4.333	3.897	4.763	3.487
极差	1.327	1.357	2.440	0.263
影响度排序	3	2	1	4

下载: 导出CSV

| 显示表格

表 7 交叉节点处截面弯矩峰值极差分析

Table 7. Range analysis of bending moment at cross section

水平	因子
水平	v_s/(m·s^-1)	φ/(°)	a_max/g	h/m
1	4.937	6.217	3.457	4.140
2	4.943	3.627	4.473	5.473
3	5.110	5.147	7.060	5.377
极差	0.173	2.590	3.603	1.333
影响度排序	4	2	1	3

下载: 导出CSV

| 显示表格

表 8 交叉节点处接头张开量峰值极差分析

Table 8. Range analysis of joint opening at cross section

水平	因子
水平	v_s/(m·s^-1)	φ/(°)	a_max/g	h/m
1	2.917	2.393	1.510	3.583
2	1.957	1.550	3.730	2.087
3	3.310	4.240	2.943	2.513
极差	1.353	2.690	2.220	1.497
影响度排序	4	1	2	3

下载: 导出CSV

| 显示表格

5. 结论

基于纵向反应位移法基本思想，提出了一种适用于T型交叉管廊纵向抗震设计方法。通过赋予地基弹簧周期性地层变形的方式，研究了T型交叉管廊的周期内变形历程，并据此讨论了交叉节点处变形和内力的最不利模式，以及计算参数对计算结果的敏感性，得到以下3点结论。

(1) 本文通过改变位移函数的方式，分为n步逐次令场地运动一个地层变形周期，可捕捉交叉型地下结构的峰值地震响应。讨论了计算步长n对计算结果精度的影响，建议计算步长n不宜小于8步。

(2) 在一个地层位移周期内，不同阶段、不同模式下T型管廊交叉节点处的变形及内力响应存在显著不同，即不同物理量的最不利模式不同。讨论了不同地震波入射角度下，T型交叉管廊的地震峰值响应。结果表明入射角对交叉节点处的轴力存在显著不利影响，随着入射角度的增大，结构的挠曲变形有所减小。给出了线弹性假定下，不同地震波入射角度时，各类最不利工况对应周期状态及地层变形，可为相关地下结构抗震设计提供参考。

(3) 考虑到场地剪切波速、入射角度、地震峰值加速度以及管廊计算埋深等四个因素的交互影响，设计了相应正交试验方案，讨论了不同参数对T型交叉管廊最不利模式的敏感性。结果表明，同一入射角度下，T型交叉节点处各类峰值响应的最不利模式一致，改变其它因素不会影响各类物理量出现峰值时所对应的周期状态，地震波入射角度对T型管廊的最不利模式起控制作用。

需要指出的是，本文计算方法并不局限于“T型”交叉管廊的地震反应计算，该计算思想经过一定拓展，也可应用于“十字型”、“L型”等交叉管廊结构，乃至地下管廊网络的抗震分析中，具体将另文讨论。

图 1 密封结构图

Figure 1. Sealing components

下载: 全尺寸图片幻灯片

图 2 小型围岩开挖卸荷系统

Figure 2. Small surrounding rock excavation unloading system

下载: 全尺寸图片幻灯片

图 3 围岩试件初始应力示意图

Figure 3. Mechanism model for rectangular roadway under excavation unloading

下载: 全尺寸图片幻灯片

图 4 应变片布置示意图

Figure 4. Arrangement of strain gauges

下载: 全尺寸图片幻灯片

图 5 应变片布置示意图

Figure 5. Arrangement of strain gauges

下载: 全尺寸图片幻灯片

图 6 岩芯试样应力-应变曲线

Figure 6. Stress-strain curves of core specimens

下载: 全尺寸图片幻灯片

图 7 预制节理试件

Figure 7. Prefabricated joint specimens

下载: 全尺寸图片幻灯片

图 8 应力路径

Figure 8. Stress paths

下载: 全尺寸图片幻灯片

图 9 T1、T2和T3测点应变-时间曲线

Figure 9. Strain-time curves of points T1, T2 and T3

下载: 全尺寸图片幻灯片

图 10 T1测点无节理试件理论值与试验值对比

Figure 10. Comparison between theoretical and experimental values of specimen without joints at point T1

下载: 全尺寸图片幻灯片

表 1 水泥砂浆的力学参数

Table 1 Mechanical parameters of cement mortar test pieces

介质	泊松比μ	弹性模量/GPa	抗压强度/MPa	密度/(g·cm^-3)
水泥砂浆	0.2	7.7	25	1.97
砂岩	0.02~0.2	4.9~78.5	20~170	2.10

下载: 导出CSV

表 2 试验方案

Table 2 Testing programs

试样编号	数量	节理尺寸/mm	模拟地应力/MPa
S-0	3	—	25
S-1	3	40×90
S-2	3	50×90
S-3	3	60×90

下载: 导出CSV

表 3 T1、T2测点径向和切向平均应变值

Table 3 Radial and sheer strain values of points T1 and T2

测点	试样编号	卸荷阶段应变峰值/10^-3	卸荷后维持阶段应变增量/10^-3
T1	S-0-R	2824.18	97.12
	S-1-R	621.39	98.98
	S-2-R	754.57	96.23
	S-3-R	1556.75	103.69
	S-0-S	-2824.18	97.12
	S-1-S	-827.57	102.43
	S-2-S	-958.32	104.94
	S-3-S	-1785.31	95.65
T2	S-0-R	1841.84	100.44
	S-1-R	352.54	96.21
	S-2-R	458.59	98.08
	S-3-R	547.83	99.79
	S-0-S	-1841.84	100.44
	S-1-S	-441.43	100.96
	S-2-S	-560.83	99.21
	S-3-S	-775.98	95.10

下载: 导出CSV

参考文献(28)

[1]	侯公羽, 梁金平, 李小瑞. 常规条件下巷道支护设计的原理与方法研究[J]. 岩石力学与工程学报, 2022, 41(4): 691-711. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202204002.htm HOU Gongyu, LIANG Jinping, LI Xiaorui, et al. Research on principles and methods of roadway support design under conventional conditions[J]. Chinese Journal of Rock Mechanics and Engineering, 2022, 41(4): 691-711. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202204002.htm
[2]	刘明星, 廖孟柯, 刘恩龙, 等. 不同结构面倾角节理岩体三轴试验力学特性[J]. 北京工业大学学报, 2018, 44(3): 336-343. https://www.cnki.com.cn/Article/CJFDTOTAL-BJGD201803003.htm LIU Mingxing, LIAO Mengke, LIU Enlong, et al. Experimental research on mechanical properties of jointed rock mass with different angles of inclination[J]. Journal of Beijing University of Technology, 2018, 44(3): 336-343. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BJGD201803003.htm
[3]	杨圣奇, 戴永浩, 韩立军, 等. 断续预制裂隙脆性大理岩变形破坏特性单轴压缩试验研究[J]. 岩石力学与工程学报, 2009, 28(12): 2391-2404. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200912005.htm YANG Shengqi, DAI Yonghao, HAN Lijun, et al. Uniaxial compression experimental research on deformation and failure properties of brittle marble specimen with pre-existing fissures[J]. Chinese Journal of Rock Mechanics and Engineering, 2009, 28(12): 2391-2404. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200912005.htm
[4]	王乐华, 柏俊磊, 孙旭曙, 等. 不同连通率节理岩体三轴加卸荷力学特性试验研究[J]. 岩石力学与工程学报, 2015, 34(12): 2500-2508. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201512014.htm WANG Lehua, BAI Junlei, SUN Xushu, et al. The triaxial loading and unloading mechanical properties of jointed rock masses with different joint connectivities[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(12): 2500-2508. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201512014.htm
[5]	杨圣奇, 陆家炜, 田文岭, 等. 不同节理粗糙度类岩石材料三轴压缩力学特性试验研究[J]. 岩土力学, 2018, 39(增刊1): 21-32. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2018S1004.htm YANG Shengqi, LU Jiawei, TIAN Wenling, et al. Experimental study of mechanical behavior of rock specimens with different joint roughness coefficient under conventional triaxial compression[J]. Rock and Soil Mechanics, 2018, 39(S1): 21-32. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2018S1004.htm
[6]	张波, 李术才, 张敦福, 等. 含充填节理岩体相似材料试件单轴压缩试验及断裂损伤研究[J]. 岩土力学, 2012, 33(6): 1647-1652. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201206007.htm ZHANG Bo, LI Shucai, ZHANG Dunfu, et al. Uniaxial compression mechanical property test, fracture and damage analysis of similar material of jointed rock mass with filled cracks[J]. Rock and Soil Mechanics, 2012, 33(6): 1647-1652. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201206007.htm
[7]	蒋爵光, 李隽蓬, 钱惠国, 等. 不同产状节理岩体隧道围岩稳定性的模型试验研究[J]. 水文地质工程地质, 1984, 5(3): 13-19. https://www.cnki.com.cn/Article/CJFDTOTAL-SWDG198405004.htm JIANG Jueguang, LI Junpeng, QIAN Huiguo, et al. Model test study on the stability of tunnel surrounding rock masses of jointed rock masses with different occurrences[J]. Hydrogeology and Engineering Geology, 1984, 5(3): 13-19. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-SWDG198405004.htm
[8]	刘刚, 赵坚, 宋宏伟, 等. 节理密度对围岩变形及破坏影响的试验研究[J]. 岩土工程学报, 2007, 29(11): 1737-1741. doi: 10.3321/j.issn:1000-4548.2007.11.025 LIU Gang, ZHAO Jian, SONG Hongwei, et al. Physical modelling of effect of joint density on deformation and failure of surrounding rocks[J]. Chinese Journal of Geotechnical Engineering, 2007, 29(11): 1737-1741. (in Chinese) doi: 10.3321/j.issn:1000-4548.2007.11.025
[9]	HUANG F, ZHU H, XU Q, et al. The effect of weak interlayer on the failure pattern of rock mass around tunnel-scaled model tests and numerical analysis[J]. Tunneling and Underground Space Technology, 2013, 35: 207-218. doi: 10.1016/j.tust.2012.06.014
[10]	朱维申, 任伟中, 张玉军, 等. 开挖条件下节理围岩锚固效应的模型试验研[J]. 岩土力学, 1997, 18(1): 1-7. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX199701000.htm ZHU Weishen, REN Weizhong, ZHANG Yujun, et al. Model testing research on stability and anchoring effect of opening excavation in jointed rocks[J]. Rock and Soil Mechanics, 1997, 18(1): 1-7. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX199701000.htm
[11]	马腾飞, 李树忱, 李术才, 等. 不同倾角多层节理深部岩体开挖变形破坏规律模型试验研究[J]. 岩土力学, 2016, 37(10): 2899-2908. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201610022.htm MA Tengfei, LI Shuchen, LI Shucai, et al. Model experimental study of deformation and failure law in excavation of deep rock mass with multi-cleftiness of different angles[J]. Rock and Soil Mechanics, 2016, 37(10): 2899-2908. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201610022.htm
[12]	钟志彬, 邓荣贵, 孙怡, 等. 隧道近距大节理硬脆性围岩破裂机理[J]. 中国公路学报, 2018, 31(5): 106-116. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGL201805014.htm ZHONG Zhibin, DENG Ronggui, SUN Yi, et al. Fracture mechanism of hard and brittle tunnel surrounding rock with nearby large joints[J]. China Journal of Highway and Transport, 2018, 31(5): 106-116. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGL201805014.htm
[13]	刘邦, 朱哲明, 周磊, 等. 贯穿隧道的节理对隧道稳定性的影响[J]. 煤炭学报, 2018, 43(5): 1296-1304. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201805012.htm LIU Bang, ZHU Zheming, ZHOU Lei, et al. Effect of the joints crossing the tunnel on tunnel stability[J]. Journal of China Coal Society, 2018, 43(5): 1296-1304. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201805012.htm
[14]	BRUNEAU G, TYLER D B, HADJIGEORGIOU J, et al. Influence of faulting on a mine shaft-a case study: Part I-background and instrumentation[J]. International Journal of Rock Mechanics and Mining Sciences, 2003, 40(1): 95-111. doi: 10.1016/S1365-1609(02)00115-6
[15]	SUORINENI F T, TANNANT D D, KAISER P K. Determination of fault-related sloughage in open stopes[J]. International Journal of Rock Mechanics and Mining Sciences, 1999, 36(7): 891-906. doi: 10.1016/S0148-9062(99)00055-8
[16]	张志强, 李宁, 陈方方, 等. 不同分布距离的软弱夹层对洞室稳定性的影响研究[J]. 岩土力学, 2007, 28(7): 1363-1368. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200707014.htm ZHANG Zhiqiang, LI Ning, CHEN Fangfang, et al. A study of influence of weak interlayers distributed at different intervals on opening stability[J]. Rock and Soil Mechanics, 2007, 28(7): 1363-1368. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200707014.htm
[17]	薄英鋆, 王华宁, 蒋明镜. 深部节理岩体开挖卸荷动力破坏机制的离散元研究[J]. 岩土工程学报 2020, 42(增刊2): 196-201. doi: 10.11779/CJGE2020S2035 BO Yingjun, WANG Huaning, JIANG Mingjing. Cracking mechanism of deep jointed rock mass induced by transient excavation via DEM[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(S2): 196-201. (in Chinese). doi: 10.11779/CJGE2020S2035
[18]	周子涵, 陈忠辉, 王建明, 等. 卸荷条件下岩石平行偏置双裂隙的扩展规律研究[J]. 岩土工程学报, 2020, 42(4): 721-730. doi: 10.11779/CJGE202004015 ZHOU Zihan, CHEN Zhonghui, WANG Jianming, et al. Propagation of offset parallel cracks in rock under unloading conditions[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(4): 721-730. (in Chinese) doi: 10.11779/CJGE202004015
[19]	苏承东, 李冰洋, 吴秋红. 石灰岩孔道试样压缩条件下变形与强度特征的试验研究[J]. 岩石力学与工程学报, 2016, 35(增刊1): 2739-2748. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2016S1018.htm SU Chengdong, LI Bingyang, WU Qiuhong. Experimental study on deformation and strength properties of hollow cylindrical specimens of limestone under compression[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(S1): 2739-2748. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2016S1018.htm
[20]	李夕兵, 吴秋红, 董陇军, 等. 含充填物的厚壁圆筒砂岩的力学性质研究[J]. 岩土工程学报, 2015, 37(12): 2149-2156. doi: 10.11779/CJGE201512003 LI Xibing, WU Qiuhong, DONG Longjun, et al. Mechanical property of thick-walled hollow cylinders of sandstone with filling[J]. Chinese Journal of Geotechnical Engineering, 2015, 37(12): 2149-2156. (in Chinese) doi: 10.11779/CJGE201512003
[21]	侯公羽, 李小瑞, 张振铎, 等. 使用小型围岩试件模拟与再现巷道围岩开挖卸荷过程的试验系统[J]. 岩石力学与工程学报, 2017, 36(9): 2136-2145. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201709007.htm HOU Gongyu, LI Xiaorui, ZHANG Zhenduo, et al. Experimental system for simulating excavation unloading process of rock around roadway by using small cylindrical hollow specimen[J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(9): 2136-2145. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201709007.htm
[22]	侯公羽. 岩石力学基础教程[M]. 北京: 机械工业出版社, 2011. HOU Gongyu. Basic Course of Rock Mechanics[M]. Beijing: China Machine Press, 2011. (in Chinese)
[23]	徐秉业, 刘信声. 应用弹塑性力学[M]. 北京: 清华大学出版社, 1995: 183-203. XU Binye, LIU Xinsheng. Applied Elastic-Plastic Mechanics[M]. Beijing: Tsinghua University Press, 1995: 183-203. (in Chinese)
[24]	李明田, 李术才, 张敦福, 等. 类岩石材料表面裂纹复合型断裂准则探讨[J]. 岩石力学与工程学报, 2011, 30(增刊1): 3326-3333. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2011S1095.htm LI Mingtian, LI Shucai, ZHANG Dunfu, et al. Investigation on mixed-mode fracture criterion of surface crack in rock-like materials[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(S1): 3326-3333. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2011S1095.htm
[25]	董晋鹏, 杨圣奇, 李斌, 等. 共面双裂隙类岩石材料抗拉强度试验研究[J]. 工程力学, 2020, 37(3): 188-201. https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX202003022.htm DONG Jinpeng, YANG Shengqi, LI Bin, et al. Experimental study on the tensile strength of rock-like materials containing two pre-existing coplanar fissures[J]. Engineering Mechanics, 2020, 37(3): 188-201. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GCLX202003022.htm
[26]	余明坤, 许国伟, 唐国栋, 等. 预制裂隙类岩石材料的分步卸载试验研究[J]. 地下空间与工程学报, 2020, 16(6): 1672-1681. https://www.cnki.com.cn/Article/CJFDTOTAL-BASE202006009.htm YU Mingkun, XU Guowei, TANG Guodong, et al. Experimental study on stepped unloading of rock-like materials with pre-existing fissure[J]. Chinese Journal of Underground Space and Engineering, 2020, 16(6): 1672-1681. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BASE202006009.htm
[27]	张波, 郭帅, 杨学英, 等. 含X型裂隙类岩石材料水力裂缝扩展研究[J]. 煤炭学报, 2019, 44(7): 2066-2073. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201907013.htm ZHANG Bo, GUO Shuai, YANG Xueying, et al. Hydraulic fracture propagation of rock-like material with X-type flaws[J]. Journal of China Coal Society, 2019, 44(7): 2066-2073. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201907013.htm
[28]	MARTIN C D. Seventeenth canadian geotechnical colloquium: The effect of cohesion loss and stress path on brittle rock strength[J]. Canadian Geotechnical Journal, 1997, 34(5): 698-725. doi: 10.1139/t97-030

施引文献(3)

期刊类型引用(2)

1.	李东桥，严晓威，梁建文，赵华. 丁字及十字交叉管廊纵向地震响应分析. 天津大学学报(自然科学与工程技术版). 2024(01): 21-30 . 百度学术
2.	梁建文，陈慧芳，李东桥，巴振宁. 地下管廊抗震研究现状综述. 天津大学学报(自然科学与工程技术版). 2024(02): 209-222 . 百度学术