列车一维碰撞能量管理的综合评价模型

张敬科; 朱涛; 王小瑞; 肖守讷; 阳光武; 杨冰

doi:10.3969/j.issn.0258-2724.20200034

列车一维碰撞能量管理的综合评价模型

doi: 10.3969/j.issn.0258-2724.20200034

西南交通大学牵引动力国家重点实验室，四川成都 610031

基金项目: 国家自然科学基金(52172409)；四川省杰出青年基金(22JCQN0156)

详细信息

作者简介:
张敬科（1992—），男，博士研究生，研究方向为机车车辆碰撞动力学，E-mail：jkzhang@my.swjtu.edu.cn

中图分类号: U270.1
计量
- 文章访问数: 590
- HTML全文浏览量: 255
- PDF下载量: 9
- 被引次数: 10
出版历程
- 收稿日期: 2020-02-18
- 修回日期: 2020-04-26
- 网络出版日期: 2020-05-14
- 刊出日期: 2020-05-14

Comprehensive Evaluation Model for One-Dimensional Crash Energy Management of Trains

State Key Laboratory of Traction Power，Southwest Jiaotong University，Chengdu 610031，China

摘要

摘要:
为了更科学全面地对列车碰撞能量方案进行配置和评价，以列车最大平均加速度、最大瞬时加速度和能量利用率为评价指标，利用序关系分析法和熵值法得到的主、客观指标权重推导出更加合理的综合权重系数；采用非线性加权综合法融合权重信息和评价指标信息，获得了列车一维碰撞能量管理综合评价模型；采用此模型对某型列车的15种能量配置方案进行综合评价. 研究结果表明：此模型综合评价结果与各指标性能变化幅度分析结果一致，具有较好的适用性和可靠性；方案4为最优方案，其最大瞬时加速度指标性能比方案8下降14.03%，但最大平均加速度指标性能和能量利用率指标性能分别比方案8提升了4.92%和19.44%；该结果既体现了对安全性影响最大的最大平均加速度指标的重要性，又综合考虑了与经济性相关的能量利用率指标.
- 碰撞能量管理 /
- 序关系分析法 /
- 熵值法 /
- 综合评价模型
Abstract:
In order to configure and evaluate the train collision energy scheme more scientifically and comprehensively, the maximum average acceleration, maximum instantaneous acceleration and energy utilization rate of trains were taken as evaluation indexes, and a more reasonable comprehensive weight coefficient was derived by using the subjective and objective index weights obtained by the order relation analysis and entropy method. A comprehensive evaluation model of train one-dimensional collision energy management was obtained by fusing weight information and evaluation index information with nonlinear weighted synthesis method, then 15 energy allocation schemes of a certain type of train were comprehensively evaluated by using this model. The results indicate that the comprehensive evaluation results are consistent with that of the change range of all indexes, and it has better applicability and reliability. Scheme 4 is the optimal scheme, in which the maximum instantaneous acceleration index decreases by 14.03%, but the maximum average acceleration index and energy utilization rate index increase by 4.92% and 19.44% compared with scheme 8 respectively. These results both reflect the importance of the maximum average acceleration index, which has the greatest impact on safety, and comprehensively considers the energy utilization rate index related to economy.
- crash energy management /
- order relationship analysis method /
- entropy method /
- comprehensive evaluation model

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. Flow chart for comprehensive evaluation of train one-dimensional collision energy management

下载: 全尺寸图片幻灯片

图 2 15种能量配置方案的各指标计算结果

Figure 2. Calculation results of each index for 15 energy allocation schemes

下载: 全尺寸图片幻灯片

图 3 各能量配置方案综合评价结果排序

Figure 3. Ranking of comprehensive evaluation results of each energy allocation scheme

下载: 全尺寸图片幻灯片

图 4 方案4、8、10、11、13各指标数值变化幅度百分比

Figure 4. Percentage changes of each indicator value of schemes 4，8，10，11 and 13

下载: 全尺寸图片幻灯片

${r_k}$ 赋值

Assignment reference table of ${r_k}$

${r_k}$	说明
1.0	指标 ${x_{k - 1}}$ 与指标 ${x_k}$ 具有同样重要性
1.2	指标 ${x_{k - 1}}$ 与指标 ${x_k}$ 稍微重要
1.4	指标 ${x_{k - 1}}$ 与指标 ${x_k}$ 明显重要
1.6	指标 ${x_{k - 1}}$ 与指标 ${x_k}$ 强烈重要
1.8	指标 ${x_{k - 1}}$ 与指标 ${x_k}$ 极端重要

下载: 导出CSV

表 2 列车碰撞能量配置方案^[10]

Table 2. Train collision energy allocation schemes

