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

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

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
计量
- 文章访问数: 280
- HTML全文浏览量: 176
- PDF下载量: 90
- 被引次数: 0
出版历程
- 收稿日期: 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全文

在城市地铁和城际铁路修建中，经常会遇到新建隧道近距离下穿既有铁路的情况^[1-2]. 新建隧道开挖施工不可避免会引起铁路路基沉降^[3-4]，诱发轨道产生几何变形，加剧轮轨动力相互作用，进而影响轨道结构服役性能和列车乘坐舒适性，严重时将造成列车脱轨和铁路停运^[5-6]. 因此，研究隧道下穿施工诱发上覆既有铁路轨面变形具有重要的工程意义.

目前，国内外学者采用理论方法分析新建隧道施工对邻近既有结构的影响时，通常采用两阶段分析方法：第1阶段，忽略既有结构的影响，计算新建隧道施工引起邻近结构所处位置土体的附加应力或自由变形；第2阶段，将既有结构视为连续地基上梁模型，并将土体附加应力或自由变形视为外部作用，应用数值方法^[7-8]或简化方法等^[9-10]分析邻近既有结构的变形和内力的变化. 两阶段分析方法由于使用简单，被广泛用于分析隧道施工对邻近管线^[11-16]、上方铁路^[17]、和邻近建筑物^[9,18]的影响.

同时，针对路基沉降引起铁路线路轨道变形开展了一定的研究：邹春华等^[19]、付龙龙等^[20]分别采用室内试验和数值方法对余弦型路基沉降与有砟轨道变形之间的互相关系进行了研究，并且，邹春华等^[21]基于弹性点支承梁模型提出了有砟铁路路基沉降引起轨面变形的计算模型；考虑无砟道床为板状多层结构，张乾等^[22]采用梁-板模型分别模拟钢轨及无砟道床，建立了无砟轨道-路基-土体的三维空间有限元模型，对地面沉降对无砟轨道系统平顺性的影响进行了分析，并指出支承层和路基表层间沉降差较大易出现离缝问题；郭宇等^[23]、Jiang等^[24]分别对双块式和整体板式无砟轨道路基沉降引起轨面变形几何特征与路基沉降波长、幅值之间的对应关系进行了研究.

以上相关理论研究中邻近既有结构均被简化为弹性地基上的单层梁模型，不能模拟地下结构与底部地基局部分离现象，不适用于计算新建隧道下穿施工引起既有无砟铁路的轨面变形. 为此，本文将无砟铁路结构视为双层地基梁模型，采用改进的Winkler弹性地基模型模拟轨道板与地基的互相作用，建立考虑轨道板底部局部脱空影响的路基沉降引起无砟铁路轨道变形的理论计算模型，进而采用解析方法推导了隧道下穿施工引起轨道变形计算式，为预测隧道下穿既有无砟铁路引起轨道变形提供理论支持.

1. 问题描述及控制方程

新建隧道下穿引起既有无砟铁路轨道变形的计算模型如所示，以隧道与铁路相交处为坐标原点，铁路线路方向和x轴建立水平坐标系. 图中： $\beta$ 为隧道与铁路的水平交角， ${R_{\text{t}}}$ 为隧道开挖半径， ${{\textit{z}}_{\text{t}}}$ 为隧道中心到铁路路基顶面的竖向距离， ${l_{\text{t}}}$ 为轨道叠合梁底部局部脱空区长度， $s\left( x \right)$ 为隧道斜交下穿引起路基沉降，（0，x_i）为铁路上任意点坐标.

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

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

下载: 全尺寸图片幻灯片

本文拟采用两阶段分析方法进行求解，并作如下基本假定：

1）钢轨视为两端自由的上层Euler梁，其竖向变形、弹性模量和截面惯性矩分别为 ${w_{\text{r}}}$ 、 ${E_{\text{r}}}$ 和 ${I_{\text{r}}}$ ；

2）轨道板和其下部支承层简化为下层弹性叠合梁，其竖向变形、弹性模量和截面惯性矩分别为 ${w_{\text{b}}}$ 、 ${E_{\text{b}}}$ 和 ${I_{\text{b}}}$ ；

