Lagrangian: singularity in inverted pendulum EoMs?

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
1 reply · 2K views
Eigenentity
Messages
7
Reaction score
0
Hi all, I'm doing some analysis of a bicycle mechanics problem and at one point the approximations I'm making mean that the problem reduces to the classic inverted pendulum. I'm very confused, as the equations I've worked out appear to have an unphysical singularity in them, and I can't see where it's coming from.

So, the problem is this: equations of motion of a point mass [tex]m[/tex] on the end of a light rigid rod (length [tex]l[/tex]), the other end of which is attached to a light, frictionless carriage, constrained to move only in the horizontal direction.

So, the problem has two degrees of freedom: [tex]x[/tex] and [tex]\theta[/tex].

By my reckoning, the position vector of the pendulum bob, counting angle clockwise from the upwards vertical, is:

[tex]\mathbf{r}=(x + l \sin\theta, l \cos\theta)[/tex]

This gives a Lagrangian of:

[tex]L = \frac{m \dot{\mathbf{r}}^2}{2} - m g \cos{\theta} = \frac{m \dot{x}^2}{2} + \frac{m l^2 \dot{\theta}^2}{2} + m \dot{x} \dot{\theta} l \cos{\theta} - m g \cos{\theta}[/tex]

Now, the Euler-Lagrange equations end up as

[tex]\ddot{x} + \ddot{\theta} l \cos{\theta} - \dot{\theta}^2 l \sin{\theta} = 0[/tex]
[tex]\ddot{\theta} + \ddot{x} \cos\theta / l - g \sin\theta/l=0[/tex]

which in turn reduce to:

[tex]\ddot{x} = \frac{\dot{\theta}^2 l - g \cos \theta }{\sin\theta}[/tex]
[tex]\ddot{\theta} = \frac{g / l - \dot{\theta}^2 \cos\theta}{\sin\theta}[/tex]

Do you spot the problem yet? There's a nasty singularity when [tex]\theta \to 0[/tex].

Can anyone tell me what's going wrong here?
 
Physics news on Phys.org
To anyone wondering, the "unphysical" singularity arises from the unphysical assumption of a light carriage. With a carriage of nonzero mass [tex]M[/tex] the EoMs become:

[tex] \ddot{x} = \frac{2 m \sin \theta \left(g \cos \theta - l \dot{\theta}^2\right)}{2 M + m (1 - \cos 2 \theta)}[/tex]
[tex] \ddot{\theta} = \frac{m l \dot{\theta}^2 \cos \theta-g (m+M)}{- m l \sin \theta - M l / \sin\theta}[/tex]

With no singularity as long as [tex]M \neq 0[/tex]!