Numerical integration of periodic orbits in the main problem of artificial satellite theory     

R. A. Broucke 《Celestial Mechanics and Dynamical Astronomy》1994,58(2):99-123

We describe a collection of results obtained by numerical integration of orbits in the main problem of artificial satellite theory (theJ ₂ problem). The periodic orbits have been classified according to their stability and the Poincaré surfaces of section computed for different values ofJ ₂ andH (whereH is thez-component of angular momentum). The problem was scaled down to a fixed value (–1/2) of the energy constant. It is found that the pseudo-circular periodic solution plays a fundamental role. They are the equivalent of the Poincaré first-kind solutions in the three-body problem. The integration of the variational equations shows that these pseudo-circular solutions are stable, except in a very narrow band near the critical inclincation. This results in a sequence of bifurcations near the critical inclination, refining therefore some known results on the critical inclination, for instance by Izsak (1963), Jupp (1975, 1980) and Cushman (1983). We also verify that the double pitchfork bifurcation around the critical inclination exists for large values ofJ ₂, as large as |J ₂|=0.2. Other secondary (higher-order) bifurcations are also described. The equations of motion were integrated in rotating meridian coordinates.  相似文献   

Quasi-critical orbits for artificial lunar satellites     

S. Tzirti  K. Tsiganis  H. Varvoglis 《Celestial Mechanics and Dynamical Astronomy》2009,104(3):227-239

We study the problem of critical inclination orbits for artificial lunar satellites, when in the lunar potential we include, besides the Keplerian term, the J ₂ and C ₂₂ terms and lunar rotation. We show that, at the fixed points of the 1-D averaged Hamiltonian, the inclination and the argument of pericenter do not remain both constant at the same time, as is the case when only the J ₂ term is taken into account. Instead, there exist quasi-critical solutions, for which the argument of pericenter librates around a constant value. These solutions are represented by smooth curves in phase space, which determine the dependence of the quasi-critical inclination on the initial nodal phase. The amplitude of libration of both argument of pericenter and inclination would be quite large for a non-rotating Moon, but is reduced to <0°.1 for both quantities, when a uniform rotation of the Moon is taken into account. The values of J ₂, C ₂₂ and the rotation rate strongly affect the quasi-critical inclination and the libration amplitude of the argument of pericenter. Examples for other celestial bodies are given, showing the dependence of the results on J ₂, C ₂₂ and rotation rate.  相似文献   

Symmetric and asymmetric librations in extrasolar planetary systems: a global view     

John D. Hadjidemetriou 《Celestial Mechanics and Dynamical Astronomy》2006,95(1-4):225-244

We present a global view of the resonant structure of the phase space of a planetary system with two planets, moving in the same plane, as obtained from the set of the families of periodic orbits. An important tool to understand the topology of the phase space is to determine the position and the stability character of the families of periodic orbits. The region of the phase space close to a stable periodic orbit corresponds to stable, quasi periodic librations. In these regions it is possible for an extrasolar planetary system to exist, or to be trapped following a migration process due to dissipative forces. The mean motion resonances are associated with periodic orbits in a rotating frame, which means that the relative configuration is repeated in space. We start the study with the family of symmetric periodic orbits with nearly circular orbits of the two planets. Along this family the ratio of the periods of the two planets varies, and passes through rational values, which correspond to resonances. At these resonant points we have bifurcations of families of resonant elliptic periodic orbits. There are three topologically different resonances: (1) the resonances (n + 1):n, (2:1, 3:2, ...), (2) the resonances (2n + 1):(2n-1), (3:1, 5:3, ...) and (3) all other resonances. The topology at each one of the above three types of resonances is studied, for different values of the sum and of the ratio of the planetary masses. Both symmetric and asymmetric resonant elliptic periodic orbits exist. In general, the symmetric elliptic families bifurcate from the circular family, and the asymmetric elliptic families bifurcate from the symmetric elliptic families. The results are compared with the position of some observed extrasolar planetary systems. In some cases (e.g., Gliese 876) the observed system lies, with a very good accuracy, on the stable part of a family of resonant periodic orbits.  相似文献   

Frozen orbits for satellites close to an Earth-like planet     

Shannon L. Coffey  André Deprit  Etienne Deprit 《Celestial Mechanics and Dynamical Astronomy》1994,59(1):37-72

