Advanced Control and Optimization with Applications to Complex Automotive SystemsView this Special Issue
Research Article | Open Access
Wuhui Yu, Xinjie Zhang, Konghui Guo, Hamid Reza Karimi, Fangwu Ma, Fumiao Zheng, "Adaptive Real-Time Estimation on Road Disturbances Properties Considering Load Variation via Vehicle Vertical Dynamics", Mathematical Problems in Engineering, vol. 2013, Article ID 283528, 9 pages, 2013. https://doi.org/10.1155/2013/283528
Adaptive Real-Time Estimation on Road Disturbances Properties Considering Load Variation via Vehicle Vertical Dynamics
Vehicle dynamics are directly dependent on tire-road contact forces and torques which are themselves dependent on the wheels’ load and tire-road friction characteristics. An acquisition of the road disturbance property is essential for the enhancement of vehicle suspension control systems. This paper focuses on designing an adaptive real-time road profile estimation observer considering load variation via vehicle vertical dynamics. Firstly, a road profile estimator based on a linear Kalman filter is proposed, which has great advantages on vehicle online control. Secondly, to minimize the estimation errors, an online identification system based on the Recursive Least-Squares Estimation is applied to estimate sprung mass, which is used to refresh the system matrix of the adaptive observer to improve the road estimation efficiency. Last, for mining road category from the estimated various road profile sequencse, a road categorizer considering road frequency and amplitude simultaneously is approached and its efficiency is validated via numerical simulations, in which the road condition is categorized into six special ranges, and this road detection strategy can provide the suspension control system with a better compromise for the vehicle ride comfort, handling, and safety performance.
The road properties have significant impact on vehicle performance because vehicle dynamics are directly dependent on tire-road contact forces and torques. Road roughness is a broad term that incorporates everything from potholes and cracks to the random deviations that exist in a profile. The analysis and estimation of a road surface are hot and challenging topics, and the researches related to road disturbances have been presented, which can be classified into three categories: road measurement, road modeling, and road estimation. Each of these categories is reviewed briefly in the sequel.
Road Measurement. It mainly focuses on measuring the road profile accurately for road serviceability, survey, and road maintenance. There are primarily three methods in use internationally. The first category consists of profilometers and profilographs [1–4]. These are manually directed and/or trailer towed mechanisms that directly contact the pavement under evaluation, such as the longitudinal profile analyser (LPA). The second category is inertial profilers such as the Mays Meter  and the GMP (general motor profilometer) , in which a vehicle mounted accelerometer is applied to collect data while traveling at normal speed conditions. González et al. also developed a road roughness estimator by the use of acceleration measurement . Furthermore, Ngwangwa et al. reconstructed road defects and road roughness classification using vehicle responses with artificial neural networks simulation . The third category, equipped with a laser (ultrasonic) transceiver , is more accurate but very expensive.
Road Modeling. It mainly focuses on providing road input source to the virtual vehicle simulation system or the Four-Post Test Rig for vehicle ride analysis. In 1970s, the PSD (power spectral density) function was used by Whitehouse and Archard  and Shinozuka and Jan  to investigate the road roughness. In 1993, Cebon proposed a method based on the IFFT (inverse fast fourier transform) to discretize PSD, which is a simple, fast, and convenient tool for generating road surfaces [12, 13]. In 1995, the contemporary international standard ISO 8608  dealt with road roughness assumes, due to classification of roads into different classes according to their unevenness, equal intensity of road unevenness in the whole range of wavelengths and a general form of the fitted PSD were given. A method based on linear filtering (autoregressive and moving average method or ARMA modelling) was proposed by Yoshimura in 1998 , which has a smaller calculation, and a faster simulation speed, but its precision is not very well. Pazooki et al. , in 2007, summarized different stochastic models of parallel road tracks and evaluated their accuracy by comparing the difference of the measured parallel tracks and the synthetic parallel tracks. In 2012, Hassan and Evans  developed a comprehensive off-road vehicle ride dynamics model considering a random roughness model of the two parallel tracks.
Road Estimation. It mainly focuses on providing real-time rough road estimation for vehicle online control system. Since road input directly affects vehicle vertical suspension dynamics, the availability of suspension sensors such as accelerator and suspension deflection sensor provide an excellent opportunity for road input estimation. Fialho and Balas developed a road adaptive active suspension using linear parameter varying gain-scheduling , in which the road estimation results are applied to the suspension control system to achieve a better balance between ride and handling performance. A new simultaneous input and state estimation algorithm were developed based on the idea of achieving minimum mean square error and minimum error variance . Two stable SISE algorithms were developed based on the minimum variance unbiased estimation technique . A Takagi-Sugeno Fuzzy observer was built for estimating both vehicle dynamics and road geometry . Vehicle sideslip and roll parameters are estimated in presence of the road bank angle and the road curvature as unknown inputs. But it is still not completely ready to be tested in experimental studies since the vehicle speed, parameters variations, and sensor noise are not considered. If the systems are subject to parameter uncertainties, the approaches proposed in [21–23] can be employed. Our previous study [24, 25] also shows the great potential in the enhancement of suspension performance by adopting road estimation in the suspension control system, even though the road estimation algorithm is very simple (just using the statistic value of the accelerator or suspension deflection sensor measurement to classify the road). Recently, a real-time estimation method based on Kalman filter is proposed to estimate the road profile, and experimental results show the accuracy and the potential of the estimation process . But in this research, the vehicle sprung mass change due to vehicle load variation is not considered, which may attenuate the robustness. Imine and Delanne developed a sliding mode observer to estimate the road profiles, which is hard for real-time implementation since a 16 degrees of freedom (DOF) full car model is too complex to be applied online [27, 28].
This paper focuses on developing a road estimation system for an online vehicle control system, which is limited by some important practical requirements. It should, for instance, be: simple enough to run in real time despite onboard processing limitations; reliable enough to operate successfully despite instrumentation failures; robust to variations in vehicle dynamics; fast enough to detect the road input changes when a car is driven on road; easy to incorporate into a control strategy.
With these requirements in mind, the objective of this research has three steps: to develop a road estimator for satisfying the above 5 practical requirements; to obtain vehicle sprung mass online to minimize the road profile estimation errors; and to classify the estimated road profile into several categories according to the main control strategy. The proposed method uses measurements from available sensors: accelerometers and suspension deflection sensors. For simplicity reasons, a quarter-car vehicle model is considered. The estimation process consisting of three blocks is shown in Figure 1.
The first block serves to calculate vehicle sprung mass online from sensor measurements, while the second block contains a Kalman filter that uses the result of the first block as a system parameter adjustor in order to improve robustness of the road estimation system. The third block serves to categorize the estimated road profile elevation into specific types.
The rest of the paper is organized as follows. Section 2 describes a road input state estimator based on a linear Kalman filter. Section 3 presents an online sprung mass estimator based on the Recursive Least-Squares Estimation. In Section 4, a novel road categorizer is proposed, and it is validated via numerical simulations. And, the paper is concluded in Section 5.
2. Road Profile Estimator
2.1. 2-DOF Quarter-Car Model
To implement the Kalman filter method, a suitable vehicle model must be developed. In order to describe the vertical dynamics of a vehicle which runs on an uneven road with a constant speed, a 2-DOF quarter-car model is represented in Figure 2. The quarter-car model does not consider the pitch and roll motions. Despite its simplicity, it captures the most basic feature of the vertical model of the vehicle . We assume that wheels are rolling without slip or contact loss. Equations (1) and (2) represent the vehicle body and the wheel motion, respectively: where and are, respectively, the sprung mass and the unsprung mass of a quarter car, represents the suspension stiffness, represents the tire vertical stiffness, is the damper damping coefficient, is the sprung mass position, is the unsprung mass position, and is the road input.
2.2. Kalman Filter for Road Input Estimation
The road input is an estimated signal and it should be a part of the system states. Hence, in this paper, state variables for the quarter car are presented as follows:
And the road roughness profile applied in the estimation observer satisfies the following equation : where and are constant real numbers. The modified quarter car has to be stable in order to provide useful result; this implies that the real parts of the eigenvalues associated with (4) have to be negative, and consequently converges to zero based on the selection of the constants and . However, since represents the road profile input, it should contain as much road input information as possible. In other words, should be damped as slowly as possible, which will affect the stability of (4). Therefore, there is a trade-off between system stability and road input information integrity. In this paper, and are set to , which are close to two times of the quarter-car tire hop frequency, is sufficient to keep (4) stable . The quarter car state equations can be written as follows: where is a state vector; is an observer vector, where is a suspension deflection measured by sensor; is a vehicle body position calculated by a twice numerical integration (trapezoidal method) of the filtered vertical acceleration signal; is a filtered vertical acceleration. and are the process and measurement noise vector, respectively, assumed to be white, zero mean, and uncorrelated.
Evolution and observation constant matrices are given as where is the identity matrix and is the sampling period. And it is easy to verify that the observability matrix is full rank:
A standard Kalman filter formulation is used for the velocity estimation [31, 32]. The time update and measurement update equations of the filter are constructed as follows: where is a process noise covariance matrix, is a measurement noise covariance matrix, is an estimation of error covariance matrix, is a Kalman gain, and is an estimated state vector. The process and measurement noise variables and determine how much should the process model and the measurements should be trusted by the filter.
2.3. Simulation Setup
Two types of road inputs were used during experimental and simulation analysis. One input is a trigonal bump with 40 mm height and 400 mm length shown in Figure 3; another input is a random road input modeled based on the inverse Fourier transform as Figure 4.
Figure 5 to Figure 8 show a comparison between the actual and the estimated road profile height. In Figures 5 and 7, the vehicle passes the bump at a constant speed of 15 km/h and 45 km/h, respectively. In Figures 6 and 8, the vehicle runs on a random road, and the speed is 30 km/h and 60 km/h, respectively. The simulation results show that, in the random road input or the low speed bump input condition, the road estimator is very efficient, while in the condition of passing the bump with a relatively high speed, as shown in Figure 7, the estimation result is not very well. This is induced by the road model defined in (4). In fact, the road model applied in the proposed estimation observer should make a compromise between system stability and the integrity of road input information. We want to keep the system stable sacrificing the integrity of road input information in this research. And this strategy decreases the road estimator performance in an impulse input condition. The higher the vehicle speed is, the sharper the road input in time domain will be and the worse the road estimator will be. But its negative effect is acceptable since most drivers will reduce the vehicle speed when they run over a bump.
3. Sprung Mass Estimation
The road estimation method presented afore is based on a hypothesis that the vehicle system parameters in Table 1 are known and invariable. But in a real vehicle control system, it is not easy to obtain the sprung mass of the equivalent 2-DOF model, because of the nonlinear properties of the suspension system and the weight distribution variation on different driving conditions. What is more, the sprung mass varies greatly between empty and full loaded, especially for a commercial vehicle. All these factors will induce a significant sprung mass change and, then, attenuate the efficiency of the road estimation system. As shown in Figure 9, 20% sprung mass error will induce an obvious road estimation error. Consequently, it is necessary to estimate the sprung mass online for improving the efficiency of road estimation.
3.1. Solution Formulation
With setting , and , (2) can be written as
The vehicle sprung mass and suspension damping coefficient are estimated simultaneously via the Recursive Least-Squares Estimation. The Recursive Least-Squares Estimation procedures are given as follows: where is a symmetric covariance matrix and is the Least-Squares estimation of ; the calculation order is described as follows.
Calculate by and according to (11) first; then, update the value of and ; and then, calculate and according to (10) and (12), respectively; last, recursive the mentioned steps in order until the error reaches the standard. As for the initial values, and , they can be calculated by the Least-Squares Estimation.
Here, the vehicle system model parameters are from Table 1; the road profile defined in Figure 4 is selected as a vertical road input and the vehicle runs with a constant speed of 20 km/h. The estimation result is shown in Figure 10, which demonstrates the estimation result converges at an acceptable region (the error is less than 5%) after 2.5 s and then shows great efficiency after 20 s. Generally, the vehicle load is a constant during a vehicle start-stop period and the convergence time for sprung mass estimation is usually far less than a vehicle start-stop period. The sprung mass estimation is a little bit time-consuming, but it is not a big issue. For example, the original sprung mass is set in normal load condition. The road estimation block will be triggered once the vehicle has been started and the vehicle velocity is greater than 10 km/h; it will obtain the sprung mass (the error is less than 5%) in 2.5 s after the road estimation block is triggered, and the sprung mass will be refreshed. And then the sprung mass will be updated in every 2.5 s until it is stable (the change is less than 1%) or the estimation time is more than 40 s.
4. Road Categorizer
In the previous sections, the road estimator performs very well in a quarter-car system. However, does this estimator have the same performance applied to an actual car online system? It cannot be denied that the estimated result, in an actual car online system, will be worse than expectation due to the effect of wheel radius, wheel contact area, and the noise signals in the control system. What is more, the estimated road profile cannot be applied to the control strategy directly since these road sequences are estimated by postestimation, which produces an unavoidable time delay. However, in many vehicle control systems, such as active/semiactive suspension system, obtaining the primary road category information is effective . Hence, it is necessary to develop a road categorizer to distinguish the road category for the main control strategy. But the most common road classification methods, estimating the road PSD (power spectral density) or the road RMS (root mean square), are not enough effective for an online control system, because these methods need to compute plenty of road height values which will consume a lot of estimation time, and this time lag usually cannot be neglected in a real-time control system.
In fact, most onboard suspension control systems, only require the mainly frequency or amplitude information of road input. In [33, 34], a road-frequency adaptive suspension is proposed, where the road surface is classified by the frequency properties according to the fact that [0–4 Hz] is the car-body frequency region and [4–8 Hz] is the human-body frequency region. [8–12 Hz] is the wheel frequency region and [12– Hz] is the harshness frequency region. Bastow et al. also categorize the road profile into four grades, which defined a very good surface with amplitudes under 5 mm, medium-quality roads with amplitudes less than 13 mm, poor-quality roads with amplitudes less than 25 mm, and off-road with the amplitudes often exceeding 25 mm . However, it is a challenging issue to design a proper filter to avoid some unimportant long waves without any phase delay, as shown in Figure 4; the road waves are not fluctuating along the horizontal plane due to the effect of some long waves which has small contribution on the vehicle vertical dynamic; consequently, we should estimate the road input frequency and amplitude avoiding the unimportant long waves. In the other hand, the road categorizer only depending on the road frequency or the road amplitude does not work very well always. And the suspension control strategy should consider the road frequency and amplitude simultaneously in some case, such that a car runs on a road twice with different speed or a car runs on different road (ISO Class A and ISO Class C) with same speed. It will be much better if the road categorizer considers road frequency and amplitude simultaneously. Hence, a novel road categorizing method considering road frequency and amplitude simultaneously is proposed in the following paragraph.
For road frequency estimation, refer to the article , a first-order zero-crossing algorithm is applied, with the road velocity state being the input, to identify the frequency components of road disturbances. The frequency estimation result of Figure 11 is shown in Figure 12, and the current estimation value is actually the last circle’s frequency.
As for road amplitude analysis, road velocity state does not contain the amplitude information, but the road amplitude estimation only via is very difficult, since some of the significant road waves are not fluctuating along the horizontal plane. As shown in Figure 11, the value of “” is one of the main amplitudes values, which is the amplitude of the wave with dominant frequency, should be estimated.
In the sequel, a novel road amplitude estimation method is proposed as shown in Figure 11, where, , , and , used to detect each complete cycle, are the three attached times when first-order zero-crossing of the road velocity sequence happens; , , and are the three road height values in the times , , , respectively. Then the road amplitude “” can be simply calculated as (13), where , are used to detect the troughs of a circle and used to detect the peak of a circles are the distance between the point (, ) to the line going through the points (, ), (, ); is two times of the amplitude ; in this paper, an approximate algorithm is used as follows:
The amplitude estimation result is shown in Figure 13; even though the algorithm is running in real time, the estimation result still has a complete cycle time lag.
As road frequency and amplitude estimation can be achieved, a road categorizer special for a semiactive suspension control system can be designed, as shown in Figure 14. A certain type of road can be detected by judging the amplitude and frequency of the estimated road profile into the specific range (I–VI), where the suggestion values of , , and are 10 mm, 4 Hz, and 8 Hz, respectively. The suspension mode selection based on the road condition in specific range (I–VI) is detailed in Table 2. Ranges I and II are the medium-low frequency and low amplitude areas, corresponding to good road surface with medium-low speed, mostly on an urban driving condition. In this case, the suspension system should be turned to soft mode to improve the vehicle ride performance. Range III is the high frequency and low amplitude area, corresponding to the expressway driving condition, in which more attention should be paid to the tire deflection, and a relative hard mode suspension system should be regulated to achieve better road holding for driving safety. Range IV is the low frequency and high amplitude area corresponding to a bad road with low speed, which is very common for an off-road. In this condition, limiting the low frequency body motion should be a key objective. Ranges V and VI correspond to poor-quality road surface with medium-high speed and the impulse road input. In these cases, the suspension system should be adjusted to attenuate the impact feeling.
An acquisition of road disturbances property is essential for the enhancement of suspensions control systems. This paper presented a method to estimate the road profile elevation based on Kalman filter. To minimize the estimation errors, an online identification system based on Recursive Least-Squares Estimation is adopted to estimate sprung mass in real time, which is applied to refresh the system matrix of the adaptive observer. And a novel road categorizer considering road frequency and amplitude simultaneously is approached to classify various road profile sequence for suspension control system. The main conclusions are as follows.(1)A road profile estimator based on linear Kalman filter is proposed, which has great advantages on practical online vehicle control. (2)An online sprung mass estimator is proposed, which demonstrates the estimation result converges at an acceptable region (the error is less than 5%) after 2.5 s and then shows great efficiency after 20 s. With this online sprung mass estimator, the accuracy road estimation result can be improved greatly.(3)A novel road amplitude estimation method is proposed. And the road condition is categorized into six special ranges according to the road frequency and amplitude estimation result simultaneously, which can provide the suspension control system with a better trade-off for the ride comfort, handling, and safety performance.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Special thanks are due to the Open Research Fund Program of the State Key Laboratory of Advanced Design and Manufacturing for Vehicle Body (31115028), the China Postdoctoral Science Foundation (2012M520028), the National Natural Science Foundation of China (51205155), and the National Basic Research Program of China (973 Program) (2011CB711201) for supporting authors’ research.
- R. P. Barre, R. T. Forbes, and S. Andrews, “The measurement and analysis of road surface roughness,” Motor Industry Research Association Report, 1970.
- M. W. Sayers, T. D. Gillespie, and W. D. O. Paterson, “Guidelines for conducting and calibrating road roughness measurements,” World Bank Technical Paper Number 46, 1986.
- ISO 8608, Mechanical Vibration Road Surface Profiles Reporting of Measured Data, International Standards Organization, 1995.
- “Standard test Method for measuring road roughness by static rod and level method,” Annual Book of ASTM Standards, vol. 4, pp. 750–755, 1996.
- “Standard test method for measuring longitudinal profile of traveled surfaces with an accelerometer: established inertial profiling Reference,” Annual Book of ASTM Standards, no. 4, pp. 702–706, 1996.
- E. B. Splanger and W. J. Kelly, “Road profilometer method for measuring road profile,” General Motors Research Publication GMR-452, 1964.
- A. González, E. J. O'Brien, Y.-Y. Li, and K. Cashell, “The use of vehicle acceleration measurements to estimate road roughness,” Vehicle System Dynamics, vol. 46, no. 6, pp. 483–499, 2008.
- H. M. Ngwangwa, P. S. Heyns, F. J. J. Labuschagne, and G. K. Kululanga, “Reconstruction of road defects and road roughness classification using vehicle responses with artificial neural networks simulation,” Journal of Terramechanics, vol. 47, no. 2, pp. 97–111, 2010.
- R. McCann and S. Nguyen, “System identification for a model-based observer of a road roughness profiler,” in Proceedings of IEEE Region 5 Technical Conference (TPS '07), pp. 336–343, April 2007.
- D. J. Whitehouse and J. F. Archard, “The properties of random surfaces of significance in their contact,” Proceedings of the Royal Society of London. Series A, vol. 316, pp. 97–121, 1970.
- M. Shinozuka and C.-M. Jan, “Digital simulation of random processes and its applications,” Journal of Sound and Vibration, vol. 25, no. 1, pp. 111–128, 1972.
- D. Cebon, Interaction between Heavy Vehicles and Roads, Society of Automotive Engineers, 1993.
- J.-J. Wu, “Simulation of rough surfaces with FFT,” Tribology International, vol. 33, no. 1, pp. 47–58, 2000.
- T. Yoshimura, “A semi-active suspension of passenger cars using fuzzy reasoning and the field testing,” International Journal of Vehicle Design, vol. 19, no. 2, pp. 150–166, 1998.
- A. Pazooki, S. Rakheja, and D. Cao, “Modeling and validation of off-road vehicle ride dynamics,” Mechanical Systems and Signal Processing, vol. 28, pp. 679–695, 2012.
- R. Hassan and R. Evans, “Road roughness characteristics in car and truck wheel tracks,” International Journal of Pavement Engineering, pp. 1–10, 2012.
- I. Fialho and G. J. Balas, “Road adaptive active suspension design using linear parameter-varying gain-scheduling,” IEEE Transactions on Control Systems Technology, vol. 10, no. 1, pp. 43–54, 2002.
- H. Fang, Y. Shi, and J. Yi, “A new algorithm for simultaneous input and state estimation,” in Proceedings of the American Control Conference (ACC 08), pp. 2421–2426, Seattle, Wash, USA, June 2008.
- H. Fang, Y. Shi, and J. Yi, “On stable simultaneous input and state estimation for discrete-time linear systems,” International Journal of Adaptive Control and Signal Processing, vol. 25, no. 8, pp. 671–686, 2011.
- H. Dahmani, M. Chadli, A. Rabhi, and A. El Hajjaji, “Vehicle dynamics and road geometry estimation using a Takagi-Sugeno fuzzy observer with unknown inputs,” in Proceedings of IEEE Intelligent Vehicles Symposium (IV '11), pp. 272–277, Baden-Baden, Germany, June 2011.
- H. Zhang, Y. Shi, and A. Saadat Mehr, “Robust static output feedback control and remote PID design for networked motor systems,” IEEE Transactions on Industrial Electronics, vol. 58, no. 12, pp. 5396–5405, 2011.
- H. Zhang, Y. Shi, and A. S. Mehr, “Robust PID control for multivariable networked control systems with disturbance/noise attenuation,” International Journal of Robust and Nonlinear Control, vol. 22, no. 2, pp. 183–204, 2012.
- H. Zhang, Y. Shi, A. Saadat Mehr, and H. Huang, “Robust FIR equalization for time-varying communication channels with intermittent observations via an LMI approach,” Signal Processing, vol. 91, no. 7, pp. 1651–1658, 2011.
- K. H. Guo and W. H. Yu, “Semi-active suspension adaptive control strategy,” Journal of Hunan University (Natural Sciences), no. 40, pp. 39–44, 2013.
- X. J. Zhang and W. H. Yu, “Semi-active suspension adaptive control strategy based on hybrid control,” in FISITA World Automotive Congress, 2012.
- M. Doumiati, A. Victorino, A. Charara, and D. Lechner, “Estimation of road profile for vehicle dynamics motion: experimental validation,” in Proceedings of the American Control Conference (ACC '11), pp. 5237–5242, San Francisco, Calif, USA, July 2011.
- H. Imine, Y. Delanne, and N. K. M'Sirdi, “Road profile input estimation in vehicle dynamics simulation,” Vehicle System Dynamics, vol. 44, no. 4, pp. 285–303, 2006.
- H. Imine, Y. Delanne, and N. K. M'Sirdi, “Road profile inputs for evaluation of the loads on the wheels,” Vehicle System Dynamics, vol. 43, pp. 359–369, 2005.
- M. M. Elmadany and Z. S. Abduljabbar, “Linear quadratic Gaussian control of a quarter-car suspension,” Vehicle System Dynamics, vol. 32, no. 6, pp. 479–497, 1999.
- S. Nguyen, System Identification for a Model Based Observer Road Roughness Profiler, University of Arkansas, 2006.
- W. Greg and B. Gary, An Introduction to the Kalman Filter, Department of Computer Science, University of North California at Chapel Hill, 2006.
- R. E. Kalman, “A newapproach to linear filtering and prediction problems,” Transactions of the ASME-Journal of Basic Engineering, no. 82, pp. 35–45, 1960.
- J. Lu, “A frequency-adaptive multi-objective suspension control strategy,” Journal of Dynamic Systems, Measurement and Control, vol. 126, no. 3, pp. 700–707, 2004.
- L. H. Nguyen, K.-S. Hong, and S. Park, “Road-frequency adaptive control for semi-active suspension systems,” International Journal of Control, Automation and Systems, vol. 8, no. 5, pp. 1029–1038, 2010.
- D. Bastow, G. Howard, and J. P. Whitehead, Car Suspension and Handling, SAE International, Warrendale, Pa, USA, 4th edition, 2004.
Copyright © 2013 Wuhui Yu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.