Transcribed Text
The MassSpring System
In this laboratory we will examine harmonic oscillation. We will model the motion of a massspring
system with differential equations.
Our objectives are as follows:
1. Determine the effect of parameters on the solutions of differential equations.
2. Determine the behavior of the massspring system from the graph of the solution.
3. Determine the effect of the parameters on the behavior of the massspring.
The primary MATLAB command used is the ode45 function.
MassSpring System without Damping
The motion of a mass suspended to a vertical spring can be described as follows. When the spring
is
not loaded it has length lo (situation (a)). When a mass m is attached to its lower end it has length l
(situation (b)). From the first principle of mechanics we then obtain
mg

(L5.1)
downward weight force upward tension force
The term g measures the gravitational acceleration (g ~ 9.8m/s² ~ 32ft/s2). The quantity k: is a spring
constant measuring its stiffness. We now pull downwards on the mass by an amount y and let the mass
go (situation (c)). We expect the mass to oscillate around the position y = 0. The second principle of
mechanics yields
mg + thiltylo)  = m dt² +y) , i.e., m + ky = 0
d2 (le
dy
(L5.2)
dt²
weight upward tension force
acceleration of mass
acceleration of mass
using (L5.1). . This ODE is secondorder.
lo
k:
V
l
V
Y
m
V
r
(a)
(b)
(c)
(d)
Equation (L5.2) is rewritten
dy
dt² + way = 0
where w? = k/m. Equation (L5.3) models simple harmonic motion. A numerical solution with ini
tial conditions y(0) = 0.1 meter and y' (()) = 0 (i.e., the mass is initially stretched downward
10cms
and
released, see setting (c) in figure) is obtained by first reducing the ODE to firstorder ODEs (see
Laboratory 4).
Let U = y'. Then v' = y" =  way = 4y. Also v(0) = y'(0) = 0. The following MATLAB program
implements the problem (with wo = 2).
function LAB05ex1
m = 1;
% mass [kg]
k = 4;
% spring constant [N/m]
omega0 = sqrt (k/m) ;
y0 = 0.1; V0 = 0;
% initial conditions
[t, Y] = ode45(Of, [0, 10], [y0,v0],[], omega0) ; % solve for 0y = Y(: , 1) ; V = Y(:,2);
% retrieve y, V from Y
figure (1) ; plot (t,y, 'b+', t,v, 'ro') ;
% time series for y and V
grid on;
%
function dYdt = f (t, Y, omega0)
y = = ;
dYdt = [ V ; omega0^2*y ];
Note that the parameter wo was passed as an argument to ode45 rather than set to its value wo = 2
directly in the function f. The advantage is that its value can easily be changed in the driver part of the
program rather than in the function, for example when multiple plots with different values of wo need
to be compared in a single MATLAB figure window.
0.2
Note that the parameter wo was passed as an argument to ode45 rather than set to its value wo = 2
directly in the function f. The advantage is that its value can easily be changed in the driver part
of
the
program rather than in the function, for example when multiple plots with different values of wo need
to be compared in a single MATLAB figure window.
0.2
0.15
0.1
0.05
0.05
O
0.1
O
O
O
0.15
O
O
O
0.2
1
2
3
4
5
6
7
8
9
10
Figure L5a: Harmonic motion
1. From the graph in Fig. L5a answer the following questions.
(a) Which curve represents y = y(t) ? How do you know?
1
1
11
/1
1
1. From the graph in Fig. L5a answer the following questions.
(a) Which curve represents y = y(t)? How do you know?
(b) What is the period of the motion? Answer this question first graphically (by reading the
period from the graph) and then analytically (by finding the period using wo).
(c) We say that the mass comes to rest if, after a certain time, the position of the mass remains
within an arbitrary small distance from the equilibrium position. Will the mass ever come to
rest? Why?
c 2011 Stefania Tracogna, SoMSS, ASU
MATLAB sessions: Laboratory 5
(d) What is the amplitude of the oscillations for y?
(e) What is the maximum velocity (in magnitude) attained by the mass, and when is it attained?
Make sure you give all the tvalues at which the velocity is maximum and the corresponding
maximum value. The tvalues can be determined by magnifying the MATLAB figure using
(d) What is the amplitude of the oscillations for y?
(e) What is the maximum velocity (in magnitude) attained by the mass, and when is it attained?
Make sure you give all the tvalues at which the velocity is maximum and the corresponding
maximum value. The tvalues can be determined by magnifying the MATLAB figure using
the magnify button
2
, and by using the periodicity of the velocity function.
(f) How does the size of the mass m and the stiffness k: of the spring affect the motion?
Support your answer first with a theoretical analysis on how wo  and therefore the period
of the oscillation  is related to m and k, and then graphically by running LAB05ex1. m first
with m = 5 and k = 4 and then with m = 1 and k: = 16. Include the corresponding graphs.
2.
The energy of the massspring system is given by the sum of the potential energy and kinetic
energy. In absence of damping, the energy is conserved.
(a) Plot the quantity E = mv²+ky² as a function of time. What do you observe? (pay close
attention to the yaxis scale and, if necessary, use ylim to get a better graph). Does the graph
confirm the fact that the energy is conserved?
(b)
Show analytically that dE = 0. (Note that this proves that the energy is constant).
(c) Plot U VS y (phase plot). Does the curve ever get close to the origin? Why or why not? What
does that mean for the massspring system?
MassSpring System with Damping
When the movement of the mass is damped due to viscous effects (e.g., the mass moves in a cylinder
containing oil, situation (d)), an additional term proportional to the velocity must be added.
The
resulting equation becomes
MassSpring System with Damping
When
the movement of the mass is damped due to viscous effects (e.g., the mass moves in a cylinder
containing oil, situation (d)), an additional term proportional to the velocity must be added. The
resulting equation becomes
m dt² + C dt + = 0 or dt² + = 0
(L5.4)
by setting p = C The program LAB05ex1 is updated by modifying the function f:
function LAB05ex1a
m = 1;
% mass [kg]
k C = 4; 1;
% spring constant [N/m]
=
% friction coefficient [Ns/m]
omega0 = sqrt (k/m) ; p = c/ (2*m);
y0 = 0.1; v0 = 0;
% initial conditions
[t,Y] = ode45(Of, , [0, 10], [y0,v0], [], omega0,p) ; % solve for 0y = Y( : , 1) ; V = Y(:,2);
% retrieve y, V from Y
figure (1) ; plot (t, y, 'ro');
% time series for y and V
grid on;
%
function dYdt = f (t,Y,omega0,p)
y = Y (1) ; V = Y (2) ;
dYdt = [ V ; ?? ] ;
% fillin dv/dt
3. Fill in LABO5ex1a.: m to reproduce Fig. L5b and then answer the following questions.
(a) For what minimal time t1 will the massspring system satisfy y(t) < 0.01 for all t > t1? You
can answer the question either by magnifying the MATLAB figure using the magnify button
3. Fill in LAB05ex1a. m to reproduce Fig. L5b and then answer the following questions.
(a) For what minimal time t1 will the massspring system satisfy y (t) < 0.01 for all t > t1? You
can answer the question either by magnifying the MATLAB figure using the magnify button
2
(include a graph that confirms your answer), or use the following MATLAB commands
(explain):
c
2011 Stefania Tracogna, SoMSS, ASU
MATLAB sessions: Laboratory 5
0.2
y(t)
v(t)=y'(t)
0.15
0.1
IV Laboravory
0.2
y(t)
v(t)=y'(t)
0.15
0.1
0.05
0.05
0.1
0.15
0.2
1
2
3
4
5
6
7
8
9
10
Figure L5b: Damped harmonic motion
for i=1 : length(y)
m(i)=max(abs(y(i:end))) ;
end
i = find(m<0.01) ; i = i(1) ;
disp (['ly/<0.01 for t>t1 with , num2str (t(i1)) , (b) What is the maximum (in magnitude) velocity attained by the mass, and when is it attained?
(b) What is the maximum (in magnitude) velocity attained by the mass, and when is it attained?
Answer by using the magnify button and include the corresponding picture.
(c) How does the size of c affect the motion? To support your answer, run the file LABO5ex1. m
for C = 2, C = 4, C = 6 and c = 8. Include the corresponding graphs with a title indicating
the value of c used.
(d) Determine analytically the smallest (critical) value of C such that no oscillation appears in
the solution.
4.
(a)
Plot the quantity E = as a function of time. What do you observe? Is the
energy conserved in this case?
(b) Show analytically that
dE
< 0 for C > 0 while dE > 0 for C < 0.
dt
dt
(c) Plot U VS y (phase plot). Comment on the behavior of the curve in the context of the motion
of the spring. Does the graph ever get close to the origin? Why or why not?
These solutions may offer stepbystep problemsolving explanations or good writing examples that include modern styles of formatting and construction
of bibliographies out of text citations and references. Students may use these solutions for personal skillbuilding and practice.
Unethical use is strictly forbidden.
function LAB05ex1
%Original values
m = 1; % mass [kg]
k = 4; % spring constant [N/m]
%Values for part 1.f) comparison for different values of k and m
% m = 5; % mass [kg]
% k = 4; % spring constant [N/m]
% m = 1; % mass [kg]
% k = 16; % spring constant [N/m]
omega0=sqrt(k/m);
y0=0.1; v0=0; % initial conditions
[t,Y]=ode45(@f,[0,10],[y0,v0],[],omega0); % solve for 0<t<10
y=Y(:,1); v=Y...