Repeat in Chapter 1, using Laplace transforms. Assume that

Derive the Laplace transform for the following time functions: [Section: 2.2] a. u(t) b. tu(t) c. \(\sin \omega t u(t)\) d. \(\cos \omega t u(t)\)

Using the Laplace transform pairs of Table 2.1 and the Laplace transform theorems of Table 2.2, derive the Laplace transforms for the following time functions: [Section: 2.2] a. \(e^{-a t} \sin \omega t u(t)\) b. \(e^{-a t} \cos \omega t u(t)\) c. \(t^3 u(t)\)

Repeat Problem 19 in Chapter 1, using Laplace transforms. Assume zero initial conditions. [Sections: 2.2 ; 2.3]

Repeat Problem 20 in Chapter 1, using Laplace transforms. Assume that the forcing functions are zero prior to t=0-. [Section: 2.2]

Repeat Problem 21 in Chapter 1, using Laplace transforms. Use the following initial conditions for each part as follows: (a) \(x(0)=4, x^{\prime}(0)=-4\); (b) \(x(0)=4, x^{\prime}(0)=1\); (c) \(x(0)=2, x^{\prime}(0)=3\), where \(x^{\prime}(0)=\frac{d x}{d t}(0)\). Assume that the forcing functions are zero prior to t = 0-. [Section: 2.2]

Use MATLAB and the Symbolic Math Toolbox to find the Laplace transform of the following time functions: [Section: 2.2] a. \(f(t)=8 t^2 \cos \left(3 t+45^{\circ}\right)\) b. \(f(t)=3 t e^{-2 t} \sin \left(4 t+60^{\circ}\right)\)

Use MATLAB and the Symbolic Math Toolbox to find the inverse Laplace transform of the following frequency functions: [Section: 2.2] a. \(G(s)=\frac{\left(s^2+3 s+10\right)(s+5)}{(s+3)(s+4)\left(s^2+2 s+100\right)}\) b. \(G(s)=\frac{s^3+4 s^2+2 s+6}{(s+8)\left(s^2+8 s+3\right)\left(s^2+5 s+7\right)}\)

A system is described by the following differential equation: \(\frac{d^3 y}{d t^3}+3 \frac{d^2 y}{d t^2}+5 \frac{d y}{d t}+y=\frac{d^3 x}{d t^3}+4 \frac{d^2 x}{d t^2}+6 \frac{d x}{d t}+8 x\) Find the expression for the transfer function of the system, Y(s)/X(s). [Section: 2.3]

For each of the following transfer functions, write the corresponding differential equation. [Section: 2.3] a. \(\frac{X(s)}{F(s)}=\frac{7}{s^2+5 s+10}\) b. \(\frac{X(s)}{F(s)}=\frac{15}{(s+10)(s+11)}\) c. \(\frac{X(s)}{F(s)}=\frac{s+3}{s^3+11 s^2+12 s+18}\)

Write the differential equation for the system shown in Figure P2.1. [Section: 2.3]

Write the differential equation that is mathematically equivalent to the block diagram shown in Figure P2.2. Assume that \(r(t)=3 t^3\). [Section: 2.3]

A system is described by the following differential equation: [Section 2.3] \(\frac{d^2 x}{d t^2}+4 \frac{d x}{d t}+5 x=1\) with the initial conditions \(x(0)=1, \dot{x}(0)=-1\). Show a block diagram of the system, giving its transfer function and all pertinent inputs and outputs. (Hint: the initial conditions will show up as added inputs to an effective system with zero initial conditions.)

Use MATLAB to generate the transfer function: [Section: 2.3] \(G(s)=\frac{5(s+15)(s+26)(s+72)}{s(s+55)\left(s^2+5 s+30\right)(s+56)\left(s^2+27 s+52\right)}\) in the following ways: a. the ratio of factors; b. the ratio of polynomials.

