3/25/2020 Control Tutorials for MATLAB and Simulink - Suspension: Simulink Modeling

Tips Effects Search Control Tutorials

TIPS ABOUT BASICS HARDWARE INDEX NEXT ►

INTRODUCTION CRUISE CONTROL MOTOR SPEED MOTOR POSITION SUSPENSION INVERTED PENDULUM

SYSTEM

MODELING Suspension: Simulink Modeling

Related Tutorial Links
ANALYSIS
Simulink Basics

Modeling w/ Simulink
CONTROL Circuit Sim Activity

PID
Related External Links
ROOT LOCUS Simulink Intro Video

FREQUENCY Simulink Modeling

Video
STATE-SPACE
Modeling Challenges
DIGITAL Video

SIMULINK

MODELING
Contents
CONTROL
Physical setup
SIMSCAPE Design requirements

Building the Model

Open-loop response

ctms.engin.umich.edu/CTMS/index.php?example=Suspension&section=SimulinkModeling 1/13
3/25/2020 Control Tutorials for MATLAB and Simulink - Suspension: Simulink Modeling

Physical setup

Designing an automotive suspension system is an interesting and challenging control problem. When the suspension system is designed,
a 1/4 model (one of the four wheels) is used to simplify the problem to a 1-D multiple spring-damper system. A diagram of this system is
shown below.

The system parameters are as follows.

(m1) body mass 2500 kg

(m2) suspension mass 320 kg

(k1) spring constant of suspension system 80,000 N/m

ctms.engin.umich.edu/CTMS/index.php?example=Suspension&section=SimulinkModeling 2/13
3/25/2020 Control Tutorials for MATLAB and Simulink - Suspension: Simulink Modeling

(k2) spring constant of wheel and tire 500,000 N/m

(b1) damping constant of suspension system 350 N.s/m

(b2) damping constant of wheel and tire 15,020 N.s/m

(u) control force = force from the controller we are going to design

Design requirements

A good automotive suspension system should have satisfactory road holding ability, while still providing comfort when riding over bumps
and holes in the road. When the vehicle is experiencing any road disturbance (i.e. pot holes, cracks, and uneven pavement),the vehicle
body should not have large oscillations, and the oscillations should dissipate quickly. Since the distance X1-W is very difficult to measure,

and the deformation of the tire (X2-W) is negligible, we will use the distance X1-X2 instead of X1-W as the output in our problem. Keep in
mind that this is an estimation.

The road disturbance (W) in this problem will be simulated by a step input. This step could represent the vehicle coming out of a pothole.
We want to design a feedback controller so that the output (X1-X2) has an overshoot less than 5% and a settling time shorter than 5

seconds. For example, when the vehicle runs onto a 10 cm high step, the vehicle body will oscillate within a range of +/- 5 mm and return
to a smooth ride within 5 seconds.

Building the Model

This system will be modeled by summing the forces acting on both masses (body and suspension) and integrating the accelerations of
each mass twice to give velocities and positions. Newton's law will be applied to each mass. Open Simulink and open a new model
window. First, we will model the integrals of the accelerations of the masses.

(1)

(2)

ctms.engin.umich.edu/CTMS/index.php?example=Suspension&section=SimulinkModeling 3/13
3/25/2020 Control Tutorials for MATLAB and Simulink - Suspension: Simulink Modeling

Insert an Integrator block (from the Continuous library) and draw lines to and from its input and output terminals.

Label the input line "a1" (for acceleration) and the output line "v1" (for velocity) To add such a label, double click in the empty space
just above the line.

Insert another Integrator block and connect it to the output of the first.

Draw a line from its output and label it "x1" (for position).

Insert a second pair of Integrators below the first with lines labeled "a2", "v2", and "x2".

Next, we will start to model Newton's law. Newton's law for each of these masses can be expressed as:

(3)

(4)

ctms.engin.umich.edu/CTMS/index.php?example=Suspension&section=SimulinkModeling 4/13
3/25/2020 Control Tutorials for MATLAB and Simulink - Suspension: Simulink Modeling

These equations can be represented with gain blocks (for 1/M1 and 1/M2) and two summation blocks.

Insert two Gain blocks, (from the Math Operations library) one attached to the inputs of each of the integrator pairs.

Edit the gain block corresponding to M1 by double-clicking it and changing its value to "1/m1".

Change the label of this Gain block to "Mass 1" by clicking on the word "Gain" underneath the block.

Similarly, edit the other Gain's value to "1/m2" and it's label to "Mass 2". (You may want to resize the gain blocks to view the contents.

To do this, single click on the block to highlight it, and drag one of the corners to the desired size.)

There are three forces acting on M1 (one spring, one damper, and the input, u) and five forces acting on M2 (two springs, two dampers,
and the input, u).

Insert two Add blocks (from the Math Operations library), one attached by a line to each of the Gain blocks.

Edit the signs of the Add block corresponding to M1 to "+--" to represent the three forces (two of which will be negative).

Edit the signs of the other Add block to "++-++" to represent the five forces, one of which will be negative.

ctms.engin.umich.edu/CTMS/index.php?example=Suspension&section=SimulinkModeling 5/13
3/25/2020 Control Tutorials for MATLAB and Simulink - Suspension: Simulink Modeling

Now, we will add in the forces acting on each mass. First, we will add in the force from Spring 1. This force is equal to a constant, k1 times

the difference X1-X2.

Insert an Add block after the upper pair of integrators.

Edit its signs to "+-" and connect the "x1" signal to the positive input and the "x2" signal to the negative input.

Draw a line leading from the output of the Add block.

Insert a Gain block above the "Mass 1" block.

