Next Article in Journal
An Enhanced Spectral Clustering Algorithm with S-Distance
Next Article in Special Issue
Dynamics of Hyperbolically Symmetric Fluids
Previous Article in Journal
Synchronization in Non-Mirror-Symmetrical Chirogenesis: Non-Helical π–Conjugated Polymers with Helical Polysilane Copolymers in Co-Colloids
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Periodic Solutions of Nonlinear Relative Motion Satellites

by
Ashok Kumar Pal
1,
Elbaz I. Abouelmagd
2,
Juan Luis García Guirao
3,4 and
Dariusz W. Brzeziński
5,*
1
Department of Mathematics and Statistics, Manipal University Jaipur, Jaipur 303007, Rajasthan, India
2
Celestial Mechanics and Space Dynamics Research Group (CMSDRG), Astronomy Department, National Research Institute of Astronomy and Geophysics (NRIAG), Helwan 11421, Cairo, Egypt
3
Departamento de Matemática Aplicada y Estadística, Universidad Politécnica de Cartagena, Hospital de Marina, 30203 Cartagena, Región de Murcia, Spain
4
Department of Mathematics, Faculty of Science, King Abdulaziz University, P.O. Box 80203, Jeddah 21589, Saudi Arabia
5
Institute of Applied Computer Science, Faculty of Electrical, Electronic, Computer and Control Engineering, Lodz University of Technology, ul. Stefanowskiego 18/22, 90-924 Łódź, Poland
*
Author to whom correspondence should be addressed.
Symmetry 2021, 13(4), 595; https://doi.org/10.3390/sym13040595
Submission received: 5 March 2021 / Revised: 30 March 2021 / Accepted: 1 April 2021 / Published: 2 April 2021
(This article belongs to the Special Issue Recent Advance in Astrophysics)

Abstract

:
The relative motion of an outline of the rendezvous problem has been studied by assuming that the chief satellite is in circular symmetric orbits. The legitimacy of perturbation techniques and nonlinear relative motion are investigated. The deputy satellite equations of motion with respect to the fixed references at the center of the chief satellite are nonlinear in the general case. We found the periodic solutions of the linear relative motion satellite and for the nonlinear relative motion satellite using the Lindstedt–Poincaré technique. Comparisons among the analytical solutions of linear and nonlinear motions and the obtained solution by the numerical integration of the explicit Euler method for both motions are investigated. We demonstrate that both analytical and numerical solutions of linear motion are symmetric periodic. However, the solutions of nonlinear motion obtained by the Lindstedt–Poincaré technique are periodic and the numerical solutions obtained by integration by using explicit Euler method are non-periodic. Thus, the Lindstedt–Poincaré technique is recommended for designing the periodic solutions. Furthermore, a comparison between linear and nonlinear analytical solutions of relative motion is investigated graphically.

1. Introduction

1.1. Importance of Relative Motion Satellites

Relative motion is considered one of the fundamental problems in orbital mechanics, it has different applications that arise in the literature of astrodynamics for rendezvous and formation of flying satellites. The first idea concerning this problem was explained by the author of [1] in the late 19th century, who proposed analyzing the Moon’s motion. His objective was to establish a more mathematically sound method for developing tables of lunar motion, this was at the time based on “practical astronomy rather than of mathematics” in his talk.
The first practical space of relative motion was in the field of intercept and space rendezvous throughout the late 1950s and extends up to the current time. In the intercept model, a chase satellite is compelled to move in such a way that its orbit intersects the trajectory of the chief satellite at a certain time. In the rendezvous problem the relative velocities of both vehicles must be pushed to zero at the intersection instant; thereby, a docking or berthing procedure or other activities may be carried out.
The authors in [2] analyzed the relative motion within the frame of developing a guidance algorithm for a rendezvous task with an assumption that the chief satellite moves in a circular orbit, and acts as a control center. It sends messages with relative location and velocity data to slave satellites. Based on that data, the satellites perform a docking maneuver and rendezvous by using an on-board propulsion system.
The comprehension of relative motion satellites is substantial for two significant applications of orbital mechanics:
  • The space rendezvous problem.
  • The formation flight problem.
The first problem investigates the designing of centralization and decentralization algorithms to dock a satellite successfully, for instance, the space shuttle with the International Space Station. These algorithms demand the modeling of relative motion satellites, without needing permanent communication with ground-based stations. Rendezvous orbital dynamics and control is one of the key elements of space flight technology for operating space rendezvous and docking missions. Rendezvous demands a specific congruent of the orbital velocities and position vectors of the two vehicles, permitting them to stay at a constant separation during orbital station-keeping. It may or may not be followed by docking or berthing, procedures which bring the spacecraft into physical connection and generate a link between them. The same rendezvous method could be used for a spacecraft “landing” on a natural celestial body with a weak gravitational field, for example, landing on one of the Martian moons will demand the same conveniency of orbital velocities, followed by a “descent” that shares some similarities with docking.
The second problem is related to the formation of a satellite system, where the dynamics of the system are exploited to place a number of satellites in space. Formation of satellites is a novel concept that distributes a large spacecraft function among several small satellites, which are less expensive and cooperative. Most applications to relative motion are in the formation satellites, which have a considerable significance, because the use of a large number of satellites with low impact satellites running, in a concrete style, could introduce a better outcome than a single high implementation satellite. Furthermore, these formations do not require a big cost, which means a great opportunity for space mission success with more resiliency. Satellite formations have good benefits for Earth observation of space missions, where the group distribution of low resolution devices, operating in conjunction with each other, may provide a higher overall information quality than a single high resolution device.

1.2. Mathematical Models of Relative Motion Satellites

