Lagrangian - rigid body problem

[Glomerular]

Homework Statement

In a uniform gravitational field, there is a uniform solid disk of of mass M and radius R. A point mass m is glued to the disk at a point that is at a distance a from the center of the disk.
The disk rolls without slipping. Find the frequency of small oscillations about the equilibrium point.

I have been told to solve the problem in two ways:
1. Considering two independent bodies, the disk and the point mass.
2. Considering only one body by using their center of mass.

Homework Equations

Euler-Lagrange equations: \frac{\partial L}{\partial q} - \frac{d}{dt}\frac{\partial L}{\partial {\dot q}}=0 (1)
Parallel axis theorem : I=I_{cm}+ MR^2 (For the second way) (2)

The Attempt at a Solution

First way:

The position of the center of mass of the disk is:

\begin{cases}<br /> x_{cm}=R\varphi \\<br /> y_{cm}=R<br /> \end{cases}
And the position of the point mass is:

\begin{cases}<br /> x_{p}=x_{cm} - a\sin(\varphi)=R\varphi - a\sin(\varphi)\\<br /> y_{p}=y_{cm} - a\cos(\varphi) =R - a\cos(\varphi)<br /> \end{cases}
So I can compute the kinetic energy K and the potential energy U as:

K=K_{cm}+K_{p} and U=U_{cm}+U_{p}

where:

K_{cm}=\frac{1}{2} M{v^2}_{cm} +\frac{1}{2}I_{disk}\dot{\varphi}^2 = \frac{3}{4}MR^2\dot{\varphi}^2

(I have used I_{disk}=\frac{1}{2}MR^2)

K_{p}=\frac{1}{2} m({\dot{x}^2}_{p}+{\dot{y}^2}_{p} )+\frac{1}{2}I_{p}\dot{\varphi}^2 = \frac{1}{2}mR^2 + \frac{1}{2}ma^2 -maR\cos(\varphi) + \frac{1}{2}ma^2\dot{\varphi}^2

(I have used I_{p}=ma^2)

U_{disk}=MgR

U_{p}=mg(R-a\cos(\varphi))

Thus, the lagrangian is:

L=\frac{1}{2}mR^2 + \frac{1}{2}ma^2 -maR\cos(\varphi) + \frac{1}{2}ma^2\dot{\varphi}^2 +\frac{3}{4}MR^2\varphi -MgR - mg(R-a\cos(\varphi))

And using the Euler-Lagrange equation (1) and approximating sin(\varphi)=\varphi I get

\ddot{\varphi}+{\varphi}[\frac{mRa + mga}{\frac{3}{2}MR^2 +ma^2}]=0

Where \omega^2 is the term next to \varphi which is not the result I should get.

Where is my error? I have spent hours looking and re-looking at this problem.

Second way:

My problem is that I don't know how to get the rolling condition for this system. The rest should be practically the same as I have done before, but computing the inertia moment by using the parallel axis theorem (2).

Thank you very much, I really appreciate some help on this.

 

[OldEngr63]

Kcm=12Mv2cm+12Idiskφ˙2=34MR2φK_{cm}=\frac{1}{2} M{v^2}_{cm} +\frac{1}{2}I_{disk}\dot{\varphi}^2 = \frac{3}{4}MR^2\varphi

Looks like you left off a dot and an exponent 2 in this expression .. just typos I presume.

(I have used Ip=ma2I_{p}=ma^2)

This is the MMOI for the point mass with respect to the center of the disk, but there is no reason to include this in the kinetic energy of the point mass. The point mass, by definition of the term point mass, has no MMOI about its own center of mass, and that would be the only term to include in addition to the point mass velocity^2 term.
The rolling constraint is easy, and you already know it. It is simply
s = R * phi
This may be differentiated as needed.

 

[Glomerular]

First of all thanks a lot for helping me with this problem.

Looks like you left off a dot and an exponent 2 in this expression .. just typos I presume.

I have already corrected that typo :)

This is the MMOI for the point mass with respect to the center of the disk, but there is no reason to include this in the kinetic energy of the point mass. The point mass, by definition of the term point mass, has no MMOI about its own center of mass, and that would be the only term to include in addition to the point mass velocity^2 term

You mean K_{p}=\frac{1}{2} m({\dot{x}^2}_{p}+{\dot{y}^2}_{p} ) with no additional rotating-motion terms right?
I feel so foolish about considering the distance with respect to the center of the disk

The rolling constraint is easy, and you already know it. It is simply
s = R * phi
This may be differentiated as needed.

I thought this would change when considering only 1 body. Then, is there any consideration I should take appart from the MMOI? It looks too easy.
Anyway, I'm going to solve it and edit #1

 

[OldEngr63]

The disk rolls without slipping.

As long as this is true, then one and the same constraint applies all the time.

 

[Glomerular]

As long as this is true, then one and the same constraint applies all the time.

You're right, I have thought about it and seems totally logic.

 

[Glomerular]

Now, the second part of the problem is solving it by using the center of mass of the system disk+point mass.
I have a doubt when computing the Inertia Moment of my new system.

The system's center of mass position is:

\begin{cases}<br /> x_{cm}=x_{0} - d\sin{\varphi}=R\varphi - d\sin{\varphi}\\<br /> y_{cm}=y_{0} - d\cos{\varphi} =R - d\cos{\varphi}<br /> \end{cases}

Where d is the distance from the center of the disk to the center of mass:

d=\frac{ma}{m+M}

Then, It is possible to write K as:

K=K_{trans}+K_{rot}

Where:

K_{trans}=\frac{1}{2} (M+m)({\dot{x}^2}_{cm}+{\dot{y}^2}_{cm})

And:

K_{rot}=\frac{1}{2}I_{cm}\dot{\varphi}^2

Here is my doubt. I get stuck when computing I_{cm}.

Should I use Steiner's theorem on each object and then sum them?
\begin {cases}<br /> I^{disk}_{cm}=I^{disk}_0 -Md^2 \\<br /> I^{p}_{cm}=I^{p}_0 -m(d-a)^2<br /> \end {cases}

And then just do I_{cm}=I^{disk}_{cm} + I^{p}_{cm}

I don't see it clear.

 

[OldEngr63]

The MMOI for the combined body involves two applications of the parallel axis theorem:

Jcombined = Jdisk + M*d^2 + Jpoint + m*(a-d)^2

Then take the equations you have for xcm and ycm for the combined system, and differentiate them with respect to time to get the velocity components to write the kinetic energy. Don't make this a hard problem; it is not.

 

[OldEngr63]

Just as a matter of curiosity, where are you from?