中国高速列车车轮多边形磨耗特征分析

王鹏; 陶功权; 杨晓璇; 谢晨希; 李伟; 温泽峰

doi:10.3969/j.issn.0258-2724.20210777

中国高速列车车轮多边形磨耗特征分析

doi: 10.3969/j.issn.0258-2724.20210777

西南交通大学轨道交通运载全国重点实验室，四川成都 610031

基金项目: 国家自然科学基金（1734201, 51805450, 52002342）；中央高校基本科研业务费专项资金（2682020CX52）

详细信息

作者简介:
王鹏（1990―），男，博士研究生，研究方向为轮轨关系，E-mail：15198285012@163.com

通讯作者:
温泽峰（1976―），男，研究员，研究方向为轮轨关系及振动噪声，E-mail：zfwen@swjtu.edu.cn

中图分类号: U270.1
计量
- 文章访问数: 1010
- HTML全文浏览量: 265
- PDF下载量: 167
- 被引次数: 7
出版历程
- 收稿日期: 2021-10-12
- 修回日期: 2021-11-23
- 网络出版日期: 2023-01-03
- 刊出日期: 2021-11-25

Analysis of Polygonal Wear Characteristics of Chinese High-Speed Train Wheels

State Key Laboratory of Rail Transit Vehicle System, Southwest Jiaotong University, Chengdu 610031, China

摘要

摘要:
车轮失圆问题广泛存在于我国高速列车，对列车乘坐舒适性和运行安全性有显著影响. 从2011年至2020年，测试了12条高铁线路中9种型号高速动车组的车轮不圆度，包括200、250、300、350 km/h 4种运营速度，共3.05万个车轮；对车轮不圆度测试数据进行特征分析，掌握我国高速动车组车轮多边形磨耗的发展规律；分析影响车轮多边形发展的关键因素，包括车辆轴距、轨道结构和研磨子修形. 结果表明：高速动车组车轮存在10～30阶多边形磨耗，多边形波长为90～288 mm，且在100～178 mm波长范围的多边形磨耗最为严重；车辆轴距、扣件类型和环境温度与车轮多边形磨耗形成密切相关，通过改善研磨子和车轮踏面匹配关系，保证踏面横向和圆周处于良好的磨耗状态，使高阶车轮多边形粗糙度水平最大下降60%.
- 高速列车 /
- 车轮 /
- 磨耗 /
- 多边形 /
- 研磨子 /
- 轴距 /
- 扣件刚度
Abstract:
Wheel out-of-roundness (OOR), which is commonly observed on the wheels of Chinese high-speed trains, has a significant influence on the vehicle ride comfort and operation safety. From 2011 to 2020, 30500 wheels of nine types of high-speed trains were selected for tests and the wheel OOR was measured. The tested high-speed trains operated on 12 high-speed railway lines with different operating speeds, including operating speeds of 200, 250, 300, and 350 km/h. The characteristics of the wheel OOR were analyzed to determine the development rules for wheel polygonal wear of high-speed EMUs in China. The key factors that affect the wheel polygon wear were analyzed, including the wheelbase, track structure, and abrasive block. The test results show that the dominant harmonic orders of the wheel polygonal wear range from 10^th to 30^th order. The corresponding wavelengths of the 10^th–30^th-order polygonization are 90–288 mm, with a wheel polygonal wear of 100‒178 mm being the most severe. The analysis of the key factors that affect the wheel polygonal wear shows that the wheelbase, fastening stiffness, and ambient temperature are closely related to the formation of wheel polygonal wear. By improving the matching relationship between the abrasive block and the wheel tread, the transverse and circumferential wear is minimized and the roughness level of the high-order polygonal wear is reduced by 60% at most.
- high-speed train /
- wheels /
- wear /
- polygon /
- abrasive block /
- wheelbase /
- fastening stiffness

HTML全文

行驶舒适性与路面的平整程度、车辆状态、运行速度等方面密切相关，被广泛认为是路面服务性能评价的重要参考指标. 路面的不平整激励是引起汽车振动并导致舒适性下降的核心因素，道路表面横向、纵向、水平向的各种波长及幅值的复杂变形组合构成了路面的不平度^[1]. 作为激励源，不平度的连续变化会引起车辆的随机振动，局部突变则会引起车辆的瞬态振动，二者共同作用，通过车身悬架系统传递至人体，最终引起其心理和生理上的舒适性感受变化.

国际上用于评估行驶舒适性的指标众多，例如舒适度指数（ride number，RN）、国际平整度指数（international roughness index，IRI）、行驶质量指数（riding quality index，RQI）以及加权加速度均方根（weighed root-mean-square acceleration，WRMSA）等. 利用这些评价指标，现有文献研究集中在宏观层面的路段行驶质量和人车振动特性进行分析^[2-4]，但对于路面局部变形（如桥头差异沉降、坑槽、窨井沉陷）导致的冲击振动问题却鲜有涉及.

已有研究表明路面局部变形引起的冲击振动对行驶舒适性具有显著影响^[5-6]，但尚未明确其中的定量关系，部分文献也提出过路面中桥头差异沉降及窨井沉陷等局部变形的控制阈值^[7-9]，但是缺乏有效的理论支撑. 另有文献提出以较短长度作为行驶质量评价单元能够更加准确地描述路面的服务性能^[10]，但是较短的单元长度可能会造成路面波长截断现象^[2]，导致评价结果失真. 由于路面局部变形随机存在于路面中，且其导致的振动具有瞬时性特点，因此，难以通过现场试验获取变形激励与振动感受之间的定量关系.

鉴于此，本文从理论上研究局部变形对舒适性的影响. 首先，建立典型路面局部变形简化模型作为振动激励源，基于两自由度车辆强迫振动模型，从理论上推导出振动模型中簧上质量位移和加速度时间历程函数，再根据ISO 2631-1^[11]中的频率加权系数分布特点设计频率加权滤波器对竖向加速度进行滤波，求解不同工况下的均方根加速度，评估人体的振动感受，最后提出不同类型局部变形的评价标准，以期为路面精细化养护管理提供依据.

1. 路面局部变形纵断面形式假定

文献[2, 12]对上海城市道路进行详细调研发现，引起道路纵断面高程局部异常波动的不平整因素主要是窨井变形、桥头差异沉降、路面修补等，各类因素的统计结果如图1所示.

图 1 路面局部变形因素统计

Figure 1. Statistical factors of local deformation in pavement

下载: 全尺寸图片幻灯片

由图1可知，窨井损坏和修补（包括路面开挖回填）是最主要的变形形式，其次是桥头差异沉降. 但是，实际路面中路面修补一般面积较大，且对行车影响十分有限；路面开挖则一般为临时性养护管理；另外，由于轮胎存在包络特性，伸缩缝问题也不具显著影响. 为进一步研究变形的影响，本文根据其余常见局部变形的形态特征，将其分为3类：

1）错台型，描述同一断面的高程突变，如桥头差异沉降；

2）双断面型，描述路面规则的凹凸变形，行进中的轮胎会在极短时间内经过变形区域的两端，如窨井沉陷或拱起；

3）曲面型，描述路面较均匀的凸起，如路面拥包等.

