Effect of Control Measures on Wheel/Rail Noise When the Vehicle Curves

: This paper developed a time domain simpliﬁed model to study the effect of control measures on wheel/rail noise when the vehicle curves. The time domain model consists of two parts, one being a vehicle-track coupling dynamic model for wheel/rail interaction, the other being a transient ﬁnite element and boundary element domain model for vibration and sound radiation. Wheel/rail noise under wheel/rail lateral creep force is predicted for a narrowly curved section of a conventional underground railway, and compared with measurement. Based on the developed model, the effect of wheel/rail friction modiﬁcation on squeal noise is investigated. In addition, effectiveness of resilient wheel and embedded track to control curve squeal noise are also assessed.


Introduction
Railway squeal noise is difficult to study because of its mechanism, wheel/rail interaction and random acoustic behaviour. Squeal generates as a railway vehicle negotiates a narrow curve, producing a tonal noise of high amplitude. Typically, the noise is radiated from the wheels and has most of its energy between 400 Hz and 10 kHz (in the frequency band where human beings are most sensitive to noise). In the process of curving, a series of railway wheel modes (particularly the wheel's axial modes) are easy to be excited, which may lead to curve squeal [1][2][3]. Known practical solutions for railway squeal in a particular curve include friction modification [4,5] and wheel damping treatments [6][7][8][9]. These solutions can reduce the occurrence of squeal for some special circumstance.
Curve squeal noise has received increased attention by means of simulations and tests. Most previous research considered that the mechanism of railway curve squeal is mainly a kind of self-excited vibration due to the nonlinear creep forces [1, 3,[10][11][12][13][14][15][16][17][18][19]. Rudd [1] found that the main excitation mechanisms of squeal were related to lateral direction stick-slip phenomena which introduced unstable wheel/rail forces and vibrations. Slope of the creep force-creepage curve becomes negative instead of a constant value when the creepages reach a high level. Many models to simulate curve squeal were developed based on Rudd's squeal theory. Most of the models begin with a characteristic of the creepage-dependent friction coefficient. Van Ruiten [10], for example, applied Rudd's model to trams. In a review article, Remington [11] described the state of the art of railway curve squeal up to 1985. He suggested that this model should include finite element models of railway wheel dynamics, numerical models for the dynamics of bogies in curves and details of the friction coefficient versus creepage. Since then, Schneider et al. [12] developed a modal model for the wheel and calculated its self-excited vibrations in the time domain through introducing Kraft's formula for the falling friction coefficient in the sliding region. Fingberg [13] included a finite element model of the wheelset and a modal model of the rail. The creep force formula is derived from Kalker's theory [20] of rolling contact, extended to include a falling regime based on Kraft's formula to resemble experimental findings more closely at large creepages. Fingberg's work has been extended further by Périard [14] who included vehicle dynamics simulations in a time-domain curve squeal calculation. Xie [21] introduced a simplified friction law into SIMPACK using a modified FASTSIM for curve squeal calculations based on Fingberg's work and TWINS. Heckl [16] developed an annular disc model and analyzed its vibration response in time domain. Heckl's work made a contribution to understanding of the disc vibration and its modal shapes when squeal noise occurs due to lateral creepage. In order to consider railway wheel geometry, Heckl's model needs to be improved. De Beer et al. [16] also presented a curve squeal model considering lateral creepages. The difference was that his model was developed in frequency domain. Noise radiation was predicted using the standard components of the TWINS software [22]. The model of de Beer was extended by Monk-Steel and Thompson [17] by including the case of wheel flange contact, which usually occurred at the leading wheel at high rail. Chiello et al. [18] presented a curve squeal model including not only tangential dynamics (friction force) but also normal dynamics on the wheel/rail contact zone. The model can predict the unstable wheel modes and the corresponding squeal level and spectrum. Huang [23] reviewed the curve squeal models, and developed a curve squeal model including a self-excited loop and relationship between wheel/rail motions and contact forces. A number of authors have proposed an alternative mechanism based on 'mode coupling', which has been explained in a simplified form by Hoffmann et al. [24,25], see also [26]. More recently Thompson comprehensively reviewed two mechanisms ('falling friction' and 'mode coupling'), experimental and theoretical work in the field of curve squeal and discussed mitigation measures in terms of these two mechanisms [27].
This paper developed a time domain simplified model to study the effect of control measures on wheel/rail noise when the vehicle curves. This model, different from the self-excited vibration model, mainly analyzed the wheel/rail noise characteristic under lateral excitation. Wheel/track vibration and sound pressure can approximate predicted for a vehicle negotiating a narrow curve. The effects of lubricant between wheel/rail interface, elastic wheel, and embedded track on wheel/rail noise reduction under lateral excitation are investigated. Figure 1 shows a three-dimensional dynamic model of a metro vehicle which coupled with a conventional ballastless track which can predict the curve squeal when a vehicle negotiates narrow curves. The dynamic model is based on the classical dynamic theory, which was also successfully applied to many research topics in rail transit field [28][29][30][31][32][33]. In this paper, the car-body and two bogies are modelled as rigid bodies. The primary and secondary suspensions are modelled as spring-damp elements. The track system is a triple-layers model as demonstrated in [30], in which the couple of rails are modelled as Timoshenko beams. The explicit integration method in [34] is used to solve dynamic equations. Wheel-rail normal force is calculated based on non-linear Hertzian contact theory. Non-linear Hertzian contact force is simulated by Equation (1). The normal force depends on the relationship between normal loading and wheel, rail deformation if the normal compression at the wheel/rail contact point is larger than zero when wheel and rail are in contact. The normal force is zero if the normal compression at the wheel/rail contact point is not larger than zero when the wheel and rail separates.

Vehicle/Track Interaction Model with Falling Friction Coefficients
where GHertz is the contact constant of wheel and rail (m/N 2/3 ). For LM wheel tread profile and CN60 rail, GHertz = 3.86 R −0.115 × 10 −8 . R is the wheel rolling radius. Zwrnc (t) is the normal compression at the wheel/rail contact point.
The tangential creep forces of wheel and rail are calculated based on Kalker's simplified creep theory, in which the falling down friction law [13] is introduced. The contact area is discretized and the tangential wheel-rail creep force is solved through iteration from the leading part of the contact area to its trailing part. Thus, the tangential wheel-rail creep forces can be written as follows: ( 1) where pxi and pyi are longitudinal and lateral tangential wheel-rail creep force. ξx, ξy, ξz are longitudinal, lateral and spin creepage, respectively. Δx is the longitudinal width of discrete grid. L1, L2, L3 are flexibility coefficients and can be written as follows: Wheel-rail normal force is calculated based on non-linear Hertzian contact theory. Non-linear Hertzian contact force is simulated by Equation (1). The normal force depends on the relationship between normal loading and wheel, rail deformation if the normal compression at the wheel/rail contact point is larger than zero when wheel and rail are in contact. The normal force is zero if the normal compression at the wheel/rail contact point is not larger than zero when the wheel and rail separates.
where G Hertz is the contact constant of wheel and rail (m/N 2/3 ). For LM wheel tread profile and CN60 rail, G Hertz = 3.86 R −0.115 × 10 −8 . R is the wheel rolling radius. Z wrnc (t) is the normal compression at the wheel/rail contact point. The tangential creep forces of wheel and rail are calculated based on Kalker's simplified creep theory, in which the falling down friction law [13] is introduced. The contact area is discretized and the tangential wheel-rail creep force is solved through iteration from the leading part of the contact area to its trailing part. Thus, the tangential wheel-rail creep forces can be written as follows: where p xi and p yi are longitudinal and lateral tangential wheel-rail creep force. ξ x , ξ y , ξ z are longitudinal, lateral and spin creepage, respectively. ∆x is the longitudinal width of discrete grid. L 1 , L 2 , L 3 are flexibility coefficients and can be written as follows: 11 , L 2 = 8a 3GC 22 , where a and b are the lengths of longitudinal and lateral semi-axis, respectively. G is the shear modulus of wheel/rail. C ij is the coefficient calculated by Kalker in 1967 [20]. The tangential wheel-rail creep forces at any point of contact area can be written as follows: where µ is the wheel/rail friction coefficient, and p is the distribution of wheel/rail normal force, the particle in contact area is in an adhesion state. The tangential wheel-rail creep forces can be calculated by Equation (2).
the particle in contact area is in a slide state. The tangential wheel-rail creep forces can be reduced to Equation (7): The wheel/rail friction coefficient µ can be modeled by a constant [35], Fingberg law (Equation (8)) [13] or piecewise friction coefficient (Equation (9)) [21].
where v(x, y) is the slip velocity between rail and wheel, µ s is the static friction coefficient.
where µ d is the sliding friction coefficient, v c the critical sliding velocity, and v(x, y) is the slip velocity in the local coordinate system. Figure 2 gives the relationship of Creep force and Creepage based on three kinds of friction coefficients. The Fingberg model and piecewise model both consider the falling down friction coefficient. The piecewise friction coefficient is simplified from Fingberg model. The falling friction coefficient needs to be taken into account in the process of curve squeal simulation. This study used the Fingberg' model.
In this study, a special type of track, i.e., embedded track, is investigated. Figure 3 shows the dynamic model of filling material (including elastomer and rail pad) in the groove. The filling material is glued to the rail and slab, and it provides elastic support and damping. There is no relative motion between filling material and rail/slab. In the dynamic simulation of embedded track, the filling material of the embedded track is approximately simulated as mass-spring-damping system. Compared with the spring-damping simulation, the effects of the inertia of elastomer can be considered in the mass-spring-damping system in this paper. The slab of the embedded track is simulated by 3D finite element in the vehicle-track coupling dynamic model.           Lateral motion equation of elastomer on the left side .
Vertical motion equation of elastomer on the left side .
Lateral motion equation of elastomer on the right side .
Vertical motion equation of elastomer on the right side .
where M filli is the mass of the i-th elastomer on one side (left or right) of the rail. ..
are lateral force between slab and elastomer, and lateral force between elastomer and rail. F l f rv(L,R)i , F l f sv(L,R)i , F l f f (L,R)i and F l f r(L,R)i are vertical force between elastomer and rail, vertical force between elastomer and slab, vertical force between the i-th and the (i − 1)-th elastomers, and vertical force between the i-th and the (i + 1)-th elastomers on the left side. F r f rv(L,R)i , F r f sv(L,R)i , F r f f (L,R)i and F r f r(L,R)i are vertical force between elastomer and rail, vertical force between elastomer and slab, vertical force between the i-th and the (i − 1)-th elastomers, and vertical force between the i-th and the (i + 1)-th elastomers on the right side.
In addition to the embedded track, the resilient wheel is also examined. Figure 5 shows the dynamic model of resilient wheel. The lateral and vertical motions of resilient wheel are written in Equations (14)- (19).
where M w1 and M w2 are the mass of the rim and the rest parts of the wheelset as shown in Figure 5.  Z 2wi are velocities of rim and the rest parts of the wheelset. Y 1wL,Ri , Y 2wi , Z 1wL,Ri and Z 2wi are displacements of rim and the rest parts of the wheelset. V is the running speed. R wi is the curve radii of i-th wheel. a 0 is the half of the distance between left and right nominal rolling circles.
.. φ sewi and φ sewi are the angular acceleration and angular of super-elevation. F eyL,Ri , F wryL,Ri , and F fyL,Ri are lateral forces due to the stiffness and damping of rubber, due to wheel/rail interaction, and due to primary suspension system. F ezL,Ri , F wrzL,Ri , and F fzL,Ri are vertical forces due to the stiffness and damping of rubber, due to wheel/rail interaction, and due to primary suspension system. Appl. Sci. 2017, 7, x;doi: www.mdpi.com/journal/applsci elastomer and rail, vertical force between elastomer and slab, vertical force between the i-th and the (i − 1)-th elastomers, and vertical force between the i-th and the (i + 1)-th elastomers on the left side.
rfr L R i F are vertical force between elastomer and rail, vertical force between elastomer and slab, vertical force between the i-th and the (i − 1)-th elastomers, and vertical force between the i-th and the (i + 1)-th elastomers on the right side.
In addition to the embedded track, the resilient wheel is also examined. Figure 5 shows the dynamic model of resilient wheel. The lateral and vertical motions of resilient wheel are written in Equations (14)- (19).

The Transient Finite and Boundary Element Model for Wheel/Rail Vibration and Sound Radiation
Three-dimensional (3D) finite element (FE) models of a standard metro wheel and a conventional slab track were developed to obtain the vibrations ( Figure 6). In the model, the normal forces and lateral creep forces obtained by the vehicle/track interaction dynamic model in Section 2.1 are applied to the tread of wheelset and rail head. Figure 6 also shows the FE models of resilient wheel and embedded track to analyze their effect on wheel/rail noise attenuation under lateral excitation. The FE model was developed in commercial software ANSYS.
where Mw1 and Mw2 are the mass of the rim and the rest parts of the wheelset as shown in Figure 5. , and fzL,Ri F are vertical forces due to the stiffness and damping of rubber, due to wheel/rail interaction, and due to primary suspension system.

The Transient Finite and Boundary Element Model for Wheel/Rail Vibration and Sound Radiation
Three-dimensional (3D) finite element (FE) models of a standard metro wheel and a conventional slab track were developed to obtain the vibrations ( Figure 6). In the model, the normal forces and lateral creep forces obtained by the vehicle/track interaction dynamic model in Section 2.1 are applied to the tread of wheelset and rail head. Figure 6 also shows the FE models of resilient wheel and embedded track to analyze their effect on wheel/rail noise attenuation under lateral excitation. The FE model was developed in commercial software ANSYS. Once the vibrations of the wheel and track are determined using by 3D FE model, sound radiations from the wheel and track can be evaluated using an appropriate vibro-acoustic prediction Once the vibrations of the wheel and track are determined using by 3D FE model, sound radiations from the wheel and track can be evaluated using an appropriate vibro-acoustic prediction method. There are three radiators: the wheels, the rails and the slabs. Equation (20) shows the 3D acoustic integral equation [36] for a sound field of frequency ω.
where p(x, y, z) denotes sound pressure amplitude in the acoustic domain, v n (x , y , z ) denotes air particle velocity amplitude in the normal direction of the domain boundary (the unit normal vector of the boundary, denoted by n, is pointing out of the domain, and finally, is the Green's function with the source at (x, y, z) and observer at (x , y , z ), and k is the wavenumber. Equation (20) is solved by 3D boundary element method (BEM). Discrete boundary surface into several sub boundary elements Ω ae . In this way, the distribution of sound pressure p and normal velocity v n on the boundary surface is approximated to the sound pressurep and normal velocitŷ v n in a series of sub boundary elements. Then, the sound pressure p and normal velocity v n can be written as: The N e i is the shape function of node i. Defined that N e i is one at the node i, and zero at other nodes. a pi and a vi are the sound pressure and normal velocity of the node i.
Defined that [N] is the (1 × n a ) matrix of shape function N e i . {p i } is the (n a × 1) matrix of sound pressure a pi . {v ni } is the (n a × 1) matrix of normal velocity a vi . Here n a is the total number of nodes. p(x , y , z ) andv n (x , y , z ) can be further written as: To use (24) and (25) solving Equation (20), the relationship ofp i andv ni needs to be determined. The relationship of them can be written as: The ith column element A bi in the coefficient matrix [A b ] can be expressed as: The ith column element B bi in the coefficient matrix [B b ] can be expressed as: where The ith column element C i in the coefficient matrix [C] can be expressed as: The ith column element D i in the coefficient matrix [D]. can be expressed as: For the BE model of wheel, the hole left in the wheel hub is closed to model the effect of axle. For the BE model of track, sound radiation only considers rail and slab surface. For a railway on the ground surface, the acoustic domain may be approximated to be a half-space. It has two acoustically hard boundaries (the sound absorption effect of the ground surface will not be considered in this paper). They are the horizontal plane beyond the top surface of the slabs. Other boundaries of the acoustic domain are radiating, including the top surface of the slabs and the surfaces of the rails. Connections between the rails and the slabs are ignored, leaving a uniform gap between the rails and the slabs. The effect of the presence of a train on the sound field will not be considered in the current paper. Thus, the acoustic domain is uniform in the track direction.

Wheel/Rail Noise Prediction and Preliminary Validation
This section gives the wheel/rail force, acceleration and sound pressure level of the wheels under the wheel/rail force excitation with and without squealing phenomenon in a sharp curve. When a vehicle curves, the leading wheel (front wheel) has a potential to produce squeal and the trailing wheel (rear wheel) should not produce squeal. Figure 7 gives the lateral creep force of the leading wheelset when a vehicle runs from a tangent line to a narrow circular curve at the speed of 45 km/h. The curve radius is 250 m, and the track cant is 96 mm. The main physical parameters of the metro vehicle used in the simulations are given in Table 1.  It can be seen from Figure 7 that the lateral creep force is very stable at tangent line. In addition, the fluctuation of lateral creep force increases gradually from transition curve to circular curve, and reaches maximum when the vehicle negotiates the circular curve. Figure 8 shows the lateral creep forces of leading and trailing wheels. It can be seen that the fluctuation of lateral creep force of leading wheel is larger than trailing wheel. Figure 9a,b give the acceleration of leading inner wheel and trailing inner wheel in time domain and frequency domain. It can be seen that the acceleration of the leading inner wheel with squeal is much larger than that of trailing inner wheel without squeal. In addition, highest peaks occur at modal frequencies corresponding to axial modes with zero-nodal-circle, showing that these modes are easy to be excited during curving. The axial mode shapes of wheel are expressed as (m, n). m represents number of nodal circle. n represents number of nodal diameter. For example, (0, 1) represents zero-nodal-circle with one-nodal-diameter. Thus, the frequencies and corresponding axle It can be seen from Figure 7 that the lateral creep force is very stable at tangent line. In addition, the fluctuation of lateral creep force increases gradually from transition curve to circular curve, and reaches maximum when the vehicle negotiates the circular curve. Figure 8 shows the lateral creep forces of leading and trailing wheels. It can be seen that the fluctuation of lateral creep force of leading wheel is larger than trailing wheel. , 7, 1144 11 of 20 Appl. Sci. 2017, 7, x; doi: www.mdpi.com/journal/applsci frequencies by simulation in this paper have a good agreement with the FRF test in the semi-anechoic room [8].   It can be seen that the acceleration of the leading inner wheel with squeal is much larger than that of trailing inner wheel without squeal. In addition, highest peaks occur at modal frequencies corresponding to axial modes with zero-nodal-circle, showing that these modes are easy to be excited during curving. The axial mode shapes of wheel are expressed as (m, n). m represents number of nodal circle. n represents number of nodal diameter. For example, (0, 1) represents zero-nodal-circle with one-nodal-diameter. Thus, the frequencies and corresponding axle mode shapes are 265 Hz (0, 1), 353 Hz (0, 0), 468 Hz (0, 2), 1167 Hz (0, 3), and 2008 Hz (0, 4). These frequencies by simulation in this paper have a good agreement with the FRF test in the semi-anechoic room [8]. Figure 10a,b give the sound pressure level (SPL) (at the point which is 30 cm away from the outside of wheel tyre) of leading wheel and trailing wheel in time domain and frequency domain. It can be seen that the SPL of the leading wheel with squeal is much higher than that of trailing wheel at frequencies of axial modes. , 7, 1144 frequencies by simulation in this paper have a good agreement with the FRF test in the semi-anechoic room [8].  Wheel/rail noise is also tested near a wheel (860 mm diameter) used in metro when the vehicle passes a sharp curve (radius = 250 m) at the speed of 45 km/h. The microphone is arranged in the outside of the leading wheel with the distance of 30 cm, which is shown in Figure 11. In the process of field test, squealing noise can be heard clearly. Figure 12 gives comparison between simulation and experimental results of SPL of the leading wheel with the same curve radius and running speed. Wheel/rail noise is also tested near a wheel (860 mm diameter) used in metro when the vehicle passes a sharp curve (radius = 250 m) at the speed of 45 km/h. The microphone is arranged in the outside of the leading wheel with the distance of 30 cm, which is shown in Figure 11. In the process of field test, squealing noise can be heard clearly. Figure 12 gives comparison between simulation and experimental results of SPL of the leading wheel with the same curve radius and running speed.
Appl. Sci. 2017, 7, 1144 13 of 20 Figure 11. Microphone arrangement.  Figure 12 shows that the SPL at the frequencies of 265-468 Hz, 1167 Hz, 2008 Hz are significant for both simulation and experimental results. These frequencies have a good agreement with the zero-nodal-circle and n-nodal-diameter axial wheel modes for simulation and experimental results. Amplitude difference percentages between simulation and experiment of 265-468 Hz, 1167 Hz, and 2008 Hz are 3%, 7%, and 6.5% of their corresponding peaks, respectively. The differences mainly come from the simplification and assumption of the actual situation from the simulation model. In other words, the current model can predict all squeal frequencies and give relatively accurate amplitudes of them. Thus, this simplified model can be used for regular investigation on wheel/rail noise when the vehicle curves.

Effect of Wheel/Rail Friction Coefficients on Wheel/Rail Noise under Lateral Excitation
This section investigates the effect of friction coefficient between wheel and rail on vibration and wheel/rail noise under lateral excitation. The friction coefficient of 0.012-0.3 are chosen for analysis because they cover the range of common friction coefficients. The value of 0.3 represents dry friction   Figure 12 shows that the SPL at the frequencies of 265-468 Hz, 1167 Hz, 2008 Hz are significant for both simulation and experimental results. These frequencies have a good agreement with the zero-nodal-circle and n-nodal-diameter axial wheel modes for simulation and experimental results. Amplitude difference percentages between simulation and experiment of 265-468 Hz, 1167 Hz, and 2008 Hz are 3%, 7%, and 6.5% of their corresponding peaks, respectively. The differences mainly come from the simplification and assumption of the actual situation from the simulation model. In other words, the current model can predict all squeal frequencies and give relatively accurate amplitudes of them. Thus, this simplified model can be used for regular investigation on wheel/rail noise when the vehicle curves.

Effect of Wheel/Rail Friction Coefficients on Wheel/Rail Noise under Lateral Excitation
This section investigates the effect of friction coefficient between wheel and rail on vibration and wheel/rail noise under lateral excitation. The friction coefficient of 0.012-0.3 are chosen for analysis  Figure 12 shows that the SPL at the frequencies of 265-468, 1167, 2008 Hz are significant for both simulation and experimental results. These frequencies have a good agreement with the zero-nodal-circle and n-nodal-diameter axial wheel modes for simulation and experimental results. Amplitude difference percentages between simulation and experiment of 265-468, 1167, and 2008 Hz are 3%, 7%, and 6.5% of their corresponding peaks, respectively. The differences mainly come from the simplification and assumption of the actual situation from the simulation model. In other words, the current model can predict all squeal frequencies and give relatively accurate amplitudes of them.
Thus, this simplified model can be used for regular investigation on wheel/rail noise when the vehicle curves.

Effect of Wheel/Rail Friction Coefficients on Wheel/Rail Noise under Lateral Excitation
This section investigates the effect of friction coefficient between wheel and rail on vibration and wheel/rail noise under lateral excitation. The friction coefficient of 0.012-0.3 are chosen for analysis because they cover the range of common friction coefficients. The value of 0.3 represents dry friction coefficient. The value of 0.012 represents the friction coefficient considering the third medium-oil between wheel/rail interface. The values of 0.1 and 0.05 represent friction coefficients between dry status and oil status. Figure 13a,b give the lateral accelerations of leading wheel with different friction coefficients in time domain and frequency domain. It can be seen that the acceleration decreases obviously at curve squeal frequencies as the friction coefficient increase. Figure 14a,b give the SPLs with different friction coefficients in time domain and frequency domain. It can be seen that SPLs also have a good reduction in significant peaks after reducing the friction coefficient. The sound reduction (SPL) of wheel considering lubricant can reduce 10-25 dB corresponding to the friction coefficient of 0.012-0.05, compared with the SPL of friction coefficient of 0.3. Therefore, it suggests that proper lubrication between the wheel and rail interface can reduce squeal noise effectively.
Appl. Sci. 2017, 7, 1144 14 of 20 coefficient. The value of 0.012 represents the friction coefficient considering the third medium-oil between wheel/rail interface. The values of 0.1 and 0.05 represent friction coefficients between dry status and oil status. Figure 13a,b give the lateral accelerations of leading wheel with different friction coefficients in time domain and frequency domain. It can be seen that the acceleration decreases obviously at curve squeal frequencies as the friction coefficient increase. Figure 14a,b give the SPLs with different friction coefficients in time domain and frequency domain. It can be seen that SPLs also have a good reduction in significant peaks after reducing the friction coefficient. The sound reduction (SPL) of wheel considering lubricant can reduce 10-25 dB corresponding to the friction coefficient of 0.012-0.05, compared with the SPL of friction coefficient of 0.3. Therefore, it suggests that proper lubrication between the wheel and rail interface can reduce squeal noise effectively.

Effect of Embedded Track on Wheel/Rail Noise under Lateral Excitation
This section investigates the effect of embedded track on wheel/rail noise under lateral excitation. The SPL of embedded track is calculated based on the embedded track dynamic model and transit BEM model in Sections 2.1 and 2.2. Figure 15 shows the SPLs of track comparison between conventional slab track and embedded track. It can be seen that the SPL of embedded track without considering wheel has a very good reduction in broad bands due to its damping of the elastic material around the rail absorbs vibration and then resulting in a low noise. Figure 16 gives the total SPLs considering standard metro wheels coupled with the conventional slab track and embedded track respectively. It can be seen that when the embedded track coupled with a squealing wheel, the reduction of sound is limited (about 2-3 dB) because the wheel is the main sound radiation source when squeal occurs, as is shown in Figure 17. In Figure 7, the SPL of a standard metro wheel is compared with the embedded track because the contribution of wheel/rail noise must be clear. The SPL of the standard metro wheel is larger than that of track in a wide frequency domain which contributes significantly to the total value of wheel-rail noise. Although the

Effect of Embedded Track on Wheel/Rail Noise under Lateral Excitation
This section investigates the effect of embedded track on wheel/rail noise under lateral excitation. The SPL of embedded track is calculated based on the embedded track dynamic model and transit BEM model in Sections 2.1 and 2.2. Figure 15 shows the SPLs of track comparison between conventional slab track and embedded track. It can be seen that the SPL of embedded track without considering wheel has a very good reduction in broad bands due to its damping of the elastic material around the rail absorbs vibration and then resulting in a low noise. Figure 16 gives the total SPLs considering standard metro wheels coupled with the conventional slab track and embedded track respectively. It can be seen that when the embedded track coupled with a squealing wheel, the reduction of sound is limited (about 2-3 dB) because the wheel is the main sound radiation source when squeal occurs, as is shown in Figure 17. In Figure 7, the SPL of a standard metro wheel is compared with the embedded track because the contribution of wheel/rail noise must be clear. The SPL of the standard metro wheel is larger than that of track in a wide frequency domain which contributes significantly to the total value of wheel-rail noise. Although the embedded track is not significant on the reduction of wheel/rail noise in the curve, it may have more effect on reducing rolling noise [37,38]. , 7, 1144 16 of 20 embedded track is not significant on the reduction of wheel/rail noise in the curve, it may have more effect on reducing rolling noise [37,38].   Appl. Sci. 2017, 7, 1144 16 of 20 embedded track is not significant on the reduction of wheel/rail noise in the curve, it may have more effect on reducing rolling noise [37,38].

Effect of Resilient Wheel on Wheel/Rail Noise under Lateral Excitation
This section investigates the effect of resilient wheel on wheel/rail noise under lateral excitation. The SPL of resilient wheel is calculated based on the resilient wheel dynamic model and transit BEM model in Sections 2.1 and 2.2. Figure 18 shows the SPLs of wheel comparison between the standard wheel and the resilient wheel. It can be seen that the resilient wheel has a very good reduction in sharp peaks, especially for the peaks of squeal frequencies, due to high damping in the resilient wheel and decoupling of the rim and web.

Effect of Resilient Wheel on Wheel/Rail Noise under Lateral Excitation
This section investigates the effect of resilient wheel on wheel/rail noise under lateral excitation. The SPL of resilient wheel is calculated based on the resilient wheel dynamic model and transit BEM model in Sections 2.1 and 2.2. Figure 18 shows the SPLs of wheel comparison between the standard wheel and the resilient wheel. It can be seen that the resilient wheel has a very good reduction in sharp peaks, especially for the peaks of squeal frequencies, due to high damping in the resilient wheel and decoupling of the rim and web.  Figure 19 gives the total SPLs of conventional slab track coupled with the standard metro wheel and resilient wheel respectively. It can be seen that when the resilient wheel coupled with the conventional slab track, the reduction of sound is significant (about 5-6 dB) because the resilient wheel can reduce squeal noise of wheel, which is main source for squeal noise, namely the main sound source is controlled. Thus, the total SPLs is reduced by using resilient wheel.   Figure 19 gives the total SPLs of conventional slab track coupled with the standard metro wheel and resilient wheel respectively. It can be seen that when the resilient wheel coupled with the conventional slab track, the reduction of sound is significant (about 5-6 dB) because the resilient wheel can reduce squeal noise of wheel, which is main source for squeal noise, namely the main sound source is controlled. Thus, the total SPLs is reduced by using resilient wheel. Figure 19 gives the total SPLs of conventional slab track coupled with the standard metro wheel and resilient wheel respectively. It can be seen that when the resilient wheel coupled with the conventional slab track, the reduction of sound is significant (about 5-6 dB) because the resilient wheel can reduce squeal noise of wheel, which is main source for squeal noise, namely the main sound source is controlled. Thus, the total SPLs is reduced by using resilient wheel.

Conclusions
A vehicle-track coupling dynamic model considering the falling wheel/rail friction coefficients was built to study the wheel/rail contact in narrow curve. In addition, a transient finite element and boundary element time domain models of vibration and sound were developed for sound radiation evaluation. Wheel/track vibration and sound pressure response are presented for a vehicle negotiating a narrow curve. The following conclusions can be drawn: (1) Different wheels have different levels of wheel/rail noise when the vehicle curves. The wheel/rail noise of the leading wheel is 20 dB larger than that of the trailing wheel. (2) Proper lubrication in wheel/rail interface can effectively reduce wheel/rail noise when the vehicle curves. The sound reduction of wheel considering lubricant (friction coefficient: 0.012~0.1) can reduce 10~25 dB compared to that of friction coefficient 0.3. (3) The application of vibration and noise control measures on wheels is better than that on track when the vehicle curves. The total sound reduction of the embedded track coupled with a squealing wheel is about only 2-3 dB. However, the total sound reduction of the elastic wheel coupled with the conventional slab track is about 5-6 dB.