CH 07
CH 07
CH 07
7.1 Re-write the first-order ODE in standard form (divide all terms by 10):
0.2 y y 0.3u
The single root is r = –1/0.2 = –5 and therefore the response is a stable exponential function. The
time constant is = 0.2 s and hence the system reaches steady state in tS = 4 = 0.8 s. The steady-
state value for a constant input u = 6 is ySS = (0.3)(6) = 1.8. A hand-drawn sketch of the response
would match the response plot (below) and include labels for the initial condition y(0) = –2,
settling time tS = 0.8 s, and steady-state response ySS = 1.8.
1.5
0.5
Output y
-0.5
-1
-1.5
-2
-2.5
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Time, s
2r 2 12r 68 0
The two roots are complex: r1, 2 3 j 5 . Hence, the homogeneous response does exhibit
oscillations. Another way to show this is to re-write the I/O equation in the standard form for a
second-order system:
4 4
tS 1.3333 s
n 0.514 34
2.5
2
Output y
1.5
0.5
-0.5
0 0.5 1 1.5 2
Time, s
2
Copyright © 2016 John Wiley & Sons
Chapter 7
The two roots are real and negative: r1 = –1 and r2 = –4.5 and hence the homogeneous response is
comprised of two decaying exponential functions (no oscillations).
b) The homogeneous response has the form yH (t ) c1et c2e4.5t and the slowest
exponential mode “dies out” at time tS = 4 s. Hence the settling time for the system is 4 s.
c) The homogeneous response will begin at the initial condition y(0) = 3 with “zero slope”
and then exponentially decay to zero in about 4 s. A sketch would match the plot below.
2.5
2
Output y
1.5
0.5
0
0 1 2 3 4 5 6
Time, s
3
Copyright © 2016 John Wiley & Sons
Chapter 7
7.4 In all cases the denominator polynomial (characteristic equation) will have the form
(s r1 )( s r2 ) 0
c
Hence the transfer function has the form G ( s )
s 2 .7 s 0 .5
2
Since the DC gain of this transfer function is G ( s 0) c / 0.5 0.5 we obtain c = 0.25 and the
transfer function is
0.25
G (s)
s 2 .7 s 0 .5
2
180
G (s) Note: G(s = 0) = 180/30 = 6
s 13s 30
2
2500
G ( s) Note: G(s = 0) = 2500/20 = 125
s 4 s 20
2
d) Roots r1, 2 0.4 j1.6 : ( s 0.4 j1.6)( s 0.4 j1.6) s 2 0.8s 2.72
0.0544 0.0544
G (s) Note: G(s=0) = 0.02
s 0.8s 2.72
2
2.72
4
Copyright © 2016 John Wiley & Sons
Chapter 7
7.5 a) We can obtain the steady-state response to a constant input by using the DC gain:
84 84
G( s) G ( s 0) 0.2593
(3s 21s 36)( s 2s 9)
2 2
(36)(9)
b) The four characteristic roots are determined by setting the denominator to zero
The three exponential decay modes for the transient response are e 3t , e 4 t (for the two real
roots) and e t (for the two complex roots). Since the slowest exponential mode is e t , the
overall settling time is tS = 4 s since e 4 is “small.”
5
Copyright © 2016 John Wiley & Sons
Chapter 7
EO ( s ) 1 1
G (s) (using RC = 0.01 s)
Ein ( s ) RCs 1 0.01s 1
Clearly the DC gain is 1 and the time constant is = RC = 0.01 s. Since the input is a constant
0.4 V, the steady-state output voltage is also 0.4 V. The output voltage will show an exponential
rise from its initial value (zero) to 0.4 V and reach steady state in settling time tS = 4 = 0.04 s. A
sketch should show eO(0) = 0, eO () 0.4 V, and settling time (i.e., match the plot below).
0.5
0.4
Output voltage, V
0.3
0.2
0.1
0
0 0.01 0.02 0.03 0.04 0.05
Time, s
6
Copyright © 2016 John Wiley & Sons
Chapter 7
First, consider the response for time 0 t 1 (Switch 1) when ein(t) = 4 V. The time constant of
the “Switch 1” model is 1 L / R1 = 0.2/1.6 = 0.125 s. Hence the settling time is 41 = 0.5 s
which is less than the switching time (1 s), so the current reaches a steady-state value. The
steady-state current for the “Switch 1” model is I LS S ein / R1 4/1.6 = 2.5 A.
When time t > 1 s, we have the “Switch 2” model with time constant 2 L /( R1 R2 ) = 0.0357
s. Hence the settling time is 42 = 0.1429 s (relative to the switch time t = 1), and the current IL
decays to zero and reaches zero at approximately time t = 1.1429 s. A sketch will match the plot
below.
2.5
2
Current, A
1.5
0.5
0
0 0.5 1 1.5
Time, s
7
Copyright © 2016 John Wiley & Sons
Chapter 7
7.8 The characteristic equation is the denominator polynomial of G(s) set to zero:
2s 2 3s 24 0 or s 2 1.5s 12 0
The roots (poles) are s1, 2 0.75 j 3.3819 . Because the roots (or poles) are complex, the
transient response does exhibit oscillations.
8
Copyright © 2016 John Wiley & Sons
Chapter 7
7.9 a) The characteristic equation is the denominator polynomial of G(s) set to zero
The roots (poles) are s1, 2 4.6667 j 28.8367 . Because the roots (or poles) are complex, the
transient response does exhibit oscillations.
Note that the standard second-order form is s 2 2 n s n2 0 . Hence the undamped natural
frequency is n 853.33 29.2119 rad/s and = 0.1598 < 1 (underdamped).
c) The period of oscillation is 2/d = 0.2178 s (= 1/4.59 Hz). The settling time for the
underdamped second-order system is tS = 4/(n) = 0.8571 s. Hence the transient response
has approximately 0.8571/0.2178 = 3.9 (~ 4) cycles before reaching steady state.
9
Copyright © 2016 John Wiley & Sons
Chapter 7
7.10 The damping ratio is the magnitude of the real part of the root divided by the distance
from the root to the origin (see Fig. 7.17). If the complex root is s a jb then we have
a
a 2 b2
10
Copyright © 2016 John Wiley & Sons
Chapter 7
7.11 We will use the log-decrement method. The peak values and peak times of the impulse
response are estimated from the impulse-response plot:
Hence the period is approximately Tperiod = 0.5 – 0.1 = 0.4 s. The log decrement is
x1 0.026
ln ln 0.7732
x2 0.012
The damping ratio is approximately
= 0.1221
4 2 2
n d / 1 2 = 15.8264 rad/s
11
Copyright © 2016 John Wiley & Sons
Chapter 7
7.12 a) The general solution form for the homogeneous response of an underdamped system
with two complex roots at r1, 2 j is
We see that the real part of the root will be the exponent of the exponential function while the
imaginary part will be the frequency of the sinusoidal function. Because we have two pairs of
complex roots we will have two damped sinusoidal functions in the form of Eq. (A). Using the
two root pairs given in the problem, the homogeneous response has the form
The unknown constants are c1, c2, 1, and 2. We could have used sine instead of cosine.
b) The two exponential functions e 1.2t and e 0.4t determine the decay rate for the two damped
sinusoidal functions. The “slowest” decaying exponential function is e 0.4t which requires about
10 s to approximately reach zero since e 4 0.0183 is “small.” Hence the settling time is 10 s.
12
Copyright © 2016 John Wiley & Sons
Chapter 7
7.13 Because the free response is x(t) = 6cos50t , we note that the undamped natural frequency
is n = 50 rad/s.
The mechanical model of an unforced mass-spring system is mx kx 0 . The standard I/O
equation form of a second-order system with no damping ( = 0) is x n2 x 0 . Hence the
undamped natural frequency of the mass-spring system is
k
n = 50 rad/s
m
13
Copyright © 2016 John Wiley & Sons
Chapter 7
2r 2 8r 6 0 or 2(r 1)( r 3) 0
b) The following transfer function can be derived from the given I/O equation:
Y (s) 3
G ( s) 2
U ( s ) 2 s 8s 6
The poles of G(s) are determined by setting the denominator polynomial to zero:
2 s 2 8s 6 0
Therefore the two poles of G(s) are s1 1 and s2 3 , which are identical to the roots in (a).
c) We can obtain a SSR for the states x1 = y and x2 y . The resulting state equation is
0 1 0
x x u
3 4 1.5
1
I A 2 4 3 0
3 4
The eigenvalues are 1 1 , 2 3 which are identical to the roots and poles in (a) and (b).
d) Because the two roots are real and negative the zero-input response consists of two decaying
exponential functions
y(t ) c1et c2e 3t
The response starts at y(0) = 5 with “zero slope” due to y (0) 0 and eventually decays to zero.
t
The slowest exponential decay mode is e which takes about 4 s to reach zero.
14
Copyright © 2016 John Wiley & Sons
Chapter 7
7.15 First, we determine the poles of the third-order transfer function G1(s)
s 3 20.2s 2 124s 24 0
The three poles are s1 0.2 and s2 10 j 4.4721 . The two underdamped poles show very
high damping ( = 0.913) and a relatively short settling time tS = 0.4 s. The single real pole
( s1 0.2 ) exhibits an exponential decay e 0.2t which takes 20 s to die out to zero. Hence the
damped sinusoidal (second-order) component of the transient response of G1(s) is barely
noticeable and the first-order component e 0.2t dominates the transient response. It is easy to see
that the single pole of the first-order transfer function G2(s) is also s 0.2 and hence its
transient response is also e 0.2t , which matches the dominant pole of G1(s).
Finally, note that the DC gain of the third-order system is G1(s = 0) = 72/24 = 3 which is also the
DC gain of the first-order system, i.e., G2(s = 0) = 0.6/0.2 = 3. So both transfer functions show
the same steady-state response to a constant input. For these reasons the first-order system G2(s)
is an excellent low-order approximation of the third-order system G1(s).
We can obtain the third-order and first-order step responses using the following MATLAB
commands:
>> sysG1 = tf(72,[1 20.2 124 24]); % define transfer function G1(s)
>> sysG2 = tf(0.6,[1 0.2]); % define transfer function G2(s)
>> t = 0:0.01:25; % define time vector
>> [y1,t] = step(sysG1,t) ; % unit-step response of G1(s)
>> [y2,t] = step(sysG2,t) ; % unit-step response of G2(s)
>> plot(t,y1,t,y2) % plot both unit-step responses
The plot is below. The two unit-step responses are almost identical, and hence G2(s) is an
excellent low-order approximate model for system G1(s).
3.5
2.5
2
Output y
0.5
0
0 5 10 15 20 25
Time, s
15
Copyright © 2016 John Wiley & Sons
Chapter 7
0.2 0.6
I A 2 4.2 2 0
2 4
b) The eigenvalues can be computed using the simple MATLAB commands; matches part (a)
c) Because the two roots are real and negative, the free (zero-input) response will consist of two
decaying exponential functions, or y (t ) c1e0.5476t c2e3.6524t . The “slow” exponential
mode is e 0.5476t which takes about 7.3 s to decay to zero.
d) MATLAB commands are probably the easiest method: use the command initial
The plot (below) verifies the free response described in part (c).
0.5
-0.5
Output y
-1
-1.5
-2
-2.5
0 2 4 6 8 10 12
Time, s
16
Copyright © 2016 John Wiley & Sons
Chapter 7
7.17 a) The MATLAB commands define state matrix A and determine its eigenvalues:
b) Because the three eigenvalues (roots) are real and negative, the free (zero-input) response will
consist of three decaying exponential functions, or y (t ) c1e t c2e 2t c3e 6t . The “slow”
exponential mode is e t which takes about 4 s to decay to zero.
c) MATLAB commands are probably the easiest method: use the command initial
The plot below verifies the free response described in part (b).
1.8
1.6
1.4
1.2
Output, y
0.8
0.6
0.4
0.2
0
0 1 2 3 4 5
Time, s
17
Copyright © 2016 John Wiley & Sons
Chapter 7
7.18 a) The MATLAB commands define state matrix A and determine its eigenvalues:
b) Because the three eigenvalues (roots) have negative real parts the free response will eventually
decay to zero. The free response will be the sum of one decaying exponential function and a
damped sinusoidal function, or y(t ) c1e 5t c2e 2t sin 8t . The “slow” exponential mode is
e 2 t which takes about 2 s to decay to zero. The free response will oscillate at frequency 8 rad/s.
c) MATLAB commands are probably the easiest method. The plot verifies part (b).
0.5
-0.5
-1
Output, y
-1.5
-2
-2.5
-3
-3.5
-4
0 1 2 3 4 5
Time, s
18
Copyright © 2016 John Wiley & Sons
Chapter 7
a) Using velocity v(t ) x (t ) as the dynamic variable the model becomes mv bv f a (t ) .
Because the force magnitude is 5 N and the very short pulse duration is 0.02 s, the input can be
represented by an ideal impulse 0.1 (t ) where (5 N)(0.02s) = 0.1 N-s is the weight or strength of
the impulse. Substituting the numerical values for m and b the system model becomes
Hence, the time constant is = 0.1667 s and the initial magnitude of the impulse response is
(0.1)(0.3333)/0.1667 = 0.2 m/s = v(0+). The impulse response is an exponential
t /
function e e . The impulse response for velocity is v(t ) 0.2e6t m/s.
6 t
b) We can integrate the velocity response in part (a) to obtain the position response:
t t
x(t ) x(0) v( )d 0.2e 6 d
0 0
Note x(0) = 0 and is the dummy integration variable. The solution is x(t ) 0.0333(1 e6t ) m
c) The Simulink model (below) will solve the mathematical model mx bx f a (t ) by
integrating acceleration twice. Both integrators have zero initial conditions. Input force is a
very short-duration pulse with a magnitude of 5 N and is created by adding two step
functions. The fixed simulation step size is 10-4 s.
Plots of velocity v(t) and position x(t) follow (next page). Note that velocity starts at rest and
exhibits a very sharp jump to nearly 0.2 m/s due to the short-duration pulse input. The simulation
nearly matches the analytic solution v(t ) 0.2e6t m/s in part (a) which was computed using an
ideal impulse. Position x(t) shows an exponential rise (from zero) to the computed steady-state of
0.0333 m.
19
Copyright © 2016 John Wiley & Sons
Chapter 7
0.2
0.18
0.16
0.14
Velocity, m/s
0.12
0.1
0.08
0.06
0.04
0.02
0
0 0.2 0.4 0.6 0.8 1
Time, s
0.035
0.03
0.025
Position, m
0.02
0.015
0.01
0.005
0
0 0.2 0.4 0.6 0.8 1
Time, s
20
Copyright © 2016 John Wiley & Sons
Chapter 7
800
G(s)
2s 24s 144s 640
3 2
The characteristic roots are the poles of the transfer function, i.e., 2s3 + 24s2 + 144s + 640 = 0.
Using MATLAB, the roots are
The three roots are –8, –2 + j6, and –2 – j6. Because the two of the roots are complex the
transient response will exhibit oscillations with frequency 6 rad/s (period = 2/6 = 1.047 s). The
real root will contribute an exponential mode e 8t which dies out in 0.5 s. The complex roots
will contribute a damped sinusoidal mode with an exponential envelop of e 2 t (dies out in 2 s).
Because the DC gain is G(s = 0) = 800/640 = 1.25 and the input is a unit step, the steady-state
response is ySS = 1.25. MATLAB’s step command can be used to obtain a numerical solution:
1.6
1.4
1.2
1
Output, y
0.8
0.6
0.4
0.2
0
0 0.5 1 1.5 2 2.5 3
Time, s
21
Copyright © 2016 John Wiley & Sons
Chapter 7
EO ( s ) 1 1
G (s)
Ein ( s ) RCs 1 0.01s 1
The Simulink model simulates the step response to a 0.4-V input. The output voltage response
eO(t) matches the sketch in Problem 7.6.
0.5
0.4
Output voltage, V
0.3
0.2
0.1
0
0 0.01 0.02 0.03 0.04 0.05
Time, s
22
Copyright © 2016 John Wiley & Sons
Chapter 7
Y (s) 0.25
G(s) 2
U ( s) s 2s 9
a) Because the system is underdamped (note the poles are 1 j 2.8284 ) we can compute the
step-response performance metrics using Table 7-4 (summarized below)
/ 1 2
Maximum overshoot M OS e = 0.3293 (or 32.9%)
peak value = (1.3293)(0.1111 m) = 0.1477 m
4
Settling time t S =4s
n
2
Period of oscillation Tperiod = 2.2214 s
n 1 2
2 1 2
Number of cycles to S-S N cycles = 1.8 cycles
An accurate sketch would show these response values - see simulation plot in part (b).
b) The below Simulink model will produce the step response - see plot of y(t) on the next page.
23
Copyright © 2016 John Wiley & Sons
Chapter 7
0.16
0.14
0.12
0.1
Position, y, m
0.08
0.06
0.04
0.02
0
0 1 2 3 4 5 6
Time, s
24
Copyright © 2016 John Wiley & Sons
Chapter 7
7.23 The I/O equation of the 1-DOF mechanical system can be derived from the transfer
function G(s) in Problem 7.22:
y 2 y 9 y 0.25u
a) The system is underdamped (note the roots are 1 j 2.8284 ) and we can use the second-
order underdamped metrics from Table 7-4 (see Problem 7.22)
Steady-state response is zero since input u = 0 and the damped sinusoid decays to zero.
4
Settling time t S =4s
n
2
Period of oscillation Tperiod = 2.2214 s
n 1 2
2 1 2
Number of cycles to S-S N cycles = 1.8 cycles
An accurate sketch would show these response values - see simulation plot in part (b) on
the next page.
b) We cannot use transfer functions since the system has non-zero initial conditions. The
below Simulink model uses integrator blocks and will produce the free response. Note that the
input is zero and the integrators must have initial conditions y (0) 0.04 m and y (0) 0 .
25
Copyright © 2016 John Wiley & Sons
Chapter 7
0.04
0.03
0.02
Position, y, m
0.01
-0.01
-0.02
0 1 2 3 4 5 6
Time, s
26
Copyright © 2016 John Wiley & Sons
Chapter 7
J k Tin (t )
( s ) 1
G (s)
Tin ( s ) 0.2 s 100
2
Hence the undamped natural frequency is n 100 / 0.2 = 22.3607 rad/s. The damping ratio
is zero since there is no friction. Hence, the homogeneous (or natural) response will be an
undamped harmonic (sinusoidal) function with the form
The forced response or particular solution (due to the sinusoidal input torque) will be a sinusoidal
function with frequency 3 rad/s (same as the input) with unknown coefficients:
Hence the total solution is the sum of sinusoidal function at two frequencies, n and 3 rad/s.
b) The Simulink model (below) will produce the desired system response.
The complete response (t) is shown on the next page. Note that the response is a summation of
two harmonic (sinusoidal) functions: the first harmonic function has a frequency of 3 rad/s
(period = 2.1 s) due to the torque input and the second harmonic function has frequency 22.36
rad/s (period = 0.28 s) due to k/J. This dual-frequency harmonic function matches the result in
part (a).
27
Copyright © 2016 John Wiley & Sons
Chapter 7
-3
x 10
8
-2
-4
-6
-8
0 2 4 6 8 10
Time, s
28
Copyright © 2016 John Wiley & Sons
Chapter 7
f a (t )
x 2 n x n2 x
m
The numerical values determined in Problem 7.11 are = 0.1221 (damping ratio) and n =
15.8264 rad/s (undamped natural frequency), and mass m = 0.2 kg. Hence the transfer function is
X (s) 1/ m 5
G (s) 2 2
Fa ( s ) s 2 n s n s 3.865s 250.47
2
The Simulink model below shows that two step functions are added to create a pulse input. One
step input is delayed by 0.01 s.
The magnitude of the pulse input was varied until the peak response of x(t) nearly matched 0.026
m (the peak response from the plot in Problem 7.11). A 10-N pulse input with pulse width of
0.01 s produced the following response (below) which nearly matches the plot in Problem 7.11.
Because the pulse magnitude is 10 N and pulse width is 0.01 s, the strength of the impulse input
is (10 N)(0.01 s) = 0.1 N-s.
0.03
0.025
0.02
0.015
Displacement, x, m
0.01
0.005
-0.005
-0.01
-0.015
-0.02
0 0.5 1 1.5 2 2.5 3
Time, s
29
Copyright © 2016 John Wiley & Sons
Chapter 7
7.26 The mathematical model of the thermal system is (see Example 4.7)
REQCT T REQqBH Ta
6
1 1 1 1 1
The equivalent thermal resistance of all surfaces is
REQ R1 R2 R6 i 1 Ri
Or, REQ = 0.018014 deg C-s/J. Hence the time constant of the first-order system is = REQC =
8,700.8 s. The settling time is tS = 4= 34,803 s = 9.67 hrs. At steady-state we have T 0 so
the first-order model becomes
T () REQ qBH Ta = (0.018014 deg C-s/J)(1000 J/s) + 10 deg C = 28.01 deg C.
The settling time and steady-state temperature match the numerical results from Problem 6.25. A
sketch of T(t) would be an exponential rise from 10 deg C to 28 deg C with settling time of 9.7
hrs (see solution of Problem 6.25).
30
Copyright © 2016 John Wiley & Sons
Chapter 7
7.27 The mathematical model of the LC tuner circuit (see Problem 3.25) is
Therefore the second-order system has zero damping ( = 0) and undamped natural frequency
n 1/ LC 4,082.5 rad/s. The voltage response eC(t) is an undamped sinusoidal function
with initial voltage eC(0) = 2.5 V and frequency 4082.5 rad/s (hence, use the cosine function):
Sketches would show that the capacitor voltage oscillates between 2.5 V at frequency 4082.5
rad/s (period = 0.0015 s) and the current oscillates between 0.2041 A at the same frequency.
Plots of capacitor voltage eC(t) and current I(t) from the numerical simulation are below. These
plots verify the sketches of the cosine and sine functions for eC(t) and I(t)
3 0.25
2.5
0.2
2
0.15
1.5
Capacitor voltage, e , V
0.1
C
1
0.5 0.05
Current, A
0 0
-0.5 -0.05
-1
-0.1
-1.5
-0.15
-2
-2.5 -0.2
-3 -0.25
0 1 2 3 4 5 0 1 2 3 4 5
Time, s -3 Time, s -3
x 10 x 10
31
Copyright © 2016 John Wiley & Sons
Chapter 7
7.28 The mathematical model of the LC tuner circuit (see Problem 3.25) is LCeC eC 0
1 1
The first state equation is x1 eC I x2
C C
1 1
The second state equation is x2 I CeC eC x1
L L
0 1/ C
x x State equation
1 / L 0
1/ C 1
det( I A) det 2 0
1 / L LC
Note that in Problem 7.27 the capacitor voltage and current of the LC tuner circuit are
eC (t ) 2.5 cos 4082.5t V and I (t ) 0.2041sin 4082.5t A. Both are undamped harmonic
functions with a frequency of 4082.5 rad/s. Because the two eigenvalues are imaginary numbers
the homogeneous response is an undamped sinusoidal function with a frequency equal to the
imaginary number ( = 4082.5 rad/s in this case).
32
Copyright © 2016 John Wiley & Sons
Chapter 7
7.29 The mathematical model of the RC circuit (see Problem 3.26) is RC eC eC 0
Using the values R = 10 and C = 185(10-6) F, the time constant is = RC = 0.00185 s. Hence
the settling time is tS = 4 = 0.0074 s. Clearly at steady state eC 0 and therefore eC () 0 .
A sketch of the capacitor voltage eC(t) would be an exponential decay from eC(0) = 200 V to zero
at approximately 0.0074 s. The mathematical solution is eC (t ) 200e t / 200e 540.54t V.
The current is I (t ) CeC which is the product of the slope (time derivative) of eC(t) and
capacitance C. The derivative is eC (t ) 1.0811(105 )e 540.54t so I (t ) CeC = 20e 540.54t A.
A sketch of the current I(t) would be an exponential rise from –10 A to zero at about tS = 0.0074
s. Sketches would match the plots below.
220
200
180
160
Capacitor voltage, e , V
C
140
120
100
80
60
40
20
0
0 0.002 0.004 0.006 0.008 0.01
Time, s
-5
Current, A
-10
-15
-20
33
Copyright © 2016 John Wiley & Sons
Chapter 7
Using R = 0.25 deg C-s/J (cardboard), C = 1237 J/deg C, and Ta = –10 deg C the model becomes
Clearly the time constant is = RC = 309.25 s and therefore the settling time is tS = 4 = 1237 s.
The steady-state temperature is Tss = Ta = –10 deg C (ambient temperature). Therefore the
temperature response T(t) starts at 80 deg C and exponentially decays to –10 deg C in 1237 s
(20.6 min). A hand-sketch would match the below plot.
100
90
80
Coffee temperature with cardboard insualtion, deg C
70
60
50
40
30
20
10
-10
-20
0 500 1000 1500
Time, s
b) Now, R = 5.5 deg C-s/J (thermos). The new time constant is = RC = 6803.5 s and
therefore the settling time is tS = 4 = 27,214 s (7.6 hr). The initial and final temperatures are
the same; only the settling time has changed. A hand-sketch would match the plot below.
100
90
80
Coffee temperature with thermos insualtion, deg C
70
60
50
40
30
20
10
-10
-20
0 0.5 1 1.5 2 2.5 3
Time, s 4
x 10
34
Copyright © 2016 John Wiley & Sons
Chapter 7
Clearly the time constant is = 0.003 s and therefore the settling time is tS = 4 = 0.012 s. The
DC gain of the amplifier transfer function is 25 and therefore the steady-state voltage for a 0.2 V
step input is 5 V. The voltage response eO(t) will exponentially reach 5 V in 0.012 s. A hand-
sketch would match the below plot.
4
Amplifier output voltage, V
0
0 0.005 0.01 0.015
Time, s
b) We can compute the steady-state spool-valve position from the individual DC gains of the
three transfer functions: power amp DC gain = 25, solenoid DC gain = 1.6, and spool-valve
DC gain = 1/1800 = 5.5556(10–4). The product of the three DC gains is 0.0222. Hence the
steady-state valve position is the product of the total DC gain (0.0222) and magnitude of the
step input voltage, ein(t) = 0.2 V, or, (0.0222)(0.2) = 0.0044 m = 4.4 mm.
The undamped natural frequency is n = 226.7787 rad/s and the damping ratio is = 200/(2n) =
0.441 and hence the spool valve step response is a damped sinusoidal that eventually settles at the
steady-state value of 0.0044 m (4.4 mm) computed in part (b). The transient step-response
characteristics can be computed using the performance equations in Table 7-4:
Peak time: t p = 0.0154 s
n 1 2
35
Copyright © 2016 John Wiley & Sons
Chapter 7
/ 1 2
Maximum overshoot: M OS e = 0.2136, or zmax = (1.2136)(0.0044) = 0.0053 m
4
Settling time: tS = 0.04 s
n
2
Period of oscillation: Tperiod = 0.0309 s
n 1 2
2 1 2
Number of cycles to steady state: N cycles ~ 1.3 cycles
A hand-sketch would show (with labels) the above calculated values and match the below plot.
-3
x 10
6
Valve position, z, m (approximate sketch)
0
0 0.01 0.02 0.03 0.04 0.05 0.06
Time, s
Finally, we can verify the approximate hand-sketch (i.e., step response of the 2nd-order spool
valve model) by a numerical simulation of the complete system using Simulink (next page):
36
Copyright © 2016 John Wiley & Sons
Chapter 7
The step response of the complete Simulink model is below. Note that the steady-state value
(0.0044 m) matches the approximate sketch. However, the addition of two first-order models
(power amp and solenoid) has delayed the response. Note that the amplifier and solenoid have
time constants of 0.003 and 0.002 s, respectively, which correspond to settling times around
0.012 and 0.008 s. Therefore, the second-order spool valve response is delayed (compare the
peak times between the approximate sketch and the full simulation). The peak response of the
full simulation is also diminished (damped) compared to the approximate sketch.
-3
x 10
6
Valve position, z, m (Simulink result)
0
0 0.01 0.02 0.03 0.04 0.05 0.06
Time, s
37
Copyright © 2016 John Wiley & Sons
Chapter 7
7.32 The measured motor characteristics are a) steady-state ss = 85 rad/s for a 2-V step
voltage input, b) settling time tS = 0.6 s, and c) the speed response from zero to 85 rad/s is an
exponential rise.
Because the response is exponential we can model the system using the “standard form” of a first-
order transfer function:
(s) a
G (s)
Ein ( s ) s 1
The time constant is = tS/4 = 0.15 s and the DC gain is G(0) = a = (85 rad/s)/(2 V) = 42.5.
Hence an appropriate transfer function model of the DC motor is
( s) 42.5
G ( s)
Ein ( s ) 0.15s 1
The following MATLAB commands simulate the step response of this transfer function model.
The plot of the simulation results (below) verifies the experimental characteristics of the DC
motor.
90
80
70
DC motor speed, rad/s
60
50
40
30
20
10
0
0 0.2 0.4 0.6 0.8 1
Time, s
38
Copyright © 2016 John Wiley & Sons
Chapter 7
b k 1 1
x x x f a or x 2 n x n2 x fa
m m m m
The mass is m = 3 kg. Maximum overshoot for a step input only depends on damping ratio
(see Table 7-4). Substitute each value of b, k and compute damping ratio:
b 90
Option A: = 0.2444
2 mk 2 (3)(11,300)
b 88
Option B: = 0.2724
2 mk 2 (3)(8,700)
Since Option B has the largest damping ratio it has the smallest percent overshoot to a step input.
39
Copyright © 2016 John Wiley & Sons
Chapter 7
7.34 a) The measured output shows a classic underdamped response which we can model
with a second-order transfer function:
( s ) an2
G ( s) 2
Tin ( s) s 2 n s n2
The steady-state output is 0.012 rad for a step input torque of 0.6 N-m. Hence the DC gain of the
transfer function is a = 0.012/0.6 = 0.02.
The peak response is 0.016 rad which is a 33.333% overshoot of the steady-state value of 0.012
rad. The maximum overshoot is
M OS e / 1 2
= 0.333333 damping ratio = 0.3301
The period of the response is 1.1 s. Therefore the damped frequency is d 2 / 1.1 = 5.712
rad/s. The undamped natural frequency is n d / 1 = 6.0512 rad/s. The second-order
2
transfer function is
( s ) 0.73234
G(s) 2
Tin ( s ) s 3.995s 36.617
b) The following MATLAB commands simulate the step response of this transfer function
model.
The simulation result (plotted below) shows a good match with Fig. P7.34; therefore the second-
order transfer function G(s) is an accurate representation of the mechanical system.
0.02
0.018
0.016
Angular displacement, rad
0.014
0.012
0.01
0.008
0.006
0.004
0.002
0
0 0.5 1 1.5 2 2.5 3
Time, s
40
Copyright © 2016 John Wiley & Sons