Flip it left-to-right by single-clicking on it and selecting Flip Block from the Rotate & Flip menu (or hit Ctrl-I).

Edit the value of this gain to "k1" and label the block "Spring 1".

Tap a line off the output of the last Add block and connect it to the input of this Gain block.

Connect the output of this Gain block (the spring force) to the second input of the Mass 1 Add block. This input should be negative
since the Spring 1 pulls down on Mass 1 when X1 > X2.

Tap a line off the spring force line and connect it to the second input of the Mass 2 Add block. This input is positive since Spring 1

pulls up on Mass 2.

ctms.engin.umich.edu/CTMS/index.php?example=Suspension&section=SimulinkModeling 6/13
3/25/2020 Control Tutorials for MATLAB and Simulink - Suspension: Simulink Modeling

Now, we will add in the force from Damper 1. This force is equal to b1 times V1-V2.

Insert an Add block below the Mass 1's first integrator.

Flip it left-to-right, and edit its signs to "+-".

Tap a line off the "v1" line and connect it to the positive input of this Add block.

Tap a line off the "v2" line and connect it to the negative input of this Add block.

Insert a Gain block to the left of this Add block and flip it left-to-right.

Edit its value to "b1" and label it "Damper 1".

Connect the output of the new Add block to the input of this gain block.

Connect the output of this gain block (the damper force) to the third input of the Mass 1 Add block. This input is negative, similar to

Spring 1's force on Mass 1.

Tap a line off Damper 1's force line and connect it to the first input (which is positive) of Mass 2's Add block.

ctms.engin.umich.edu/CTMS/index.php?example=Suspension&section=SimulinkModeling 7/13
3/25/2020 Control Tutorials for MATLAB and Simulink - Suspension: Simulink Modeling

Now we will add in the force from Spring 2. This force acts only on Mass 2, but depends on the ground profile, W. Spring 2's force is equal

to X2-W.

Insert a Step block in the lower left area of your model window. Label it "W".

Edit it's Step Time to "0" and it's Final Value to "0". (We will assume a flat road surface for now).

Insert an Add block to the right of the W Step block and edit its signs to "-+".

Connect the output of the Step block to the positive input of this Add block.

Tap a line off the "x2" signal and connect it to the negative input of the new Add block.

Insert a Gain block to the right of this Add block and connect the Add's output to the new Gain's input.

Change the value of the gain to "k2" and label it "Spring 2".

Connect the output of this block (Spring 2's force) to the fourth input of Mass 2's Add block. This force adds in in the positive sense.

ctms.engin.umich.edu/CTMS/index.php?example=Suspension&section=SimulinkModeling 8/13
3/25/2020 Control Tutorials for MATLAB and Simulink - Suspension: Simulink Modeling

Next, we will add in the force from Damper 2. This force is equal to b2 times V2-d/dt(W). Since there is no existing signal representing the

derivative of W we will need to generate this signal.

Insert a Derivative block (from the Continuous library) to the right of the W step block.

Tap a line of the Step's output and connect it to the input of the Derivative block.

Insert an Add block after the Derivative block and edit it's signs to "+-".

Connect the Derivative's output to the positive input of the new Add block.

Tap a line off the "v2" line and connect it to the negative input of this Add block.

Connect the output of this Add block (Damper 2's force) to the fifth input of Mass 2's Add block. This force also adds in with positive

sign.

ctms.engin.umich.edu/CTMS/index.php?example=Suspension&section=SimulinkModeling 9/13
3/25/2020 Control Tutorials for MATLAB and Simulink - Suspension: Simulink Modeling

The last force is the input U acting between the two masses.

Insert a Step block (from the Sources library) in the upper left of the model window.

Connect its output to the remaining input of Mass 1's Add block (with positive sign).

Tap a line off this signal and connect it to the remaining input of Mass 2's Add block (with negative sign).

Edit this Step block's Step Time to "0" and leave its Final Value "1".

Label this Step block "U".

Finally, to view the output (X1-X2) insert a Scope block (from the Sinks library) and connect it to the output of the rightmost Add block.

ctms.engin.umich.edu/CTMS/index.php?example=Suspension&section=SimulinkModeling 10/13
3/25/2020 Control Tutorials for MATLAB and Simulink - Suspension: Simulink Modeling

You can download a model file for the complete system by right clicking here and then selecting Save link as ....

Open-loop response

To simulate this system, first, an appropriate simulation time must be set. Select Model Configuration Parameters from the Simulation

menu and enter "50" in the Stop Time field. 50 seconds is long enough to view the open-loop response. The physical parameters must

now be set. Run the following commands at the MATLAB prompt:

m1 = 2500;

m2 = 320;

k1 = 80000;

k2 = 500000;

ctms.engin.umich.edu/CTMS/index.php?example=Suspension&section=SimulinkModeling 11/13
3/25/2020 Control Tutorials for MATLAB and Simulink - Suspension: Simulink Modeling

b1 = 350;

b2 = 15020;

Run the simulation (Ctrl-T or Run from the Simulation menu). When the simulation is finished, if you double-click on the Scope block you

should see the following output.

Observing the above, we would like to improve the response of the suspension through the control of the suspension control force U. The

model created here will be employed for controller design and analysis within Simulink in the Suspension: Simulink Controller Design
page.

Published with MATLAB® 9.2

All contents licensed under a Creative Commons Attribution-ShareAlike 4.0 International License.

ctms.engin.umich.edu/CTMS/index.php?example=Suspension&section=SimulinkModeling 12/13
3/25/2020 Control Tutorials for MATLAB and Simulink - Suspension: Simulink Modeling

ctms.engin.umich.edu/CTMS/index.php?example=Suspension&section=SimulinkModeling 13/13