Problems on Simple Harmonic Motion

1.

        Note:  I showed this one in class! Here it is again in Maple.

>    equ1:=c1*sin(omega*t)+c2*cos(omega*t);

equ1 := c1*sin(omega*t)+c2*cos(omega*t)

>    first_derivative:=diff(equ1,t);second_derivative:=diff(first_derivative,t);

first_derivative := c1*cos(omega*t)*omega-c2*sin(omega*t)*omega

second_derivative := -c1*sin(omega*t)*omega^2-c2*cos(omega*t)*omega^2

Notice:   second_derivative = -omega^2 (equ1) and so it follows that c1*sin(omega*t)+c2*cos(omega*t)  satisfies the given differential equation.

2.

          Note:  This I also did in class. Here it is again in Maple.

>    equ2:=M*cos(omega*t+delta);

equ2 := M*cos(omega*t+delta)

>    expand(equ2,t);

M*cos(omega*t)*cos(delta)-M*sin(omega*t)*sin(delta)

Just let c[1] = -M*sin(delta)  and c[2] = M*cos(delta)  and the above becomes equ1. (Remember, delta  is a constant so that sin(delta)  and cos(delta)  are both constants.)

3.

        This one is a straight forward calculation and algebraic re-arrangement. Do it by hand!! Here it is in Maple Syrup:

>    position:=M*cos(omega*t+delta);velocity:=diff(position,t);acceleration:=diff(velocity,t);

position := M*cos(omega*t+delta)

velocity := -M*sin(omega*t+delta)*omega

acceleration := -M*cos(omega*t+delta)*omega^2

>    acceleration/position;

-omega^2

Notice:   acceleration = -omega^2*position  so acceleration is directly proportional to displacement.

4.

     This was done in class (a long  time ago -- before calling it the "integrative activity") using calculations. Another way to think about it, without  calculations, is as follows. We know the velocity and acceleration curves are ordinary sine and cosine curves respectively. Therefore, we know that these have their maximum and minimum values only where they have horizontal tangent lines  (this is not true in general -- can you think of reasons why?). Since speed is the absolute value of velocity, we know the speed is at a maximum when the velocity curve is at a max or a min. These are precisely where the tangent lines to its graph are horizontal. The tangent lines are horizontal exactly where the derivative of the velocity function is zero and the derivative of the velocity function is the acceleration function! That is, the speed is at a max when the acceleration is zero. Conversely, when the acceleration is zero, the velocity curve has a horizontal tangent line and this only happens when the velocity is at a maximum or minimum value (since the velocity curve is a sine curve).

5.

        a)

>    position:=1/5*cos(8*t);

position := 1/5*cos(8*t)

>    velocity:=diff(position,t);acceleration:=diff(velocity,t);

velocity := -8/5*sin(8*t)

acceleration := -64/5*cos(8*t)

>    acceleration/position;

-64

Note:  aceleration = -64(position) so we have that x'' = -k*x/m  where -k/m  = -64  ( so that k/m = omega^2 ) .  i.e. the position equation obeys the differential equation  for Hooke's Law.

      b)         Displacement s  is at a maximum from equilibrium when the velocity is zero: i.e. when Diff(s,t) = 0  and   -8/5*sin(8*t) = 0  when 8*t = 0+k*Pi . This means the first two time are t[1] = 0  and t[2] = Pi/8  seconds. The plot below of position vs. time confirms this.

>    plot(position,t=0..2,title=`position vs time`,thickness=2);

[Maple Plot]

               The Speed (recall: s = abs(v)  ) is at a maximum when the acceleration is zero (which is where the tangent line to the velocity graph is horizontal). This happens when Diff(v,t) = 0   which has the solution 8*t = Pi/2+k*Pi  or t = Pi/16+k*Pi/8  so that the first two times are t[1] = Pi/16  and t[2] = 3*Pi/16  and these are confirmeds in the plot below of velocity vs. time.

>    plot(velocity,t=0..2,title=`Velocity vs time`,thickness=2,colour=blue);

[Maple Plot]

6.

>    s:=sin(Pi/2*t)+cos(Pi/2*t);
v:=diff(s,t);
a:=diff(v,t);

s := sin(1/2*Pi*t)+cos(1/2*Pi*t)

v := 1/2*cos(1/2*Pi*t)*Pi-1/2*sin(1/2*Pi*t)*Pi

a := -1/4*sin(1/2*Pi*t)*Pi^2-1/4*cos(1/2*Pi*t)*Pi^2

a)

Note:  The particle changes direction when its velocity  is zero.

>    solve(v=0,t);

1/2

To get the REAL  solution you must change 1/2*cos(1/2*Pi*t)*Pi-1/2*sin(1/2*Pi*t)*Pi = 0  into the form tan(Pi/2*t) = 1  which has the general solution Pi*t/2 = arctan(1)+k*Pi    or   t = 2/Pi*(Pi/4+k*Pi)  and so the first two solutions are t[1] = 1/2  and t[2] = 5/2  and these are confrimed in the plot below.

>    plot(v,t=0..4,thickness=2,title=`velocity vs time`);

[Maple Plot]

b)

To show the particle is in simple harmonic motion it is enough to show that the equation Diff(`s(t)`,t,t) = -omega^2*s(t)  is satisfied. The unfamiliar notation is just Maple's way of saying "the acceleration (second derivative of position) is directly  proportional to the position." An easy way to show this is to divide the acceleration equation by the position equation and verify that the ratio is a constant. See below.

>    simplify(a/s);

-1/4*Pi^2

This confirms that, indeed, a(t) = -1*Pi^2/4 s(t)  and so the motion is simple harmonic . Therefore, -omega^2 = -Pi^2/4  so omega = Pi/2  (which we already knew)

c)

The amplitude  of the motion is abs(M)  where M  is obtained from the equivalent equation of the motion s(t) = M*cos(omega*t+delta) . To get this a different way  from what we did in class (or in the notes above) all you have to realize is that amplitude of s(t) corresponds to the maximum  displacement of the particle. This occurs any time the velocity is zero. We already know that v(t) = 0  when t = 1/2  so we simply "plug" this value into the postion function to find the amplitude:

>    subs(t=1/2,s);amplitude=simplify(%);

sin(1/4*Pi)+cos(1/4*Pi)

amplitude = 2^(1/2)

Therefore, it follows that there are two different ways to find the amplitude. The way we did it in class (and in the notes above) is way cooler .