Seismic Responses of RC Frame Structures with Special-Shaped
Columns and Girder-Transfer Story

ZHANGXin-pei; SHUTao; ZHAO Yong; HUHong-wei

Abstract

High rising buildings with a large bay supermarket on their two lowest stories and dwelling houses on the other stories were taken as the discussed object. Based on the relative national design codes and their engineering conditions, seismic responses of 72 high rising buildings of RC frame structures with special-shaped columns and a girder-transfer story were calculated respectively using the modal analysis method and the elastic time-history method, and in the calculation, different combinations of fortification intensities 7 and 8 and site soil typesⅠtoⅣwere considered. From the calculated result, the changes in their story drift, storyshear andmomenton structural bottomwere analyzed underthe condition of frequently occurred earthquakes, and some conclusions helpful for the design of high rising buildings were obtained.

Keywords:

随着我国高速铁路的迅猛发展，轮轨间动作用力不断增大，车轮多边形磨耗问题日益突出. 对2013年以前运营的数百个车轮进行了非圆化磨耗测试，统计结果显示，96%的车轮发生了偏心磨损，其次是低阶谐波磨损. 近些年来（主要指2014 年以后），高速列车车轮普遍出现了高阶多边形磨损，而车轮因多边形引起的轮轨冲击力达到车轮静载荷的1倍～3倍，轴箱加速度达到20g～40g^[1]. 金学松等^[1-3]通过长期的现场跟踪测试提出了我国最早发现的直线电机地铁车轮9边形磨耗的机理是轮对一阶弯曲振动所致，也正在通过测试给出多边形的一些机理解释. Morys^[4]将车轮和制动盘视为三维旋转弹簧-阻尼单元连接的刚体，构建了车辆-轨道动力学模型，结合轮对长期磨耗模型研究了车轮多边形的发展规律发现，多边形磨耗的发展过程很大程度上取决于轨道特征和激励频率. Meywerk^[5]建立单个轮对-轨道动力学模型，研究了轮对柔性和轨道弹性对车轮多边形磨耗的影响. Johansson等^[6]基于滚动接触理论，建立三维车辆-轨道动力学模型，研究认为车轮多边形磨耗与车辆/轨道耦合系统的垂直共振有关. Wu等^[7-8]研究认为同一转向架两轮对之间钢轨的3阶弯曲振动是导致车轮高阶多边形磨耗的主要因素，然而在运行中能否被激发还有待验证. Meinke等^[9]提出高速轮对不同于普速轮对，其动力学是由陀螺力矩和惯性力矩决定的，研究了静、动不平衡对车轮多边形磨耗的影响. 胡晓依等^[10]将车辆/轨道耦合动力学模型与Archard磨耗模型相结合，建立长期磨损迭代模型模拟运行时车轮多边形的整个发展历程. 宋志坤等^[11]通过对比分析车轮粗糙度对车轮多边形的影响发现，提高镟修质量可以减缓多边形的增长速度.

上述研究主要是从现场跟踪或是车辆-轨道大系统动力学仿真入手，获得了一些可以解释车轮多边形现象的原因. 对于高速动车组车轮高阶多边形磨耗的成因，业界提出其受多重因素影响，很难给出统一的解释. 而现场目前不得不采用的方式是镟轮，这不仅缩短了车轮的使用寿命，也不符合低碳减排的国家要求. 所以研究车轮多边形磨耗形成和发展机理，明确其主要影响因素，是当前车轮多边形磨耗问题研究的重点. 为此，本文基于转子动力学理论建立轮轨动力学模型，结合磨耗模型研究多边形磨耗的发生和演化.

1. 车轮多边形磨耗发生与演化模型

1.1 轮对-轨道转子系统动力学模型

研究表明，轮对蠕滑是车轮多边形磨耗发生和发展的主要原因^[12]. 轨道车辆结构庞大，自由度较多，其非线性因素复杂，通常根据研究需要对重点关注点精细建模，对于其他部分做适当简化. 因此，本文选取自由度较少的轮对轴箱装置为研究对象，如图1所示，将其简化为轮对-轨道转子系统，其中包括轴承、车轮、车轴及模拟轨道的两个大直径圆盘. 图1中：m₁、m₄分别为左、右轴箱质量；m₂、m₃分别为左、右车轮的质量；e₁、e₂分别为左右车轮的偏心量；k_bx1、k_bx2（k_bz1、k_bz2）分别为左、右轴承的纵向（垂向）刚度；c_bx1、c_bx2（c_bz1、c_bz2）分别为左、右轴承的纵向（垂向）的阻尼；k₁、c₁、L₁分别为左、右轴段的刚度、阻尼和长度；k₂、c₂、L₂分别为中间轴段的刚度、阻尼和长度； k_R、c_R、f_st分别为扣件的刚度、阻尼和预压缩量； ω 为圆盘绕 O₁自转的角速度.

图 1 轮对转子系统动力学模型

Figure 1. Dynamic model of wheelset rotor system

下载: 全尺寸图片幻灯片

为了分析转子的基本特征，假设：

1）车轮为刚性圆盘；

2）车轴为等直圆轴，具有一定的弯曲刚度和无限大的扭转刚度，两端相同的轴箱弹性支承；

3）将簧上质量等效为轨道预压力，忽略轮轨接触斑处相对变形，假设轮轨一体支撑在扣件上.

