Name ME6 Final. I certify that I upheld the Stanford Honor code during this exam Monday December 2, 2005 3:30-6:30 p.m. ffl Print your name and sign the honor code statement ffl You may use your course notes, homework, books, etc. ffl Write your answers on this handout ffl Where space is provided, show your work to get credit ffl If necessary, attach extra pages for scratch work ffl Best wishes for a fun vacation. Merry Christmas and Happy New Year! Page Value Score 2 2 3 3 4 7 5 4 6 4 7 0 8 2 9 7 Total 00 F. (2 pts.) Miscellaneous (a) ( pt.) What is your Meyer-Briggs (human metrics) personality type? (b) ( pt.) Suppose you plan on teaching or TAing Dynamic Systems. Provide suggestions for improving the course. Specifically, what you would add, remove, or change in the course (e.g., content, homework, book, lectures, etc.)?
F.2 (3 pts.) Coupled motions of Wilbur-force pendulum N N o Shown to the right is a rigid body B that is attached to a spring at B c, the mass center of B. The other end of the spring is attached to point N o which is fixed in a Newtonian reference frame N. The equations governing this model of the system are B m ẍ + k x x + k c = 0 I + k c x + k = 0 B c Quantity Symbol Type Mass of B m constant Central moment of inertia of B about vertical I constant Linear spring constant modeling extensional flexibility k x constant Linear spring constant modeling torsional flexibility k constant Linear spring constant modeling coupled flexibility k c constant Translational stretch of spring from equilibrium x dependent variable Rotational stretch of spring from equilibrium dependent variable (a) (2 pts.) Write the ODEs in the matrix form M Ẍ + B _X + KX = 0, where X = " # ẍ _x x 0 + + = _ 0 (b) (9 pts.) The solution to the previous set of ODEs has the form X(t)=U Λe pt where p is a constant (to-be-determined) and U is a non-zero 2 matrix of constants (to-be-determined). Fill in the two blanks in the following polynomial equation that governs the values of = - p 2. Note: the blanks only involve m, I, k x, k, k c. x. 2 + Λ + = 0 2
(c) (6 pts.) For certain values of m, I, k x, k, and k c, the matrix A = M - K = Find the eigenvalues and corresponding eigenvectors of A. 0 0. = U = 2 = U 2 = Imaginary Axis (d) (2 pts.) Calculate the values of p, p 2, p 3, p 4 and draw their locations in the complex plane. p ;2 = Real Axis p 3;4 = (e) ( pt.) The solution is stable/neutrally stable/unstable. (f) (4 pts.) Write the solution for X(t) in terms of the yet-to-be-determined constants c, c 2, c 3, c 4, the sine and cosine functions, and t. x(t) (t) = n c sin( t) + c 2 cos( t)o 3 + n c 3 sin( t) + c 4 cos( t)o
(g) (2 pts.) Determine c, c 2, c 3, c 4 when x(0) = 0:2, (0)=0, _x(0) = 0, and _ (0) = 0. c = c 2 = c 3 = c 4 = (h) (2 pts.) Using the aforementioned initial values, write explicit solutions for x(t) and (t). x(t) = (t) = (i) ( pt.) Using trigonometric identities, equation (2.9), and cos(a) =-cos(a+ß), show x(t) = 0:2 sin(-0:6 t + ß 2 ) sin(3:6 t + ß 2 ) = 0:2 cos(0:6 t) cos(3:6 t) 0.2 0.2 (t) = 0:2 sin(-0:6 t) sin(3:6 t) = 0:2 sin(0:6 t + ß) sin(3:6 t) 0.5 0.5 Displacement 0. 0.05 0-0.05-0. Rotation angle (radians) 0. 0.05 0-0.05-0. -0.5-0.5-0.2 0 5 0 5 20 25 30 35 40-0.2 0 5 0 5 20 25 30 35 40 (j) (2 pts.) Give an interpretation of the time-behavior of x(t) and (t). 4
F.3 (38 pts.) Linearization of equations of motion for a swinging spring A straight, massless, linear spring connects a particle Q to a point N o which is fixed in a Newtonian reference frame N. The system identifiers and governing equations are shown below. N Quantity Identifier Type Local gravitational constant g constant Mass of Q m constant Natural spring length L n constant Linear spring constant k constant Spring stretch y variable Pendulum angle" variable m (L n + y) + 2 m _ _y + mg sin( ) = 0 m ÿ m (L n +y) _ 2 + ky mg cos( ) = 0 N o b y k θ L n n z n y n x y b x (a) ( pt.) Classify the previous equations by picking the relevant qualifiers from the following list. Uncoupled Linear Homogeneous Constant-coefficient st-order Algebraic Coupled Nonlinear Inhomogeneous Variable-coefficient 2nd-order Differential (b) (3 pts.) Find the condition on y nom such that y = y nom (a constant) and = nom =0 satisfy the nonlinear ODEs. y nom = (c) (0 pts.) Using a Taylor series, linearize the differential equations in perturbations about this nominal solution. In other words, define perturbations of y and as ey = y y nom (y nom is constant) e = nom ( nom =0) and form a set of differential equations which arelineariney, _ ey, ëy, e, _ e, ë. = 0 5 = 0
(d) (7 pts.) Replace y and in the original nonlinear ODEs with y = ey + y nom (y nom is the constant determined earlier) = e + nom ( nom =0) Assuming ey and e are small and using the small-approximations technique of Chapter 20, form a set of ODEs that are linear in ey, _ ey, ëy, e, _ e, ë. = 0 = 0 (e) ( pt.) Classify your equations in part (3d) by picking the relevant qualifiers from the following list. Uncoupled Linear Homogeneous Constant-coefficient st-order Algebraic Coupled Nonlinear Inhomogeneous Variable-coefficient 2nd-order Differential (f) ( pt.) The linearized ODEs resulting from the Taylor series technique are identical to those obtained from the small-approximations technique. True/False. (g) (5 pts.) Find analytical solutions for ey(t) and e (t) when g =9:8 m=sec 2, m = kg, L n =0:5 m, k = 00 n=m, ey(0) = 0: m, e (0) = ffi,and _ ey(0) = _ e (0) = 0 ey(t) = e (t) = 6
(h) (4 pts.) Make a rough sketch of your linearized solutions for y(t)and (t)for 0 t 6 sec. y(t) ß (t) ß Spring stretch (meters) Pendulum angle (degrees) 0 2 4 6 8 0 2 4 6 0 2 4 6 8 0 2 4 6 (i) (6 pts.) The computer solution for ey and e to the original nonlinear ODEs for 0 t 6 sec are plotted below. Comment on the similarities and differences between the plots you produced with the ones below. Give a reason for the similarities and differences. 0.2 30 Spring stretch (meters) 0.5 0. 0.05 0 0 2 4 6 8 0 2 4 6 Similarities: Differences: Reason for differences: Pendulum angle (degrees) 20 0 0-0 -20-30 0 2 4 6 8 0 2 4 6 Similarities: Differences: Reason for differences: 7
F.4 (29 pts.) Dynamic response for an air-conditioner on a building An air conditioner is bolted to the roof of a one story building. The air conditioner's motor is unbalanced and its eccentricity is modeled as a particle of mass m attached to the distal end of a rigid rod of length r. When the motor spins with angular speed Ω, it causes the building's roof of mass M to vibrate. The stiffness and material damping in each column that supports the roof is modeled as a linear horizontal spring (k) and linear horizontal damper (b). Homework 4.5 showed that the ODE governing the horizontal displacement x of the building's roof is (M +m)ẍ + 2 b _x + 2 kx = mrω 2 sin(ωt) (a) (3 pts.) Express! n,, A in terms of M, m, b, k, r so the ODE is x k,b M Ωt r m k,b ẍ + 2! n _x +! 2 n x = AΩ2 sin(ωt)! n = = A = (b) (2 pts.) For certain values of M, m, b, k, andr, this ODE simplifies to Calculate numerical values for! n and. ẍ + 3 _x + 900 x = x0-4 Ω 2 sin(ωt)! n = rad sec = (c) (7 pts.) Fillinnumerical valuesinthefollowing expressions for x ss (t), the steady-state part of x(t). Ω ( rad sec ) x ss(t) 20 Λ sin( t + ) 30 Λ sin( t + ) 40 Λ sin( t + ) 8
(d) (4 pts.) The building's occupants complain that the roof shakes too much. Comment on Ω the effect small variations of M, m, b, k, and r have on:,! n, and the magnitude of x ss (t). Fill in each element in the table by writing (decreases), 0 (no effect), + (increases), or? (if it may decrease or increase). For the 2 nd -to-last column, assume the air conditioner's normal operating speed is Ω = 28 rad rad, and for last column, assume Ω = 32. sec sec Ω ß 28 rad sec Ω jx!n ss (t)j Balancing the motor (r! 0) 0 0 Increasing the motor speed Ω Decreasing the motor speed Ω Adding mass to the roof (increasing M) Removing mass from the roof Stiffening the support columns (increasing k) Adding damping to the columns (increasing b) Ω ß 32 rad sec jx ss (t)j (e) (3 pts.) List three ways to change the motor and minimize the roof shaking. 9