基于双层梁理论的隧道下穿诱发无砟铁路变形计算方法

雷鸣; 张丙强; 刘海; 黄志斌

doi:10.3969/j.issn.0258-2724.20230033

基于双层梁理论的隧道下穿诱发无砟铁路变形计算方法

doi: 10.3969/j.issn.0258-2724.20230033

雷鸣^1,,
张丙强^2, ,,
刘海²,
黄志斌³

1.
长沙学院土木工程学院，湖南长沙 410022
2.
福建理工大学土木工程学院，福建福州 350118
3.
中铁二十四局集团福建铁路建设有限公司，福建福州 350013

基金项目: 福建省自然科学基金（2020J01882，2023J01344）；湖南省教育厅重点项目（19A044）

详细信息

作者简介:
雷鸣（1977—），男，副教授，博士，研究方向为岩土工程，E-mail：540584998@qq.com

通讯作者:
张丙强（1979—），男，教授，博士，研究方向为地下工程，E−mail：zbq@fjut.edu.cn

中图分类号: U211.1
计量
- 文章访问数: 278
- HTML全文浏览量: 173
- PDF下载量: 90
- 被引次数: 1
出版历程
- 收稿日期: 2023-02-07
- 修回日期: 2023-10-01
- 网络出版日期: 2024-01-19
- 刊出日期: 2023-10-31

Theoretical Method for Calculating Rail Deformation of Ballastless Railway Caused by Tunnel Undercrossing Based on Dual Beam Model

1.
School of Civil Engineering, Changsha University, Changsha 410022, China
2.
School of Civil Engineering, Fujian University of Technology, Fuzhou 350118, China
3.
Fujian Railway Construction Co., Ltd., of China Railway 24th Bureau Group, Fuzhou 350013, China

摘要

摘要:
为研究隧道下穿无砟铁路轨道时板底部脱空对轨道变形的影响，提出一种改进的隧道下穿引起无砟铁路变形的计算方法. 首先，将无砟铁路轨道结构简化为双层地基梁模型，建立隧道下穿施工引起无砟铁路变形的控制方程；然后，将无砟铁路分为中间脱空段和两端接地段共三部分，推导隧道下穿施工诱发无砟铁路变形的计算式，对比路基不均匀沉降引起无砟轨道变形的理论计算值与数值模拟结果，验证了理论计算方法的正确性；最后，探讨了新建隧道埋深、隧道下穿施工引起周围地层的损失率，以及既有铁路与隧道间的水平夹角对无砟铁路轨道变形的影响. 研究结果表明：当隧道从铁路下方6 m处垂直穿越既有铁路时，轨道中点变形将达到最大值；隧道施工引起周围地层损失率从0.25%增大到2.50%时，轨道中点变形和轨道板底部脱空区宽度将分别增大4.0倍和2.2倍.
- 轨道变形 /
- 隧道下穿 /
- 双层梁模型 /
- Winkler地基 /
- 理论方法
Abstract:
To analyze the effect of a void under the track slab on the rail deformation of a ballastless railway caused by tunnel undercrossing, an improved calculation method was proposed to predicate the tunnel undercrossing-induced rail deformation of the ballastless railway. Firstly, the rail of the ballastless railway was simplified as a dual subgrade beam model, and the governing equation was established for the rail deformation of ballastless railway caused by tunnel undercrossing. Then, the ballastless railway was divided into three parts including the middle section above a void and the two sections connecting with the subgrade, and the formulas for the rail deformation of the ballastless railway caused by tunnel undercrossing were derived. The proposed theoretical calculation method was verified by comparing the theoretically calculated and numerically simulated results of rail deformation of the ballastless railway caused by uneven settlement of the subgrade. Finally, the influences of parameters on the rail deformation of the ballastless railway were discussed, including the burial depth of the new tunnel, the ground loss rate caused by tunnel undercrossing construction, and the intersection angle between the railway and the tunnel. The results show that when the undercrossing tunnel is 6 m below the railway, the deformation at the midpoint of the rail reaches the maximum. When the ground loss rate caused by tunnel undercrossing construction increases from 0.25% to 2.50%, the deformation at the midpoint of the rail and the width of the void under the track slab increase by 4.0 and 2.2 times, respectively.
- rail deformation /
- tunnel undercrossing /
- dual beam model /
- Winkle foundation /
- theoretical method

HTML全文

双制式列车能实现市域铁路与市内城轨交通的互联互通、贯通运营，减少换乘. 然而，双制式列车存在车体—钢轨回流、车体电位抬升和轴承电蚀等问题，车体回流过大、车体电位过高会对车载弱电设备造成影响，并恶化列车的电磁环境，而轴承的损伤直接影响列车的行车安全. 这些问题与列车整车接地系统的设计密切相关，因此，研究双制式列车接地系统配置具有重要意义.

目前，针对列车接地系统的研究主要集中在单一制式列车方面，如直流（DC）制或交流（AC）制列车. 文献[1]量化了直流制地铁列车轴承电流的大小及其与牵引回流的关系，并提出一种抑制轴承电流的接地方式. 文献[2]建立了动态工况下的地铁列车车网联合仿真模型，对列车接地回流分布进行仿真分析. 文献[3-5]建立高速动车组接地系统的仿真模型，通过优化列车接地系统，降低车体环流和车体—轴端电位.

对于双制式牵引供电系统，为限制不同供电制式系统之间发生牵引回流交换，建议在双制式系统转换段使用绝缘节与特殊馈电装置和开关装置连接^[6]. 针对过渡区段设置方案，文献[7-8]研究了双制式牵引供电系统复合钢轨电位问题，提出应在交、直转换段钢轨处设置绝缘节；文献[9-10]分析了当双制式线路的直、交流区段由同一主所供电时，交流区段产生的谐波问题会影响直流区段；文献[11]建立了双制式列车通过锚段关节式电分相过程的等值电路模型，并分析了列车通过锚段关节式电分相时的暂态特性.

双制式列车接地系统需兼容交、直牵引回流. 目前，对于双制式列车接地系统的研究较少. 文献[12]介绍了双制式列车牵引传动系统. 文献[13]提出一种双制式列车接地系统方案设计，提出在各节车厢车体与汇流排之间均应设置保护电阻.

