A Single-Degree-of-Freedom Vehicle Model The mass–spring–damper system shown in Figure 5.101 represents a vehicle traveling on a rough road. Assume that the surface of the road can be approximated as a sine wave z = Z0sin(ωt), where Z0 = 0.01 m and ω = 3.5 rad/s. The mathematical model of the

Question

A Single-Degree-of-Freedom Vehicle Model

The mass-spring-damper system shown in Figure 5.101 represents a vehicle traveling on a rough road. Assume that the surface of the road can be approximated as a sine wave [latex]z=Z_{0} \sin (\omega t)[/latex], where [latex]Z_{0}=0.01 \mathrm{~m}[/latex] and [latex]\omega=3.5 \mathrm{rad} / \mathrm{s}[/latex]. The mathematical model of the system is given by an ordinary differential equation

[latex]m \ddot{x}+b \dot{x}+k x=b \dot{z}+k z[/latex]

where [latex]m=3000 \mathrm{~kg}, b=2000 \mathrm{~N} \cdot \mathrm{s} / \mathrm{m}[/latex], and [latex]k=50 \mathrm{kN} / \mathrm{m}[/latex].

a. Build a Simulink model of the system based on the mathematical representation and find the displacement output [latex]x(t)[/latex].

b. Convert the ordinary differential equation to a transfer function and repeat Part (a). Assume zero initial conditions.

c. Build a Simscape model of the physical system and find the displacement output [latex]x(t)[/latex].

Accepted Answer

a. Solving for the highest derivative of the output [latex]x[/latex] gives

[latex]\ddot{x}=\frac{1}{m}(k z+b \dot{z}-k x-b \dot{x}) .[/latex]

The corresponding Simulink block diagram is shown in Figure 5.102. Note that the displacement input [latex]z(t)[/latex] is a sine function, which can be defined using a Sine Wave block available in the library of Sources. Double-click on the block and type 0.01 for the Amplitude and 3.5 for the Frequency to define the input [latex]z(t)=0.01 \sin (3.5 t)[/latex].

b. The transfer function relating the input [latex]z(t)[/latex] to the output [latex]x(t)[/latex] is

[latex]\frac{X(s)}{Z(s)}=\frac{b s+k}{m s^{2}+b s+k} ext {. }[/latex]

The Simulink block diagram built based on the transfer function is shown in Figure 5.103, where a Transfer Fcn block is used to represent the vehicle system. Double-click on the block and type [ [latex]\left[\begin{array}{ll}\mathrm{b} & \mathrm{k}\end{array} ight][/latex] for the Numerator coefficient and [latex]\left[\begin{array}{lll}m & \mathrm{~b} & \mathrm{k}\end{array} ight][/latex] for the Denominator coefficient to define the transfer function [latex]X(s) / Z(s)[/latex].

c. The Simscape block diagram corresponding to the physical system is shown in Figure 5.104, which can be created by following these steps:

1. Type ssc_new at the MATLAB Command window to open the main Simscape library and create a new model.

2. Open the library of Simscape/Foundation Library/Mechanical/Translational Elements and drag the Mass, Translational Damper, and Translational spring into the model window. Double-click on these blocks to define the parameters Mass, Damping coefficient, and spring rate as [latex]\mathrm{m}, \mathrm{b}[/latex], and [latex]\mathrm{k}[/latex].

3. To add the representation of the displacement input, open the library of simscape/ Foundation Library/Mechanical/Mechanical Sources and drag the Ideal Translational Velocity Source into the model window. Note that two types of inputs are available in Simscape for translational mechanical systems, and they are used to define either a force or a velocity input. Therefore, the displacement input function in this example must be converted to the velocity by taking the time derivative.

4. To add the sensor to measure the displacement of the mass, open the library of Simscape/Foundation Library/Mechanical/Mechanical Sensors and drag the Ideal Translational Motion Sensor into the model window.

5. Now open Simulink libraries to add the source and the scope. As mentioned in Step 3, taking the time derivative of [latex]z(t)[/latex] gives the velocity input function [latex]\dot{z}(t)=0.035 \cos (3.5 t)[/latex]. Open the library of Simulink/Sources to drag the Clock block and open the library of Simulink/User-Defined Functions to drag the Fun block. Double-click on the Fun block and type [latex]0.035^{*} \cos (3.5 * u)[/latex] for the Expression, where [latex]u[/latex] is the default name of the input to the Fun block, and here it represents the time [latex]t[/latex]. Note that the Simulink-PS Converter and PS-Simulink Converter blocks are used to convert Simulink signals into physical signals or vice versa.

6. Orient the blocks and connect them as shown in Figure 5.104.

Define the values of the parameters [latex]m, b[/latex], and [latex]k[/latex] in the MATLAB Command window. Run all simulations and the same curve as shown in Figure 5.105 can be obtained, which is the resulting displacement output [latex]x(t)[/latex] of the vehicle due to the roughness of the road.

Question 5.18: A Single-Degree-of-Freedom Vehicle Model The mass–spring–dam...

Related Answered Questions

Verified Answer:

Verified Answer:

A Two-Degree-of-Freedom Mass–Spring System Consider the two-degree-of-freedom mass–spring system shown in Figure 5.106. The mathematical model of the system is given by a set of ordinary differential equations m1x1 + (k1 + k2)x1 – k2x2 = 0, m2x2 – k2x1 + k2x2 = 0, where m1 = m2 = 5 kg, ...

Verified Answer:

Verified Answer:

A Single-Degree-of-Freedom Gear–Train System For the gear–train system shown in Figure 5.95a, derive the differential equation of motion. The mass moments of inertia of the two gears about their respective fixed centers are IC1 and IC2. ...

Verified Answer:

Verified Answer:

Verified Answer:

A Cart–Inverted-Pendulum System Consider the mechanical system shown in Figure 5.79, in which a uniform rod of mass M and length L is pivoted on a cart of mass m. An external force f is applied to the cart. Assume that the pendulum is constrained to move in a vertical plane, and the cart moves ...

Verified Answer:

Verified Answer:

A Pure Rolling Disk Consider the system shown in Figure 5.60a, in which a uniform disk of mass m and radius r rolls on a horizontal surface. A translational spring of stiffness k is attached to the disk. Assuming that there is no slipping between the disk and the surface, derive the differential ...

Verified Answer: