Homework and Reading guide



ME 111, Fall 2008


Homework guidelines: Homework is due at midnight Fri. in the box outside 122F Votey Hall. On every homework assignment please do the things listed below. After a few weeks of getting used to these you will have to follow all of them to get any credit on homework.

a) On the top right corner neatly print the following, making appropriate substitutions as appropriate:
    Sally Rogers
     HW 2 Due Oct. 5, 2007
     ME 111
b) STAPLE your homework in the top left corner.
c) Please: use only the front side of each sheet.
d) You are encouraged to work with others on homework; you will gain more by explaining to and learning from other students. At the top clearly acknowledge all help you got from TAs, faculty, students, or ANY other source (but for lecture, text, and office hours only). Examples could be "Mary Jones pointed out to me that I needed to draw the second FBD in problem 2." or "Nadia Chow showed me how to do problem 3 from start to finish." or "I basically copied this solution from a Tau Beta Sigma frat file." etc. If your instructor or grader thinks you are taking too much from other sources he/she will tell you. In the mean time, don't violate academic integrity rules by being being unclear about what of your presentation you worked out on your own.
e) Every use of force, moment, momentum, or angular momentum balance must be associated with a clear correct free body diagram . Every use of electrical laws must be accompanied by a clear correct circuit diagram, etc.
f) Your vector notation must be clear and correct.
g) Every line of every calculation must be dimensionally correct. (Carry your units.)
h) Your work should be laid out neatly enough to read. Part of your job as an engineer is not just to get the right answer, but to communicate its justification clearly. So, that is part of your job on the homework as well.
i) All computer code and output must be adequately documented and commented; your name should be clearly visible, as printed by the computer (e.g., title plots with your name, put your name in a comment in the first line of any *.m files, etc.) All computer plots must be titled and have axes labeled.
j) Write equations in symbolic form before plugging in numerical values and solving. It makes it easier for you and us to check your work. This is also a much desired practice for your exam solutions.


Homework grading and solutions: We will grade carefully two of the homework problems each week counting for 50% of the total points; the rest will be checked for an honest attempt to complete them and count for the remaining 50%. The carefully graded problems will not be announced in advance. We will penalize late homework and will not accept homework after the solutions are posted. Solutions will be posted at the course web site on the Monday after the homework is due.
Homework and reading assignments: (O) indicates Ogata reading assignments. Notes: (1) homework assignments are subject to change each week depending on how the lecture is progressing; and (2) the problems from Ogata are available below for each hw assignment by download because they are from the 3rd edition (whose problems I prefer). The text in the 4th ed. is equivalent.
Homework 1, Due Fri., Sept. 19:
Topics: Introduction to system dynamics; MATLAB review
Reading guide:
(O) Ch. 1 (Introduction to mathematical modeling, analysis and design of dynamical systems).
Problems:
(0) Please read the homework guidelines at the top of this WWW page. In particular, for this first assignment (and all subsequent MATLAB assignments), please pay attention to item i)
Please do and turn in one hw in groups of two. I have assigned groups randomly. Download the list here: MATLAB groups.

