# TheJ2 Relative Perturbation Analysis of Satellite Formation under the Requirement of Relative Position Maintenance with Millimeter-Level Accuracy

Authors: Jingji Wang; Chunyang Liu
Journal: International Journal of Aerospace Engineering
 (2021)

---

## Abstract

With high-precision DEM (Digital Elevation Model) and GMTI (Ground Moving Target Indicator) as the demand background, the influence ofJ2 zonal harmonic term perturbation on the relative motion of the millimeter-level short-range leader-follower satellites in near-circular orbit is studied through the relative perturbation method. An equation of motion that can describe the motion of the leader-follower satellites under the influence of J2 perturbation in near-circular orbit is derived, and the characteristics of the trajectory of in-plane periodic motion are analyzed. A study shows that under the influence of the relative perturbation of the J2 term, the in-plane periodic motion of the leader-follower satellites in near-circular orbit is a symmetrical closed “drop-shaped” trajectory with a period of 2π/n. By comparing with the results of numerical simulations, the correctness of the conclusions obtained in this paper is verified. According to the research results, it can be known that only using a thruster as the actuator to maintain the relative position can no longer meet the requirements of the long-term mm-level relative position maintenance. In the future, a new technical approach needs to be explored to achieve the long-term relative position maintenance with millimeter-level control accuracy.


---

## Body

## 1. Introduction