Repeat Problem 13 for the following transfer function: [Section: 2.3] \(G(s)=\frac{s^4+25 s^3+20 s^2+15 s+42}{s^5+13 s^4+9 s^3+37 s^2+35 s+50}\)

Use MATLAB to generate the partial fraction expansion of the following function: [Section: 2.3] \(F(s)=\frac{10^4(s+5)(s+70)}{s(s+45)(s+55)\left(s^2+7 s+110\right)\left(s^2+6 s+95\right)}\)

Use MATLAB and the Symbolic Math Toolbox to input and form LTI objects in polynomial and factored form for the following frequency functions: [Section: 2.3] a. \(G(s)=\frac{45\left(s^2+37 s+74\right)\left(s^3+28 s^2+32 s+16\right)}{(s+39)(s+47)\left(s^2+2 s+100\right)\left(s^3+27 s^2+18 s+15\right)}\) b. \(G(s)=\frac{56(s+14)\left(s^3+49 s^2+62 s+53\right)}{\left(s^3+81 s^2+76 s+65\right)\left(s^2+88 s+33\right)\left(s^2+56 s+77\right)}\)

Find the transfer function, \(G(s)=V_o(s) / V_i(s)\), for each network shown in Figure P2.3. [Section: 2.4]

Find the transfer function, \(G(s)=V_L(s) / V(s)\), for each network shown in Figure P2.4. [Section: 2.4]

Find the transfer function, \(G(s)=V_o(s) / V_i(s)\), for each network shown in Figure P2.5. Solve the problem using mesh analysis. [Section: 2.4]

Repeat Problem 19 using nodal equations. [Section: 2.4]

(a) Write, but do not solve, the mesh and nodal equations for the network of Figure P2.6. [Section: 2.4] (b) Use MATLAB, the Symbolic Math Toolbox, and the equations found in part a to solve for the transfer function, \(G(s)=V_o(s) / V(s)\). Use both the mesh and nodal equations and show that either set yields the same transfer function. [Section: 2.4]

Find the transfer function, \(G(s)=V_o(s) / V_i(s)\), for each operational amplifier circuit shown in Figure P2.7. [Section: 2.4]

Find the transfer function, \(G(s)=V_o(s) / V_i(s)\), for each operational amplifier circuit shown in Figure P2.8. [Section: 2.4]

Find the transfer function, \(G(s)=X_1(s) / F(s)\), for the translational mechanical system shown in Figure P2.9. [Section: 2.5]

Find the transfer function, \(G(s)=X_2(s) / F(s)\), for the translational mechanical network shown in Figure P2.10. [Section: 2.5]

Find the transfer function, \(G(s)=X_2(s) / F(s)\), for the translational mechanical system shown in Figure P2.11. (Hint: place a zero mass at \(x_2(t)\).) [Section: 2.5]

For the system of Figure P2.12 find the transfer function, \(G(s)=X_1(s) / F(s)\). [Section: 2.5]

Find the transfer function, \(G(s)=X_3(s) / F(s)\), for the translational mechanical system shown in Figure P2.13. [Section: 2.5]

Find the transfer function, \(X_3(s) / F(s)\), for each system shown in Figure P2.14. [Section: 2.5]

Write, but do not solve, the equations of motion for the translational mechanical system shown in Figure P2.15. [Section: 2.5]

For the unexcited (no external force applied) system of Figure P2.16, do the following: (a) Write the differential equation that describes the system. (b) Assuming initial conditions \(x(0)=x_0\) and \(\dot{x}(0)=x_1\), write a Laplace transform expression for X(s). (c) Find x(t) by obtaining the inverse Laplace transform from the result in Part c. d. What will be the oscillation frequency in \(\mathrm {H z}\) for this system?

For each of the rotational mechanical systems shown in Figure P2.17, write, but do not solve, the equations of motion. [Section: 2.6]

For the rotational mechanical system shown in Figure P2.18, find the transfer function \(G(s)= \theta_2(s) / T(s)\) [Section: 2.6]

