n-body problem
The n-body problem is an ancient, classical problem[1] of predicting the individual motions, and forces on same, of a group of celestial objects interacting with each other gravitationally. Solving this problem -- from the time of the Greeks and on -- has been motivated by the need to understand the motions of the Sun, planets and the visible stars. In the 20th century, globular-cluster star-structures became an important n-body problem too.[2] The n-body problem in general relativity is considerably more difficult.
The classical physical problem can be informally stated as: given only the present positions and velocities of a group of celestial bodies, otherwise their orbits, predict their forces owning to their interactions; and as a consequence, predict their motions for all past and future times.
To this purpose the two-body problem has been completely solved. For n = 3, solutions exist for special cases, and also for the so-called Restricted 3-body Problem discussed below. A general, classical solution in terms of first integrals is known to be impossible. An exact theoretical solution for arbitrary n can be approximated via Taylor series, but in practice such an infinite series must be truncated, giving at best only an approximate solution. In addition, many solutions by numerical integration exist, but these too are approximate solutions (see below).
This page develops a general (statics), closed-form solution for determining the reactive forces for the general n-body problem case, owing to a concentrated applied body force (an external force and moment). Rigid-body (i.e., where A 's = 1, see below) applications include distribution of useful loads in a finite element model (FEM); 3D rigid-body rivet analyses; and Astronomy problems. Spatial (Astrodynamic) problems and the like are particular solutions to this new general n-body problem solution; and once the initial positions and forces of the bodies, {mζ}, ζ = 1, 2, ... N, are known, velocities and accelerations, i.e. their motions, may be determined. Determining equations and numerical values for Astronomy problems is beyond this page's scope (see Astrodynamics). Structural analyses or soil analyses elastic solutions (A 's < 1) are possible too. Again, the latter applications are beyond the scope of this page. Note especially: rigid-body solutions for elastic structures are not valid.
Informal version of the Newton n-body problem
Some textbooks state the n-body problem problem the following way: consider point-masses , ..., in three-dimensional (physical) space. Suppose the attractive force experienced between each pair of particles is Newtonian (i.e., conform to Sir Isaac Newton's (1643-1727) Laws of Motion). Then determine the spatial position of each particle at every future (or past) moment of time if initial positions and velocities in space are specified for every particle at some instant, of time.
The above declaration is misleading because Newton, by using three orbital positions of a planet's orbit—obtained from John Flamsteed [3]—was able to produce an equation predicting the planet's motion, i.e., gave its orbital properties. Having done so they soon discovered those equations did not predict correctly some orbits and Newton realized it was because of the gravitational interactive forces amongst all the planets effecting all their orbits. Thus the need and rise of the n-body problem in the early 17th century. So, strictly speaking, the attractive forces do not conform to Newton's Laws of Motion completely; they conform to Jean Le Rond D'Alembert's non-Newtonian 1st and 2nd Principles and to the n-body problem physics (and equations) developed below (see Planet Forces and Motions Within Our Solar System below). An aside: these mathematically undefined perturbations still exist even today and the planetary orbits have to be constantly updated.
The above goes right to the heart of the matter as to what exactly the n-body problem is all about. It is not sufficient to just specify the initial position and velocity to determine a planets orbit—the gravitational interactive forces have to be known too.
Newton does not say it directly but implies in his Principia the n-body problem is unsolvable[4],.[5] Newton said (ref Cohen) in his Principia, paragraph 21: "And hence it is that the attractive force is found in both bodies. The Sun attracts Jupiter and the other planets, Jupiter attracts its satellites and similarly the satellites act on one another. And although the actions of each of a pair of planets on the other can be distinguished from each other and can be considered as two actions by which each attracts the other, yet inasmuch as they are between the same, two bodies they are not two but a simple operation between two termini. Two bodies can be drawn to each other by the contraction of rope between them. The cause of the action is twofold, namely the disposition of each of the two bodies; the action is likewise twofold, insofar as it is upon two bodies; but insofar as it is between two bodies it is single and one"...Newton concluded via his 3rd Law that "according to this Law all bodies must attract each other."
Perhaps the best way to say what the n-body problem is, is to say what was initially stated: given only the present positions and velocities of a group of celestial bodies, otherwise their orbits, predict their forces owning to their interactions; and as a consequence, predict their motions for all past and future times.
King Oscar II Prize: Historical Perspective
The problem of finding the general solution of the n-body problem was considered very important and challenging. Indeed in the late 19th century King Oscar II of Sweden, advised by Gösta Mittag-Leffler, established a prize for anyone who could find the solution to the problem. The announcement was quite specific:
Given a system of arbitrarily many mass points that attract each according to Newton's law, under the assumption that no two points ever collide, try to find a representation of the coordinates of each point as a series in a variable that is some known function of time and for all of whose values the series converges uniformly.
In case the problem could not be solved, any other important contribution to classical mechanics would then be considered to be prize-worthy.
The prize was finally awarded to Poincaré, even though he did not solve the original problem. (The first version of his contribution even contained a serious error; for details see the article by Diacu). The version finally printed contained many important ideas which led to the development of chaos theory. The problem as stated originally was finally solved by Karl Fritiof Sundman for n = 3 (see below).
Classical Mathematical formulation
Since gravity is responsible Newton reasoned for the motions of planets and stars, he expressed their gravitational interactions in terms of his Law of Universal Gravitation, but soon found, and like everyone else since, it was not a general solution to the n-body problem (see below). However, his Law was the catalysis which started the modern quest for the solution to the n-body problem.
The Law of Universal Gravitation[6] is expressed in scalar form as:
where F, a scalar, is, by Newton's 3rd Law, the force between the two masses; m1 and m2 are two mass bodies attracted to each other, assumed reduced to point-masses; and r, a scalar too, is the distance between the masses. Newton proved in his Principia a spherically-symmetric body can be modeled as a point-mass.
Notice the difference between the scalar form and the vector form:
where both F and r are vectors giving both magnitude and direction in some pre-defined Euclidean coordinate system.
The constant G is special, it is the gravitational constant. G, the constant of proportionality, is the parameter making Newton's Law emperialistic (i.e., capable of being verified or disproved by observation or experiment). To this purpose see Lord Cavendish's experiment measuring of G. Then what is G physically? G allows mass (matter phenomenon) when combined with distance (spatial phenomenon) to be linked to force ("force" phenomenon). All three phenomena have different units of measurement. Newton's equation is empirical because G is empirical, but (m1m2/r3)r is not a force and so G(m1m2/r3)r is not equal to Fg (gravitational force) -- but only equivalent to it (and that only by assuming we know what "force" is: in this case it is a field force). His equation establishes only an equivalence relationship via G: the equation maps the mass and spatial phenomena -to- force. As will be explain below, force is a functional.
G allows Newton's 2nd Law and all others too, to not require incorporating any constants of proportionality.
Either the scalar or vector forms can be set equal to Newton's 2nd Law, F = m a , where a is the directed acceleration. Notice the absents of a constant of proportional. Continuing with this purpose of equating and using the vector form, let F ij = f(mi,mj,rij), where i and j are integers, 1, 2, 3,... Otherwise, F is a function of mass and geometry. Then to employ Newton's Law Universal, let F12 = (Gm1m2/r3) r be the force on m1 owning to mass m2; likewise, let F 21 =(Gm1m2/r3)r be the force on m2 owning to m1. This way one arrives at two separate equations, and by Newton's 3rd Law, F 12 = F 21 (see the Two-body Problem below). Imagine if there were three masses involved: there would be six separate equations, F12, F21, F13, F31, F23, F 32. If four masses (or more) were involved there would be a mathematical mess.
Newton's Law of Universal Gravitation is generally considered applying only to celestial objects, but is equally useful for Earth science (calculating Ocean tides, calculating the 13-mile high ellipsoidal budge at the equator, calculating the fuel requirements verses weight for rockets, etc.). For those applications there is the "Laws of Weights".[7] Bodies weigh most at the surface of the Earth of course. Below the surface the weight decreases as the distance to the center decreases linearly; above the surface the weight decrease as to the square of the distance increases (Newton's Law).
Newton thought in terms of a change in momenta ∆m''v; where as Euler, for changing parameters as a function of other variables, thought in calculus terms, and momenta became a rate of change between variables, like F = m''a = m d2r/dt2. This was not just a change in notation started by G. W. Leibniz (1646-1716) but a change in concept. In calculus terms Newton's notation indicated a rate of change in variables—with variables wearing hat-like dots (for, say, a change in velocity with one dot), double dots for acceleration a. Newton's momentum became Euler's Eulerian acceleration with a = d2r/dt2. In what follows Newton's notation is used.
Most textbooks present the following general solution approach for the n-body problem: it means in mathematical terms finding a global solution describing the -body problem of celestial mechanics, which is an initial-value problem for ordinary, second order, differential equations. Given initial values for the positions and velocities of n particles (j = 1,...,n) with for all mutually distinct j and k , find the solution of the second order system:
Eq.
where are constants representing the masses of n point-masses, , are 3-dimensional vector functions of the time variable t, describing the positions of the point masses, and G is the constants of proportionality; i.e., gravitational constant.
This equation is Newton's second law of motion on the (left side) combined with his Law of Universal Gravitation (right side). The left-hand side is the mass times acceleration for a particular jth particle (F = ma in vector form); whereas the right-hand side is the sum of the forces on that particle relative to another particle. The forces are assumed here to be gravitational as given by Newton's law of universal gravitation; thus, they are proportional to the masses involved, and vary as the inverse square of the distance between the masses. The power in the denominator is three instead of two to balance the vector difference in the numerator, which is used to specify the direction of the force. An aside, Euler actually transformed Newton's 2nd Law equation — of changing momentum, assumed equivalent to net force - into an ordinary, second order differential equation.
A general, closed-form, n-body problem solution (i.,e., n-formulas) can not be derived employing Newton's Law of Universal Gravitation equation as a solution approach.
see General considerations for the Classical Solution below for further discussion.
Two-body problem
The two-body problem, n = 2, was completely solved by Johann Bernoulli (1667-1748) by classical theory (and not by Newton) by assuming the main point-mass was fixed, is outlined here[8] Consider then the motion of two bodies, say Sun-Earth, with the Sun fixed, then:
The equation describing the motion of mass m2 relative to mass m1 is readily obtained from the differences between these two equations and after canceling common terms gives: α + (η/r3)r = 0, where
- α is the Eulerian acceleration d2r/dt2;
- r = r2 - r1 is the vector position of m2 relative to m1;
- and η = G(m1 + m2).
This is the fundamental differential equation for the two-body problem Bernoulli solved in 1734. Notice for this approach forces have to be determined first, then the equation of motion resolved. This differential equation has elliptic , or parabolic or hyperbolic solutions[9], [10], [11]
It is incorrect to think of m1 (the Sun) as fixed in space when applying Newton's Law of Universal Gravitation, and to do so leads to erroneous results. Dr. Clarence Cleminshaw calculated the approximate position of the Solar System's true barycenter, a result achieved mainly by combining only the masses of Jupiter and the Sun. Science Program stated in reference to his work: "The Sun contains 98 per cent of the mass in the solar system, with the superior planets beyond Mars accounting for most of the rest. On the average, the center of the mass of the Sun-Jupiter system, when the two most massive objects are considered alone, lies 462,000 miles from the Sun's center, or some 30,000 miles above the solar surface! Other large planets also influence the center of mass of the solar system, however. In 1951, for example, the systems' center of mass was not far from the Sun's center because Jupiter was on the opposite side from Saturn, Uranus and Neptune. In the late 1950s, when all four of these planets were on the same side of the Sun, the system's center of mass was more than 330,000 miles form the solar surface, Dr. C. H. Cleminshaw of Griffith Observatory in Los Angeles has calculated."[12] This means the Sun wobbles and Sun-spots are possibly caused via the movement of the barycenter, owing to Jupiter's 11-year cycles, producing Sun-spots every 22 years. It further needs to be pointed out the total mass orbiting the Sun is probably equal to the Sun's own mass.
![](http://upload.wikimedia.org/wikipedia/commons/thumb/1/13/Kipler%27s_Error.jpg/220px-Kipler%27s_Error.jpg)
The Sun wobbles as it rotates around the galactic center, dragging the Solar System and Earth along with it. What mathematician Kepler did in arriving at his three famous equations was curve-fit the apparent motions of the planets using Tycho Brahe's data, and not curve-fitting their true circular motions about the Sun (see Figure). Both Robert Hooke and Newton were well aware Newton's Law of Universal Gravitation did not hold for the forces associated with elliptical orbits.[13] In fact, Newton's Universal Law doesn't account for the orbit of Mercury, the Asteroid Belt's gravitational behavior, or Saturn's Rings.[14] Newton stated (in the 11th Section of the Principia) the main reason however for failing to predict the forces for elliptical orbits was his math model was for a body confined to a situation that "hardly exist in the real world," namely, the "motions of bodies attracted toward an unmoving center." Some present physics and astronomy textbooks don't emphasize the negative significance of Newton's assumption and end up teaching that his math model is in effect reality. It is to be understood the classical two-body problem solution above is a mathematical toy. See also Kepler's first law of planetary motion.
Newton conveniently fixed the Sun so he could do simple calculations (in effect he cheated), but all following after him have also made the same mistake by analytically fixing the Sun too. They perpetuated the mathematical toy (see Truesdell's Essays in the History of Mechanics referenced below). An aside: Newtonian physics doesn't include (among other things) relative motion and may be the root of the reason Newton "fixed" the Sun.[15][16] The Sun's wobbling means the n-body problem's solution realistically is much more complicated than maybe previously thought (see below).
Three-body problem
In the past not much was known about the n-body problem for n equal to or greater than three.[17] The case for n = 3 has been the most studied. Many earlier attempts to understand the Three-body problem were quantitative, aiming at finding explicit solutions for special situations.
- In 1687 Isaac Newton published in the Principia the first steps in the study of the problem of the movements of three bodies subject to their mutual gravitational attractions, but his efforts resulted in verbal descriptions and geometrical sketches; see especially Book 1, Proposition 66 and its corollaries (Newton, 1687 and 1999 (transl.), see also Tisserand, 1894).
- In 1767 Euler found collinear motions, in which three bodies of any masses move proportionately along a fixed straight line. The circular restricted three-body problem is the special case in which two of the bodies are in circular orbits (approximated by the Sun-Earth-Moon system and many others).
- In 1772 Lagrange discovered two classes of periodic solution, each for three bodies of any masses. In one class, the bodies lie on a rotating straight line. In the other class, the bodies lie at the vertices of a rotating equilateral triangle. In either case, the paths of the bodies will be conic sections. Those solutions led to the study of central configurations , for which for some constant k>0 .
- A major study of the Earth-Moon-Sun system was undertaken by Charles-Eugène Delaunay, who published two volumes on the topic, each of 900 pages in length, in 1860 and 1867. Among many other accomplishments, the work already hints at chaos, and clearly demonstrates the problem of so-called "small denominators" in perturbation theory.
- In 1917 Forest Ray Moulton published his now classic, An Introduction to Celestial Mechanics (see references) with its plot of the Three body restricted Solution (see below).[18][19]
![](http://upload.wikimedia.org/wikipedia/commons/f/f9/N-body_problem_%283%29.gif)
Specific solutions to the Three-body problem result in chaotic motion with no obvious sign of a repetitious path.[citation needed]
The restricted problem (both circular and elliptical) was worked on extensively by many famous mathematicians and physicists, most notably by Poincaré at the end of the 19th century. Poincaré's work on the Restricted Three-body Problem was the foundation of deterministic chaos theory.[citation needed] In the restricted problem, there exist five equilibrium points. Three are collinear with the masses (in the rotating frame) and are unstable. The remaining two are located on the third vertex of both equilateral triangles of which the two bodies are the first and second vertices. See Figure below.
This may be easier to visualize if one considers the more massive body (e.g., Sun) to be "stationary" in space, and the less massive body (e.g., Jupiter) to orbit around it, with the equilibrium points maintaining the 60 degree-spacing ahead of, and behind, the less massive body almost in its orbit (although in reality neither of the bodies are truly stationary, as they both orbit the center of mass of the whole system—about the barycenter). For sufficiently small mass ratio of the primaries, these triangular equilibrium points are stable, such that (nearly) massless particles will orbit about these points as they orbit around the larger primary (Sun). The five equilibrium points of the circular problem are known as the Lagrange points. See figure below:
![](http://upload.wikimedia.org/wikipedia/commons/9/97/Restricted_3-Body_1.jpg)
In the Restricted 3-Body Problem math model figure above (ref. Moulton), Lagrangian points L4 and L5 are where the Trojan planetoids resided; m1 is the Sun and m2 is Jupiter. L2 is where the asteroid belt is. It has to be realized for this model, this whole Sun-Jupiter diagram is rotating about its barycenter. The h-circles and closed loops echo the electromagnetic fluxes issued from the Sun and Jupiter. (An aside, see Meirovitch, pages 414 and 413 for the Three-body Problem.)
The Restricted Three-body Problem assumes the mass of one of the bodies is negligible[citation needed]; For a discussion of the case where the negligible body is a satellite of the body of lesser mass, see Hill sphere; for binary systems, see Roche lobe; for another stable system, see Lagrangian point.
Practical n-body Problem Theoretical Solution
This solution occurred while trying to distribute useful applied loads (i.e., human body loads, luggage, attached extra hardware, etc.) to the grids (spatial points) belonging to surrounding basic structures of finite element models. How else would one distribute applied loads if not by a n-body problem type solution? The similarity for determining the reactive loads for a 3D-rivit analyses also became apparent.
Turns out any solution to the n-body problem requires using Jean Le Rond D'Alembert's 1st and 2nd Principles, is key.[20]
Hookean Switches and Springs, Point-Masses and the Inertia Tensor
Index notation is used in this section.[21][22] Indexes and scalars are plain italic text; vectors are in bold, italic text; Cartesian tensors[23] are in plain italic text.
It is convenient to define a vector operator, a functional, to be employed in the solution of the n-body Problem.
Hookean Switches: Let there be a vector function, 0 ≤ Aζi ≤ 1, where i := x, y, z or 1, 2, 3, are the component directions, ζ := 1,...N, is the point-masses' mζ coordinate number. The Aζi's are functionally the feasible, translational reaction component's allowable degrees of freedom (DOF) at Pζ⇒ (xζ, yζ, zζ), specified for any and all point-masses mζ, where the N point-masses mζ are defined in a right-handed, Euclidean coordinate system in File:Reals R.jpg3 with (x, y, z) coordinates. Coupling forces are reacted at these point-masses too and are subject to the same feasible translational constraints.
Applying Hookean switches A 's at the supports or reaction points for every feasible translational reaction, i.e., for each Pζ, would be the same as showing those feasible or allowable translational vectors visually on a free-body diagram; but instead of showing by default all three DOF possible translational vectors, say ( P ζ )i, with i = 1,2,3. There must be at least three point-masses for stability, having a total of six DOF , meaning there must be at least six feasible reactions for these three points; i.e., six Aiζ's. (This latter requirement is the same as balancing an aircraft's FEM on three points and then applying Jean Le Rond D'Alembert's two Principles, resulting in no load on the three points.)
(Definition continues.) Let the radius vector be:
Eq.
where i, j, k with the "hats" are the unit directional vectors, and the A 's are the Hookean switches or springs defined above. For example, for mζ at Pζ ⇒ (xζ, yζ, zζ): if all Aζi = 0 then r ζ = 0; if all Aζi = 1 then there are three reaction components at ζ. These conditions describe Hookean switches. All other values of Aζi ≤ 1, but not zero, describe Hookean springs. For these conditions, if the reaction points' material is structural; i.e., elastic ( σ = E ε ), then it will be further assumed the structure will conform to a similar linear law, say Hooke's Law, ( F = k x ) of load-deformation proportionality, where k is the spring rate, and wherein the A 's are analogous to k (this assumption is not entirely necessary as F = f(x) can be any functional without its exact definition given).
(Definition continues.) Then for each point Pζ there may be 0, or up to 3 DOF, meaning mζ may be, depending on Aζi's components for any particular DOF in that direction, elastic or rigid body. If r ζ is defined with Hookean Switches then r ζ are said to be a Hookean vector.
The concept is: Hookean switches allow point-masses to respond, or not respond, as reaction forces in selected directions, by setting selected Aiζ's components to 0; and vice-verse. When all the feasible Aζi = 1 the solution is a rigid body solution. In this way Hookean switches allows for the solution of the n-body Problem. Hookean springs allows a type of finite element model (FEM) to be constructed (i.e., 3D rivet analyses).
Let M := m1 + m2 + m3 + ...+ mN, for ζ = 1, 2, ..., N point-masses. Of course for each mζ there is a spatial point Pζ ⇒ (xζ, yζ, zζ).
Let r ζi at Pζ, where ζ = 1,...N mass points, be a Hookean vector for each mζ in all (x, y, z) space, for a right-hand, rectangular Cartesian coordinate system. The A 's here are the components of A. Then the six massless barycenters (six naturally develop when using index notation) for {M} ≤ 3 are, by simple statics:
Eq.s
If the reaction or support pattern is semi-symmetrical then the centroid set will be less than six, as, say when for example, xij = xji; if the pattern is entirely symmetrical then these six centroids will revert back to the usual three. Continuing, the massless centroidal mass moments of inertia for N masses, mς, again by simple statics are:
Eq.s
where, say x barxy, etc, are the barycenters calculated above.[24][25]
But before that...qualification for a right-hand coordinate system arises because some USA Aerospace firms use the left-hand system, making their signs for the cross products positive. If the inertia tensor if symmetric and positive-definite (via the value of the its determinant), then the coordinate system is right-handed ― always check this, and if not then your signs are wrong.
Applied forces and moments may be Hookean too. These Hookean definitions (and concept) are missing from physics textbooks, but are somewhat similar in concept to grid-point releases in a finite elements analyses. Hookean points are not Lagrangian coordinates.
The Force
Interestingly, if empirical force laws and Euclidean geometry are true, and the Principle of Equilibrium is true of course, then any logical combination therein is also true. In all that follows is naught but simple physics (statics).
If, for example, when the vector si = Є ijk rj θ k (analogous to s = rθ) from simple geometry is combined with D'Alembert's two Principles R = 0 and M = 0 , then a formula called The Force, can be developed giving the reactive forces for a pattern of points subject to an external resultant force F i and moment M0i. rj may be a Hookean vector, and si covers all space. It will be shown The Force can be used to completely solve any n-body Problem for any n ≥ 3.
Let the radius vectors to the locations of the ζ = 1, 2, ...N reaction points or supports points be the set r ⇒{rζi)|i = 1, 2, 3} for mζ ≥ 3 mass-points, react to an applied resultant force F i and moment ( M 0)i, i = 1, 2, 3 components, at File:C g s.jpg. Then reaction displacements (∆sζ)i of any point (r ζ)i owing to angular rotations (∆θζ)i are:
Eq.
where File:Delta S 01.jpg is a vector. Assume only infinitesimal angular rotations (∆θζ)i are associated with reaction point motions File:Delta S 01.jpg in any instant of time. Further, since this is an instant of time (time is fixed) the solution becomes one of statics (i.e., quasi-steady). Remember, ζ := 1,...N, is the point-mass mζ coordinate number, where m ≥ 3.
The load center is defined to be at File:C g s.jpg, who coordinates are basic input. File:Centroids 01.jpg are use later to calculate the sum of the moments — are equal to the {Pζ} coordinate set for ζ = 1, ... N reaction points, and are the centers of gravity of those point-masses. Real masses may be reduced to point-masses the same way Newton did it. In this analysis physical real mass is not a variable; only forces, reactions, moments, and geometry are real variables.
The Hookean File:Delta S 01.jpg equations in expanded matrix form are exactly analogous to M = r x F or bisor (second order tensor) Mi = ЄijkrjFk.
Assume infinitesimal moments (torques) M ζi are used in lieu of ∆θζi. Expressed mathematically:
P ζi are the three force components at the support or reaction points. File:Rbar 01.jpgζi are the distances from the support patterns centroid to the support points, which is a detail explained more fully below. This angular rotations -to- moments assumption will be incorporated below. Interestingly, this same assumption is employed in the displacement method in finite element analysis, but there in a different context. (Other assumption relationships may be employed too.) The left hand side expanded in matrix form is:
The [f(File:Rbar 01.jpg)]3,3 matrix is actually the inertia tensor of the reaction or support pattern, a mathematical detail given below. What follows next is to determine the set of distances File:Rbar 01.jpgςi from the centroids of the reaction or support pattern to the individual mass-points. To that purpose let:
The radius vector rςi are rays to point-masses mς at B; the radius vector r o i connects to the tail of the position vector File:Rbar 01.jpgζi at A; and the position vector in turn connects to the point-masses at B; and closes the loop (see figure), where:
![](http://upload.wikimedia.org/wikipedia/commons/thumb/b/b6/R_bar_equatons.jpg/220px-R_bar_equatons.jpg)
- File:Rbar 01.jpgζi, i = 1, 2, 3, are the position vectors from the six centroids of the reaction pattern -to- the reaction points or supports points (or planets, etc.). The three centroids really consist of the six centroids as developed in the Section directly above: File:Rbar Converted.jpg.
- r0i are the radius vectors (rays from the origin to...) to the reaction or support pattern's mass-points at {Pζi}.
- rζi are the distances from the origin to the reaction pattern's point-masses.
These equations in expanded form are:
Also note, mς can be the centroids of reduced point-masses, but usually are just a mζ 's coordinate.
Suffice it to say when the reactions P ςi are multiplied by File:Rbar 01.jpgςi, theFile:Rbar 01.jpgςi takes the form of a skewed Hermitian matrix:
Eq.
And that product times the inverse of the inertia tensor [Iς]−1, is equal to the internal moment Mςi:; or re-written in terms of reactions:
P ζi are the individual coupling reactions from the reduced moments, plus the additional translational reactions distributed about the reaction points mζi owing to an externally applied force F ζi and Moi. The resultant internal moments M ζi as coupling reactions at each point are equal to the sum of resultant external moments at the centroids of the reaction pattern, by equilibrium considerations; i.e., by D'Alembert's Second Principle (and not by Newton's Third Law, which would be incorrect reasoning).[26]
(Explanatory note.) Some explanation is in order before preceding: First, referring to the cantilever beam with the off-set load, if the load F1 is moved to the cantilever's tip's end along the shaft's axis, then (to calculate loads at the root) if the original force, F1, is reduplicated and transferred to the tip as F 2, then it has to be balanced by an up-load, F3, for the two to remain in equilibrium. This up-load F3, combined with the original down-load F1, is a couple and equivalent to a moment M (torque). Second, at the root of the cantilever the reaction forces at the attachment points (not shown) not only react with the vertical force F2 in pure upward vertical shear divided by the number of attachment points, but by horizontal reaction couples in- and out- of the wall too. The new couple (moment M ) is reacted by only reaction couples at the root parallel to the wall. In this way (symbolically) the external moment Mexternal = M internal = rFreaction; i.e., the external moment(s) equal the in-plane and out-of-plane static reaction couples at the attachments in the wall (D'Alembert's Second Principle). Looking at the big picture, the total external forces and moments are equal to the sum of all the internal coupling reactions (moments) plus the direct shear reactions (translational reactions) divided by the number of attachments. D'Alembert's two Principles combined.
(Continuing.) Let the reactions Pςi above owe to a known external applied force Fi and known applied external moment Moi. The Fi when moved from its initial position to the centroids of the reaction or support pattern, creates moments and these moments are of the form [xL]3,3 F i. The [xL]3,3 is developed now:
and rpi ⇒File:C g s.jpg is the position of the applied force and moment. Then in expanded scalar form:
Eq.
[xL]3,3 is a skewed Hermitian.
Putting all the above together results in an equation call The Force. The Force formula gives a set of balanced reactions at the reaction points or support points (Lb.). Its proof is by calculating equilibrium again, once those reactions have been found (see Equilibrium of Forces and Moments section below).
Eq.
where i = x, y, z; ζ = 1, 2, ..., N.
This solution takes an external force and moment and moves it over to the centroids of the reactions or supporting pattern (system); and by incorporating the massless inertial properties of that pattern, beams the applied loads to those reaction points (or support points). Although somewhat messy algebraically, it's that simple. It's analogous to a 3D rigid body, rivet analysis.
The above in expanded matrix form, except for the inverse of the massless inertia tensor [ I ς]−1, for the three orthogonal reaction components at a ζ-reaction point is:
Eq.
In the expanded version (except for [ I ς]−1 terms) one can clearly see the right side of the equation is composed of two parts: the translational forces and the coupling forces.
In the two above equations:
- ζ = 1, 2, 3, ... N point-masses mζ coordinate number; brackets [ ] and { } represent matrices and -1 is for inverse; and subscript L stands for Load position, s stands for supports;
- ( P ζ )i are the vector reactions at each support point Pi⇒ (xi, yi, zi);
- ( F L)i are the three components of the applied force (lb.), transposed; ( M 0)i are the the three components of the applied moment (lb.-in.), transposed; F and M o are basic inputs;
- The load centers at File:C g s.jpg is basic input (in.).
- [ I ς]−1 is the inverse of the massless inertia tensor of the reaction pattern (via Crout's Method), is a symmetrical 3x3 matrix (in.−2);
- [xL]i is the position of the applied force and moment, 3x3 matrix; [xL]i( F L) are the resulting couples in the support pattern caused by moving the applied force over to the pattern's centroids (in.).
- [xs]i are the centroids of the reaction pattern, 3x3 matrix (in.);
- [xs]i[ I ς]−1[[xL]i{ F L}iT + { M o}iT]3,3] are the unbalanced moments reacted as couples, 3x3 matrix (lb.);
- {( F L)i/(Σ(Aζ)i)}T1,3 are the symmetrical (translational) reactions (lb.).
The [xs]i[Iς]−1[xL] term is a congruent transformation (Meirovitch) and for all physical entities having like congruent transformations, like stress, strain, they are always positive definite if the right-hand rule is used: a fundamental property of all matter.
{Aζ}i are the Hookean switches defined above. If a particular motion is planar then out-of-plane reaction forces (P ζ)i will be zero by default. ( P ζ)i has been computerized, see below.
Equilibrium of Forces and Moments
The sum of the forces for the entire support or reaction point set, {Pζ} set, for ζ = 1, ... N reaction points, are:
Eq.
The applied load center is at File:C g s.jpg and is basic input. The c.g.'s (centers of gravity) of the reaction points (the mζ's reduced to a support-point set of point-masses) are at File:Centroids 01.jpg, i = 1, 2, ...,N; are in most cases the same set as the {Pζ} set for ζ = 1, ... N reaction points. Then the sum of the moments for the entire supports or reactions are:
Eq.
Numerical Example
Example problem: a large video game monitor weighing 179 pounds is mounted into a cantilevered (from the floor) flying (airborne) console, is supported internally by the console by 20 support (attachment) points. The external applied force, has a load c.g. of File:C g s.jpg⇒(3.96442, 13.5625, 44.89365) inches, is close to these support points. The problem is to determine three rigid-body, unit-load sets for the entire support pattern. In this example, only (as an example) a unit load of 1.0 Lb. in the vertical direction is applied to determine its unit reactions. (Two more unit load runs would ordinarily be needed to determine sideways and forward unit reactions too. Otherwise, three sets of unit-load reactions are needed to obtain a complete loading profile for, say, a FEM.) Support coordinates in inches are given in Tables 1 and 2; from which the six support c.g.'s as well as the inertia components are calculated (only once) in the Tables below.
Extracting values from these tables gives:
And the inertia tensor is:
And in matrix form:
The load transfer matrix is:
![](http://upload.wikimedia.org/wikipedia/commons/3/3c/Transfer_Load_XL.jpg)
And the support transfer matrix which transfer the applied moment as couples to the support points, is:
All internal reactions at the N-support points are given by the following three equations:
Reaction moments equations:
Planet Forces and Motions Within Our Solar System
If we define the Solar System's domains as having expanded as far as the gravitational fields of the Sun, Jupiter, and the other outer-planets, then the Solar System, beyond the planets Neptune and Pluto, extends perhaps a thousand times further than we have yet explored. But out there we know there are billions of small objects in that zone known as Oort Cloud. If all the known mass except the Sun — that is, the planets, asteroids, moons, Trojan planetoids, etc. — is subtracted from the Sun's known mass, then it is conjectured: the remaining mass of the Solar System would be equal to all that mass in the zone, such that, that mass plus the known body masses equals the Sun's mass. This is not only reasonable gravitationally, but application of the general n-body Problem indicates it is true too, because without this additional mass — the planets, etc., has too little mass to absorb the Sun's gravitational forces — result in higher interactive forces via the n-body solution calculations, than those via Newton's Law of Universal Gravitation. These higher forces are reduced to Newtonian forces by assuming the zone 's mass is almost equal to the Sun's.
Calculating Solar Forces: suppose high above the Solar plane a picture was taken by a satellite, instantly freezing all motion of all bodies in the Solar System. This would result in a quasi-steady (static) picture of fixed positions, zero velocities and accelerations. The positions become variables (coordinates) for input into the general solution of the n-body Problem; and the main masses are already known. With this knowledge, a suggested Solar System general n-body type solution may be effected by:
- Calculate all forces of known bodies via Newton's Law of Universal Gravitation, whose values are to be used one at a time discretely as applied external forces; and as a comparisons of the n-body solution forces -to- Newtonian forces.
- At the inner diameter of Oort's Cloud establish, say, a hundred (or more) mass-points uniformly distributed around the inside of its sphere, whose total mass sum is equal to the Sun's mass, less the known masses.
- Use the quasi-steady positions of the known masses combined with the mass-points of Oort's Cloud mass-points, as input to the n-body program.
- Run the n-body computer program first using the Sun as the known force; and the planets, etc., as unknown reactions; and record the position of the Solar System's barycenter.
- Repeat a similar run for each known mass-point, with one mass-point selected as the known external applied force; (use Newton's Law to determine the applied force).
- Repeat computer runs time-wise incrementally until the orbit of the Solar System's barycenter is known (shades of Cleminshaw).
If these forces are known then the eccentricities of the planet's orbits can be determined. This task is way to large for it to be presented here.
n-body particle simulations can be effected via computer.[27] This program demonstrates body-to-body interactions, i.e., almost head-on collisions, almost always results in bodies orbiting bodies rather than collisions.
More General considerations for the Classical Solution
For every solution of the problem, not only applying an isometry or a time shift but also a reversal of time (unlike in the case of friction) gives a solution as well.
In the physical literature about the -body problem ( ≥ 3), sometimes reference is made to the impossibility of solving the -body problem (via employing the above approach). However, care must be taken when discussing the 'impossibility' of a solution, as this refers only to the method of first integrals (compare the theorems by Abel and Galois about the impossibility of solving algebraic equations of degree five or higher by means of formulas only involving roots).
Employing Newton's generalized Law of Universal Gravitation formula above, the -body problem contains 6 variables, since each point particle is represented by three space (displacement) and three momentum components. First integrals (for ordinary differential equations) are functions that remain constant along any given solution of the system, the constant depending on the solution. In other words, integrals provide relations between the variables of the system, so each scalar integral would normally allow the reduction of the system's dimension by one unit. Of course, this reduction can take place only if the integral is an algebraic function not very complicated with respect to its variables. If the integral is transcendent the reduction cannot be performed.
The -body problem has 10 independent algebraic integrals:
- three for the center of mass;
- three for the linear momentum;
- three for the angular momentum;
- one for the energy.
This allows the reduction of variables to 6 − 10. The question at that time was whether there exist other integrals besides these 10. The answer was given in the negative in 1887 by Heinrich Bruns:
Theorem (First integrals of the -body problem) The only linearly independent integrals of the -body problem, which are algebraic with respect to , and are the 10 described above.
(This theorem was later generalized by Poincaré). These results however do not imply there does not exist a general solution to the -body problem, or that the perturbation series (Lindstedt series) diverges. Indeed Sundman provided such a solution by means of convergent series. (See Sundman's theorem for the 3-body problem).
(Classical) Power Series Solution
This is the most elemental way classically of solving the n-body Problem. The theoretical expression is often called "The n-body problem by Taylor series", which is an implementation of the Power series solution of differential equations.
We start by defining the differential equations system[citation needed]:
where (in ) the upper index indicates the second derivative with respect to time , represents the number of each body and the coordinate.
Because and are given as initial conditions, then every are known. Doing implicit derivation over every results in which at are known because each depends on known precalculated and given constants and then the Taylor series are constructed theoretically in such way, performing this process infinitely.
Numerical integration
n-body problems can be solved by numerically integrating the differential equations of motion. Many different ways to do this to varying degrees of accuracy and speed exist.[28]
The simplest integrator is the Euler method, but this is only first order. A second order method is leapfrog integration, but higher-order integration methods such as the Runge–Kutta methods can be employed. Symplectic integrators are often used for n-body problems.
Numerical integration has a time complexity of O(n2), but tree structured methods, such as Barnes-Hut simulation, can improve this to O(n log n), or even to O(n) such as with the fast multipole method.
Sundman's Theorem for the 3-body Problem
In 1912, the Finnish mathematician Karl Fritiof Sundman proved there existed a series solution in powers of for the 3-body Problem. This series is convergent for all real t, except initial data which correspond to zero angular momentum. However these initial data are not generic since they have Lebesgue measure zero.
An important issue in proving this result is the fact that the radius of convergence for this series is determined by the distance to the nearest singularity. Therefore it is necessary to study the possible singularities of the 3-body Problems. As it will be briefly discussed below, the only singularities in the 3-body Problem are:
- binary collisions,
- triple collisions.
Now collisions, whether binary or triple (in fact of arbitrary order), are somehow improbable since it has been shown they correspond to a set of initial data of measure zero. However there is no criterion known to be put on the initial state in order to avoid collisions for the corresponding solution. So Sundman's strategy consisted of the following steps:
- He first was able, using an appropriate change of variables, to continue analytically the solution beyond the binary collision, in a process known as regularization.
- He then proved triple collisions only occur when the angular momentum L vanishes. By restricting the initial data to he removed all real singularities from the transformed equations for the 3-body problem.
- The next step consisted in showing that if , then not only can there be no triple collision, but the system is strictly bounded away from a triple collision. This implies, by using Cauchy's existence theorem for differential equations, there are no complex singularities in a strip (depending on the value of L) in the complex plane centered around the real axis (shades of Kovalevskaya).
- The last step is then to find a conformal transformation which maps this strip into the unit disc. For example if (the new variable after the regularization) and if [clarification needed] then this map is given by:
This finishes the proof of Sundman's Theorem. Unfortunately the corresponding convergent series converges very slowly. That is, getting the value to any useful precision requires so many terms, that his solution is of little practical use.
A Generalized Sundman Global Solution
In order to generalize Sundman's result for the case n > 3 (or n = 3 and c = 0[clarification needed]) one has to face two obstacles:
- As it has been shown by Siegel, collisions which involve more than two bodies cannot be regularized analytically, hence Sundman's regularization cannot be generalized.
- The structure of singularities is more complicated in this case: other types of singularities may occur (see below).
Lastly, Sundman's result was generalized to the case of n > 3 bodies by Q. Wang in the 1990s. Since the structure of singularities is more complicated, Wang had to leave out completely the questions of singularities. The central point of his approach is to transform, in an appropriate manner, the equations to a new system, such that the interval of existence for the solutions of this new system is .
Singularities of the n-body problem
There can be two types of singularities of the n-body problem:
- collisions of two or more bodies, but for which q(t) (the bodies' positions) remains finite. (In this mathematical sense, a "collision" means that two point-like bodies have identical positions in space.)
- singularities in which a collision does not occur, but q(t) does not remain finite. In this scenario, bodies diverge to infinity in a finite time, while at the same time tending towards zero separation (an imaginary collision occurs "at infinity").
The latter ones are called Painlevé's conjecture (no-collisions singularities). Their existence has been conjectured for n > 3 by Painlevé (see Painlevé's conjecture). Examples of this behavior have been constructed by Xia[29] and Gerver.
See also
- Einstein–Infeld–Hoffmann equations
- Few-body systems
- Natural units
- n-body choreography
- Virial theorem
Notes
- ^ Leimanis and Minorsky: Our interest is with Leimanis, who first discusses some history about the n-body problem, especially Ms. Kovalevskaya's ~1868-1888, twenty-year complex-variables approach, failure; Section 1: The Dynamics of Rigid Bodies and Mathematical Exterior Ballistics (Chapter 1, the motion of a rigid body about a fixed point (Euler and Poisson equations); Chapter 2, Mathematical Exterior Ballistics), good precursor background to the n-body problem; Section 2: Celestial Mechanics (Chapter 1, The Uniformization of the Three-body Problem (Restricted Three-body Problem); Chapter 2, Capture in the Three-Body Problem; Chapter 3, Generalized n-body Problem).
- ^ Heggie and Hut have many tables and charts, etc.
- ^ See David H. and Stephen P. H. Clark's The Suppressed Scientific Discoveries of Stephen Gray and John Flamsteed, Newton's Tyranny, W. H. Freeman and Co., 2001. A popularization of the historical events and bickering between those parties, but more importantly about the results they produced.
- ^ See Principia, Book Three, System of the World, "General Scholium," page 372, last paragraph. Newton was well aware his math model did not reflect physical reality. This edition referenced is from the Great Books of the Western World, Volume 34, which was translated by Andrew Motte and revised by Florian Cajori. This same paragraph is on page 1160 in Stephen Hawkins' huge On the Shoulders of Giants, 2002 edition; is a copy from Daniel Adee's 1848 addition. Cohen also has translated new editions: Introduction to Newton's 'Principia' , 1970; and Isaac Newton's Principia, with Varian Readings, 1972. Cajori also wrote a History of Science, which is on the Internet.
- ^ See I. Bernard Cohen's Scientific American article listed in the References.
- ^ Newton's Law of Universal Gravitation can be derived via the potential function V(r) = m/r, where m is mass and r is the spatial vector. See Brouwer and Clemence book: their Chapter III (the heart of the book), Gravitational Attraction Between Bodies of Finite Dimensions, assumes Newton's Law of Universal Gravitation in order to prove it yields the potential function V(r) = m/r. Going the other way (derivation not in their book) and deriving Newton's Law raises big questions as to which is more fundamental, Newton's Law; or the potential function and LaPlace's equation?
- ^ The Elements of Mechanical And Electrical Engineering, Volume I, Chapter Elementary Mechanics, pp. 318-319, The Colliery Engineer Co., 1898.
- ^ See Bate, Mueller, and White: Chapter 1, "Two-Body Orbital Mechanics," pp 1-49. These authors were from the Dept. of Astronautics and Computer Science, United States Air Force Academy. See Chapter 1. Their textbook is not filled with advanced mathematics.
- ^ For the classical approach, if the common center of mass (i.e., the barycenter) of the two bodies is considered to be at rest, then each body travels along a conic section which has a focus at the barycenter of the system. In the case of a hyperbola it has the branch at the side of that focus. The two conics will be in the same plane. The type of conic (circle, ellipse, parabola or hyperbola) is determined by finding the sum of the combined kinetic energy of two bodies and the potential energy when the bodies are far apart. (This potential energy is always a negative value; energy of rotation of the bodies about their axes is not counted here)
- If the sum of the energies is negative, then they both trace out ellipses.
- If the sum of both energies is zero, then they both trace out parabolas. As the distance between the bodies tends to infinity, their relative speed tends to zero.
- If the sum of both energies is positive, then they both trace out hyperbolas. As the distance between the bodies tends to infinity, their relative speed tends to some positive number.
- ^ For this approach see Lindsay's Physical Mechanics, Chapter 3, "Curvilinear Motion in a Plane," and specifically paragraph 3-9, "Planetary Motion"; and continue reading on to the Chapter's end, pp. 83-96. Lindsay presentation goes a long way in explaining these latter comments for the fixed 2-body problem; i.e., when the Sun is assumed fixed.
- ^ Note: The fact a parabolic orbit has zero energy arises from the assumption the gravitational potential energy goes to zero as the bodies get infinitely far apart. One could assign any value to the potential energy in the state of infinite separation. That state is assumed to have zero potential energy by convention.
- ^ Science Program's “ The Nature of the Universe " states Clarence Cleminshaw (1902-1985) served as Assistant Director of Griffith Observatory from 1938-1958 and as Director from 1958-1969. Some publications by Cleminshaw, C. H.: “Celestial Speeds,” 4 1953, equation, Kepler, orbit, comet, Saturn, Mars, velocity; Cleminshaw, C. H.: “The Coming Conjunction of Jupiter and Saturn,” 7 1960, Saturn, Jupiter, observe, conjunction; Cleminshaw, C. H.: “The Scale of The Solar System,” 7 1959, Solar system, scale, Jupiter, sun, size, light.
- ^ See. I. Bernard Cohen's Scientific American article.
- ^ Brush, Stephen G. Editor: Maxwell on Saturn's Rings, MIT Press, 1983.
- ^ See J. Bronowski and Bruce Mazlish's The Western Intellectual Tradition, Dorset Press, 1986, for a discussion of this apparent lack of understanding by Newton. Also see Truesdell's Essays in the History of Mechanics for additional background about Newton accomplishments or lack therein.
- ^ As Hufbauer points out, Newton miscalculated and published unfortunately the wrong value for the Sun's mass twice before he got it correct in his third attempt.
- ^ See Leimanis and Minorsky's historical comments.
- ^ See Moulton's Restricted Three-body Problem 's analytical and graphical solution.
- ^ See Meirovitch's book: Chapters 11, Problems in Celestial Mechanics; 12, Problem in Spacecraft Dynamics; and Appendix A: Dyadics.
- ^ See Gallian and Wilson: this short, highly technical physics paper is a full-blown example of a real-world application of D'Alembert's 1st and 2nd Principles applied to a body floating in space just above the Earth's surface (and not in outer-space). The equations developed there are similar to some equations developed here. D'Alembert theory is: if a body is supported via three points, for a total of six-degrees of freedom -- for complete initial stability -- then, when the applied external loads are applied and equal to the reacting internal inertial loads, there are no forces on the three points (the acid test). All forces, as well as couples, are balanced this way: external = internal.
- ^ See Koreneu: Supplement, Elements of Tensor Algebra and Indicial Notation in Mechanics.
- ^ See also Eisele and Mason's Applied Matrix and Tensor Analysis.
- ^ MYKLESTAD, Nils O.: Cartesian Tensors, D. Van Nostran Co., 1967.
- ^ These equations, except for the addition of the A 's, can be found in any standard dynamics textbooks; say one such as Meriam's Engineering Mechanics, Volume 2.
- ^ . The above equations are computerized (see below) so if the first, front, upper quadrant is used there is no need to worry about getting the signs wrong. If however, the left, front, upper quadrant is used too, then a matrix, and its inverse for an unsymmetrical matrix from the right side is not equal to its inverse and transpose on the left side. Each side must be independently developed. This problem occurs when an aircraft's main coordinate system splits it from the nose -to- the tail down the middle. Be careful using the n-body solution in other quadrants. Ref Gelman.
- ^ It is further conjectured some resultant internal moment sets M ζi may not be unique owning to the excessive number of reaction or support points associated with some n-body problems -- is a subject beyond the present paper's scope and development; nevertheless, if true it is contrary to classical mechanics.
- ^ See: " Gravitational N-Body Simulation " written by Henley Quadling, is a simple 16 bit DOS application allowing simulation of N particles subject to their gravitational interactions. Initial velocities and different masses can be entered. There is no limit on the number of particles; the only limiting factor is the speed of your cpu. The calculation is a true 3D calculation and there is an explicit accuracy control.
- ^ N-Body/Particle Simulation Methods
- ^ Xia, Zhihong (1992). "The Existence of Noncollision Singularities in Newtonian Systems". Annals Math. 135 (3): 411–468. JSTOR 2946572.
References
- Heggie, Douglas and Hut, Piet: The Gravitational Million-Body Problem, A Multidisciplinary Approach to Star Cluster Dynamics, Cambridge University Press, 357 pages, 2003.
- Leimanis, E., and Minorsky, N.: Dynamics and Nonlinear Mechanics, Part I: Some Recent Advances in the Dynamics of Rigid Bodies and Celestial Mechanics (Leimanis), Part II: The Theory of Oscillations (Minorsky), John Wiley & Sons, Inc., 1958.
- Moulton, Forest Ray: An Introduction to Celestial Mechanics, Dover, 1970.
- Meirovitch, Leonard: Methods of Analytical Dynamics, McGraw-Hill Book Co., 1970.
- Brouwer, Dirk and Clemence, Gerald M.: Methods of Celestial Mechanics, Academic Press, 1961.
- Cohen, Bernard I.: "Newton's Discovery of Gravity," Scientific American, pp. 167–179, Vol. 244, No. 3, Mar. 1980.
- Cohen, Bernard I.: The Birth of a New Physics, Revised and Updated, W.W. Norton & Co., 1985.
- Science Program's “ The Nature of the Universe ,” booklet, published by Nelson Doubleday, Inc., in 1968:
- Bate, Roger R.; Mueller, Donald D.; and White, Jerry: Fundamentals of Astrodynamics, Dover, 1971.
- Batin, Richard H.: An Introduction to The Mathematics and Methods of Astrodynamics, AIAA, 1987.
- Gallian, Dave A. and Wilson, Henry E.: "The Integration of NASTRAN Into Helicopter Airframe Design/Analysis," American Helicopter Society Pub., Reprint No. 780, May 1973. This is a paper, not a book.
- Lindsay, Robert Bruce: Physical Mechanics, 3rd Ed., D. Van Nostrand Co., Inc., 1961.
- Gelman, Harry: Part I: The second orthogonality conditions in the theory of proper and improper rotations: Derivation of the conditions and of their main consequences, J. Res. NBS 72B (Math. Sci.)No. 3, 1968. Part II: The intrinsic vector; Part III: The Conjugacy Theorem,J. Res. NBS 72B (Math. Sci.) No. 2, 1969. A Note on the time dependence of the effective axis and angle of a rotation, J. Res. NBS 72B (Math. Sci.)No. 3&4, Oct. 1971. These papers are on the Internet.
- Meriam, J. L.: Engineering Mechanics, Volume 1 Statics, Volume 2 Dynamics, John Wiley & Sons, 1978.
- Quadling, Henley: "Gravitational N-Body Simulation: 16 bit DOS version," June 1994. nbody*.zip is available at the http://www.ftp.cica.indiana.edu: see external links.
- Korenev, G. V.: The Mechanics of Guided Bodies, CRC Press, 1967.
- Eisele, John A. and Mason, Robert M.: Applied Matrix and Tensor Analysis, John Wiley & Sons, 1970.
- Murray, Carl D. and Dermott, Stanley F.: Solar System Dynamics, Cambridge University Press, 606 pages, 2000.
- Bronowski, Jacob and Mazlish, Bruce: The Western Intellectual Tradition, from Leonardo to Hegel, Dorsey Press. 1986.
- Truesdell, Clifford: Essays in the History of Mechanics, Springer-Verlay, 1968.
- Hufbauer, Dr. Karl C. (History of Science): Exploring the Sun, Solar Science since Galileo, The Johns Hopkins University Press, 1991. This book was sponsored by the NASA History Office.
- Diacu, F.: The solution of the n-body problem, The Mathematical Intelligencer,1996,18,p. 66–70
- Mittag-Leffler, G.: The n-body problem (Price Announcement), Acta Matematica, 1885/1886,7
- Saari, D.: A visit to the Newtonian n-body problem via Elementary Complex Variables, American Mathematical Monthly, 1990, 89, 105–119
- Newton, I.: Philosophiae Naturalis Principia Mathematica, London, 1687: also English translation of 3rd (1726) edition by I. Bernard Cohen and Anne Whitman (Berkeley, CA, 1999).
- Wang, Qiudong (1991). "The global solution of the n-body problem". Celestial Mechanics and Dynamical Astronomy. 50 (1): 73–88. Bibcode:1991CeMDA..50...73W. ISSN 0923-2958.
{{cite journal}}
:|access-date=
requires|url=
(help) - Sundman, K. F.: Memoire sur le probleme de trois corps, Acta Mathematica 36 (1912): 105–179.
- Tisserand, F-F.: Mecanique Celeste, tome III (Paris, 1894), ch.III, at p. 27.
- Hagihara, Y: Celestial Mechanics. (Vol I and Vol II pt 1 and Vol II pt 2.) MIT Press, 1970.
- Boccaletti, D. and Pucacco, G.: Theory of Orbits (two volumes). Springer-Verlag, 1998.
- Havel, Karel. N-Body Gravitational Problem: Unrestricted Solution (ISBN 978-09689120-5-8). Brampton: Grevyt Press, 2008. http://www.grevytpress.com
- Saari, D. G.; Hulkower, N. D. (1981). "On the Manifolds of Total Collapse Orbits and of Completely Parabolic Orbits for the n-Body Problem". Journal of Differential Equations. 41 (1): 27–43. doi:10.1016/0022-0396(81)90051-6.
External links
- Three-Body Problem at Scholarpedia
- More detailed information on the three-body problem
- Regular Keplerian motions in classical many-body systems
- Applet demonstrating chaos in restricted three-body problem
- Applets demonstrating many different three-body motions
- On the integration of the n-body equations
- Java applet simulating Solar System
- Java applet simulating a ring of bodies orbiting a large central mass
- Java applet simulating dust in the Solar System
- Java applet simulating a stable solution to the equi-mass 3-body problem
- Java applet simulating choreographies and other interesting n-body solutions
- A java applet to simulate the 3-d movement of set of particles under gravitational interaction
- Javascript Simulation of our Solar System
- The Lagrange Points - with links to the original papers of Euler and Lagrange, and to translations, with discussion
- http://ftp.math.utah.edu/pub/tex/bib/toc/sciam1980.html#242%281%29:January:1980
- http://www.ftp.cica.indiana.edu