The satellite formation system [1] refers to a space system composed of multiple physically disconnected satellites in order to complete the same mission. Satellite formation can not only complete complex tasks that cannot be completed by a single satellite but also have more advantages than a single satellite in terms of reliability, flexibility, and economy. Therefore, satellite formation systems have become a research hotspot in the aerospace field and are widely used in many other fields, such as earth observation, astronomical observation, and deep space exploration. This article takes the high-precision DEM (Digital Elevation Model) and GMTI (Ground Moving Target Indicator) tasks as the background. As shown in Figure 1, the leader-follower satellite means that the deputy satellite and the chief satellite operate in the same orbit in a certain sequence, and the orbit elements of leader-follower satellites are only different in argument of latitude. Therefore, the leader-follower satellite is equivalent to a sparse antenna array, using ground echo interferometric imaging principles and formation flying satellite control technology to complete DEM and GMTI tasks.Figure 1
Diagram of DEM and GMTI tasks.Whether the satellite formation system can meet the DEM mission measurement accuracy and the minimum detectable speed of the GMTI mission not only depends on the performance of payload, length of baseline, and data processing but also depends on the control accuracy of relative position maintenance of the payload of the leader-follower satellites. In low-precision observation tasks, it only needs to control the relative position accuracy of the interference baseline at the meter or decimeter level. Therefore, in the study of relative dynamics, only the secular and long-periodic effects of various perturbations on the relative position need to be paid attention to. As users have higher and higher requirements for performance and accuracy, the accuracy requirements for the relative position maintenance of formation flying satellite have also been increased from the meter level and decimeter level to the millimeter level, which raises the higher control accuracy of the relative position maintenance of the satellite formation. Therefore, it is necessary not only to offset the effects of the secular and long-periodic contributions due to perturbation, but also to compensate for the short-periodic impact of perturbation. Hence, it is important to conduct a detailed analysis of the shape of the trajectory of formation flying satellites under the impact of perturbation. Among the many perturbations, the influence of theJ2 term perturbation occupies a major position. So this paper’s aim is mainly to conduct a detailed analysis of the influence of the J2 perturbation on the relative position of the leader-follower satellites.A lot of research efforts had been done in the field of satellite relative motion under the influence ofJ2 perturbation based on state transition matrices (STMs) and relative orbital elements (ROEs) [2]. Schweighart and Sedwick [3] derived a linear model by incorporating the first-order term of the J2 perturbation into the C-W equation with a meter-level position error when compared with a numerical integration. Izzo et al. [4] considered the atmospheric drag and J2 perturbation in the extension of the C-W model and used the Floquet theory to predict the relative motion under the perturbations with a higher accuracy than that of the C-W equation. Stringer et al. [5] used the quadratic Volterra method developing a nonlinear model for including the J2 perturbation with submeter-level accuracy.Unlike the above research, Schaub and Alfriend [6] and Gim and Alfriend [7] considered that the ROEs vary slowly with time and developed the ROE-based equations for relative motion under J2 perturbation. Schaub and Alfriend [6] used the J2-perturbed Hamiltonian deriving a relative state equation defined by Delaunay element differences and explored the impact of J2 perturbation on relative orbits. Gim and Alfriend [7] developed a state transition matrix (STM) including the effect of secular, long, and short-periodic contributions due to J2 to the first order. Gim and Alfriend [8] then developed a model that can avoid the singularity on equatorial orbits by using differential equinoctial elements. D’Amico [9] developed a linear relative motion model which captures the time evolution of quasinonsingular ROE subjects (relative semimajor axis, the relative mean longitude, the relative eccentricity, and inclination vectors) to first-order J2 perturbations in near-circular orbits. By using the same ROE state, Gaias et al. [10] developed a more complete relative motion model which corrects the original J2 formulation shortcomings of D’Amico’s [9] work and includes the effects of atmospheric drag, with a higher accuracy than that of D’Amico’s [9] model. In addition, Mahajan et al. [11] developed a STM for relative motion under the complete zonal perturbations based on the studies of Gim and Alfriend [7]. More recent work by Russell et al. [12] developed an analytical STM for relative motion including J2 and J3 perturbations by using the Vinti theory. Burnett et al. [13] derived a new linearized differential equation model including zonal harmonic perturbations and corresponding perturbation solutions for relative motion without averaging the perturbing acceleration or its gradient. The perturbation solution of their model depends only on time, the initial chief orbital elements, and the deputy initial conditions in the chief Hill frame and has a similar error to the GA-STM for zero initial chief orbit eccentricity.As a nonlinear extension to the ROE-based linear relative motion, Alfriend [14] and Alfriend and Yan [15] derived a nonlinear model for propagating the mean OED in arbitrarily eccentric orbits while including second-order J2 effects. Based on the Gim-Alfriend’s studies, Sengupta et al. [16] developed a second-order state transition tensor (STT) and a nonlinear solution including the effects of J2 perturbation by combining the Keplerian STT with Gim-Alfriend’s [7] linear model. Then, Yang et al. [17] derived a more complete second-order STT for relative motion under the J2-perturbed elliptic orbits based on the Sengupta et al.’s [16] work. The two studies above [16, 17] did not consider the effect of J2 perturbation of the second-order transformation from the rectilinear relative state to the osculating ROEs and the second-order propagation of the mean ROEs. Similar to the methods in [16, 17], Zhen et al. [18] developed a second-order analytical STT for relative motion under the J2-perturbed elliptic orbits by using the geometric method, and the result is more accurate than that of previous linear or nonlinear analytic methods. Moreover, Mahajan et al. [19] considered the J2 zonal harmonic and sectorial and tesseral harmonics in their work, presenting a STM for perturbed satellite relative motion with the position errors below 1% of the formation size after two days in both the lowly and highly eccentric reference orbits. Recently, Gaias et al. [20] presented an analytical framework for the precise modeling of the relative motion in low earth orbits. Their model includes the first-order expansion of the effects due to any even zonal harmonics and the second-order expansion of the unperturbed and J2 terms and can achieve high-precision mean/osculating orbital element conversions.In addition, in order to achieve high-precision relative orbit maintenance, some scholars studied the hovering formation and orbit determination on the basis of the above dynamic research in recent years. Rao et al. [21] obtained a “teardrop” hovering formation by designing a set of relative orbit elements and proposed a new impulsive control strategy to keep the deputy satellite in the hovering pattern for a long time. Sun et al. [22] proposed an innovative orbit determination method which makes use of gravity gradients for low earth-orbiting satellites. Bai et al. [23] proposed a teardrop hovering formation design for the elliptical reference orbit with a J2 perturbation.In response to the above-mentioned background and research status, this paper studies the shape of relative motion of the leader-follower satellites in near-circular orbit under the effect ofJ2 perturbation. On the basis of the predecessors, an equation of motion that can describe the relative motion of the leader-follower satellites in near-circular orbit under the influence of J2 perturbation is derived, and the characteristics of in-plane periodic motion of the leader-follower satellites under the influence of J2 term perturbation are analyzed. The equation of motion derived in this paper has a very simple form and is more desirable and achievable in onboard computers. By comparing with the result of numerical simulations under the same conditions, the relative position error of this equation is less than 20 mm in the radial direction and 60 mm in the tangential direction in a day. The major contributions of this study are summarized as follows:(1)
We derive a high-precision equation of motion that can describe the relative motion of the leader-follower satellites in near-circular orbit under the influence ofJ2 perturbation. The equation has a very simple form and is more desirable and achievable in onboard computers, so it can be used as a state equation for high-precision relative navigation(2)
We analyze the characteristics of relative motion of leader-follower satellites under the influence ofJ2 perturbation, which can provide a theoretical basis for the design of controllers of mm-level relative orbit maintenanceThe outline of this paper is as follows: the coordinate system and the mathematical description of formation flying are introduced in Section2. An equation of motion that can describe the relative motion of the leader-follower satellites in near-circular orbit under the influence of J2 perturbation is derived, and the characteristics of in-plane periodic motion of the leader-follower satellites under the influence of J2 term perturbation are analyzed in Section 3. Numerical simulations used to validate the results in this paper are shown in Section 4. In Section 5, we summarize the study and propose a new method for future relative orbit maintenance.