Find the transfer function, \(\frac{\theta_1(s)}{T(s)}\), for the system shown in Figure P2.19.

For the rotational mechanical system with gears shown in Figure P2.20, find the transfer function, \(G(s)=\theta_3(s) / T(s)\). The gears have inertia and bearing friction as shown. [Section: 2.7]

For the rotational system shown in Figure P2.21, find the transfer function, \(G(s)=\theta_2(s) / T(s)\). [Section: 2.7]

Find the transfer function, \(G(s)=\theta_2(s) / T(s)\), for the rotational mechanical system shown in Figure P2.22. [Section: 2.7]

Find the transfer function, \(G(s)=\theta_4(s) / T(s)\), for the rotational system shown in Figure P2.23. [Section: 2.7]

For the rotational system shown in Figure P2.24, find the transfer function, \(G(s)=\theta_L(s) / T(s)\). [Section: 2.7]

For the rotational system shown in Figure P2.25, write the equations of motion from which the transfer function, \(G(s)=\theta_1(s) / T(s)\), can be found. [Section: 2.7]

Given the rotational system shown in Figure P2.26, find the transfer function, \(G(s)=\theta_6(s) / \theta_1(s)\). [Section: 2.7]

In the system shown in Figure P2.27, the inertia, J, of radius, r, is constrained to move only about the stationary axis A. A viscous damping force of translational value \(f_v\) exists between the bodies J and M. If an external force, f(t), is applied to the mass, find the transfer function, \(G(s)=\theta(s) / F(s)\). [Sections: 2.5; 2.6]

For the combined translational and rotational system shown in Figure P2.28, find the transfer function, G(s)=X(s)/T(s). [Sections: 2.5; 2.6; 2.7]

Given the combined translational and rotational system shown in Figure P2.29, find the transfer function, G(s)=X(s)/T(s). [Sections: 2.5 ; 2.6 ]

For the motor, load, and torque-speed curve shown in Figure P2.30, find the transfer function, \(G(s)=\theta_L(s) / E_a(s)\). [Section: 2.8]

The motor whose torque-speed characteristics are shown in Figure P2.31 drives the load shown in the diagram. Some of the gears have inertia. Find the transfer function, \(G(s)=\theta_2(s) / E_a(s)\). [Section: 2.8]

A dc motor develops \(55 \mathrm{~N}-\mathrm{m}\) of torque at a speed of \(600 \mathrm{rad} / \mathrm{s}\) when 12 volts are applied. It stalls out at this voltage with \(100 \mathrm{~N}-\mathrm{m}\) of torque. If the inertia and damping of the armature are \(7 \mathrm{~kg}-\mathrm{m}^2\) and \(3 \mathrm{~N}-\mathrm{m}-\mathrm{s} / \mathrm{rad}\), respectively, find the transfer function, \(G(s)=\theta_L(s) / E_a(s)\), of this motor if it drives an inertia load of \(105 \mathrm{~kg}-\mathrm{m}^2\) through a gear train, as shown in Figure P2.32. [Section: 2.8]

In this chapter, we derived the transfer function of a dc motor relating the angular displacement output to the armature voltage input. Often we want to control the output torque rather than the displacement. Derive the transfer function of motor that relates output torque to input armature voltage. [Section: 2.8]

Find the transfer function, \(G(s)=X(s) / E_a(s)\), for the system shown in Figure P2.33. [Sections: 2.5-2.8]

Find the series and parallel analogs for the translational mechanical system shown in Figure 2.20 in the text. [Section: 2.9]

Find the series and parallel analogs for the rotational mechanical systems shown in Figure P2.17(b) in the problems. [Section: 2.9]

A system's output, c, is related to the system's input, r, by the straight-line relationship, c=5r+7 Is the system linear? [Section: 2.10]

Consider the differential equation \(\frac{d^2 x}{d t^2}+3 \frac{d x}{d t}+2 x=f(x)\) where f(x) is the input and is a function of the output, x. If f(x)=sin x, linearize the differential equation for small excursions. [Section: 2.10] a. x=0 b. \(x=\pi\)

