Report on the derivation of the solid pendulum with mass-spring system demonstration

Nasser M. Abbasi, updated August 28, 2012

Physical description of the problem

The problem is described in the following diagram

report_1.gif

A solid pendulum has a mass-less tube at the end of it with spring around it that has a bob at its end. Hence the mass-spring can only move in a horizontal direction inside the tube. There are 2 degrees of freedom for this system x(t) and θ(t) . The Lagrangian is derived for the case if a spring is present of not.

Derivation for the solid pendulum with the case when the spring is present

The kinetic energy KE is given by

report_2.gif

report_3.gif

and the potential energy is given by (assuming the pivot at zero potential, and negative potential is below that)

report_4.gif

report_5.gif

The Lagrangian is found and the differential equations derived

report_6.gif

report_7.gif

report_8.gif

report_9.gif

Parameters are given values, NDSolve is used to generate numerical solution

report_10.gif

report_11.gif

report_12.gif

report_13.gif report_14.gif

Derivation for the solid pendulum with the case when spring is missing

The kinetic energy remain the same and is given by

report_15.gif

report_16.gif

and the potential energy is given by (assuming the pivot at zero potential, and negative potential is below that)

report_17.gif

report_18.gif

The Lagrangian is found and the differential equations derived

report_19.gif

report_20.gif

report_21.gif

report_22.gif

Parameters are given values, NDSolve is used to generate numerical solution

report_23.gif

report_24.gif

report_25.gif

report_26.gif report_27.gif

Derivation for the massless pendulum with the spring present

This case was not part of the demo. But I include the model here for illustration. In this case, the  pendulum is assumed massless. As follows

report_28.gif

In this case, for some reason, singularity can be generated by NDSolve under some very specific conditions.

Kinetic energy

report_29.gif

report_30.gif

Potential energy

report_31.gif

report_32.gif

Lagrnagian is found and the 2 equations derived

report_33.gif

report_34.gif

report_35.gif

report_36.gif

NDSolve is used to solve and solution is plotted

report_37.gif

report_38.gif

report_39.gif

report_40.gif report_41.gif

Singularity exist under some specific conditions. Needs more investigation

report_42.gif

report_43.gif

NDSolve::ndsz:

Spikey Created with Wolfram Mathematica 8.0