IOSR Journal of Computer Engineering (IOSR-JCE) e-issn: 2278-0661,p-ISSN: 2278-8727, Volume 17, Issue 2, Ver. 1 (Mar Apr. 2015), PP 65-74 www.iosrjournals.org Car Dynamics using Quarter Model and Passive Suspension, Part VI: Sprung-mass Step Response Galal Ali Hassaan (Emeritus Professor, Department of Mechanical Design & Production Engineering, Faculty of Engineering/Cairo University, Giza, Egypt) Abstract: The objective of the paper is to investigate the step response of a 2 DOF quarter-car model with passive suspension. The mathematical models of the sprung-mass displacement and acceleration as response to the step road disturbance are derived. The step response of the system is plotted using MATLAB for a 100 mm step amplitude. The effect of the suspension parameters on the step response of the sprung-mass is investigated in terms of motion displacement and acceleration. The maximum percentage overshoot and settling time of the step response are evaluated as time-based characteristics of the dynamic system for a specific range of the suspension parameters. The maximum acceleration of the sprung-mass due to the step road disturbance did not conform with the ride comfort requirements according to ISO 2631-1 standard. Keywords: Car dynamics; quarter car model; passive suspension; step response, sprung-mass acceleration I. Introduction Car dynamics can be studied in various ways including time response to speed humps, frequency response and step response. The latest is borrowed from automatic control engineering to investigate the dynamic system in the time domain using the most known input type. This type of input is used by a large number of researchers dealing with car dynamics. Masi (2001) used computer simulation to examine the effect of control methods on the performance of semiactive dampers in controlling the dynamics of a quarter-car model. One of the inputs he used was a step input of 12.7 mm amplitude [1]. Askerda et. al. (2003) developed a methodology for analyzing the impact of data errors on control system dependability. Their applied their technique on a car slip-control brake-system. They used a step input command signal in terms of the normalized car-slip [2]. Toshio and Atsushi (2004) presented the construction of a pneumatic active suspension system for a one-wheel car model using fuzzy reasoning and disturbance observer. They used a step input of 1 V amplitude to the actuator [3]. Sakman, Guclu and Yagiz (2005) investigated the performance of fuzzy logic controlled active suspension on a nonlinear model with 4 DOF. They demonstrated the effectiveness of the fuzzy logic controller for active suspension systems. They used a half-car model with 4 DOF. They used a step road surface input of 35 mm amplitude resulting in chassis maximum acceleration of about 3.5 m/s 2 [4]. Faheem, Alam and Thomson (2006) investigated the response of quarter and half-car models to road excitation. They studied the time response to a step road input and the frequency response of the sprung and unsprung masses [5]. Maila and Priyandoko (2007) presented the design of a control technique for active suspension system of a quarter-car model using adaptive fuzzy logic and active force control. They tested their control strategy using a 50 mm step road input resulting in a 75 mm maximum sprung-mass displacement and 20 m/s2 maximum acceleration using passive elements. The active force control techniques reduced the maximum acceleration to about 10 m/s 2 [6]. Hanafi (2009) designed a fuzzy logic control system for a quartercar model. He used a step road disturbance of 100 mm amplitude [7]. Alexandru and Alexandru (2010) presented a virtual prototype of an automobile suspension system with force generator actuator. They used a quarter-car model and a step input of 100 mm amplitude. The time response of the sprung-mass exhibited an overshoots and undershoots [8]. Unaune, Pawar and Mohite (2011) simulated the time response of a quarter-car model to a step input with passive elements suspension. They investigated the effect of mass ratio, stiffness ratio and damping coefficient ratio on the sprung-mass dynamics. They didn t show the amplitude of the step input used in the simulation [9]. Fayyad (2012) presented a control system for active suspension system to improve ride comfort even at resonance. He used a step road input of 80 mm amplitude to a quarter-car model. The passive elements resulted in a sprung-mass overshoot of 59 % while with the active control it was 22.5 %. The maximum acceleration to the step input used was 17 m/s2 with passive elements and 1.3 m/s 2 with active suspension [10]. Chikhale and Deshmukh (2013) performed a vibration analysis of a quarter-car model to a step input using MATLAB and MSc. ADAMS. They used the system state space representation in deriving the mathematical model of the system. They used a step input of 25 mm amplitude. The sprung-mass had 62 % overshoot in its displacement response to the step input [11]. Faruk, Bature, Babani and Dankadai (2014) used PID conventional DOI: 10.9790/0661-17216574 www.iosrjournals.org 65 Page
controller and integrated fuzzy logic control schemes for good driving comfort and effective isolation of road disturbance. They simulated the quarter-car step input response for a 1 m step amplitude!!. The time response without control had a 75 % overshoot reduced to zero with PID and fuzzy controllers [12]. Patole and Swant (2015) studied the efforts of researchers about nonlinearity of passive suspension parameters. They referred to ride comfort through vibration amplitude and ride comfort level for model inputs including a step type [13]. II. The Quarter-Car Model Researchers use quarter car model in studying the vehicle dynamics. This model depending on the degree of simplification used by the researcher may be considered as a single (1DOF), two (2DOF) or three degree of freedom (3DOF) [14,15,16]. A quarter-car model with 2 DOF will be considered in this work because of the complexity of its mathematical model if 3 DOF one is considered. The line diagram of the 2 DOF quartercar model is shown in Fig.1 [17]. Fig. 1 2 DOF quarter-car model [17]. The model of Fig.1 consists of: - One tire and wheel and its attachments having the parameters: mass (m 1 ), damping coefficient (c 1 ) and stiffness (k 1 ). - Sprung mass which is one quarter of the car chassis having the parameters: mass (m 2 ), damping coefficient (c 2 ) and stiffness (k 2 ). - c 2 and k 2 are the suspension parameters. III. Model Parameters The 2 DOF quarter-car model parameters used in this analysis work are given in Table 1. The suspension parameters k 2 and c 2 are set in the range given in the table to investigate their effect on the step response of the vehicle. The tire and wheel parameters (m 1, k 1 and c 1 ) and sprung mass (m 2 ) are borrowed from reference [17]. Table 1 2 DOF model parameters Parameters Description Value k 1 (N/m) Tire stiffness 135000 c 1 (Ns/m) Tire damping coefficient 1400 m 1 (kg) Un-sprung mass 49.8 k 2 (N/m) Suspension stiffness 1000 3000 c 2 (Ns/m) Suspension damping coefficient 1000 4000 m 2 (kg) Sprung mass 466.5 IV. Sprung-Mass Displacement The mathematical model of the 3 DOF quarter-car model is defined by three differential equations obtained by applying Newton's third law to the free body diagram of the three lumped masses of Fig.1. The dynamic motions of the system are: x 1, x 2 and x 3. The irregularities of the road is defined by the dynamic motion y which is assumed sinusoidal. Of course, the irregularities are random. The frequency contents of this irregularity can be defined using FFT (fast Fourier transform). Thus, the sinusoidal input motion y represents one component in this transform. All the parameters of the system are passive and assumed having the constant values presented in section 3. The differential equations are written in standard form as follows: For the tire and wheel mass, m 1 : m 1 x 1 '' + (c 1 +c 2 )x 1 ' c 2 x 2 ' + (k 1 +k 2 )x 1 k 2 x 2 = c 1 y' +k 1 y (1) For the sprung mass, m 2 : DOI: 10.9790/0661-17216574 www.iosrjournals.org 66 Page
m 2 x 2 '' - c 2 x 1 ' + c 2 x 2 ' - k 2 x 1 + k 2 x 2 = 0 (2) The Laplace transform of Eqs.1 and 2 gives: [m 1 s 2 + (c 1 +c 2 )s + k 1 +k 2 ]X 1 (s) (c 2 s+k 2 )X 2 (s) = (c 1 s+k 1 )Y(s) (3) And (m 2 s 2 +c 2 s+k 2 )X 2 (s) (c 2 s+k 2 )X 1 (s) = 0 (4) Where X 1 (s) and X 2 (s) are the Laplace transform of the motions x 1 and x 2 respectively, and s is the Laplace operator. Solving Eqs.3 and 4 for X2(s)/Y(s) gives: X 2 (s)/y(s) = (c 1 s+k 1 )(c 2 s+k 2 ) / {[m 1 s 2 + (c 1 +c 2 )s + k 1 +k 2 ] (m 2 s 2 +c 2 s+k 2 ) (c 2 s+k 2 ) 2 } (5) Eq.5 is the transfer function of the dynamic system. Writing Eq.5 in a standard form gives: X 2 (s)/y(s) = N 1 (s) / D 1 (s) (6) Where: N 1 (s) = c 1 c 2 s 2 + (c 1 k 2 +c 2 k 1 )s + k 1 k 2 and D 1 (s) = m 1 m 2 s 4 + [m 1 c 2 +m 2 (c 1 +c 2 )]s 3 + [m 1 k 2 +c 2 (c 1 +c 2 )+m 2 (k 1 +k 2 )-c 2 2 ]s 2 + [k 2 (c 1 +c 2 )+c 2 (k 1 +k 2 )-2c 2 k 2 ]s + [k 2 (k 1 +k 2 )-k 2 2 ] V. Sprung-Mass Acceleration - In the derivation of the sprung mass displacement, the variables x, x, x and y are used. - For sake of the mathematical model of the sprung mass acceleration, the acceleration variables a 1, a 2, velocity variables v 1, v 2 and displacement variables x 1, x 2 will be used. - The differential equations of the dynamic system (Eqs.1 and 2) in terms of the acceleration, velocity and displacement variables will take the form: m 1 a 1 + (c 1 +c 2 )v 1 c 2 v 2 + (k 1 +k 2 )x 1 k 2 x 2 = c 1 y + k 1 y (7) And m 2 a 2 - c 2 v 1 + c 2 v 2 k 2 x 1 + k 2 x 2 = 0 (8) - The acceleration, velocity and displacement variables are related through: a 1 = v 1 (9) a 2 = v 2 (10) v 1 = x 1 (11) v 2 = x 2 (12) - Taking the Laplace transform of Eqs.7 through 12 gives the Laplace transform ratio A2(s)/Y(s) as: A 2 (s)/y(s) = N 2 (s) / D 2 (s) (13) Where: N 2 (s) = c 1 c 2 s 4 + (k 1 c 2 + k 2 c 1 ) s 3 + k 1 k 2 s 2 (14) D 2 (s) = m 1 m 2 s 4 + [m 2 (c 1 +c 2 ) + m 1 c 2 ] s 3 + [m 2 (k 1 +k 2 ) + c 2 (c 1 + c 2 ) + m 1 k 2 c 2 2 ] s 2 + [c 2 (k 1 +k 2 ) + k 2 (c 1 +c 2 ) 2k 2 c 2 ] s + k 2 (k 1 +k 2 ) k 2 2 (15) VI. Sprung-Mass Step Displacement Response The sprung mass displacement is evaluated using the transfer function of the spring-mass displacement relative to the ground motion (Eq.6). MATLAB command 'step' is used for this purpose for an input step of 100 mm amplitude [18]. The dynamics of the quarter-car model is evaluated in response to a step road input for suspension parameters in the range: 1000 c 2 4000 Ns/m, 1000 k 2 3000 N/m (16) The step response of the sprung mass for the suspension parameters in the range of Eq.16 is shown in Figs.2 through 6. DOI: 10.9790/0661-17216574 www.iosrjournals.org 67 Page
Fig.2 Sprung-mass step response for k 2 = 1000 N/m. Fig.3 Sprung-mass step response for k 2 = 1500 N/m. Fig.4 Sprung-mass step response for k 2 = 2000 N/m. DOI: 10.9790/0661-17216574 www.iosrjournals.org 68 Page
Fig.5 Sprung-mass step response for k 2 = 2500 N/m. Fig.6 Sprung-mass step response for k 2 = 3000 N/m. - The step response of the sprung-mass exhibits an overshoot and undershoot depending of the level of the suspension damping coefficient. - The maximum overshoot and maximum undershoot decreases as the damping coefficient increases. - All the step response curves of the sprung-mass displacement intersect in more than one point. - The effect of the suspension parameters c 2 and k 2 on the maximum percentage overshoot of the sprung-mass step response is shown graphically in Fig.7. Fig.7 Sprung-mass maximum displacement overshoot. DOI: 10.9790/0661-17216574 www.iosrjournals.org 69 Page
- The maximum percentage overshoot decreases as the suspension damping coefficient increases and the suspension stiffness decreases. - The settling time of the sprung-mass displacement is shown in Fig.8. Fig.8 Sprung-mass displacement settling time. - The settling time decreases as the suspension damping coefficient increases. - The effect of the suspension stiffness is not uniform. - versed between the first and second intersections. VII. Sprung-Mass Step Acceleration Response - The step response of the sprung-mass acceleration is simulated using Eq.13 and the 'step' command of MATLAB. - The step response in terms of the sprung-mass acceleration is in Figs.9-13 for suspension stiffness of 1000, 1500, 2000, 2500 and 3000 N/m respectively. Fig.9 Sprung-mass step response (acceleration) for k 2 = 1000 N/m. DOI: 10.9790/0661-17216574 www.iosrjournals.org 70 Page
Fig.10 Sprung-mass step response (acceleration) for k 2 = 1500 N/m. Fig.11 Sprung-mass step response (acceleration) for k 2 = 2000 N/m. Fig.12 Sprung-mass step response (acceleration) for k 2 = 2500 N/m. DOI: 10.9790/0661-17216574 www.iosrjournals.org 71 Page
Fig.13 Sprung-mass step response (acceleration) for k 2 = 3000 N/m. - The sprung-mass acceleration exhibits overshoot and undershoot characteristics relative to the zero acceleration line. - The effect of the suspension parameters on the maximum acceleration of the sprung-mass is shown in Fig.14. Fig.14 Maximum acceleration of the sprung-mass due to the step input. - The maximum acceleration of the sprung mass increases with increased damping. - The maximum acceleration increases almost linearly with the suspension damping coefficient. - The suspension stiffness has almost no effect on the maximum acceleration of the sprung-mass. - Fig.15 shows the effect of the suspension parameters on the acceleration response settling time as evaluated by MATLAB. Fig.15 Settling time of the sprung-mass due to the step input. DOI: 10.9790/0661-17216574 www.iosrjournals.org 72 Page
- The settling time increases to a maximum and then decreases again. - For all values of the suspension stiffness, the settling time is maximum at a damping coefficient of about 2000 Ns/m. - All the characteristic curves of the settling time intersect at a damping coefficient of about 1700 Ns/m. VIII. Conclusion - The dynamic model of the 2 DOF quarter-car dynamic system was derived assuming passive parameters for all dampers and springs in the system. - The step response of the system was presented in terms of the suspension parameters (damping coefficient and stiffness). - Using MATLAB, the step response of the sprung-mass was presented for suspension parameters of damping coefficient in the range 1000 to 4000 Ns/m and stiffness in the range 1000 to 3000 N/m. - The maximum overshoot of the sprung-mass displacement has decreased from 20.30 to 2.73 % when the damping coefficient of the suspension increased from 1000 to 4000 Ns/m at 1000 N/m suspension stiffness. - The settling time of the sprung-mass displacement has decreased from 2.95 to 0.248 s when the damping coefficient of the suspension increased from 1000 to 4000 Ns/m at 1000 N/m suspension stiffness. - The suspension stiffness did not affect the maximum acceleration of the sprung-mass during the step response. - With suspension damping coefficient of 1700 Ns/m, increasing the suspension stiffness decreased the settling time of the sprung-mass acceleration. With damping coefficient > 1700 Ns/m, this action was reversed. - The sprung-mass settling time of its acceleration step response was increased with increased suspension damping coefficient up to 2000 Ns/m after which it started decreasing with less rate. - The maximum acceleration due to the 100 mm step road disturbance was completely extremely uncomfortable according to ISO 2631-1 [19]. - Such sudden road disturbance has to be avoided because of its harmful effect on passengers and vehicles. - The alternative is using more smooth speed control humps such as simple harmonic motion hump [15] or polynomial hump [20]. References [1]. J. Masi, Effect of control techniques on the performance of semiactive dampers, M. Sc. Thesis, Faculty of Virginia Polytechnic Institute and State University, December 2001. [2]. O. Askerdal, M. Gafvert M. Hiller and N. Suri, Analyzing the impact of data errors in safety critical control systems, IEICE Transaction on Information and Systems, E86-D, No.12, 2623-2633, 2003. [3]. Y. Toshio and T. Atsushi, Pneumatic active suspension system for a one wheel car model using fuzzy reasoning and disturbance observer, Journal of Zhejiang University Science, 5 (9), 1060-1068, 2004. [4]. L. Sakman, R. Guclu and N. Yagiz, Fuzzy logic control of vehicle suspension with dry friction nonlinearity, Sadhana, 30 (5), 649-659, 2005. [5]. A. Faheem, F. Alam and V. Thomas, The suspension dynamic for a quarter-car model and half-car model, 3rd BSME-ASME International Conference on Thermal Engineering, 20-22 December, Dhaka, 2006. [6]. M. Maila and G. Priyandoko, Simulation of a suspension system with adaptive fuzzy active force control, International Journal of Simulation Modeling, 6 (1), 25-36, 2007. [7]. D. Hanafi, The quarter car fuzzy controller design based on model from intelligent system identification, IEEE Symposium on Industrial Electronics and Application, Kuala Lumpur, Malaysia, October 4-6, 930-933, 2009. [8]. C. Alexandru and P. Alexandru, The virtual prototype of a mechatronic suspension system with active force control, WSEAS Transactions on Systems, 9 (9), 927-936, 2010. [9]. D. Unaune, M. Pawar and S. Mohite, Ride analysis of quarter vehicle model, Proceedings of the First International Conference on Modern Trends in Industrial Engineering, Surat, Gujarat, India, November, 17-19, 2011. [10]. S. Fayyad, Constructing control system for active suspension system, Contemporary Engineering Sciences, 5 (6), 189-200, 2012. [11]. S. Chikhale and S. Deshm, Comparative analysis of vehicle suspension system in matlab-simulink and MSc-ADAMS with the help of quarter-car model, International Journal of Innovations Research in Science, 2 (10), 5452-5459, 2013. [12]. M. Faruk, A. Bature, S. Batani and N. Dankadai, Conventional and intelligent controller for quadratic car suspension system, International Journal of Technical Research and Applications. 2, Special Issue, 24-27, 2014. [13]. S. Patole and S. Swant, An overview of disarray inside performance analysis of half car model passive vehicle dynamic system subjected to different road profiles with wheel base delay and nonlinear parameters, International Journal of Innovative Research in Advanced Engineering, 2 (1), 63-66, 2015. [14]. G. Litak, M. Borowiec, M. Ali, L. Saha and M. Friswell, Pulsive feedback control of a quarter car model forced by a road profile, Chaos Solutions and Fractals, 33, 1672-1676, 2007. [15]. G. A. Hassaan, Car dynamics using quarter model and passive suspension, Part II: A novel simple harmonic hump, Journal of Mechanical and Civil Engineering, 12(1), 93-100, 2015. [16]. G. A. Hassaan, Car dynamics using quarter model and passive suspension, Part V: Frequency response considering driver-seat, International Journal of Applied Sciences and Engineering Research, 2(1), Accepted for Publication, 2015. [17]. A. Florin, M. Joan-Cosmin and P. Lilian, Passive suspension modeling using MATLAB, quarter-car model, input signal step type, New Technologies and Products in Machine Manufacturing, January, 258-263, 2013. [18]. C. Houpis and S. Sheldon, Linear control systems analysis and design with MATLAB, CRC Press, 2013. [19]. Y. Marjanen, Validation and improvement of the ISO 2631-1 standard method for evaluating discomfort from whole-body vibration in a multi-axis environment, Ph. D. Thesis, Loughborough University, UK, January 2010. DOI: 10.9790/0661-17216574 www.iosrjournals.org 73 Page
[20]. G. A. Hassaan, Dynamics using quarter model and passive suspension, Part III: A novel polynomial hump, Journal of Mechanical and Civil Engineering, 12(1), Version III, 51-57, 2015. Biography Galal Ali Hassaan Emeritus Professor of System Dynamics and Automatic Control. Has got his B.Sc. and M.Sc. from Cairo University in 1970 and 1974. Has got his Ph.D. in 1979 from Bradford University, UK under the supervision of Late Prof. John Parnaby. Now with the Faculty of Engineering, Cairo University, EGYPT. Research on Automatic Control, Mechanical Vibrations, Mechanism Synthesis and History of Mechanical Engineering. Published 10 s of research papers in international journal and conferences. Author of books on Experimental Systems Control, Experimental Vibrations and Evolution of Mechanical Engineering. Reviewer and member of editorial boards of some international journals. DOI: 10.9790/0661-17216574 www.iosrjournals.org 74 Page