# Motion equations for Woodpecker toy (multibody system)

I am trying to understand a rigid multibody model of a Woodpecker toy (see figure below). Now I am not going to go into details about the model or justify this approach, I am just trying to understand the equations of motion for the woodpecker during different stages of its path down the road. More specifically I am trying to understand how they are derived: Is it through the Lagrangian? The Hamiltonian? Or Newtonian mechanics? (The last one is most likely).

For the record, the Woodpecker toy consists of a rod (static) on which a sleeve (ring) is connected. The sleeve is then connected to the Woodpecker body joined by a spring. The equation(s) of motion for the two stages of the Woodpecker toy are (using this model):

Eq. I : ($$\theta$$ is the only degree of freedom (DOF))

$$(I_2+m_2b^2)\ddot{\theta}=-c\theta+m_2bg_{gr}\tag{1}$$

Eq. II : (two DOF, vertical motion $$z$$ included)

$$(I_2+m_2b^2(1-\frac{m_2}{m_1+m_2}))\ddot{\theta}=-c\theta\tag{2}$$ $$(m_1+m_2)\ddot{z}+m_2b\ddot{\theta}=(m_1+m_2)g_{gr}\tag{3}$$

Now, the first equation corresponds to the Woodpecker being jammed (not moving vertically), thus only DOF is $$\theta$$. The equation is easily derived from the law of angular momentum: $$\frac{dL}{dt}=I\boldsymbol\alpha+2rp_{||}\boldsymbol\omega.$$ Since it is planar motion it reduces to one equation (the equation above).

The second equation corresponds to the Woodpecker (free) falling. It moves in vertical direction as well as rotates. We have two (coupled) equations for this motion, but how are they derived exactly? I recognize the first equation $$(2)$$ as the Law of angular momentum, but inside it there is this new factor/coefficient $$(1-\frac{m_2}{m_1+m_2})$$, that I believe comes from the motion of the sleeve. But why this term exactly? I cannot see the connection. Equation $$(3)$$ is Newton's second law, but why is there an angular acceleration component in there? How was this equation derived?

• Note that all three formulations of mechanics will lead to the same equations of motion, so your first series of questions is a bit confusing. Feb 7, 2019 at 16:24
• I think it's reasonable in absolute terms (though perhaps not by the Physics SE homework guidelines) to ask something like "how they are derived? Newtonian? Lagrangian? Hamiltonian?" because though you should end up with equivalent solutions from all methods, it's frequently found that one of them gives you a particularly easy procedure. However, I'm not aware if there's any concrete rule for knowing which method is the best beforehand, since the difficulty of methods is subjective in the first place.
– user191954
Feb 7, 2019 at 16:27
• Thanks for the replies! Kyle Kanos and Chair: Firstly I am intrested in how to derive them, that was the question, because now I really don't have an idea why the equations hold. Secondly I am intrested in what formalism, but as I understood it (from my intense research on the subject this week) the number of equations differ between Lagrange equations and Hamiltons equations (the order of the ODE:s are different). Also they have different variables: Lagrange only generalized coordinates/velocities and the Hamiltonian a generalized momentum variable. I am aware of transformation between them. Feb 7, 2019 at 18:41

Your equations are incomplete, because you haven't included the angle $$\theta_k$$ of the sleeve relative to the pole. For simplicity, you could take the two "jammed" angles as being symmetric, i.e. $$\theta_{k1}$$ and $$-\theta_{k1}$$.

When the toy is jammed, in equation 1 you should have a term $$c(\theta - \theta_k)$$ not your $$c\theta$$ term.

When the toy is not jammed, you have a 3-DOF system, not two, with the variables $$\theta$$, $$\theta_k$$, and $$z$$.

You also need a criterion to tell when the toy is jammed and when it is not. A simple way to model that would be to assume the pole is "infinitely rough," and if $$|\theta| = \theta_{k1}$$ the toy is jammed. A more realistic model might include using Coulomb friction to model the vertical friction force on the pole as a function of the reaction forces on the ring.

When the toy is jammed, you can calculate the moment that the ring is applying to the pole. When that moment becomes zero, the ring will "unjam" and lose contact with the pole. The motion is then governed by the "free" equations of motion until $$|\theta| = \theta_{k1}$$ and the ring makes contact with the pole again.

One interesting question is whether in the limit the motion is periodic - i.e. there is a cycle of jammed and unjammed motions that repeats, with the same time interval and the same vertical drop down the pole for each repetition. The complete cycle might contain more than one jammed section, depending on the amplitude of the bird's motion.

BUT … all of the above is an over-simplification, because it ignores the relative vertical motion of the bird and the ring!

At this point, I'm inclined to say "including that is left as an exercise for the reader" :)

• Thank you for the reply! You made a intresting point there about the friction. I was aware that this model was quite simplified, so I posted a new question with a more sophisticated model using three DOF. Since friction, slip and stick are hard to model, I am thinking of starting with a rod that only has stiction, and then move on to the more realistic case you pointed out. Please take a look at my new question; the model is better. But I still have trouble understanding how to derive the equations of motion. The reports make a similar comment to yours: "We leave it to the reader!" Feb 7, 2019 at 20:02