up

Spring mass system on a rotating disk

Nasser M. Abbasi

June 28, 2015
June 28, 2015 this page in PDF

1 Physical description of the problem

The following diagram shows the physical layout that illustrates the dynamics of a spring mass system on a rotating table or a disk. Mass (the bob) is attached to the end of a spring.

Mass vibrates moving back and forth at the end of a spring that is laid out along the radius of a spinning disk. The bob is considered a point mass. The spring is fixed to the center of the disk which is the origin of the inertial coordinates system. The bob is considered part of the overall disk mass for the purpose of calculation of the mass moment of inertia, therefore the overall mass moment of inertia around the axis of rotation will change as the bob moves along the radius.

There is no friction between the bob and the table. The spring mass is negligible. The spring is assumed to be infinitely rigid against any rotational movement. Hence the spring will move only in the radial direction (as observed in the inerial frame).

In other words, the spring will remain straight all the time no matter how fast the disk spins and how heavy the bob is.

This is important because it means the bob will have no speed component perpendicular to the radius which would normally be    β€²
xπœƒ and also it will have no Euler acceleration term perpendicular to the radius which would normally be xπœƒβ€²β€² . However, there will be Coriolis acceleration due to the point mass changing position along the radius as the disk spins.

This causes a Coriolis force which generates a torque on the disk. This torque however will change in value as the bob moves and the net effect is that the angular momentum of the overall system remain constant as will be shown below.

PIC

The equations of motion for the bob and the disk are derived using both Newton's and Lagrangian methods. It is a good practice to start the analysis by making a velocity, acceleration and force diagrams. The velocity and acceleration diagrams of the system are given below

PIC

The force diagram is

PIC

The force k(x βˆ’ x0)  is the spring force on the bob mass.

There will be a reaction force at the center of the disk equal to this force where the spring is attached. But this reaction force produces no torque on the disk since its arm length is zero.

2 Derivations using Newton's laws

For the bob, applying F  = ma  along the only degree of freedom the bob can move along, which is the radius of the disk gives

pict

Now       β€²β€²
Ο„ = Iπœƒ is applied to the disk (and remembering to include the bob mass as part of the disk) in order to obtain the equation of motion for the πœƒ  degree of freedom.

The torque is due to the Coriolis acceleration. The Coriolis force βˆ’ m (2xβ€²πœƒβ€²)  and the distance from the origin is x  , hence the toque is βˆ’ 2mxx  β€²πœƒβ€² (the minus sign due to the direction of the torque being clockwise). Therefore        β€²β€²
Ο„ =  Iπœƒ gives

            ( M  R2       )
βˆ’ 2mxx β€²πœƒβ€² =  ------+ mx2    πœƒβ€²β€²
                2

Hence

              β€² β€²
πœƒβ€²β€² = --βˆ’-4mxx--πœƒ---
     M  R2 + 2mx2
(2)

The above 2 equations are now solved numerically for x(t)  and πœƒ (t)

3 Derivations using Lagrangian method

The Lagrangian L  = T βˆ’ V  is first determined, where T  is the overall kinetic energy of the system and     V  is the overall potential energy of the system.

     1          2
V =  -k (x βˆ’ x0)
     2

and

                             (      )
     1   (  β€² 2     β€² 2)   1   M R2     β€²2
T =  -m   (x ) + (xπœƒ )   + --  ------ (πœƒ )
     2                     2     2

Hence

              1   (             )   1 ( M R2 )         1
L =  T βˆ’ V =  -m   (xβ€²)2 + (xπœƒ β€²)2 + --  ------ (πœƒβ€²)2 βˆ’ -k(x βˆ’ x0 )2
              2                     2     2            2

The x  motion gives

pict

Hence

pict

which is the same as equation (1) found using Newton's method. Now for the πœƒ  motion

pict

Hence

pict

which is the same as equation (2) found using Newton's method. Now for the πœƒ  motion.

4 Verification that angular momentum remain constant

It is important to verify that the angular momentum remain constant based on the above model. The angular momentum must be constant since there are no dissipative forces and no external forces. The initial angular momentum is given by       (               )
I πœƒβ€² = M  R2 + mx  (0 )2  πœƒβ€²(0)  where x (0)  is the initial position of the bob and  β€²
πœƒ (0)  is the initial angular velocity given to the disk. Let Ξ»  be the value of this initial angular momentum. Hence at any point of time in the future the following relation must hold

pict

Now, taking derivative with time of the above gives

pict

Which is the same equation of motion for πœƒ  found above using Newton's and Lagrangian methods. Hence the assumption that the angular momentum remain constant leads to the equation of motion. This shows that the angular momentum being constant does not contradicts the assumptions used to derive the equations of motion.

5 Solving and simulation of the system

The equations of motion are solved numerically. At the end of integration step, the new angular momentum is calculated and was verified that it remained equal to L (0 )

Using Mathematica, it is not required to convert the equations to first order, and equations (1) and (2) can be given to NDSolve as is. The following is an example showing how this was done. In this example, the following initial conditions and parameters are used

R = 2  meter,M  =  10  kg,x0 = R βˆ•2,k =  1  N/m,m =  5  kg,x (0) = 1.3  meter,xβ€²(0) = 0,πœƒ (0) = Ο€βˆ•4,πœƒβ€²(0) = 0.007  rad/sec.

The equations are solved and x (t)  and  β€²
πœƒ (t)  are plotted. This was done for t = 120  seconds


PIC


The plots are

PIC

and

PIC

To verify that the angular momentum is conserved, it is calculated and plotted as well. The result shows it is constant. This verifies that the numerical solution is correct.

PIC