3）采用改进Winkler地基模拟叠合梁下的路基作用，地基的弹性系数为 ${k_{\text{b}}}$ ；

4）采用线性弹簧模拟扣件系统，钢轨扣件系统单位长度的等效弹性系数为 ${k_{\text{r}}}$ ；

5）隧道施工引起周围土体沉降曲线符合高斯分布.

根据Euler梁挠曲理论，叠合梁变形诱发钢轨（上层梁）变形的控制方程为

${E_{\text{r}}}{I_{\text{r}}}\frac{{{{\rm{d}}^4}{w_{\text{r}}}\left( x \right)}}{{{\rm{d}}{x^4}}} = - {k_{\text{r}}}\left[ {{w_{\text{r}}}\left( x \right) - {w_{\text{b}}}\left( x \right)} \right].$

(1)

新建隧道施工前，叠合梁在轨道系统及自重作用下产生均匀沉降. 新建隧道下穿施工引起路基沉降，改变轨道板叠合梁与地基土间的互相作用，导致叠合梁产生变形. 若叠合梁变形大于路基顶面土体沉降，则叠合梁与地基土体间保持接触，两者之间仍然存在互相作用；否则，叠合梁与地基土体将发生分离（即局部脱空），两者互相作用力为0. 轨道板叠合梁（下层梁）变形的控制方程为

