Comment on First Results for the Milankovitch Cycles

This is a very good method which uses spherical polar coordinates for the two masses of a dumbbell. In far infra red rotational spectroscopy it would be a diatomic molecule such as nitrogen or oxygen, a linear rotating molecule with specific moments of inertia with well defined rotational kinetic energy. The extra ingredient in your calculation is that the gravitational file on both masses is considered. The earth is a symmetric top, so gravitational fields on points on the earth’s surface can be considered. For clarity, I will write out your method in a note.

Sent: 13/02/2017 13:49:00 GMT Standard Time
Subj: Re: First Results for the Milankovitch Cycles

See sketch: The model is two dumbbell masses moving freely around the middle point (centre of mass) of their connection (angles theta and phi). Their distance to the middle point is h=const. The middle point moves with coordinates r, theta_1, phi_1 around the central mass M (at rest). A So this is a two-body problem with vectorial coordinates

r_1 = R + h_1,
r_2 = R + h_2.


h_1 = – h_2

the kinetic energies of both masses are the same and the total kinetic energy takes a simple form. If the correct distances to the centre are taken for the potential energy, the calculation becomes unhandable, therefore I took the simplified form

U = – m M G / r

for both masses. In total we have 5 degrees of freedom or independent Lagrange equations. They are compiled in section 6 of the protocol. The subsystems (theta, phi) and (r, theta_1, phi_1) decouple. The solutions are oscillatory in the parameter range where the conic section of the centre of mass is an ellipse. The exact term of U would introduce a coupling between both subsystems.
The dumbbell without central motion should be equivalent to a rigid body with special moments of inertia.


Am 13.02.2017 um 13:29 schrieb EMyrone:

This looks very interesting, the problem has been essentially solved. Can you send over the lagrangian used for this solution? It will be very interesting. The result is right, the nutations and precessions are very slow compared with the orbital interval. Congratulations on this solution!

To: EMyrone
Sent: 12/02/2017 19:36:52 GMT Standard Time
Subj: dumbbell orbit

I plotted the orbit of the centre of the dumbbell (blue) and one mass of
the dumbbell (red). It can be seen that the centre of mass moves on a
plane elliptic orbit. The dumbbell shows nutation and precession. This
is the most simple model without a rigid body, requiring two mass points
with constraints. These mass points represent a ring which gives the
same Lagrange equations. This may also be a method for computing the
Milankovitch cycles where the parameters are different:
nutation/precession is very slow compared to the elliptic orbit.



  1. No trackbacks yet.

Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s

%d bloggers like this: