About the preliminary design of the suspension spring and shock absorber

Similar documents
Aspects Concerning Modeling and Simulation of a Car Suspension with Multi-Body Dynamics and Finite Element Analysis Software Packages

Study on Dynamic Behaviour of Wishbone Suspension System

Researches regarding a pressure pulse generator as a segment of model for a weighing in motion system

Comparative study between double wish-bone and macpherson suspension system

The influence of thermal regime on gasoline direct injection engine performance and emissions

Chapter 2 Dynamic Analysis of a Heavy Vehicle Using Lumped Parameter Model

Hardware structures of hydronic systems for speed control

Influence of Parameter Variations on System Identification of Full Car Model

KINEMATICAL SUSPENSION OPTIMIZATION USING DESIGN OF EXPERIMENT METHOD

Study concerning the loads over driver's chests in car crashes with cars of the same or different generation

Mathematical Modelling and Simulation Of Semi- Active Suspension System For An 8 8 Armoured Wheeled Vehicle With 11 DOF

ENERGY RECOVERY SYSTEM FROM THE VEHICLE DAMPERS AND THE INFLUENCE OF THE TANK PRESSURE

Simulation and Analysis of Vehicle Suspension System for Different Road Profile

SUSPENSION OF A MOUNTAIN BIKE SVOČ FST Bc. Vít Prošek University of West Bohemia Univerzitni 8, Pilsen Czech Republic

NOVEL DAMPER FOR PASSIVE SECURITY INCREASING

INTERCONNECTION POSSIBILITIES FOR THE WORKING VOLUMES OF THE ALTERNATING HYDRAULIC MOTORS

ISSN: SIMULATION AND ANALYSIS OF PASSIVE SUSPENSION SYSTEM FOR DIFFERENT ROAD PROFILES WITH VARIABLE DAMPING AND STIFFNESS PARAMETERS S.

Comparison between Optimized Passive Vehicle Suspension System and Semi Active Fuzzy Logic Controlled Suspension System Regarding Ride and Handling

THE INFLUENCE OF THE WHEEL CONICITY ON THE HUNTING MOTION CRITICAL SPEED OF THE HIGH SPEED RAILWAY WHEELSET WITH ELASTIC JOINTS

FLUID FLOW MODELLING OF A FLUID DAMPER WITH SHIM LOADED RELIEF VALVE

Experimental Investigation of Effects of Shock Absorber Mounting Angle on Damping Characterstics

Constructive Influences of the Energy Recovery System in the Vehicle Dampers

Effect of plus sizing on driving comfort and safety of users

Comparing PID and Fuzzy Logic Control a Quarter Car Suspension System

Scroll Compressor Oil Pump Analysis

Comparison Of Multibody Dynamic Analysis Of Double Wishbone Suspension Using Simmechanics And FEA Approach

QUARTER CAR SUSPENSION SYSTEM WITH ONE DEGREE OF FREEDOM SIMULATED USING SIMULINK. L. Bereteu, A. Perescu

MODELS FOR THE DYNAMIC ANALYSIS OF THE SUSPENSION SYSTEM OF THE VEHICLES REAR AXLE

Design and Analysis of Hydraulic Chassis with Obstacle Avoidance Function

Physical Modelling of Mine Blast Impact on Armoured Vehicles

Design of Damping Base and Dynamic Analysis of Whole Vehicle Transportation based on Filtered White-Noise GongXue Zhang1,a and Ning Chen2,b,*

Impact of air conditioning system operation on increasing gases emissions from automobile

Chapter 2. Background

CHAPTER 4: EXPERIMENTAL WORK 4-1

Design of pneumatic proportional flow valve type 5/3

Investigation of Semi-Active Hydro-Pneumatic Suspension for a Heavy Vehicle Based on Electro-Hydraulic Proportional Valve

A REVIEW OF TWO WHEELER VEHICLES REAR SHOCK ABSORBER

The Shock Absorber Handbook Second Edition

Modeling of 17-DOF Tractor Semi- Trailer Vehicle

Active Systems Design: Hardware-In-the-Loop Simulation

THE FREE PISTON ENGINE AS PRIME MOVER FOR CONVERTER BASED VEHICLE PROPULSION SYSTEMS

Modeling, Design and Simulation of Active Suspension System Frequency Response Controller using Automated Tuning Technique

