On the Artificial Equilibrium Points in a Generalized Restricted Problem of Three Bodies

The present article studies the stability conditions of central control artificial equilibrium generalized restricted problem of three bodies. It is generalized in the sense that here we have taken the larger primary body to be in shape of an oblate spheroid. The equilibrium points are sought by the application of the propellant for which it would just balance the gravitational forces. The launching flight of such a satellite is seen to be applicable for having arbitrary space stations for these different missions. Specialty of the result of the investigation lies in the fact that an arbitrary space station can be formed to attain any specified mission.


Introduction
An important recent publication on total sailing is a book by Mclnnes [1] which provides a large number of references on the topic.Although solar sailing has been considered as a practical means of space craft propulsion only relatively recently, but the fundamental ideas are by no means new.In 1920s the Soviet father of astronautics Konstantin Tsioslkovsky and his coworker Tsander wrote of using tremendous mirror of very thin sheets [2,3] and using the pressure of sun light to attain cosmic velocities.The concept of solar sailing was reinvented much later by Richard Darwin at the IBM Watson laboratory in New Jersey and he published the first paper on solar sail in journal "Jet Propulsion" in 1958 [4], and he was the first man to coin the word solar sailing.In the early 1970s at the Battelle Laboratories in Ohio [5] discovered a trajectory that could allow a solar sail to comet Hally.
At the initial stage an extremely large (800 m × 800 m) three axis stabilized square solar sail with four deployable booms was considered but was dropped in 1977 due to the perceived risks associated with boom.Nowadays a square solar sail configuration is seen as optimum for these smaller solar sails.
The mission using solar-electric propulsion was initially strongly supported but later on due to escalating cost the propulsion was dropped [6].The Russian com-pany NPO Energia deployed a spinning 20 m reflector in February 1993 and further refinement of the reflector mission continued.Recently European space agency (ESA) and German aerospace agency (DLR) cofounded the fabrication of a 20 m × 20 m solar sail which was successfully ground tested in December 1999 [7].After the ground test the agency funded a series of mission studies and later it is being funded in USA and other centers of cosmic research.Looking to future it is being planned to set up a mission to study the various inner and outer planets by means of establishing artificial stations in their neighborhoods.
Solar sail is a proposed form of spacecraft propulsion that takes advantage of the radiation pressure to propel a spacecraft by means of a large membrane mirror.The impact of the photons emitted by the sun on the surface of the sail and its further reflection accelerates the spacecraft.Although this acceleration produced by the solar radiation pressure is smaller than the one achieved by the traditional propulsion systems, this one is continuous and unlimited.This makes long-time missions more accessible [1].It also opens a wide range of possible mission applications that cannot be achieved by a traditional spacecraft [8].
Sun-sail and other hybrid sails have been recently proposed to explore various cosmic research programmes.
An electric sail was proposed to be designed which could be capable of guarantying the fulfillment of trajectories that would be otherwise unfeasible through conventional system.In particular, the proposal [9] was to analyze the electric sail capabilities of generating a class of displaced non Keplerian orbits, useful for the observation of Polar Regions.Secondly earth-escape strategies happen to be an important problem.With growing interest in solar sailing comes the requirement to provide a basis for future detailed planetary escape mission [10] analysis by drawing together the prior work, clarifying and explaining previously anomalies.
The Lagrangian points of the circular restricted problem of three bodies (CR3BP) are known to be the five positions in an orbital system rotating with the two massive bodies where a small object (e.g. a spacecraft) affected by their gravity and centripetal forces can be stationary relative to the two larger bodies.Now if some propellant forces such as solar sail or electrodynamics or magnetic are introduced, then we have some other or more equilibrium points [11] which are artificial equilibrium points or Lagrangian points [AEPS].The problem of describing the location of AEPS and of investigating their stability property has been studied by several authors.In this connection, the works [12,13] may be specially mentioned, which investigated the effects of the thrust due to the radiation pressure and showed that there are seven equilibrium points and also linear stability was studied.Some scholars [14] by introducing a propellant force acting radially have also shown the existence of seven and more equilibrium points and have examined the stability numerically.Really different studies regarding the effect of the solar sail and the corresponding existence of equilibrium points have been made and in this connection the mention may be made [1,8,10].The low thrust systems were studied by some scholars [15].These studies have shown the existence of infinite equilibrium points depending on the magnitude of the propulsive acceleration.However only a subset of the potentially achievable AEPS turns out to be stable and could be used by a spacecraft without the use of a suitable control system.The first proposed the use of a stationary satellite hovering above the earth poles with solar radiation pressure balancing gravitational forces [16].In 1994 only continuous surfaces of unstable AEPS were derived [10], when low thrust control acceleration was added.It was proved only in 2007 that [15] a marginal stability can be found in a region of space far enough from the second primary.Really the topology of such a subset of stable AEPS is strictly dependant on the propulsion system employed by the spacecraft.It has recently been pointed [17] that for a low propulsive acceleration, the stable AEPS are confined to a very restricted region around the classical Lagrangian points.These artificially generated equilibrium points offer the possibility of considering very interesting mission applications.For example we may take Geostorm warning mission and the Polar observer mission [8].The Geostorm is a mission concept where a modest sail is placed sunwards of the classical Earth-Sun 1 point.Then with a magnetometer we can detect the solar wind polarity and give enhanced warning of the geomagnetic storms, doubling the time of alert of conventional 1 Halo orbiter such as SOHO.The aim of the polar observer mission is to provide constant viewing of the polar region and could be useful to image the Polar Regions or carry out studies on the climate evolution in the Arctic or Antarctic zone.The present paper presents a generalization of the work [17] in the sense that the formers have taken the shape of the earth to be of a sphere, whereas different from them, we have taken the shape of the earth to be an oblate spheroid.We have followed the same process as in the referred work [17].Our process differs from that of [17][18][19][20][21][22] only having mathematical expressions and in our final results.Our work consists of five sections.It starts with a general introduction giving a background of the problem and some progress that have come to our notice and also some applications.The second deals with differential equations of motion and the propellant forces leading to the equilibrium positions, i.e.AEPS.The third section deals with the linearization of the equation of motion in Hamiltonian form.In the fourth section we find out the stability conditions restricting to linear stability alone.In the fifth section we find out minimum controlled artificial equilibrium points and we finally draw our conclusion of the investigation.

L L
The paper has been concluded with the appendices giving the calculations of the various terms occurring in the paper and in the last section a final conclusion of our studies has been given.

Equations of Motion
Using dimensionless variables and a system   2 ; , M x y the equations of motion referred to the solution [5,6] when the origin is taken at the smaller mass M 2 may be written as: 2 , 2 .
x y x ny y nx where, for M fo Y-axis and the axes are taken to be rotating along a normal to r xyne with the angular velocity n].pla The inertial kinetic energy per unit mass is given by and the total potential due to the gravitational and the control acceleration is: where & x y a a are the x and y components of the control acceleration.Thus the Equations (2.1) may be written as: Then the equilibrium point  0 0 , x y for the control acceleration will be given by where   0 0 , a   are the radial distance of the equilibrium point from the primaries.

Linearized Equations of Motion Conditions
In this section we shall find out the conditions for the linear stability of the equilibrium points.We have for the Lagrangian function as: -L T V  and the canonical variables p x & p y for Hamiltonian variables Hence the Hamiltonian H may be written as Hence the transformed Hamiltonian may be written after eliminating the constant term as where,

Linear Stability Conditions
The corresponding characteristic equation will be Hence the Linearized equations of motion may be written as where f and g may be written as   be the two oscillation frequencies, then they will be given as where g and for linear stability we must have

Minimum Control Artificial Equilibrium Points
To minimize the objective function, let us consider 2 x y J a a   and writing the value for a x & a y , we get and equating the corresponding coefficients of μ, μ 2 and the value of  corresponding to these δ 0 , δ 1 , δ 2 & δ 3 will be the maximum value of  .
Hence we have,   x y a a a where, we have taken    will be the corresponding dimensional acceleration, being taken to be the distance between the primaries.It may be verified that when l 1   and 1 , the result (5.1) & (5.2) are constant giving the control acceleration to be zero.

A 
Let us look for an analytical estimation of minimum distance from the second primary (ρ) allowing the linear stability conditions.We have We have ignored the term of the   ,since 0 From the 1 st two conditions we obtain our stability boundary limit, i.e.;     ] For earth-moon and 2.3 10 and so it may be checked as Thus the value coincides with the value attained by [Bombardelli] and the corresponding control acceleration at the boundary will be In dimensional units the minimum distance from the second primary in order to have stability will be   .
So it shows that with the presence of obliquity the amount of acceleration decreases and a lesser power of thrust will be needed.

Conclusion
We have investigated the properties of minimum control artificial equilibrium points in the planar circular restricted three-body problem while the effect of the oblateness of the bigger primary body is also taken into account.We have found the analytical expressions which characterize their location, control acceleration and stability properties.It is seen that due to the presence of oblateness in the expression all the properties are likely to be affected and the disturbance is more natural since the study of effect of oblateness is quite necessary for exhaustive study of the effect.The specialty of the presence of the oblateness lies in the fact that the amount of acceleration is less and consequently a less power of thrust is required for the mission.

Appendix I (a) Expression for a, b and c:
Let the suffix (0) denote the terms corresponding to the equilibrium position   0 0 , x y and the suffix   1 A denote the contribution due to the presence of the oblateness of the larger primary body.
   in the expression for Now are the coefficients of , ,

Appendix II
Expression for f and g : We shall now onwards ignore second and higher powers of A 1 and also multiple μA 1 and so on.Thus We have  where we have dropped the suffix (0) with x 0 , y 0 , δ 0 , 0 Neglecting the term of Appendix IV We shall ignore the terms of the order of 2  and

A f a c A
Ignoring the terms of the order of 2  and also taking into account   Similarly, we shall have the expression for 2 f g  as follows: Thus neglecting the terms as mentioned above, we have Now let us transfer the origin to the equilibrium point the new coordinates be written as   ,


,  and they are given by (Refer to Appendix Section I(a))

pression for f and g and
dimensional control acceleration will be