Effect of vehicletrack vertical coupling vibrations on framemounted traction motor dynamics
Yuewei Yu^{1} , Leilei Zhao^{2} , Changcheng Zhou^{3}
^{1, 2, 3}School of Transportation and Vehicle Engineering, Shandong University of Technology, Zibo, China
^{3}Corresponding author
Journal of Vibroengineering, Vol. 22, Issue 1, 2020, p. 184196.
https://doi.org/10.21595/jve.2019.20900
Received 10 July 2019; received in revised form 22 September 2019; accepted 3 November 2019; published 15 February 2020
JVE Conferences
In order to reveal the effect of the vehicletrack vertical coupling vibrations on the framemounted traction motor dynamics, a vehicletrack vertical coupling dynamic model with considering the influence of the framemounted traction motor is established, and the correctness of the model is verified by real vehicle test. In case of the investigated vehicle model, the influences of the vehicletrack vertical coupling vibrations and the suspension parameters on the framemounted traction motor dynamics are discussed. The results show that the traction motor is significantly affected by the train system, when the motor is equivalent to the bogie frame mass, the phenomenon of underestimation exists to evaluate the vibration of the motor. In addition, suspension parameters have a great impact on the traction motor dynamics, rational selection suspension parameters can help to attenuate the vibration of the traction motor, and alleviate uneven the distribution of the air gap magnetic field of the traction motor.
 A vehicletrack vertical coupling dynamic model with considering the influence of the framemounted traction motor is established.
 The influences of the vehicletrack vertical coupling vibrations on the framemounted traction motor dynamics are discussed.
 The influences of the suspension parameters on the vertical vibration responses of the framemounted traction motor dynamics are discussed.
Keywords: train, traction motor, mounted on the bogie frame, electromechanical coupling, dynamic characteristics.
1. Introduction
Finding out the vibration mechanism of each subsystem is the key to ensure the safety and smooth operation of trains, and is also the premise of the design of railway vehicles’ main components [1, 2]. Many scholars have carried out a lot of research around this. For example, Dimitriu [3] compared the vertical vibration of trains under the rigid and elastic carbody from the comfort perspective by establishing a vertical dynamic model of railway vehicles with considering the elastic deformation of the car body. Lu et al. [4] proposed a solution method of the vehicletrack vertical coupling dynamic model based on the symplectic method, and discussed the influence of the track damping and vehicle speed on the train vertical dynamic responses. Nagai et al. [5] studied the influence of the passenger number and the passenger distribution on the railway vehicle’s car body vibration by establishing a passengerbody coupling vibration analysis model. However, for trains driven by traction motors, due to the limitation of computer performance, most of the early studies on train’s dynamic performance based on multibody dynamics always simplified the traction motor to the springup and/or springdown mass. In fact, for railway vehicles with power bogie system, the traction motor and the train system are a coupling system [6, 7], therefore, when ignoring the coupling between the traction motor and the train system, it will inevitably affect the analysis results of the train system’s vibration response, and also cannot analyze the dynamic characteristics of the traction motor during train operation.
In recent years, with the continuous deepening of research, researchers in various countries have carried out many related studies on trains and traction motors [2]. Alfi et al. [8] analyzed the influence of the motor connection stiffness on the critical speed of railway vehicles, and used dynamic absorbers to explain the impact of the drive system elastic suspension on vehicle lateral stability. Yao et al. [9] discussed the mechanism of the impact of elastic suspension on locomotive dynamic performance. The results show that when the lateral suspension frequency of the drive system is close to the wheelset rotational frequency, the locomotive has the best lateral dynamic performance. Besides, the small lateral suspension damping helps to improve the lateral stability of the locomotive and reduce the force of the wheel and the rail. Luo et al. [10] analyzed the influence of the suspension parameters of traction drive device on locomotive lateral dynamic performance. The results indicate that the mass of the traction drive device and the lateral stiffness of primary suspension have a great influence on the rigid suspension locomotive, but little influence on the elastic suspension locomotive. When locomotive’s running speed is above 160 km/h, the suspension mode of the traction device should adopt elastic suspension type. In addition, on the basis of considering the coupling vibration of the train system and the traction motor synthetically, Senini et al. [11] used computer to establish a dynamic simulation method for the wheelrail interactions, and studied the impact of the wheel slip/rotation effect and the wheel diameter on the load of the motor. Zhai [12, 13] studied the vibration characteristics of the traction motor and the dynamic interaction between the locomotive and the track by establishing a locomotivetrack vertical coupling dynamic model with considering the influence of the axlemounted traction motor, and then, he used the model to calculate and compare the dynamic characteristics of the locomotive to the track under the rigid and elastic suspension modes of the axlemounted traction motor, and studied the influence of the vehicle vibration on the dynamics of pantograph mounting system [14]. Chen et al. [15, 16] studied the effect of the locomotive dynamic characteristics on the root crack of the transmission system under traction drive conditions by establishing a heavyhaul electric locomotivetrack vertical longitudinal coupling dynamic model with considering the dynamic effect of gear transmission. Liu et al. [17] analyzed the dynamic force of the wheel and the rail under curve driving conditions by using the locomotivetrack vertical coupling dynamic model.
As stated above, although many scholars have done a lot of research on the dynamic characteristics of the train system and the traction motor under the coupling of the train and the motor, most of them are aimed at axlemounted traction motor. Therefore, it is impossible to accurately describe the actual vibration of the train system and the traction motor under the influence of the framemounted traction motor. It shows that, the vertical coupling vibration mechanism between the framemounted traction motor and the train system needs to be further studied. In this paper, according to the structural characteristics of a train with framemounted traction motor, a vehicletrack vertical coupling dynamic model with considering the influence of the framemounted traction motor is established. Then, by considering both the track irregularity excitation and the traction motor vibration characteristics, the vertical vibration responses of the train with framemounted traction motor are simulated, and the correctness of the dynamic model established is verified by the real vehicle test. In case of the investigated vehicle model, the influences of the vehicletrack vertical coupling vibrations and the suspension parameters on the framemounted traction motor dynamics are discussed. The research made in this paper can provide an effective model reference for the analysis of the dynamic characteristics of the framemounted traction motor and the optimization design of the motor suspension parameters.
2. Vehicletrack vertical coupling dynamic model with considering the framemounted traction motor
2.1. Physical model of the vertical coupling system
According to the actual structure of a train, the influence of the gear meshing in the transmission system and the radial clearance in the motor bearing are neglected, the effect of the lubricating oil film between the bearing inner/outer rings and the rolling body is not considered, and in the vertical direction, the bearing is equivalent to an elastic element, the motor is equivalent to a single mass springdamping vibration system, a vehicletrack vertical coupling vibration dynamic model with considering the influence of the framemounted traction motor is established, as shown in Fig. 1. Here, ${M}_{c}$, ${M}_{t}$, ${M}_{w}$, and ${M}_{m}$ are the half of the car body mass, the bogie frame mass, the wheelset and motor rotor total mass, and the motor stator and the shell total mass, respectively; ${J}_{c}$, ${J}_{t}$, and ${J}_{m}$ are the half of the car body moment inertia, the frame moment inertia, and the motor stator and the shell total moment inertia, respectively; ${C}_{pz}$ and ${C}_{sz}$ are the vertical equivalent damping of one side of the primary suspension and secondary suspension, respectively; ${K}_{pz}$ and ${K}_{sz}$ are the vertical equivalent stiffness of one side of the primary suspension and secondary suspension, respectively; ${C}_{m1}$ and ${C}_{m2}$ are the damping between the frame and driving system, respectively; ${K}_{m1}$ and ${K}_{m2}$ are the stiffness between the frame and driving system, respectively; ${K}_{mb}$ is the equivalent stiffness of the bearing; ${m}_{r}$, ${M}_{si}$, and ${M}_{bi}$ are the rail mass per unit length, the half of the sleeper mass, and the single side ballast mass, respectively; ${I}_{Y}$ is the rail section inertia, $E$ is the rail elastic modulus; ${K}_{pi}$, ${K}_{bi}$, ${K}_{wi}$, and ${K}_{fi}$ are the rail pad stiffness, the ballast support stiffness, the ballast shear stiffness, and the subgrade support stiffness, respectively; ${C}_{pi}$, ${C}_{bi}$, ${C}_{wi}$, and ${C}_{fi}$ are the rail pad damping, the ballast damping, the ballast shear damping, and the subgrade damping, respectively. In addition, ${z}_{c}$, ${z}_{t1}$${z}_{t2}$, ${z}_{m1}$${z}_{m4}$, and ${z}_{w1}$${z}_{w4}$ are the vertical displacements of the car body, the bogie frame, the traction motor, and the wheelset, respectively; ${\beta}_{c}$, ${\beta}_{t1}$${\beta}_{t2}$, and ${\beta}_{m1}$${\beta}_{m4}$ are the pitch displacements of the car body, the bogie frame, and the traction motor, respectively; ${z}_{r}$, ${z}_{s}$, and ${z}_{b}$ are the vertical displacements of the rail, the sleeper, and the ballast, respectively. ${P}_{1}$${P}_{4}$ are the wheelrail forces; ${z}_{v1}$${z}_{v4}$ are the track irregularities of each wheelset, and ${z}_{v2}={z}_{v1}\left(t{\tau}_{1}\right)$, ${z}_{v3}={z}_{v1}\left(t{\tau}_{2}\right)$, ${z}_{v4}={z}_{v1}\left(t{\tau}_{3}\right)$, where, $t$ denotes time, $\tau $ is the time lag, ${\tau}_{1}=2{L}_{t}/v$, ${\tau}_{2}=2{L}_{c}/v$, ${\tau}_{3}=2\left({L}_{t}+{L}_{c}\right)/v$, ${L}_{c}$ and ${L}_{t}$ are the half distances of the vehicle spacing and the bogie wheelbase, respectively, $v$ is the train running speed.
Fig. 1. Vehicletrack vertical coupling dynamic model with considering the influence of the framemounted traction motor
2.2. Mathematical model of the vertical coupling system
According to the Newton’s second law and the Elastic vibration theory, the vibration differential equations of each subsystem in the model shown in Fig. 1 can be derived as follows:
– Vertical vibration of the car body:
${K}_{sz}\left({z}_{c}+{L}_{c}{\beta}_{c}{z}_{t2}\right){C}_{sz}\left({\dot{z}}_{c}+{L}_{c}{\dot{\beta}}_{c}{\dot{z}}_{t2}\right).$
– Pitching vibration of the car body:
${K}_{sz}({z}_{c}+{L}_{c}{\beta}_{c}{z}_{t2}){C}_{sz}({\dot{z}}_{c}+{L}_{c}{\dot{\beta}}_{c}{\dot{z}}_{t2})]{L}_{c}.$
Vertical vibration of the front bogie frame:
${K}_{m1}({z}_{t1}{L}_{b}{\beta}_{t1}{z}_{m1}+{L}_{1}{\beta}_{m1}){K}_{pz}({z}_{t1}{L}_{t}{\beta}_{t1}{z}_{w1})$
${C}_{m1}({\dot{z}}_{t1}{L}_{b}{\dot{\beta}}_{t1}{\dot{z}}_{m1}+{L}_{1}{\dot{\beta}}_{m1}){C}_{pz}({\dot{z}}_{t1}{L}_{t}{\dot{\beta}}_{t1}{\dot{z}}_{w1})$
${C}_{sz}({\dot{z}}_{t1}{\dot{z}}_{c}+{L}_{c}{\dot{\beta}}_{c}){K}_{m2}[{z}_{t1}({L}_{b}{L}_{1}{L}_{2}){\beta}_{t1}{z}_{m1}{L}_{2}{\beta}_{m1}]$
${C}_{m2}[{\dot{z}}_{t1}({L}_{b}{L}_{1}{L}_{2}){\dot{\beta}}_{t1}{\dot{z}}_{m1}{L}_{2}{\dot{\beta}}_{m1}]$
${K}_{m2}[{z}_{t1}+({L}_{b}{L}_{1}{L}_{2}){\beta}_{t1}{z}_{m2}+{L}_{2}{\beta}_{m2}]$
${C}_{m2}[{\dot{z}}_{t1}+({L}_{b}{L}_{1}{L}_{2}){\dot{\beta}}_{t1}{\dot{z}}_{m2}+{L}_{2}{\dot{\beta}}_{m2}]$
${K}_{m1}({z}_{t1}+{L}_{b}{\beta}_{t1}{z}_{m2}{L}_{1}{\beta}_{m2}){K}_{sz}({z}_{t1}{z}_{c}+{L}_{c}{\beta}_{c})$
${C}_{m1}\left({\dot{z}}_{t1}+{L}_{b}{\dot{\beta}}_{t1}{\dot{z}}_{m2}{L}_{1}{\dot{\beta}}_{m2}\right).$
– Pitching vibration of the front bogie frame:
${K}_{pz}({z}_{t1}+{L}_{t}{\beta}_{t1}{z}_{w2}){L}_{t}{C}_{pz}({\dot{z}}_{t1}+{L}_{t}{\dot{\beta}}_{t1}{\dot{z}}_{w2}){L}_{t}$
$+{K}_{m1}({z}_{t1}{L}_{b}{\beta}_{t1}{z}_{m1}+{L}_{1}{\beta}_{m1}){L}_{b}+{C}_{m1}({\dot{z}}_{t1}{L}_{b}{\dot{\beta}}_{t1}{\dot{z}}_{m1}+{L}_{1}{\dot{\beta}}_{m1}){L}_{b}$
${K}_{m1}({z}_{t1}+{L}_{b}{\beta}_{t1}{z}_{m2}{L}_{1}{\beta}_{m2}){L}_{b}{C}_{m1}({\dot{z}}_{t1}+{L}_{b}{\dot{\beta}}_{t1}{\dot{z}}_{m2}{L}_{1}{\dot{\beta}}_{m2}){L}_{b}$
$+{K}_{m2}[{z}_{t1}({L}_{b}{L}_{1}{L}_{2}){\beta}_{t1}{z}_{m1}{L}_{2}{\beta}_{m1}]({L}_{b}{L}_{1}{L}_{2})$
$+{C}_{m2}[{\dot{z}}_{t1}({L}_{b}{L}_{1}{L}_{2}){\dot{\beta}}_{t1}{\dot{z}}_{m1}{L}_{2}{\dot{\beta}}_{m1}]({L}_{b}{L}_{1}{L}_{2})$
${K}_{m2}[{z}_{t1}+({L}_{b}{L}_{1}{L}_{2}){\beta}_{t1}{z}_{m2}+{L}_{2}{\beta}_{m2}]({L}_{b}{L}_{1}{L}_{2})$
${C}_{m2}\left[{\dot{z}}_{t1}+\left({L}_{b}{L}_{1}{L}_{2}\right){\dot{\beta}}_{t1}{\dot{z}}_{m2}+{L}_{2}{\dot{\beta}}_{m2}\right]\left({L}_{b}{L}_{1}{L}_{2}\right).$
– Vertical vibration of the rear bogie frame:
${K}_{m1}({z}_{t2}{L}_{b}{\beta}_{t2}{z}_{m3}+{L}_{1}{\beta}_{m3}){K}_{pz}({z}_{t2}+{L}_{t}{\beta}_{t2}{z}_{w4})$
${C}_{m1}({\dot{z}}_{t2}{L}_{b}{\dot{\beta}}_{t2}{\dot{z}}_{m3}+{L}_{1}{\dot{\beta}}_{m3}){C}_{pz}({\dot{z}}_{t2}+{L}_{t}{\dot{\beta}}_{t2}{\dot{z}}_{w4})$
${K}_{m2}[{z}_{t2}({L}_{b}{L}_{1}{L}_{2}){\beta}_{t2}{z}_{m3}{L}_{2}{\beta}_{m3}]$
${C}_{m2}[{\dot{z}}_{t2}({L}_{b}{L}_{1}{L}_{2}){\dot{\beta}}_{t2}{\dot{z}}_{m3}{L}_{2}{\dot{\beta}}_{m3}]$
${K}_{m2}[{z}_{t2}+({L}_{b}{L}_{1}{L}_{2}){\beta}_{t2}{z}_{m4}+{L}_{2}{\beta}_{m4}]$
${C}_{m2}[{\dot{z}}_{t2}+({L}_{b}{L}_{1}{L}_{2}){\dot{\beta}}_{t2}{\dot{z}}_{m4}+{L}_{2}{\dot{\beta}}_{m4}]$
${K}_{m1}({z}_{t2}+{L}_{b}{\beta}_{t2}{z}_{m4}{L}_{1}{\beta}_{m4}){K}_{sz}({z}_{t2}{z}_{c}{L}_{c}{\beta}_{c})$
${C}_{sz}\left({\dot{z}}_{t2}{\dot{z}}_{c}{L}_{c}{\dot{\beta}}_{c}\right){C}_{m1}\left({\dot{z}}_{t2}+{L}_{b}{\dot{\beta}}_{t2}{\dot{z}}_{m4}{L}_{1}{\dot{\beta}}_{m4}\right).$
– Pitching vibration of the rear bogie frame:
${K}_{pz}({z}_{t2}+{L}_{t}{\beta}_{t2}{z}_{w4}){L}_{t}{C}_{pz}({\dot{z}}_{t2}+{L}_{t}{\dot{\beta}}_{t2}{\dot{z}}_{w4}){L}_{t}$
$+{K}_{m1}({z}_{t2}{L}_{b}{\beta}_{t2}{z}_{m3}+{L}_{1}{\beta}_{m3}){L}_{b}+{C}_{m1}({\dot{z}}_{t2}{L}_{b}{\dot{\beta}}_{t2}{\dot{z}}_{m3}+{L}_{1}{\dot{\beta}}_{m3}){L}_{b}$
${K}_{m1}({z}_{t2}+{L}_{b}{\beta}_{t2}{z}_{m4}{L}_{1}{\beta}_{m4}){L}_{b}{C}_{m1}({\dot{z}}_{t2}+{L}_{b}{\dot{\beta}}_{t2}{\dot{z}}_{m4}{L}_{1}{\dot{\beta}}_{m4}){L}_{b}$
$+{K}_{m2}[{z}_{t2}({L}_{b}{L}_{1}{L}_{2}){\beta}_{t2}{z}_{m3}{L}_{2}{\beta}_{m3}]({L}_{b}{L}_{1}{L}_{2})$
$+{C}_{m2}[{\dot{z}}_{t2}({L}_{b}{L}_{1}{L}_{2}){\dot{\beta}}_{t2}{\dot{z}}_{m3}{L}_{2}{\dot{\beta}}_{m3}]({L}_{b}{L}_{1}{L}_{2})$
${K}_{m2}[{z}_{t2}+({L}_{b}{L}_{1}{L}_{2}){\beta}_{t2}{z}_{m4}+{L}_{2}{\beta}_{m4}]({L}_{b}{L}_{1}{L}_{2})$
– Vertical vibration of the first wheelset:
${K}_{mb}\left({z}_{w1}{z}_{m1}\right){K}_{H}\left({z}_{w1}{z}_{v1}\right){P}_{1}\left(t\right)+{P}_{0}.$
– Vertical vibration of the second wheelset:
${K}_{mb}\left({z}_{w2}{z}_{m2}\right){K}_{H}\left({z}_{w2}{z}_{v2}\right){P}_{2}\left(t\right)+{P}_{0}.$
– Vertical vibration of the third wheelset:
${K}_{mb}\left({z}_{w3}{z}_{m3}\right){K}_{H}\left({z}_{w3}{z}_{v3}\right){P}_{3}\left(t\right)+{P}_{0}.$
– Vertical vibration of the forth wheelset:
${K}_{mb}\left({z}_{w4}{z}_{m4}\right){K}_{H}\left({z}_{w4}{z}_{v4}\right){P}_{4}\left(t\right)+{P}_{0}.$
– Vertical vibration of the first traction motor:
${K}_{m2}[{z}_{m1}+{L}_{2}{\beta}_{m1}{z}_{t1}+({L}_{b}{L}_{1}{L}_{2}\left){\beta}_{t1}\right]{K}_{mb}({z}_{m1}{z}_{w1})$
${C}_{m2}[{\dot{z}}_{m1}+{L}_{2}{\dot{\beta}}_{m1}{\dot{z}}_{t1}+({L}_{b}{L}_{1}{L}_{2}\left){\dot{\beta}}_{t1}\right]$
${C}_{m1}({\dot{z}}_{m1}{L}_{1}{\dot{\beta}}_{m1}{\dot{z}}_{t1}+{L}_{b}{\dot{\beta}}_{t1})+{M}_{m}e{\omega}^{2}\mathrm{c}\mathrm{o}\mathrm{s}\left(\omega t\right).$
– Pitching vibration of the first traction motor:
${K}_{m2}[{z}_{m1}+{L}_{2}{\beta}_{m1}{z}_{t1}+({L}_{b}{L}_{1}{L}_{2}\left){\beta}_{t1}\right]{L}_{2}$
${C}_{m2}\left[{\dot{z}}_{m1}+{L}_{2}{\dot{\beta}}_{m1}{\dot{z}}_{t1}+\left({L}_{b}{L}_{1}{L}_{2}\right){\dot{\beta}}_{t1}\right]{L}_{2}.$
– Vertical vibration of the second traction motor:
${K}_{m2}[{z}_{m2}{L}_{2}{\beta}_{m2}{z}_{t1}({L}_{b}{L}_{1}{L}_{2}){\beta}_{t1}{K}_{mb}({z}_{m2}{z}_{w2}\left)\right]$
${C}_{m2}[{\dot{z}}_{m2}{L}_{2}{\dot{\beta}}_{m2}{\dot{z}}_{t1}({L}_{b}{L}_{1}{L}_{2}\left){\dot{\beta}}_{t1}\right]$
${C}_{m1}({\dot{z}}_{m2}+{L}_{1}{\dot{\beta}}_{m2}{\dot{z}}_{t1}{L}_{b}{\dot{\beta}}_{t1})+{M}_{m}e{\omega}^{2}\mathrm{c}\mathrm{o}\mathrm{s}\left(\omega t\right).$
– Pitching vibration of the second traction motor:
$+{C}_{m2}[{\dot{z}}_{m2}{L}_{2}{\dot{\beta}}_{m2}{\dot{z}}_{t1}({L}_{b}{L}_{1}{L}_{2}\left){\dot{\beta}}_{t1}\right]{L}_{2}$
$+{K}_{m2}[{z}_{m2}{L}_{2}{\beta}_{m2}{z}_{t1}({L}_{b}{L}_{1}{L}_{2}\left){\beta}_{t1}\right]{L}_{2}$
${C}_{m1}\left({\dot{z}}_{m2}+{L}_{1}{\dot{\beta}}_{m2}{\dot{z}}_{t1}{L}_{b}{\dot{\beta}}_{t1}\right){L}_{1}.$
– Vertical vibration of the third traction motor:
${K}_{m2}[{z}_{m3}+{L}_{2}{\beta}_{m3}{z}_{t2}+({L}_{b}{L}_{1}{L}_{2}\left){\beta}_{t2}\right]{K}_{mb}({z}_{m3}{z}_{w1})$
${C}_{m2}[{\dot{z}}_{m3}+{L}_{2}{\dot{\beta}}_{m3}{\dot{z}}_{t2}+({L}_{b}{L}_{1}{L}_{2}\left){\dot{\beta}}_{t2}\right]$
${C}_{m1}\left({\dot{z}}_{m3}{L}_{1}{\dot{\beta}}_{m3}{\dot{z}}_{t2}+{L}_{b}{\dot{\beta}}_{t2}\right)+{M}_{m}e{\omega}^{2}\mathrm{c}\mathrm{o}\mathrm{s}\left(\omega t\right).$
– Pitching vibration of the third traction motor:
${C}_{m2}[{\dot{z}}_{m3}+{L}_{2}{\dot{\beta}}_{m3}{\dot{z}}_{t2}+({L}_{b}{L}_{1}{L}_{2}\left){\dot{\beta}}_{t2}\right]{L}_{2}$
${K}_{m2}[{z}_{m3}+{L}_{2}{\beta}_{m3}{z}_{t2}+({L}_{b}{L}_{1}{L}_{2}\left){\beta}_{t2}\right]{L}_{2}$
$+{C}_{m1}\left({\dot{z}}_{m3}{L}_{1}{\dot{\beta}}_{m3}{\dot{z}}_{t2}+{L}_{b}{\dot{\beta}}_{t2}\right){L}_{1}.$
– Vertical vibration of the forth traction motor:
${K}_{m2}[{z}_{m4}{L}_{2}{\beta}_{m4}{z}_{t2}({L}_{b}{L}_{1}{L}_{2}\left){\beta}_{t2}\right]{K}_{mb}({z}_{m4}{z}_{w4})$
${C}_{m2}[{\dot{z}}_{m4}{L}_{2}{\dot{\beta}}_{m4}{\dot{z}}_{t2}({L}_{b}{L}_{1}{L}_{2}\left){\dot{\beta}}_{t2}\right]$
${C}_{m1}({\dot{z}}_{m4}+{L}_{1}{\dot{\beta}}_{m4}{\dot{z}}_{t2}{L}_{b}{\dot{\beta}}_{t2})+{M}_{m}e{\omega}^{2}\mathrm{c}\mathrm{o}\mathrm{s}\left(\omega t\right).$
– Pitching vibration of the forth traction motor:
$+{K}_{m2}[{z}_{m4}{L}_{2}{\beta}_{m4}{z}_{t2}({L}_{b}{L}_{1}{L}_{2}\left){\beta}_{t2}\right]{L}_{2}$
${C}_{m1}\left({\dot{z}}_{m4}+{L}_{1}{\dot{\beta}}_{m4}{\dot{z}}_{t2}{L}_{b}{\dot{\beta}}_{t2}\right){L}_{1}+{C}_{m2}\left[{\dot{z}}_{m4}{L}_{2}{\dot{\beta}}_{m4}{\dot{z}}_{t2}\left({L}_{b}{L}_{1}{L}_{2}\right){\dot{\beta}}_{t2}\right]{L}_{2}.$
– Vertical vibration of the rail:
${\sum}_{i=1}^{N}{K}_{pi}{z}_{k}\left({x}_{i}\right){\sum}_{h=1}^{NM}{z}_{h}\left({x}_{i}\right){q}_{h}\left(t\right)+{\sum}_{i=1}^{N}{C}_{pi}{z}_{k}\left({x}_{i}\right){\dot{z}}_{si}\left(t\right)$
$+{\sum}_{i=1}^{N}{K}_{pi}{z}_{k}\left({x}_{i}\right){z}_{si}\left(t\right)+{\sum}_{n=1}^{4}{P}_{n}\left(t\right){z}_{k}\left({x}_{wn}\right),\uff08k=1,\dots ,NM\uff09.$
– Vertical vibration of the sleeper:
$+{C}_{pi}{\sum}_{h=1}^{NM}{z}_{h}\left({x}_{i}\right){\dot{q}}_{h}\left(t\right)+{K}_{pi}{\sum}_{h=1}^{NM}{z}_{h}\left({x}_{i}\right){q}_{h}\left(t\right),\uff08i=1,\dots ,N\uff09.$
– Vertical vibration of the ballast:
${K}_{wi}{z}_{b(i+1)}\left(t\right)+{C}_{wi}{\dot{z}}_{b(i1)}\left(t\right)({K}_{bi}+{K}_{fi}+2{K}_{wi}){z}_{bi}\left(t\right)+{K}_{wi}{z}_{b(i1)}\left(t\right),$
$\uff08i=1,\dots ,N\uff09.$
Here, ${L}_{1}$ is the distance between the driving system mass center and the end beam suspension point, ${L}_{2}$ is the distance between the driving system mass center and the crossbeam suspension point, ${L}_{b}$ is the half of the frame length, $e$ is the rotor eccentricity, $\omega $ is the speed of the rotor. ${P}_{0}$ is the static force of the wheel and the rail, ${P}_{n}\left(t\right)$ is the vertical force of the wheel and the rail, $NM$ is the rail modal order, $N$ is the sleeper number; ${q}_{k}\left(t\right)$ is the regular vibration mode coordinate of the rail, in which $k$ is the mode order number; ${z}_{h}\left({x}_{i}\right)$ is the transposed matrix of ${z}_{k}\left({x}_{i}\right)$; the expression of ${z}_{k}\left({x}_{i}\right)$ and ${P}_{n}\left(t\right)$ can be written as [18]:
${P}_{n}\left(t\right)=\left\{\begin{array}{l}{\left\{\frac{1}{G}\left[{z}_{wn}\right(t){z}_{r}({x}_{wn},t){z}_{vn}(t\left)\right]\right\}}^{3/2},{z}_{wn}\left(t\right){z}_{r}({x}_{wn},t){z}_{vn}\left(t\right)\ge 0,\\ \\ 0,{z}_{wn}\left(t\right){z}_{r}({x}_{wn},t){z}_{vn}\left(t\right)0,\end{array}\right.$
where, $G$ is the contact constant of the wheel and the rail, for the taper tread wheel $G=$4.57$R$^{0.149}×10^{8} (m/N^{2/3}), for the wear type tread wheel $G=$3.86$R$^{0.115}×10^{8} (m/N^{2/3}), $R$ is the wheel radius.
2.3. Numerical solution of the system dynamic response
By analyzing the mathematical model of the vertical coupling system, it can be seen that, the equation group is a large, complex nonlinear dynamic differential equation group (with 18+$NM$+2$N$ DOFs). It is difficult to solve it by the analytic method. Therefore, we need to use numerical integration method to solve it. According to Eqs. (121), they can be expressed as the following matrix and vector form:
where, $\mathbf{M}$, $\mathbf{C}$, and $\mathbf{K}$ are the mass, damping, and stiffness matrices of the vertical coupling system, respectively; $\left\{\ddot{\mathbf{X}}\right\}$, $\left\{\dot{\mathbf{X}}\right\}$, and {$\mathbf{X}$} are the generalized acceleration, velocity, and displacement vectors, respectively; {$\mathbf{P}$} is the force vector.
According to Eq. (22), using the Newmark$\beta $ implicit integration algorithm, the timedomain responses of the vehicletrack vertical coupling vibration dynamic model with considering the influence of the framemounted traction motor can be calculated. The iterative process of solving the equation is as follows:
where, $\mathrm{\Delta}t$ is the time integration step, $m$ represents the iterative step; $\beta $ and $\gamma $ are the parameters to control the accuracy and stability of the integration method.
3. Test verification of the vertical coupling dynamic model
In order to verify the correctness of the vertical coupling dynamic model, the field measurements have been made by Shandong Zhiheng Vibration Damping Tech Co., Ltd, on a highspeed train which has been widely employed in China, and the analysis and test results are compared and analyzed. The wheel tread type of the train is wear type. The parameters of the traction motor and the train system can be found in Table 1, and the parameters of the ballasted track employed in this paper can be referenced to Zhai et al. [19]. The train running speed is 200 km/h, the corresponding motor rotor speed is 392 r/s, the eccentricity of the rotor is 1.5×10^{4} m, and the sampling frequency is 1 000 Hz.
Table 1. System parameters values
Symbol

Unit

Value

Symbol

Unit

Value

${M}_{c}$

kg

17 000

${K}_{pz}$

N/m

275 000

${M}_{t}$

kg

1 500

${K}_{sz}$

N/m

200 000

${M}_{w}$

kg

700

${K}_{mb}$

N/m

1.0×10^{4}

${M}_{m}$

kg

220

${K}_{m1}$

N/m

1.0×10^{6}

${J}_{c}$

kg·m^{2}

1 138 500

${K}_{m2}$

N/m

1.0×10^{6}

${J}_{t}$

kg·m^{2}

1 355

${L}_{c}$

m

9.0

${J}_{m}$

kg·m^{2}

13.7

${L}_{t}$

m

1.2

${C}_{pz}$

N·s/m

30 000

${L}_{1}$

m

0.3

${C}_{sz}$

N·s/m

40 000

${L}_{2}$

m

0.5

${C}_{m1}$

N·s/m

100

${L}_{b}$

m

1.8

${C}_{m2}$

N·s/m

100

$R$

m

0.457 5

Track irregularity is the excitation source of the wheelrail system and the main cause of the vehicle vibration and wheelrail force. It has an important impact on the dynamic characteristics of railway vehicles [20, 21]. Among them, the track vertical irregularity is the main cause of the trains’ vertical vibration. In order to make the simulation conditions as consistent as possible with the line conditions in the real vehicle test, the German track spectrum is adopted for the track irregularity. Since the track irregularity is usually expressed in the form of power spectral density, the timefrequency conversion method provided by Zhai et al. [19] is used to transform the track irregularity into a timedomain sample as the input of the system, and the vertical vibration responses of the train system and the traction motor are analyzed. Here, the spatial frequency analytical expression of the power spectral density of the track vertical irregularity is [22]:
where, ${S}_{v}\left(\mathrm{\Omega}\right)$ is the power spectral density of the track vertical irregularity, $\mathrm{\Omega}$ is the spatial frequency of the track irregularity, ${A}_{v}$ is the roughness coefficient of the rail; ${\mathrm{\Omega}}_{c}$ and ${\mathrm{\Omega}}_{r}$ are the truncated space frequencies. Each known parameter value is shown in Table 2. Here, the low interference spectrum is suitable for the vehicle speed of greater than or equal to 250 km/h, the high interference spectrum is suitable for the vehicle speed below 250 km/h.
Table 2. Parameters of the track vertical irregularity
Track level

Parameter


${A}_{v}$
_{}/(m^{2}·rad/m)

${\mathrm{\Omega}}_{c}$
_{}/(rad/m)

${\mathrm{\Omega}}_{r}$
_{}/(rad/m)


Low interference

4.032×10^{7}

0.824 6

0.020 6

High interference

10.8×10^{7}

0.824 6

0.020 6

Fig. 2 shows the simulation analysis and test results of the vertical vibration accelerations of the traction motor and the car body, respectively.
From Fig. 2, it can be seen that, the vertical vibration acceleration responses of the traction motor and the car body calculated by simulation are in good agreement with the vibration responses measured by the actual vehicle test, but the peak values are slightly smaller than the measured results. In which, the maximum values of the vertical vibration accelerations of the traction motor and the car body simulated are 41.157 m/s^{2} and 1.429 m/s^{2}, respectively, and the measured maximum values are 43.661 m/s^{2} and 1.479 m/s^{2}, respectively. Thus, the vertical vibration accelerations of the traction motor and the car body calculated by the dynamic model established in this paper are in good agreement with the measured data, which shows that the vehicletrack vertical coupling vibration dynamic model with considering the influence of the framemounted traction motor proposed is correct.
Fig. 2. Vertical random vibration response of the railway vehicle: a) the vertical vibration acceleration of the traction motor, b) the vertical vibration acceleration of the car body
4. The influence of the vehicletrack vertical coupling vibrations on the framemounted traction motor dynamics
In order to find out the influences of the vehicletrack vertical coupling vibrations on the framemounted traction motor dynamics, by using the dynamic model established in Section 2, taking the highspeed train described in Section 3 as an example, the vertical vibration responses of the train with or without framemounted traction motor under the excitation of track irregularity shown in Eq. (23) are analyzed. In which, the simulation time step is 0.002 s, the simulation time is 10 s, and the train speed is 250 km/h. In addition, in the process of simulation calculation, except the traction motor, other dynamic parameters and line parameters of the train are identical.
4.1. Influence of train system on the vibration response of the traction motor
Figs. 34 show the power spectral density curves of the vertical and pitching accelerations of the traction motor, and the relative clearance between the stator and rotor of the motor with or without the train system obtained by simulation. Since the similarity of the vibration laws of the four motors, only the analysis results of the first motor are given here. Note that, the influence of the random track irregularity disturbance is neglected when the train system is not considered.
Fig. 3. Power spectral density curve of the traction motor vibration acceleration: a) the vertical acceleration of the traction motor, b) the pitching acceleration of the traction motor
Fig. 3 and Fig. 4 show that, the vertical vibration responses of the traction motor increase significantly after considering the train system. When considering the train system, the root mean square values of the vertical and pitching vibration accelerations of the traction motor are 23.095 m/s^{2} and 23.121 rad/s^{2}, respectively; the root mean square value of the relative clearance between the stator and rotor of the motor is 7.650×10^{3} m. When the train system is not considered, the root mean square values of the vertical and pitching vibration accelerations of the traction motor are 17.185 m/s^{2} and 17.185 rad/s^{2}, respectively; the root mean square value of the relative clearance between the stator and rotor of the motor is 1.013×10^{4} m. The results show that the traction motor is significantly affected by the train system. Therefore, the influence of the train system should be taken into account when analyzing the vibration characteristics of the traction motor.
Fig. 4. Power spectral density curve of the relative clearance between the stator and rotor
4.2. Influence of electromechanical coupling on the vibration response of the traction motor
Fig. 5 and Fig. 6 show the power spectral density curves of the vertical and pitching acceleration of the traction motor, and the relative clearance between the stator and rotor of the motor with or without the electromechanical coupling obtained by simulation. Since the similarity of the vibration laws of the four motors, only the analysis results of the first motor are given here. Here, when the electromechanical coupling is not considered, the motor is equivalent to the frame mass.
Fig. 5. Power spectral density curve of the traction motor vibration acceleration: a) the vertical acceleration of the traction motor, b) the pitching acceleration of the traction motor
From Fig. 5 and Fig. 6, it can be seen that, the vertical vibration responses of the traction motor increases with considering the electromechanical coupling. When considering the electromechanical coupling, the root mean square values of the vertical and pitching vibration accelerations of the traction motor are 23.095 m/s^{2} and 23.121 rad/s^{2}, respectively; the root mean square value of the relative clearance between the stator and rotor of the motor is 7.650×10^{3} m. When the electromechanical coupling is not considered, the root mean square values of the vertical and pitching vibration accelerations of the traction motor are 22.579 m/s^{2} and 21.761 rad/s^{2}, respectively; the root mean square value of the relative clearance between the stator and rotor of the motor is 0 m. The results show that when the motor is equivalent to the bogie frame mass, the phenomenon of underestimation exists to evaluate the vibration of the motor. Therefore, the coupling effect between the train system and the traction motor should be considered in the analysis of the vibration characteristics of the motor.
Fig. 6. Power spectral density curve of the relative clearance between the stator and rotor
5. The influence of the suspension parameters on the framemounted traction motor dynamics
In order to find out the influences of the suspension parameters on the framemounted traction motor dynamics, taking the highspeed train described in Section 3 as an example, the influences of the stiffness parameters ${K}_{m1}$, ${K}_{m2}$ and damping parameters ${C}_{m1}$ and ${C}_{m2}$ on the vibration responses of the traction motor are studied by using the coupling system model.
Fig. 7. The influences of the suspension parameters on the vibration responses of the traction motor: a) the influence of ${K}_{m1}$, b) the influence of ${K}_{m2}$, c) the influence of ${C}_{m1}$, d) the influence of ${C}_{m2}$
Fig. 7 shows the change of the root mean square values of the vertical and pitching vibration accelerations of the traction motor, and the relative clearance between the motor stator and rotor with the stiffness parameters ${K}_{m1}$, ${K}_{m2}$ and damping parameters ${C}_{m1}$ and ${C}_{m2}$, respectively. In order to represent each response in a same figure, ${\ddot{z}}_{m1}$, ${\ddot{\beta}}_{m1}$, and ${z}_{m1}$${z}_{w1}$ are treated as follows: ${\sigma}_{1}={\sigma}_{{\ddot{z}}_{m1}}/{\sigma}_{{\ddot{z}}_{m1b}}$, ${\sigma}_{2}={\sigma}_{{\ddot{\beta}}_{m1}}/{\sigma}_{{\ddot{\beta}}_{m1b}}$, ${\sigma}_{3}={\sigma}_{({z}_{m1}{z}_{w1})}/{\sigma}_{({z}_{m1b}{z}_{w1b})}$, where, $b$ represents the vibration responses under the original vehicle parameters. Since the similarity of the vibration laws of the four motors, only the analysis results of the first motor are given here.
As can be seen from Fig. 7, with the increase of ${K}_{m1}$ and ${K}_{m2}$, ${\sigma}_{1}$ and ${\sigma}_{2}$ decrease first, then increase and finally stabilize, ${\sigma}_{3}$ decreases first and then stabilizes; with the increase of ${C}_{m1}$ and ${C}_{m2}$, ${\sigma}_{1}$ and ${\sigma}_{2}$ decrease first and then increase, and ${\sigma}_{3}$ hardly changes. It can be seen that, under the condition of a given mass of the motor, there are optimal ${K}_{m1}$, ${K}_{m2}$, ${C}_{m1}$, and ${C}_{m2}$ values which minimize the vertical and pitching vibration accelerations of the motor, and the values of ${K}_{m1}$ and ${K}_{m2}$ in a certain range of stiffness can help to reduce the relative clearance between the stator and the rotor of the motor, that is, to improve the uneven distribution of the air gap magnetic field of the motor. Among them, the sensitivity degree of ${\sigma}_{1}$ and ${\sigma}_{2}$ to each suspension parameter is in the order of ${K}_{m1}$, ${K}_{m2}$, ${C}_{m1}$, and ${C}_{m2}$.
6. Conclusions
In this paper, a vehicletrack vertical coupling dynamic model with considering the influence of the framemounted traction motor is established. And then, by considering both the track irregularity excitation and the traction motor vibration characteristics, the vertical vibration responses of the train with framemounted traction motor are simulated, and the correctness of the dynamic model established is verified by the real vehicle test. According to the results, the following conclusions can be drawn:
1) Traction motor is significantly affected by the train system, when the motor is equivalent to the bogie frame mass, the phenomenon of underestimation exists to evaluate the vibration of the motor, therefore, when analyzing the vibration characteristics of the traction motor, the coupling effect between the train system and the traction motor should be considered.
2) Suspension parameters have a great impact on the vibration responses of the traction motor. There are optimal ${K}_{m1}$, ${K}_{m2}$, ${C}_{m1}$, and ${C}_{m2}$ values which minimize the vertical and pitching vibration accelerations of the motor. Rational selection of suspension parameters will help to attenuate the vibration of the traction motor, and alleviate uneven the distribution of air gap magnetic field of the traction motor. The dynamic model can be further used to optimize the stiffness and damping parameters of the traction motor suspension system.
References
 Evans J., Berg M. Challenges in simulation of rail vehicle dynamics. Vehicle System Dynamics, Vol. 47, Issue 8, 2009, p. 10231048. [Publisher]
 Zhai W. M. VehicleTrack Coupled Dynamics. Fourth Edition, Science Press, Beijing, 2015. [CrossRef]
 Dimitriu M. On the critical points of vertical vibration in a railway vehicle. Archive of Mechanical Engineering, Vol. 61, Issue 4, 2014, p. 609625. [Publisher]
 Lu F., Kennedy D., Williams F. W., et al. Symplectic analysis of vertical random vibration for coupled vehicle track systems. Journal of Sound and Vibration, Vol. 317, Issues 12, 2008, p. 236249. [Publisher]
 Nagai M., Yoshida H., Tohtake T., et al. Coupled vibration of passenger and lightweight carbody in consideration of humanbody biomechanics. Vehicle System Dynamics, Vol. 44, 2006, p. 601611. [Publisher]
 Wu X. W., Chi M. R., Gao H. The study of postderailment dynamic behavior of railway vehicle based on running tests. Engineering Failure Analysis, Vol. 44, Issue 5, 2014, p. 382399. [Publisher]
 Wu X. W., Chi M. R., Gao H. Postderailment dynamic behaviour of a highspeed train under earthquake excitations. Engineering Failure Analysis, Vol. 64, 2016, p. 97110. [Publisher]
 Alfi S., Mazzola L., Bruni S. Effect of motor connection on the critical speed of highspeed railway vehicles. Vehicle System Dynamics, Vol. 46, 2008, p. 201214. [Publisher]
 Yao Y., Zhang K. L., Luo S. H. The mechanism of dynamic effect of the driving system elastically mounted on the bogie frame of locomotive. Journal of Vibration Engineering, Vol. 25, Issue 5, 2012, p. 481487. [CrossRef]
 Luo Y., Jin D. C. Research on the rules of suspension parameters to driving equipments suspended in bogie frames. China Railway Science, Vol. 4, 2007, p. 7882. [CrossRef]
 Senini F., Flinders S., Oghanna W. Dynamic Simulation of WheelRail Interaction for Locomotive Traction Studies. Railroad Conference, Pittsburgh, 1993. [CrossRef]
 Zhai W. M. Locomotiveandtrack system coupling dynamics and its application to the study of hoeomotwe performance. China Railway Science, Vol. 2, 1996, p. 5873. [CrossRef]
 Zhai W. M. Study of wheel/rail dynamic performance of electric drive locomotive. Electric Drive for Locomotives, Vol. 6, 1996, p. 1013. [CrossRef]
 Zhai W. M., Cai C. B. Effect of locomotive vibrations on pantographcatenary system dynamics. Vehicle System Dynamics, Vol. 29, 1998, p. 4758. [Publisher]
 Chen Z. G., Zhai W. M., Wang K. Y. Dynamic investigation of a locomotive with effect of gear transmissions under tractive conditions. Journal of Sound and Vibration, Vol. 408, 2017, p. 220233. [Publisher]
 Chen Z. G., Zhai W. M., Wang K. Y. Vibration feature evolution of locomotive with tooth root crack propagation of gear transmission system. Mechanical Systems and Signal Processing, Vol. 115, 2019, p. 2944. [Publisher]
 Liu P. F., Zhai W. M., Wang K. Y., et al. Theoretical and experimental study on vertical dynamic characteristics of sixaxle heavyhaul locomotive on curve. Transport, Vol. 33, 2018, p. 291301. [Publisher]
 Jin X. S., Wen Z. F., Wang K. Y., et al. Threedimensional train–track model for study of rail corrugation. Journal of Sound and Vibration, Vol. 293, Issues 35, 2006, p. 830855. [Publisher]
 Zhai W. M., Wang K. Y., Cai C. B. Fundamentals of vehicletrack coupled dynamics. Vehicle System Dynamics, Vol. 47, Issue 11, 2009, p. 13491376. [Publisher]
 Sowiński B. Interrelation between wavelengths of track geometry irregularities and rail vehicle dynamic properties. The Archives of Transport, Vols. 2526, Issue 12, 2013, p. 97108. [CrossRef]
 Sowiński B., Kostrzewski M. Railway vehicle transient vibration as a factor of track irregularities formation. Proceedings of the 13th Mini Conference on Vehicle System Dynamics, Identification and Anomalies, 2012. [CrossRef]
 Mao J. F., Yu Z. W., Xiao Y. J. Random dynamic analysis of a trainbridge coupled system involving random system parameters based on probability density evolution method. Probabilistic Engineering Mechanics, Vol. 46, 2016, p. 4861. [Publisher]