The relative motion between two satellites is described by several mathematical models [3,4,5,6]. The “Hill–Clohessy–Wiltshire model (HCW)” is as a classical approach described by a set of linearized, time-invariant ordinary differential equations. However, this model assumes a circular chief orbit and it is valid only for small initial separations between the satellites.
An assumption of an elliptical shape of the chief satellite orbit results in time-varying differential equations. In this work such a model is represented by the Tschauner–Hempel (TH) equations. However, since the applicability of this model is also limited to small initial separations between the satellites, this investigation deploys it only for comparative purposes.
In cases where a chief satellite (target or passive vehicle) is moving in a circular orbit, the orbital relative motion can be described using HCW equations [2]. A principal assumption for the application of HCW equations is a distance of separation between chief and deputy (chase or active vehicle) satellites of less than 1 km.
The authors of [7] were the first to find a closed-form solution for linearized relative motion within the frame of elliptic orbits. The authors of [8] found a singularity in Lawden’s solution and then constructed a state-transition matrix that depends explicitly on the true and the eccentric anomaly, overcoming the singularity. We find a closed-form solution for a linearized relative motion within the frame of elliptic orbits in [7] for the first time. The author of [8] found a singularity in Lawden’s solution. Afterwards, they proposed a new state-transition matrix without a singularity by explicitly using the true and the eccentric anomaly data. They derived the state transition matrix (STM) of the (TH) equations for orbits with small eccentricity, based mainly on the solutions of the HCW equations using a perturbation theory. In the framework of improving the obtained solution, approximation solutions of the second–order relative motion equations set in spherical coordinates are constructed [9].
Several versions of the relative motion dynamics have been developed to involve the perturbations effect [10,11,12,13]. The inhomogeneous gravitational field generated by the presence of the non-sphericity of objects, such as the Earth’s oblateness, is based on an assumption that the gravitational field is a sum of zonal harmonics coefficients; the authors of [10] derived a set of relative dynamics time-variant equations under the effect of the zonal harmonic parameter. An analytical solution of relative motion subject to Lorentz-force perturbations is constructed in [14]. Within the frame of developing the analysis of relative motion using state transition, matrices for the elliptical case are found, by the authors in [15,16], to be under the perturbation effect.
Furthermore, an outstanding piece of work studied the reformulation of relative motion based on non-linear system dynamics; see [17]. Another excellent piece of work proposed to analyze relative motion dynamics, control techniques and specific applications for formation flying, deployment, station keeping, and reconfiguration; see [18,19,20,21,22].
The aim of the proposed paper is to find the Periodic solutions of nonlinear relative motion using the perturbation technique of the Lindstedt–Poincaré method, which enable us to remove the secular terms, which will be growth with increasing time and will go to infinity with increasing time to infinity. However, there are some numerical methods, which can be used to estimate such solutions; for example: The Verrlet method, the step size control algorithm and the explicit Runge–Kutta method; for details, comparisons and effective methods, see [23,24]. However, we emphasize that the numerical methods are not used to verify the analytical solutions in most cases, and vice versa. In addition, there is no full warranty for convergence of numerical solutions in most cases. Thus, we are interested to find the periodic solution of non-linear relative motion, which is more realistic in practical applications by using the perturbation technique of Lindstedt–Poincaré.

2. Formulation of Motion Model

We consider a set of two satellites model, the first is called target satellite (or chief or passive vehicle) and the second is called chaser (or deputy or active vehicle). In general, this model is implemented to accomplish some space mission such as space rendezvous. The mission between two space crafts or satellites is successful when the two satellites have the same velocity and position vectors at a certain time. In the initiates of rendezvous time the two parts in a rendezvous sequence are:
  • Phasing of maneuvers to accomplish a rendezvous in the timing sequence, which will bring the two satellites into close proximity.
  • Maneuvers of finalizing rendezvous that involve the relative motion between the two satellites for rendezvous and docking.
Now we will assume that the inertial references frame with origin, which coincides with the center of the main body, say Earth’s center; the orthogonal unit vector is denoted by { i x , i y , i z } . The vector i x in the direction of the equinoxes and the vectors i y is perpendicular to it, while i z passes through the North Pole. To describe a relative motion satellite in Local-Vertical-Local-Horizontal (LVLH) coordinates, we proposed that the orientation of these frame by the unit vector triad { i r , i θ , i h } , where the vector i r lies in the chief radial direction, i h lies in the direction of the chief angular momentum, and i θ completes the right-handed orthogonal triad, such that i θ = i h i r ; see Figure 1.
According to Figure 1, the deputy satellite relative position vector in LVLH coordinates can be written as:
R = ξ i r + ϑ i θ + ζ i h ,
where ξ , ζ , and ϑ are the components of the relative position vector along the radial, long track, and out of plane, respectively.
Since r d = r c + r , then
r d = ( r + ξ ) i r + ϑ i θ + ζ i h ,
where r c = r i r represents the current radius of the chief’s satellite orbit. The LVLH frame rotation with angular velocity ω is given by
ω = ω r i r + ω θ i θ + ω h i h ,
where ω r = ω θ = 0 and ω = ω h i h , while the attached frame rotates with the angular velocity ω h = θ ˙ = f ˙ . We have to note that θ = Ω + f in which Ω and f are the argument of perigee and the true anomaly, respectively, within the frame of two-body motion. The angular momentum magnitude of the chief satellite is expressed by
h = r 2 θ ˙ , θ ¨ = 2 r ˙ θ ˙ r ,
where r is the trajectory of the chief (target) satellite, which may be circular or elliptical. On the other hand, the expression of the chief’s true anomaly rate θ ˙ can be written in terms of a semi-latus rectum as
θ ˙ = μ p r 4 ,
wherein the elliptical case p = a ( 1 e 2 ) , a are the semi-major axes and e is the eccentricity, μ is the gravitational constant and r = p / ( 1 + e cos f ) . The angular velocity and the acceleration of the the LVLH reference frame are given by
ω = θ ˙ i h = h r 2 i h , ω ˙ = θ ¨ i h = 2 r ˙ θ ˙ r i h .
Now, the general vectorial equations that model the relative motion in the LVLH reference frame are
R ¨ + 2 ω × R ˙ + ω ˙ × R + ω × ( ω × R ) = T ,
where T is the sum of total external forces acting on the deputy and chief, respectively; ( . ) and ( . . ) denote the first and second derivative with respect to time.
Substituting Equations (1) and (2) into Equation (3), we obtain
ξ ¨ 2 θ ¨ ϑ θ ˙ ϑ ˙ θ ˙ 2 ξ i r + ϑ ¨ + θ ¨ ξ + 2 θ ˙ ξ ˙ θ ˙ 2 ϑ i θ + ζ ¨ i h = T ,
Now we can write the right-hand side of Equation (4) in the form
T = T d T c ,
where T d is the external forces acting on the deputy and T c is the external forces acting on the chief satellite, which can be written in the case of the gravitational field being generated by a spherical object, such as
T d = μ R d r d 3 , T c = μ R c r c 3 ,
Hence,
T = μ ( r + ξ ) i r + ϑ i θ + ζ i h ( r + ξ ) 2 + ϑ 2 + ζ 2 3 / 2 + μ 1 r 2 i r ,
Substitution Equation (5) into Equation (4); then, the vector form or relative motion equations will be controlled by
ξ ¨ 2 θ ¨ ϑ θ ˙ ϑ ˙ θ ˙ 2 ξ i r + ϑ ¨ + θ ¨ ξ + 2 θ ˙ ξ ˙ θ ˙ 2 ϑ i θ + ζ ¨ i h = μ ( r + ξ ) i r + ϑ i θ + ζ i h ( r + ξ ) 2 + ϑ 2 + ζ 2 3 / 2 + μ 1 r 2 i r ,
The vector Equation (6) can be degraded to
ξ ¨ 2 θ ˙ ϑ ˙ θ ˙ 2 ξ θ ¨ ϑ = μ ( r + ξ ) ( r + ξ ) 2 + ϑ 2 + ζ 2 3 / 2 + μ r 2 , ϑ ¨ + 2 θ ˙ ξ ˙ θ ˙ 2 ϑ + θ ¨ ξ = μ ϑ ( r + ξ ) 2 + ϑ 2 + ζ 2 3 / 2 , ζ ¨ = μ ζ ( r + ξ ) 2 + ϑ 2 + ζ 2 3 / 2 .
Equation (7) represents the general equations of the satellite relative motion in LVLH coordinates.

