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?