Consider the differential equation \(\frac{d^3 x}{d t^3}+10 \frac{d^2 x}{d t^2}+20 \frac{d x}{d t}+15 x=f(x)\) where f(x) is the input and is a function of the output, x. If \(f(x)=3 e^{-5 x}\), linearize the differential equation for x near 0. [Section: 2.10]

Many systems are piecewise linear. That is, over a large range of variable values, the system can be described linearly. A system with amplifier saturation is one such example. Given the differential equation \(\frac{d^2 x}{d t^2}+17 \frac{d x}{d t}+50 x=f(x)\) assume that f(x) is as shown in Figure P2.34. Write the differential equation for each of the following ranges of x : [Section: 2.10] a. \(-\infty<x<-3\) b. -3<x<3 c. \(3<x<\infty\)

For the translational mechanical system with a nonlinear spring shown in Figure P2.35, find the transfer function, G(s)=X(s)/F(s), for small excursions around f(t)=1. The spring is defined by \(x_s(t)=1-e^{-f_s(t)}\), where \(x_s(t)\) is the spring displacement and \(f_s(t)\) is the spring force. [Section: 2.10]

Enzymes are large proteins that biological systems use to increase the rate at which reactions occur.For Example, food is usually composed of large molecules that are hard to digest; enzymes break down the large molecules into small nutrient as part of the digestive process. One Such Enzyme Is amylase, contained in human saliva. It is commonly known that if you place a piece of uncooked pasta in your mouth its taste will change from paper-like to sweet as amylase breaks down the carbohydrates into sugars. Enzyme breakdown is often expressed by the following relation: \(S+E \underset{k_1}{\stackrel{k_s}{\rightleftarrows}} C \stackrel{k_1}{\rightarrow} P\) In this expression a substrate (S) interacts with an enzyme (E) to form a combined product (C) at a rate \(k_1\). The intermediate compound is reversible and gets disassociated at a rate \(k_{-1}\). Simultaneously some of the compound is transformed into the final product (P) at a rate \(k_2\). The kinetics describing this reaction are known as the Michaelis-Menten equations and consist of four nonlinear differential equations. However, under some conditions these equations can be simplified. Let \(E_0\) and \(S_0\) be the initial concentrations of enzyme and substrate, respectively. It is generally accepted that under some energetic conditions or when the enzyme concentration is very big \(\left(E_0 \gg S_0\right)\), the kinetics for this reaction are given by \(\begin{aligned} \frac{d S}{d t} & =k_\psi\left(\tilde{K}_s C-S\right) \\ \frac{d C}{d t} & =k_\psi\left(S-\tilde{K}_M C\right) \\ \frac{d P}{d t} & =k_2 C \end{aligned}\) where the following constant terms are used (Schnell, 2004): \(\begin{aligned} k_\psi & =k_1 E_0 \\ \tilde{K}_s & =\frac{k-1}{k_\psi} \end{aligned}\) and \(\tilde{K}_M=\tilde{K}_s+\frac{k_2}{k_\psi}\) a. Assuming the initial conditions for the reaction are \(S(0)=S_0, E(0)=E_0, C(0)=P(0)=0\), find the Laplace transform expressions for S, C, and P : \(\mathscr{L}\{S\}, \mathscr{L}\{C\}\), and \(\mathscr{L}\{P\}\), respectively. b. Use the final theorem to find \(S(\infty), C(\infty)\), and \(P(\infty)\).