取O-xyz为固定坐标系，圆盘所在平面与弹性轴两端支承点连线的交点O为固定坐标系原点，y轴沿转子轴线，圆盘所在平面为xOz坐标参考平面. 圆盘瞬时位置和受力情况见图2，图中：O₁为圆盘形心，C为圆盘质心，两者之间的距离O₁C，即偏心量e；r为动挠度；F为轴的弹性恢复力；R为黏性阻尼力；W为重力；F_x和F_z分别为轮轨接触纵向蠕滑力和法向力. 选取（x(t)，z(t)）为圆盘形心O₁广义坐标；t为时间.

图 2 圆盘瞬时位置及受力

Figure 2. Disk instantaneous position and force

下载: 全尺寸图片幻灯片

1.2 轮轨滚动接触力模型

围绕轮轨滚动接触问题的求解，研究发展了许多理论模型，较著名的理论有Carter 理论、Kalker理论、沈氏理论等，其中Kalker 简化理论在铁路领域应用最广泛^[13]. 蠕滑力除了与蠕滑率有关以外还与接触区的形状、表面润滑状态、正压力的分布等有关. Kalker公式^[14-15]为

$\left\{\begin{array}{l}{F}_{x}=-Gab{C}_{11}{\gamma }_{x},\\ {F}_{y}=-Gab{C}_{22}{\gamma }_{y}-G(ab{)}^{3\text{/}2}{C}_{23}{\xi }_{{\textit{z}}},\end{array} \right.$

(1)

式中：F_y为横向蠕滑力；a、b分别为接触椭圆的长、短半轴长；G为剪切模量；C₁₁、C₂₂ 、C₂₃为蠕滑系数；γ_x、 γ_y、 ξ_z分别为纵向、横向和自旋蠕滑率.

根据接触椭圆偏心率及第一、二类完整椭圆积分可确定出接触斑长、短半轴长度^[13]分别为

$\left\{ \begin{gathered} a = {m_{\rm{c}}}{\left(\frac{{3{\text{π}} ({\delta _1} + {\delta _2})}}{{2\beta }}{F_{\textit{z}}}\right)^{1/3}} , \\ b = {n_{\rm{c}}}{\left(\frac{{3{\text{π}} ({\delta _1} + {\delta _2})}}{{2\beta }}{F_{\textit{z}}}\right)^{1/3}}, \\ \end{gathered} \right.$

(2)

$\left\{ \begin{gathered} {m_{\rm{c}}} = {\left(\frac{{2{E_{\rm{p}}}}}{{{\text{π}} (1 - e_{\rm{c}}^2)}}\right)^{1/3}}, \\ {n_{\rm{c}}} = {\left(\frac{{2{E_{\rm{p}}}\sqrt {1 - e_{\rm{c}}^2} }}{{\text{π}} }\right)^{1/3}} , \\ {\delta _i} = \frac{{1 - {\upsilon _i}^2}}{{{\text{π}} {E_i}}}, \quad i=1, 2, \\ \beta = 1{\text{/}}{R_{11}} + 1{\text{/}}{R_{12}} + 1{\text{/}}{R_{21}} + 1{\text{/}}{R_{22}}, \end{gathered} \right.$

(3)

式中：E_p为第二类完整椭圆积分；e_c为接触椭圆偏心率；υ₁、υ₂分别为接触体1、2的泊松比；E₁、E₂分别为接触体1、2的弹性模量；R₁₁、R₁₂和R₂₁、R₂₂分别为车轮和钢轨沿接触坐标系纵、横向的主曲率半径.

合成的蠕滑力F_rr为

$F_{{\rm{rr}}} = {(F_x^2 + F_y^2)^{\tfrac{1}{2}}} .$

(4)

修正蠕滑力F_r为

${F_{\text{r}}} = \left\{ \begin{array}{*{20}{l}} \mu {F_{\textit{z}}}\left[\left(\dfrac{F_{\rm{rr}}}{{\mu {F_{\textit{z}}}}}\right) - \dfrac{1}{3}{{\left(\dfrac{{F_{\rm{rr}}}}{{\mu {F_{\textit{z}}}}}\right)}^2} + \dfrac{1}{{27}}{{\left(\dfrac{{F_{\rm{rr}}}}{{\mu {F_{\textit{z}}}}}\right)}^3}\right],\;\\ \quad F_{\rm{rr}} \leqslant 3\mu {F_{\textit{z}}}, \\ \mu {F_{\textit{z}}},\quad F_{\rm{rr}} \gt 3\mu {F_{\textit{z}}}, \end{array} \right.$

(5)

式中：μ为摩擦系数；${F_{\textit{z}}} = - {k_{\rm{R}}}({f_{{\rm{st}}}} + {\textit{z}})$.

1.3 轮对转子系统的动力学方程

参照多跨不平衡轴系的非线性动力学建模方法^[16]，结合拉格朗日方程建立轮对-轨道转子系统的动力学方程，如式（6）.

$\left\{ \begin{array}{l} {m_1}{{\ddot x}_1} + {c_1}{{\dot x}_1} + {k_1}{x_1} - {k_1}{x_2} = {F_{{\rm{b}}x1}}, \\ {m_1}{{\ddot {\textit{z}}}_1} + {c_1}{{\dot {\textit{z}}}_1} + {k_1}{{\textit{z}}_1} - {k_1}{{\textit{z}}_2} = {F_{{\rm{b}}{\textit{z}}1}} - {m_1}g, \\ {m_2}{{\ddot x}_2} + {c_2}{{\dot x}_2} - {k_1}{x_1} + ({k_1} + {k_2}){x_2} - {k_2}{x_3} = \\ \quad {m_2}{e_1}{\omega ^2}\cos \;\omega t + {F_{x1}}, \\ {m_2}{{\ddot {\textit{z}}}_2} + {c_2}{{\dot {\textit{z}}}_2} - {k_1}{{\textit{z}}_1} + ({k_1} + {k_2}){{\textit{z}}_2} - {k_2}{{\textit{z}}_3} =\\ \quad {m_2}{e_1}{\omega ^2}\sin \;\omega t - {m_2}g + {F_{{\textit{z}}1}}, \\ {m_3}{{\ddot x}_3} + {c_3}{{\dot x}_3} - {k_2}{x_2} + ({k_1} + {k_2}){x_3} - {k_1}{x_4} =\\ \quad {m_3}{e_2}{\omega ^2}\cos\; \omega t + {F_{x2}}, \\ {m_3}{{\ddot {\textit{z}}}_3} + {c_3}{{\dot {\textit{z}}}_3} - {k_2}{{\textit{z}}_2} + ({k_1} + {k_2}){{\textit{z}}_3} - {k_1}{{\textit{z}}_4} =\\ \quad {m_4}{e_2}{\omega ^2}\sin \;\omega t - {m_3}g + {F_{{\textit{z}}2}}, \\ {m_4}{{\ddot x}_4} + {c_4}{{\dot x}_4} + {k_1}{x_4} - {k_1}{x_3} = {F_{{\rm{b}}x2}}, \\ {m_4}{{\ddot {\textit{z}}}_4} + {c_4}{{\dot {\textit{z}}}_4} + {k_1}{{\textit{z}}_4} - {k_1}{{\textit{z}}_3} = {F_{{\rm{b}}{\textit{z}}2}} - {m_4}g . \end{array} \right.$

(6)

表示成矩阵形式为

${\boldsymbol{M\ddot u}} + {\boldsymbol{C\dot u}} + {\boldsymbol{Ku}} = {{\boldsymbol{F}}_{\text{b}}} + {{\boldsymbol{F}}_{\text{e}}} + {{\boldsymbol{F}}_{\text{r}}} + {{\boldsymbol{F}}_{\text{g}}},$

(7)

式中：M、C、K、F_b、F_e、F_r和F_g分别为转子系统的质量矩阵、阻尼矩阵、刚度矩阵、轴承力向量、不平衡力向量、轮轨接触力向量和重力向量；$ \ddot{{\boldsymbol{u}}}、\dot{{\boldsymbol{u}}}、{\boldsymbol{u}} $分别为模型中质点的加速度向量、速度向量和位移向量.

车轴各轴段的刚度如式（8）.

${k_i} = \frac{{12{E_{{\rm{a}}i}}{I_i}}}{{L_i^3}},$

(8)

式中：E_ai为车轴材料的弹性模量；I_i为车轴截面惯性矩，如式（9），d_aoi和d_aii为第i轴段的车轴外径和内径；L_i为轴段长度.

${I_i} = \frac{{\text{π}} }{{64}}(d_{{\rm{ao}}i}^4 - d_{{\rm{ai}}i}^4).$

(9)

1.4 车轮圆周磨耗深度模型

在车轮多边形磨耗预测中，通常采用的材料磨耗模型有Archard 模型和磨耗功模型. 研究表明，两者对应的车轮多边形演变趋势基本相同，但后者所得结果与现场实测数据更为接近，且计算效率较高^[17]，因此本文采用磨耗功模型，如式（10）.

$\Delta m = {K_{\rm{w}}} {W_{\rm{w}}} ,$

(10)

式中：Δm为磨耗质量；K_w为磨耗功系数；W_w为接触斑内的磨耗功.

假设车轮在磨耗过程中钢轨平顺，横向磨耗均匀，仅考虑周向失圆. 结合式（10）推导可得，车轮圆周磨耗深度Δr为

$\begin{split} & \Delta r = R(\theta - 2{\text{π}} ) - R(\theta ) = \frac{{\Delta m}}{{\rho {A_{\rm{w}}}}} = \frac{{{K_{\rm{w}}}{W_{\rm{w}}}}}{{\rho {A_{\rm{w}}}}} = \\ &\quad {K_{\rm{w}}}\frac{{{F_{{x}}} {{{s}}_{\rm{w}}}}}{{\rho {b_{\rm{w}}}{{{s}}_{\rm{w}}}}} = {K_{\rm{w}}}\frac{{{F_x}}}{{\rho {b_{\rm{w}}}}}, \end{split}$

(11)

式中：R(θ−2π)为车轮滚动圆前一次磨耗后的半径，θ为车轮滚动圆的圆周某处对应的角度；R(θ)为后一次磨耗后的半径；A_w为磨耗面积；ρ为材料密度；s_w为轮轨相对滑动距离；b_w为磨耗区的平均宽度.

由式（11）可见，车轮磨耗与纵向蠕滑力有关，当纵向蠕滑力按照一定的频率周期性变化时，导致磨耗深度周期性变化，多边形磨耗纵向蠕滑力最大时，将出现磨耗峰值，且在车轮纵向振动的一个周期内只出现一次. 设车轮纵向振动频率为f_x ,则相邻两个磨耗峰值的时间间隔为1/ f_x ,轮对转频为 f₂，则车轮转动周期为1/ f₂. 可得车轮转动一圈内出现的磨耗峰值数（阶次）n为

$n = \frac{{{f_x}}}{{{f_{\text{2}}}}}.$

(12)

轮对转频 f₂与列车运行速度v的关系为$ {f_{\text{2}}} = v/({\text{π}} D) $，其中，D为车轮名义直径. 代入式（12）可知，车轮纵向振动频率为

${f_x} = \frac{{nv}}{{{\text{π}} D}}.$

(13)

2. 多边形磨耗仿真分析

选取某型高速动车组轮对参数：车轴内径为30 mm，车轴外径为190 mm，轴颈中心距为1956 mm，滚动圆间距为1493 mm，寿命周期内车轮直径D_w=920～830 mm，ρ=7850 kg/m³，E₁=E₂=2.06 × 10¹¹ Pa，v₁=v₂=0.3，e=0.1 mm，k_R=20 MN/m，（轴承刚度）k_b= 1×10¹⁰ N/m，μ = 0.25，v=300 km/h.

2.1 车轮多边形磨耗特征

分析e=0.1 mm、v=300 km/h、D_w=920 mm时，左侧车轮的纵向振动位移时域图和频域图、轴心轨迹及车轮圆周磨耗图以表征车轮的振动特征和磨耗特征，如图3.

图 3 轮对系统的振动和磨耗特征

Figure 3. Vibration characteristics of wheelset rotor system

下载: 全尺寸图片幻灯片

由图3（a）、（b）可见：纵向振动位移为含有两种频率成分的周期信号，其低频频率等于轮对转频，高频频率接近于600 Hz，频谱图可以有效表征轮对-轨道转子系统的振动特征. 由图3（c）、（d）可见：轴心轨迹呈波浪状椭圆，车轮圆周磨耗形状为偏心多阶多边形；车轮多边形磨耗与轮对纵向周期振动有关，车轮磨耗形状图可有效表征车轮多边形的磨耗特征.

2.2 偏心量的影响

偏心量e在0.01～0.10 mm时对轮对-轨道转子系统振动特征的影响如图4. 图4（a）可见，不同偏心量下，振动位移包含两种频率成分：低频频率等于轮对转频，其幅值A随着偏心量的增大而增大，说明偏心量引起的惯性力增加，使得轮对振动幅值增大；高频频率均固定在580 Hz附近，说明高频频率与偏心量无关.

图 4 不同偏心量的幅频

Figure 4. Amplitude-frequency diagram under different eccentricities

下载: 全尺寸图片幻灯片

e=0.01，0.10，0.50，1.00 mm时，偏心量对车轮磨耗特征的影响如图5. 由图5可见：随着偏心量的增大，车轮出现了明显的椭圆化和高阶多边形磨耗现象，高阶多边形磨耗幅值略有减小，但阶次恒定. 此时的车轮磨耗为一阶叠加了多阶的多边形磨耗，且多阶多边形的阶次与偏心量的变化无关.

图 5 不同偏心量的车轮磨耗形状

Figure 5. Wheel wear shapes under different eccentricity

下载: 全尺寸图片幻灯片

2.3 运行速度和车轮直径的影响

已有研究发现，车轮多边形磨耗的发生与演化取决于运行速度、车轮直径和轮轨固有振动^[18]. 高铁为300 km/h恒速运行方式，由于目前动车组采用镟轮检修模式，车轮寿命期内的车轮直径在920～830 mm变化，对应取轮对角速度为180.8～205.9 rad/s，相当于轮对转频为25.2～58.0 Hz.

图6为不同角速度的频谱图，由图6可见，不同轮对角速度下，轮对振动位移也包含两种频率成分：低频频率随着转速的增大而增大，大小等于该转速对应的转频；高频频率仍然为一个固定频率，大小约为580 Hz.

图 6 不同角速度ω的频谱

Figure 6. Spectrum diagram under different angular velocities

下载: 全尺寸图片幻灯片

为了进一步研究轮对转频对车轮多边形磨耗特征的影响，定义频率比N为固定频率f₁和轮对转频f₂之比. 取频率比N=18.0～20.5，绘制车轮圆周磨耗图，如图7. 由图7可见：当固定频率与轮对转频比为整数倍时，即N=18.0，19.0，20.0时，车轮圆周方向上发生并演化出了多边形，多边形阶数等于频率比，且车轮转速越高，多边形阶数越小；当固定频率与轮对转频比不为整数，即N=18.5，19.5，20.5时，未演化成带有明显阶次的车轮多边形，车轮磨耗为均匀磨耗. 可见，车轮多边形磨耗的演化呈现明显的“频率整分”特征.

图 7 不同频率比下的车轮磨耗形状

Figure 7. Wheel wear shapes under different frequency ratios

下载: 全尺寸图片幻灯片

车轮纵向振动存在一个固定频率，当该频率接近于车轮转频的整数倍时，车轮磨耗快速演化为阶数等于该整数倍的多边形.

2.4 现场跟踪实测验证

现场实测发现，列车运行速度为 300 km/h，在车轮轮径寿命周期从920～830 mm 变化过程中，存在3个车轮多边形磨耗高速发展时期，车轮直径为830、875 mm和915 mm，对应18阶、19阶和20阶车轮多边形磨耗，而在其他直径未发生明显阶次的多边形磨耗^[19].

将现场实测数据代入式（12）和式（13），可计算出对应的轮对转频f₂和纵向振动频率f_x，见表1.

表 1 车轮多边形磨耗参数分析

Table 1. Analysis of wheel polygonal wear parameters

n	D/mm	ω/（rad·s⁻¹）	f₂/Hz	f_x/Hz
20	915	182.1	29.0	580.1
19	875	190.5	30.3	576.3
18	830	200.8	32.0	575.6

下载: 导出CSV

| 显示表格

由表1可见：v=300 km/h，车轮直径分别为830、875 mm和915 mm的3种情况均对应一个大小约为580 Hz的固定频率，该“固定频率”分别等于对应的轮对转频的整数倍，同时也分别发生了与相应整数倍相等的18阶、19阶和20阶车轮多边形磨耗；在车轮轮径寿命周期内，对应的轮对转频不能整除固定频率的车轮直径时，却未检测到明显的多边形磨耗. 综上，现场实测数据分析表明，车轮边形磨耗呈“频率固定和转频整分”特征，此磨耗特征与本文建立的轮对-轨道转子动力学系统车轮多边形演化模型所得结论一致.

3. 固定频率的来源与影响

为了探究上述研究中的固定频率来源，利用本文建立的轮对转子动力学模型进行模态分析，并与建立的轮轨接触实体有限元模型模态分析结果进行对比验证. 进而研究车轴直径等系统参数对固定频率的影响，并确定各参数的灵敏度.

3.1 集中质量模型模态分析

令系统主振动$ u = \varphi {\text{sin}}(\omega t + \phi ) $，代入式（7）对应的系统${\boldsymbol{M}}\ddot {\boldsymbol{u}} + {\boldsymbol{Ku}} = {\boldsymbol{0}},\;{{\boldsymbol{u}} \in {{\boldsymbol{R}}_{8 \times 1}}}$，可得系统的特征方程为

$({\boldsymbol{K}} - {\omega ^2}{\boldsymbol{M}}){\boldsymbol{\phi}} = {\boldsymbol{0}},$

(14)

式中：${\boldsymbol{ \phi}} = {({\phi _1},{\phi _2},\cdots,{\phi _8})^{\text{T}}} $，为系统振动幅值列向量.

由式（14）得

$\left\{ {\begin{array}{*{20}{l}} {{\boldsymbol{\phi}} _{i}^{\rm{T}}}{\boldsymbol{M}}{{\boldsymbol{\phi}} _{j}} = 0,\; {{\boldsymbol{\phi}} _{i}^{\text{T}}}{\boldsymbol{K}}{{\boldsymbol{\phi}} _{j}} = 0,\quad{i \ne j} , \\ {{\boldsymbol{\phi}} _{i}^{\text{T}}}{\boldsymbol{M}}{{\boldsymbol{\phi}} _{j}} = {{m_{{\rm{p}}i}}},\; {{\boldsymbol{\phi}} _{i}^{\text{T}}}{\boldsymbol{K}}{{\boldsymbol{\phi}} _{j}} = {{k_{{\rm{p}}i}}},\quad{i = j} , \end{array}} \right.$

式中：${\boldsymbol{\phi}} _{i}$、${\boldsymbol{\phi}} _{j} $为系统第i、j阶模态的振幅向量；m_pi、k_pi为当i=j时系统的模态质量和模态刚度.

利用Kronecker符号：

$\left\{ {\begin{array}{*{20}{l}} {{{\boldsymbol{\phi}} _{i}^{\rm{T}}}{{M}}{{\boldsymbol{\phi}} _{j}} = {\delta _{{{ij}}}}{m_{{\rm{p}}i}}}, \\ {{{\boldsymbol{\phi}} _{i}^{\rm{T}}}{{K}}{{\boldsymbol{\phi}} _{j}} = {\delta _{{{ij}}}}{k_{{\rm{p}}i}}} , \end{array}} \right.$

(15)

式中：i=j时${\delta _{{{ij}}}}=1 $，${i \ne j} $时$ {\delta _{{{ij}}}}=0$.

由$ {{\boldsymbol{\phi }}_{i}^{\text{T}}}{{K}}{{\boldsymbol{\phi}} _{i}} = {\omega _i}^2{{\boldsymbol{\phi}} _{i}^{\text{T}}}{{M}}{{\boldsymbol{\phi}} _{i}} \Rightarrow {k_{{\rm{p}}i}} = {\omega _i}^2{m_{{\rm{p}}i}} $，得第i阶固有频率为

${\omega _i} = \sqrt {\frac{{{k_{{\rm{p}}i}}}}{{{m_{{\rm{p}}i}}}}} .$

(16)

由式（7）可得

$\left\{ \begin{gathered} {m_1}{{\ddot x}_1} + {c_1}{{\dot x}_1} + ({k_1} - {k_{\rm{b}}}){x_1} - {k_1}{x_2} = 0 , \\ {m_1}{{\ddot {\textit{z}}}_1} + {c_1}{{\dot {\textit{z}}}_1} + ({k_1} - {k_{\rm{b}}}){{\textit{z}}_1} - {k_1}{{\textit{z}}_2} = 0 , \\ {m_2}{{\ddot x}_2} + {c_2}{{\dot x}_2} - {k_1}{x_1} + ({k_1} + {k_2}){x_2} - {k_2}{x_3} + \mu {k_{\rm{R}}}{{\textit{z}}_2} = 0, \\ {m_2}{{\ddot {\textit{z}}}_2} + {c_2}{{\dot {\textit{z}}}_2} - {k_1}{{\textit{z}}_1} + ({k_1} + {k_2} + {k_{\rm{R}}}){{\textit{z}}_2} - {k_2}{{\textit{z}}_3} = 0 , \\ {m_3}{{\ddot x}_3} + {c_3}{{\dot x}_3} - {k_2}{x_2} + ({k_1} + {k_2}){x_3} - {k_1}{x_4} + \mu {k_{\rm{R}}}{{\textit{z}}_3} = 0, \\ {m_3}{{\ddot {\textit{z}}}_3} + {c_3}{{\dot {\textit{z}}}_3} - {k_2}{{\textit{z}}_2} + ({k_1} + {k_2} + {k_{\rm{R}}}){{\textit{z}}_3} - {k_1}{{\textit{z}}_4} = 0 , \\ {m_4}{{\ddot x}_4} + {c_4}{{\dot x}_4} + ({k_1} - {k_{\rm{b}}}){x_4} - {k_1}{x_3} = 0 , \\ {m_4}{{\ddot {\textit{z}}}_4} + {c_4}{{\dot {\textit{z}}}_4} + ({k_1} - {k_{\rm{b}}}){{\textit{z}}_4} - {k_1}{{\textit{z}}_3} = 0 . \\ \end{gathered} \right.$

(17)

式（17）结合特征方程式（14）可得：

$\begin{gathered} \left( \begin{gathered} \left[ {\begin{array}{*{20}{c}} {{k_1} - {k_{\rm{b}}}} & 0 & { - {k_1}} & 0 & 0 & 0 & 0 & 0 \\ 0 & {{k_1} - {k_{\rm{b}}}} & 0 & { - {k_1}} & 0 & 0 & 0 & 0 \\ { - {k_1}} & 0 & {{k_1} + {k_2}} & {\mu {k_{\rm{R}}}} & { - {k_2}} & 0 & 0 & 0 \\ 0 & { - {k_1}} & 0 & {{k_1} + {k_2} + {k_{\rm{R}}}} & 0 & { - {k_2}} & 0 & 0 \\ 0 & 0 & { - {k_2}} & 0 & {{k_1} + {k_2}} & {\mu {k_{\rm{R}}}} & { - {k_1}} & 0 \\ 0 & 0 & 0 & { - {k_2}} & 0 & {{k_1} + {k_2} + {k_{\rm{R}}}} & 0 & { - {k_1}} \\ 0 & 0 & 0 & 0 & { - {k_1}} & 0 & {{k_1} - {k_{\rm{b}}}} & 0 \\ 0 & 0 & 0 & 0 & 0 & { - {k_1}} & 0 & {{k_1} - {k_{\rm{b}}}} \end{array}} \right] - {\omega ^2}\left[ {\begin{array}{*{20}{c}} {{m_1}} & {} & {} & {} & {} & {} & {} & {} \\ 0 & {{m_1}} & {} & {} & {} & {} & {} & {} \\ 0 & 0 & {{m_2}} & {} & {} & {} & {} & {} \\ 0 & 0 & 0 & {{m_2}} & {} & {} & {} & {} \\ 0 & 0 & 0 & 0 & {{m_3}} & {} & {} & {} \\ 0 & 0 & 0 & 0 & 0 & {{m_3}} & {} & {} \\ 0 & 0 & 0 & 0 & 0 & 0 & {{m_4}} & {} \\ 0 & 0 & 0 & 0 & 0 & 0 & 0 & {{m_4}} \end{array}} \right] \\ \end{gathered} \right)\left[ {\begin{array}{*{20}{c}} {{\phi _1}} \\ {{\phi _2}} \\ {{\phi _3}} \\ {{\phi _4}} \\ {{\phi _5}} \\ {{\phi _6}} \\ {{\phi _7}} \\ {{\phi _8}} \end{array}} \right] = {\boldsymbol{0}} \\ \end{gathered}.$

(18)

将系统质量参数m₁、m₂、m₃、m₄和刚度参数k₁、k₂、k_b、k_R代入求解，获得了式（18）中的8阶模态，其中固定频率附近的各阶模态见图8. 由图8可见：固定频率等于2.1节中发现的580 Hz附近时，轮对转子系统的振动为二阶弯曲振动.

图 8 不同固有频率下的振型

Figure 8. Mode shapes at different natural frequencies

下载: 全尺寸图片幻灯片

3.2 实体有限元模型模态分析

为了进一步研究固定频率与钢轨、轮对及耦合的关系，除上述轮对的尺寸参数外，取60 kg钢轨断面，钢轨长度6.0 m，轨枕距625 mm，扣件宽度170 mm，扣件垂向刚度20 MN/m，转向架固定轴距2.5 m，建立轮轨耦合实体有限元模型进行模态分析，研究转动效应对轮轨耦合系统固有频率的影响.

由表2可见：不考虑转动效应的刚轮柔轨、柔轮刚轨、柔轮柔轨，分别对应有1、4、5个固有频率接近“固定频率”580 Hz，而不同转速下轮轨耦合系统的固有频率均出现了5个接近固定频率的值. 说明转速和钢轨柔性对固定频率附近的固有频率影响不大，可以认为轮对对固定频率的贡献比钢轨大，虽然钢轨柔性对固定频率的贡献不大，但是将轮对和钢轨都考虑成柔性更符合实际情况. 所以在研究不同转速对轮轨耦合系统的固有频率影响时，将轮轨考虑成柔性.

表 2 轮轨耦合系统的固有频率

Table 2. Natural frequency of rotor system of wheel set Hz

不考虑转动效应			考虑转动效应
刚轮柔轨	柔轮刚轨	柔轮柔轨	柔轮柔轨
519.40	567.16	565.42	565.42
559.49	567.19	580.73	580.61
574.77	577.74	583.07	582.97
593.11	578.03	584.48	584.37
605.90	585.43	585.47	585.65
608.93	585.45	586.78	586.94
622.27	748.22	591.33	591.35
648.97	769.92	621.80	621.81

下载: 导出CSV

| 显示表格

考虑转动效应下柔性轮对柔性钢轨的振动模态如图9所示，由图9可知：固定频率附近的柔性轮轨振动均对应轮对的2阶弯曲，而钢轨则分别呈现8阶横向弯曲、8阶横向弯扭耦合和8阶垂向弯曲；轮对的2阶弯曲振动是固定频率的主要来源，这一结论验证了式（18）中的集中质量模型的分析结果，也再次验证了本文所建立的轮对轨道转子动力学模型的可行性. 综上所述，出现580 Hz左右的“频率固定”现象的原因，主要由轮对的2阶弯曲所致，钢轨柔性对固定频率的影响不大.

图 9 考虑转动效应下轮轨耦合系统振动模态

Figure 9. Vibration modes of wheel-rail coupling system considering rotation effect

下载: 全尺寸图片幻灯片

3.3 扣件刚度对固定频率的影响

由式（8）、（18）可见，车轮多边形磨耗的固定频率主要与车轴直径、轮对质量、轴承刚度和扣件刚度有关. 分别研究扣件刚度为20、50、100、150 MN/m时，扣件刚度对固定频率的影响.

图10为同扣件刚度幅频图，由图10可见：随着扣件刚度的增加，固定频率及其对应的振动幅值均略有增大. 图11为不同扣件刚度的磨耗特征，由图11可见：随着扣件刚度的增大，由于固定频率的增大，与转频的比值逐渐远离整数倍，车轮磨耗形状由k_R=20 MN/m的明显多阶多边形逐渐变为k_R=150 MN/m均匀磨耗. 可见，增大扣件刚度对车轮多边形磨耗有一定的抑制作用.

图 10 不同扣件刚度幅频图

Figure 10. Amplitude-frequency diagram under different fastener stiffnesses

下载: 全尺寸图片幻灯片

图 11 不同扣件刚度的磨耗特征

Figure 11. Wear characteristics under different fastener stiffnesses

下载: 全尺寸图片幻灯片

3.4 车轴外径对固定频率的影响

转速一定，取车轴外径分别为171、190 mm和209 mm，研究车轴外径对固定频率的影响.

图12为不同车轴直径（d_ao）的幅频图，由图12可见：随着车轴外径的增加，由于轮对刚度增加，固定频率增大，转频和固定频率的幅值均减小. 图13为不同车轴直径的磨耗特征，由图13可见：增大和减小车轴直径，均能因固定频率不再等于转频的整数倍而避免多阶多边形的发展；减小车轴外径会加速车轮的均匀磨耗. 可见，通过增大车轴外径可有效抑制车轮多边形.

图 12 不同车轴直径的幅频图

Figure 12. Amplitude-frequency diagram under different diameters

下载: 全尺寸图片幻灯片

图 13 不同车轴直径的磨耗特征

Figure 13. Wear characteristics under different diameters

下载: 全尺寸图片幻灯片

3.5 车轮质量对固定频率的影响

取车轮质量（m）分别为360、400 kg和440 kg，研究车轮质量对固定频率的影响. 由图14的振动特征可见：随着车轮质量的减小，固定频率增大，转频和固定频率的幅值均减小. 由图15的磨耗特征可见：增大和减小车轮质量，均能因固定频率不再等于转频的整数倍而避免多阶多边形的发展；增大车轮质量会加速车轮的均匀磨耗量. 因此，可以通过车轮轻量化抑制车轮多边形磨耗.

图 14 不同车轮质量的幅频图

Figure 14. Amplitude-frequency diagram under different wheel masses

下载: 全尺寸图片幻灯片

图 15 不同车轮质量的磨耗特征

Figure 15. Wear characteristics under different wheel masses

下载: 全尺寸图片幻灯片

3.6 支承刚度对固定频率的影响

取支承刚度分别为5000、10000、15000 MN/m，研究支承刚度对固定频率的影响.

由图16中的振动特征可见：随着支承刚度的减小，固定频率减小，转频和固定频率的幅值均增大. 由图17中的磨耗特征可见：增大和减小支承刚度，均能因固定频率不再等于转频的整数倍而避免多阶多边形的发展；减小支承刚度会加速车轮的均匀磨耗量. 因此，可以通过增大支承刚度抑制车轮多边形磨耗.

图 16 不同支承刚度的幅频图

Figure 16. Amplitude-frequency diagram under different support stiffnesses

下载: 全尺寸图片幻灯片

图 17 不同支承刚度的磨耗特征

Figure 17. Wear characteristics under different support stiffnesses

下载: 全尺寸图片幻灯片

3.7 灵敏度分析

改变车轴外径等系统参数可一定程度抑制车轮多边形磨耗的演化，为了给现场抑制车轮多边形提供理论依据，分别研究各系统参数对固定频率的灵敏度.

由式（14）可得$ K{\boldsymbol{\phi}} = {\omega ^2}M{\boldsymbol{\phi}} $，则设计变量α对固定频率的灵敏度为

$\lambda = \frac{{\partial {\omega ^2}}}{{\partial \alpha }} = {{\boldsymbol{\phi}} ^{\text{T}}}\left(\frac{{\partial K}}{{\partial \alpha }} - {\omega ^2}\frac{{\partial M}}{{\partial \alpha }}\right){\boldsymbol{\phi}} .$

(19)

设$ \Delta K $、$ \Delta M $是由参数变化$ \Delta \alpha $而产生的刚度矩阵和质量矩阵的增量，则由式（19）可得固定频率灵敏度的摄动公式为

$\lambda = \frac{{\Delta {\omega ^2}}}{{\Delta \alpha }} = {{\boldsymbol{\phi}} ^{\text{T}}}\left(\frac{{\Delta K}}{{\Delta \alpha }} - {\omega ^2}\frac{{\Delta M}}{{\Delta \alpha }}\right){\boldsymbol{\phi}} .$

(20)

图18为固定频率灵敏度分析结果，由图18可见：随着车轮质量的增大，固定频率灵敏度快速下降；随着车轴外径增加、支承刚度增大，固定频率灵敏度快速增大；随着扣件刚度增大，固定频率灵敏度几乎不变. 上述系统参数对固定频率的灵敏度从大到小的顺序依次为车轴外径>车轮质量>支承刚度>扣件刚度. 可见，轮对自身特性是车轮多边形磨耗的主要影响因素，而轨道特性的影响较小.

图 18 固定频率灵敏度分析结果

Figure 18. Fixed frequency sensitivity analysis results

下载: 全尺寸图片幻灯片

结合3.3节～3.6节的研究结果可见：通过调整系统参数改变固定频率的大小，使其避开被轮对转频整分，可避免车轮磨耗演化为车轮多边形. 具体参数调整可依据本文确定的灵敏度，也可以通过改变运行速度以改变轮对转频，来达到同样的目的.

4. 结　论

本文建立了车轮多边形磨耗发生和演化模型，表征了轮对-轨道转子动力学系统的振动特征和车轮多边形磨耗特征，揭示了车轮多边形发生和演化的规律. 主要结论如下：

1）基于轮对-轨道转子动力学理论建立轮轨接触动力学模型，结合磨耗功模型研究多边形磨耗，并用现场跟踪实测数据验证了模型的可靠性.

2）通过频谱分析和磨耗特征分析发现车轮多边形磨耗遵循“定频整分”的规律.

3）通过模态分析发现固定频率与车轮转动激发的轮对2阶弯曲振动的固有频率一致，而钢轨的柔性对固定频率的贡献不大.

4）通过破坏整分条件即可抑制多边形发展的思路，确定了改变固定频率的灵敏度顺序，为抑制车轮多边形磨耗提供了理论依据.

本文建立的高速车轮多边形磨耗模型中，未考虑动载荷及悬挂等的影响，而磨耗演化全过程预测研究也将是下一步进行的工作.

References

Proportional views

Seismic Responses of RC Frame Structures with Special-Shaped Columns and Girder-Transfer Story

Abstract

1. 车轮多边形磨耗发生与演化模型

1.1 轮对-轨道转子系统动力学模型

1.2 轮轨滚动接触力模型

1.3 轮对转子系统的动力学方程

1.4 车轮圆周磨耗深度模型

2. 多边形磨耗仿真分析

2.1 车轮多边形磨耗特征

2.2 偏心量的影响

2.3 运行速度和车轮直径的影响

2.4 现场跟踪实测验证

3. 固定频率的来源与影响

3.1 集中质量模型模态分析

3.2 实体有限元模型模态分析

3.3 扣件刚度对固定频率的影响

3.4 车轴外径对固定频率的影响

3.5 车轮质量对固定频率的影响

3.6 支承刚度对固定频率的影响

3.7 灵敏度分析

4. 结　论

References

Proportional views

Catalog

通讯作者: 陈斌, bchen63@163.com

Proportional views

Related

Seismic Responses of RC Frame Structures with Special-Shaped Columns and Girder-Transfer Story

Abstract

1. 车轮多边形磨耗发生与演化模型

1.1 轮对-轨道转子系统动力学模型

1.2 轮轨滚动接触力模型

1.3 轮对转子系统的动力学方程

1.4 车轮圆周磨耗深度模型

2. 多边形磨耗仿真分析

2.1 车轮多边形磨耗特征

2.2 偏心量的影响

2.3 运行速度和车轮直径的影响

2.4 现场跟踪实测验证

3. 固定频率的来源与影响

3.1 集中质量模型模态分析

3.2 实体有限元模型模态分析

3.3 扣件刚度对固定频率的影响

3.4 车轴外径对固定频率的影响

3.5 车轮质量对固定频率的影响

3.6 支承刚度对固定频率的影响

3.7 灵敏度分析

4. 结 论

References

Proportional views

Catalog

通讯作者: 陈斌, bchen63@163.com

Proportional views

Related

Export File

Citation

Format

Content

4. 结　论