方案	头车钩缓装置		中间车钩缓装置		吸能装置
方案	阻抗力/kN	额定容量/kJ	阻抗力/kN	额定容量/kJ	阻抗力/kN	额定容量/kJ
1	1000	387.5	1000	743	1400	970
2		387.5	1200	883
3		387.5	1400	1023
4	1200	457.5	1000	743	1400	970
5			1200	883
6			1400	1023
7	1000	387.5	1000	743	1600	1110
8			1200	883
9			1400	1023
10	1200	457.5	1000	743	1600	1110
11			1200	883
12			1400	1023
13	1400	527.5	1000	743	1600	1110
14			1200	883
15			1400	1023

下载: 导出CSV

表 3 15种能量配置方案的各指标数据一致化和无量纲化处理结果

Table 3. Data consistency and dimensionless processing results of 15 energy allocation schemes

方案	能量利用率/%	最大平均加速度/（m•s⁻²）	最大瞬时加速度/（m•s⁻²）
1	68.5178	18.7772	36.2090
2	54.1003	37.5183	42.3262
3	42.0498	45.7863	33.6498
4	64.1945	59.7502	56.4960
5	51.3224	50.1960	53.4373
6	40.9935	48.3586	38.5186
7	46.0210	54.0545	55.4972
8	51.7137	56.8104	64.4235
9	35.2421	47.0725	41.2653
10	62.3947	55.5243	54.8106
11	49.7770	56.4430	64.2363
12	38.2938	52.0334	43.3251
13	60.1842	55.7080	54.4985
14	47.7034	54.9731	63.6745
15	37.4918	56.9942	47.6321

下载: 导出CSV

指标序关系及 ${r_k}$ 赋值表

Index order relation and assignment table of ${r_k}$

类别	取值
指标序关系	${x_2} \succ {x_3} \succ {x_1} (x_1^* \succ x_2^* \succ x_3^*)$
${r_2}$	1.8
${r_3}$	1.2

下载: 导出CSV

表 5 G1-法计算的各评价指标权重

Table 5. Weight of each evaluation index by G1 method

指标类型	权重
最大平均加速度	0.4954
最大瞬时加速度	0.2752
能量利用率	0.2294

下载: 导出CSV

表 6 熵值法计算的各评价指标权重

Table 6. Weight of each evaluation index by entropy method

指标类型	权重
最大平均加速度	0.3328
最大瞬时加速度	0.3336
能量利用率	0.3336

下载: 导出CSV

表 7 各评价指标综合权重

Table 7. Comprehensive weight of each evaluation index

指标类型	权重
最大平均加速度	0.4948
最大瞬时加速度	0.2755
能量利用率	0.2297

下载: 导出CSV

参考文献(20)