2.1. General Circular Relative Motion Satellites

Now, we denote v c and v d as the potential functions for gravitational forces T c and T d , acting on the chief and deputy satellite, respectively; hence,
v c = μ r , v d = μ r 1 + R 2 r 2 + 2 ξ r 1 / 2 ,
Note that ( ξ / r ) = ( ξ / R ) . ( R / r ) ; Equation (8) is a generating function for Legendre polynomials with argument ( ξ / R ) ; consequently,
v d = μ r n = 0 ( 1 ) n R r n P n ξ / R ,
where P n is the Legendre polynomial with degree n, where
P 0 ξ / R = 1 , P 1 ξ / R = ξ R , P 2 ξ / R = 1 2 ξ R 2 1 ,
Then,
v d = μ R 1 ξ r + ( 2 ξ 2 ϑ 2 ζ 2 ) 2 r 2 + v h ,
where v h is the higher-order of a Legendre polynomial given by
v h = μ r n = 3 ( 1 ) n R r n P n ξ / R .
Then the total external forces can be rewritten as
T = Δ v ,
where
Δ v = Δ v d Δ v c .
Since Δ v d is the gravitational acceleration acting on the deputy and Δ v c is the gravitational acceleration acting on the chief’s satellite, while Δ v is the total relative gravitational acceleration, then
Δ v = 2 μ ξ r 3 + v h ξ i r + μ ϑ r 3 + v h ϑ i θ + μ ζ r 3 + v h ζ i h ,
and
T = 2 μ ξ r 3 v h ξ i r μ ϑ r 3 + v h ϑ i θ μ ζ r 3 + v h ζ i h ,
Utilizing Equations (9) and (10) with Equation (4), the vector equation of relative motion can be rewritten as
ξ ¨ 2 θ ¨ ϑ θ ˙ ϑ ˙ θ ˙ 2 ξ i r + ϑ ¨ + θ ¨ ξ + 2 θ ˙ ξ ˙ θ ˙ 2 ϑ i θ + ζ ¨ i h = 2 μ ξ r 3 v h ξ i r + μ ϑ r 3 v h ϑ i θ + μ ζ r 3 v h ζ i h ,
Hence, Equation (11) can be written as three differential equations of second order
ξ ¨ 2 θ ˙ ϑ ˙ θ ˙ 2 + 2 μ r 3 ξ θ ¨ ϑ = v h ξ , ϑ ¨ + 2 θ ˙ ξ ˙ θ ˙ 2 μ r 3 ϑ + θ ¨ ξ = v h ϑ , ζ ¨ + μ ζ r 3 = v h ζ .
Equation (12) represents the general dynamical system of the relative motion satellite. However, we have to observe that the quantities in the right-hand sides, v h / ξ , v h / ϑ , v h / ζ , represent the components of perturbing gravitational accelerations, which are considered the main contribution for higher order of nonlinearity to a dynamical system of relative motion.

2.2. Linear Circular Relative Motion Satellites

In the case where the chief’s orbit is moving in circular orbits, then r = a and θ ˙ = n are constants. Thereby, the rotation frame with these properties is referred to as Hill’s frame and rotates with angular velocity n = μ / a 3 , and the resulting equations are referred to as HCW equations, as in [2]. In addition, if the perturbing acceleration or the higher order on nonlinearity terms ( v h = 0 ) are neglected. Then the linearized relative motion equations will be governed by
ξ ¨ 2 n ϑ ˙ 3 n 2 ξ = 0 , ϑ ¨ + 2 n ξ ˙ = 0 , ζ ¨ + n 2 ζ = 0 .
Let ξ = x , ϑ = y , ζ = z ; thereby, the above equations can be written as
x ¨ 2 n y ˙ 3 n 2 x = 0 , y ¨ + 2 n x ˙ = 0 , z ¨ + n 2 z = 0 .
The above dynamical system is famous and is called the Hill–Clohessy–Wiltshire equations, and represents the linearity of relative satellite motion. Furthermore, the dynamical system of relative motion (HCW) is invariant under the symmetry ( x , y , z ) ( x , y , z ) .

3. Solutions of HCW Equations

