[Math] Pendulum with moving support

geometryphysicspolar coordinates

Consider a mass $m$ connected to an oscillating pivot by a massless rod of fixed length $l$. The pivot is located at $(x_0(t), 0)$, where $x_0(t)=A_0\cos(\omega_0t)$. The $y$-axis is positive downwards. The angle between $y$ and the rod is $\theta(t)$. $\theta$ is the only generalized coordinate.

I have a hard time deriving the kinetic energy, which should be $T=\frac{1}{2}m\left(l^2\dot{\theta}^2+2l\dot{\theta}\dot{x_0}\cos\theta+\dot{x}_0^2\right)$.

Here's my thinking:

$$v_{rod}=l\dot{\theta}$$
$$v_{pivot}=\dot{r}_{pivot}$$

$r_{pivot}=x_0\implies\dot{r}_{pivot}=\dot{x}_0$
So $$v_{total}^2=\left(l\dot{\theta}+\dot{x}_0\right)^2$$

But something doesn't add up. I think that something is wrong with the velocity equation for the pivot, because if $\vec{r}_{pivot}=r\hat{r}$ then it somehow depends on $\theta$, but the angle of the moving pivot with the $y$-axis is always fixed. I'm quite confused. Would appreciate some help.

Best Answer

For me, the first step is to make a good sketch, then write the position vector and check it against the sketch. Carefully checking it for various angles is so important. When $\vec{r}(t)$ is correct, then take the derivative.
$$\vec{r} = \vec{r}_{pivot} + \vec{r}_{rod} \\ \vec{r}(t) = x_0(t) \hat{i} + (L\sin\theta(t) \hat{i} + L\cos\theta(t) \hat{j}) \\ \vec{v}(t) = [\dot{x}_0(t) + L\dot{\theta}\cos\theta(t)] \hat{i} - L\dot{\theta}\sin\theta(t) \hat{j} \\ \vec{v}\cdot\vec{v} = (\dot{x}_0 + L\dot{\theta}\cos\theta)^2 + L^2\dot{\theta}^2\sin^2\theta \\ v^2 = \dot{x}_0^2 + 2L\dot{x}_0\dot{\theta}\cos\theta + L^2\dot{\theta}^2$$

Having a good expression for the position vector also comes in handy for writing the potential energy function, $PE=-\vec{F}_{gravity}\cdot\vec{r}$.

Clarifications:

The position of the mass $m$ is given by $\vec{r}$, which always points from the origin to the mass at the end of the pendulum rod. There is also a vector $\vec{r}_{pivot}$ that always points from the origin to the pivot point of the pendulum. The position of the pivot is given as always on the $x$-axis at a distance $x_0(t)$ from the origin, so $\vec{r}_{pivot}=x_0(t)\hat{i}$

The vector $\vec{r}_{rod}$ always points along the pendulum rod from the pivot point to the mass no matter what angle, $\theta$, the rod makes with the vertical. The components of $\vec{r}_{rod}$ are found by projecting $\vec{r}_{rod}$ onto the horizontal and vertical axes. The length of $\vec{r}_{rod}$ is always $L$. By our choice of $\theta$, the vertical projection is $L\cos\theta$ and the horizontal projection is $L\sin\theta$.

$\theta$ is the angle that the rod makes with the vertical line through the pivot point. This is an important point for two reasons. First, $\theta$ is often used as the polar coordinate about the origin, but here it is a generalized coordinate that indicates the rotation of the pendulum rod about its pivot point. Second, we have chosen $\theta$ so that it is zero when the pendulum is in equilibrium. We could have used a different angle, such as the angle $\phi$ between the rod and horizontal, say, but then we would have a different expression for $v^2$.

Thus, the vector from the origin to the mass is the sum of the vector from the origin to the pivot and the vector from the pivot to the mass.

As for the formula for the potential energy, let $\vec{F}$ be the gravitational force on mass $m$ at position $\vec{r}$ in a uniform gravitational field. We can calculate the force two ways. One is simply $\vec{F}=mg\hat{j}$, with the unit vector $\hat{j}$ pointing down along the positive $y$-axis. The other is $\vec{F}=-\vec{\nabla}{V(\vec{r})}$, where $V(\vec{r})$ is the potential energy function. Since the gravitational field is uniform, $\vec{F}$ has no $\vec{r}$ dependence. It is not difficult to let $V(\vec{r})=-\vec{F}\cdot\vec{r}$ and verify that $$-\vec{\nabla}{V(\vec{r})} = -\vec{\nabla}{(-\vec{F} \cdot \vec{r})}=(\vec{F} \cdot \vec{\nabla})\vec{r} = \vec{F} $$ Therefore $V(\vec{r}) = -\vec{F} \cdot \vec{r}$ satisfies our one requirement for it to be the potential energy function. When we have a complicated expression for $\vec{r}$ and we need the potential energy, this is a very handy formula. Two warnings, however: First, our force $\vec{F}$ must be the same everywhere, no $\vec{r}$ dependence, which is often the case, for this formula to be valid. Second, the potential energy is not a unique quantity, so we might not get the exact same function as someone else, the difference usually being an additive constant. We are free to add or subtract any constant we want to our potential energy.

Related Question