目前，对列车接地系统的研究主要以仿真建模为主，缺乏统一的数学模型，这限制了列车接地系统在各区段运行过程性能的深入研究. 为此，本文提出列车接地系统的车—轨链式等效电路模型，并结合牵引供电系统进行车-地一体化联合计算，分析双制式列车运行在直流区段、交流区段的车体环流分布，并对列车接地系统提出改进建议. 该方法也适用于单供电制式或多供电制式列车接地系统的分析.

1. 双制式列车接地系统

双制式列车可在交、直流区段贯通运行. 交、直流区段之间需通过交—直转换段进行分离，其结构如图1所示.

图 1 双制式牵引供电系统

Figure 1. Dual-system traction power supply system

下载: 全尺寸图片幻灯片

为适应直流和交流两种供电制式，双制式列车除了和普通地铁列车一样装有牵引逆变器以外，也装有与牵引变压器和牵引变流器. 列车在交—直转换段完成不同制式的切换.

国内某型双制式列车为6节编组，其接地系统如图2所示. 列车牵引传动和接地系统均首尾对称. 在直流牵引工况下，列车工作接地与保护接地共用同一接地汇流排. 在交流牵引工况下，列车工作接地设置专用汇流排与保护接地的汇流排分开.

图 2 线路列车既有接地系统

Figure 2. Existing train grounding system

下载: 全尺寸图片幻灯片

1.1 车—轨链式电路

列车接地系统可以分为车体、接地汇流排和钢轨三层结构，如图3所示. 图中：Z_ctlj和Z_ct分别为列车车间连接线阻抗和列车车体阻抗，Z_hlp为汇流排间连接线阻抗，R₀为车体与接地汇流排之间的保护电阻，Z_jd1和Z_jd2为接地汇流排与轮对之间连接线阻抗，Z_r1为同一转向架轮对间的钢轨阻抗，Z_r2为不同一转向架间的钢轨阻抗，Z_r3为不同车厢间的钢轨阻抗. 接地汇流排切面包含4个节点（图3中 ①～④）.

图 3 列车接地系统电路模型

Figure 3. Circuit model of train grounding system

下载: 全尺寸图片幻灯片

单个接地汇流排的节点导纳矩阵 ${{\boldsymbol{Y}}_{{\text{cb}}}}$ 如式（1）所示.

$\qquad\qquad\qquad\quad {{\boldsymbol{Y}}_{{\text{cb}}}} = \left( {\begin{array}{*{20}{c}} {1/{R_0}}&{ - 1/{R_0}}&0&0 \\ { - 1/{R_0}}&{1/{R_0} + 1/{Z_{{\text{jd}}2}} + 1/{Z_{{\text{jd}}1}}}&{ - 1/{Z_{{\text{jd}}2}}}&{ - 1/{Z_{{\text{jd}}1}}} \\ 0&{ - 1/{Z_{{\text{jd}}1}}}&{1/{Z_{{\text{jd}}1}} + 1/{Z_{{\text{r}}1}}}&{ - 1/{Z_{{\text{r}}1}}} \\ 0&{ - 1/{Z_{{\text{jd}}2}}}&{ - 1/{Z_{{\text{r}}1}}}&{1/{Z_{{\text{jd}}2}} + 1/{Z_{{\text{r}}1}}} \end{array}} \right).$

(1)

各接地汇流排之间还存在车间连接线、汇流排间连接线等横连线，其中，车间连接线的横连线导纳矩阵Y_lj如式（2）所示.

${{\boldsymbol{Y}}_{{\text{lj}}}} = \left( {\begin{array}{*{20}{c}} {1/{Z_{{\text{ctlj}}}}}&{}&{ - 1/{Z_{{\text{ctlj}}}}}&{} \\ {}&{\boldsymbol{O}}&{}&{\boldsymbol{O}} \\ { - 1/{Z_{{\text{ctlj}}}}}&{}&{1/{Z_{{\text{ctlj}}}}}&{} \\ {}&{\boldsymbol{O}}&{}&{\boldsymbol{O}} \end{array}} \right),$

(2)

式中：O为3阶零矩阵.

由于列车接地汇流排的布置具有对称性，可以根据列车接地汇流排划分切面，形成车—轨链式等效电路，如图4所示. 图中，Y_cbi为第i个接地汇流排节点导纳矩阵，Z_hli为横连线阻抗矩阵，I_cti为注入第i个接地汇流排的牵引电流向量， $i = 1,2, \cdots ,m$ .

图 4 车—轨链式电路

Figure 4. Train–rail chain circuit

下载: 全尺寸图片幻灯片

形成的列车车—轨导纳矩阵为

${{\boldsymbol{Y}}_{{\text{CT}}}} = \left( {\begin{array}{*{20}{c}} {{\boldsymbol{Z}}_{{\text{hl1}}}^{ - 1} + {{\boldsymbol{Y}}_{{\text{cb1}}}}}&{ - {\boldsymbol{Z}}_{{\text{hl1}}}^{ - 1}}&{}&{}&{}&{} \\ { - {\boldsymbol{Z}}_{{\text{hl1}}}^{ - 1}}&{{\boldsymbol{Z}}_{{\text{hl1}}}^{ - 1} + {{\boldsymbol{Y}}_{{\text{cb2}}}} + {\boldsymbol{Z}}_{{\text{hl2}}}^{ - 1}}&{ - {\boldsymbol{Z}}_{{\text{hl2}}}^{ - 1}}&{}&{}&{} \\ {}&{ - {\boldsymbol{Z}}_{{\text{hl2}}}^{ - 1}}&{{\boldsymbol{Z}}_{{\text{hl2}}}^{ - 1} + {{\boldsymbol{Y}}_{{\text{cb3}}}} + {\boldsymbol{Z}}_{{\text{hl3}}}^{ - 1}}&{ - {\boldsymbol{Z}}_{{\text{hl3}}}^{ - 1}}&{}&{} \\ {}&{}& & \vdots & &{} \\ {}&{}&{}&{ - {\boldsymbol{Z}}_{{\text{hl}}m{{-2}}}^{ - 1}}&{{\boldsymbol{Z}}_{{\text{hl}}m{{-2}}}^{ - 1} + {{\boldsymbol{Y}}_{{\text{cb}}m{{-}}1}} + {\boldsymbol{Z}}_{{\text{hl}}m{{-1}}}^{ - 1}}&{ - {\boldsymbol{Z}}_{{\text{hl}}m{{-1}}}^{ - 1}} \\ {}&{}&{}&{}&{ - {\boldsymbol{Z}}_{{\text{hl}}m{{-1}}}^{ - 1}}&{{{\boldsymbol{Y}}_{{\text{cb}}m}} + {\boldsymbol{Z}}_{{\text{hl}}m{{-1}}}^{ - 1}} \end{array}} \right).$

