File Exchange

image thumbnail

Animation of a Quarter-Car Automotive Suspension

version 1.0.0.1 (90.7 KB) by James Allison
A standard linear quarter-car model is simulated and then animated.

93 Downloads

Updated 14 Mar 2019

View License

This model was developed for use in teaching modeling, simulation, and optimization in graduate engineering courses. The simulation is particularly useful for demonstrating different simulation designs that are optimized for different applications (e.g., maximum handling capability, maximum comfort, etc.). The model is based on the electric vehicle case study found in this dissertation:

http://deepblue.lib.umich.edu/handle/2027.42/58449

This code was used in the creation of this corresponding video, used for teaching:

https://youtu.be/okiScBW7exQ

Note: this submission does not include code for the multi-objective optimization highlighted in the youtube video.

Cite As

James Allison (2019). Animation of a Quarter-Car Automotive Suspension (https://www.mathworks.com/matlabcentral/fileexchange/35478-animation-of-a-quarter-car-automotive-suspension), MATLAB Central File Exchange. Retrieved .

Comments and Ratings (35)

I use this example for assignments in some courses, so I have not made the optimization part of the code public. If you would like to see a simple example of multi-objective optimization, I would recommend:

https://www.mathworks.com/matlabcentral/fileexchange/64787-demonstration-of-two-multi-objective-optimization-strategies

For a more detailed treatment of suspension design optimization, please refer to:

https://www.mathworks.com/matlabcentral/fileexchange/40504-co-design-of-an-active-automotive-suspension-using-dt

The code in the following HCS design optimization example may be useful for learning how to use fmincon for design optimization problems:

https://www.mathworks.com/matlabcentral/fileexchange/64651-design-optimization-of-a-helical-compression-spring?s_tid=prof_contriblnk

Sir,
Thanks for your reply.
I have one more doubt.
From where can I get the the Optimization file (opt_qcar.m) which is shown in your video
Thanks

I checked the zip file from MLC FX, and it does contain the mat file with the data. I suspect that something is filtering out this mat file when you attempt to download it or open it.

I will also update the YouTube link.

sir,
Hello Mr, James Allison sir,
i am facing the same problem, it seems some thing missing in this part of the program,
Error using load
Unable to read file IRI_737b: No such file or directory.

Error in Untitled5 (line 43)
load IRI_737b % road profile data

plz can u halp me in details sir..

little

Mike Duke

Change the fresh in the

I have the same problem as ben and BS CE. The error code is:
??? Error: File: plotsusp.m Line: 42 Column: 3
Expression or statement is incorrect--possibly unbalanced (, {, or [.

Error in ==> sim_qcar at 79
plotsusp([z0(i), z1(i)*zmf, z2(i)*zmf, t(i)],road_x,road_z,x(i),umf);

Any idea how to solve the problem?

Thank you

What the unit of the z1 amd z2 vector(suspension travel)? The unit, according to the code appears to be in centimeters.

Sankar: You could change the road profile data by replacing the .mat file with one you generate, as long as the units, spatial resolution, and variable names are identical.

Lucas: If I remember correctly, the unit is meters. The profile was generated using filtered noise, as described in Section 8.3.3 of: http://systemdesign.illinois.edu/publications/All08b.pdf

What is the unit of the pavement profile? And by what method did you get the profile? Thank you

BS CE

Mr. Allison, I have found following while running the function as in video....

??? Error: File: plotsusp.m Line: 42 Column: 3
Expression or statement is incorrect--possibly unbalanced (, {, or [.

Error in ==> sim_qcar at 79
plotsusp([z0(i), z1(i)*zmf, z2(i)*zmf, t(i)],road_x,road_z,x(i),umf);

Hi.. Is there any way to change the road profile data..

Henry Piz

Amritam:

I would first point you to the following journal paper:
http://systemdesign.illinois.edu/publications/All14b.pdf

and then my dissertation:
http://systemdesign.illinois.edu/publications/All08b.pdf
(chapter on EV design)

Both discuss these metrics. If I remember correctly my dissertation cites at least one reference that compares the simple metrics I use in these models to much more sophisticated metrics (e.g., frequency-based constraints for comfort). The reference then concludes that the simple integral of acceleration metric is very accurate even though it is not very complex. There should be other references that address both rattle space and handling.

Regarding multi-objective optimization, I would suggest the following review article:

http://link.springer.com/article/10.1007%2Fs00158-003-0368-6

AMRITAM DAS

Hi, Prof.James Allison. You have mentioned that you can't upload the optimization code. But can you at least provide a document or slides which explain the theory behind defining maximum handling capability, maximum comfort and design a multi-objective optimization problem. I am doing my thesis on similar topic. I believe it will be of great help.

Thank you Patrick for pointing out a fix for newer versions of MATLAB.

Works on R2015b if you change 'refresh' to 'drawnow'

I've not tested this code on more recent versions of MATLAB. The animation worked well with R2011b, but perhaps the animation no longer works with more recent releases of MATLAB.

Xaris

I can only see the plots and the picture of the suspension. I can't see the kinematic vision

Naveen: IRI_737b is a .mat file included with this submission. You need to make sure it is on your path to run this code.

Hello Mr, James Allison sir,
it seems some thing missing in this part of the program,
Error using load
Unable to read file IRI_737b: No such file or directory.

Error in Untitled5 (line 43)
load IRI_737b % road profile data

plz can u halp me sir..

ben

Hello Mr, James Allison

it seems some thing missing in this part of the program, [~, ....

[~,istart] = min(abs(xstart-road_x));
xend = curr_x + fw;
[~,iend] = min(abs(xend-road_x));
xpstart = xstart-curr_x;

can you help please? because i cant run it .. it gives me this error

??? Error: File: plotsusp.m Line: 42 Column: 3
Expression or statement is incorrect--possibly unbalanced (, {, or [.

Saman: Thank you for your interest in the optimization code. I decided not to distribute it because I use it in optimization classes that I teach. You might be interested in this submission:

http://www.mathworks.com/matlabcentral/fileexchange/40504

It includes the optimization code, but is a more complex example.

To answer your other question, f1 is an interface function between qcar (the simulation) and fmincon.

Saman: Thank you for your *interest* in the optimization code. ...

Hi James
Could you please add (opt_qcar) or send it to me (saman.abdanan@gmail.com). I saw your simulation on youtube for optimization. Could you please explain what was f1 in your m.file line 17?

Alexandru: Perhaps there is another function 'ss' on your path that is superseding the built-in ss function:

http://www.mathworks.com/help/control/ref/ss.html

You can use the 'which' command to find out if this is the issue.

Another possibility is the ss is part of the control system toolbox. If you don't have that toolbox, then ss will not work.

Alexandru

Hello,

I am interested in your app. But I have problems when running it. I get this error:
Undefined function 'ss' for input arguments of type 'double'.

Error in sim_qcar (line 48)
qcar = ss(Aqcar,Bqcar,Cqcar,Dqcar);

What am I missing?
Thank you!
Alex

umf is a multiplication factor that allows quick exploration of road profiles with different levels of roughness. It should be set to 1 for the output to correspond to the road profile in the .mat file.

Nice catch Louis. Yes, dt should be changed to dt2 on line 62.

Louis

Very helpful model, I have two comments though:
1. I think umf should be set to 1 (or removed)
2. on line 62 the division by dt should be changed to a division by dt2

These two details made the script output inaccurate results; the script otherwise gives consistent results.

Mayur

James,

Thanks for this model. I am trying to understand what umf = 3 is, and why you have used it. Could you please explain.

Thanks,
Mayur

very nice !!!

Updates

1.0.0.1

Updated description

MATLAB Release Compatibility
Created with R2011b
Compatible with any release
Platform Compatibility
Windows macOS Linux