Humans are able to stand on two legs through a complex feedback system that includes several sensory inputs— equilibrium and visual along with muscle actuation. In order to gain a better understanding of the workings of the postural feedback mechanism, an individual is asked to stand on a platform to which sensors are attached at the base. Vibration actuators are attached with straps to the individual’s calves. As the vibration actuators are stimulated, the individual sways and movements are recorded. It was hypothesized that the human postural dynamics are analogous to those of a cart with a balancing standing pole attached (inverted pendulum). In that case, the dynamics can be described by the following two equations: \(\begin{aligned} J \frac{d^2 \theta}{d t^2} & =m g l \sin \theta(t)+T_{\text {bal }}+T_d(t) \\ T_{\text {bal }}(t) & =-m g l \sin \theta(t)+k J \theta(t)-\eta J \dot{\theta}(t) \\ & -\rho J \int_0^t \theta(t) d t \end{aligned}\) where m is the individual's mass; l is the height of the individual's center of gravity; g is the gravitational constant; J is the individual's equivalent moment of inertia; \(\eta, \rho\), and k are constants given by the body's postural control system; \(\theta(t)\) is the individual's angle with respect to a vertical line; \(T_{\text {bal }}(t)\) is the torque generated by the body muscles to maintain balance; and \(T_d(t)\) is the external torque input disturbance. Find the transfer function \(\frac{\Theta(s)}{T_d(s)}\) (Johansson, 1988).

Figure P2.36 shows a crane hoisting a load. Although the actual system’s model is highly nonlinear, if the rope is considered to be stiff with a fixed length L, the system can be modeled using the following equations: \(\begin{aligned} m_L \ddot{x}_{L a} & =m_L g \phi \\ m_T \ddot{x}_T & =f_T-m_L g \phi \\ x_{L a} & =x_T-x_L \\ x_L & =L \phi \end{aligned}\) where \(m_L\) is the mass of the load, \(m_T\) is the mass of the cart, \(x_T\) and \(x_L\) are displacements as defined in the figure, \(\phi\) is the rope angle with respect to the vertical, and \(f_T\) is the force applied to the cart (Marttinen, 1990). (a) Obtain the transfer function from cart velocity to rope angle \(\frac{\Phi(s)}{V_T(s)}\) (b) Assume that the cart is driven at a constant velocity \(V_0\) and obtain an expression for the resulting \(\phi(t)\). Show that under this condition, the load will sway with a frequency \(\omega_0=\sqrt{\frac{g}{L}}\). (c) Find the transfer function from the applied force to the cart's position, \(\frac{X_T(s)}{F_T(s)}\). (d) Show that if a constant force is applied to the cart, its velocity will increase without bound as \( \rightarrow \infty\).

In 1978, Malthus developed a model for human growth population that is also commonly used to model bacterial growth as follows. Let N(t) be the population density observed at time t. Let K be the rate of reproduction per unit time. Neglecting population deaths, the population density at a time \(t+\Delta t\) (with small \(\Delta t\) ) is given by \(N(t+\Delta t) \approx N(t)+K N(t) \Delta t\) which also can be written as \(\frac{N(t+\Delta t)-N(t)}{\Delta t}=K N(t)\) Since N(t) can be considered to be a very large number, letting \(\Delta t \rightarrow 0\) gives the following differential equation (Edelstein-Keshet, 2005): \(\frac{d N(t)}{d t}=K N(t)\) a. Assuming an initial population \(N(0)=N_0\), solve the differential equation by finding N(t). b. Find the time at which the population is double the initial population.

In order to design an underwater vehicle that has the characteristics of both a long-range transit vehicle (torpedo-like) and a highly maneuverable low-speed vehicle (boxlike), researchers have developed a thruster that mimics that of squid jet locomotion (Krieg, 2008). It has been demonstrated there that the average normalized thrust due to a command step input, \(U(s)=\frac{T_{r e f}}{s}\) is given by: \(T(t)=T_{r e f}\left(1-e^{-\lambda t}\right)+a \sin (2 \pi f t)\) where \(T_{\text {ref }}\) is the reference or desired thrust, \(\lambda\) is the system's damping constant, a is the amplitude of the oscillation caused by the pumping action of the actuator, f is the actuator frequency, and T(t) is the average resulting normalized thrust. Find the thruster's transfer function \(\frac{T(s)}{U(s)}\). Show all steps.