建立激励模型需首先假定各类局部变形的函数形式，由于路面局部变形处纵断面高程一般具有突变特性，传统简谐函数难以描述其断面特点，因此对3类变形进行简化（如图2所示），以高差 $\Delta g$ 作为变形衡量指标，纵断面高程y(t)如式（1）～（3）所示.

1）错台型

$y\left( t \right) = \left\{ {\begin{array}{*{20}{l}} {0,\;\;t = 0,} \\ \Delta g,\;\;0 \lt t \leqslant l/v, \end{array}} \right.$

(1)

式中：l 为变形部分长度； t 为时间； v 为车速.

2）双断面型

$y\left( t \right) = \left\{ {\begin{array}{*{20}{l}} \begin{gathered} 0,\;\;t = 0, \hfill \\ \Delta g,\;\;\;0 \lt t \leqslant d/v, \hfill \\ \end{gathered} \\ {0,\;\;t \gt d/v,} \end{array}} \right.$

(2)

式中： $d$ 为凹凸部分直径或长度.

3）曲面型

$y\left( t \right) = \left\{ {\begin{array}{*{20}{l}} \begin{gathered} 0,\;\;t = 0,\;\;\;\;\;\;\; \hfill \\ - \dfrac{{4\Delta gv}}{{{l^2}}}\left( {v{t^2} - lt} \right),\;\;\;0 \lt t \leqslant l/v, \hfill \\ \end{gathered} \\ {0,\;\;t \gt l/ v\;.\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;} \end{array}} \right.$

(3)

图 2 典型路面局部变形纵断面示意

Figure 2. Profile diagrams of typical local deformation

下载: 全尺寸图片幻灯片

2. 人体竖向位移及竖向加速度推导

振动舒适性评价一般以加速度为基础指标，已有研究中加速度多为现场实测值，不同车辆之间难以对比. 本文通过两自由度车辆强迫振动模型的求解，从理论上推导人体的竖向位移和加速度响应.

2.1 两自由度车辆强迫振动模型

为统一各国的平整度度量标准，世界银行曾在巴西进行了大量平整度试验（international road roughness experiment），试验中首次构想出两自由度1/4理想车模型^[13]，由于理想车的运动方程恒定，固定车辆参数后，相当于利用同一辆车在任意地方行驶，因此，其测试方法具有很好的移植性. 1/4车模型如图3所示，图中： ${m_{\text{s}}}$ 为簧上质量， ${m_{\text{t}}}$ 为簧下质量， ${K_{\text{s}}}$ 为悬架刚度系数， ${K_{\text{t}}}$ 为轮胎刚度系数， ${C_{\text{s}}}$ 为悬架阻尼系数， ${Z_{\text{s}}}$ 和 ${Z_{\text{t}}}$ 分别为簧上质量和簧下质量的绝对位移.

图 3 1/4车模型

Figure 3. Quarter car model

下载: 全尺寸图片幻灯片

运动方程为

${m_{\text{s}}}{{{\mathop Z\limits^{ \bullet \bullet }}_{\rm{s}}}} + {C_{\text{s}}}\left( {{{\mathop Z\limits^{ \bullet } }_{\rm{s}}} - {{\mathop Z\limits^{ \bullet }}_{\rm{t}}}} \right) + {K_{\text{s}}}\left( {{Z_{\text{s}}} - {Z_{\text{t}}}} \right) = 0 \text{，}$

(4)

${m_{\rm{t}}}{{\mathop Z\limits^{ \bullet \bullet }}_{\rm{t}}} + {C_{\rm{s}}}\left( {{{\mathop Z\limits^{ \bullet }}_{\rm{t}}} - {{\mathop Z\limits^{ \bullet}}_{\rm{s}}}} \right) + {K_{\rm{s}}}\left( {{Z_{\rm{t}}} - {Z_{\rm{s}}}} \right) + {K_{\rm{t}}}({Z_{\rm{t}}} - y) = 0.$

(5)

将式（4）、（5）等式左右两侧同时除以m_s，得

${{{\mathop Z\limits^{ \bullet \bullet} }_{\rm{s}}}} + c\left( {{{\mathop Z\limits^{ \bullet } }_{\text{s}}} - {{\mathop Z\limits^ \bullet }_{\text{t}}}} \right) + {k_{\text{2}}}\left( {{Z_{\text{s}}} - {Z_{\text{t}}}} \right) = 0 \text{，}$

(6)

$u{{{\mathop Z\limits^{ \bullet \bullet }}_{\rm{t}}}} + c\left( {{{\mathop Z\limits^ \bullet }_{\text{t}}} - {{\mathop Z\limits^ \bullet }_{\text{s}}}} \right) + {k_{\text{2}}}\left( {{Z_{\text{t}}} - {Z_{\text{s}}}} \right) + {k_{\text{1}}}({Z_t} - y) = 0 \text{，}$

(7)

式中：系数 $c = {C_{\text{s}}}/{m_{\text{s}}} = 6.00\;{{\rm{s}}^{ - 1}}；$ 系数 $u = {m_{\text{t}}}/{m_{\text{s}}} = 0.15$ ；系数 ${k_1} = {K_{\text{t}}}/{m_{\text{s}}} = 653\;{{\rm{s}}^{ - 2}}；$ 系数 ${k_2} = {K_{\text{s}}}/{m_{\text{s}}} = 63.3\;{{\rm{s}}^{ - 2}}$ .

2.2 任意路面变形激励下人体竖向位移的时域函数

路面局部变形对两自由度车辆振动系统的激励不具有简谐特性和周期特性，在这种激振作用下，系统通常仅有瞬态振动，故为求得加速度的时域函数，可以采用卷积积分法^[14]求解簧上质量的加速度响应，具体过程如下：

1）求解1/4车振动系统对于单位脉冲激励 $\delta \left( t \right)$ 的响应函数；

2）利用卷积积分法求解路面局部变形激励作用下的系统响应.

簧上质量的竖向位移可近似认为是人体的竖向位移. 由于运动方程存在耦合，故需要对方程进行解耦，运动方程如式（8）所示.

${\boldsymbol{M}}\mathop {\boldsymbol{Z}}\limits^{ \bullet \bullet } + {\boldsymbol{C}}\mathop {\boldsymbol{Z}}\limits^ \bullet + {\boldsymbol{KZ}} = {\boldsymbol{Y}}，$

(8)

式中：M、C、K、Z、Y分别为质量矩阵、阻尼矩阵、刚度矩阵、位移矩阵、输入矩阵，

$\begin{split}& {\boldsymbol{M}} = \left[ {\begin{array}{*{20}{c}} {{m_{\text{s}}}}&0 \\ 0&{{m_{\text{t}}}} \end{array}} \right],\;\;\; \\ & {\boldsymbol{C}} = \left[ {\begin{array}{*{20}{c}} {{C_{\text{s}}}}&{ - {C_{\text{s}}}} \\ { - {C_{\text{s}}}}&{{C_{\text{s}}}} \end{array}} \right], \\& {\boldsymbol{K}} = \left[ {\begin{array}{*{20}{c}} {{K_{\text{s}}}}&{ - {K_{\text{s}}}} \\ { - {K_{\text{s}}}}&{{K_{\text{s}}} + {K_{\text{t}}}} \end{array}} \right],\;\; \\& {\boldsymbol{Z}} = \left[ {\begin{array}{*{20}{c}} {{Z_{\text{s}}}} & {{Z_{\text{t}}}} \end{array}} \right]^{\text{T}},\; \\& \;{\boldsymbol{Y}} = \left[ {\begin{array}{*{20}{c}} 0 & {{K_t}y} \end{array}} \right]^{\text{T}}. \\ \end{split}$

运动方程中阻尼矩阵无法通过实模态变换实现对角化，因此采用复模态变换求解^[15]. 首先，令式（8） ${\boldsymbol{Y}} =\bf{ 0}$ 得到自由振动方程，设特解为 ${\boldsymbol{\mu }}{{\rm{e}}^{{\boldsymbol{\lambda }}t}}$ ， ${\boldsymbol{\mu }}$ 、λ分别为特征向量和特征值，可求解振动系统的特征方程为

$\left| {{{\boldsymbol{\lambda }}^2}{\boldsymbol{M}} + {\boldsymbol{\lambda C}} + {\boldsymbol{K}}} \right| = {\boldsymbol{0}}.$

(9)

二自由度系统存在4个特征值，且为2对共轭复数. 将1/4车模型参数代入特征方程可得4个特征值分别为

$\left\{ {\begin{array}{*{20}{l}} {{\lambda _1} = - 20.4 + 64.4{\rm{i}},\;}\\ {{\lambda _2} = - 2.6 + 7.3{\rm{i}},}\\ {{\lambda _3} = - 20.4 - 64.4{\rm{i}},\;}\\ {{\lambda _4} = - 2.6 - 7.3{\rm{i}}.} \end{array}} \right.$

(10)

代入自由振动方程，对应的特征向量构成复模态矩阵 ${\boldsymbol{\mu }} = \left[ {\begin{array}{*{1}{c}}{{{\boldsymbol{\mu }}_1}}\;\;{{{\boldsymbol{\mu }}_2}}\;\;{{{\boldsymbol{\mu }}_3}}\;\;{{{\boldsymbol{\mu }}_4}}\end{array}} \right]$ ，其中：

$\left\{ \begin{array}{*{20}{l}} {{\boldsymbol{\mu }}_1} = {\left[ {\begin{array}{*{20}{c}} {{\mu _{11}}}\;\;{{\mu _{12}}} \end{array}} \right]^{\rm{T}}} = {\left[ {\begin{array}{*{20}{c}} 1\;\;{ - 4.2 + 10.5{\rm{i}}} \end{array}} \right]^{\rm{T}}}, \hfill \\ {{\boldsymbol{\mu }}_2} = {\left[ {\begin{array}{*{20}{c}} {{\mu _{21}}}\;\;{{\mu _{22}}} \end{array}} \right]^{\rm{T}}} = {\left[ {\begin{array}{*{20}{c}} 1\;\;{0.074 + 0.054{\rm{i}}} \end{array}} \right]^{\rm{T}}}, \hfill \\ {{\boldsymbol{\mu }}_3} = {\left[ {\begin{array}{*{20}{c}} {{\mu _{31}}}\;\;{{\mu _{32}}} \end{array}} \right]^{\rm{T}}} = {\left[ {\begin{array}{*{20}{c}} 1\;\;{ - 4.2 - 10.5{\rm{i}}} \end{array}} \right]^{\rm{T}}}, \hfill \\ {{\boldsymbol{\mu }}_4} = {\left[ {\begin{array}{*{20}{c}} {{\mu _{41}}}\;\;{{\mu _{42}}} \end{array}} \right]^{\rm{T}}} = {\left[ {\begin{array}{*{20}{c}} 1\;\;{0.074 - 0.054{\rm{i}}} \end{array}} \right]^{\rm{T}}}, \hfill \\ \end{array} \right.$

(11)

由于上述复模态矩阵无法直接变换，因此这里引入状态变量 ${\boldsymbol{x}}$ 对方程解耦^[15]，补充恒等式 ${\boldsymbol{M}}\mathop {\boldsymbol{Z}}\limits^ \bullet - {\boldsymbol{M}}\mathop {\boldsymbol{Z}}\limits^ \bullet = {\bf{0}}$ ，将运动方程转换为状态方程：

$\mathop {\widetilde{M}}\limits \mathop {\boldsymbol{x}}\limits^ \bullet + \mathop {\widetilde{K}}\limits {\boldsymbol{x}} = \mathop {\widetilde{Y}}\limits \text{，}$

(12)

式中：

${\boldsymbol{x}} = {\left[ {\begin{array}{*{20}{c}} {\mathop {\boldsymbol{Z}}\limits^ \bullet }&{\boldsymbol{Z}} \end{array}} \right]^{\rm{T}}};\;\;\mathop {\boldsymbol{\widetilde{M}}}\limits = \left[ {\begin{array}{*{20}{c}} {\bf{0}}&{\boldsymbol{M}} \\ {\boldsymbol{M}}&{\boldsymbol{C}} \end{array}} \right]；$

$\mathop {\boldsymbol{\widetilde{K}}}\limits = \left[ {\begin{array}{*{20}{c}} { - {\boldsymbol{M}}}&{\bf{0}} \\ {\bf{0}}&{\boldsymbol{K}} \end{array}} \right];\;\;\mathop {\boldsymbol{\widetilde{Y}}}\limits = {\left[ {\begin{array}{*{20}{c}} {\bf{0}}&{\boldsymbol{Y}} \end{array}} \right]^{\rm{T}}}.$

特征向量构成复模态矩阵：

${\boldsymbol{V}} = \left[ {\begin{array}{*{20}{c}} {\begin{array}{*{20}{c}} {{\lambda _1}{{\boldsymbol{\mu }}_1}} \\ {{{\boldsymbol{\mu }}_1}} \end{array}}&{\begin{array}{*{20}{c}} {{\lambda _2}{{\boldsymbol{\mu }}_2}} \\ {{{\boldsymbol{\mu }}_2}} \end{array}}&{\begin{array}{*{20}{c}} {{\lambda _3}{{\boldsymbol{\mu }}_3}} \\ {{{\boldsymbol{\mu }}_3}} \end{array}}&{\begin{array}{*{20}{c}} {{\lambda _4}{{\boldsymbol{\mu }}_4}} \\ {{{\boldsymbol{\mu }}_4}} \end{array}} \end{array}} \right] \text{.}$

(13)

${\boldsymbol{V}}$ 具有关于 ${\boldsymbol{\widetilde{M}}}$ 和 $\mathop {\boldsymbol{\widetilde{K}}}\limits$ 的加权正交性，满足

${\lambda _i}\mathop {\boldsymbol{\widetilde{M}}}\limits {{\boldsymbol{V}}_i} + \mathop {\boldsymbol{\widetilde{K}}}\limits {{\boldsymbol{V}}_i} = 0,\;\;\;\;i = 1,2,3,4,$

(14)

式中： ${{\boldsymbol{V}}_i}$ 为矩阵V的第i个列向量.

代入 $\mathop {\boldsymbol{\widetilde{M}}}\limits$ 、 $\mathop {\boldsymbol{\widetilde{K}}}\limits$ ，并作复模态变换 ${\boldsymbol{x}} = {\boldsymbol{Vp}}$ ，其中p为模态变换响应列阵，代入式（12）可对振动方程解耦：

$\begin{split}& {{\rm{diag}}\left( {{m_1},{m_2},{m_3},{m_4}} \right)\mathop {p}\limits^ \bullet + } \;{\rm{diag}} ( - {\lambda _1}{m_1}, - {\lambda _2}{m_2},\\&\quad - {\lambda _3}{m_3}, - {\lambda _4}{m_4} ) {\boldsymbol{p}} = {{\boldsymbol{V}}^{\rm{T}}}{\tilde{\boldsymbol Y}} \text{，}\end{split}$

(15)

式中： ${m_i} = {{\boldsymbol{V}}_i}^{\rm{T}}\mathop {\boldsymbol{\widetilde{M}}}\limits {{\boldsymbol{V}}_i}$

将式（10）、（11）代入式（15），求得

$\left\{ {\begin{array}{*{20}{l}} {{m_1} = \left( {1\;730 - 1\;774{\rm{i}}} \right){m_{\rm{s}}}},\\ {{m_2} = {m_{\rm{s}}}14{\rm{i}}},\\ {{m_3} = \left( {1\;730 + 1\;774{\rm{i}}} \right){m_{\rm{s}}}},\\ {{m_4} = - {m_{\rm{s}}}14{\rm{i}}}. \end{array}} \right.$

(16)

当初始输入为单位脉冲信号，式（15）化简为

$\mathop {{p_i}}\limits^ \bullet - {\lambda _i}{p_i} = \frac{{{\mu _{i2}}}}{{{m_i}}}\delta \left( t \right)\; \text{.}$

(17)

式中： ${p_i}$ 为响应列阵中第i个值.

式(17)等式两边取跨零积分，得

$\int_{0 - \varepsilon }^{0 + \varepsilon } {\left( {\mathop {{p_i}}\limits^ \bullet - {\lambda _i}{p_i}} \right)} {\rm{d}}t = \frac{{{\mu _{i2}}}}{{{m_i}}}\int_{0 - \varepsilon }^{0 + \varepsilon } {\delta \left( t \right){\rm{d}}t} \text{，}$

(18)

式中： $\varepsilon$ 为极小量.

得到系统的单位脉冲激励响应为

${h_i}(t) = \frac{{{\mu _{i2}}}}{{{m_i}}}{{\rm{e}}^{{\lambda _i}t}} \text{.}$

(19)

利用卷积积分方法，可求得在任意激励 ${K_{\text{t}}}y(t)$ 下的响应为

${p_i} = \int_0^t {h\left( {t - \tau } \right)} {K_{\text{t}}}y\left( \tau \right){\rm{d}}\tau = \dfrac{{{K_{\text{t}}}{\mu _{i2}}}}{{{m_i}}}\int_0^t {{{\rm{e}}^{{\lambda _i}\left( {t - \tau } \right)}}y\left( \tau \right){\rm{d}}\tau } ,$

(20)

代入 ${\boldsymbol{x}} = {\boldsymbol{Vp}}$ 得到

$\begin{split}& \left[ {\begin{array}{*{20}{c}} {\mathop {{\rm{ }}{\boldsymbol{Z}}}\limits^ \bullet }\\ {\boldsymbol{Z}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {\begin{array}{*{20}{c}} {{\lambda _1}{{\boldsymbol{\mu }}_1}}\\ {{{\boldsymbol{\mu }}_1}} \end{array}\;\;\begin{array}{*{20}{c}} {{\lambda _2}{{\boldsymbol{\mu }}_2}}\\ {{{\boldsymbol{\mu }}_2}} \end{array}\;\;\begin{array}{*{20}{c}} {{\lambda _3}{{\boldsymbol{\mu }}_3}}\\ {{{\boldsymbol{\mu }}_3}} \end{array}\;\;\begin{array}{*{20}{c}} {{\lambda _4}{{\boldsymbol{\mu }}_4}} \\ {{{\boldsymbol{\mu }}_4}} \end{array}} \end{array}} \right]\;\times\\&\quad {\left[ {\begin{array}{*{20}{c}} {{p_1}\;\;{p_2}\;\;{p_3}\;\;{p_4}} \end{array}} \right]^{\rm{T}}}. \end{split}$

(21)

可得

${\boldsymbol{Z}} = \left[ {\begin{array}{*{20}{c}} {\displaystyle \sum\limits_{i = 1}^4 {{\mu _{i1}}{p_i}} }\\ {\displaystyle \sum\limits_{i = 1}^4 {{\mu _{i2}}{p_i}} } \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {\displaystyle \sum\limits_{i = 1}^4 {\dfrac{{{K_{\rm{t}}}{\mu _{i1}}{\mu _{i2}}}}{{{m_i}}}\int\nolimits_0^t {{{\rm{e}}^{{\lambda _i}\left( {t - \tau } \right)}}y\left( \tau \right){\rm{d}}\tau } } }\\ {\displaystyle \sum\limits_{i = 1}^4 {\dfrac{{{K_{\rm{t}}}{\mu _{i2}}^2}}{{{m_i}}}\int\nolimits_0^t {{{\rm{e}}^{{\lambda _i}\left( {t - \tau } \right)}}y\left( \tau \right){\rm{d}}\tau } } } \end{array}} \right]，$

(22)

${Z_{\text{s}}}$ 即为 $y\left( t \right)$ 激励下的人体竖向位移.

2.3 局部变形激励下的人体竖向加速度

由式（22）可知：人体竖向位移的变化振动与系统的参数及路面激励相关，1/4车模型参数固定后，位移则是路面激励的函数. 对位移求二阶导数可以得到加速度函数^[16]，下面分别推导3类典型局部变形激励下的加速度响应.

2.3.1 错台型

车辆在经过错台型变形会经历两个阶段，以桥头跳车为例，第一阶段为驶上但是未离开桥面，第二阶段为驶离桥面. 车辆经过两处错台，引起两次跳车. 但是，实际桥梁具有一定长度，汽车经过第一处错台导致的加速度响应会在短时间内衰减为0，而此时车辆仍未驶下桥面，故每次错台视作为一处变形. 考虑驶上错台的加速度时间历程，即 $0 \leqslant t \leqslant l/ v$ ，联立式（1）、（22）求解得

${Z_{\text{s}}} = \left\{ \begin{gathered} 0,\;\;\;t = 0, \hfill \\ {K_{\text{t}}}\Delta g\sum\limits_{i = 1}^4 {\frac{{{\mu _{i2}}}}{{{\lambda _i}{m_i}}}\left( {{{\rm{e}}^{{\lambda _i}t}} - 1} \right)} ,\;\;\;0 \lt t \leqslant l /v, \hfill \\ \end{gathered} \right.$

(23)

对其求二阶导数，得到竖向加速度为

${a_{\text{s}}} = \left\{ \begin{gathered} 0,\;\;\;t = 0, \hfill \\ {K_{\text{t}}}\Delta g\sum\limits_{i = 1}^4 {\frac{{{\lambda _i}{\mu _{i2}}}}{{{m_i}}}{{\rm{e}}^{{\lambda _i}t}}} ,\;\;\;0 \lt t \leqslant l/ v. \hfill \\ \end{gathered} \right.$

(24)

2.3.2 双断面型

双断面型变形的长度一般较小，例如窨井的直径一般仅为0.5 m左右，导致车辆驶过瞬间会连续经过变形两端，且间隔时间极短，两次振动响应会发生叠加，所以车轮驶上驶下变形的两次冲击都需要考虑在内，加速度响应为分段函数，联立式（2）、（22）求解，结果如式（25）所示.

${Z_{\text{s}}} = \left\{ {\begin{array}{*{20}{l}} \begin{gathered} 0,\;\;\;t = 0, \hfill \\ {K_{\text{t}}}\Delta g \displaystyle \sum\limits_{i = 1}^4 {\dfrac{{{\mu _{i2}}}}{{{\lambda _i}{m_i}}}\left( {{{\rm{e}}^{{\lambda _i}t}} - 1} \right),\;\;\;0 \lt t \lt d/v,} \hfill \\ \end{gathered} \\ {{K_{\text{t}}}\Delta g \displaystyle \sum\limits_{i = 1}^4 {\dfrac{{{\mu _{i2}}}}{{{\lambda _i}{m_i}}}\left( {{{\rm{e}}^{{\lambda _i}t}} - {{\rm{e}}^{{\lambda _i}\left( {t - d/v} \right)}}} \right),\;\;\;t \geqslant d/v.} } \end{array}} \right.$

(25)

对式（25）求二阶导，得到竖向加速度为

${a_{\text{s}}} = \left\{ {\begin{array}{*{20}{l}} \begin{gathered} \;0,\;\;\;t = 0, \hfill \\ {K_{\text{t}}}\Delta g \displaystyle \sum\limits_{i = 1}^4 {\dfrac{{{\lambda _i}{\mu _{i2}}}}{{{m_i}}}{{\rm{e}}^{{\lambda _i}t}},\;\;\;0 \lt t \lt d/v,} \hfill \\ \end{gathered} \\ {{K_{\text{t}}}\Delta g \displaystyle \sum\limits_{i = 1}^4 {\dfrac{{{\lambda _i}{\mu _{i2}}}}{{{m_i}}}\left( {{{\rm{e}}^{{\lambda _i}t}} - {{\rm{e}}^{{\lambda _i}\left( {t - d/v} \right)}}} \right),\;\;\;t \geqslant d/v.} } \end{array}} \right.$

(26)

2.3.3 曲面型

此类变形断面为平滑曲线，在变形长度较大时，振动响应会减弱. 联立式（3）、（22）求解，结果如式（27）所示. 同理可得竖向加速度结果如式（28）所示.

$\qquad \qquad \begin{split} {Z_{\rm{s}}} = \left\{ \begin{array}{l} \;0\;,\;\;\;t = 0,\;\;\\ \displaystyle \sum\limits_{i = 1}^4 {\dfrac{{4\Delta gv{K_{\rm{t}}}{\mu _{i2}}}}{{{\lambda _i}^3{m_i}{l^2}}}\left( \begin{array}{l} v{\lambda _i}^2{t^2} - 2v{{\rm{e}}^{{\lambda _i}t}} + {\lambda _i}l{{\rm{e}}^{{\lambda _i}t}} - {\lambda _i}^2lt + 2v{\lambda _i}t - {\lambda _i}l + 2v \end{array} \right),} \;\;\;0 < t < l/v,\\ \displaystyle \sum\limits_{i = 1}^4 {\dfrac{{4\Delta gv{K_{\rm{t}}}{\mu _{i2}}}}{{{\lambda _i}^3{m_i}{l^2}}}{{\rm{e}}^{{\lambda _i}\left( {t - l/v} \right)}}\left( \begin{array}{l} {\lambda _i}l{{\rm{e}}^{{\lambda _i}l/v}} - 2v{{\rm{e}}^{{\lambda _i}l/v}} + {\lambda _i}l + 2v \end{array} \right),} \;\;\;t \geqslant l/v. \end{array} \right. \end{split}$

(27)

${a_{\text{s}}} = \left\{ {\begin{array}{*{20}{l}} \begin{gathered} \;0,\;\;\;t = 0,\;\; \hfill \\ \sum\limits_{i = 1}^4 {\frac{{4\Delta gv{K_{\text{t}}}{\mu _{i2}}}}{{{\lambda _i}{m_i}{l^2}}}\left( {\left( {{\lambda _i}l - 2v} \right){{\rm{e}}^{{\lambda _i}t}} + 2v} \right),\;\;\;0 \lt t \lt l/v,} \hfill \\ \end{gathered} \\ \begin{gathered} \sum\limits_{i = 1}^4 {\frac{{4\Delta gv{K_{\text{t}}}{\mu _{i2}}}}{{{\lambda _i}{m_i}{l^2}}}{{\rm{e}}^{{\lambda _i}\left( {t - l/v} \right)}}\left( {{\lambda _i}l{{\rm{e}}^{{\lambda _i}l/v}} + {\lambda _i}l - 2v{{\rm{e}}^{{\lambda _i}l/v}} + 2v} \right),} \hfill \\ \;\;\;\;\;\;\;t \geqslant l/v.\;\; \hfill \\ \end{gathered} \end{array}} \right.$

(28)

3. 基于频率加权的振动舒适性评价

局部变形激励导致的人车振动为瞬时冲击，不具周期性特征，ISO 2631-1^[11]推荐采用运行均方根法来评价这种路面偶然激励导致的振动感受，其形式为

${a_w}\left( {{t_0}} \right) = \sqrt {\frac{1}{{{\tau _0}}}\int_{{t_0} - {\tau _0}}^{{t_0}} {a_{\rm{w}}^2\left( t \right){\rm{d}}t} } \text{，}$

(29)

式中： ${t_0}$ 为瞬时激励发生时刻； ${\tau _0}$ 为持续时间； ${a_{\rm{w}}}\left( t \right)$ 为瞬时频率加权加速度.

欲求 ${a_{\rm{w}}}\left( t \right)$ ，需要对竖向加速度的时域值进行频率加权滤波. 由第2节可知，竖向加速度时域值与行车速度、变形长度、变形高差相关，因本文主要研究不同类型局部变形引起的人体振动，故先统一车速和变形长度两个变量. 根据现行城市道路规范^[17]，道路设计速度为20～60 km/h，本文车速取40 km/h；实际路面中双断面型变形多表现为窨井的凸起或沉陷损坏，窨井的直径一般为0.6～0.7 m^[2]，本文取0.6 m.

3.1 频率加权滤波器设计

人体对不同频率的振动敏感程度不同，一般认为0.5～80.0 Hz是振动舒适性评价的考虑范围，竖向振动频率加权函数^[18]为

${W_{\text{k}}}\left( f \right) = \left\{ \begin{gathered} 0.5,\;\;\;0.5 \lt f \leqslant 2.0, \hfill \\ {f/4},\;\;\;2.0 \lt f \leqslant 4.0, \hfill \\ 1.0,\;\;\;4.0 \lt f \leqslant 12.5, \hfill \\ {12.5/f},\;\;\;12.5 \lt f \leqslant 80.0, \hfill \\ \end{gathered} \right.$

(30)

式中：f 为竖向振动频率.

${W_{\text{k}}}$ 在低频段和高频段均呈现抑制性，而在中间频段则具有较好的通过性，所以可以设计带通滤波器进行频率加权滤波. 因运行均方根计算对加速度相位要求不高，故采用IIR （infinite impulse response）数字滤波器设计，并借助MATLAB信号处理工具箱完成. 巴特沃斯滤波器原型具有通带内最大平坦和随频率升高而单调递减特性，与标准中的 ${W_{\text{k}}}$ 具有较高契合度，因此本文采用巴特沃斯带通滤波器.

数字滤波器通过间接法设计，步骤如下：

步骤1　确定数字滤波器技术指标；

步骤2　将数字滤波器技术指标转换为模拟低通滤波器技术指标；

步骤3　设计模拟低通滤波器（滤波器原型）；

步骤4　通过频率转换，将模拟低通滤波器转变为模拟带通等滤波器；

步骤5　模拟滤波器数字化.

由MATLAB中buttord函数先确定滤波器阶数，再由buttap函数设计出滤波器原型，并借助lp2bp函数实现模拟低通至模拟带通的转换，最后由bilinear函数完成数字化.

根据 ${W_{\text{k}}}$ 渐近线形式，频率加权函数在4.0～12.5 Hz，权重保持为1，考虑巴特沃斯滤波器的过渡特性，试算后取通带范围为2.5～18.0 Hz，通过上述过程可最终设计出频率加权函数滤波器，本文同时设计了一个对比滤波器，两者的区别在于阻带最小衰减量不同. 设计滤波器的幅频特性与ISO 2631-1^[11]中的频率加权函数对比如图4所示.

图 4 竖向加速度的频率加权曲线

Figure 4. Frequency weighting curves of vertical acceleration

下载: 全尺寸图片幻灯片

图4可见，设计滤波器在高频段具有一定误差，但在25.0 Hz以下人体最敏感的低频段吻合度高，且在80.0 Hz之后迅速衰减为0，而对比滤波器则与之相反. 考虑到路面局部变形的激振频率主要为低频振动（如图5），高频段的权重实际意义不大，即保证通带频率范围的特性更为关键，因此，本文设计的带通滤波器适用性更高.

图 5 竖向加速度频谱

Figure 5. Spectrum of vertical acceleration

下载: 全尺寸图片幻灯片

3.2 3类局部变形激励作用下振动舒适性评价

因局部变形高差与时间变化无关，可视为固定系数，因此，将代表性车速和变形长度代入式（24）、（26）和式（28）可得竖向加速度的时域函数，通过带通滤波器后可得到瞬时频率加权加速度，频率记权前后的加速度变化情况如图6所示.

根据运行均方根的定义，通过数值积分方法计算得到3类局部变形作用下的均方根加速度，如式（31）～（33）所示. 其中，加速度持续时间取为1 s^[11].

图 6 滤波前后的竖向加速度对比

Figure 6. Vertical acceleration comparison before and after filtering

下载: 全尺寸图片幻灯片

1）错台型

${a_{{\rm{w}}}}= {K_{\text{t}}}\Delta g{\left[ {{{\int_0^1 {\left( {\sum\limits_{i = 1}^4 {\frac{{{\lambda _i}{\mu _{i2}}}}{{{m_i}}}{{\rm{e}}^{{\lambda _i}t}}} } \right)} }^2}{\text{d}}t} \right]^{0.5}} = 42.04\Delta g \text{.}$

(31)

2）双断面型

$\begin{split}& {a_{{\rm{w}}}} = {K_{\text{t}}}\Delta g \;\times \hfill \\&\;\; {\left[ {\int_0^{\tfrac{l}{v}} {{{\left( {\sum\limits_{i = 1}^4 {\frac{{{\lambda _i}{\mu _{i2}}}}{{{m_i}}}{{\rm{e}}^{{\lambda _i}t}}} } \right)}^2}} {\rm{d}}t + \int_{\tfrac{l}{v}}^1 {{{\left( {\sum\limits_{i = 1}^4 {\frac{{{\lambda _i}{\mu _{i2}}}}{{{m_i}}}\left( {{{\rm{e}}^{{\lambda _i}t}} - {{\rm{e}}^{{\lambda _i}\left( {t - \tfrac{l}{v}} \right)}}} \right)} } \right)}^2} {\rm{d}}t} } \right]^{0.5}} \hfill \\&\;\; =72.07\Delta g\text{.} \\ \end{split}$

(32)

3）曲面型

${a_{{\rm{w}}}} = \frac{{4\Delta gv{K_{\text{t}}}}}{l}{\left( {A + B} \right)^{0.5}}= 61.28\Delta g \text{，}$

(33)

式中：

$\begin{gathered} A = \int_0^{\tfrac{l}{v}} {{{\left( {\sum\limits_{i = 1}^4 {\frac{{{\mu _{i2}}}}{{{\lambda _i}{m_i}}}\left( {\left( {{\lambda _i}l - 2v} \right){{\rm{e}}^{{\lambda _i}t}} + 2v} \right)} } \right)}^2}} {\rm{d}}t ；\hfill \\ B = \int_{\tfrac{l}{v}}^1 {{{\left( {\sum\limits_{i = 1}^4 {\frac{{{\mu _{i2}}}}{{{\lambda _i}{m_i}}}{{\rm{e}}^{{\lambda _i}\left( {t - l/v} \right)}}\left( {{\lambda _i}l{{\rm{e}}^{{\lambda _i}l/v}} + {\lambda _i}l - 2v{{\rm{e}}^{{\lambda _i}l/v}} + 2v} \right)\;} } \right)}^2} {\rm{d}}t}. \hfill \\ \end{gathered}$

变形高差与均方根加速度之间的理论关系如图7所示. 结合ISO 2631-1^[11]给出的舒适度评价标准，本文计算出的局部变形评价标准如表1所示. 当错台型、双断面型及曲面型局部变形分别达到 4.76、2.78、3.26 cm时，行驶质量均达到极不舒适.

综合上述分析，3类局部变形激励作用时的均方根加速度均和变形高差呈正比关系；在典型行驶车速40 km/h和变形长度0.6 m情况下，同样的变形高差作用，双断面型对人体舒适性影响最强烈，曲面型次之，错台型影响相对最小. 值得注意的是，实际路面中，曲面型变形多表现为路面拥包等损坏，长度不一. 当式（33）中局部变形长度l较大时，均方根加速度a_w的比例系数会减小，如图8所示. 本文所取的0.6 m变形长度几乎为曲面型变形的最不利长度，当长度约为0.3 ～ 1.5 m，同等高差时曲面型影响大于错台型，路面养护管理中应当重点考虑该长度范围.

图 7 均方根加速度与局部变形高差的关系

Figure 7. Relationship between root-mean-square acceleration and height difference of local deformation

下载: 全尺寸图片幻灯片

表 1 局部变形高差评价标准

Table 1. Evaluation criteria for height difference of local deformation

均方根加速度/ (m • s⁻²)	振动感受	△g 范围/cm
均方根加速度/ (m • s⁻²)	振动感受	错台型	双断面型	曲面型
< 0.315	没有不舒适	< 0.75	< 0.45	< 0.51
0.315～0.630	有些不舒适	0.75～ 1.50	0.45～ 0.87	0.51～ 1.03
0.500～1.000	比较不舒适	1.19～ 2.38	0.69～ 1.39	0.82～ 1.63
0.800～1.600	不舒适	1.90～ 3.81	1.11～ 2.22	1.31～ 2.61
1.250～2.500	很不舒适	2.97～ 5.95	1.73～ 3.47	2.04～ 4.08
> 2.000	极不舒适	> 4.76	> 2.78	> 3.26

下载: 导出CSV

| 显示表格

图 8 变形长度与比例系数的关系

Figure 8. Relationship between length of deformation and evaluation coefficient

下载: 全尺寸图片幻灯片

4. 结　论

1）对路面中常见的局部变形形式进行分析总结，建立了错台型、双断面型和曲面型3类典型局部变形的纵断面模型.

2）求解两自由度汽车强迫振动模型，推导出任意路面变形激励作用时人体竖向位移的时域变化函数；根据3类局部变形的纵断面特征从理论上得出人体竖向加速度的时域解，并对函数式中变形长度和行驶车速两个参数作出合理假定.

3）依据ISO 2631-1评价标准中竖向加速度的频率加权系数分布，设计了频率加权函数数字带通滤波器；对3类局部变形激励下的竖向加速度信号进行滤波，得到瞬时频率加权加速度并求解出均方根加速度；最后根据标准中舒适度划分方法确定了基于舒适性的3类局部变形评价标准. 结果表明：3类局部变形激励作用下，均方根加速度均与局部变形高差呈正比关系，同等变形程度下，双断面型变形的行驶不适度最高.

4）路面养护管理中，可结合高精度检测设备及时获取路面严重错台型及双断面型局部变形的位置，做好养护和修复工作，提高路面服务性能和寿命.

5）文中研究结果适用于城市道路沥青路面，城市道路水泥路面可以参考；对于公路路面中的类似问题，需要做进一步的研究.

图 1 车轮多边形磨耗照片

Figure 1. Polygonal wheel wear photo

下载: 全尺寸图片幻灯片

图 2 车轮失圆阶次表示

Figure 2. Wheel OOR described in the harmonic order levels

下载: 全尺寸图片幻灯片

图 3 A厂列车车轮多边形磨耗阶次

Figure 3. Harmonic order levels of wheel polygonal wear for manufacturer A

下载: 全尺寸图片幻灯片

图 4 A厂列车车轮典型多边形磨耗

Figure 4. Typical wheel polygonal wear for manufacturer A

下载: 全尺寸图片幻灯片

图 5 B厂列车车轮多边形磨耗阶次

Figure 5. Harmonic order levels of wheel polygonal wear for manufacturer B

下载: 全尺寸图片幻灯片

图 6 B厂车轮典型多边形磨耗

Figure 6. Typical wheel polygonal wear for manufacturer B

下载: 全尺寸图片幻灯片

图 7 高速动车组多边形磨耗特征

Figure 7. Wheel polygonal wear characteristic of high-speed EMU

下载: 全尺寸图片幻灯片

图 8 车轮径跳

Figure 8. Radial run-out of wheels

下载: 全尺寸图片幻灯片

图 9 不同线路条件下轨头垂向敲击垂向响应

Figure 9. Vertical frequency response function of railhead vertical excitation from impact hammer test with different railway lines

下载: 全尺寸图片幻灯片

图 10 多边形磨耗的极坐标表示

Figure 10. Polar-coordinate system of wheel OOR

下载: 全尺寸图片幻灯片

图 12 车轮踏面不同横向位置圆周磨耗

Figure 12. OOR wear at different transverse positions of wheel tread

下载: 全尺寸图片幻灯片

图 11 研磨子对车轮多边形磨耗影响

Figure 11. Influence of the abrasive block on roughness level of wheel OOR

下载: 全尺寸图片幻灯片

图 13 车轮踏面外形

Figure 13. Wheel tread profile

下载: 全尺寸图片幻灯片

表 1 高速动车组车轮不圆测试信息

Table 1. Measurement information of wheel OOR of high-speed EMUs

厂家	车型	设计速度/ （km·h⁻¹）	车轮初始直径/mm	运行线路	列次/列	测试车轮数/个
A	A1	200	860	兰州—乌鲁木齐	16	4122
	A2	250	860	大同—西安、贵阳—广州、南宁—广州、南昌—福州	130	4080
	A3	250	920	成都—贵阳	2	1608
	A4	300	860	大同—西安、北京—广州	53	6088
B	B1	200	890	兰州—乌鲁木齐	40	1800
	B2	250	860	大同—西安、上海—昆明	2	440
	B3	300	920	武汉—广州、北京—广州、北京—上海	502	8456
	B4	300	920	哈尔滨—大连	78	1902
	B5	350	920	大同—西安、郑州—徐州、北京—广州、北京—上海、北京—天津	12	2004

下载: 导出CSV

表 2 A厂列车车轮多边形磨耗波长

Table 2. Wavelengths of wheel polygonal wear for manufacturer A

车型	直径/mm	主要阶次/阶	波长/mm
A1	860	26～28	96～102
A2	860	14 22～24 29～30	192 112～122 90～92
A3	920	15 22～26	192 111～131
A4	860	10 17	250 150

下载: 导出CSV

表 3 B厂列车车轮多边形磨耗波长

Table 3. Wavelengths of wheel polygonal wear for manufacturer B

车型	直径/mm	主要阶次	波长/mm
B1	890	14～16 22～28	174～200 100～126
B2	860	14 22～24	192 112～122
B3	920	12 18～20	240 143～160
B4	920	21～23	124～136
B5	920	10 16～17	288 167～178

下载: 导出CSV

表 4 2种轴距下多边形磨耗对比

Table 4. Comparison of wheel polygonal wear under two wheelbases

轴距/m	车型	阶次	波长/mm	频率/Hz
2.5	A1～A4、 B2～B5	10～15	192、 240～288	330～354
		16～30	90～178	520～613
		21～23	124～136	613～672
2.7	B1	14～16	174～200	266～306
2.7	B1	21～28	100～126	423～533

下载: 导出CSV

表 5 高铁线路调查结果

Table 5. Investigation results for high-speed railway lines

线路	主要运行车型	轨道板类型	扣件类型	扣件间距/mm	扣件刚度/（kN·mm⁻¹）	频响主要频率/Hz
大西线	A2、A4、B2、B5	CRTS Ⅰ	WJ-8	629	20～40
贵广线	A2	CRTS Ⅰ	WJ-8	629	20～40	165
成贵线	A3	CRTS Ⅰ/CRTS Ⅲ	WJ-8	629	20～40	203
京广线	A4、B3、B5	CRTS Ⅰ	WJ-8/Vossloh 300	650	20～40	142
京沪线	B3、B5	CRTS Ⅱ	WJ-8/Vossloh 300	650	20～40	120
兰新线	A1、B1	CRTS Ⅰ	Vossloh	650	20～40	160
沪昆线	B2	CRTS Ⅰ	WJ-8	650	20～40	155
哈大线	B4	CRTS Ⅰ	WJ-7	629	30～40	252

下载: 导出CSV

参考文献(18)

[1]	陶功权,周小江,周业明,等. B型地铁车轮失圆问题分析[J]. 机械工程学报,2020,56(14): 152-160. doi: 10.3901/JME.2020.14.152 TAO Gongquan, ZHOU Xiaojiang, ZHOU Yeming, et al. Analysis of the wheel out-of-roundness of type B metro train[J]. Journal of Mechanical Engineering, 2020, 56(14): 152-160. doi: 10.3901/JME.2020.14.152
[2]	TAO G Q, WANG L F, WEN Z F, et al. Measurement and assessment of out-of-round electric locomotive wheels[J]. Proceedings of the Institution of Mechanical Engineers, Part F: Journal of Rail and Rapid Transit, 2018, 232(1): 275-287. doi: 10.1177/0954409716668210
[3]	金学松,吴越,梁树林,等. 车轮非圆化磨耗问题研究进展[J]. 西南交通大学学报,2018,53(1): 1-14. doi: 10.3969/j.issn.0258-2724.2018.01.001 JIN Xuesong, WU Yue, LIANG Shulin, et al. Mechanisms and countermeasures of out-of-roundness wear on railway vehicle wheels[J]. Journal of Southwest Jiaotong University, 2018, 53(1): 1-14. doi: 10.3969/j.issn.0258-2724.2018.01.001
[4]	TAO G Q, WEN Z F, JIN X S, et al. Polygonisation of railway wheels: a critical review[J]. Railway Engineering Science, 2020, 28(4): 317-345. doi: 10.1007/s40534-020-00222-x
[5]	刘佳,韩健,肖新标,等. 高速车轮非圆化磨耗对轴箱端盖异常振动影响初探[J]. 机械工程学报,2017,53(20): 98-105. doi: 10.3901/JME.2017.20.098 LIU Jia, HAN Jian, XIAO Xinbiao, et al. Influence of wheel non-circular wear on axle box cover abnormal vibration in high-speed train[J]. Journal of Mechanical Engineering, 2017, 53(20): 98-105. doi: 10.3901/JME.2017.20.098
[6]	ZHANG J, HAN G X, XIAO X B, et al. Influence of wheel polygonal wear on interior noise of high-speed trains[J]. Journal of Zhejiang University—Science A (Applied Physics & Engineering), 2014, 15(12): 1002-1018.
[7]	HU W G, LIU Z M, LIU D K, et al. Fatigue failure analysis of high speed train gearbox housings[J]. Engineering Failure Analysis, 2017, 73: 57-71. doi: 10.1016/j.engfailanal.2016.12.008
[8]	彭来先,韩健,初东博,等. 高速动车组垂向止挡异常振动特性及成因分析[J]. 机械工程学报,2019,55(12): 121-127. doi: 10.3901/JME.2019.12.121 PENG Laixian, HAN Jian, CHU Dongbo, et al. Analysis of abnormal vibration characteristics and causes of vertical block in high-speed EMU[J]. Journal of Mechanical Engineering, 2019, 55(12): 121-127. doi: 10.3901/JME.2019.12.121
[9]	肖俊恒,闫子权,涂英辉,等. 轮轨振动对高铁扣件伤损的影响分析[J]. 中国铁路,2017(11): 10-14. doi: 10.19549/j.issn.1001-683x.2017.11.010 XIAO Junheng, YAN Ziquan, TU Yinghui, et al. On impact of wheel-rail vibration on HSR fastener damage[J]. China Railway, 2017(11): 10-14. doi: 10.19549/j.issn.1001-683x.2017.11.010
[10]	PALLGEN G. Unrunde rader an eisenbahnfahraeugen[J]. Der Eisenbahningenieur, 1998, 49(1): 56-60.
[11]	JOHANSSON A. Out-of-round railway wheels—assessment of wheel tread irregularities in train traffic[J]. Journal of Sound and Vibration, 2006, 293(3/4/5): 795-806.
[12]	QU S, ZHU B, ZENG J, et al. Experimental investigation for wheel polygonisation of high-speed trains[J]. Vehicle System Dynamics, 2021, 59(10): 1573-1586. doi: 10.1080/00423114.2020.1772984
[13]	WU Y, DU X, ZHANG H J, et al. Experimental analysis of the mechanism of high-order polygonal wear of wheels of a high-speed train[J]. Journal of Zhejiang University—Science A (Applied Physics & Engineering), 2017, 18(8): 579-592. doi: 10.1631/jzus.A1600741
[14]	CAI W B, CHI M R, WU X W, et al. Experimental and numerical analysis of the polygonal wear of high-speed trains[J]. Wear, 2019, 440/441: 1-12.
[15]	金学松,吴越,梁树林,等. 高速列车车轮多边形磨耗、机理、影响和对策分析[J]. 机械工程学报,2020,56(16): 118-136. doi: 10.3901/JME.2020.16.118 JIN Xuesong, WU Yue, LIANG Shulin, et al. Characteristics, mechanism, influences and countermeasures of polygonal wear of high-speed train wheels[J]. Journal of Mechanical Engineering, 2020, 56(16): 118-136. doi: 10.3901/JME.2020.16.118
[16]	European Committee for Standardization. Railway applications-acoustics-rail and wheel roughness measurement related to noise generation: EN 15610[S]. Brussels: CEN-CENELEC Management Centre, 2019.
[17]	翟志浩,和振兴. 极寒环境下铁路扣件新型网孔式弹性垫板动力性能及影响研究[J]. 中国科学:技术科学,2020,50(2): 235-244. ZHAI Zhihao, HE Zhenxing. Research on dynamic performance and influence of new mesh-hole elastic pads for railway fasteners in an extremely cold environment[J]. Scientia Sinica (Technologica), 2020, 50(2): 235-244.
[18]	乔青峰,李明星,赵晓男,等. 研磨子抑制高速列车车轮多边形磨耗的机理研究[J]. 摩擦学学报,2020,40(2): 234-239. doi: 10.16078/j.tribology.2019037 QIAO Qingfeng, LI Mingxing, ZHAO Xiaonan, et al. Mechanism of suppression of polygonal wear of wheel on high-speed trains by abrasive block[J]. Tribology, 2020, 40(2): 234-239. doi: 10.16078/j.tribology.2019037

施引文献

期刊类型引用(3)

1.	蔡崇庆，姜涛，屈裕坤，姜学鹏. 新型深井立式车库通风设计优化. 暖通空调. 2025(02): 26-31 . 百度学术
2.	单勇，赵维豪. 基于通风与空气净化控制的动车车厢空气污染物的模拟研究. 轨道交通装备与技术. 2024(S1): 74-77+81 . 百度学术
3.	朱奋田，郑琦耀，许良中，昌超，牛纪强. 风环境下明线运行列车动力包通风和烟气扩散特性. 交通运输工程学报. 2024(05): 248-258 . 百度学术