THEORETICAL INVESTIGATIONS ABOUT THE INFLUENCES OF THE PISTON PIN STRESS ON THE COMBUSTION ENGINE PERFORMANCES

The Modeling and Simulation of DC Traction Power Supply Network for Urban Rail Transit Based on Simulink

Semi-Active Suspension for an Automobile

MAGNETORHEOLOGICAL FLUID DAMPER WITH EXTERNAL MAGNETIC CIRCUIT. Alexandru Boltosi, Adrian Chiriac, Liviu Bereteu, Ramona Nagy, Voichiţa Hule

The analysis of the accuracy of the wheel alignment inspection method on the side-slip plate stand

Steering Dynamics of Tilting Narrow Track Vehicle with Passive Front Wheel Design

A CONTRIBUTION TO PASSIVE SHOCK ABSORBER FORCE INVESTIGATION

Technical elements for minimising of vibration effects in special vehicles

Determination of anti pitch geometry. acceleration [1/3]

Test Rig Design for Measurement of Shock Absorber Characteristics

Active Suspensions For Tracked Vehicles

Study of intake manifold for Universiti Malaysia Perlis automotive racing team formula student race car

The Application of Simulink for Vibration Simulation of Suspension Dual-mass System

Influence of shock absorber temperature on vehicle ride comfort and road holding

Development and Control of a Prototype Hydraulic Active Suspension System for Road Vehicles

FEASIBILITY STYDY OF CHAIN DRIVE IN WATER HYDRAULIC ROTARY JOINT

SIMULATION OF A SPARK IGNITION ENGINE WITH CYLINDERS DEACTIVATION

NUMERICAL ANALYSIS OF LOAD DISTRIBUTION IN RAILWAY TRACK UNDER WHEELSET

EFFECTIVE SOLUTIONS FOR SHOCK AND VIBRATION CONTROL

Development of analytical process to reduce side load in strut-type suspension

Clamping Force Effects on the Behaviour of Asymmetrical Friction Connections (AFC)

MARINE FOUR-STROKE DIESEL ENGINE CRANKSHAFT MAIN BEARING OIL FILM LUBRICATION CHARACTERISTIC ANALYSIS

STIFFNESS CHARACTERISTICS OF MAIN BEARINGS FOUNDATION OF MARINE ENGINE

Plastic Ball Bearing Design Improvement Using Finite Element Method

Optimization of Seat Displacement and Settling Time of Quarter Car Model Vehicle Dynamic System Subjected to Speed Bump

COMPUTATIONAL MODELING OF HEAVY DUTY TRUCK DRIVESHAFT

The operating principle and experimental verification of the hydraulic electromagnetic energy-regenerative shock absorber

Design and Analysis of suspension system components

International Journal of Scientific & Engineering Research, Volume 7, Issue 3, March ISSN DESIGN AND ANALYSIS OF A SHOCK ABSORBER

Relevant friction effects on walking machines

Control and operating conditions and hydrokinetic converter slip in the vehicle s transmission system

Modeling, Design and Simulation of Active Suspension System Root Locus Controller using Automated Tuning Technique.

Registration of the diagnostic signals of the starting system for selected faults

MODELING OF SUSPENSION SYSTEM OF A LIGHT TRACKED VEHICLE

Control and Simulation of Semi-Active Suspension System using PID Controller for Automobiles under LABVIEW Simulink

Particular bi-fuel application of spark ignition engines

MECA0494 : Braking systems

Vibration Measurement and Noise Control in Planetary Gear Train

Experimental Characterization of Gas Filled Hydraulic Damper Using Ramp Excitation

2. Write the expression for estimation of the natural frequency of free torsional vibration of a shaft. (N/D 15)

Shimmy Identification Caused by Self-Excitation Components at Vehicle High Speed

Simulation of Dynamics of System with Hydraulic Lines and Linear Hydraulic Motor with Mass Load

INVESTIGATION OF FRICTION COEFFICIENTS OF ADDITIVATED ENGINE LUBRICANTS IN FALEX TESTER

VALIDATION OF A HUMAN-AND-HARDWARE-IN-THE- LOOP CONTROL ALGORITHM

Vibration Analysis of an All-Terrain Vehicle

Faculty Code: AU13. Faculty Name: RAJESH. M. Designation: LECTURER