We say that a planet is Earth-like if the coefficient of the second order zonal harmonic dominates all other coefficients in the gravity field. This paper concerns the zonal problem for satellites around an Earth-like planet, all other perturbations excluded. The potential contains all zonal coefficientsJ ₂ throughJ ₉. The model problem is averaged over the mean anomaly by a Lie transformation to the second order; we produce the resulting Hamiltonian as a Fourier series in the argument of perigee whose coefficients are algebraic functions of the eccentricity — not truncated power series. We then proceed to a global exploration of the equilibria in the averaged problem. These singularities which aerospace engineers know by the name of frozen orbits are located by solving the equilibria equations in two ways, (1) analytically in the neighborhood of either the zero eccentricity or the critical inclination, and (2) numerically by a Newton-Raphson iteration applied to an approximate position read from the color map of the phase flow. The analytical solutions we supply in full to assist space engineers in designing survey missions. We pay special attention to the manner in which additional zonal coefficients affect the evolution of bifurcations we had traced earlier in the main problem (J ₂ only). In particular, we examine the manner in which the odd zonalJ ₃ breaks the discrete symmetry inherent to the even zonal problem. In the even case, we find that Vinti's problem (J ₄+J ₂ ² =0) presents a degeneracy in the form of non-isolated equilibria; we surmise that the degeneracy is a reflection of the fact that Vinti's problem is separable. By numerical continuation we have discovered three families of frozen orbits in the full zonal problem under consideration; (1) a family of stable equilibria starting from the equatorial plane and tending to the critical inclination; (2) an unstable family arising from the bifurcation at the critical inclination; (3) a stable family also arising from that bifurcation and terminating with a polar orbit. Except in the neighborhood of the critical inclination, orbits in the stable families have very small eccentricities, and are thus well suited for survey missions.  相似文献   

The manifolds of families of 3D periodic orbits associated to Sitnikov motions in the restricted three-body problem     

E. A. Perdios 《Celestial Mechanics and Dynamical Astronomy》2007,99(2):85-104

This paper deals with the Sitnikov family of straight-line motions of the circular restricted three-body problem, viewed as generator of families of three-dimensional periodic orbits. We study the linear stability of the family, determine several new critical orbits at which families of three dimensional periodic orbits of the same or double period bifurcate and present an extensive numerical exploration of the bifurcating families. In the case of the same period bifurcations, 44 families are determined. All these families are computed for equal as well as for nearly equal primaries (μ = 0.5, μ = 0.4995). Some of the bifurcating families are determined for all values of the mass parameter μ for which they exist. Examples of families of three dimensional periodic orbits bifurcating from the Sitnikov family at double period bifurcations are also given. These are the only families of three-dimensional periodic orbits presented in the paper which do not terminate with coplanar orbits and some of them contain stable parts. By contrast, all families bifurcating at single-period bifurcations consist entirely of unstable orbits and terminate with coplanar orbits.  相似文献   

Periodic orbits and bifurcations in the Sitnikov four-body problem     

P. S. Soulis  K. E. Papadakis  T. Bountis 《Celestial Mechanics and Dynamical Astronomy》2008,100(4):251-266

We study the existence, linear stability and bifurcations of what we call the Sitnikov family of straight line periodic orbits in the case of the restricted four-body problem, where the three equal mass primary bodies are rotating on a circle and the fourth (small body) is moving in the direction vertical to the center mass of the other three. In contrast to the restricted three-body Sitnikov problem, where the Sitnikov family has infinitely many stability intervals (hence infinitely many Sitnikov critical orbits), as the “family parameter” ż₀ varies within a finite interval (while z ₀ tends to infinity), in the four-body problem this family has only one stability interval and only twelve 3-dimensional (3D) families of symmetric periodic orbits exist which bifurcate from twelve corresponding critical Sitnikov periodic orbits. We also calculate the evolution of the characteristic curves of these 3D branch-families and determine their stability. More importantly, we study the phase space dynamics in the vicinity of these orbits in two ways: First, we use the SALI index to investigate the extent of bounded motion of the small particle off the z-axis along its interval of stable Sitnikov orbits, and secondly, through suitably chosen Poincaré maps, we chart the motion near one of the 3D families of plane-symmetric periodic orbits. Our study reveals in both cases a fascinating structure of ordered motion surrounded by “sticky” and chaotic orbits as well as orbits which rapidly escape to infinity.  相似文献   

Distance in the space of energetically bounded Keplerian orbits     

Jared M. Maruskin 《Celestial Mechanics and Dynamical Astronomy》2010,108(3):265-274