$\qquad\qquad\qquad {E_{\text{b}}}{I_{\text{b}}}\dfrac{{{{\rm{d}}^4}{w_{\text{b}}}\left( x \right)}}{{{\rm{d}}{x^4}}} =\left\{ {\begin{array}{*{20}{l}} { {k_{\text{r}}}\left[ {{w_{\text{r}}}\left( x \right) - {w_{\text{b}}}\left( x \right)} \right] + \gamma h},\quad{\left| x \right| \leqslant {l_{\text{t}}}/2}, \\ {{k_{\text{r}}}\left[ {{w_{\text{r}}}\left( x \right) - {w_{\text{b}}}\left( x \right)} \right] + \gamma h - {k_{\text{b}}}\left[ {{w_{\text{b}}}\left( x \right) - s\left( x \right)} \right]},\quad{\left| x \right| > {l_{\text{t}}}/2} , \end{array}} \right.$

(2)

式中： $\gamma$ 、 $h$ 分别为叠合梁重度和高度； $s\left( x \right)$ 计算式^[7,13-14]如式（3）.

$s\left( x \right) = \frac{{{\text{π}} R_{\text{t}}^2{V_{\text{t}}}}}{{\sqrt {2{\text{π}} } {i_{\text{s}}}}}\exp \left[ { - 0.5{{\left( {\frac{{x\sin \;\beta }}{{{i_{\text{s}}}}}} \right)}^2}} \right],$

(3)

式中： ${V_{\rm{t}}}$ 为隧道施工引起的平均地层损失比，与施工条件相关； ${i_{\text{s}}}$ 为路基顶面沉降槽宽度半径，采用建议公式 ${i_{\text{s}}} = \kappa {{\textit{z}}_{\rm{t}}}$ 进行计算， $\kappa$ 为地层沉降槽宽度参数，与土体条件相关.

2. 解析方法

将式（1）变形可得

${w_{\text{b}}}\left( x \right) = \frac{{{E_{\text{r}}}{I_{\text{r}}}}}{{{k_{\text{r}}}}}\frac{{{{\rm{d}}^4}{w_{\text{r}}}\left( x \right)}}{{{\rm{d}}{x^4}}} + {w_{\text{r}}}\left( x \right).$

(4)

对式（4）求4次导数，可得

$\frac{{{{\rm{d}}^4}{w_{\text{b}}}\left( x \right)}}{{{\rm{d}}{x^4}}} = \frac{{{E_{\text{r}}}{I_{\text{r}}}}}{{{k_{\text{r}}}}}\frac{{{{\rm{d}}^8}{w_{\text{r}}}\left( x \right)}}{{{\rm{d}}{x^8}}} + \frac{{{{\rm{d}}^4}{w_{\text{r}}}\left( x \right)}}{{{\rm{d}}{x^4}}}.$

(5)

将式（4）、（5）依次代入式（2）中各分段方程，整理后可得关于钢轨变形的控制方程.

1）当 $\left| x \right| \leqslant {l_{\text{t}}}/2$ 时

脱空区上方钢轨的变形控制方程为

$\frac{{{{\rm{d}}^8}{w_{\text{r}}}\left( x \right)}}{{{\rm{d}}{x^8}}} + C\frac{{{{\rm{d}}^4}{w_{\text{r}}}\left( x \right)}}{{{\rm{d}}{x^4}}} = \frac{{{k_{\text{r}}} \gamma h}}{{{E_{\text{r}}}{I_{\text{r}}} {E_{\text{b}}}{I_{\text{b}}}}},$

(6)

式中： $C = \dfrac{{{k_{\text{r}}}}}{{{E_{\text{b}}}{I_{\text{b}}}}} + \dfrac{{{k_{\text{r}}}}}{{{E_{\text{r}}}{I_{\text{r}}}}}$ .

求解式（6），可得钢轨中间脱空段的变形表达式为

$\begin{split} \\[-10pt] & {w_{\text{r}}}\left( x \right) = {{\rm{e}}^{{\beta _1}x}}\left( {{c_1}\cos\; {\beta _1}x + {c_2}\sin \;{\beta _1}x} \right) + \\ &\quad {{\rm{e}}^{ - {\beta _1}x}}\left( {{c_3}\cos \;{\beta _1}x + {c_4}\sin\; {\beta _1}x} \right) + {c_5} + {c_6}x +\\ &\quad {c_7}{x^2} + {c_8}{x^3} + \frac{{\gamma h}}{{24\left( {{E_{\text{r}}}{I_{\text{r}}} + {E_{\text{b}}}{I_{\text{b}}}} \right)}}{x^4}, \end{split}$

(7)

式中： ${\beta _1} = \sqrt[\uproot{13}{4}]{{\dfrac{{{k_{\text{r}}}}}{{4{E_{\text{b}}}{I_{\text{b}}}}} + \dfrac{{{k_{\text{r}}}}}{{4{E_{\text{r}}}{I_{\text{r}}}}}}}$ ；c₁～c₈为待定参数，可用边界条件求出.

将式（7）代入式（4），可求得中间脱空区上方的轨道板叠合梁变形表达式为

$\begin{split} & {w_{\text{b}}}\left( x \right) = {D_1}\left[ {{{\rm{e}}^{{\beta _1}x}}\left( {{c_1}\cos \;{\beta _1}x + {c_2}\sin \;{\beta _1}x} \right)}+ \right. \\ &\quad \left. { {{\rm{e}}^{ - {\beta _1}x}}\left( {{c_3}\cos \;{\beta _1}x + {c_4}\sin \;{\beta _1}x} \right)} \right] + \\ &\quad {c_5} + {c_6}x + {c_7}{x^2} + {c_8}{x^3} +\\ &\quad \frac{{\gamma h}}{{24\left( {{E_{\text{r}}}{I_{\text{r}}} + {E_{\text{b}}}{I_{\text{b}}}} \right)}}{x^4} + \frac{{\gamma h {E_{\text{r}}}{I_{\text{r}}}}}{{{k_{\text{r}}}\left( {{E_{\text{r}}}{I_{\text{r}}} + {E_{\text{b}}}{I_{\text{b}}}} \right)}}, \end{split}$

(8)

式中： ${D_1} = 1 - 4{\beta _1^4}\dfrac{{{E_{\text{r}}}{I_{\text{r}}}}}{{{k_{\text{r}}}}}$ .

2）当 $\left| x \right| > {l_{\text{t}}}/2$ 时

脱空区外侧钢轨变形控制方程为

$\begin{split} & \frac{{{{\rm{d}}^8}{w_{\text{r}}}\left( x \right)}}{{{\rm{d}}{x^8}}} + E\frac{{{{\rm{d}}^4}{w_{\text{r}}}\left( x \right)}}{{{\rm{d}}{x^4}}} + F {w_{\text{r}}}\left( x \right) = \\ &\quad F s\left( x \right) + \frac{{{k_{\text{r}}} \gamma h}}{{{E_{\text{r}}}{I_{\text{r}}} {E_{\text{b}}}{I_{\text{b}}}}}, \end{split}$

(9)

式中： $E = \dfrac{{{k_{\text{r}}}\left( {{E_{\text{r}}}{I_{\text{r}}} + {E_{\text{b}}}{I_{\text{b}}}} \right)}}{{{E_{\text{r}}}{I_{\text{r}}} {E_{\text{b}}}{I_{\text{b}}}}} + \dfrac{{{k_{\text{b}}}}}{{{E_{\text{b}}}{I_{\text{b}}}}}$ ， $F = \dfrac{{{k_{\text{r}}}{k_{\text{b}}}}}{{{E_{\text{r}}}{I_{\text{r}}} {E_{\text{b}}}{I_{\text{b}}}}}$ .

求解式（9），并结合边界条件 ${\left. {{w_{\text{r}}}} \right|_{x \to \infty }} = 0$ ，可得脱空区外侧钢轨变形表达式为

$\begin{split} & {w_{\text{r}}}\left( x \right) = {\bar w_{\text{r}}}\left( x \right) + \frac{{\gamma h}}{{{k_{\text{b}}}}} + {{\rm{e}}^{ - {\beta _2}x}}\left( {{c_9}\cos\; {\beta _2}x + {c_{10}}\sin \;{\beta _2}x} \right) + \\ &\quad {{\rm{e}}^{ - {\beta _3}x}}\left( {{c_{1{\text{1}}}}\cos \;{\beta _3}x + {c_{1{\text{2}}}}\sin \;{\beta _3}x} \right), \end{split}$

(10)

式中：c₉～c₁₂为待定系数，可由边界条件确定. ${\bar w_{\text{r}}}\left( x \right)$ 为路基沉降部分对应的方程特解，可采用展开级数法进行求解， $\beta _2$ 、 $\beta _3$ 取值如式（11）.

$\qquad\qquad\qquad \left\{\begin{array}{l} \beta _2=\sqrt[\uproot{13}{{4}}]{{\dfrac{{E - \sqrt {{E^2} - 4F} }}{8}}},\;{\beta _3} = \sqrt[\uproot{13}{{4}}]{{\dfrac{{E + \sqrt {{E^2} - 4F} }}{8}}},\quad\;{E^2} - 4F > 0,\\ {\beta _2} = \dfrac{{\sqrt {2\sqrt[4]{F} - \sqrt {2\sqrt F } - E} }}{2},\;{\beta _3} = \dfrac{{\sqrt {2\sqrt[4]{F} + \sqrt {2\sqrt F } - E} }}{2},\quad\;{E^2} - 4F < 0. \end{array}\right.$

(11)

将 $s\left( x \right)$ 采用三角函数进行展开，记为

$s\left( x \right) = {a_0} + \sum\limits_{n = 1}^{\infty } {\left( {{a_n}\cos {\frac{{n{\text{π}} x}}{{{l_{\text{s}}}}}} } \right)},$

(12)

式中： ${a_0} = \dfrac{1}{{{l_{\text{s}}}}}\displaystyle\int\limits_0^{ + {{{l_{\text{s}}}} / 2}} {s\left( x \right){\rm{d}}x}$ ， ${a_n} = \dfrac{2}{{{l_{\text{s}}}}}\displaystyle\int\limits_0^{ + {{{l_{\text{s}}}} / 2}} {s\left( x \right)\cos {\frac{{n{\text{π}} x}}{{{l_{\text{s}}}}}}{\rm{d}}x}$ ， ${l_{\text{s}}}$ 为线路计算长度，n为三角级数项的编号.

路基沉降对应的轨道变形特解为

${\bar w_{\text{r}}}\left( x \right) = {a_0} + \sum\limits_{n = 1}^{\infty } {{\dfrac{{F\cos \left( {\dfrac{{n{\text{π}} x}}{{{l_{\text{s}}}}}} \right)}}{{{{\left( {\dfrac{{n{\text{π}} }}{{{l_{\text{s}}}}}} \right)}^8} + E{{\left( {\dfrac{{n{\text{π}} }}{{{l_{\text{s}}}}}} \right)}^4} + F}}} }.$

(13)

同理，将式（10）代入式（4），可求得脱空区外侧轨道板叠合梁变形表达式为

$\begin{split} & {w_{\text{b}}}\left( x \right) = {\bar w_{\text{r}}}\left( x \right) + \frac{{{E_{\text{r}}}{I_{\text{r}}}}}{{{k_{\text{r}}}}}\frac{{{{\rm{d}}^4}{{\bar w}_{\text{r}}}\left( x \right)}}{{{\rm{d}}{x^4}}} + \frac{{\gamma h}}{{{k_{\text{b}}}}} +\\ &\quad {D_2}{{\rm{e}}^{ - {\beta _2}x}}\left( {{c_9}\cos \;{\beta _2}x + {c_{10}}\sin \;{\beta _2}x} \right) +\\ &\quad {D_3}{{\rm{e}}^{ - {\beta _3}x}}\left( {{c_{11}}\cos \;{\beta _3}x + {c_{12}}\sin \;{\beta _3}x} \right), \end{split}$

(14)

式中： ${D_{\text{2}}} = 1 - 4\dfrac{{{E_{\text{r}}}{I_{\text{r}}}}}{{{k_{\text{r}}}}}{\beta _{\text{2}}^4}$ ， ${D_{\text{3}}} = 1 - 4\dfrac{{{E_{\text{r}}}{I_{\text{r}}}}}{{{k_{\text{r}}}}}{\beta _{\text{3}}^4}$ .

根据Euler梁变形理论，可求得钢轨和轨道板叠合梁的转角 ${\theta _j}$ 、弯矩 ${M_j}$ 和剪力 ${Q_j}$ （ $j$ =r，b分别表示钢轨和叠合梁）：

$\left\{ {\begin{array}{*{20}{l}} {{\theta _j} = \dfrac{{{\rm{d}}{w_j}}}{{{\rm{d}}x}}}, \\ {{M_j} = {{ - }}{E_j}{I_j}\dfrac{{{{\rm{d}}^2}{w_j}}}{{{\rm{d}}{x^2}}}}, \\ {{Q_j} = {{ - }}{E_j}{I_j}\dfrac{{{{\rm{d}}^3}{w_j}}}{{{\rm{d}}{x^3}}}}. \end{array}} \right.$

(15)

1）由于模型对称性，当 $x = 0$ 时，钢轨和叠合梁的转角与剪力都为 ${\text{0}}$ ，即

$\left\{ {\begin{array}{*{20}{c}} {{{\left. {\dfrac{{{\rm{d}}{w_j}}}{{{\rm{d}}x}}} \right|}_{x = 0}} = 0}, \\ {{{\left. {\dfrac{{{{\rm{d}}^3}{w_j}}}{{{\rm{d}}{x^3}}}} \right|}_{x = 0}} = 0} . \end{array}} \right.$

(16)

2）当 $x = {l_{\rm{t}}}/2$ 时，钢轨和叠合梁的位移、转角、弯矩与剪力都保持连续，即

$\left\{ {\begin{array}{*{20}{l}} {{{\left. {{w_j}} \right|}_{{x^ - } = {{{l_{\text{t}}}} / 2}}} = {{\left. {{w_j}} \right|}_{{x^ + } = {{{l_{\text{t}}}} / 2}}}}, \\ {{{\left. {\dfrac{{{\rm{d}}{w_j}}}{{{\rm{d}}x}}} \right|}_{{x^ - } = {{{l_{\text{t}}}} / 2}}} = {{\left. {\dfrac{{{\rm{d}}{w_j}}}{{{\rm{d}}x}}} \right|}_{{x^ + } = {{{l_{\text{t}}}} / 2}}}} ,\\ {{{\left. {\dfrac{{{{\rm{d}}^2}{w_j}}}{{{\rm{d}}{x^2}}}} \right|}_{{x^ - } = {{{l_{\text{t}}}} / 2}}} = {{\left. {\dfrac{{{{\rm{d}}^2}{w_j}}}{{{\rm{d}}{x^2}}}} \right|}_{{x^ + } = {{{l_{\text{t}}}} / 2}}}}, \\ {{{\left. {\dfrac{{{{\rm{d}}^3}{w_j}}}{{{\rm{d}}{x^3}}}} \right|}_{{x^ - } = {{{l_{\text{t}}}} / 2}}} = {{\left. {\dfrac{{{{\rm{d}}^3}{w_j}}}{{{\rm{d}}{x^3}}}} \right|}_{{x^ + } = {{{l_{\text{t}}}} / 2}}}} . \end{array}} \right.$

(17)

根据上述推导的钢轨和轨道板叠合梁变形计算式，利用边界条件（式（16）、（17））可建立起一个关于待定系数的12 × 12阶线性方程组，编制MATLAB程序可求出待定系数. 回代入式（7）、（10）即为隧道下穿施工诱发无砟铁路钢轨变形. 由于脱空区长度 ${l_{\rm{t}}}$ 未知，需要采用多次迭代法进行计算.

3. 数值验证

郭宇等^[23]基于通用有限元软件ABAQUS，对余弦型路基沉降引起双块式无砟轨道的几何变形规律进行了模拟分析. 模型中采用空间梁单元模拟钢轨，采用弹簧－阻尼器模拟扣件系统，道床板、支承层和路基结构均采用实体单元模拟，结构材料参数见表1. 在轨枕与混凝土道床板、道床板与支承层、支承层与路基表面均进行接触设置，并充分考虑轨道结构和路基之间可能出现离缝甚至空吊现象.

表 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

| 显示表格

为验证本文理论方法的正确性，采用本文理论方法对文献[23]中不同路基沉降幅值和波长条件下的轨面变形进行分析，模型中钢轨抗弯刚度为6.624 MN·m²，叠合梁抗弯刚度为1.225 GN·m²，扣件刚度为30 MN·m⁻¹，路基弹性系数为200 MN·m⁻¹. 本文计算结果与文献[23]的分析数据对比如图2所示，图2（a）中路基沉降波长取20 m，幅值分别取5 mm和20 mm；图2（b）中路基沉降幅值10 mm，波长分别取10 m和40 m.

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

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

下载: 全尺寸图片幻灯片

从图2中可看出：当轨道结构未产生脱空（5 mm/20 m、10 mm/40 m）时，理论解析计算结果与数值模拟结果几乎一致；当轨道结构产生脱空（20 mm/20 m、10 mm/10 m）时，考虑轨道结构脱空的理论解析计算结果与数值模拟结果基本一致，小于未考虑轨道结构脱空的理论计算结果；当路基沉降为10 mm/10 m时，未考虑轨道结构脱空的理论方法求解出的钢轨变形最大值（7.9 mm）约为考虑轨道结构脱空的理论计算结果（2.3 mm）的3.4倍. 分析表明，当路基不均匀沉降引起无砟轨道结构产生离缝或空吊时，采用本文提出的理论方法仍能较为准确计算出轨面各点处变形值.

4. 参数分析

为控制新建隧道下穿既有铁路引起轨道变形，实际工程中通常采用减小隧道施工引起周围地层损失率或选择不同穿越角度和隧道埋深等措施. 为此，本节将针对上述3个参数对隧道下穿施工诱发轨面变形的影响进行分析. 算例中：开挖半径 ${R_{\text{t}}}$ =3.1 m，地层损失率 ${V_{\text{t}}}$ =0.5%，隧道埋深 ${{\textit{z}}_{\text{t}}}$ =6 m，地层沉降槽宽度参数 $\kappa$ =0.5，隧道与铁路线路垂直相交；无砟铁路结构参数取值与验证算例相同.

4.1 隧道埋深的影响

为分析隧道埋深对下穿施工引起钢轨变形的影响，采用本文理论方法分别计算不同隧道埋深（ ${{\textit{z}}}_{\text{t}}$ =2，4，6，10，20，40 m）条件下轨道的变形，其中轨道中点变形值及轨道板叠合梁脱空区宽度随隧道埋深的变化曲线如图3所示.

图 3 隧道埋深的影响

Figure 3. Influence of buried depth of new tunnel

下载: 全尺寸图片幻灯片

从图3中可看出，随着隧道埋深的逐步增大，轨道中点变形值及叠合梁底部脱空区宽度均先增大而后减小，当隧道埋深为6 m时达到峰值. 这是由于随着隧道埋深的增加，隧道施工引起路基顶部的沉降槽宽度将逐渐增大，但沉降槽最大沉降值却逐渐减小. 当新建隧道埋深较小时，隧道施工引起路基顶部的沉降槽的宽度较小，轨道叠合梁底部脱空区宽度和最大变形值较小；随着隧道埋深的增大，路基顶部沉降槽的宽度将增大，导致叠合梁底部脱空区宽度和钢轨中心点变形逐渐增大；当隧道埋深增大到一定值时，路基顶部的沉降逐渐平缓，导致叠合梁脱空区宽度逐渐减小，直至随路基协同变形. 由此可见，选择适当的隧道埋深穿越既有铁路可以有效控制隧道下穿施工引起上方轨道变形.

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

当隧道施工中地层损失率 ${V_{\text{t}}}$ 分别取0.25%、0.50%、1.00%、1.50%、2.00%和2.50%时，隧道下穿施工诱发钢轨中点变形值及叠合梁底部脱空区宽度变化如图4所示. 由图可知：随着隧道施工过程中周围地层损失率的逐渐增大，钢轨中点变形值及叠合梁底部脱空区宽度均逐渐增大；当隧道施工中地层损失率从0.25%增大到2.50%时，钢轨中点变形值和叠合梁脱空区宽度分别从9.7 mm、4 m增大到48.6 mm和13 m，依次增长了约4.00倍、2.25倍. 由此可见，当减小隧道施工中地层损失率是控制隧道下穿施工引起钢轨变形的有效措施之一.

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

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

下载: 全尺寸图片幻灯片

4.3 隧道与铁路水平夹角的影响

当新建隧道与铁路间的水平夹角 $\;\beta$ 分别取15°、30°、45°、60°、75° 和90° 时，新建隧道下穿施工诱发轨道中点变形值及叠合梁底部脱空区宽度变化如图5所示. 由图可知：随着隧道与铁路线路间水平夹角的逐渐增大，钢轨中点变形值逐渐减小，而叠合梁底部脱空区宽度逐渐先增大、后基本保持不变. 这是由于当新建隧道与既有铁路线路间的水平夹角较小时，隧道施工引起上方铁路路基顶部沉降较为平缓，轨道叠合梁与路基之间保持跟随变形；随着隧道与铁路线路间水平夹角的增大，路基顶部沉降槽宽度逐渐变窄，轨道叠合梁与路基间产生变形差，轨道叠合梁底部将产生局部脱空，轨道中点变形将减小. 由此可见，选择较大角度穿越既有铁路有利于控制隧道下穿施工引起轨道的变形.

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

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

下载: 全尺寸图片幻灯片

5. 小　结

1）将无砟铁路轨道结构简化为双层地基梁，采用无拉力弹性地基模拟轨道板叠合梁与路基的互相作用，并根据叠合梁与路基间的变形差判别两者的接触状态，建立隧道下穿施工引起无砟铁路轨道变形的控制方程.

2）将双层梁分为中间脱空段和两端接地段共三部分，采用解析方法分别推导出上、下层梁（钢轨和轨道板叠合梁）变形的通解和特解表达式，并结合边界条件和分段处连续条件求出变形表达式中的待定系数.

3）随着隧道埋深的增大，轨道最大变形值及叠合梁底部脱空区宽度均呈先增大而后减小趋势. 随着隧道施工中周围地层损失率的增大，钢轨最大变形值及叠合梁底部脱空区宽度均逐渐增大. 随着隧道与铁路线路间水平夹角的增大，钢轨最大变形值逐渐减小，而叠合梁底部脱空区宽度逐渐增大、后保持不变.

图 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

施引文献

附加材料(0)

访问统计

点击查看大图

图(5) / 表(1)

计量

文章访问数: 280
HTML全文浏览量: 176
PDF下载量: 90
被引次数: 0

1. 问题描述及控制方程
2. 解析方法
3. 数值验证
4. 参数分析
4.1 隧道埋深的影响
4.2 隧道施工引起地层损失率的影响
4.3 隧道与铁路水平夹角的影响
5. 小　结

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

doi: 10.3969/j.issn.0258-2724.20230033

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

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

计量

出版历程