(3)

对于牵引电流向量，需分别考虑列车处于直流和交流两种工况. 含有牵引电流注入的车—轨切面电流向量如式（4）所示.

$\begin{gathered} {{\boldsymbol{I}}_{{\text{ct}}i}}(t) = \left\{ {\begin{array}{*{20}{c}} {{{\left[ {\begin{array}{*{20}{c}} 0&{{{\dot I}_{{\text{tract}}i}}(t)}&0&0 \end{array}} \right]}^{\mathrm{T}}}\begin{array}{*{20}{c}} , \;\; &{{T_{\mathrm{r}}}(t) = {T_{{\mathrm{r}}0}}}, \end{array}} \\ {{{\left[ {\begin{array}{*{20}{c}} 0&0&{{{\dot I}_{{\text{tract}}i}}(t)}&0 \end{array}} \right]}^{\mathrm{T}}}\begin{array}{*{20}{c}} , \;\; &{{T_{\mathrm{r}}}(t) = {T_{{\mathrm{r}}1}}}, \end{array}} \end{array}} \right. \end{gathered}$

(4)

式中：T_r(t)为时刻t列车的供电制式状态；T_r0为列车处于直流制牵引供电；T_r1为列车处于交流制牵引供电； $\dot{I}_{{\mathrm{tract}}i}(t)$ 为时刻t注入第i个接地汇流排的牵引电流，直流制式下 $\dot{I}_{{\mathrm{tract}}i}(t)$ 仅取实部.

对列车车—轨首、末切面的电流向量进行修正，修正后的车—轨首、末切面的电流向量I_ct1new、I_ctmnew如式（5）所示.

$\begin{gathered} \left\{ {\begin{array}{*{20}{c}} {{{\boldsymbol{I}}_{{\text{ct1new}}}} = {{\boldsymbol{I}}_{{\text{ct1}}}} + {{\left[ {\begin{array}{*{20}{c}} 0&0&{ - {{\dot I}_{{\text{GL}}}}(t)}&0 \end{array}} \right]}^{\mathrm{T}}}}, \\ {{{\boldsymbol{I}}_{{\text{ct}}{{m}}{\text{new}}}} = {{\boldsymbol{I}}_{{\text{ct}}{{m}}}} + {{\left[ {\begin{array}{*{20}{c}} 0&0&0&{ - {{\dot I}_{{\text{GR}}}}(t)} \end{array}} \right]}^{\mathrm{T}}}}, \end{array}} \right. \\ \end{gathered}$

(5)

$\sum\limits_{i = 1}^m {{{\dot I}_{{\text{tract}}i}}(t) = {{\dot I}_{{\text{tr}}}}(t) = {{\dot I}_{{\text{GL}}}}(t) + {{\dot I}_{{\text{GR}}}}(t)},$

(6)

式中： $\dot{I}_{\mathrm{tr}}(t)$ 为列车取流， $\dot{I}_{\mathrm{GL}}(t)$ 为车—轨首切面钢轨电流， $\dot{I}_{\mathrm{GR}}(t)$ 为车—轨末切面钢轨电流. 直流制式下 $\dot{I}_{\mathrm{GL}}(t)$ 、 $\dot{I}_{\mathrm{GR}}(t)$ 仅取实部.

列车车—轨电流向量为

${{\boldsymbol{I}}_{{\text{CT}}}} = {({{\boldsymbol{I}}_{{\text{ct1new}}}},{{\boldsymbol{I}}_{{\text{ct2}}}},{{\boldsymbol{I}}_{{\text{ct3}}}}, \cdots ,{{\boldsymbol{I}}_{{\text{ct}}(m - 2)}},{{\boldsymbol{I}}_{{\text{ct}}(m - 1)}},{{\boldsymbol{I}}_{{\text{ct}}m{\text{new}}}})^{\mathrm{T}}}.$

(7)

形成的车—轨链式等效电路节点电压向量为

${{\boldsymbol{U}}_{{\text{CT}}}} = {\boldsymbol{Y}}_{{\text{CT}}}^{ - 1}{{\boldsymbol{I}}_{{\text{CT}}}}.$

(8)

2. 双制式牵引供电系统建模

2.1 交、直流区段牵引供电系统建模

交、直流区段的牵引网均满足平行多导体结构，故可采用链式电路模型表示. 交、直流区段l长度的牵引网阻抗矩阵 ${{\boldsymbol{Z}}_l}$ 、导纳矩阵 ${{\boldsymbol{Y}}_l}$ 如式（9）所示^[14].

$\left\{ {\begin{array}{*{20}{l}} {{\boldsymbol{Z}}_l} = \sinh (\sqrt {{\boldsymbol{ZY}}} l){{({\boldsymbol{ZY}})}^{ - \frac{1}{2}}}{\boldsymbol{Z}} =\\ \quad {\boldsymbol{Z}}{\displaystyle\sum\limits_{n = 0}^\infty {\dfrac{{{l^{2n + 1}}}}{{(2n + 1)!}}{{({\boldsymbol{ZY}})}^n}} }, \\ {{\boldsymbol{Y}}_l}/2 = {\boldsymbol{Z}}{{({\boldsymbol{ZY}})}^{ - \frac{1}{2}}}\tanh (\sqrt {{\boldsymbol{ZY}}} l/2) = \\ \quad {{\boldsymbol{Z}}_l}^{ - 1}\displaystyle\sum\limits_{n = 1}^\infty {\dfrac{{{l^{2n}}}}{{(2n)!}}{{({\boldsymbol{ZY}})}^n}}, \end{array}} \right.$

(9)

式中：Z、Y分别为牵引网单位阻抗矩阵、导纳矩阵.

在直流区段中，只考虑Z、Y的电阻部分.

牵引变电所可采用诺顿等效电路模型. 交流牵引变电所α、β供电臂的等效电流源 $\dot{I}_{{\mathrm{equ}}\alpha}$ 、 $\dot{I}_{{\mathrm{equ}}\beta }$ 可用式（10）求得^[15].

$\left\{ {\begin{array}{*{20}{l}} {\left[ {\begin{array}{*{20}{c}} {{{\dot{I}}_{{\text{equ}}\alpha }}}&{} \\ {}&{{{\dot{I}}_{{\text{equ}}\beta }}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{Z_\alpha }} \\ {{Z_\beta }} \end{array}} \right] = {\boldsymbol{N}}{{\boldsymbol{E}}_{{\mathrm{ABC}}}} - {Z_\gamma }\left[ {\begin{array}{*{20}{c}} {{{\dot{I}}_\beta }} \\ {{{\dot{I}}_\alpha }} \end{array}} \right], } \\ {{{\boldsymbol{Z}}_{\alpha \beta }} = \left[ {\begin{array}{*{20}{c}} {{ Z_\alpha }}&{{ Z_\gamma }} \\ {{ Z_\gamma }}&{{ Z_\beta }} \end{array}} \right],} \end{array}} \right.$

(10)

式中：N为电压变化矩阵，E_ABC为牵引变压器一次侧电压向量，Z_α、Z_β、 $Z_\gamma$ 为变电所T型等效电路等值阻抗，Z_αβ为变电所等效阻抗矩阵， $\dot{I}_\alpha$ 、 $\dot{I}_\beta$ 分别为变电所α、β供电臂端口电流.

直流牵引变电所使用24脉波整流机组（24 MPR），其等效电流源I_d和端口等效电阻r₀为^[16]

$\left\{ {\begin{array}{*{20}{l}} {{I_{\text{d}}} = \dfrac{{{U_{\text{d}}}}}{{{r_0}}} = \dfrac{1}{{{r_0}}}\dfrac{{24\sqrt 2 }}{{\text{π}} }\sin \dfrac{{\text{π}} }{{24}}{U_2} = 1.41\dfrac{{{U_2}}}{{{r_0}}}}, \\ {{r_0} = \dfrac{3}{{2{\text{π}} }}{X_{\text{c}}}}, \end{array}} \right.$

(11)

式中：U_d为24 MPR等效电压源，U₂为24 MPR二次侧线电压额定值，X_c为24 MPR换流阻抗.

2.2 交、直流区段联合牵引供电计算

在双制式线路上存在不同供电制式的列车运行，故需对交、直流区段进行联合牵引供电计算.

对于牵引网的建模，若线路均为高架区段，交流区段采用直供带回流方式，其牵引网由回流线、接触网、钢轨、综合地线构成，单位牵引网阻抗、导纳为4阶矩阵；直流区段牵引网由接触网、钢轨、排流网构成，单位牵引网阻抗、导纳为3阶矩阵. 为形成统一的牵引网阻抗、导纳矩阵，直流区段需添加虚拟导线使其单位牵引网阻抗、导纳变为4阶矩阵. 为使得系统节点导纳矩阵非奇异，虚拟导线与钢轨之间可设置大电阻横向连接. 双制式线路牵引网模型如所示. 图中， $\dot{I}_{\mathrm{trdc}}$ 、 $\dot{I}_{\mathrm{trac}}$ 分别为列车在直流区段、交流区段取流.

图 5 双制式牵引网模型

Figure 5. Dual-system traction network model

下载: 全尺寸图片幻灯片

根据叠加定理，交、直流区段统一牵引供电计算中节点电压方程组为

$\left\{ {\begin{array}{*{20}{c}} {{{\boldsymbol{U}}_{{\text{Sdc}}}} = {\boldsymbol{G}}_{{\text{S}}j(a + b) \times j(a + b)}^{ - 1}{{\boldsymbol{I}}_{{\mathrm{dcs}}}}}, \\ {{{\boldsymbol{U}}_{{\text{Sac}}}} = {\boldsymbol{Y}}_{{\text{S}}j(a + b) \times j(a + b)}^{ - 1}{{\boldsymbol{I}}_{{\mathrm{acs}}}}}, \end{array}} \right.$

(12)

式中：U_Sdc为系统节点电压向量直流分量，U_Sac为系统节点电压向量交流分量，I_dcs为直流电流向量，I_acs为交流电流向量，G_S为仅直流源作用时系统牵引网节点导纳矩阵，Y_S为仅交流源作用时系统牵引网节点导纳矩阵，j为单位牵引网阻抗导纳矩阵阶数，a、b分别为直流区段、交流区段切面数.

3. 车-地一体的交、直流区段联合牵引供电计算

由双制式牵引系统供电计算可以得到各时刻列车的取流信息以及列车位置相邻切面的钢轨电流信息. 根据这些数据，结合车—轨链式等效电路模型，可以得出列车在各时刻的车体环流分布以及各节车厢车体—轴端电位.

车-地一体的交、直流区段联合牵引供电计算电路模型如图6所示，图6（a）中列车接地系统等值电路见图6（b）. 图中：Z_c为接触网阻抗，Z_r为钢轨阻抗，Z_p为排流网阻抗，Z_fb为回流线阻抗， Z_r-p为钢轨对排流网阻抗，Z_p-e为排流网对地阻抗，Z_gw为贯通地线阻抗，Z_e为贯通地线接地阻抗.

图 6 车—地一体的交、直流区段联合牵引供电计算电路模型

Figure 6. Circuit model of united traction power supply calculation for vehicle-ground integration in AC and DC sections

下载: 全尺寸图片幻灯片

该联合计算流程如图7所示. 该方法不仅适用于双流制接地系统的研究，也可用于单一供电制式或者其他多供电制式列车接地系统的分析.

图 7 供电系统与列车接地系统联合计算流程

Figure 7. Flowchart of united calculation of power supply system and train grounding system

下载: 全尺寸图片幻灯片

4. 列车接地保护电阻优化模型

4.1 目标函数

从车体—轴端电位角度，建立列车接地保护电阻优化模型，尽可能地抑制车体—轴端电位大小.

选择列车运行过程中，列车各车车体—轴端电位最大值作为目标函数f₁(X)，如式（13）所示.

$\left\{ {\begin{array}{*{20}{c}} {\min {f_1}(X) = {U_{{\text{CZ,}}k}}}, \\ {X = \{ {x_k}\} ,\quad {x_k} \in \varepsilon } , \end{array}} \right.$

(13)

式中：U_CZ,k为列车第k节车厢车体—轴端电位最大值，X为列车各车保护电阻值配置集合，x_k为列车第k节车厢的保护电阻值，ε为以10 mΩ为步长的随机离散变量组合.

4.2 约束条件

求解单目标函数式（13），其约束条件为列车保护电阻上、下限值以及通过保护电阻的车体电流上限值，如式（14）所示.

$\left\{ {\begin{array}{*{20}{c}} {{x_{\min }} \leqslant {x_k} \leqslant {x_{\max }}}, \\ {I_{{\text{CT}},k} \leqslant {I_{{\text{CT}},\max }}}, \end{array}} \right.$

(14)

式中：x_max、x_min分别为保护电阻设置的上、下限，I_CT,k为列车第k节车厢保护电阻的电流最大值，I_CT,max为列车保护电阻的电流上限值.

5. 算例分析

为验证算法的有效性，以国内某条双制式轨道交通线路（图8）为例，分析列车处于直流区段和交流区段2种情况下的车体环流分布，并与实测结果进行对比.

图 8 双制式轨道交通线路

Figure 8. Dual-system rail transit line

下载: 全尺寸图片幻灯片

5.1 仿真条件

该线路全长28.22 km，各牵混所以及主所位置见表1. 全线共设置2个直流车站和5个交流车站，所有车站均为高架站. T1和T2为区间牵引所，TS1为带牵引所的直流区段车站，S2～S7为普通车站，K为线路设计起点，TPS1和TPS2为交流区段牵引主变电所. 交—直转换段NS设置在S2和S3站点之间. K—NS为直流区段，采用1500 V柔性接触网授流，通过钢轨回流；采用地面制动电阻吸收列车再生制动能量；TS1处的钢轨电位限制器（OVPD）接地. NS—S7为交流区段，采用25 kV柔性接触网授流，通过回流线、贯通地线以及钢轨回流. 定义K—S7为上行方向， S7—K为下行方向.

表 1 牵引所以及主所位置

Table 1. Location of traction station and main substation km

牵混所编号	位置	主所编号	位置
T1	0.15	TPS1	19.81
TS1	1.74	TPS2	28.21
T2	3.70

下载: 导出CSV

| 显示表格

该型列车车辆参数见表2所示，其接地系统如图2所示. 定义其中1车为头车，方向朝S7；6车为尾车，方向朝TS1. 全天采用K—S7单一运行路线，发车间隔设置为500 s.

表 2 车辆参数

Table 2. Train parameters

参数	数值
列车编组	6A
车重/t	223.88
结构速度/（km·h⁻¹）	120
最大加速度/（m·s⁻²）	1.1 （DC），0.9 （AC）
最大减速度/（m·s⁻²）	1.2

下载: 导出CSV

| 显示表格

列车接地系统中各连接线阻抗以及车体阻抗由现场实际测量列车得出，对应参数见表3.

表 3 列车接地系统模型参数

Table 3. Model parameters of train grounding system

参数	AC	DC
Z_ct/mΩ	1.500	0.027
Z_ctlj/mΩ	0.82	0.82
R₀/mΩ	50.6	50.0
Z_jd1/mΩ	0.89	0.89
Z_jd2/mΩ	1.21	1.21
Z_hlp/mΩ	3.9	3.9
Z_r/（mΩ·km⁻¹）	191.0^[17]	36.4^[18]

下载: 导出CSV

| 显示表格

牵引计算模拟得到的上下行列车功率-时间曲线如图9所示.

图 9 上、下行列车功率曲线

Figure 9. Power curves of up and down trains

下载: 全尺寸图片幻灯片

5.2 实测验证

对该线路某辆列车的一组动力单元的接地回流系统的电流进行同步检测，在列车保护电阻处加装电流传感器，规定流出车体电流方向为正. 监测点的详细布置方案如图10所示，图中，k-p为第k节车厢第p个接地汇流排处的保护电阻，余图同. 不同位置设备的检测信号通过GPS同步授时.

图 10 监测点布置方案

Figure 10. Monitoring point layout

下载: 全尺寸图片幻灯片

5.2.1 直流区段

在直流工况下，列车牵引电流通过各车接地汇流排分散注入至钢轨. 选取列车在上行方向TS1—K区间运行过程，计算并分析列车车体环流分布情况. 该区间列车取流以及左、右切面钢轨电流计算值如图11所示. 列车一组动力单元的车体环流分布计算值和实测值如图12所示.

图 11 列车取流以及钢轨电流（直流）

Figure 11. Train’s current and rail current (DC)

下载: 全尺寸图片幻灯片

在该直流区间中，列车各保护电阻的电流计算值如图13所示. 由图可知，通过列车保护电阻电流的最大值为22.3 A. 该时刻列车处于牵引工况，车体环流路径为4～6车流入车体，1～3车流出车体.

图 12 保护电阻电流计算与实测对比（直流）

Figure 12. Comparison of protection resistance current calculation and measured value (DC)

下载: 全尺寸图片幻灯片

图 13 列车各保护电阻电流（直流）

Figure 13. Protection resistance current of train (DC)

下载: 全尺寸图片幻灯片

5.2.2 交流区段

列车在交流工况下，牵引电流通过3、4车牵引变压器一次侧专用汇流排注入至钢轨. 选取列车在上行方向S5—S4区间运行过程计算列车取流以及左、右切面钢轨电流计算值，如图14所示. 列车一组动力单元的车体环流分布计算值和实测值如图15所示.

图 14 列车取流以及钢轨电流（交流）

Figure 14. Train’s current and rail current (AC)

下载: 全尺寸图片幻灯片

在该交流区间中，列车各保护电阻的电流计算值如图16所示. 通过列车保护电阻电流的最大峰值为9.6A，该时刻列车处于牵引工况，车体环流路径为3～6车流入车体，1～2车流出车体.

图 16 列车各保护电阻电流（交流）

Figure 16. Protection resistance current of train (AC)

下载: 全尺寸图片幻灯片

计算与实测存在误差的原因是由于列车车体环流分布受列车司机操纵设置以及钢轨均流线等横联线位置的影响. 列车1～3车监测位置的保护电阻电流实测值与计算值的总体变化趋势大致相同，说明该车-地一体联合牵引供电计算模型可用于分析列车车体环流分布，验证了列车在交、直流区段时模型的准确性.

图 15 保护电阻电流计算与实测对比（交流）

Figure 15. Comparison of protection resistance current calculation and measured value (AC)

下载: 全尺寸图片幻灯片

5.3 对列车接地系统配置的讨论与分析

在正常工况下，列车接地系统配置需能抑制车体环流以及车体—轴端电位. 对于交流制动车组，车体—轴端电位的最大值不应超过1 V^[19]. 地铁列车的车体—轴端电位和车体环流没有明确的限值. 关于列车保护电阻，TB/T 2977—2016^[20]规定，列车车体与轨道之间的保护电阻值不应大于50 mΩ^[20].

由图12、15可见，在交、直流区段中该型双制式列车的车体环流分布总体上是头尾车（1、6车）最大，中间车（3、4车）最小. 为使列车在交、直流区段运行过程中降低头尾车车体—轴端电位的同时，尽量抑制车体环流，本文对既有的列车接地系统方案提出2种改进方案：方案1为更改1、6车的保护电阻值，方案2为在1、6车增设1处接地汇流排，并更改1、6车的保护电阻值. 改进方案如图17所示.

图 17 列车接地系统改进方案

Figure 17. Improved scheme of train grounding system

下载: 全尺寸图片幻灯片

设置x_min=10 mΩ，x_max=50 mΩ，I_CT,max=60 A，使用2种改进方案重新计算，可得其在交、直流区段运行中的车体电流和车体—轴端电位分布.

方案1的车体—轴端电位以及车体环流分布最大值如图18所示. 当1、6车的保护电阻从50 mΩ降为10 mΩ，列车在直流区段的车体—轴端电位最大值从1.25 V降至1.04 V，而车体电流最大值从22.33 A增至88.01 A. 列车在交流区段的车体—轴端电位最大值从0.89 V降至0.67 V，而车体电流最大值从16.86 A增至45.98 A.

图 18 列车的车体—轴端电位和车体环流分布（方案1）

Figure 18. Distribution of current circulation and potential of train’s body–axle end (scheme 1)

下载: 全尺寸图片幻灯片

方案2的车体—轴端电位以及车体环流分布最大值如图19所示. 在1、6车各增设1处汇流排，1、6车的保护电阻从50 mΩ降为10 mΩ，列车在直流区段的车体—轴端电位最大值从2.12 V降至1.64 V，而车体电流最大值从28.99 A增至107.97 A. 列车在交流区段的车体—轴端电位最大值从0.79 V降至0.50 V，而车体电流最大值从10.95 A增至25.11 A.

图 19 列车的车体—轴端分布和车体环流分布（方案2）

Figure 19. Distribution of current circulation and potential of train’s body–axle end (scheme 2)

下载: 全尺寸图片幻灯片

方案1相较方案2，列车的车体—轴端电位最大值降低36.58%～41.04%，车体环流最大值降低18.49%～22.97%；并且在方案1中，1、6车采用20 mΩ保护电阻接地可使列车的车体—轴端电位最大值为1.15 V，车体电流最大值为50.30 A，达到抑制车体—轴端电位的最优效果.

6. 结　论

1）建立了车地一体化的交、直流区段联合牵引供电计算模型，并将模型计算值与实测数据进行比较，验证模型的正确性.

2）列车车体环流除受自身牵引取流影响外，也受邻车位置和工况、牵引网横联线设置的影响. 在交、直流区段中，列车头尾车的车体电流和车体—轴端电位值最大，中间车的车体环流和车体—轴端电位值最小.

3）建立了列车接地保护电阻优化模型，并以既有的双制式列车接地系统为例，讨论分析了2种改进的接地系统方案. 方案1中头尾车采用20 mΩ保护电阻接地，能达到抑制车体—轴端电位的最优效果.

本文仅对稳态情况下列车接地保护电阻的设置进行了探讨，车体浪涌的暂态模型还需进行进一步研究.

图 1 隧道下穿既有铁路力学分析示意

Figure 1. Mechanical analysis of existing railway caused by tunnel undercrossing

下载: 全尺寸图片幻灯片

图 2 本文结果与文献[23]数据的对比

Figure 2. Comparisons between results of this study and data from literature [23]

下载: 全尺寸图片幻灯片

图 3 隧道埋深的影响

Figure 3. Influence of buried depth of new tunnel

下载: 全尺寸图片幻灯片

图 4 隧道施工引起地层损失率的影响

Figure 4. Influence of ground loss rate caused by tunnel undercrossing construction

下载: 全尺寸图片幻灯片

图 5 地铁线路与隧道水平夹角的影响

Figure 5. Influence of intersection angle between between railway and tunnel

下载: 全尺寸图片幻灯片

表 1 无砟轨道材料参数

Table 1. Parameters of ballastless rail material

结构弹性模量/Pa 泊松比说明

钢轨 2.10 × 10¹¹ 0.30 T60 轨
道床板 3.25 × 10¹⁰ 0.17 C40 混凝土
支承层 2.55 × 10¹⁰ 0.17 C20 混凝土
路基 180 0.25 密度 300 kg/m³

下载: 导出CSV

参考文献(24)

[1]	刘建友,赵勇,吕刚,等. 隧道下穿高速铁路路基风险分级及风险评价方法研究[J]. 现代隧道技术,2020,57(6): 8-16,54. LIU Jianyou, ZHAO Yong, LV Gang, et al. Research on risk classification and assessment method for tunnels crossing under high-speed railway subgrade[J]. Modern Tunnelling Technology, 2020, 57(6): 8-16,54.
[2]	仇文革,万世付,高刚刚,等. 砂卵石地层盾构隧道下穿铁路咽喉区地表沉降控制研究[J]. 现代隧道技术,2021,58(5): 37-45. QIU Wenge, WAN Shifu, GAO Ganggang, et al. Study on the control measures against ground settlement induced by shield tunnel construction underneath railway throat in sandy cobble stratum[J]. Modern Tunnelling Technology, 2021, 58(5): 37-45.
[3]	吕培林,周顺华. 软土地区盾构隧道下穿铁路干线引起的线路沉降规律分析[J]. 中国铁道科学,2007,28(2): 12-16. doi: 10.3321/j.issn:1001-4632.2007.02.003 LV Peilin, ZHOU Shunhua. Analysis on upper rail settlement in soft ground resulting from shield tunnelling across main railway line[J]. China Railway Science, 2007, 28(2): 12-16. doi: 10.3321/j.issn:1001-4632.2007.02.003
[4]	孙连勇,黄永亮,王启民,等. 地铁盾构隧道下穿既有铁路变形控制研究[J]. 现代隧道技术,2018,55(5): 140-145,173. SUN Lianyong, HUANG Yongliang, WANG Qimin, et al. Study of controlling measures for existing railway deformation caused by metro shield tunnel under-passing construction[J]. Modern Tunnelling Technology, 2018, 55(5): 140-145,173.
[5]	张乾,蔡小培,钟阳龙,等. 无砟轨道路基不均匀沉降区高速列车动力特性[J]. 西南交通大学学报,2023,58(1): 133-140. ZHANG Qian, CAI Xiaopei, ZHONG Yanglong, et al. Dynamic characteristics of high-speed trains in differential subgrade settlement zone of ballastless track[J]. Journal of Southwest Jiaotong University, 2023, 58(1): 133-140.
[6]	唐钱龙,李永恒,彭立敏,等. 盾构下穿施工变形条件下既有铁路行车限速标准研究[J]. 铁道科学与工程学报,2022,19(9): 2700-2708. TANG Qianlong, LI Yongheng, PENG Limin, et al. Study on speed limit standard of existing railway traffic under the deformation condition of shield tunneling construction[J]. Journal of Railway Science and Engineering, 2022, 19(9): 2700-2708.
[7]	KLAR A, VORSTER T E, SOGA K, et al. Elastoplastic solution for soil-pipe-tunnel interaction[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2007, 133(7): 782-792. doi: 10.1061/(ASCE)1090-0241(2007)133:7(782)
[8]	ZHANG Z G, HUANG M S, ZHANG M X. Deformation analysis of tunnel excavation below existing pipelines in multi-layered soils based on displacement controlled coupling numerical method[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2012, 36(11): 1440-1460. doi: 10.1002/nag.2098
[9]	杨成永,马文辉,费腾,等. 考虑轴力和剪切效应的盾构隧道纵向变形分析[J]. 西南交通大学学报,2022,57(1): 139-147. YANG Chengyong, MA Wenhui, FEI Teng, et al. Analysis of longitudinal deformation of shield tunnel structures with consideration of axial force and shear effect[J]. Journal of Southwest Jiaotong University, 2022, 57(1): 139-147.
[10]	LI P, DU S J, WANG Y H, et al. Timoshenko beam solution for the response of existing tunnels because of tunneling underneath[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2016, 40(5): 766-784. doi: 10.1002/nag.2426
[11]	KLAR A, VORSTER T E B, SOGA K, et al. Soil-pipe interaction due to tunnelling: comparison between Winkler and elastic continuum solutions[J]. Géotechnique, 2005, 55(6): 461-466.
[12]	LIU X, FANG Q, ZHANG D L, et al. Behaviour of existing tunnel due to new tunnel construction below[J]. Computers and Geotechnics, 2019, 110: 71-81. doi: 10.1016/j.compgeo.2019.02.013
[13]	LIN C G, HUANG M S, NADIM F, et al. Tunnelling-induced response of buried pipelines and their effects on ground settlements[J]. Tunnelling and Underground Space Technology, 2020, 96: 103193.1-103193.17. doi: 10.1016/j.tust.2019.103193
[14]	LIU X, JIANG A N, HAI L, et al. Study on soil gap formation beneath existing underground structures due to new excavation below[J]. Computers and Geotechnics, 2021, 139: 104379.1-104379.12. doi: 10.1016/j.compgeo.2021.104379
[15]	可文海,管凌霄,刘东海,等. 盾构隧道下穿管道施工引起的管-土相互作用研究[J]. 岩土力学,2020,41(1): 221-228,234. KE Wenhai, GUAN Lingxiao, LIU Donghai, et al. Research on upper pipeline-soil interaction induced by shield tunnelling[J]. Rock and Soil Mechanics, 2020, 41(1): 221-228,234.
[16]	冯国辉,徐长节,郑茗旺,等. 新建隧道下穿既有隧道引起的隧-土相互作用研究[J]. 工程力学,2023,40(5): 59-68. FENG Guohui, XU Changjie, TEY Mingwang, et al. Study of tunnel-soil interaction induced by tunneling underlying[J]. Engineering Mechanics, 2023, 40(5): 59-68.
[17]	黄晓林,周志刚,刘齐建. 浅埋隧道开挖引起的高铁轨道变形与内力分析[J]. 湖南科技大学学报(自然科学版),2019,34(1): 42-48. HUANG Xiaolin, ZHOU Zhigang, LIU Qijian. Study on internal forces and deformation of the high-speed rail by the excavation of the shallow tunnel[J]. Journal of Hunan University of Science & Technology (Natural Science Edition), 2019, 34(1): 42-48.
[18]	DECK O, SINGH A. Analytical model for the prediction of building deflections induced by ground movements[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2012, 36(1): 62-84. doi: 10.1002/nag.993
[19]	邹春华,周顺华,王炳龙,等. 路基不均匀沉降对有砟轨道沉降影响的模型试验[J]. 同济大学学报(自然科学版),2011,39(6): 862-869. doi: 10.3969/j.issn.0253-374x.2011.06.014 ZOU Chunhua, ZHOU Shunhua, WANG Binglong, et al. Model test study of influence of differential subgrade settlement on ballasted track settlement[J]. Journal of Tongji University (Natural Science), 2011, 39(6): 862-869. doi: 10.3969/j.issn.0253-374x.2011.06.014
[20]	付龙龙,宫全美,周顺华,等. 列车荷载作用下有砟轨道轨面沉降与路基不均匀沉降间的相关关系[J]. 振动与冲击,2013,32(14): 23-28,39. FU Longlong, GONG Quanmei, ZHOU Shunhua, et al. Relationships between ballasted track settlement and differential subgrade settlement under train load[J]. Journal of Vibration and Shock, 2013, 32(14): 23-28,39.
[21]	邹春华,周顺华,王长丹,等. 路基不均匀沉降引起有砟轨道沉降的计算方法[J]. 同济大学学报(自然科学版),2013,41(8): 1237-1242. ZOU Chunhua, ZHOU Shunhua, WANG Changdan, et al. Calculation method of ballasted track settlement caused by subgrade differential settlement[J]. Journal of Tongji University (Natural Science), 2013, 41(8): 1237-1242.
[22]	张乾,蔡小培,钟阳龙,等. 路基沉降区双块式无砟轨道损伤行为与轨面不平顺分析[J]. 中南大学学报(自然科学版),2022,53(7): 2573-2581. ZHANG Qian, CAI Xiaopei, ZHONG Yanglong, et al. Analysis of damage and rail irregularity of double-block ballastless track in subgrade settlement area[J]. Journal of Central South University (Science and Technology), 2022, 53(7): 2573-2581.
[23]	郭宇,高建敏,孙宇,等. 路基沉降与双块式无砟轨道轨面几何变形的映射关系[J]. 铁道学报,2016,38(9): 92-100. GUO Yu, GAO Jianmin, SUN Yu, et al. Mapping relationship between subgrade settlement and rail deflection of the double-block ballastless track[J]. Journal of the China Railway Society, 2016, 38(9): 92-100.
[24]	JIANG H G, LI X L, XIN G F, et al. Geometry mapping and additional stresses of ballastless track structure caused by subgrade differential settlement under self-weight loads in high-speed railways[J]. Transportation Geotechnics, 2019, 18: 103-110. doi: 10.1016/j.trgeo.2018.10.007

施引文献

期刊类型引用(1)
1. 武小年，吴庭，黄昭文，张润莲. 基于混沌映射和NFSR的16比特动态S盒构造方法. 信息安全与通信保密. 2024(11): 10-19 . 百度学术
其他类型引用(0)

附加材料(0)

访问统计

点击查看大图

图(5) / 表(1)

计量

文章访问数: 278
HTML全文浏览量: 173
PDF下载量: 90
被引次数: 1

1. 双制式列车接地系统
1.1 车—轨链式电路
2. 双制式牵引供电系统建模
2.1 交、直流区段牵引供电系统建模
2.2 交、直流区段联合牵引供电计算
3. 车-地一体的交、直流区段联合牵引供电计算
4. 列车接地保护电阻优化模型
4.1 目标函数
4.2 约束条件
5. 算例分析
5.1 仿真条件
5.2 实测验证
5.3 对列车接地系统配置的讨论与分析
6. 结　论

1. 双制式列车接地系统
1.1 车—轨链式电路
2. 双制式牵引供电系统建模
2.1 交、直流区段牵引供电系统建模
2.2 交、直流区段联合牵引供电计算
3. 车-地一体的交、直流区段联合牵引供电计算
4. 列车接地保护电阻优化模型
4.1 目标函数
4.2 约束条件
5. 算例分析
5.1 仿真条件
5.2 实测验证
5.3 对列车接地系统配置的讨论与分析
6. 结　论

参考文献(24)

施引文献

附加材料(0)

访问统计

结构	弹性模量/Pa	泊松比	说明
钢轨	2.10 × 10¹¹	0.30	T60 轨
道床板	3.25 × 10¹⁰	0.17	C40 混凝土
支承层	2.55 × 10¹⁰	0.17	C20 混凝土
路基	180	0.25	密度 300 kg/m³

基于双层梁理论的隧道下穿诱发无砟铁路变形计算方法

doi: 10.3969/j.issn.0258-2724.20230033

作者简介: 雷鸣（1977—），男，副教授，博士，研究方向为岩土工程，E-mail：540584998@qq.com

通讯作者: 张丙强（1979—），男，教授，博士，研究方向为地下工程，E−mail：zbq@fjut.edu.cn

计量

出版历程

Theoretical Method for Calculating Rail Deformation of Ballastless Railway Caused by Tunnel Undercrossing Based on Dual Beam Model

1. 双制式列车接地系统

1.1 车—轨链式电路

2. 双制式牵引供电系统建模

2.1 交、直流区段牵引供电系统建模

2.2 交、直流区段联合牵引供电计算

3. 车-地一体的交、直流区段联合牵引供电计算

4. 列车接地保护电阻优化模型

4.1 目标函数

4.2 约束条件

5. 算例分析

5.1 仿真条件

5.2 实测验证

5.2.1 直流区段

5.2.2 交流区段

5.3 对列车接地系统配置的讨论与分析

6. 结 论

期刊类型引用(1)

其他类型引用(0)

计量

出版历程

目录

1. 双制式列车接地系统

1.1 车—轨链式电路

2. 双制式牵引供电系统建模

2.1 交、直流区段牵引供电系统建模

2.2 交、直流区段联合牵引供电计算

3. 车-地一体的交、直流区段联合牵引供电计算

4. 列车接地保护电阻优化模型

4.1 目标函数

4.2 约束条件

5. 算例分析

5.1 仿真条件

5.2 实测验证

5.3 对列车接地系统配置的讨论与分析

6. 结 论

作者简介:
雷鸣（1977—），男，副教授，博士，研究方向为岩土工程，E-mail：540584998@qq.com

通讯作者:
张丙强（1979—），男，教授，博士，研究方向为地下工程，E−mail：zbq@fjut.edu.cn

6. 结　论

6. 结　论