In order to construct the solutions of the dynamical system which represent the linear relative motion satellite, we will assume that the deputy satellite starts its motion under the following set of the initial conditions
x ( 0 ) = x 0 0 , x ˙ ( 0 ) = x ˙ 0 0 , y ( 0 ) = y 0 0 , y ˙ ( 0 ) = y ˙ 0 0 , z ( 0 ) = z 0 0 , z ˙ ( 0 ) = z ˙ 0 0 ,
Thereby, after integration of the second equation in System (13), we obtain
y ˙ + 2 n x = X 0 ,
where X 0 is a constant which can be determined from the initial conditions, so one obtains
X 0 = y ˙ 0 0 + 2 n x 0 0 ,
and Equation (15) will take the following form
y ˙ = X 0 2 n x .
Substituting y ˙ into the first equation in System (13), we get
x ¨ + n 2 x 2 X 0 n = 0 .
Let ψ = x 2 X 0 / n , then d ψ / d t = d X / d t and d 2 ψ / d t 2 = d 2 X / d t 2 and, the Equation (17) becomes
d 2 ψ d t 2 + n 2 ψ = 0 ,
Hence, Equation (18) has the solution
ψ = X sin [ n t + θ 0 ] .
Then the final form solution is
x = X sin [ n t + θ 0 ] + 2 n y ˙ 0 0 + 2 n x 0 0 .
Substituting Equation (19) into Equation (16) and integrating it, we get
y = 2 X cos [ n t + θ 0 ] 3 y ˙ 0 0 + 2 n x 0 0 t + Y 0 ,
where Y 0 the constant of the integration.
To complete our steps, we have to integrate the third equation in System (13), which represents the out-plane motion; hence, we get
z = Z sin ( n t + φ 0 ) .
Using the obtained solutions of Equations (19)–(21), the solutions of System (13) are summarized in the following formulas
x = X sin [ n t + θ 0 ] + 2 n y ˙ 0 0 + 2 n x 0 0 , y = 2 X cos [ n t + θ 0 ] 3 y ˙ 0 0 + 2 n x 0 0 t + Y 0 , z = Z sin ( n t + φ 0 ) ,
where X , Y , Z , θ 0 , and φ 0 are arbitrary constants that must match the initial conditions.
Here, we observe that the second formula in Solutions (22) contains an unbounded term, which will grow to infinity with increasing time. To get bounded solutions, in particular, periodic solutions, the coefficient of the unbounded term must equal zero; these will give an additional condition to get periodic and bounded solutions; thereby, we will obtain the new conditions on the initial conditions as
y ˙ 0 0 + 2 n x 0 0 = 0 .
So we rewrite the periodic solutions of System (13) in the following formulas
x = X sin [ n t + θ 0 ] , y = Y 0 + 2 X cos [ n t + θ 0 ] , z = Z sin ( n t + φ 0 ) .
Applying the first derivative to Equations (24), we obtain the relative velocities in the form
x ˙ = n X cos [ n t + θ 0 ] , y ˙ = 2 n X sin [ n t + θ 0 ] , z ˙ = n Z cos ( n t + φ 0 ) ,
If we apply the initial conditions on Equation (24) and the condition in Equation (23), then the constants of integration can be determined by
Y 0 = y 0 0 2 n x ˙ 0 0 , X 2 = x 0 0 2 + 1 n 2 x ˙ 0 0 2 , Z 2 = z 0 0 2 + 1 n 2 z ˙ 0 0 2 ,
and
tan θ 0 = n x 0 0 x ˙ 0 0 , tan φ 0 = n z 0 0 z ˙ 0 0 .
Equations (24) and (25) represent the relative positions and velocities of the chaser satellite for space rendezvous within the framework of linearized motion.
With the help of Equations (24)–(27), the initial amplitude of the chaser satellite for space rendezvous, i.e., ( x 0 0 , y 0 0 , z 0 0 ) is taken as ( 0.13 , 0.12 , 0.10 ) , with different values of phase angle. The projection of the periodic solution in the xy-plane is shown in Figure 2, whereas the orbit in three dimensions is shown in Figure 3 at different values of the phase angle.
It is clear from Figure 2a that the analytical solution of linear relative motion in two dimensions is periodic, and the obtained solution by a numerical integration of the explicit Euler method using step size 10 3 at a set of initial condition x 0 0 , y 0 0 , x ˙ 0 0 , y ˙ 0 0 as 0 , 1.2 , 1.0 , 0 is also periodic; this solution is shown in Figure 2b. In addition, the two solutions are symmetric due to the y-axis ( x = 0 ) .
The linear motion in three dimensions is also periodic due two both solutions. We remark that both analytical and numerical solutions of linear relative motion in three dimensions are also periodic. The analytical solution is given in Figure 3a, while the solution by a numerical integration of the explicit Euler method using step size 10 3 at a set of initial condition x 0 0 , y 0 0 , x ˙ 0 0 , y ˙ 0 0 as 0 , 1.2 , 1.0 , 0 is given in Figure 3b.

Nonlinear Circular Relative Motion Satellites

In the case of considering the perturbed potential v h 0 , and keeping the nonlinearity up to second order terms, the relative motion satellite equations of System (12) are approximated to the following equations
ξ ¨ 2 n ϑ ˙ 3 n 2 ξ = 3 n 2 2 r ϑ 2 + ζ 2 2 ξ 2 , ϑ ¨ + 2 n 2 ξ ˙ = 3 n 2 ξ ϑ r , ζ ¨ + n 2 ζ = 3 n 2 ξ ζ r .
The above dynamical system consists of the three second-order differential equations with second degree, which represent the relative motion of the deputy satellite.
Now we intend to find the periodic solutions of these equations to find the relative position and velocity of the deputy satellite to accomplish space rendezvous. The solutions of these equations will be generalized or extended for the linear solutions; then, we have to find the solutions under the previous conditions of linear relative motion. Unfortunately, there are no exact analytical solutions for these equations, and the exact solutions are not possible. Thereby, we must introduce techniques or perturbation methods, which can be used to evaluate approximate or semi-analytical solutions; in particular, we must construct periodic solutions. So, the next section will be devoted to analysis of a perturbation technique, which enables us to find an appropriate approximated analytical solution to perform a space rendezvous mission.

4. Periodic Solution of Nonlinear Relative Motion Satellites

This section is devoted to constructing the periodic solutions of the nonlinear relative motion satellite. From physical point of view, it is sufficient to find the solutions of the dynamical system within the frame of second degree or third in maximum. This view can be applied for the model of relative motion satellite, because the relative distance (distance between chief and deputy satellites) is very small compared to the initial position of the target satellite.

4.1. Legitimacy of Perturbation Techniques

The classical perturbation theory, which is called the straightforward expansion technique, the Lindstedt–poincaré method, and the multiple scales method are used when the dynamical systems are either weakly nonlinear or weakly non-autonomous; this means that the effects of either nonlinearity or non-autonomy are very small. Alternatively, the dynamical systems with these properties can be considered either linear or quasi-linear. In most cases, the systems of differential equations have linear expressions and small nonlinear or non-autonomous expressions in such a way that they are separated from each other.
In the context of either weak nonlinearity or weak non-autonomy, more approximations may still be required by finding a parameter in the dynamical system, which can be used as a very small perturbation effect, but its value will be replaced by one at the end. This technique is reasonable and accepted when the terms of the dynamical system are quasi-linear or small in themselves. The obtained parameters in this case are called “place keeping parameters”.

4.2. Legitimacy of Relative Motion Equations

