Modelling Car Suspension with ODE's: Damped Free Oscillations Part 2

in #steemstem6 years ago (edited)

In the mathematical modelling of the mass-damper-spring system (Figure 1), we've explored the underdamped case in the previous post.

s26.png
Figure 1. Schematic of the mass-damper-spring system (adapted from Wikimedia Commons)

We seen that with an under damped suspension system (case I), a car can wallow up and down whenever it goes over a speed bump, making it very difficult for the driver to control the vehicle.

Before we explore the case II, critical damping, let's recap the equation of motion for the mass-damper-spring system...

s1.png

where:

  • s23.png = sprung mass [kg]
  • s24.png = damping coefficient [Ns/m]
  • s25.png = spring constant [N/m]

...we can have 3 different cases...

Case I < 0 Under damping
Case II = 0 Critical damping
Case III > 0 Over damping

Case II: Critically Damped

Critical damping is the border case between oscillatory (under damped) and non-oscillatory (over damped) motions.

Let's take explore our car suspension model with the damper tuned to be critically damped. The parameters from the previous post...

  • sprung mass m = 300kg
  • spring constant k = 1200N/m
  • initial displacement q13.png = 0.15m
  • initial velocity q27.png = 0m/s

The system is critically damped when...

s10.png

Substituting above parameters, we need to tune the damper to...

s11.png

A diagram of the system is shown in Figure 2 below...

s12.png
Figure 2. Suspension schematic with parameters required for critical damping

Now that we have the damping coefficient for the critically damped case, substituting these into equation (1), the equation of motion becomes...

s2.png

The characteristic equation of (2)...

s3.png

The only root of equation (3) is...

s4.png

Now according to equation (8) of Post #3, the general solution is...

s5.png

Applying the first initial condition (initial displacement)...

s6.png

The first derivative of the general solution, equation (4) is...

s7.png

Now we can apply the second initial condition (initial velocity)...

s8.png

Therefore, the particular solution is...

s9.png

The behaviour of the particular solution is shown in Figure 3...

s13.png
Figure 3. Behaviour of the suspension system with s28.png and s27.png

As you can see, the car pretty much returns to its equilibrium position in less than 4 seconds, which is quite a slow response for most real-world applications.

Exploring the initial value problem

Let's now explore the solution with a few different permutations for the initial conditions...

  • initial displacement q13.png = 0.15m, therefore A = 0.15
  • initial velocity q27.png = -0.6m/s

Substituting this values into equation (5)...

s14.png

This gives rise to the particular solution is...

s15.png


What about the following?...

  • initial displacement q13.png = 0.15m, therefore A = 0.15
  • initial velocity q27.png = -0.3m/s

This gives rise to the particular solution...

s16.png


And the following permutation...

  • initial displacement q13.png = 0.15m, therefore A = 0.15
  • initial velocity q27.png = 0.3m/s

...gives rise to the particular solution...

s17.png


All of the above solutions are superimposed in Figure 4...

s22.png
Figure 4. Various particular solutions to the critically damped car suspension problem.


Finally, let's see how the suspension behaves with zero initial displacement, but a non-zero initial velocity. For instance...

  • initial displacement q13.png = 0.0m, therefore A = 0.0
  • initial velocity q27.png = 0.3m/s

Again, from equation (5), we have...

s30.png

...giving rise to the particular solution...

s31.png

A few solutions showing the behaviour of the mechanism with varying values of s32.png...

s29.png
Figure 5. Various particular solutions to the critically damped car suspension problem.

Discussion about the above case...

Critical damping is what car manufacturers aim for in most cases. As we have seen with under damping, the residual vibrations that occur after hitting a bump can be very dangerous in terms of steering and car control. It can also cause passengers who are susceptible to motion sickness to become ill.

With over damping, which we will explore in the next post, the shock absorber is stiff and sluggish, which means it is slow to respond to bumps. This means that bumps are not absorbed as much and car returns to its equilibrium position slowly, making the ride quite uncomfortable. It feels like a car with no suspension, and may compromise the car's grip on rougher roads, because the wheels are in less contact with the ground at certain times (the car might do a "jump" when it hits a bump).

As I have mentioned previously, the parameters used for the this particular suspension set up is quite "soft" and "squishy" for most real-world applications. We have a soft spring at 1200N/m, and the critical damping coefficient for this softly sprung car will dissipate a bump of 15cm in about 4 seconds.


Credits:

All equations in this tutorial were created with QuickLatex

All graphs were created with www.desmos.com/calculator


First Order Differential Equations

  1. Introduction to Differential Equations - Part 1
  2. Differential Equations: Order and Linearity
  3. First-Order Differential Equations with Separable Variables - Example 1
  4. Separable Differential Equations - Example 2
  5. Modelling Exponential Growth of Bacteria with dy/dx = ky
  6. Modelling the Decay of Nuclear Medicine with dy/dx = -ky
  7. Exponential Decay: The mathematics behind your Camping Torch with dy/dx = -ky
  8. Mixing Salt & Water with Separable Differential Equations
  9. How Newton's Law of Cooling cools your Champagne
  10. The Logistic Model for Population Growth
  11. Predicting World Population Growth with the Logistic Model - Part 1
  12. Predicting World Population Growth with the Logistic Model - Part 2
  13. What's faster? Going up or Coming Down?

First order Non-linear Differential Equations

  1. There's a hole in my bucket! Let's turn it into a cool Math problem!
  2. The Calculus of Hot Chocolate Pouring!
  3. Foxes hunting Bunnies: Population Modelling with the Predator-Prey Equations

Second Order Differential Equations

  1. Introduction to Second Order Differential Equations
  2. Finding a Basis for solutions of Second Order ODE's
  3. Roots of Homogeneous Second Order ODE's and the Nature their Solutions
  4. Modelling with Second Order ODE's: Undamped Free Oscillations
  5. Modelling Car Suspension with ODE's: Damped Free Oscillations Part 1
  6. Modelling Car Suspension with ODE's: Damped Free Oscillations Part 2

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

Sort:  

Excellent Mathematical description. Often hard to explain, but you did so with ease. Keep it up :)

Coin Marketplace

STEEM 0.29
TRX 0.12
JST 0.034
BTC 63755.07
ETH 3328.24
USDT 1.00
SBD 3.91