[1]	JACOBSEN K, TYRELL D. Impact test of a crash-energy management passenger rail car[C]//2004 ASME/IEEE Joint Rail Conference. Baltimore: ASME, 2004: 19-26.
[2]	TYRELL D, JACOBSEN K, PARENT D. Preparations for a train-to-train impact test of crash-energy management passenger rail equipment[C]//2005 Joint Rail Conference. Pueblo: ASME, 2005: 107- 116.
[3]	CAROLAN M, MUHLANGER M, PERLMAN B, et al. Occupied volume integrity testing: elastic test results and analyses[C]//ASME 2011 Rail Transportation Division Fall Technical Conference. Minneapolis: ASME, 2011: 47-57.
[4]	CAROLAN M, PERLMAN B, TYRELL D. Crippling test of a budd pioneer passenger car[C]//ASME 2012 Joint Rail Conference. Philadelphia: ASME, 2012: 225-235.
[5]	WALDECK H. ICE动力车碰撞仿真[J]. 电力牵引快报,1996(4): 26-32. WALDECK H. ICE power car crash simulation[J]. Electric Traction Report, 1996(4): 26-32.
[6]	TYRELL D,张钟浩. 机车抗碰撞能力改进研究[J]. 国外内燃机车,2002(3): 23-31. TYRELL D, ZHANG Zhonghao. Research on improvement of anti-collision ability of locomotive[J]. Foreign Diesel Locomotive, 2002(3): 23-31.
[7]	Rail Track PLC. Structural requirements for railway vehicle: GM/RT—2100[S]. London: Safety and Standards Directo Rate, 1994.
[8]	LU G. Energy absorption requirement for crashworthy vehicles[J]. Proceedings of the Institution of Mechanical Engineers,Part F:Journal of Rail and Rapid Transit, 2002, 216(1): 31-39.
[9]	杨宝柱,肖守讷,杨超. 能量吸收方案对列车碰撞响应的影响[J]. 城市轨道交通研究,2018(3): 54-57. YANG Baozhu, XIAO Shoune, YANG Chao. The effect of energy absorbing schemes on train crash responses[J]. Urban Mass Transit, 2018(3): 54-57.
[10]	袁成标. 基于列车碰撞动力学的纵向碰撞建模及能量管理研究[D]. 成都: 西南交通大学, 2018.
[11]	郭亚军. 综合评价理论、方法及应用[M]. 北京: 科学出版社, 2007: 16-17, 45-47, 73.
[12]	王文华,黄仁乐,皱保平,等. 基于序关系分析法的配电网风险评估方法[J]. 计算机与数字工程,2018,46(8): 1516-1519,1552. doi: 10.3969/j.issn.1672-9722.2018.08.007 WANG Wenhua, HUANG Renle, ZHOU Baoping, et al. Based on-sequence relation analysis risk assessment method of distribution network[J]. Computer and Digital Engineering, 2018, 46(8): 1516-1519,1552. doi: 10.3969/j.issn.1672-9722.2018.08.007
[13]	李航,李根,张帅阳,等. 基于序关系分析法的车辆主动悬架控制器仿真研究[J]. 中国工程机械学报,2018,16(3): 231-236. LI Hang, LI Gen, ZHANG Shuaiyang, et al. Simulation research on vehicle active suspension controller based on the order relation analysis method[J]. Chinese Journal of Construction Machinery, 2018, 16(3): 231-236.
[14]	梁东,李远富,樊敏. 基于区间数距离的铁路线路方案决策模型与方法[J]. 西南交通大学学报,2019,54(4): 823-830. LIANG Dong , LI Yuanfu , FAN Min. Decision making model and method based on distance measure between interval numbers in railway location[J]. Journal of Southwest Jiaotong University, 2019, 54(4): 823-830.
[15]	徐健,杜贞栋,林洪孝,等. 基于序关系分析法的节水型社会评价指标权重的确定[J]. 水电能源科学,2014,32(10): 132-134. XU Jian, DU Zhendong, LIN Hongxiao, et al. Confirmation of index weight of water-saving social construction based on order relation analysis method[J]. Water Resources and Power, 2014, 32(10): 132-134.
[16]	高瑾瑾,郑源,李涧鸣. 基于改进序关系分析法的抽水蓄能电站经济效益综合评价[J]. 长江科学院院报,2018,35(4): 137-142. doi: 10.11988/ckyyb.20161372 GAO Jinjin, ZHENG Yuan, LI Jianming. Comprehensive evaluation of economic benefits of pumped storage power plants by improved order relation analysis method[J]. Journal of Yangtze River Scientific Research Institute, 2018, 35(4): 137-142. doi: 10.11988/ckyyb.20161372
[17]	梁晨. 子午线轮胎综合接地性能评价体系与方法研究[D]. 南京: 江苏大学, 2013.
[18]	余建英, 何旭宏. 数据统计分析与SPSS应用[M]. 北京: 人民邮电出版社, 2003: 87.
[19]	王学军,郭亚军,兰天. 构造一致性判断矩阵的序关系分析法[J]. 东北大学学报,2006,27(1): 115-118. WANG Xuejun, GUO Yajun, LAN Tian. Rank correlation analysis of formation of consistent judgment matrix[J]. Journal of Northeastern University, 2006, 27(1): 115-118.
[20]	刘祎. 多功能包装机设计质量综合评价系统研究[D]. 太原: 中北大学, 2016.

施引文献

期刊类型引用(3)

1.	冯国江，齐红瑞，范晓鹏，朱涛，谷鹏，肖守讷. 基于某特殊铁路标准的电力机车耐撞性设计方法. 铁道车辆. 2024(03): 117-123 . 百度学术
2.	Yanwen Liu，Bing Yang，Shoune Xiao，Tao Zhu，Guangwu Yang，Ruixian Xiu. Research on the Influence of Multiple Parameters on the Responses of a B-type Subway Train. Chinese Journal of Mechanical Engineering. 2022(04): 329-346 . 必应学术
3.	吕天一，肖守讷，朱涛，张敬科，王小瑞，李玉如，阳光武，杨冰. 采用多目标进化算法的列车碰撞能量配置方案优化. 中国机械工程. 2021(18): 2262-2267 . 百度学术

其他类型引用(7)

附加材料(0)

访问统计

点击查看大图

图(4) / 表(7)

计量

文章访问数: 590
HTML全文浏览量: 255
PDF下载量: 9
被引次数: 10

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

列车一维碰撞能量管理的综合评价模型

doi: 10.3969/j.issn.0258-2724.20200034

作者简介: 张敬科（1992—），男，博士研究生，研究方向为机车车辆碰撞动力学，E-mail：jkzhang@my.swjtu.edu.cn

计量

出版历程