## 2. The Mathematical Description of Formation Flying

In order to derive the equation of the motion of leader-follower satellites under the influence of theJ2 relative perturbation, it is first necessary to establish a coordinate system and a mathematical model to describe the relative motion [24] of the satellite formation. As shown in Figure 2, denote the chief satellite as S and the deputy satellite as C. The orbital coordinate system OS−xyz of the chief satellite S is defined as follows: the origin OS is located at the center of the chief satellite and moves with it, the x-axis takes the radial direction of the reference satellite and points from the center of the earth to S, the z-axis points to the positive normal direction of the chief satellite S orbital plane, and the y-axis is determined by the right-hand rule.Figure 2
Relative motion coordinate system.Assuming that the chief satelliteS is moving in a near-circular orbit, the following C-W equation [25] can be obtained to describe the relative motion while ignoring high-order small quantities (since the distance between the chief satellite and the deputy satellite ρ is much smaller than the geocentric distance of the chief satellite r, we can ignore the ρ/r2 and higher power terms):(1)d2xdt2−2ndydt−3n2x=∆fx,d2ydt2+2ndxdt=∆fy,d2zdt2+n2z=∆fz,where x,y,z is the relative position of the deputy satellite C in the orbital coordinate system of the chief satellite S, n=μ/a3 is the average angular velocity of the chief satellite S, μ is the geocentric gravitational constant, and a is the semimajor axis of the chief satellite. ∆fx, ∆fy, and ∆fz are the differences between the perturbation acting on the chief satellite and the deputy satellite, that is, the relative perturbation.Taking the initial moment as the moment when the chief satellite is located at the RAAN (Right Ascension of Ascending Node), in a near-circular orbit, the argument of latitude of the chief satellite is(2)ut=nt.Then, there is(3)ddt∙=nddu∙,d2dt2∙=n2d2du∙.In the context of this paper, the effect of the perturbation on the argument of latitude of the satellite is secondary, so we can neglect it and substitute formula (3) into formula (1); then, the independent variable t of the C-W equation can be changed to u, and we have(4)d2xdu2−2dydu−3x=1n2∆fx,d2ydu2+2dxdu=1n2∆fy,d2zdu2+z=1n2∆fz.In order to facilitate the analysis of the motion deviation of the deputy satellite in the relative motion orbit under the influence of relative perturbation, let∆x,∆y,∆z replace x,y,z in the above formula; then, equation (4) can be changed to(5)d2∆xdu2−2d∆ydu−3∆x=1n2∆fx,d2∆ydu2+2d∆xdu=1n2∆fy,d2∆zdu2+∆z=1n2∆fz,where ∆x=x−xp, ∆y=y−yp, ∆z=z−zp, x,y,z is the position of the deputy satellite in the relative motion orbit only under the influence of two-body gravity, and xp,yp,zp is the actual position of the deputy satellite in the relative motion orbit under the influence of relative perturbation and two-body gravity. The solution ∆x,∆y,∆z of equation (5) can describe the relative motion only under the influence of relative perturbation.


## 3.J2 Relative Perturbation Analysis

In this section, the relative perturbation method is used to derive the equation of relative motion of the leader-follower satellites under the influence ofJ2 perturbation. The gravitational potential function of the zonal harmonic term of the earth [26] is(6)U=μr1−J2Re22r23sin2φ−1−J3Re32r35sin3φ−3sinφ−J4Re48r435sin4φ−30sin2φ+3,where Re is the radius of the earth, r,φ is the geocentric distance and geocentric latitude of the spacecraft in the ground-fixed coordinate system, and J2, J3, and J4 are the coefficients of zonal harmonic term perturbation.Then, the perturbation potential function ofJ2 (the central gravitational term is omitted) is(7)∆U=−μJ2Re22r33sin2φ−1.The components of the perturbation force in spherical coordinates are(8)fr=∂∆U∂r=32J2μRe2r43sin2φ−1fα=∂∆U∂α1rcosφ=0fφ=∂∆U∂φ=−32J2μRe2r4sin2φ.From the relationship between satellite spherical coordinates and orbital coordinates, it can be seen that the radial perturbation forces under satellite spherical coordinates and orbital coordinates are the same, and the tangential and normal perturbation forcesfy and fz have the following conversion relationship(9)fy=fφcosβfz=fφsinβ,where β is the angle between the orbital plane and the meridian plane of the point where the satellite is located. From the spherical triangle relationship, we have(10)cosi=cosφsinβsinφ=sinisinucosβ=tanφcotu.From equations (6), (7), (8), (9), and (10), we have the perturbation force of the J2 zonal harmonic term in orbital coordinates:(11)fx=−k1−3sin2isin2ufy=−ksin2isin2ufz=−ksin2isinu,where k=3/2J2μRe2/r4 and i and u are the inclination and argument of latitude of the spacecraft in the ground-fixed coordinate system.In equation (11), by taking partial differentiation of r, i, and u, we can obtain the following equation:(12)∆fx=k41−3sin2isin2ur∙∆r+3sin2isin2u∙∆i+3sin2isin2u∙∆u,∆fy=k4sin2isin2ur∙∆r−sin2isin2u∙∆i−2sin2icos2u∙∆u,∆fz=k4sin2isinur∙∆r−2cos2isinu∙∆i−sin2icosu∙∆u.In the above formula,∆r is the difference of the geocentric distance between the deputy satellite and the chief satellite, ∆i is the difference of the inclination between the deputy satellite and the chief satellite, and ∆u is the difference of the argument of latitude between the deputy satellite and the chief satellite.As we all know, the orbit elements of leader-follower satellites are only different in argument of latitude. Therefore, in a near-circular orbit, there are no differences in the geocentric distance and inclination between the chief and deputy satellite, that is,∆r=∆i=0. Then, equation (12) can be changed to(13)∆fx=3ksin2isin2u∙∆u,∆fy=−2ksin2icos2u∙∆u,∆fz=−ksin2icosu∙∆u.Assuming that the deputy satellite has no relative position and velocity deviation between the influence of relative perturbation and two-body gravity at the initial moment, the initial value is(14)∆x0=∆y0=∆z0=0,(15)d∆xduu=0=d∆yduu=0=d∆zduu=0=0.Substituting equations (14) and (15) into equation (5), we can get the following analytical solution:(16)∆x=2k3n2sin2isinu1−cosu∙∆u,(17)∆y=−k3n2sin2i5+cosu1−cosu∙∆u,(18)∆z=−k2n2usin2isinu∙∆u.In formulas (16), (17), and (18), ∆x,∆y is the in-plane motion trajectory of the leader-follower satellites in near-circular orbit only under the influence of J2 perturbation, and ∆z reflects the out-of-plane motion only under the influence of J2 perturbation. It can be seen from the above formula that the relative position deviation in the three directions under the influence of J2 perturbation is positively correlated with the difference of argument of latitude ∆u between the deputy satellite and the chief satellite. Therefore, the greater the distance between the leader-follower satellites, the greater the relative position error of the formation flying satellite.Because the relative motion of short-range leader-follower satellites in thez direction under the influence of J2 perturbation is less than 1 mm, we mainly discuss in-plane motion. For in-plane motion, let(19)A=k3n2∆usin2i.Then, formulas (16) and (17) can be simplified as(20)∆x=2Asinu1−cosu,∆y=−A5+cosu1−cosu.The parametric equation represented by equation (20) describes the in-plane motion of leader-follower satellites in near-circular orbit only under the influence of J2 perturbation. In equation (20), by finding the first-order and second-order derivatives for t, we can obtain the velocity and acceleration of the in-plane motion only under the influence of J2 perturbation as follows:(21)∆ẋ=2nAcosu−cos2u,∆ẏ=−nA4sinu+sin2u,(22)∆x¨=2n2A2sin2u−sinu,∆y¨=−n2A4cosu+2cos2u.By analyzing equations (20), (21), and (22), it can be obtained that the in-plane motion of leader-follower satellites in near-circular orbit only under the influence of J2 perturbation has the following law:(a)
The curve described by the parametric equation (20) is a symmetrically closed “drop-shaped” line as shown in Figure 3, and ∆x=0 is the symmetry axis of the curve. Let t=0 be the starting time of the in-plane motion, which corresponds to O0,0 in the figure. The deputy satellite moves periodically in the clockwise direction from point O, and the period is 2π/n. The detailed parameters at the special points in the motion are shown in Table 1(b)
In formula (20), the parameter “A” determines the magnitude of the in-plane motion only under the influence of J2 perturbation. When the distance between the leader-follower satellites is certain, the magnitude of “A” is also related to the orbital altitude and the inclination of the chief satellite. As the orbital altitude increases, “A” continues to decrease, and the in-plane relative position deviation of the leader-follower satellites also decreases. When 0°<i<90°, “A” increases with the increase in inclination i, and the in-plane relative position deviation of leader-follower satellites also increases. When i=90°, the in-plane relative position deviation reaches its maximum value. When 90°<i<180°, “A” decreases as the inclination i increases, and the in-plane relative position deviation of the leader-follower satellites also decreases(c)
Figures4 and 5 show the velocity and acceleration of the deputy satellite’s in-plane motion in near-circular orbit only under the influence of J2 relative perturbation at an orbital altitude of 500 km. It can be seen from the figure that in order to control this in-plane motion within the mm level, the required acceleration is of the order of 10−7m/s2, and the thruster is required to provide a continuous thrust of about 0.1 mN. In order to eliminate this movement caused by the relative perturbation of J2, more frequent thruster control must be relied on. The long-term frequent jets of the thruster will increase fuel consumption and greatly reduce the service life of the satellites. Therefore, at this stage, simply using the thruster as the actuator cannot achieve mm-level long-term relative position maintaining controlFigure 3
In-plane motion trajectory of leader-follower satellites under the influence ofJ2 relative perturbation.Table 1
Special point parameters in the process of the in-plane motion only under the influence ofJ2 relative perturbation.
ORPQu (rad)2kπ, k=0,1,2⋯2π/3+2kπ, k=0,1,2⋯π+2kπ, k=0,1,2⋯4π/3+2kπ, k=0,1,2⋯t (s)2kπnπn23+2k2k+1πnπn43+2k∆x (km)033/2A (positive maximum)0−33/2A (negative maximum)∆y (km)0−274A-8A (maximum)−274A∆ẋ (km·s-1)00-4 nA0∆ẏ (km·s-1)0−532nA0−332nAFigure 4
In-plane motion velocity of leader-follower satellites under the influence ofJ2 relative perturbation.Figure 5
In-plane motion acceleration of leader-follower satellites under the influence ofJ2 relative perturbation.


## 4. Simulation

According to the background, the simulation parameters select the commonly used orbital parameters for earth observation using synthetic aperture radar satellites. The two satellites are, respectively, selected for simulation in the form of forward and backward flying formation with a distance of 100 m and 200 m. Select a low-earth circular orbit with an orbit altitude of 500 km, and the remaining initial orbit parameters are shown in Tables2 and 3. The numerical simulator takes the initial analytic equations of motion as described in equation (23) and integrates the absolute motion of each satellite in time. The motion of each satellite is then differenced and converted to a coordinate system that was mentioned in Section 2.(23)r¨=−μr3r+FJ2.Table 2
Initial orbit parameters when the two satellites are 200 m apart.
Parameter nameSatellite 1Satellite 2Semimajor axis (km)6878.146878.14Eccentricity00Inclination (°)97.406597.4065RAAN (°)00Argument of perigee (°)00True anomaly (°)00.00166603Table 3
Initial orbit parameters when the two satellites are 100 m apart.
Parameter nameSatellite 1Satellite 2Semimajor axis (km)6878.146878.14Eccentricity00Inclination (°)97.406597.4065RAAN (°)00Argument of perigee (°)00True anomaly (°)00.000833015In order to eliminate the error of integration, we take the initial analytic equations of motion as described in equation (24) as a numerical simulator and integrate the absolute motion of each satellite in time once more and perform the same processing as above. By subtracting the above two results, the relative motion result of the satellites only under the J2 perturbation can be obtained, which can be compared with the analytical solution result in this paper. By comparing the conclusions obtained in this paper with the results of the numerical simulation, the correctness of the results obtained in this paper is verified.(24)r¨=−μr3r.The simulation results shown in Figures6(b) and 7(b) show that under the influence of the relative perturbation of J2, the in-plane periodic motion of the leader-follower satellites is a symmetrical closed “drop-shaped” trajectory, which is coincident with the results of numerical simulations (the results shown in Figures 6(a) and 7(a)), and the relative position error of this equation is less than 20 mm in the radial direction and 60 mm in the tangential direction in a day, as shown in Figure 8.Figure 6
In-plane motion only under the influence ofJ2 perturbation when two satellites are 200 m apart.
(a)
The result of numerical simulation(b)
Analytical solution resultFigure 7
In-plane motion only under the influence ofJ2 perturbation when two satellites are 100 m apart.
(a)
The result of numerical simulation(b)
Analytical solution resultFigure 8
The error between the result of the analytical solution and the result of the numerical simulations.


## 5. Conclusions

In this paper, the relative perturbation method is used to study the influence of the relative perturbation ofJ2 zonal harmonic term perturbation on the relative motion of leader-follower satellites under the requirement of millimeter-level accuracy in near-circular orbit. Studies have shown that under the influence of the J2 relative perturbation, the in-plane periodic motion of the leader-follower satellites in near-circular orbit is a symmetrical closed “drop-shaped” trajectory with a period of 2π/n. In order to eliminate this motion caused by the J2 relative perturbation, more frequent thruster control must be relied on. The long-term frequent jets of the thruster will increase fuel consumption and greatly reduce the service life of the satellite. Therefore, only using a thruster as the actuator to maintain the relative position can no longer meet the requirements of the long-term mm-level relative position maintenance. In the future, we can follow the principle of flywheel angular momentum exchange and use the principle of linear momentum exchange for relative position control. By installing a two-dimensional sliding mechanism on the two satellites and controlling the sliding mechanism to move, the linear momentum exchange caused by the continuous motion of the sliding mechanism can change the relative position of the formation satellites. Since the relative motion under the influence of J2 perturbation is a continuous closed periodic trajectory and the motion amplitude is limited, this method can well offset this “drop-shaped” motion. Moreover, this method does not need to carry additional fuel and can provide continuous control power only by relying on electric energy. It has advantages over the use of thrusters in terms of accuracy and fuel saving.



---

*Source: 6687468-2021-02-26.xml*