Modelling Car Suspension with ODE's: Damped Free Oscillations Part 1
In my previous post, I stated that all real world free oscillating systems are damped because their vibrations eventually attenuate away, due to energy dissipation to the surrounding environment.
However, for many engineering applications, the natural damping may not be enough. One way to damp oscillations quicker is to have the vibrating mass move through a viscous fluid, thus the energy of the motion is dissipated as heat through the thick fluid (as it causes the fluid to move).
Indeed, many real-life engineering applications use the concept of viscous damping, such as shock absorbers on mountain bike, motor cycle and car suspension systems.
![256px-Lotus_12,_rear_suspension.jpg](https://steemitimages.com/DQmNtetYtuS2HsDg7SkzxJUBUh2PWzFxmqhcytepQV1onFk/256px-Lotus_12%2C_rear_suspension.jpg)
Figure 1. Rear suspension on a historic Lotus Formula 1 car.
Source: Wikimedia Commons
In Figure 1, we can see that the engineers have coiled the spring over a shock absorber which is filled with oil. This is a common arrangement in many modern cars to make mounting this type of suspension easier.
A schematic of a mass-spring mechanism undergoing viscous damping is shown in Figure 2...
![r18.png](https://steemitimages.com/640x0/https://steemitimages.com/DQmamxpmXXVJ2qmckto7wzJ7UohpQgkEFEkRfZg3srkZY6D/r18.png)
Figure 2. Schematic of the mass-spring-damper system. Adapted from Calculus 5th Edition, James Stewart.
Like the spring, the damping force is always opposing the motion of the mass. A good approximation (confirmed by physical experiments) is the damping force is linearly proportional to the velocity of the mass, hence .
Thus with a damped system, the mass now has 2 opposing forces. Firstly...
![q3.png](https://steemitimages.com/640x0/https://steemitimages.com/DQmQRTWYwRJkQpNkq1tzKzY5jnEE2k7czMNkevqipZToiKV/q3.png)
And secondly...
![r1.png](https://steemitimages.com/DQmUa2MYUzD8QpJe4vJKpyX5SwLjDMgCxfxCfteBJSkmbRi/r1.png)
These forces are shown in blue, in Figure 2. Again, by Newton's 2nd Law, the resultant force on an object is equal to its mass by its acceleration...
![r2.png](https://steemitimages.com/640x0/https://steemitimages.com/DQmXBKVuJn8L2c31mJFtHxYrYGMUh9vNviB6GRahFbedEVM/r2.png)
Expressing acceleration as the second time derivative of displacement, we can rearrange equation (1) into a Second Order, homogeneous, linear ODE with constant coefficients...
![r3.png](https://steemitimages.com/DQmQ5Ue9io2fc8LMy4DYqLZvFHuYS6yA2xmdu5iTdLGt7sz/r3.png)
The characteristic equation of (2) is...
![r4.png](https://steemitimages.com/640x0/https://steemitimages.com/DQmR8c2T9X55Uff1iPH7bHP18enH6nPGWCj6agQSb1RAVaC/r4.png)
...and the roots of (3) are...
![r5.png](https://steemitimages.com/DQmddYcqFkKJp66Xkyx1mPsp95fG2XPwv58oeQdAANJvxNu/r5.png)
Now let the discriminant Δ equal the term inside the square-root and we have the same situation as that of Post #3 in this series: Roots of Homogeneous Second Order ODE's and the Nature their Solutions, where the behaviour of the solution will depend on the value of Δ.
So with...
![r6.png](https://steemitimages.com/640x0/https://steemitimages.com/DQmXFAooK6yFwLMiWppoqY3Ts1xAJ8nED4oUsRb6BMkd37e/r6.png)
...we can have 3 different cases...
Case I | ![]() |
![]() |
Under damping |
Case II | ![]() |
![]() |
Critical damping |
Case III | ![]() |
![]() |
Over damping |
Case I: Under damping
Suppose that we are investigating the suitability of a new car suspension system. If we focus our attention on just one corner of the car, we can simplify the suspension to our mathematical mass-spring-damper model.
Let's assign the system system with the following parameters...
- mass m = 300kg (total car body mass of 1200kg)
- spring constant k = 1200N/m (note, this is a very soft spring compared to real world applications)
- damping constant c = 600Ns/m
- initial displacement
= 0.15m (let's say this wheel has been driven off a curb of 15cm height, and therefore the spring is stretched by this amount as our starting position)
- initial velocity
= 0m/s
A schematic of this system is shown in Figure 3 below...
![r32.png](https://steemitimages.com/DQmdc2af6jwaAR1ZXckDQnYRJAqXfnuXPcpL1CoVhP11ZJq/r32.png)
Figure 3. Schematic of a simplified car suspension
(adapted from Wikimedia Commons)
Substituting this values for m, c and k in to equation (2), we have...
![r10.png](https://steemitimages.com/640x0/https://steemitimages.com/DQmXATRL3mPWjnw7AcPG7N9GqXmKHvZ3up7Txg6Ey7TaTWY/r10.png)
Characteristic equation for (4) is...
![r11.png](https://steemitimages.com/DQmUJxpNqk2ELwHoiNXjXWUUdEA91tNiBzQYWxpMCYQvCpQ/r11.png)
The roots of (5) are...
![r12.png](https://steemitimages.com/640x0/https://steemitimages.com/DQmaSHLCv3qmbPBiiD1GhTadu8vzi7jF8Jpor7mxRYNTzde/r12.png)
As we have complex conjugate roots, the general solution to equation (4), according equation (11) of Post #3 is...
![r13.png](https://steemitimages.com/DQmQMhnXRSbgVSihsgy25pGRXzKEzrWFDkd2nXsQwGTJX4x/r13.png)
Let's find the particular solution of (6) by applying the initial conditions. First, applying the initial displacement condition...
![r14.png](https://steemitimages.com/640x0/https://steemitimages.com/DQmR8y2wrgZ2v31xbFj8iFc1CcR1kA8dojwFRFEpL3sSrsA/r14.png)
The first derivative of (6) is...
![r15.png](https://steemitimages.com/DQmbYaa61XbLYo1LximxawEhqmu1tRCyGKsHbRtLNrhn6r4/r15.png)
Now we can apply the initial velocity to equation (7)...
![r16.png](https://steemitimages.com/640x0/https://steemitimages.com/DQmVsc2UsskpN1QmxeD1rV8A84XpbiPrLTxjU2icLXPTNL5/r16.png)
Therefore, the particular solution is...
![r17.png](https://steemitimages.com/DQmVC9hZZ9qXRPyQWEpo77NVCK62rLsA8heLUjYxwaEXQe2/r17.png)
The "amplitude phase-shift" equivalent solution is...
![r18.png](https://steemitimages.com/640x0/https://steemitimages.com/DQmVWZp2MEsnZG1gsQuCovpe4oGrSMRqyCF6Ucj9Nf4nAu6/r18.png)
The behaviour of the solution of equation (9) shown in Figure 3 below...
![r21.png](https://steemitimages.com/DQmY8e16U3Tjx1hzRTkmxRz8Bv3rN8fv57VeZEJHmpU8pPR/r21.png)
Figure 3. Free oscillations with damping coefficient c = 600Ns/m
As you can see, the motion attenuates quickly after 2 oscillations (including the initial displacement), with the peak of the second oscillation a mere fraction of the amplitude of the first. Note, although the oscillations are damped out quickly, this is still a case of under damping according to the mathematics.
So what happens when we lower the damping coefficient (or tune the damper down) to say, c = 150Ns/m?
Equation (2) becomes...
![r26.png](https://steemitimages.com/640x0/https://steemitimages.com/DQmSYUp7xEUMvGBz7vVQVLxSdHuo6atvHaoGJvRLG1oj994/r26.png)
And the characteristic equation for (10) is...
![r27.png](https://steemitimages.com/DQmW6XiNGohb7dC4CY8DiTSfHiiwuGn4NKE2vYSGL3qte9b/r27.png)
The roots of (11) are...
![r28.png](https://steemitimages.com/640x0/https://steemitimages.com/DQmZk2EHStvMPz1u61gZdEH2g88R1EPqDVFFs8QJaPiSF1A/r28.png)
As previously, from these roots, the general solution is...
![r29.png](https://steemitimages.com/DQmdZEz2iRBJpNC6XJBaL1tat4RfYrpCtgYfbSAWjGkgxd5/r29.png)
And as previously, applying both initial conditions, the particular solution is...
![r30.png](https://steemitimages.com/640x0/https://steemitimages.com/DQmV62ngGr4mvzbpsQqfo2Zym4ezVCoBmzTwpYuRUwxrTT4/r30.png)
The "amplitude phase-shift" equivalent solution...
![r31.png](https://steemitimages.com/DQmVNRNLsSqbuK3N6qiYA5S4yM82TW5eH6mrsfxXntuwqpX/r31.png)
The behaviour of (14) shown in Figure 4 below...
![r25.png](https://steemitimages.com/640x0/https://steemitimages.com/DQmeERpmdSdHkdLfC7a1wcwPBLkHD7qHX4VXh2g7wt9pMZ9/r25.png)
Figure 4. Free oscillations with damping coefficient c = 150Ns/m
Now there are at least 6 oscillations before the motion peters out. Do you think this is a good performance from a car's suspension? This car's going to bump up and down 6 times with every pothole it drives over!
Take a look at the video below (published by Open University) which excellently explains what happens to a car when it is severely underdamped. Enjoy.
So you can see, setting the correct level damping is very important for car control when it comes to suspension tuning.
We'll look into overdamped and critically damped cases in the next post.
Credits:
All equations in this tutorial were created with QuickLatex
All graphs were created with www.desmos.com/calculator
First Order Differential Equations
- Introduction to Differential Equations - Part 1
- Differential Equations: Order and Linearity
- First-Order Differential Equations with Separable Variables - Example 1
- Separable Differential Equations - Example 2
- Modelling Exponential Growth of Bacteria with dy/dx = ky
- Modelling the Decay of Nuclear Medicine with dy/dx = -ky
- Exponential Decay: The mathematics behind your Camping Torch with dy/dx = -ky
- Mixing Salt & Water with Separable Differential Equations
- How Newton's Law of Cooling cools your Champagne
- The Logistic Model for Population Growth
- Predicting World Population Growth with the Logistic Model - Part 1
- Predicting World Population Growth with the Logistic Model - Part 2
- What's faster? Going up or Coming Down?
First order Non-linear Differential Equations
- There's a hole in my bucket! Let's turn it into a cool Math problem!
- The Calculus of Hot Chocolate Pouring!
- Foxes hunting Bunnies: Population Modelling with the Predator-Prey Equations
Second Order Differential Equations
- Introduction to Second Order Differential Equations
- Finding a Basis for solutions of Second Order ODE's
- Roots of Homogeneous Second Order ODE's and the Nature their Solutions
- Modelling with Second Order ODE's: Undamped Free Oscillations
- Modelling car suspension with Second Order ODE's: Damped Free Oscillations
Please give me an Upvote and Resteem if you have found this tutorial helpful.
Please ask me a maths question by commenting below and I will try to help you in future posts.
Tip me some DogeCoin: A4f3URZSWDoJCkWhVttbR3RjGHRSuLpaP3
Tip me at PayPal: https://paypal.me/MasterWu
Content type: long, technical
![clar.png](https://steemitimages.com/DQmfH2dRuSi1jUNra9ebptNK5xxkQsHpM2boFFZT7BbFhX6/clar.png)
![spel.png](https://steemitimages.com/DQmSdZZAMgcwtaN86uYm2Q39ALBT4AcJbw3UQjmJW7ZgvrD/spel.png)
Awarded 5 out of 6 owls:
Details: The originality owl was not awarded since it requires that the math is explained in creative/novel way.