Input parameter for rail car system modeling.
Abstract
The PID classic control systems are often employed for rail car systems to reduce the vibrations and disturbance rate during movement. In this study, the dynamic modeling and simulation of PID controls for rail car systems were carried out. Using 9 degrees of freedom, the modeling process comprises the representation of the rail car system and the rail track followed by the generation of equations of motion as well as differential equations for the rail car body, wheel sets and bogie. The represented systems are simulated in the MATLAB Simulink 2018 environment based on the equations of motion generated, and subsequently vibration analysis was carried out. The PID control system tuned according to the NicholsZiegler rule was introduced to minimize the vibrations and disturbance rate. The performance of the control and the rail car system in terms of the input step response, bandwidth, frequency, phase margin, frequency and input and output rejections was evaluated. The control system demonstrated significant robustness in providing the required active control for the system, while there was improved stability and reduction in noise and vibration under control action of the PID, thus improving ride comfort.
Keywords
 advanced controls
 modeling
 rail car
 simulation
 vibration
1. Introduction
The rail car system uses steel wheels moving on steel rails, and its suspension system comprises the bogie, primary and secondary suspension elements as well as springs and dampers. The rail car has the power bogies in the front and rear positions each having four wheels arranged in pairs and rigidly connected through an axle. The sets of wheel, which rotates at the same speed, are connected to the bogie through the primary suspension system, which is harder and stiffer in order to minimize load disturbances and uneven weight distributions. This is to maintain a good balance of the rail car while moving along its track. The stiff primary suspension system connects the wheel set to the bogie, while the soft secondary suspension system connects the bogie to the rail car body in order to isolate it from disturbances that stem from rail irregularities, uneven track profiles and its associated vibrations. While the primary suspension system is designed to provide guidance and rail car stability amidst load and weight variations, the secondary suspension system is to enhance comfortable ride by isolating rail car body from rail track irregularities. These design requirements of both suspension systems are to dampen the effect of vibration and increase the overall system’s performance [1, 2]. For a rail car, when vibrations are not kept within the permissible limits, the comfort, safety and health of passengers are at risk, and the time interval for maintenance activities for the components will likely increase. Hence, undesirable vibration in the rail car system is the cause of noise, considerable energy loss, reduction in the system’s performance, fatigue or fracture of some component, instability of a moving rail car and displacement of the rail track, amongst others. Over the years, there are three main types of suspension systems often employed to check vibration in the rail car system, namely, the passive, semiactive and fully active suspension systems [3, 4]. The design requirements of the suspension system of the rail car are to provide support against the dynamic loads and weight of the railcar; prevent load and rail disturbances; isolate the rail car body disturbances that could offset its balance; prevent irregular motions such as bouncing, yawing, etc.; provide guidance along the rail track; and optimize the curving, braking and other car maneuvering performances. The passive suspension system, which consists of springs, absorbs and stores the energy absorbed, while the dampers mounted on each wheel act as the shock absorber that dissipates the energy stored in the spring and reduces the vibrations from the rail transmitted to the vehicle ([5, 6, 7]. The passive suspension system is costeffective and driven by simple technology, but its demerit lies in the fact that it cannot measure some critical parameters such as the velocity, displacement, acceleration, etc. of the rail car system in real time in order to effect the needed adjustments to stabilize the rail car system; hence, it is rigid and cannot reach the compromise between alternate hard primary and soft secondary suspension systems amidst load and rail irregularities. On the other hand, the active suspension system uses sensors, comparators and actuator for measuring and monitoring some critical parameters that influence rail car stability in real time [8, 9]. The measured system’s parameters are compared with the threshold already preset on the controller, which can guarantee comfortable ride. The steadystate errors generated are eliminated through adjustment and compensation for such errors, and as such, there is highlevel damping without compromising the system’s performance. The third category of suspension system is the semiactive suspension system, which employs a spring and controllable damper. The spring element stores the energy, while the controllable damper dissipates the energy stored. The semiactive suspension systems combine the features of the passive and active suspension systems such as the use of passive damper and an actively controlled spring. The merit of the semiactive suspension is that it is cost and energy effective [9, 10, 11].
Over the years, researchers have employed several approaches ranging from classic to advance systems to control the suspension system in order to minimize the effect of vibrations. Such approaches include the use of proportionalintegralderivative (PID) control, Fuzzy PID, Linear Quadratic Regulator vibration controller, Adaptive NeuroFuzzy Inference System control and magnetorheological dampers, amongst others [12, 13, 14, 15, 16].
2. Materials and method
The analysis of vibration and control in the rail car system starts with the schematic representation of the rail car and track system including their degrees of freedoms and subsequent generation of equations of motion. The modeling is done with the masses of the system (rail car, body and wheel sets) having 6 degrees of freedom in the longitudinal, lateral and vertical directions as well as in the roll, pitch and yaw rotational motions. Figures 1 and 2 illustrate the free body diagram for the rail car and track system in their degrees of freedom.
In order to give a complete vibration analysis and control, the model design and simulation of the suspension systems are considered as both linear and nonlinear systems. There are three forces acting on the rail car suspension system, namely, the spring force, rolling resistance and forces due to the wheeltrack interactions. The spring force and the rolling resistance act on the rail car body with mass
Eq. (3) expresses the summation of the forces
The spring deflection is expressed as Eq. (7):
where
The transfer function was used for representing the linear systems, and the inputs are the load changes, applied forces as well as the uneven track profiles, while the output of the system is the acceleration and displacement of the rail car body as well as the deflection of the suspension systems.
According to Sezer and Atalay [19], the vectors for the displacement, rail car disturbance and control forces are expressed as Eqs. (12)–(14), respectively:
where.
The mathematical model can be obtained in either the state space or using the transfer function for necessary control actions. The control action is to improve the speed of response, system’s balance and stability and to reduce the steadystate error as well as the amplitude of oscillations. Eqs. (15) and (16) express the form of the secondorder systems:
where.
The performance of the control system is measured by the settling time, delay time rise time, percent overshoot and the steadystate error.
The settling time
On the other hand, the rise time is the time it takes the signal to change from a low value to a high value (say 10–90% or 0–100%). The time it takes the peak value to occur known as the time is expressed as Eq. (23):
The steadystate error
Similarly, the delay time is the time required for the response to reach half the final value for the first time, while the percent overshoot is the percent by which the step response of the system exceeds the final steadystate value. It is a parameter that defines the instability of a system (Eq. (25)):
Eq. (26) relates the percent overshoot to the damping ratio:
The input parameters of the system are presented in Table 1.
2.1 The proportionalintegralderivative (PID) control
The PID control represents the proportional, integral and derivative controls. It is a form of classic control designed to automatically reduce rise and settling times, steadystate errors and percent overshoot. The block diagram of the PID controller is illustrated in Figure 4. The threshold value otherwise referred to as the reference point is preset on the controller while realtime measurement using sensors is taken. The error generated, which is the difference between the threshold value and the actual measurement, represents the deviation from the ideal process; thereafter, the actuator effects realtime control to adjust process variables. The output signal of the PID controllers often responds to changes over time with respect to the actual measurements and set points. The PID variables are iteratively adjusted until the steadystate error from the output signal is eliminated. This control action is done by adjusting the controller gain with resulting decrease in the rise time and percent increase in the overshoot, which makes the system go unstable. This rise time is further reduced with the integral control action. Finally, the derivative action is introduced to compensate for the offset. This reduces the percent overshoot and settling time, thus making the system stable over time.
Eq. (27) gives the expression for the control action of the PID controller:
The NicholsZiegler tuning rules employed for tuning the PID control as well as the summary of the effects of its control action on the PID are presented in Tables 2 and 3, respectively.
S/N  Parameter  Notation  Value  Unit 

1  Average mass of the rail car 

50,500  kg 
2  Average mass of bogie 

2410  kg 
3  Mass of primary suspension system 

30,000  kg 
4  Mass of secondary suspension system 

30,000  kg 
5  Moments of inertia 

56,900  kgm^{2} 
6  Rail car roll inertia 

68,200  kgm^{2} 
7  Rail car pitch inertia 

71,000  kgm^{2} 
8  Average mass of first wheel set and axle 

1300  kg 
9  Average mass of second wheel 

1300  kg 
10  Distance between the centre of gravity and the front position of the rail car 

6  m 
11  Distance between the centre of gravity and the middle position of the rail car 

6  m 
12  Distance between the centre of gravity and the rear position of the rail car 

6  m 
13  Spring constant of the primary suspension system 

2.4 × 10^{6}  N/m 
14  Spring constant of the secondary suspension system 

5.6 × 10^{5}  N/m 
15  Spring constant of the wheel 

4.0 × 10^{5}  N/m 
16  Damping constant of the primary suspension system 

1.2 × 10^{3}  Ns/m 
17  Damping constant of the secondary suspension system 

2.95 × 10^{4}  Ns/m 
18  Damping constant of the wheel 

5.0 × 10^{4}  Ns/m 
S/N  Type of controller  K_{p}  K_{i}  K_{d} 

1  P  0.5 Kcr 

0 
2  PI  0.45 Kcr  0.83 Pcr  0 
3  PID  0.6 Kcr  0.5 Pcr  0.125 Pcr 
S/N  Controller response  Rise time  Overshoot  Setting time 

1  K_{p}  Decrease  Increase  Small change 
2  K_{i}  Decrease  Increase  Increase 
3  K_{d}  Small change  Decrease  Decrease 
The signal (U) which passes through the controller computes the derivative and integral of error signal. The signal error is thereafter sent to the system in order to obtain the system’s output (Y). The PID controller was designed in the MATLAB Simulink 2018 environment to generate a continuous time control. Using the NicholsZiegler rules, the tuning of PID controller was done by generating the system’s transfer function and subsequent importation of the parameters obtained into the linear timeinvariant system.
The PID control system and its connection to the rail car system are shown in Figures 5 and 6, respectively.
The aim of the control design is to keep the system variables close to the reference in order to compensate for the effect of load and rail disturbances. The system requirement is to check unpleasant motion and ensure rail car stability by reaching a compromise between the stiff primary suspension and soft secondary suspension system.
The actively controlled suspension system can be activated via the use of solenoid, hydraulic, electromagnetic means or through a magnetorheological damper. This system is designed to use the solenoid actuators because of its lightweight, simplicity in structure, ease of installation and short response time, which makes it highly sensitive to disturbances.
3. Results and discussion
Figure 7 shows the step response before the iterative adjustment of the PID control. The amplitude of oscillation, which is a function of the percent overshoot, is 2 mm, and the system could not return to the equilibrium position after 3 s. The shape of the plot represents a system that is underdamped, which signifies the need for damping to minimize unwanted motion. The system whose step response is depicted in Figure 7 is relatively unstable as vibration will reduce the system’s and ride performance.
Figure 8 shows the step response from the controlled system. When compared to Figure 1, the amplitude of oscillation has reduced to 1.15 mm and settling time 0.5 s under the effect of the PID control action. The system is relatively stable as the plot gives an indication of a critically damped system. This implies that under the effect of the PID control system, the system could regain its stability and equilibrium position after encountering some level of disturbances.
Figure 9 shows the response of the PID controller in order to determine whether the control system meets the design requirements. The tuned response represents the actual measurement of the process variables, while the baseline response represents the reference or set point. The differences between the tuned response and the baseline response give the steadystate error. The nature of the plot also indicates that the system is critically damped; thus, the damping ratio
Figure 9 shows the reference tracking of the tuned response and the reference (baseline). The design requirement of the control system includes setpoint tracking; hence, this plot shows how the closedloop system responds to a step change in set point. From the plot, the steadystate error is minimal, thus indicating the effectiveness of the control system.
Figure 10 shows the closedloop step response to a step disturbance at the system’s output. This is important in analyzing the sensitivity of the control system to noise measurement. The negligible value of the steadystate error as shown by the degree of agreement between the tuned output response of the system and the baseline (reference) indicates the high sensitivity and compensation of the system’s control to the noise measurement.
Figure 11 shows the closedloop system response to load disturbance. The plot represents a step disturbance from the system’s input. The differences between the tuned response and the baseline response, which is a function of the steadystate error, are significantly large. The implication of this is that the control system is insensitive to the rejection of load and other input disturbances, which is capable of offsetting the balance of the system. The large steadystate error resulting from input disturbance rejection stems from the fact that a single PID may not be able to satisfy all the design requirements at the same time; hence, there is always a performance tradeoff amongst the reference tracking, percent overshoot and input disturbance rejection. However, the use of Fuzzy PID or ISAPID controller can be used to meet the design requirements significantly. This will improve the response of the reference tracking with the provision of an additional tuning parameter, which allows for independent control of the effect of the reference signal on the proportional action.
The robustness and performance of the PID control are directly proportional to the degree of stability of the rail car, ride comfort and performance of the rail car. The summary of the performance of the PID control system is presented in Table 4.
Parameter  Tuned response  Baseline (reference) response 

Rise time  0.07 s  0.0847 s 
Settling time  0.499  1.94 s 
Overshoot  4.76%  10.3% 
Peak  1.05  1.1 mm 
Phase margin  83^{o}  60.4^{o} 
Frequency  16.6  17.8 rad 
Closed loop  Stable  Stable 
From Table 4, the tuned response has better performance and robustness compared to the baseline due to periodic iterative adjustment to eliminate the steadystate error for each step input. The rise time of the tuned response indicates a fast response time of the control system to disturbances or changes when compared to the baseline. This signifies some degree of delay in the time it takes the system to respond to fluctuations. In addition, comparing the settling time for both responses, the tune response settles faster after some disturbances compared to the baseline response. This explains the ability of the rail car to regain its stability after encountering some level of disturbances with the active control system. Also, the percent overshoot for the tuned response is still within the range of the permissible oscillation (5%) for critically damped system which indicates that the rail car system is relatively stable amidst load and rail disturbances. An increase in the phase margins implies a significant reduction in the percent overshoot and bandwidth.
Figure 12 shows the bode response of the control system. This is the plot of the frequency and phase response of the control system. The phase margin was found to be 83^{o} at a frequency of 16.6 rad/s for the tuned response compared to the baseline response with a phase margin and frequency of 60.4^{o} and 17.8 rad/s, respectively. The relationship between the phase margin and percent overshoot is inversely proportional. Hence, the high value of the phase margin results in significant decrease in the percent overshoot, thus improving the rail car stability. Also, the higher closedloop bandwidth results in faster rise time. The rise time was found to be 0.07 s for the tuned response and 0.0847 s for the baseline response. This implies that for the baseline response, the percent overshoot is still significant to offset the stability of the rail car system.
Figures 13 and 14 show the result of the linearization of the rail car system. This is to determine the dynamics of the system in real time and within time and frequency domains. Figure 13 shows the impulse response of the rail car system, which is the degree of rail car body displacement as a result of load or rail disturbances. The maximum amplitude oscillation is 0.02 mm, which is negligible and insufficient to offset the rail car balance. In addition, the settling time (less than 1 s) is also within the permissible range indicating the ability of the system to regain its stability within a second as a result of load or rail disturbances.
Figure 14 shows the bode plot of the rail car system; the phase margin was found to be 48^{o} at a frequency of 20.2 rad/s which signifies significant reduction in the percent overshoot due to the compensation for steadystate error by the derivative action of the PID.
4. Conclusion
The control of undesirable vibrations of the rail car system in 6 degrees of freedom was carried out using the PID control system. The initial displacement and vertical acceleration that characterize the performance of the system on encountering rail and load disturbances were minimized with the iterative adjustment of the PID controller according to NicholsZiegler tuning rules. In addition, the performance of the control and the rail car system in terms of the input step response, bandwidth, frequency, phase margin and input and output rejections was within the acceptable range. Hence, the PID control system shows significant robustness in providing the required active control for the system, while the rail car system shows improved stability and reduction in vibration under control action of the PID, thus improving ride comfort. However, a single PID may not sufficiently satisfy all the design requirements at the same time resulting in performance tradeoff. However, Fuzzy PID or ISAPID controller can be used to meet the design requirements significantly. This will further improve the performance and robustness of the rail car system.
References
 1.
AlZughaibi A, Davies H. Controller design for active suspension system of quarter car with unknown mass and timedelay. International Journal of Mechanical, Aerospace, Industrial, Mechatronic and Manufacturing Engineering. 2015; 9 (8):14841489  2.
Sharma RC, Sharma SK. Sensitivity analysis of threewheel vehicle’s suspension parameters influencing ride behaviour. Noise & Vibration Worldwide. 2018; 49 (7–8):272280  3.
Bideleh SMM, Mei TX, Berbyuk V. Robust control and actuator dynamics compensation for railway vehicles. Vehicle System Dynamics. 2016; 54 (12):17621784  4.
Yao JL, Shi WK. Development of a sliding mode controller for semiactive vehicle suspensions. Journal of Vibration and Control. 2013; 19 :11521160  5.
Arefsoliman MA. Adaptive LQR control strategy for active suspension system. 2011. SAE Technical Paper, 2011010430  6.
Herbst G. A simulative study on active disturbance rejection control (ADRC) as a control tool for practitioners. Electronics. 2013; 2 (3):246279  7.
Colombo EF, Gialleonardo ED, Facchinetti A, Bruni S. Active car body roll control in railway vehicles using hydraulic actuation. Control Engineering Practice. 2014; 31 :2434  8.
Gohrle C, Schindler A, Wagner A, Sawodny O. Design and vehicle implementation of preview active suspension controllers. IEEE Transactions on Control Systems Technology. 2014; 22 (3):11351142  9.
Jin XS, Xiao XB, Ling L, Zhou L, Xiong JY. Study on safety boundary for highspeed trains running in severe environments. International Journal of Rail Transportation. 2013; 1 (1–2):87108. DOI: 10.1080/23248378.2013.790138  10.
Nguyen SD, Choi SB, Nguyen QH. An optimal design of interval type2 fuzzy logic system with various experiments including magnetorheological fluid damper. IMechE Part C: Journal of Mechanical Engineering Science. 2014; 228 (17):30903106  11.
Nguyen SD, Nguyen QH. Design of active suspension controller for train cars based on sliding mode control, uncertainty observer and neurofuzzy system. Journal of Vibration and Control. 2017; 23 (8):13341353  12.
Kim HC, Choi SB, Lee GS. Performance analysis of a semiactive railway vehicle suspension featuring MR dampers. In: SPIE Smart Structures and Materials + Nondestructive Evaluation and Health Monitoring, San Diego, CA. Bellingham, WA: International Society for Optics and Photonics (SPIE); 2014. p. 905711  13.
Kumar PS, Sivakumar K, Kanagarajan R, Kuberan S. 2799 adaptive neuro fuzzy inference system control of active suspension system with actuator dynamics. Journal of Vibroengineering. 2018; 20 (1):541549  14.
MatamorosSanchez AZ, Goodall RM. Novel mechatronic solutions incorporating inerters for railway vehicle vertical secondary suspensions. Vehicle System Dynamics. 2015; 53 (2):13136  15.
Oh JS, Shin YJ, Koo HW, Kim HC, Park J, Choi SB. Vibration control of a semiactive railway vehicle suspension with magnetorheological dampers. Advances in Mechanical Engineering. 2016; 8 (4):113  16.
Suarez B, Felez J, Maroto J, Rodriguez P. Sensitivity analysis to assess the influence of the inertial properties of railway vehicle bodies on the vehicle’s dynamic behaviour. Vehicle System Dynamics. 2013; 51 (2):251279  17.
Zolotas AG, Goodall RM. Modelling and control of railway vehicle suspensions. In: Turner MC, Bates DG, editors. Mathematical Methods for Robust and Nonlinear Control, Lecture Notes in Control and Information Sciences. New York: Springer; 2007. pp. 373412  18.
Liang L, Xiao XB, Xiong JY, Zhou L, Wen ZF, Jin XS. A threedimensional model for coupling dynamics analysis of high speed traintrack system. Journal of Zhejiang University Science A (Applied Physics and Engineering). 2014:121  19.
Sezer S, Atalay AE. Application of fuzzy logic based control algorithms on a railway vehicle considering random track irregularities. Journal of Vibration and Control. 2011; 18 (8):11771198  20.
Ulum Z, Affaf M, Salmah JCE, Suparwanto A. Active suspension systems design of a light rail vehicle using MPC with preview information disturbance. In: 5th International Conference on Instrumentation, Control, and Automation (ICA) Yogyakarta, Indonesia, August 9–11, 2017. pp. 1823  21.
Sharma SK, Kumar A. Ride performance of a highspeed rail vehicle using controlled semi active suspension system. Smart Materials and Structures. 2017; 26 (1–19):055026