The test bench for testing torsional stiffness of active anti-roll bar made of extended profiles with rectangular cross-section

EFFECTIVENESS OF THE ACTIVE PNEUMATIC SUSPENSION OF THE OPERATOR S SEAT OF THE MOBILE MACHINE IN DEPEND OF THE VIBRATION REDUCTION STRATEGIES

1036. Thermal-hydraulic modelling and analysis of hydraulic damper for impact cylinder with large flow

2764. Outer characteristic simulation and performance analysis of variable shock absorber

Research on Damping Characteristics of Magneto-rheological Damper Used in Vehicle Seat Suspension

Selection criteria of the addendum modification coefficients of spur gear pairs with smaller number of pinion teeth

MODELLING AND SIMULATION OF HYDROPNEUMATIC SUSPENSION FOR A CAR

STUDY REGARDING THE MODELING AND SIMULATION ON THE INFLUENCE OF AUTOMOBILE BRAKE SYSTEMS ON ACTIVE SAFETY

A STUDY OF HYDRAULIC RESISTANCE OF VISCOUS BYPASS GAP IN MAGNETORHEOLOGICAL DAMPER

Dynamic Simulation of Vehicle Suspension Systems for Durability Analysis

DESIGN AND ANALYSIS OF SPRING SUSPENSION SYSTEM

Magneto Rheological Semi-Active Damper with External By-pass Circuit in Modular Structure

Transcription:

IOP Conference Series: Materials Science and Engineering PAPER OPEN ACCESS About the preliminary design of the suspension spring and shock absorber To cite this article: I Preda 2016 IOP Conf. Ser.: Mater. Sci. Eng. 147 012128 View the article online for updates and enhancements. Related content - The preliminary design of the optical Thomson scattering diagnostic for the National Ignition Facility P Datte, J S Ross, D Froula et al. - Preliminary Design of Aerial Spraying System for Microlight Aircraft Zamri Omar, Nurfazliawati Idris and M.Zulafif Rahim - Preliminary design of a supercritical CO2 wind tunnel B Re, A Rurale, A. Spinelli et al. This content was downloaded from IP address 148.251.232.83 on 12/07/2018 at 08:52

About the preliminary design of the suspension spring and shock absorber I Preda 1 1 Automotive and Transportation Department, Transilvania University of Brasov, Brasov, Romania E-mail: pion@unitbv.ro Abstract. The aim of this paper is to give some recommendation for the design of main-spring and shock absorber of motor vehicle suspensions. Starting from a 2DoF model, the suspension parameters are transferred on the real vehicle on the base of planar schemes for the linkage. For the coil spring, the equations that must be fulfilled simultaneously permit to calculate three geometrical parameters. The indications presented for the shock absorber permit to obtain the damping coefficients in the compression and rebound strokes and to calculate the power dissipated during the vehicle oscillatory movement. 1. Introduction The aim of this work is to give the main directions to be followed during the earliest design stages of passive-suspension s spring and shock absorber, or, in other words, to indicate the main steps needed to be done when the design starts from scratch. The design of a vehicle suspension system starts with very few input parameters. Simple models are used during initial simulations in order to ensure the desired compromise between comfort and dynamic performance qualities, at different vehicle speeds and loads. That stage leads to the setup of the needed suspension parameters on the model and vehicle. 2. The suspension model and its correspondence with the real vehicle To design the vehicle suspension, a set of input data are necessary. To obtain this information, the automotive engineers imagine and use models to simulate the ride behaviour of the suspension. Corresponding to the needs, the models can have different complexities. For basic studies of vehicle s ride behaviour and quality, the most used model is the so-called quarter-vehicle model or vehicle-corner model [17], [2], [4], [19]. As presented in the left side of figure 1, this includes only two elements with concentrated inertial properties: one for a wheel and the other for the part of the body supported by that wheel. Because normally the up-and-down movement is of most interest, only the mass of those inertial elements will be considered (neglecting the moments of inertia involved in the oscillation). That means the model will have only two degrees of freedom, the vertical translation of the masses m W of the wheel (unsprung mass) and m B of the vehicle quarter (sprung mass). In the model, the vehicle body is linked to the wheel through massless spring and damper, while the connection between the wheel and the ground is realised by the tyre (represented here as a spring-damper combination). Suited for frequencies up to 30 50 Hz (bigger of 2 to 4 times as the natural frequency of the unsprung mass), this model has only six parameters: the masses, the stiffness k S of the suspension spring, the damping coefficient c D of the shock absorber, and the Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI. Published under licence by Ltd 1