In this paper we derive an explicit, analytic formula for the geodesic distance between two points in the space of bounded Keplerian orbits in a particular topology. The specific topology we use is that of a cone passing through the direct product of two spheres. The two spheres constitute the configuration manifold for the space of bounded orbits of constant energy. We scale these spheres by a factor equal to the semi-major axis of the orbit, forming a linear cone. This five-dimensional manifold inherits a Riemannian metric, which is induced from the Euclidean metric on \mathbbR⁶{\mathbb{R}^6}, the space in which it is embedded. We derive an explicit formula for the geodesic distance between any two points in this space, each point representing a physical, gravitationally bound Keplerian orbit. Finally we derive an expression for the Riemannian metric that we used in terms of classical orbital elements, which may be thought of as local coordinates on our configuration manifold.  相似文献   

Balanced earth satellite orbits     

V. Kudielka 《Celestial Mechanics and Dynamical Astronomy》1994,60(4):455-470

An analytic model for third-body perturbations and for the second zonal harmonic of the central body's gravitational field is presented. A simplified version of this model applied to the Earth-Moon-Sun system indicates the existence of high-altitude and highly-inclined orbits with their apsides in the equator plane, for which the apsidal as well as the nodal motion ceases. For special positions of the node, secular changes of eccentricity and inclination disappear too (balanced orbits). For an ascending node at vernal equinox, the inclination of balanced orbits is 94.56°, for a node at autumnal equinox 85.44°, independent of the eccentricity of the orbit. For a node perpendicular to the equinox, there exist circular balanced orbits at 90° inclination. By slightly adjusting the initial inclination as suggested by the simplified model, orbits can be found — calculated by the full model or by different methods — that show only minor variations in eccentricity, inclination, argument of perigee, and longitude of the ascending node for 10⁵ revolutions and more. Orbits near the unstable equilibria at 94.56° and 85.44° inclination show very long periodic librations and oscillations between retrogade and prograde motion.Retired from IBM Vienna Software Development Laboratory.  相似文献   

Resonant Periodic Orbits of Trans-Neptunian Objects     

Kotoulas  Thomas A.  Hadjidemetriou  John D. 《Earth, Moon, and Planets》2002,91(2):63-93

We study two and three-dimensional resonant periodic orbits, usingthe model of the restricted three-body problem with the Sun andNeptune as primaries. The position and the stability character ofthe periodic orbits determine the structure of the phase space andthis will provide useful information on the stability and longterm evolution of trans-Neptunian objects. The circular planarmodel is used as the starting point. Families of periodic orbitsare computed at the exterior resonances 1/2, 2/3 and 3/4 withNeptune and these are used as a guide to select the energy levelsfor the computation of the Poincaré maps, so that all basicresonances are included in the study. Using the circular planarmodel as the basic model, we extend our study to more realisticmodels by considering an elliptic orbit of Neptune and introducingthe inclination of the orbit. Families of symmetric periodicorbits of the planar elliptic restricted three-body problem andthe three-dimensional problem are found. All these orbitsbifurcate from the families of periodic orbits of the planarcircular problem. The stability of all orbits is studied. Althoughthe resonant structure in the circular problem is similar for allresonances, the situation changes if the eccentricity of Neptuneor the inclination of the orbit is taken into account. All theseresults are combined to explain why in some resonances there aremany bodies and other resonances are empty.  相似文献   

Analytical investigations of quasi-circular frozen orbits in the Martian gravity field     

Xiaodong Liu  Hexi Baoyin  Xingrui Ma 《Celestial Mechanics and Dynamical Astronomy》2011,109(3):303-320

Frozen orbits are always important foci of orbit design because of their valuable characteristics that their eccentricity and argument of pericentre remain constant on average. This study investigates quasi-circular frozen orbits and examines their basic nature analytically using two different methods. First, an analytical method based on Lagrangian formulations is applied to obtain constraint conditions for Martian frozen orbits. Second, Lie transforms are employed to locate these orbits accurately, and draw the contours of the Hamiltonian to show evolutions of the equilibria. Both methods are verified by numerical integrations in an 80 × 80 Mars gravity field. The simulations demonstrate that these two analytical methods can provide accurate enough results. By comparison, the two methods are found well consistent with each other, and both discover four families of Martian frozen orbits: three families with small eccentricities and one family near the critical inclination. The results also show some valuable conclusions: for the majority of Martian frozen orbits, argument of pericentre is kept at 270° because J ₃ has the same sign as J ₂; while for a minority of ones with low altitude and low inclination, argument of pericentre can be kept at 90° because of the effect of the higher degree odd zonals; for the critical inclination cases, argument of pericentre can also be kept at 90°. It is worthwhile to note that there exist some special frozen orbits with extremely small eccentricity, which could provide much convenience for reconnaissance. Finally, the stability of Martian frozen orbits is estimated based on the trace of the monodromy matrix. The analytical investigations can provide good initial conditions for numerical correction methods in the more complex models.  相似文献   