The Gompertz growth model is commonly used to model tumor cell growth. Let v(t) be the tumor's volume, then \(\frac{d v(t)}{d t}=\lambda e^{-\alpha t} v(t)\) where \(\lambda\) and \(\alpha\) are two appropriate constants (EdelsteinKeshet, 2005). (a) Verify that the solution to this equation is given by \(v(t)=v_0 e^{\lambda / \alpha\left(1-e^{-\alpha t}\right)}\), where \(v_0\) is the initial tumor volume. (b) This model takes into account the fact that when nutrients and oxygen are scarce at the tumor's core, its growth is impaired. Find the final predicted tumor volume (let \(t \rightarrow \infty\) ). (c) For a specific mouse tumor, it was experimentally found that \(\lambda=2.5\) days, \(\alpha=0.1\) days with \(v_0= 50 \times 10^{-3} \mathrm{~mm}^3\) (Chignola, 2005). Use any method available to make a plot of v(t) vs. t. (d) Check the result obtained in Part b with the results from the graph found in Part c. A muscle hanging from a beam is shown in Figure P2.37(a) (Lessard, 2009). The \(\alpha\)-motor neuron can be used to electrically stimulate the muscle to contract and pull the mass, m, which under static conditions causes the muscle to stretch. An equivalent mechanical system to this setup is shown in Figure P2.37(b). The force \(F_{i s o}\) will be exerted when the muscle contracts. Find an expression for the displacement \(X_1(s)\) in terms of \(F_1(s)\) and \(F_{\text {iso }}(s)\).

A muscle hanging from a beam is shown in Figure P2.37(a) (Lessard, 2009). The \(\alpha\)-motor neuron can be used to electrically stimulate the muscle to contract and pull the mass, m, which under static conditions causes the muscle to stretch. An equivalent mechanical system to this setup is shown in Figure P2.37(b). The force \(F_{i s o}\) will be exerted when the muscle contracts. Find an expression for the displacement \(X_1(s)\) in terms of \(F_1(s)\) and \(F_{\text {iso }}(s)\).

A three-phase ac/dc converter supplies dc to a battery charging system or dc motor (Graovac, 2001). Each phase has an ac filter represented by the equivalent circuit in Figure P2.38. Derive that the inductor current in terms of the two active sources is \(I_{a c F}(s)=\frac{1+R C s}{L C s^2+R C s+1} I_{a c R}(s)+\frac{C s}{L C s^2+R C s+1} V_a(s)\)

A photovoltaic system is used to capture solar energy to be converted to electrical energy. A control system is used to pivot the solar platform to track the sun's movements in order to maximize the captured energy. The system consists of a motor and load similar to that discussed in Section 2.8. A model has been proposed (Agee, 2012) that is different from the model developed in the chapter in the following ways: (1) the motor inductance was not neglected and (2) the load, in addition to having inertia and damping, has a spring. Find the transfer function, \(\theta_m(s) / E_a(s)\), for this augmented system assuming all load impedances have already been reflected to the motor shaft.

In a paint mixing plant, two tanks supply fluids to a mixing cistern. The height, h, of the fluid in the cistern is dependent upon the difference between the input mass flow rate, q, and the output flow rate, \(q_e\). A nonlinear differential equation describing this dependency is given by (Schiop, 2010) \(\frac{d h}{d t}+\frac{A_e}{A} \sqrt{2 g h}=\frac{q}{\rho A}\) where A= cross-sectional area of the cistern, \(A_e=\) cross-sectional area of the exit pipe, g= acceleration due to gravity, and \(\rho=\) liquid density. (a) Linearize the nonlinear equation about the equilibrium point \(\left(h_0, q_0\right)\) and find the transfer function relating the output cistern fluid level, H(s), to the input mass flow rate, Q(s). (b) The color of the liquid in the cistern can be kept constant by adjusting the input flow rate, q, assuming the input flow's color is specifically controlled. Assuming an average height, \(h_{a v}\), of the liquid in the cistern, the following equation relates the net flow of color to the cistern to the color in the cistern. \(e_1 q-e q_e=\frac{d}{d t}\left(\rho A e h_{a v}\right)\) where \(e_1=\) fractional part of flow representing color into the cistern, and e= fractional part of the cistern representing color in the cistern. Assume that the flow out of the cistern is constant and use the relationship, \(q_e=\rho A_e \sqrt{2 g h_{a v}}\), along with the given equation above to find the transfer function, E(s)/Q(s), that relates the color in the cistern to the input flow rate.

