Use the model from Section 4.4.2 in the book for a 1/4-car model. I looked up the values for a Mazda Miata:
- k=6000 kg/m (note the weird units, k = 6000*9.81 N/m)
- m = 250 kg
- b = 5370 N s/m.
Write a computer program (Matlab is fine as long as you trust the answer) to compute how the car moves up and down for different driving velocities and different road types. Once you have that, submit plots for
- Ride height vs x
- Suspension compression vs x (the suspension will "bottom out" at 4cm deflection)
- The net acceleration vs x (if the acceleration is more negative than -9.81 the wheel will leave the road).
Added October 24, 8:00pm If you use ode45() in Matlab, it needs a function that gives the right hand side of the differential equation. For some reason, it really, really wants that function to be a function of x and t only. But in order for us to compute y and dydt we need to send it other variables too, like wavenos, phi and sk0. The way to deal with this is with an anonymous function. If my right hand side function is called carrhs.m then use the following syntax in ode45():
Code: Select all
[t z] = ode45(@(t,x) carrhs(t,x,sk0,wavenos,phi,v,m,b,k),[0 timetodrive1km], [icz iczdot])
Code: Select all
@(t,x) carrhs(t,x,sk0,wavenos,phi,v,m,b,k)