stiffness and damping coefficient of the tyre, k T and c T. Due to the mechanical properties of the tyre rubber and of the low frequencies of interest (0 25 Hz) in the case of the ride study, the tyre damping may be neglected (c T 0). Accordingly with the destination of the vehicle, the design engineer will set-up the best values for the model parameters k T, k S and c D. In the literature there are indications about how to do that by simulation in time or frequency domain in the case of deterministic of stochastic excitation. A possible way to do that it is presented in [9]. The next step after the tuning of the quarter-car model is to transfer or translate the values of the design parameters from the model to the real vehicle, i.e. the values k S of the suspension main spring stiffness, c D of the shock absorber damping and k T of the tyre stiffness. To realize the correspondence between all the parameters of the model and the real vehicle it is necessary to obtain almost similar dynamic behaviours for the both systems: the mechanical energy (both, kinetic and potential components) and theirs rates of change must be equal at any moment for the dynamic model and the real vehicle. The results of parameters translation from the model to the vehicle depend on the type, layout and dimensions of the suspension linkage. How the values of the coefficients k S and c D are changed during this transfer (while the tyre stiffness k T remains the same) it is presented in [8]. In that paper, the demonstration is based on the generally adopted (even though it is not entirely true) hypothesis of planar mechanisms approximating some of the most spread independent suspension linkages. Such a scheme of calculus it is presented in the right side of figure 1. Figure 1. Calculus scheme of the short-long-arm suspension linkage and its 2DoF model. Starting from the values of the real vehicle s tyre stiffness and shock absorber damping, some recommendations regarding the design of the coil spring and shock absorber will be made further. 3. The coil spring The suspension spring design can be realized relying on machine elements publications such as [18], [6], or automotive design publications, such as [17], [2], [1]. Accordingly to the literature, the next three equations must be simultaneously fulfilled: D m = R dia d; τ = F [8 c cor (R dia ) D m )]/( d 3 ); k = (G d 4 )/(8 N c D m 3 ) (1) The notations in this equation system represent: d wire diameter; D m mean coil (wounding) diameter; E modulus of elasticity (Young s modulus); ν Poisson s ratio; G=E/[2 (1+ν)] shear modulus; τ torsion stress; R dia ratio of coil and wire diameters; c cor stress correction factor (as 2

function of R dia ) considering the uneven distribution of torsion stress in the case of variable spring load (which is the case of vehicle suspension); N c number of active coils; F compression force; k spring stiffness. The spring dimensions (d, D m ) and the force F are illustrated also in figure 2. Figure 2. Scheme with the calculus elements of the coil spring The stress correction factor c cor (also named Whal s factor) depends on R dia =D m /d and can be calculated with the theoretical derived equation [18] c cor = (4R dia - 1)/ (4R dia - 4) + 0.615/R dia (2) or can be taken from tables with experimentally determined data [1]. From the equations (1), three unknowns from the previous list of magnitudes can be determined, while the others must be adopted. For example, the three unknowns may be: d, D m and N c. In that case, the material properties must be adopted (the shear modulus G and the admissible torsion stress τ adm ), while the stiffness coefficient k (the spring rate) is already known (as presented previously). The usual maximum dynamic load on the tyre is indicated in the literature: about 2 times [17] up to 2.5 times [11] bigger than the static load. That means the maximal compression force F max acting on the coil spring will not exceed 2 times the static load, because at the extremities of the linkage stroke, the main spring is helped (protected) by rubber or polyurethane spring buffers. For example, in the case of the MacPherson suspension presented in [5], the maximum force acting on the spring is approximately 1.7 F st. For other spring types, as leaf and torsion bar springs, the equation system (1) will be replaced by other specific equations, but involving also dimensional and material (strength) parameters. 4. The shock absorber At low frequencies, the shock absorber is the suspension main element reducing the energy of oscillations. It serves to dissipate the kinetic and potential energy existing in the suspension during relative wheel-body movements. The shock absorber is an essential element to realize a good compromise between a low level of the vibrations transmitted from the ground to the vehicle body (comfortability) and an adequate control of the unsprung and sprung masses, in order to keep good dynamic performances (active safety and handling). The working principle of the modern shock absorbers consists in generating friction forces opposing to the relative movements of the damper ends by forcing a quantity of oil to pass through calibrated orifices, as illustrated in figure 3. Many of the automotive specialists recommend, and motor-vehicle manufacturers adopt, designs with asymmetric characteristics that offer, for the same relative speed, smaller damping force on the compression stroke as in the rebound (extension) stroke. That behaviour will reduce the summative force transmitted together by the spring and shock absorber to the body when the wheel escalates a bump, while at the rebound the damper will consume most of the potential energy accumulated in the spring. If c c and c r are the damping coefficients of the shock absorber in the compression and rebound strokes, theirs ratio R f =c r /c c (3) 3