Control of HIV/AIDS. HIV inflicts its damage by infecting healthy CD4 + T cells (a type of white blood cell) that are necessary to fight infection. As the virus embeds in a T cell and the immune system produces more of these cells to fight the infection, the virus propagates in an opportunistic fashion. As we now develop a simple HIV model, refer to Figure P2.39. Normally T cells are produced at a rate s and die at a rate d. The HIV virus is present in the bloodstream in the infected individual. These viruses in the bloodstream, called free viruses, infect healthy T cells at a rate \(\beta\). Also, the viruses reproduce through the T cell multiplication process or otherwise at a rate k. Free viruses die at a rate c. Infected T cells die at a rate \(\mu\). A simple mathematical model that illustrates these interactions is given by the following equations (Craig, 2004): \(\begin{aligned} \frac{d T}{d t} & =s-d T-\beta T v \\ \frac{d T^*}{d t} & =\beta T v-\mu T^* \\ \frac{d v}{d t} & =k T^*-c v \end{aligned}\) where T = number of healthy T cells \(T^*\) = number of infected T cells v = number of free viruses a. The system is nonlinear; thus linearization is necessary to find transfer functions as you will do in subsequent chapters. The nonlinear nature of this model can be seen from the above equations. Determine which of these equations are linear, which are nonlinear, and explain why. b. The system has two equilibrium points. Show that these are given by \(\left(T_0, T_0^*, v_0\right)=\left(\frac{s}{d}, 0,0\right)\) and \(\left(T_0, T_0^*, v_0\right)=\left(\frac{c \mu}{\beta k}, \frac{s}{\mu}-\frac{c d}{\beta k}, \frac{s k}{c \mu}-\frac{d}{\beta}\right)\)

