Anlayesha Paper
Anlayesha Paper
Suspension and vibration analysis on half car model with and without active controller
Nemmadi Anlaysha , Amit Kumar Mehar
Raghu Engineering College, Dakamarri, Bheemunipatnam, Visakhapatnam.
energy and do not require external power. Semi-active suspensions incorporate
A MATLAB Simulink environment is utilized to numerically simulate a half signal or power supply. Full-active suspensions employ actuators, typically
car model with independent front and rear vertical passive suspensions. The hydraulic cylinders, to generate desired forces, necessitating external power for
purpose of this simulation is to analyze the reaction forces exerted by the front operation.
and rear wheels in response to the pitch and bounce degrees of freedom of the The design and development of automotive suspension systems have been of
car. By implementing a conventional passive suspension between the car body great interest for nearly a century. Early systems were adapted from horseless
and wheel assembly, consisting of a spring and a damper, the model aims to carriage practices. With the increase in vehicle speeds, complex vibration
provide insights into the ride characteristics. The chosen spring-damper problems affecting ride comfort and safety have emerged. Solving these
characteristics of the suspension are pre-selected with the intention of problems involves reducing excitation levels from road surface irregularities
emphasizing specific objectives, such as passenger comfort or suspension and designing effective suspension systems to maintain acceptable comfort
deflection. These characteristics play a vital role in determining the behavior of levels and ensure safety on various road conditions. Scientific approaches to
the suspension system. The simulation considers various conflicting objectives, suspension problems have gained momentum over time.
allowing for a comprehensive evaluation of the model's ride characteristics. The primary objective of a suspension system is to provide a high level of ride
Overall, the numerical simulation within the MATLAB Simulink environment quality and protect the vehicle structure by effectively isolating it from road
offers valuable insights into the behavior and ride characteristics of the half car surface irregularities. This necessitates a soft suspension. Additionally, the
model, shedding light on the reaction forces, suspension performance. suspension should ensure lateral stability and controllability under different
Keywords: half car model, vibration, analytical and numerical analysis. running conditions, support variable static loads, and maintain a stiff
Introduction suspension.
The car suspension system is responsible for separating the car body from its Despite significant theoretical and practical advancements in suspension
wheels. As vehicle capabilities have increased, there has been a significant systems, most current road vehicles still face the fundamental conflict of
improvement in the performance of suspension systems. To achieve an optimal compromising suspension parameters. Further improvements using
suspension system, several performance characteristics must be taken into conventional elements alone appear challenging. However, active and semi-
account. active suspension systems offer new possibilities for enhancing vehicle
These characteristics primarily deal with regulating body and suspension behavior, although their adoption in the road vehicle industry remains cautious
movement, as well as distributing forces. Ideally, the suspension should isolate due to cost and complexity concerns. Advances in computer technology and
the car body from road disturbances, as well as inertial disturbances during mathematical analysis provide opportunities for theoretical analysis to support
cornering, braking, or acceleration. Additionally, it should minimize the practical system design. Utilizing a simple model with few parameters, as long
transmission of vertical forces to ensure passenger comfort, achieved by as it reasonably represents the real system, can provide economical and efficient
reducing vertical car body acceleration. Excessive wheel travel can lead to guidelines for system design.
suboptimal tire-to-road contact, resulting in poor handling and adhesion. In the figure, a general representation of the car suspension system is depicted.
Therefore, maintaining the optimal tire-to-road contact on all four wheels is By employing accurate representations of vehicle dynamics, guide way inputs,
crucial for good handling characteristics. and appropriate performance criteria, these facilities enable the prediction and
In traditional suspension systems, these characteristics often conflict and fail to evaluation of suspension system performance.
0
⎡ ⎤ 0
⎡ 0 ⎤
⎢ ⎥
⎢ 0 ⎥ ⎢ 0 ⎥
( ) ⎢ 0 ⎥
⎢ ⎥
⎢ 0 ⎥
𝐵=⎢ 0 ⎥ 𝐶=⎢ ⎥
⎢ ⎥ 0
⎢ ⎥
⎢ ⎥
⎢ ⎥
⎢ 0 ⎥
Fig 2: car model. ⎢ ⎥
⎢ ⎥
⎣ ⎦
⎣ ⎦
quality, reduced vibrations, and enhanced vehicle stability. Continuous research solution and under which the steady state closed loop system is stable are
and development in this area contribute to the ongoing advancements in active The pair (𝐴, 𝐵) defines a stabilized system
suspension control strategies. The pair (𝐴, 𝐶) defines a detectable system
Linear Quadratic Regulator Controller Design of Active Suspension System for the 4-DOF vehicle model
Linear quadratic regulator control (LQR) is a powerful concept of optimally Let us consider the state space representation of linear time invariant system of
controlling linear systems whose states are available to be fed back, i.e. the half car (4DOF) model. The governing equations of motion for the half car
the states can be measured online. In this section mathematics of LQR model written in the following state space form
controller was discussed for the design of active suspension system to improve ẋ = Ax + Bu + Dr
the road handling and ride comfort for the full car and half car, quarter car Where A, B and D are state spaces matrices which are formulated in the
model. The LQR approach of vehicle suspension control is widely used in the mathematical modeling. The control algorithm adopted for the 2DOF vehicle
studies of controller design for the active suspension system. The LQR same as for 4DOF model but only the difference is state space representation
optimum control provides flexibility that that factors of the performance index changes and gain matrix changes. The selection of weighting parameters should
can be weighted according to the designer’s goal and other any constraints. The be taken to minimize the sprung body acceleration, pitch acceleration and
theory of optimal control is concerned with operating a dynamic system at suspension deflection. The following Q and R matrices chosen to find out the
minimum cost and the case where the system dynamics are described by a set control force (𝑡).
of linear differential equations and the cost is described by a quadraticfunction 1 0 0 0 0 0 0 0 0
⎡0 1 0 0 0 0 0 0 0⎤
is called the LQ problem. The LQR control design based on the automotive ⎢0
⎢ 0 1 0 0 0 0 0 ⎥⎥
performance criteria specified for the full car, half car and quarter car models 𝑄 = 1 × 10 ⎢0 0 0 1 0 0 0 0 ⎥ 𝐴𝑛𝑑 𝑅 =
⎢0 0 0 0 0 0 0 0 ⎥
will demonstrate the effect of LQR control to passenger cars with active ⎢0 0 0 0 0 0 0 0 ⎥
suspension.Fig 6.1 represents the application of LQR controller the plant to ⎢0 0 0 0 0 0 0 0 ⎥
⎣0 0 0 0 0 0 0 0 ⎦
estimate the optimal control force u to minimize the performance index.
1 0
Linear quadratic regulator (LQR) is an optimal control problem made the 0 1
following assumptions:
The state equation of plant is linear
The cost function is quadratic
Consider time linear time invariant system
𝑥̇ (𝑡) = 𝐴𝑥(𝑡) + 𝐵𝑢(𝑡) + 𝐷𝑟
The optimal solution of the linear regulator problem is to determine the control
law
𝑢(𝑡) = −𝐾𝑥(𝑡)
When substituting in the equation then the new system state becomes
𝑥̇ (𝑡) = (𝐴 − 𝐵𝐾)𝑥(𝑡) + 𝐷𝑟(𝑡) Fig 7: simulink model for LQR controller.
The procedure of optimization consists of finding control input U, In the development of a Simulink model for an active suspension system, the
which minimizes the performance index 𝐽 constitutes the performance specifications of the vehicle, in this case, the Toyota Avalon XLE 2015, are
characteristics requirement as well as the controller input limitation. The incorporated. These specifications are typically provided in the form of
optimal controller of the given system is defined as to find a control input 𝑢(𝑡) variables and are uploaded to the workspace of MATLAB from an M-file
minimizing the following quadratic performance index. editor.
Once the specifications are available in the MATLAB workspace, they can be
𝐽= (𝑥(𝑡) 𝑄 𝑥(𝑡) + 𝑢(𝑡) 𝑅𝑢(𝑡))𝑑𝑡
utilized within the Simulink model. The Simulink model is designed to simulate
Where 𝑢(𝑡) is state input vector, 𝑥(𝑡) is the state variable vector and 𝑄 and 𝑅 the behavior of the active suspension system based on the provided
are positive definite weighting matrices. The state feedback gain matrix 𝐾 specifications and control algorithms.
(controller) given by During the simulation, the variables representing the vehicle specifications,
𝐾=𝑅 𝐵 𝑃 such as mass, spring constants, damping coefficients, and other relevant
And P is the unique non negative definite solution of the algebraic Riccati parameters, are used in the model to define the characteristics of the suspension
equation (ARE) system.
𝐴 𝑃 + 𝑃𝐴 − 𝑃𝐵𝑅 𝐵 𝑃 + 𝑄 = 0 The simulation results are then analyzed and visualized using various tools
The necessary conditions under which riccati equation has a steady-state available in Simulink. One commonly used tool is the Scope block, which
International Journal Of Innovative Trends In Multi Disciplinary Scientific Research(Ijtmsr)
ISSN (Online): 2451 – 5632
(www.ijtmsr.com) Volume 8, Issue 4, 2023
allows for the observation and plotting of relevant signals and variables during
0.0035
the simulation. The Scope block provides a graphical representation of the 0.003
system's response, enabling engineers to evaluate the performance of the active 0.0025
suspension system based on the specified specifications. 0.002
By incorporating the Toyota Avalon XLE 2015 specifications into the Simulink 0.0015 30 kmph
0.001
model and visualizing the simulation results using tools like the Scope block, 40 kmph
0.0005
engineers can assess the effectiveness of the active suspension system and make 0 45 kmph
necessary adjustments or improvements to optimize the vehicle's ride comfort -0.0005 0 5 10 15
and handling characteristics. -0.001
-0.0015
0.006 -0.002
0.005
0.004 Graph 4: vehicle pitch at different velocities.
0.003 30 kmph The graph above depicts the pitch of a vehicle at various speeds. The preceding
0.002 graph clearly shows that if the vehicle goes at the slowest possible speed, it will
40 kmph
0.001
45 kmph have the greatest displacement.
0
-0.001 0 5 10 15 0.02
-0.002
-0.003 0.015
0.01
Graph 1: vehicle displacement at different velocities.
passive
The graph above depicts vehicle displacement as a function of vehicle velocity. 0.005
active
The preceding graph clearly shows that if the vehicle goes at the slowest
0
possible speed, it will have the greatest displacement. 0 5 10 15
-0.005
0.05
0.04 -0.01
0.03 30 kmph
Graph 5: vehicle displacement.
0.02 40 kmph
0.01 The graph above depicts the vehicle displacement suspension systems, both
45 kmph
0 passive and active. External force is exerted to the vehicle's suspension stem in
-0.01 0 5 10 15 this case. Then we may minimise the quantity of waves and enhance the
vehicle's damping. As shown in the preceding shot.
Graph 2: front wheel displacement at different velocities.
The graph above depicts the front wheel displacement as the vehicle's velocity 0.05
changes. The preceding graph clearly shows that if the vehicle goes at the
0.04
slowest possible speed, it will have the greatest displacement.
0.03
0.0001 passive
0.02
active
0.00005 30 kmph 0.01
40 kmph
0
0 45 kmph 0 5 10 15
0 5 10 15 -0.01
Graph 3: Rear wheel displacement at different velocities. The graph above depicts the passive and active suspension systems of the front
The graph above depicts the rear wheel displacement as the vehicle's velocity wheel displacement. External force is being thrust to the vehicle's suspension
changes. The preceding graph clearly shows that if the vehicle goes at the system in this case. Then we may minimise the quantity of waves and enhance
slowest possible speed, it will have the greatest displacement. the vehicle's damping. As shown in the preceding shot.
International Journal Of Innovative Trends In Multi Disciplinary Scientific Research(Ijtmsr)
ISSN (Online): 2451 – 5632
(www.ijtmsr.com) Volume 8, Issue 4, 2023
-0.004
References
-0.006
1.Shpetim Lajqi and StanislavPehan, “Designs and Optimizations of Active and
Graph 7: vehicle pitch behavior. Semi-Active Non-linear Suspension Systems for a Terrain Vehicle”, Journal of
The graph above depicts the passive and active car pitch suspension systems. Mechanical Engineering 58(2012)12, 732-743.
External force is applied to the vehicle's suspension system in this case. Then 2. Anirudh R. Iyer, R. Venkatachalam and A. Balaraju, “Analysis of Optimum
we may minimise the quantity of waves and enhance the vehicle's damping. As Suspension Parameters of a Semi-Independently Suspended Automobile”,
shown in the preceding shot. International Journal of Innovative Research in Science, Engineering and
natural frequency matlab ansys Technology, Vol. 3, Issue 5, May 2014, ISSN: 2319-8753. Bamankar,
mass 0.965 0.9927 International Journal of Advanced Engineering Research and Studies E-
ISSN2249–8974 Int. J. Adv. Engg. Res. Studies/V/I/Oct.-Dec,2015/01-03
inertia 3.24 1.955
3. Florin Mailat, Stefania Donescuand VeturiaChiroiu, “On the Automotive
frone wheel 11.373 11.31
Semi-Active Suspensions”, Proceedings of The Romanian Academy, Series A,
rear wheel 18.92 19.08
Volume 5, Number 1/2004.
Table 2: natural frequency of half car model.
4. S. Prabhakar, Dr.K.Arunachalam, “Simulation And Analysis Of Passive
25 Suspension System For Different Road Profiles With Variable Damping And
Stiffness Parameters”, National Conference On Recent Trends And
20
Developments In Sustainable Green Technologies,JCHPS Special Issue 7:
15 2015,ISSN: 0974-2115.
matlab
5. CatalinAlexandru, and PetreAlexandru, “A comparative analysis between the
10 ansys
vehicles passive and active suspensions”, International Journal of Mechanics,
5 Issue 4, Volume 5, 2011.