Microsoft Word - 22-2420_s Engineering, Technology & Applied Science Research Vol. 8, No. 6, 2018, 3603-3608 3603 www.etasr.com Pehlivan et al: Modeling and Validation of 2-DOF Rail Vehicle Model Based on Electro–Mechanical… Modeling and Validation of 2-DOF Rail Vehicle Model Based on Electro–Mechanical Analogy Theory Using Theoretical and Experimental Methods Fatih Pehlivan Faculty of Engineering Department of Mechanical Engineering Karabuk University Karabuk, Turkey fatihpehlivan@karabuk.edu.tr Cihan Mizrak Faculty of Engineering Department of Mechatronic Engineering Karabuk University Karabuk, Turkey cihanmizrak@karabuk.edu.tr Ismail Esen Faculty of Engineering Department of Mechanical Engineering Karabuk University Karabuk, Turkey iesen@karabuk.edu.tr Abstract—This paper presents theoretical and experimental results on modeling and simulation of two degrees of freedom rail vehicle by using electro-mechanical similarity theory. In this study, the equations of motion were derived using Newton’s second law of motion and then mechanical and equivalent electrical circuits were obtained with the help of a free body diagram. A schema in Simulink allowing analyzing of the behavior of the primary and secondary suspension was created. In order to verify the electrical model, transfer function and schema were developed in Simulink. The simulation results were compared with the experimental data and the comparison showed that the results of the mechanical experiments were close to the simulation results, but the electrical results showed better periodic behavior. Keywords-electro-mechanical analogy; equivalent electrical circuit; modeling; simulation; simulink; vehicle; suspension I. INTRODUCTION Urban and intercity transportation have been provided by highways in many countries and so issues regarding traffic congestion, traffic accidents and increasing transportation cost are emerging. In developed countries, railway transport systems are preferred for either urban or intercity transportation because of the fact that they are one of the safest and the most economical transport modes [1–3]. Recently the demand for railway transport systems has increased globally. The most crucial effect determining the driving safety and ride comfort is vehicle vibrations caused by dynamic forces occurring during movement. It was shown that during movement, vibrations occur due to dynamic forces usually caused by wheel-rail interaction [4–7]. This vertical oscillation of the rail vehicle can be minimized by optimizing spring stiffness and damping coefficient of the shock absorber as shown in [8, 9]. Models with various degrees of freedom (DOF) were examined to determine the dynamic behavior of the rail vehicles in [10]. Authors in [7] realized dynamic solutions of rail vehicles from single-DOF to multi-DOF and also worked on rail-wheel contact problems. Authors in [11] studied the dynamic behavior of two degrees of freedom (2-DOF) wagon models under random vibrations. Authors in [12] designed the vehicle suspension system of the quarter vehicle model in Matlab/Simulink and MSC-ADAMS and compared the results. The mathematical model of the system was created using the state space model and the vibration of the system was investigated by the unit step input function. Authors in [13] studied the mathematical modeling and simulation of a quarter- car with two DOF. The state space mathematical model of the system was derived using Newton's second law of motion and free-body diagram. Afterwards, the performance of the system was examined with Matlab/Simulink. Passive, semi-active and active suspension systems have been tested with unit step input function. Authors in [14] studied the effect of spring and damping coefficients on acceleration and deflection using Matlab/Simulink. Authors in [15] aimed to improve driving comfort and road handling. Authors in [16] modeled a wagon moving on a bridge considered to be an elastic beam and done some calculations in terms of vibration and acoustics. They also calculated the response of the harmonic force that comes from the bridge and the rail. Authors in [17] studied the responses of the contact force and displacement caused by wheel-rail interaction. Analysis of the 5-DOF quarter-rail model supported by a 3-layer rail system was examined by using the equations of motion of vehicle and rail. Most of the studies performed regarding rail vehicles have been related to keeping the vehicle vibrations at optimum level. In this paper, vehicle vibrations were studied, however, results were achieved by creating not only mechanical prototypes but also their equivalent electrical circuits. Experimental data obtained from the mechanical mechanism or the equivalent electrical circuit were compared with the simulation results. II. MODEL OF THE SYSTEM A. Mathematical Model of the 2-DOF Rail Vehicle In order to constitute the mathematical model of the rail vehicle, it is beneficial to understand the quarter, half and full car models which are widely used in the literature due to their similarity with the models of rail vehicles. The modeled 2-DOF Engineering, Technology & Applied Science Research Vol. 8, No. 6, 2018, 3603-3608 3604 www.etasr.com Pehlivan et al: Modeling and Validation of 2-DOF Rail Vehicle Model Based on Electro–Mechanical… rail vehicle shown in Figure 1(a), consists of two masses: m1 is the mass representing the rail vehicle body or wagon and m2 is the mass representing the bogie system. k1 and b1 are the stiffness and damping coefficients of the secondary suspension respectively, k2 and b2 are the stiffness and damping coefficients of the primary suspension respectively, x1 and x2 are the vertical displacements of the rail vehicle body and the bogie system respectively and y is rail excitation as shown in [18–20]. D’Alembert’s principle was used to draw the free- body diagram of this system shown in Figure 1(b). (a) (b) Fig. 1. (a) The 2-DOF rail vehicle model. (b) Vehicle free body diagram. Accordingly, the equations for the vertical motion of the rail vehicle body and the bogie system are obtained by (1) and (2): ����� � ���� ��� � � � �� � ����� 0 (1) ���� � ���� ��� � � �������� � ��� � ��� � � � �� � ����� 0 (2) B. Mechanical and Equivalent Electrical Circuit Substituting values in (3), (4) and (5), the force representation for the 2-DOF rail vehicle can be written as: ��� � ��� � ��� 0 (3) ��� � ��� � ��� � ��� � ��� 0 (4) �� � ��� � ��� 0 (5) The mechanical circuit of the system is obtained using (3), (4) and (5) as shown in Figure 2. Fig. 2. The mechanical circuit of the modeled rail vehicle with 2 DOF. Mechanical systems have one excitation source (force), while electrical systems have two excitation sources (voltage difference and current). In this sense, mechanical systems have two different electrical counterparts. The first similarity theory developed to find the first electrical charge is obtained by using the force-voltage difference similarity and the Kirchoff voltage law (KVL). To find the other electrical counterpart, the second theory of similarity is established, the force-flow similarity and the Kirchoff current law (KCL) are used [21]. The electrical responses of the mechanical elements in both theories are obtained by (6)-(8) which are generated using Kirchoff's laws and are shown in Table I: ���� �� ��� � ���� � � � ������ (6) ���� � ������� � !��� � � " � !����� (7) !��� # �$����� � $��� % � � & � ������ (8) TABLE I. ANALOGOUS QUANTITIES Mechanical Quantity Electrical Quantity (1st similarity theory) Electrical Quantity (2nd similarity theory) Force (F) Voltage (V) Current (I) Velocity (v) Current (I) Voltage (V) Mass (m) Inductance (L) Capacitance (C) Friction (b) Resistance (R) Inverse Resistance (1/R) Spring constant (k) Inverse Capacitance (1/C) Inverse Inductance (1/L) Here, the second theory is selected to be applied to the modeled 2-DOF rail vehicle. According to the velocity-voltage similarity or the second similarity theory, the connections in parallel or serial of the mechanical circuit will be the same for the electrical system and so the electrical equivalent circuit of the rail vehicle is obtained by using the mechanical circuit of the system shown in Figure 2. The obtained circuit illustrated in Figure 3, and consists of two capacitances: C1 is mass representing rail vehicle body or wagon and C2 is mass representing the bogie system. L1 and R1 are the inductance and resistance stiffness and damping coefficients of secondary suspension respectively, L2 and R2 are the stiffness and damping coefficients of primary suspension respectively, and V is the excitation voltage. Fig. 3. The mechanical circuit of the modeled 2-DOF rail vehicle. C. Theoretical Methods In order to check the theoretical implementations of the system, a crank-connecting rod mechanism is considered for a harmonic rail excitation in the design of the mechanical system because it is simple and the most fundamental mechanism to transform the rotary motion to translational motion or vice versa [22]. The trigonometric method is used for the kinematic analysis and the velocity of the wheel is calculated for the rail excitation as shown in the following equations: � cos * �+ cos , (9) Engineering, Technology & Applied Science Research Vol. 8, No. 6, 2018, 3603-3608 3605 www.etasr.com Pehlivan et al: Modeling and Validation of 2-DOF Rail Vehicle Model Based on Electro–Mechanical… cos , -1 � sin� , (10) 123 4 % 123 5 6 (11) cos , 71 � 8%6 9 � sin� * (12) � cos * �+71 � 8%6 9 � sin� * (13) � cos * �7+� � +� 8%6 9 � sin� * (14) � � * sin * � ��� :-+� � ; � � sin� *< (15) � � * sin * � :=5 >�%� 123 �5<�-?�=>�%� 123� 5 (16) The 2-DOF rail vehicle consists of three main parts representing rail actuator, rail vehicle body and bogie system. They are separately designed and assembled together with other parts such as suspension, step motor, linear bearing etc. Figure 4 shows the assembly view of the designed parts of the objective rail vehicle. Fig. 4. The designed 2-DOF rail vehicle. D. Transfer Function of the Designed System Taking the Laplace transformation of (1) and (2), either the connection between the rail vehicle body and bogie system or the connection between the bogie system and the rail actuator are found as in the following equations. Equation (23) gives the transfer function of the 2-DOF rail vehicle: ��@��A� � ��@��A� � �A@��A� � �A@��A� ��A�@��A� (17) @��A� B��C �D��CD�� ��CD�� E @��A� (18) ��FG�A� � @��A�H � �AFG�A� � @��A�H��F@��A� � @��A�H � �AF@� �A� � @��A�H ��A�@��A� (19) � �A � ���G�A� � � �A � ���@��A� F��A� �� � � ��A � ��� � ���H@� �A� (20) � �A � ���G�A� � � �A � ���@��A� F��A� � � � � ��A � ��� � ���H B��C �D��CD�� ��CD�� E @��A� (21) J��C� K�C� ����C�D����� D�����CD�� �� ������CLD�M�CND�O�C�D�"�CD���� (22) where: P �� � � �� � � �� �, Q ���� � � � � ���� � ����, and # ��� � ��� R�A� J��C�K�C� J� �C� K �C� S (23) E. Mathematical Model of the System in MATLAB Simulink Simulink offers a graphical solution for modeling, simulating and analyzing dynamic systems over customizable blocks and tool library [23]. A Simulink model of the system is created using the mathematical model obtained through (1) and (2) as shown in Figure 6. F. Result of the Theoretical Methods in MATLAB Simulink In order to verify the electrical model, a transfer function and a model were developed in Simulink as shown in Figure 5. The results were obtained by using (16) for the rail excitation and the system parameters shown in Table II. Figure 7 shows the vertical displacement of the rail vehicle body with time and it is seen that the transfer function and the Simulink model of the mechanical system have one-to-one correspondence with the electrical model of the system. III. EXPERIMENTAL STUDY Experiments were conducted to validate both mechanical and electrical approaches. The 2-DOF rail vehicle model prototype was created. The theoretical study and the parameters shown in Table II were used in simulations to compare the experimental results with the theoretical ones. Fig. 5. Matlab representation of the electrical model, the transfer function and the Simulink model of the system Engineering, Technology & Applied Science Research Vol. 8, No. 6, 2018, 3603-3608 3606 www.etasr.com Pehlivan et al: Modeling and Validation of 2-DOF Rail Vehicle Model Based on Electro–Mechanical… Fig. 6. Mathematical model of the system in Simulink. TABLE II. PARAMETER VALUES USED IN THE 2-DOF RAIL VEHICLE Parameters Values Mass of the rail vehicle body or wagon, m1 0,794kg Damping coefficients of secondary suspension, b1 57,53N.s/m Stiffness of secondary suspension, k1 400N/m Mass of the bogie system, m2 0,748kg Damping coefficients of primary suspension, b2 57,53N.s/m Stiffness of primary suspension, k2 400N/m Crank length, R 0,025m Connecting rod length, l 0,084m Engine speed, * 74rpm Fig. 7. Verification of the electrical model of the system by using the transfer function and the Simulink model. A. Mechanical Test Rig The design of the experimental prototype, which comprises the primary and secondary suspension mechanisms and the excitation mechanism, is shown in Figure 8. Aluminum shock absorbers, stepper motor and its drive module, control card, chromium rods, linear bearings, rod holders, power supply and vibration datalogger were used in the design. The parts representing the rail excitation, the bogie and the wagon were made up of particle board. It was preferred because it is lighter than all other metals, including aluminum and is more rigid than composite or plastic. Mechanical vibration experiments were performed on the prototype rail vehicle test rig shown in Figure 9. In the mechanical test rig, a stepper motor controls the rail excitation by using a crank-connecting rod mechanism and the bottom particle board which is below the primary suspension and provides the excitation as velocity to the system. In this way, both primary and secondary suspensions start to move and acceleration of the wagon which is the top particle board is taken with vibration datalogger mounted to the wagon. (a) (b) Fig. 8. (a) Crank-connecting rod mechanism, (b) kinematic analysis. Fig. 9. Mechanical test rig Engineering, Technology & Applied Science Research Vol. 8, No. 6, 2018, 3603-3608 3607 www.etasr.com Pehlivan et al: Modeling and Validation of 2-DOF Rail Vehicle Model Based on Electro–Mechanical… B. Electrical Test Rig To verify the theoretical study and the designed and manufactured prototype of the 2-DOF mechanical rail vehicle model, it is necessary to establish the equivalent electrical test rig of the system. For this purpose, the first step is to create voltage representing the rail excitation. An arbitrary voltage generator can supply the existing excitation because of the fact that it can produce periodic signals as the rail excitation in the mechanical test rig. An AA Tech AWG-1020 function generator was used in the experiments. The current excitation signal was generated as shown in Figure 10 via Easywave software and was sent to the voltage generator. Fig. 10. The electrical excitation signal For convenience, the electrical parameters were selected such that the ratio of electrical parameters and their representing mechanical parameters be 1:100, which means 1kg corresponds to 0.01F etc. The parameter values are shown in Table III. Except for the inductance element, all values shown in Table III were obtained with elements available in the market. In this study, two inductance elements having values of 250mH were needed. An inductor of this value is difficult to find on the market and it is obtained by wrapping copper wire on ferrite [24]. The number of turns is calculated by using (24): S 7 &MT (24) where L is the desired inductance value and AL is the magnetic conductivity. The magnetic conductivity of the ferrite used to wrap the copper wire is 24700x10 -9 H. The electrical experiments were performed on the created electrical circuit shown in Figure 11. Using the GW INSTEK GDS-1102-U digital oscilloscope, the voltage values on the capacitor (C1), which is the electrical equivalent of the wagon, were measured and transferred to the Freewave software and experimental data were obtained. TABLE III. PARAMETER VALUES USED IN THE 2-DOF RAIL VEHICLE Equivalent Parameter Used Values Wagon mass, C1 Capacitor 0,008F Damping coefficients of secondary suspension, R1 Resistor 1,7424Ω Stiffness of secondary suspension, L1 Inductor 0,25H Mass of the bogie system, C2 Capacitor 0,0075F Damping coefficients of primary suspension, R2 Resistor 1,7424Ω Stiffness of primary suspension, L2 Inductor 0,25H Fig. 11. Electrical equivalent circuit set up IV. RESULTS To compare the experimental results with the theoretical mechanical and electrical approach, a series of simulations were conducted in Matlab. The vertical vibrations obtained from the mechanical experiment (Figure 12) and the vertical vibrations obtained from the electrical experiment (Figure 13) were compared with the simulation results. Fig. 12. Mechanical experiment and theoretical results comparison of vibration values of the rail vehicle model (a) first and (b) second experiment. Fig. 13. Electrical experiment and theoretical results comparison of vibration values of the rail vehicle model (a) first and (b) second experiment. Engineering, Technology & Applied Science Research Vol. 8, No. 6, 2018, 3603-3608 3608 www.etasr.com Pehlivan et al: Modeling and Validation of 2-DOF Rail Vehicle Model Based on Electro–Mechanical… Both experimental methods were applied two times and the sampling duration was 6s in each experiment. When the root- mean-square (RMS) values of the vertical vibrations obtained from the mechanical test rig, the equivalent electrical test rig and the theoretical model were compared to each other, it was seen that the mechanical test rig and the theoretical methods had the same results by a maximum error of 9.23% as shown in Table IV. On the other hand, the maximum error was 25.71% when the equivalent electrical test rig was compared with the theoretical methods as shown in Table V. TABLE IV. RMS-THEORITICALVALUE RESULT COMPARISON - MECHANICAL TEST RIG RMS Values Error Mechanical test rig (1st experiment) 1,084 6,47% Mechanical test rig (2nd experiment) 1,266 9,23% Theoretical method 1,159 TABLE V. RMS-THEORITICAL VALUE COMPARISON-EQUIVALENT ELECTRICAL TEST RIG RMS Values Error Equivalent electrical test rig (1st experiment) 0,11717 25,71% Equivalent electrical test rig (2nd experiment) 0,12749 19,17% Theoretical method 0,15773 V. CONCLUSION Two-DOF rail vehicle has been studied and modeled by mechanical test rig, its equivalent electrical circuit and theoretical methods. Both the mechanical test rig and its equivalent electrical circuit were validated by simulations and both methods were found to be capable of determining the vertical vibrations of both bogie and wagon. The mechanical test rig RMS results are closer to the simulation results in comparison with the electrical circuit, but the electrical results show better periodic behavior. REFERENCES [1] S. Bhattacharjee, A. R. Goetz, “Impact of light rail on traffic congestion in Denver”, Journal Of Transport Geography, Vol. 22, pp. 262-270, 2012 [2] R. Lalive, S. Luechinger, A. Schmutzler, Does Supporting Passenger Railways Reduce Road Traffic Externalities?, University of Zurich, Working Paper Series, No. 110, University of Zurich, 2013 [3] M. Metin, Rayli Sistem Araclarının Modellenmesi ve Titresimlerinin Kontrolu, Yildiz Technical University, 2007 (in Turkish) [4] A. Wijeweera, H. To, M. B. Charles, K. Sloan, “A time series analysis of passenger rail demand in major Australian cities”, Economic Analysis And Policy, Vol. 44, No. 3, pp. 301-309, 2014 [5] E. Matsika, S. Ricci, P. Mortimer, N. Georgiev, C. O’Neill, “Rail vehicles, environment, safety and security”, Research In Transportation Economics, Vol. 41, No. 1, pp. 43-58, 2013 [6] M. Graa, M. Nejlaoui, A. Houidi, Z. Affi, L. Romdhane, “Modeling and Simulation for Vertical Rail Vehicle Dynamic Vibration with Comfort Evaluation”, in: Multiphysics Modelling and Simulation for Systems Design and Monitoring, Springer, 2015 [7] R. V. Dukkipati, V. K. Garg, Dynamics of Railway Vehicle Systems, Academic Press, 1984 [8] A. Orvnas, “Methods for Reducing Vertical Carbody Vibrations of a Rail Vehicle - A Literature Survey”, in: Report in Railway Technology, KTH Engineering Sciences, 2010 [9] V. Barethiye, G. Pohit, A. Mitra, “A combined nonlinear and hysteresis model of shock absorber for quarter car simulation on the basis of experimental data”, Engineering Science And Technology, An International Journal, Vol. 20, No. 6, pp. 1610-1622, 2017 [10] C. Mizrak, Orta Hizli Calışan Rayli Sistem Bojilerinin Titresimlerinin Modellenmesi ve Optimizasyonu, Karabuk University, 2015 (in Turkish) [11] G. R. M. Mastinu, M. Gobbi, G. D. Pace, “Analytical formulae for the design of a railway vehicle suspension system”, Journal Of Mechanical Engineering Science, Vol. 215, No. 6, pp. 683-698, 2001 [12] S. J. Chikhale, D. S. P. Deshmukh, “Comparative Analysis Of Vehicle Suspension System in Matlab-Simulink and MSc- Adams with the help of Quarter Car Model”, International Journal Of Innovative Research In Science, Engineering And Technology, Vol. 2, No. 8, pp. 4074-4081, 2013 [13] P. Sathishkumar, J. Jancirani, D. John, S. Manikandan, “Mathematical modelling and simulation quarter car vehicle suspension”, International Journal of Innovative Research in Science, Engineering and Technology, Vol. 3, No. 1, pp. 1280-1283, 2014 [14] K. K. Jagtap, D. R. Dolas, “Simulation of Quarter Car Model”, International Journal Of Engineering Research And General Science, Vol. 3, No. 6, pp. 242-249, 2015 [15] A. Florin, M. R. Ioan-Cozmin, P. Liliana, “Pasive Suspension Modeling Using MATLAB, Quarter Car Model, Imput Signal Step Type”, New Technologies And Products In Machine Manufacturing Technologies, Romania, 2013 [16] N. Ouelaa, A. Rezaiguia, B. Laulagnet, “Vibro-acoustic modelling of a railway bridge crossed by a train”, Applied Acoustics, Vol. 67, No. 5, pp. 461-475, 2006 [17] R. U. A. Uzzal, W. Ahmed, S. Rakheja, “Dynamic Analysis of Railway Vehicle-Track Interactions due to Wheel Flat with a Pitch-Plane Vehicle Model”, Journal of Mechanical Engineering, Vol. 39, No. 2, pp. 86-94, 2008 [18] M. Metin, R. Guclu, “Vibrations control of light rail transportation vehicle via PID type fuzzy controller using parameters adaptive method”, Turkish Journal of Electrical Engineering and Computer Sciences, Vol. 19, No. 5, pp. 807-816, 2011 [19] E. J. O’Brien, C. Bowe, P. Quirke, D. Cantero, “Determination of longitudinal profile of railway track using vehicle-based inertial readings”, Journal of Rail and Rapid Transit, Vol. 231, No. 5, pp. 518- 534, 2016 [20] C. Mizrak, I. Esen, “Determining Effects of Wagon Mass and Vehicle Velocity on Vertical Vibrations of a Rail Vehicle Moving with a Constant Acceleration on a Bridge Using Experimental and Numerical Methods”, Shock And Vibration, Vol. 2015, Article ID183450, 2015 [21] X. Xu, H. Jiang, M. Gao, “Modeling and Validation of Air Suspension with Auxiliary Chamber Based on Electromechanical Analogy Theory”, Applied Mechanics And Materials, Vol. 437, No. 1, pp. 190-193, 2013 [22] D. Hroncova, M. Binda, P. Sarga, F. Kicak, “Kinematical analysis of crank slider mechanism using MSC Adams/View”, Procedia Engineering, Vol. 48, pp. 213-222, 2012 [23] A. Tandel, A. R. Deshpande, S. P. Deshmukh, K. R. Jagtap, “Modeling, Analysis and PID Controller Implementation on Double Wishbone Suspension Using SimMechanics and Simulink”, Procedia Engineering, Vol. 97, pp. 1274-1281, 2014 [24] A. Goldman, Magnetic Components for Power Electronics, Springer, 2002