In line with the reasons outlined in the previous subsection, we can neither apply the multiple scales technique nor the Lindstedt–Poincaré method to System (28). Thereby, we will use the place keeping parameters method, after carrying out the relative equations of motion in dimensionless variables.
In order to set the equations in dimensionless variables, we assume that ( ξ , ϑ , ζ ) = r ( x , y , z ) ; then, System (28) can be written as
x ¨ 2 n y ˙ 3 n 2 x = 3 2 n 2 ( y 2 + z 2 2 x 2 ) , y ¨ + 2 n x ˙ = 3 n 2 x y , z ¨ + n 2 z = 3 n 2 x z .
Since R 0 < < r , we can write ϵ = R 0 r < < 1 , where R 0 is the measure of relative orbit size. Regarding this idea, we can use the place keeping parameters method, and replace the coordinate ( x , y , z ) by ( ϵ x , ϵ y , ϵ z ) . Then System (29) can be written as
x ¨ 2 n y ˙ 3 n 2 x = 3 2 n 2 ϵ ( y 2 + z 2 2 x 2 ) , y ¨ + 2 n x ˙ = 3 ϵ n 2 x y , z ¨ + n 2 z = 3 ϵ n 2 x z .
The above equations represent the nonlinear relative motion satellite in dimensionless variables, which are convenient for applying some perturbation techniques, such as multiple scales technique [25,26,27] and KBM methods or Lindstedt–Poincaré [28,29,30,31]. However, we will examine the provided solution by the latter method with the initial conditions of the linear relative motion satellite; these conditions are describe in Equation (14).

5. Periodic Solutions by the Lindstedt–Poincar é Technique

Now we will apply the Lindstedt–Poincaré technique on the relative satellite motion to find the periodic solution of the equations of motion of the second degree, say System (30). For simplicity, we suppose that θ = n t ; then, System (30) can be written in the form
d 2 x d θ 2 2 d y d θ 3 x = 3 2 ϵ ( y 2 + z 2 2 x 2 ) , d 2 y d θ 2 + 2 d x d θ = 3 ϵ x y , d 2 z d θ 2 + z = 3 ϵ x z .
To account for the nonlinear dependence of the frequency, we will introduce the stretch variable, τ = ω θ ; we denote d Γ / d τ = Γ and d 2 Γ / d τ 2 = Γ ; then, System (31) can be written in the form
ω 2 x 2 ω y 3 x = 3 2 ϵ ( y 2 + z 2 2 x 2 ) , ω 2 y + 2 ω x = 3 ϵ x y , ω 2 z + z = 3 ϵ x z .
Now we look for approximated solutions in the form of a power series as
x ( τ , ϵ ) = n = 0 ϵ n x n ( τ ) , y ( τ , ϵ ) = n = 0 ϵ n y n ( τ ) , z ( τ , ϵ ) = n = 0 ϵ n z n ( τ ) ,
and
ω ( ϵ ) = n = 0 ϵ n ω n , ω 0 = 1 ,
The initial conditions of position components can be rewritten in the form
x ( 0 , ϵ ) = n = 0 ϵ n x n ( 0 ) = x 0 0 , y ( 0 , ϵ ) = n = 0 ϵ n y n ( 0 ) = y 0 0 , z ( 0 , ϵ ) = n = 0 ϵ n z n ( 0 ) = z 0 0 ,
while the initial conditions of velocity components can also be rewritten in the form
x ( 0 , ϵ ) = n = 0 ϵ n x n ( 0 ) = x 0 0 , y ( 0 , ϵ ) = n = 0 ϵ n y n ( 0 ) = y 0 0 , z ( 0 , ϵ ) = n = 0 ϵ n z n ( 0 ) = z 0 0 .
Substituting Equation (33) into Equation (32) with the help of Equation (34) and keeping the terms up to first order in ϵ , we obtain
x 0 2 y 0 3 x 0 + ϵ x 1 2 y 1 3 x 1 + 2 ω 1 x 0 2 ω 1 y 0 3 2 ( y 0 2 + z 0 2 2 x 0 2 ) = 0 , y 0 + 2 x 0 + ϵ y 1 + 2 x 1 + 2 ω 1 y 0 + 2 ω 1 x 0 3 x 0 y 0 = 0 , z 0 + z 0 + ϵ z 1 + z 1 + 2 ϵ ω 1 Z 0 3 x 0 z 0 = 0
where the frequency is
ω ( ϵ ) = 1 + ϵ ω 1 ,
The set of initial conditions in Equations (35) and (36) may be rewritten as
x ( 0 , ϵ ) = x 0 ( 0 ) + ϵ x 1 ( 0 ) = x 0 0 , y ( 0 , ϵ ) = y 0 ( 0 ) + ϵ y 1 ( 0 ) = y 0 0 , z ( 0 , ϵ ) = z 0 ( 0 ) + ϵ z 1 ( 0 ) = z 0 0 ,
and
x ( 0 , ϵ ) = x 0 ( 0 ) + ϵ x 1 ( 0 ) = x 0 0 , y ( 0 , ϵ ) = y 0 ( 0 ) + ϵ y 1 ( 0 ) = y 0 0 , z ( 0 , ϵ ) = z 0 ( 0 ) + ϵ z 1 ( 0 ) = z 0 0 ,
Equating each of the coefficients of ϵ that have the same power in System (37), using Conditions (38) and (39), we get the unperturbed equations system in the form
x 0 2 y 0 3 x 0 = 0 , y 0 + 2 x 0 = 0 , z 0 + z 0 = 0 ,
and the set of initial conditions is
x 0 ( 0 ) = x 0 0 , x 0 ( 0 ) = x ˙ 0 0 , y 0 ( 0 ) = y 0 0 , y 0 ( 0 ) = y ˙ 0 0 , z 0 ( 0 ) = z 0 0 , z 0 ( 0 ) = z ˙ 0 0 ,
while the equations which represent the first corrections are
x 1 2 y 1 3 x 1 = 3 2 y 0 2 + z 0 2 2 x 0 2 2 ω 1 x 0 + 2 ω 1 y 0 , y 1 + 2 x 1 = 3 x 0 y 0 2 ω 1 y 0 2 ω 1 x 0 , z 1 + z 1 = 3 x 0 z 0 2 ω 1 z 0 ,
with the initial conditions
x 1 ( 0 ) = 0 , x 1 ( 0 ) = 0 , y 1 ( 0 ) = 0 , y 1 ( 0 ) = 0 , z 1 ( 0 ) = 0 , z 1 ( 0 ) = 0 ,
To accomplish our procedure and find the approximated periodic solutions, we have to find solutions for both the dynamical systems in Equations (40) and (42)
Zero order solutions
According to the solutions of linear motion, the general periodic solutions of the unperturbed relative motion will be controlled by
x 0 = X sin [ τ + θ 0 ] , y 0 = Y 0 + 2 X cos [ τ + θ 0 ] , z 0 = Z sin [ τ + φ 0 ] ,
With a help of Conditions (41), the constants of integration are given by
Y 0 = y 0 0 2 x ˙ 0 0 , X 2 = x 0 0 2 + x ˙ 0 0 2 , Z 2 = z 0 0 2 + z ˙ 0 0 2 ,
and
tan θ 0 = x 0 0 x ˙ 0 0 , tan φ 0 = z 0 0 z ˙ 0 0 .
First order solutions
Substituting Equation (44) into Equation (42), we obtain
x 1 2 y 1 3 x 1 = X 11 + X 12 cos ( τ + θ 0 ) + X 13 sin ( τ + θ 0 ) + X 14 sin ( τ + φ 0 ) + X 15 cos 2 ( τ + θ 0 ) + X 16 cos 2 ( τ + φ 0 ) , y 1 + 2 x 1 = Y 11 sin ( τ + θ 0 ) + Y 12 cos ( τ + θ 0 ) + Y 13 sin 2 ( τ + θ 0 ) , z 1 + z 1 = Z 11 + Z 12 sin ( τ + φ 0 ) + Z 13 cos ( 2 τ + θ 0 + φ 0 ) ,
where X 1 i , i = 1 , 2 , 3 , 6 , Y i j and Z 1 j , j = 1 , 2 , 3 are constants, which are calculated by the following relations
X 11 = 3 4 2 X 2 + 2 Y 0 2 + Z 2 , X 12 = 6 X Y 0 , X 13 = 4 ω 1 X , X 14 = 4 ω 1 Z , X 15 = 9 2 X 2 , X 16 = 3 4 Z 2 , Y 11 = 3 X Y 0 , Y 12 = 2 ω 1 X , Y 13 = 3 X 2 , Z 11 = 3 2 X Z cos ( θ 0 φ 0 ) , Z 12 = 2 ω 1 Z , Z 13 = 3 2 X Z .
Equation (45) consists of a dynamical system of second-order non-homogeneous differential equations. After integrating this system with the help of Conditions (43), we observe that the terms with coefficients X 1 i = 1 , 2 , 3 , 4 , Y 11 , Y 12 , and Z 12 represent secular terms, as well as the extra term with coefficient S = 2 X 2 + Z 2 + 3 X 2 cos 2 θ 0 . Therefore, if these coefficients are nonzero, then the correction solutions will involve secular terms in the scale τ . However, that is exactly what we need to avoid. Hence we will remove these terms by equating their coefficients by zeros. These procedures ensure that the solutions will not include secular terms and at least no secularities appear in the first two terms in the perturbation series. Furthermore, the terms with coefficients X 1 i = 1 , 2 , 3 , 4 , Y 11 , Y 12 , and Z 12 represent solutions for the associated homogeneity of Equation (45); hence, we set the coefficients to zeros before integration.
In this context, the removal of the secular terms yields
ω 1 = 0 , Y 0 = 0 , S = 0 .
Then the first correction will be controlled by
x 1 ( τ ) = 1 2 3 2 X 2 + Z 2 + 2 X 2 7 cos 2 θ 0 + 2 cos ( 2 θ + τ ) Z 2 cos 2 ϕ 0 + 2 cos ( τ + ϕ 0 ) sin 2 1 2 τ , y 1 ( τ ) = 1 4 9 X 2 sin 2 θ 0 14 X 2 sin ( 2 θ 0 τ ) + Z 2 8 cos 1 2 ( τ + 4 ϕ 0 ) sin 3 1 2 τ + 6 sin τ + X 2 12 sin τ sin 2 ( θ 0 + τ ) + 6 sin ( 2 θ 0 + τ ) , z 1 ( τ ) = X Z 3 cos ( θ 0 ϕ 0 ) cos ( θ 0 + ϕ 0 ) 2 cos ( τ + θ 0 ϕ 0 ) sin 2 1 2 τ .
Therefore, the periodic solutions of the nonlinear relative motion satellite can be written in the form
x ( τ , ϵ ) = x 0 ( τ ) + ϵ x 1 ( τ ) , y ( τ , ϵ ) = y 0 ( τ ) + ϵ y 1 ( τ ) , z ( τ , ϵ ) = z 0 ( τ ) + ϵ z 1 ( τ ) ,
where τ = n ω t .
Since we use the method of place keeping parameters, we have to take ϵ = 1 ; then, Solutions (47) will be rewritten as
x ( τ , ϵ ) = x 0 ( τ ) + x 1 ( τ ) , y ( τ , ϵ ) = y 0 ( τ ) + y 1 ( τ ) , z ( τ , ϵ ) = z 0 ( τ ) + z 1 ( τ ) .
Substituting Equations (44) and (46) into Equation (48), the periodic solutions of the nonlinear relative motion satellite, with the first corrections, can be estimated approximately by using the Lindstedt–Poincaré technique.
With the help of Equation (48), the initial amplitude of the chaser satellite for space rendezvous, i.e., ( x 0 0 , y 0 0 , z 0 0 ) , works similar to in the previous section ( 0.13 , 0.12 , 0.10 ) with different values for the phase angle. The projection of the periodic solution in the xy-plane is shown in Figure 4, whereas the orbit in three dimensions is shown in Figure 5 with different values of phase angle.
A comparison between the nonlinear relative motion solutions of the Lindstedt–Poincaré technique obtained by the numerical integration of the explicit Euler method are shown in Figure 4, in two dimensions, and in Figure 5, in three-dimensional motion. While numerical solutions for linear motion are periodic, the numericality of nonlinear motion is not periodic; this can be observed in Figure 4b and Figure 5b.
The projections in the xy-plane of the periodic solution of linear and nonlinear motion are shown in Figure 6 at different phases, whereas the differences between orbits in three dimensions of the linear and nonlinear motion at the same different phases angle with φ 0 = π / 4 are shown in Figure 7.
The graphical investigations in both Figure 6 and Figure 7 show that the analytical solutions of linear and nonlinear motions are periodic; however, there are differences between the style of the two motions. The periodicity of the nonlinear relative motion solutions gives a considerable verification due to the solution obtained with the Lindstedt–Poincaré technique.

6. Conclusions

In this paper, we studied the relative motion of an outline of the space rendezvous problem, assuming that the chief satellite was in circular symmetric orbits. It was found that the equations that describe the motion of the deputy satellite with respect to the references fixed at the center of the chief satellite is nonlinear in the general case. Then the periodic solutions first for linear relative motion satellite were evaluated. Here a key role is the symmetry of the periodic solution joined with the fact that the system HCW is invariant under negative symmetry.
The Lindstedt–Poincaré technique was used to get the approximated periodic solutions for nonlinear differential equations of the second degree. This method is referred to by many sources as a method of successive approximations. In our study, we use it to find the first corrections of a relative motion satellite. In this context, the obtained solutions are more general compared to the linear solutions. The analysis of a linear relative motion satellite and its periodic solutions are constructed. The approximate periodic solutions of the nonlinear relative motion satellite are constructed using the Lindstedt–Poincaré technique.
Comparisons among the analytical solution of linear motion and the obtained solution by the numerical integration of the explicit Euler method are investigated, and we show that bth solutions are periodic in two and three dimensions. Moreover, comparisons among the obtained solutions of nonlinear relative motion using the Lindstedt–Poincaré technique and the numerical solution by integration of the explicit Euler method were studied. However, the numerical solutions are not periodic in both motions in two and three dimensions. Thus, the Lindstedt–Poincaré technique is recommended for designing the periodic solutions.

Author Contributions

Conceptualization: E.I.A.; methodology: E.I.A., D.W.B., J.L.G.G. and A.K.P.; software: D.W.B.; validation: E.I.A., D.W.B., J.L.G.G. and A.K.P.; formal analysis: J.L.G.G.; investigation: E.I.A. and A.K.P.; resources: A.K.P.; writing—original draft preparation: E.I.A., D.W.B. and A.K.P.; writing—review and editing: E.I.A., D.W.B., J.L.G.G. and A.K.P.; visualization: D.W.B.; supervision: E.I.A.; project administration: J.L.G.G.; funding acquisition: D.W.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research is co-financed by the Institute of Applied Computer Science, Lodz University of Technology as a part of its research activity nr 501/2-24-1-2.

Data Availability Statement

Not applicable.

Acknowledgments

Publication of this research is financed by the Institute of Applied Computer Science, Lodz University of Technology, Poland as a part of its statutory activity nr 501/2-24-1-2.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hill, G.W. Researches in the lunar theory. Am. J. Math. 1878, 1, 5–26. [Google Scholar] [CrossRef]
  2. Clohessy, W.; Wiltshire, R. Terminal guidance system for satellite rendezvous. J. Aerosp. Sci. 1960, 27, 653–658. [Google Scholar] [CrossRef]
  3. Sengupta, P.; Vadali, S.R.; Alfriend, K.T. Modeling and control of satellite formations in high eccentricity orbits. J. Astronaut. Sci. 2004, 52, 149–168. [Google Scholar] [CrossRef]
  4. Lee, W.-K.; Yoo, S.-M.; Park, S.-Y.; Choi, K.-H.; Chang, Y.-K. Determination of initial conditions for satellite formation elying in elliptical orbits. J. Astron. Space Sci. 2005, 22, 21–34. [Google Scholar] [CrossRef] [Green Version]
  5. Park, H.-E.; Park, S.-Y.; Lee, S.-J.; Choi, K.-H. Analysis of linear and nonlinear relative orbit dynamics for satellite formation flying. J. Astron. Space Sci. 2009, 26, 317–328. [Google Scholar] [CrossRef]
  6. Leomanni, M.; Garulli, A.; Giannitrapani, A.; Quartullo, R. Satellite relative motion modeling and estimation via nodal elements. J. Guid. Control Dyn. 2020, 43, 1904–1914. [Google Scholar] [CrossRef]
  7. Lawden, D.F. Optimal Trajectories for Space Navigation; Butterworths: London, UK, 1963; Volume 3. [Google Scholar]
  8. Carter, T.E. State transition matrices for terminal rendezvous studies: Brief survey and new example. J. Guid. Control Dyn. 1998, 21, 148–155. [Google Scholar] [CrossRef]
  9. Karlgaard, C.D.; Lutze, F.H. Second-order relative motion equations. J. Guid. Control Dyn. 2003, 26, 41–49. [Google Scholar] [CrossRef] [Green Version]
  10. Ross, I.M. Linearized dynamic equations for spacecraft subject to j perturbations. J. Guid. Control Dyn. 2003, 26, 657–659. [Google Scholar] [CrossRef]
  11. Vadali, S.R. Model for linearized satellite relative motion about a J2-perturbed mean circular orbit. J. Guid. Control Dyn. 2009, 32, 1687–1691. [Google Scholar] [CrossRef]
  12. Guffanti, T.; D’Amico, S. Linear models for spacecraft relative motion perturbed by solar radiation pressure. J. Guid. Control Dyn. 2019, 42, 1962–1981. [Google Scholar] [CrossRef]
  13. Franzini, G.; Innocenti, M. Relative Motion Dynamics with Arbitrary Perturbations in the Local-Vertical Local-Horizon Reference Frame. J. Astronaut. Sci. 2020, 67, 98–112. [Google Scholar] [CrossRef]
  14. Pollock, G.E.; Gangestad, J.W.; Longuski, J.M. Analytical solutions for the relative motion of spacecraft subject to lorentz–force perturbations. Acta Astronaut. 2011, 68, 204–217. [Google Scholar] [CrossRef]
  15. Yamada, K.; Kimura, M.; Shima, T.; Yoshikawa, S. New state transition matrix for formation flying in j2-perturbed elliptic orbits. J. Guid. Control Dyn. 2012, 35, 536–547. [Google Scholar] [CrossRef]
  16. Jiang, W.; Li, J.; Jiang, F.; Bernelli-Zazzera, F. A simple method to design non-collision relative orbits for close spacecraft formation flying. Sci. China Phys. Mech. Astron. 2018, 61, 054511. [Google Scholar] [CrossRef]
  17. Bando, M.; Ichikawa, A. Periodic orbits of nonlinear relative dynamics and satellite formation. J. Guid. Control Dyn. 2009, 32, 1200–1208. [Google Scholar] [CrossRef]
  18. Alfriend, K.; Yan, H. Evaluation and comparison of relative motion theories. J. Guid. Control Dyn. 2005, 28, 254–261. [Google Scholar] [CrossRef]
  19. Glassmeier, K.-H.; Boehnhardt, H.; Koschny, D.; Kührt, E.; Richter, I. The rosetta mission: Flying towards the origin of the solar system. Space Sci. Rev. 2007, 128, 1–21. [Google Scholar] [CrossRef]
  20. Capó-Lugo, P.A.; Bainum, P.M. Orbital Mechanics and Formation Flying: A Digital Control Perspective; Elsevier: Amsterdam, The Netherlands, 2011. [Google Scholar]
  21. Sullivan, J.; Grimberg, S.; D’Amico, S. Comprehensive survey and assessment of spacecraft relative motion dynamics models. J. Guid. Control Dyn. 2017, 40, 1837–1859. [Google Scholar] [CrossRef]
  22. Canuto, E.; Novara, C.; Carlucci, D.; Montenegro, C.P.; Massotti, L. Spacecraft Dynamics and Control: The Embedded Model Control Approach; Butterworth-Heinemann: Oxford, UK, 2018. [Google Scholar]
  23. Andreev, V.; Goryainov, S.; Krasilnikov, A. Six–body problem solution using symplectic integrators. In Proceedings of the 2016 IEEE NW Russia Young Researchers in Electrical and Electronic Engineering Conference (EIConRusNW), St. Petersburg, Russia, 2–3 February 2016; pp. 116–119. [Google Scholar]
  24. Andreev, V.; Goryainov, S.; Krasilnikov, A. N-body problem solution with composition numerical integration methods. In Proceedings of the 2016 XIX IEEE International Conference on Soft Computing and Measurements (SCM), St. Petersburg, Russia, 25–27 May 2016; pp. 196–198. [Google Scholar]
  25. Nayfeh, A.H. Introduction to Perturbation Techniques; John Wiley & Sons: Hoboken, NJ, USA, 2011. [Google Scholar]
  26. Abouelmagd, E.I.; Elshaboury, S.; Selim, H. Numerical integration of a relativistic two-body problem via a multiple scales method. Astrophys. Space Sci. 2016, 361, 38. [Google Scholar] [CrossRef]
  27. Abouelmagd, E.I.; Guirao, J.L.G.; Pal, A.K. Periodic solution of the nonlinear Sitnikov restricted three-body problem. New Astron. 2020, 75, 101319. [Google Scholar] [CrossRef]
  28. Lü, J.; Wang, Q.; Wang, S. Analytic approach on geometric structure of invariant manifolds of the collinear lagrange points. Sci. China Phys. Mech. Astron. 2012, 55, 1703–1712. [Google Scholar] [CrossRef]
  29. Abouelmagd, E.I.; Mortari, D.; Selim, H.H. Analytical study of periodic solutions on perturbed equatorial two-body problem. Int. J. Bifurc. Chaos 2015, 25, 1540040. [Google Scholar] [CrossRef]
  30. Abouelmagd, E.I.; Guirao, J.L.; Vera, J.A. Dynamics of a dumbbell satellite under the zonal harmonic effect of an oblate body. Commun. Nonlinear Sci. Numer. Simul. 2015, 20, 1057–1069. [Google Scholar] [CrossRef] [Green Version]
  31. Pal, A.K.; Kushvah, B.S. Geometry of halo and Lissajous orbits in the circular restricted three–body problem with drag forces. Mon. Not. R. Astron. Soc. 2015, 446, 959–972. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Configuration of relative motion satellite.
Figure 1. Configuration of relative motion satellite.
Symmetry 13 00595 g001
Figure 2. Solutions of linear relative motion in two dimensions. (a) Analytical solution of linear relative motion at different phases. (b) Numerical solutions by using numerical integration of explicit Euler method.
Figure 2. Solutions of linear relative motion in two dimensions. (a) Analytical solution of linear relative motion at different phases. (b) Numerical solutions by using numerical integration of explicit Euler method.
Symmetry 13 00595 g002
Figure 3. Solutions of linear relative motion in three dimensions (a) Analytical solutions at different phases when φ 0 = π / 4 . (b) Numerical solutions by using numerical integration of the explicit Euler method.
Figure 3. Solutions of linear relative motion in three dimensions (a) Analytical solutions at different phases when φ 0 = π / 4 . (b) Numerical solutions by using numerical integration of the explicit Euler method.
Symmetry 13 00595 g003
Figure 4. Solution of nonlinear relative motion in two dimensions. (a) Solutions of the Lindstedt–Poincaré technique at different phases. (b) Numerical solutions by integration of the explicit Euler method.
Figure 4. Solution of nonlinear relative motion in two dimensions. (a) Solutions of the Lindstedt–Poincaré technique at different phases. (b) Numerical solutions by integration of the explicit Euler method.
Symmetry 13 00595 g004
Figure 5. Solution of nonlinear relative motion in three dimensions; (a) Solutions of the Lindstedt–Poincaré technique at different phases when φ 0 = π / 4 . (b) Numerical solutions by integration of the explicit Euler method.
Figure 5. Solution of nonlinear relative motion in three dimensions; (a) Solutions of the Lindstedt–Poincaré technique at different phases when φ 0 = π / 4 . (b) Numerical solutions by integration of the explicit Euler method.
Symmetry 13 00595 g005
Figure 6. Comparison between analytical solutions of linear and nonlinear relative motion in two dimensions at different phases.
Figure 6. Comparison between analytical solutions of linear and nonlinear relative motion in two dimensions at different phases.
Symmetry 13 00595 g006
Figure 7. Comparison between analytical solutions of linear and nonlinear relative motion in three dimensions at different phases when φ 0 = π / 4 .
Figure 7. Comparison between analytical solutions of linear and nonlinear relative motion in three dimensions at different phases when φ 0 = π / 4 .
Symmetry 13 00595 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pal, A.K.; Abouelmagd, E.I.; García Guirao, J.L.; Brzeziński, D.W. Periodic Solutions of Nonlinear Relative Motion Satellites. Symmetry 2021, 13, 595. https://doi.org/10.3390/sym13040595

AMA Style

Pal AK, Abouelmagd EI, García Guirao JL, Brzeziński DW. Periodic Solutions of Nonlinear Relative Motion Satellites. Symmetry. 2021; 13(4):595. https://doi.org/10.3390/sym13040595

Chicago/Turabian Style

Pal, Ashok Kumar, Elbaz I. Abouelmagd, Juan Luis García Guirao, and Dariusz W. Brzeziński. 2021. "Periodic Solutions of Nonlinear Relative Motion Satellites" Symmetry 13, no. 4: 595. https://doi.org/10.3390/sym13040595

APA Style

Pal, A. K., Abouelmagd, E. I., García Guirao, J. L., & Brzeziński, D. W. (2021). Periodic Solutions of Nonlinear Relative Motion Satellites. Symmetry, 13(4), 595. https://doi.org/10.3390/sym13040595

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop