Eccentrically hollow cylinder - Lagrangian (external forces)

[widelec]

Homework Statement

I am providing a solution up to the point when I'm having a little issue with defining the generalized force.

An eccentrically hollow cylinder of radius r rolls down a plane of inclination angle \alpha. Inside the cylinder, there is a cylinder-shaped hole of radius \frac{r}{2} and center located \frac{r}{2} from the cylinder's center. Using the Lagrange method, formulate the differential equations of the system, assuming the density of the cylinder is \rho and its length is equal to 1 \ m.

Assume the potential energy is zero and there is an external force acting on the system.

Variables and data:

l,r - length (= 1 m) and radius of the cylinder
d - distance between the cylinder's center of rotation and its center of mass (calculated below)
m=l \rho \pi r^{2} - mass of the full cylinder
m_{s}=l \rho \pi \left [ r^{2} - \left ( \frac{r}{2} \right )^{2} \right ] = \frac{3}{4}m - mass of the hollowed cylinder
g - gravitational acceleration
x_{1},y_{1} - coordinates of the center of rotation
x_{s},y_{s} - coordinates of the center of mass
\varphi_{1} - rotation angle of the cylinder

(h_s is not considered in this approach - no potential energy)

Homework Equations

The constraints are:

1. y_{1}=0
2. x_{1}=C_{1}-r \varphi_{1} (C1 being the initial position)
3. x_{s}=x_{1}+d \cos(\varphi_1), \ y_{s}=d \sin(\varphi_1)

\varphi_{1} is chosen as the generalized coordinate.

The Attempt at a Solution

The differentials are:

\dot{x}_{1}=-r \dot{\varphi}_{1}
\dot{x}_{s}=-r \dot{\varphi}_{1}-d \dot{\varphi}_{1} \sin(\varphi_{1})
\dot{y}_{s}=d \dot{\varphi}_{1} \cos(\varphi_{1})

The distance between centers of rotation and mass is calculated:

d=\frac{1}{m_{s}}\sum _{i=1}^{2} m_{i}d_{i}=\frac{4}{3m} \left (\frac{1}{2} m \cdot 0 +\frac{1}{4}m \cdot \frac{1}{2}r\right )=\frac{1}{6}r

Moment of inertia of the hollow cylinder with regard to the center of rotation (Steiner's theorem):

J_{0}=\frac{1}{2}mr^2 - \left [ \frac{1}{2} \left ( \frac{1}{4}m\right ) \left ( \frac{r}{2}\right )^2 + \left ( \frac{1}{4}m\right ) \left ( \frac{r}{2}\right )^2 \right ] = \frac{13}{32}mr^2

and in regard to the center of mass:

J_{s}=J_{o}-m_{s}d^2=\frac{37}{96}mr^2

Thus, the kinetic energy (and the Lagrangian, since potential is zero) of the system is equal to:

\mathcal{L}=T=\frac{1}{2}J_s \dot{\varphi}_1^2 + \frac{1}{2}m_s \left ( \dot{x}_s^2+\dot{y}_s^2\right )

Skipping the calculations, I've come to the following form:

\mathcal{L}=\frac{37}{64}mr^2\dot{\varphi}_1^2+\frac{1}{8}mr^2 \dot{\varphi}_1^2 \sin\varphi_1

So the respective differentials are:

\frac{\partial \mathcal{L}}{\partial \dot{\varphi}_1}=\frac{37}{32}mr^2\dot{\varphi}_1+\frac{1}{4}mr^2 \dot{\varphi}_1\sin\varphi_1
\frac{d}{dt} \left [ \frac{\partial \mathcal{L}}{\partial \dot{\varphi}_1} \right ]=\frac{37}{32}mr^2\ddot{\varphi}_1+\frac{1}{4}mr^2\left ( \ddot{\varphi}_1 \sin{\varphi_1} + \dot{\varphi}_1^2\cos{\varphi_1}\right )
\frac{\partial \mathcal{L}}{\partial \varphi}_1=\frac{1}{8}mr^2\dot{\varphi}_1^2\cos\varphi_1

Now, we need to formulate the Lagrange equation:

\frac{d}{dt} \left [ \frac{\partial \mathcal{L}}{\partial \dot{\varphi}_1} \right ]-\frac{\partial \mathcal{L}}{\partial \varphi}_1=\tilde{P}_{\varphi_1}

where \tilde{P}_{\varphi_1} is the generalized force acting in the direction of \varphi_1.

This is the point where I'm having trouble defining the generalized force.

Since we assume the cylinder does not slip along the plane, we should consider only the frictional force that causes the cylinder to roll.

Because frictional force acts in the same direction as the X axis, I write it out as:

F_{x_1}=-m_sg \sin{\alpha}=-\frac{3}{4}mg\sin{\alpha}

Now, to convert it to act along \varphi_1, I believe it needs to be done this way:

\tilde{P}_{\varphi_1}=\frac{\partial x_1}{\partial \varphi_1}F_{x_1}=(-r) \cdot (-\frac{3}{4}mg\sin{\alpha})=\frac{3}{4}mgr\sin{\alpha}

Is this going to be the only external force acting on the system? Does hollowing the cylinder eccentrically make any difference in terms of forces when compared to a "standard" full cylinder rolling down a plane?

Or perhaps I'm making a mistake somewhere? Any help would be greatly appreciated.

 

[paisiello2]

1) Assumption about no potential energy doesn't seem correct. Not only does the cylinder have potential energy but the potential energy changes as it rolls down the ramp.

2) However, maybe if you assume that gravity is an external force in this situation then you can ignore the potential energy.

3) When you introduced the constraint x_{1}=C_{1}-r \varphi_{1} then I think this is effectively your friction force. So no need to add this as an external force.

 

[widelec]

paisiello2, thank you very much for your reply.

However, I did come to a conclusion after an extensive discussion with my professor. I am now 100% sure this is correct.

We have to consider virtual work as follows (minus signs due to directions opposite to axes):

\delta A=(-m_{s}g\sin{\alpha}) \delta x_s+(-m_{s}g\cos{\alpha})\delta y_s

where:

\delta x_s=\frac{\partial x_s}{\partial \varphi_1}, \ \delta y_s=\frac{\partial y_s}{\partial \varphi_1}

thus:

\delta A=\frac{3}{4}mgr\left[ \sin{\alpha}-\frac{1}{6}\cos(\alpha+\varphi_1)\right ]=\tilde{P}_{\varphi_1}

Substituting the above to the aforementioned Lagrange equation along with the previously calculated derivatives gives a result consistent with the one given in my coursebook where potential energy is considered (and no external forces acting, in turn).