(1) Using MATLAB, write a script file to solve the system of 6 linear equations for the 6 unknown reactions for the plate of Problem 1 of ME 012, Prelim 1, 2004. Click here to download a copy of this exam.
Print out and turn in: (a) the script file you write and (b) the output results. Solve the 6 equations for 6 unkowns by hand to check the MATLAB solutions - just kidding!;P (Actually, not hard for this problem.)
(2) Derive the equation of motion for the simple pendulum with point mass m , length l, linear torsional damping at the pivot with constant c, and gravitational constant g (as we did in ME 012: FBD, linear or angular momentum balance, polar coordinates, etc.). Convert it to a system of first order equations. Then, solve the equations using MATLAB's ode45 using values you choose for m and l and over a time interval that you choose. Graph the angle and angular rate vs. time on separate plots. Print and turn in: (a) the derivative function file containing the pendulum equations, (b) the script file you write to solve them, (3) and the plots.
(3) Write a MATLAB function file to plot a `perfect' circle of any radius. Print and turn in the function file and a plot.
MATLAB file 1, MATLAB file 2, MATLAB file 3
Homework 2, Due Mon., Sept. 29.
Topic: Mechanical systems modeling review -- mechanical elements (mass, moment of inertia and parallel axis theorem, springs, dampers, friction); linear and angular momentum balance laws (force, torque); kinematics (straight-line and angular displacement, velocity, acceleration); energy methods (power balance, conservation of energy, work, kinetic and potential energy); equations of motion.
Reading guide: (O) Ch. 3
Problems: (1) B-3-4 (2) B-3-8(a) (3) B-3-10 (4) B-3-13 (5) B-3-21.
Problems for download
Solutions
Homework 3, Due Tues. Oct. 7.
Topic: Mechanical systems review continued.
Reading guide: (O) Ch. 3.
Problems: (1) B-3-14 (2) B-3-22. Hint: compare your answer to that obtained by energy methods in Ex. 3-4 on p. 79 and in class by power balance. (3) B-3-29. Modification: Change second sentence to read: " Find the force necessary to accelerate the mass upwards at 1 m/s^2." (4) B-3-30. Modification: Change the first sentence to read: " For the hoist shown in Figure 3-69, given that the plate AB and the block of mass m are accelerating at 1 m/s^2 due to the force F = 5N, find the mass m . (5) B-3-2 (harder). Hints and modification: For background, see (i) my notes on torsional pendula from class and (ii) the section "Experimental determination of moment of inertia" on p. 69 which is accomplished using a shaft of known stiffness. (a) In this problem, however, you are to find J without knowing the stiffness of the shaft. In the problem figure, you need to assume that the rotating body is suspended at its center of mass and that the mass is axisymmetric with respect to the rotation axis. (b) In addition, find the unknown stiffness of the shaft. Assuming that the mass is suspended from a thin circular steel shaft of given length and diameter, use your knowledge from mechanics of materials (see torsion chapter in Beer and Johnston) to propose a way to compare the stiffness obtained from the measurements to the theoretical shaft stiffness.
Problems for download
Solutions
Homework 4, Due Tues., Oct. 14.
Topic: Behavior of systems -- solving (1st and 2nd order) governing linear ordinary differential equations, using various methods; free (or natural) and forced (or driven) systems; step, ramp, and impulsive forcing; solving with MATLAB. How do these systems behave as functions of time? How do you change the system parameters to achieve desired behavior?
Reading guide: (O) Ch. 8. (corresponds to Ch. 6 in 3rd ed.)
Problems: (1) B-6-1 (simple 1st order electrical system). Solve by hand and numerically by using MATLAB with electrical parameters of your choice. Turn in: a script file that solves and plots the true AND numerical solutions for e_0(t); a derivative file; and one plot (e_0(t) vs. time). (2) B-6-7. Notes: there is gravity and 'x' is measured from the rest position of the mass. (3) B-6-8. Notes/suggestions/hints. There is gravity. Designate the masses as m_1 = 20 kg and m_2 = 2 kg. The position 'x' is measured from the rest position of m_1. I will give more suggestions and hints in class. (4) B-6-9 Notes/hints. We derived the equations of motion for this problem in class: two coupled second order nonlinear o.d.e.'s. (Please use the coordinate system choice we made in class!) Convert these to a system of four first order o.d.e.'s. Using MATLAB, solve the system of equations numerically for theta, thetadot, y, and ydot using the parameters and initial conditions given. Turn in: a script file that solves and plots the solution; a derivative file; and 2 plots (One with y and ydot vs. time and the other with theta and thetadot vs. time).
Problems for download
Solutions
Homework 5, Due Fri. Oct. 24.
Topic: Behavior of systems (continued) --
Reading guide: (O) Ch. 8.
Problems: (1) B-6-4. Ignore completely the problem instructions and follow mine here. The end of spring k_1 is given a constant displacement of x_i(t) = U at time t=0. (a) Find x_0(t). (b) What are x_0, y, x_i, and F (force in the system) as t -> infinity. Hints will be given in class. (2) B-6-12. Instead of delta(t), replace it simply with an impulsive force F(t). Follow the analogous problem solved in class. (3) B-6-17. Ignore completely the problem instructions and follow mine here. (a) Solve for x1(t) and x_2(t) by hand via the method discussed in class for the analogous problem. (b) Obtain their solutions using MATLAB. (c) Find the long-time solutions x1(inf), x2(inf), x1dot(inf), and x2dot(inf) from the hand solutions. Check to see that the t=0+ and long-time (t->inf) solutions agree with the plots of the numerical solution. Turn in: the hand solutions plus a script driver file, derivative file and plots of x1, x1dot, x2, and x2dot.

Hints : download the analogous example from class plus its MATLAB solutions here: Massless multiple spring-damper system example

Problems for download
Solutions
Homework 6, Due Wed., Nov. 19.
Topic: Behavior of mechanical and electrical systems in response to sinusoidal (or harmonic) forcing; multi-degree-of-freedom translational and rotational systems response; solving with MATLAB.
Reading guide: (O) Ch. 7.
Problems: [NOTE: THESE PROBLEMS ARE IN YOUR BOOK (4th edition)!]
(1)
B-9-4. Ignore gravity. Notes: steady-state solutions are the long-time or persisting solutions. (The transient (homogeneous) solutions become neglibible in some finite time and the steady-state (particular) solutions persist forever.) (2) B-9-14, Notes and additional assignments on B-9-14: (1) When the problem asks to "find the first and second modes of vibration", it means find the two natural frequencies and corresponding eigenvectors analytically as in lecture; (2) in addition to the hand solutions, find the eigenvalues and eigenvectors analytically with Mathematica (or MathCad or Maple) and numerically with MATLAB; and, (3) simulate the motions of the masses numerically and plot the mass positions vs. time using both (a) MATLAB and (b) Working Model 2D (WM2D). Export the position data from WM2D and import it to MATLAB and plot both on the same graph. For the numerical simulation, please use, m = 2 kg and k = 10 N/m. Use x1_dot(0)=x2_dot(0)=0 with two different initial conditions on the positions: (1) (x1(0)/x2(0))_1 corresponding to normal mode 1 and (2) (x1(0)/x2(0))_2 corresponding to normal mode 2. Use the numerical plots to estimate the two normal mode frequencies. Verify that there is agreement between the analytical solutions and numerical estimates of the frequencies.
Extra credit: redo the problem analytically and numerically with the middle spring having stiffness n*k instead of 2*k. What happens to the normal mode motions and frequencies when n is small and n is big? Do the numerical solutions agree with what behavior that the analytical solutions predict in each case? Explain what is happening physically.
Solutions
Homework 7, Due Wed., Dec. 3.
Topic: Equilibrium, linear and linearized stability, and control.
Reading guide: Class notes.
Problems: (1)Do the WM part from the last assignment where you are to compare the WM results to the MATLAB ones.

(2) My First Robotics Problem: Balance the inverted Pendulum.

In class we studied the stability of the two equilibriums of the simple pendulum (no friction at hinge or air drag, point mass m, massless rod of length ell). Using a LINEARIZED STABILITY ANALYSIS, we found that theta=0 (inverted or upright pos.) is an UNSTABLE equilibrium and theta=pi (hanging) is also an UNSTABLE (or, at best 'NEUTRALLY STABLE) equilibrium, in response to small disturbances or PERTURBATIONS.

In this problem, we replace the pendulum hinge with a motor that can supply a torque T_m = T_M(theta, theta_dot); that is, the motor torque is determined by an instantaneous measurement of the angle theta and its time rate of change, theta_dot. The aim is to design the motor to provide a torque that can stabilize the unstable inverted pendulum; that is, by using sensors to measure angle and angle rate, the goal for the automatic motor measurement and control system to stabilize the pendulum near theta=0 by providing appropriate torques to 'dampen' small disturbances to the perfect system (the system at rest, at theta=0). To picture what is going on, imagine that a small handle is rigidly attached to the pendulum at its base at right angles to its plane of motion. While holding the pendulum by the handle, imagine trying to restore it to the vertical as it tips, by twisting the handle in response to the tipping, where your response is proportional to the 'tip' and 'rate of tip' as you 'see' and 'feel' those quantities. (Of course, in the idealized model, we are ignoring the significant effect of the twisting friction between your hand and the handle.)

In particular, I have suggested a combination of LINEAR PROPORTIONAL and DERIVATIVE control torques, T_m= b*theta + c*theta_dot; many other functions of angle and angle rate will work, too.

Steps:

(A) To start, created a new FBD and re-derive the equations of motion using the same configuration from class, but this time with the motor torque acting on the end of the pendulum rod.

(B) Next, your job is to find constants b and c as functions of m, ell, and g so that the equilibrium q^* = {theta^*, theta_dot^*}={0,0} is stable; that is, so that perturbations to the system at this equilibrium die away exponentially. In mathematical terms, equilibriums are stable if the REAL PARTS OF ALL OF THE EIGENVALUES of the JACOBIAN (or 'derivative') MATRIX evaluated at q^*={0,0} ARE NEGATIVE. Use the linearized stability analysis procedure we employed in several examples in lecture.

(C) Then, once you have determined stabilizing values for b and c, test them on the FULL NONLINEAR equation of motion. That is, in the usual manner, do the followin:

(i) convert the new torque-driven equation of motion into a system of first order o.d.e's;

(ii) pick values for m and ell and find b and c in terms of those (and g); and,

(iii) integrate the system of equations using ode45.m and show that small initial disturbances to the system (i.e., small theta_0 and theta_dot_0), decay over time, thus demonstrating that the motor control torque can stabilize the pendulum ion the neighborhood of theta=0. Produce plots that show this behavior in the angle and angle rate.

(D) What is biggest disturbance that the systems is stable with respect to?