Effect of 3rd-degree gravity harmonics and Earth perturbations on lunar artificial satellite orbits     

S. Tzirti  K. Tsiganis  H. Varvoglis 《Celestial Mechanics and Dynamical Astronomy》2010,108(4):389-404

In a previous work we studied the effects of (I) the J ₂ and C ₂₂ terms of the lunar potential and (II) the rotation of the primary on the critical inclination orbits of artificial satellites. Here, we show that, when 3rd-degree gravity harmonics are taken into account, the long-term orbital behavior and stability are strongly affected, especially for a non-rotating central body, where chaotic or collision orbits dominate the phase space. In the rotating case these phenomena are strongly weakened and the motion is mostly regular. When the averaged effect of the Earth’s perturbation is added, chaotic regions appear again for some inclination ranges. These are more important for higher values of semi-major axes. We compute the main families of periodic orbits, which are shown to emanate from the ‘frozen eccentricity’ and ‘critical inclination’ solutions of the axisymmetric problem (‘J ₂ + J ₃’). Although the geometrical properties of the orbits are not preserved, we find that the variations in e, I and g can be quite small, so that they can be of practical importance to mission planning.  相似文献   

Gyldén systems: Rotation of pericenters     

André Deprit  Bruce Miller  Carol A. Williams 《Astrophysics and Space Science》1989,159(2):239-270

A canonical transformation in phase space and a rescaling of time are proposed to reduce a Keplerian system with a time-dependent Gaussian parameter to a perturbed Keplerian system with a constant Gaussian parameter. When the time variation is slow, the perturbation through second order in the reduced problem is conservative, and, as a result, the orbital space of the averaged system is a sphere on which the phase flow causes a differential rotation representing a circulation of the line of apsides. The flow presents two isolated singularities corresponding to circular orbits travelled respectively in the direct and in the retrograde sense, and a degenerate manifold of fixed points corresponding to the collision orbits. Normalization beyond order two does not break the degeneracy. Adiabatic invariants, which are conservative functions, may be computed from the normalized Hamiltonian evaluated here to the fourth order. Nonetheless so high an approximation gives little information because the normalizing Lie transformation depends explicitly on the time through mixed secular-periodic terms. As an application, an estimate is offered for the apsidal rotation that a second order time derivative in the mass of the sun would induce on planetary orbits. This suggests an observational method for determining the latter parameter for the solar wind, but the predicted motions are too slow for the current level of observational precision.  相似文献   

Planar periodic orbits in exterior resonances with Neptune     

George Voyatzis  Thomas Kotoulas 《Planetary and Space Science》2005,53(11):1189-1199

In the framework of the planar restricted three-body problem we study a considerable number of resonances associated to the basic dynamical features of Kuiper belt and located between 30 and 48 a.u. Our study is based on the computation of resonant periodic orbits and their stability. Stable periodic orbits are surrounded by regular librations in phase space and in such domains the capture of trans-Neptunian object is possible. All the periodic orbits found are symmetric and there is an indication of the existence of asymmetric ones only in a few cases. In the present work first, second and third order resonances are under consideration. In the planar circular case we found that most of the periodic orbits are stable. The families of periodic orbits are temporarily interrupted by collisions but they continue up to relatively large values of the Jacobi constant and highly eccentric regular motion exists for all cases. In the elliptic problem and for a particular eccentricity value of the primary bodies, the periodic orbits are isolated. The corresponding families, where they belong to, bifurcate from specific periodic orbits of the circular problem and seem to continue up to the rectilinear problem. Both stable and unstable orbits are obtained for each case. In the elliptic problem, the unstable orbits found are associated with narrow chaotic domains in phase space. The evolution of the orbits, which are located in such chaotic domains, seems to be practically regular and bounded for long time intervals.  相似文献   

Secular perturbations of resonant asteroids     

Yoshihide Kozai 《Celestial Mechanics and Dynamical Astronomy》1985,36(1):47-69

Secular perturbations of asteroids are derived for mean motion resonance cases under the assumptions that the disturbing planets are moving along circular orbits on the same plane and that critical arguments are fixed at stable equilibrium points. Under these assumptions the equations of motion are reduced to those of one degree of freedom with the energy integral. Then equi-energycurves on (2g-X) plane (g and X being, respectively, the argument of perihelion and (1–e²)^1/2) are derived for given values of the two constant parameters, the semi-major axis and =(1–e²)^1/2 cos i, and the variations of the eccentricity and the inclination as functions of the argument of perihelion are graphically estimated. In fact this method is applied to numbered asteroids with commensurable mean motions to estimate the ranges of the variations of orbital elements.The same method is also applied to the Pluto-Neptune system and the results are found to agree with those of numerical integrations and show that the argument of perihelion of Pluto librates around 90°.  相似文献   

Orbits near critical inclination,including lunisolar perturbations     

Michael E. Hough 《Celestial Mechanics and Dynamical Astronomy》1981,25(2):111-136

An improved theory is presented of long period perigee motion for orbits near the critical inclinations 63.4° and 116.6°. Inclusion of lunisolar perturbations andall measured zonal harmonic coefficients from a recent Earth model are significant improvements over existing theories. Phase portraits are used to depict the interaction between eccentricity magnitude and argument of perigee. The Hamiltonian constant can be chosen as the parameter to display a family of phase plane trajectories consisting of libration, circulation, and asymptotic motion along separatrices near equilibrium points. A two parameter family of phase portraits is defined by the other two integrals, the average semimajor axis and component of angular momentum resolved along the Earth's polar axis. There are regions of the parameter space where the stability and total number of equilibria can change, or two separatrices can coalesce. These phenomena signal large qualitative changes in phase portrait topology. Numerical studies show that lunisolar perturbations control stability of equilibria for orbits with semimajor axes exceeding 1.4 Earth radii. Moreover, a theory which includes lunisolar perturbations predicts larger maximum fluctuations in eccentricity and faster oscillations near stable equilibria compared to a theory which models only the zonal harmonics.  相似文献   

Periodic orbits around an oblate spheroid     

I. Stellmacher 《Celestial Mechanics and Dynamical Astronomy》1981,23(2):145-158

The general properties of certain differential systems are used to prove the existence of periodic orbits for a particle around an oblate spheroid.In a fixed frame, there are periodic orbits only fori=0 andi near /2. Furthermore, the generating orbits are circles.In a rotating frame, there are three families of orbits: first a family of periodic orbits in the vicinity of the critical inclination; secondly a family of periodic orbits in the equatorial plane with 0<e<1; thirdly a family of periodic orbits for any value of the inclination ife=0.  相似文献   

Extension of the critical inclination     

Xiaodong Liu  Hexi Baoyin  Xingrui Ma 《Astrophysics and Space Science》2011,334(1):115-124

The critical inclination is of special interest in artificial satellite theory. The critical inclination can maintain minimal deviations of eccentricity and argument of pericentre from the initial values, and orbits at this inclination have been applied to some space missions. Most previous researches about the critical inclination were made under the assumption that the oblateness term J ₂ is dominant among the harmonic coefficients. This paper investigates the extension of the critical inclination where the concept of the critical inclination is different from that of the traditional sense. First, the study takes the case of Venus for instance, and provides some preliminary results. Then for general cases, given the values of argument of pericentre and eccentricity, the relationship between the multiplicity of the solutions for the critical inclination and the values of J ₂ and J ₄ is analyzed. Besides, when given certain values of J ₂ and J ₄, the relationship between the multiplicity of the solutions for the critical inclination and the values of semimajor axis and eccentricity is studied. The results show that for some cases, the value of the critical inclination is far away from that of the traditional sense or even has multiple solutions. The analysis in this paper could be used as starters of correction methods in the full gravity field of celestial bodies.  相似文献   

Evolution of a Keplerian disk of colliding and fragmenting particles: a kinetic model with application to the Edgeworth-Kuiper belt     

Alexander V. Krivov  Miodrag Srem?evi?  Frank Spahn 《Icarus》2005,174(1):105-134

We present a kinetic model of a disk of solid particles, orbiting a primary and experiencing inelastic collisions. In distinction to other collisional models that use a 2D (mass-semimajor axis) binning and perform a separate analysis of the velocity (eccentricity, inclination) evolution, we choose mass and orbital elements as independent variables of a phase space. The distribution function in this space contains full information on the combined mass, spatial, and velocity distributions of particles. A general kinetic equation for the distribution function is derived, valid for any set of orbital elements and for any collisional outcome, specified by a single kernel function. The first implementation of the model utilizes a 3D phase space (mass-semimajor axis-eccentricity) and involves averages over the inclination and all angular elements. We assume collisions to be destructive, simulate them with available material- and size-dependent scaling laws, and include collisional damping. A closed set of kinetic equations for a mass-semimajor axis-eccentricity distribution is written and transformation rules to usual mass and spatial distributions of the disk material are obtained. The kinetic “core” of our approach is generic. It is possible to add inclination as an additional phase space variable, to include cratering collisions and agglomeration, dynamical friction and viscous stirring, gravity of large perturbers, drag forces, and other effects into the model. As a specific application, we address the collisional evolution of the classical population in the Edgeworth-Kuiper belt (EKB). We run the model for different initial disk's masses and radial profiles and different impact strengths of objects. Our results for the size distribution, collisional timescales, and mass loss are in agreement with previous studies. In particular, collisional evolution is found to be most substantial in the inner part of the EKB, where the separation size between the survivors over EKB's age and fragments of earlier collisions lies between a few and several tens of km. The size distribution in the EKB is not a single Dohnanyi-type power law, reflecting the size dependence of the critical specific energy in both strength and gravity regimes. The net mass loss rate of an evolved disk is nearly constant and is dominated by disruption of larger objects. Finally, assuming an initially uniform distribution of orbital eccentricities, we show that an evolved disk contains more objects in orbits with intermediate eccentricities than in nearly circular or more eccentric orbits. This property holds for objects of any size and is explained in terms of collisional probabilities. The effect should modulate the eccentricity distribution shaped by dynamical mechanisms, such as resonances and truncation of perihelia by Neptune.  相似文献   

Periodic orbits and their bifurcations in a 3-D system     

G. Contopoulos  B. Barbanis 《Celestial Mechanics and Dynamical Astronomy》1994,59(3):279-300

We study some simple periodic orbits and their bifurcations in the Hamiltonian . We give the forms of the orbits, the characteristics of the main families, and some existence diagrams and stability diagrams. The existence diagram of the family 1a contains regions that are stable (S), simply unstable (U), doubly unstable (DU) and complex unstable (). In the regionsS andU there are lines of equal rotation numberm/n. Along these lines we have bifurcations of families of periodic orbits of multiplicityn. When these lines reach the boundary of the complex unstable region, they are tangent to it. Inside the region there are linesm/n, along which the orbits 1a, describedn-times, are doubly unstable; however, along these lines there are no bifurcations ofn-ple periodic orbits. The families bifurcating from 1a exist only in certain regions of the parameter space (, ). The limiting lines of these regions join at particular points representing collisions of bifurcations. These collisions of bifurcations produce a nonuniqueness of the various families of periodic orbits. The complicated structure of the various bifurcations can be understood by constructing appropriate stability diagrams.  相似文献   

Straight-line oscillations generating three-dimensional motions in the photogravitational restricted three-body problem     

E. A. Perdios  V. S. Kalantonis  C. N. Douskos 《Astrophysics and Space Science》2008,314(1-3):199-208

The Sitnikov configuration is a special case of the restricted three-body problem where the two primaries are of equal masses and the third body of a negligible mass moves along a straight line perpendicular to the orbital plane of the primaries and passes through their center of mass. It may serve as a toy model in dynamical astronomy, and can be used to study the three-dimensional orbits in more applicable cases of the classical three-body problem. The present paper concerns the straight-line oscillations of the Sitnikov family of the photogravitational circular restricted three-body problem as well as the associated families of three-dimensional periodic orbits. From the stability analysis of the Sitnikov family and by using appropriate correctors we have computed accurately 49 critical orbits at which families of 3D periodic orbits of the same period bifurcate. All these families have been computed in both cases of equal and non-equal primaries, and consist entirely of unstable orbits. They all terminate with coplanar periodic orbits. We have also found 35 critical orbits at which period doubling bifurcations occur. Several families of 3D periodic orbits bifurcating at these critical Sitnikov orbits have also been given. These families contain stable parts and close upon themselves containing no coplanar orbits.  相似文献