Hybrid vehicle. Problem 23 in Chapter 1 discusses the cruise control of serial, parallel, and split-power hybrid electric vehicles (HEVs). The functional block diagrams developed for these HEVs indicated that the speed of a vehicle depends upon the balance between the motive forces (developed by the gasoline engine and/or the electric motor) and running resistive forces. The resistive forces include the aerodynamic drag, rolling resistance, and climbing resistance. Figure P2.40 illustrates the running resistances for a car moving uphill (Bosch, 2007). The total running resistance, \(F_w\), is calculated as \(F_w=F_{R o}+F_L+F_{S t}\), where \(F_{R o}\) is the rolling resistance, \(F_L\) is the aerodynamic drag, and \(F_{S t}\) is the climbing resistance. The aerodynamic drag is proportional to the square of the sum of car velocity, v, and the head-wind velocity, \(v_{h w}\), or \(v+v_{h w}\). The other two resistances are functions of car weight, G, and the gradient of the road (given by the gradient angle \(\alpha\) ), as seen from the following equations: \(F_{R o}=f G \cos \alpha=f m g \cos \alpha\) where \(\begin{aligned} f & =\text { coefficient of rolling resistance } \\ m & =\text { car mass, in } \mathrm{kg} \\ g & =\text { gravitational acceleration, in } \mathrm{m} / \mathrm{s}^2 \\ F_L & =0.5 \rho C_w A\left(v+v_{\mathrm{hw}}\right)^2 . \end{aligned}\) and \(\begin{aligned} \rho & =\text { air density, in } \mathrm{kg} / \mathrm{m}^3 \\ C_w & =\text { coefficient of aerodynamic drag } \\ A & =\text { largest cross-section of the car, in } \mathrm{kg} / \mathrm{m}^2 \\ F_{S t} & =G \sin \alpha=m g \sin \alpha . \end{aligned}\) The motive force, F, available at the drive wheels is: \(F=\frac{T i_{\text {tot }}}{r} \eta_{t o t}=\frac{P \eta_{t o t}}{v}\) where \(\begin{aligned} T & =\text { motive torque } \\ P & =\text { motive power } \\ i_{\text {tot }} & =\text { total transmission ratio } \\ r & =\text { tire radius } \\ \eta_{\text {tot }} & =\text { total drive-train efficiency. } \end{aligned}\) The surplus force, \(F-F_w\), accelerates the vehicle (or retards it when \(F_w>F\) ). Letting \(a=\frac{F-F_w}{k_m \cdot m}\), where a is the acceleration and \(k_m\) is a coefficient that compensates for the apparent increase in vehicle mass due to rotating masses (wheels, flywheel, crankshaft, etc.): (a) Show that car acceleration, \({ }^{21} a\), may be determined from the equation: \(F=f m g \cos \alpha+m g \sin \alpha+0.5 \rho C_{\mathrm{w}} A\left(v+v_{h w}\right)^2+k_m m a\) (b) Assuming constant acceleration and using the average value for speed, find the average motive force, \(F_{a v}\) (in N), and power, \(P_{a v}\) (in \(\mathrm{kW}\) ) the car needs to accelerate from 40 to \(60 \mathrm{~km} / \mathrm{h}\) in 4 seconds on a level road, \(\left(\alpha=0^{\circ}\right)\), under windless conditions, where \(v_{h w}=0\). You are given the following parameters: \(\mathrm{m}=1590 \mathrm{~kg}, \mathrm{~A}=2 \mathrm{~m}^2, f=0.011, \rho=1.2 \mathrm{~kg} / \mathrm{m}^3\), \(C_w=0.3, \eta_{\text {tot }}=0.9, k_m=1.2\). Furthermore, calculate the additional power, \(P_{\text {add }}\), the car needs after reaching \(60 \mathrm{~km} / \mathrm{h}\) to maintain its speed while climbing a hill with a gradient \(\alpha=5^{\circ}\). (c) The equation derived in Part a describes the nonlinear car motion dynamics where F(t) is the input to the system, and v(t) the resulting output. Given that the aerodynamic drag is proportional to \(v^2\) under windless conditions, linearize the resulting equation of motion around an average speed, \(v_{\mathrm{o}}=50 \mathrm{~km} / \mathrm{h}\), when the car travels on a level road, \({ }^{22}\) where \(\alpha=0^{\circ}\). (Hint: Expand \(v^2-v_0^2\) in a truncated Taylor series). Write that equation of motion and represent it with a block diagram in which the block \(G_v\) represents the vehicle dynamics. The output of that block is the car speed, v(t), and the input is the excess motive force, \(F_e(t)\), defined as: \(F_e=F-F_{S t}-F_{R o}+F_o\), where \(F_o\) is the constant component of the linearized aerodynamic drag. (d) Use the equation in Part c to find the vehicle transfer function: \(G_v(s)=V(s) / F_e(s)\).

Parabolic trough collector. In a significant number of cases, the open-loop transfer function from fluid flow to fluid temperature in a parabolic trough collector can be approximated (Camacho, 2012) by: \(P(s)=\frac{K}{1+\tau s)} e^{-s T}\) (a) Write an analytic expression for the unit step response of the open-loop system assuming that h(t) represents the output temperature and q(t) the input fluid flow. (b) Make a sketch of the unit step response of the open-loop system. Indicate on your figure the time delay, the settling time, the initial and final values of the response, and the value of the response when \(t=\tau+T\). (c) Call the output temperature h(t) and the input fluid flow q(t). Find the differential equation that represents the open-loop system.