may have values in the interval 2 5 [17], [7], [19], with a tendency at the modern vehicles to limit that value to 2 3, in order to avoid the body sinking when the vehicle starts passing over obstacles [4] and not to impede the tyre to maintain the contact with the ground. Figure 3. Simple principle scheme presenting the oil circulation into a shock absorber When the oil passes through orifices with constant flowing section, the damping force is practically proportional with the square of the relative speed, F=sgn(v) c v 2. (4) But, with modern shock absorber designs, that parabolic force-speed characteristic can be transformed due to supplementary flowing sections opened by complex control valves when the oil pressure difference in the working chambers attains certain values. Therefore, using valves, the actual characteristic shape F=F(v) may be very different as the natural parabolic one from equation (4). Today, the characteristics are preferred almost linear (F=c v), figure 4, or even regressive, in order to reduce the possibility of excessive forces apparition. Figure 4. Schematic typical work diagram (left) and characteristic diagram (right) of a modern bi-linear shock absorber There are many specialized publications presenting models for calculating different types of shock absorbers: [17], [7], [10], [14], [15], [3], [13], [16], that may be used for mono- or bi-tubular designs; without or with separation piston; with normal- or magnetorheological oil. These theoretical models are then verified in practice. At some designs with expansion reservoir, the influence of the pressured gas may be considered in a more complex model of the shock absorber or as a supplementary spring added to the main one. To experimentally obtain the characteristics, the damper is excited with sinusoidal displacement u and speed v [12]: u =a sin(ω t) =a sin(2 f t); v =du/dt =a ω cos(ω t) =a 2 f cos(2 f t) (5) where a is a constant displacement amplitude and f and ω are constant frequency and pulsation. By varying the frequency in several steps, different force-displacement curves F=F(u), named work diagrams, are obtained, as shown in figure 4. Then, by taking the peak values of the damper force at the displacement u=±a, which corresponds with the velocity v=±2 f a, the curve F=F(v), named 4

characteristic diagram, is generated. Normally, the half-cycle representing the rebound stroke is associated with negative damper velocities and forces. In [4] it is demonstrated that an individual characteristic diagram of any shock absorber can be replaced for simulation purpose with an equivalent bilinear one: F c = c c v if v 0; F r = c r v if v < 0; R f =c r /c c =constant (6) At that point, knowing the mass m B of the vehicle body, the spring stiffness k S and adopting the damping ratio (7) it can be calculated the equivalent damping coefficient c D, considered constant for the both strokes of the shock absorber. The instantaneous power dissipated by the shock absorber is P D (t) = F(t) v(t) = c v 2 (t). (8) that, for sinusoidal displacement, will become P Ds (t) = c a 2 ω 2 cos 2 (ω t). (9) Further, integrating this power P Ds (t) over a half cycle it obtains the dissipated energy L s1/2 = P Ds (t) = c a 2 ω /4. (10) The problem is now to realize the equivalence of the constant coefficient c D with the coefficients c c and c r in the equation (6). The logical possibility to make this approximation must take into account the equality of the energy consumed in an oscillation cycle by the equivalent and by the real shock absorbers: 2 c D a 2 ω /4 = (c c + c r ) a 2 ω /4 which means c D = (c c + c r )/2, (11) i.e. the equivalence equation indicated in the literature. Now, with the equations (6) and (11), the true compression and rebound damping coefficients (c c and c r ), of the shock absorber can be calculated. For a sinusoidal excitation with the frequency f [cycles/s], the mean power dissipated by the suspension (transformed in heat) is P Dmean = f (c c + c r ) a 2 ω /4 = f c D a 2 ω /8 (12) and may be used for the calculus of the thermal balance. Because the damping coefficients c c and c r in the compression and rebound strokes are known now, the design of the shock absorber may continue with the dimensioning of the piston orifices (as presented in [17] and [7]) and also with the strength calculus op the parts, because the mean and maximal forces produced by the shock absorber (F c and F r ) can be estimated. 5. Conclusions For the study of the ground vehicles independent suspensions, simple dynamic and mathematic models are used, mainly in the first stages of the design process. After the tuning of the suspension parameters realized on the model, the suspension stiffness and the damping coefficient are transferred (with bigger values) to the real vehicle s main-spring and shock absorber. Based on these necessary values, in the paper are presented some recommendations for the design process of the suspension spring and shock absorber. These may be useful to establish both the spring s wire and mean coil wounding diameters and the shock absorber s damping coefficients and mean power dissipated at different working regimes. 5

6. References [1] Bosch Automotive Handbook 2011 8th Edition (Plochingen: Robert Bosch GmbH). [2] Câmpian V, Vulpe V, Ciolan G, Enache V, Preda I and Câmpian O V 1989 Automobile (Brasov: Universitatea din Brasov) [3] Dixon J C 2007 The Shock Absorber Handbook Second Edition (Chichester: Wiley) [4] Genta G and Morello L 2009 The Automotive Chassis Vol. 2 - System Design. (Springer) [5] Merel J and Deschaume A 2003 Filtrage d'une jambe de force de type Mac Pherson Ingénieurs de l'automobile Mars 2003 pp 44-48 [6] Mogan G L and Butnariu S L 2013 Organe de masini. Teorie-Proiectare-Aplicatii - sistem integrat Second Edition (Brasov: Universitatea Transilvania din Brasov) [7] Pham A T 1987 La dynamique automobile. Calcule des elements de suspension. La fonction elastique. La fonction d'amortissement (Paris: CLESIA, C87-02); [8] Preda I. 2016 Independent Suspension The Equivalence of Model and Vehicle Parameters. Proceedings of the European Automotive Congress EAEC-ESFA 2015 (Springer International Publishing) Chapter 60, pp 651-664 [9] Preda I 2016 Model Based Algorithm for the Study of the Vehicle Suspension Applied Mechanics and Materials 823 pp 247-252 [10] Radulescu C. 2002 Contributii privind studiul suspensiei semiactive a autoturismelor Ph.D. Thesis (Bucharest: Military Technical Academy) [11] Reimpell J, Stoll H and Betzler J W 2001 The Automotive Chassis. Engineering Principles. Second Edition (Oxford: Butterworth-Heinemann) [12] Rill G. 2007 Vehicle Dynamics - Fundamentals and Modeling Aspects. (Regensburg: Fachhochschule Regensburg) [13] Sandu C Southward S and Richards R 2010 Comparison of Linear, Nonlinear, Hysteretic, and Probabilistic Models for Magnetorheological Fluid Dampers ASME J of Dynamic Systems Measurement and Control, Special Issue on Physical System Modeling 132(6) pp 164-172 [14] Simms A and Crolla D 2002 The influence of damper properties on vehicle dynamic behaviour SAE Technical Papers 2002-01-0319 [15] Sireteanu T and Ghita G 2004 Controlul inteligent al vibratiilor mecanice prin fluide magnetoreologice cu implementarea automatizarilor complexe cu logica programabila (Bucharest: Institutul de Mecanica Solidelor al Academiei Romane, Chapter 2003-2004) [16] Tabacu C. 2013 Contributii la îmbunatatirea confortului si stabilitatii vehiculelor rutiere, prin optimizarea caracteristicilor amortizoarelor de suspensie Ph.D. Thesis (Bucharest: Politehnica University of Bucharest). [17] Untaru M Fratila G Potincu G Seitz N Tabacu I Peres G and Tiberiu M. 1982 Calculul si constructia automobilelor. (Bucharest: Editura Didactica si Pedagogica) [18] Wahl A M 1944 Mechanical springs (Cleveland: Penton Publishing) [19] Wong J Y 2001 Theory of Ground Vehicles Third Edition (Ottawa: Wiley) 6