Next Article in Journal
Experimental Investigation of Flame Anchoring Behavior in a LOX/LNG Rocket Combustor
Next Article in Special Issue
Design of Aerospace Vehicles’ Thermal Protection Based on Heat-Insulating Materials with Optimal Structure
Previous Article in Journal
A Coaxial Pulsed Plasma Thruster Model with Efficient Flyback Converter Approaches for Small Satellites
Previous Article in Special Issue
Low-Dimensional Models for Aerofoil Icing Predictions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Tether System at the L1, L2 Collinear Libration Points of the Mars–Phobos System: Analytical Solutions

by
Vladimir S. Aslanov
* and
Daria V. Neryadovskaya
Department of Theoretical Mechanics, Samara National Research University, 34, Moskovskoye Shosse, Samara 443086, Russia
*
Author to whom correspondence should be addressed.
Aerospace 2023, 10(6), 541; https://doi.org/10.3390/aerospace10060541
Submission received: 5 May 2023 / Revised: 1 June 2023 / Accepted: 3 June 2023 / Published: 5 June 2023

Abstract

:
This paper is dedicated to identifying stable equilibrium positions of the tether systems attached to the L1 or L2 libration points of the Mars–Phobos system. The orbiting spacecraft deploying the tether is at the L1 or L2 libration point and is held at one of these unstable points by the low thrust of its engines. In this paper, the analysis is performed assuming that the tether length is constant. The equation of motion for the system in the polar reference frame is obtained. The stable equilibrium positions are found and the dependence of the tether angular oscillation period on the tether length is determined. An analytical solution in the vicinity of the stable equilibrium positions for small angles of deflection of the tether from the local vertical is obtained in Jacobi elliptic functions. The comparison of the numerical and analytical solutions for small angles of deflection is performed. The results show that the dependencies of the oscillation period on the length of the tether are fundamentally different for L1 and L2 points. Analytical expressions for the tether tension are derived, and the influence of system parameters on this force is investigated for static and dynamic cases.

1. Introduction

The space tether is a type of tether that is made of high-strength fiber used to connect spacecraft to each other or to other masses. The space tether systems allow us to perform missions that are impossible, impractical or unprofitable to accomplish with the help of other space equipment. For example, tether systems can be used for docking between spacecraft [1], as a space elevator [2,3,4,5,6], for payload orbital transfer [7,8], for exploring deep space [9], the atmosphere and the surface of the planets and their moons [10], as well as asteroids. For instance, Mashayekhi and Misra studied the effect of attaching a tether and ballast mass to an asteroid with subsequent cutting of the tether [11]. In recent decades, the use of space tethers near the collinear Lagrangian points has received considerable attention [12,13,14,15]. Ref. [14] focused on the development of a new mission to explore Phobos using a tether system anchored below the L1 Mars–Phobos libration point and deployed toward Mars at a length slightly greater than the distance from Phobos. Paper [15] showed the maintenance of an L1-type artificial equilibrium point in the Sun (Earth + Moon) circular restricted three-body problem by means of an electric solar wind sail. The tether capture system is also a promising method for removing space debris [16,17,18,19,20,21]. The topic of dynamics and control of tether systems has received substantial attention [22,23,24,25,26,27,28,29,30,31,32,33,34,35,36]. Huang et al. examined several new applications for the space tether during operation in orbit, focusing on the structure, dynamics and control [23]. Paper [24] discussed the diversity of tether modelling that has been undertaken recently, and showed that dynamics and control are the two fundamentally important aspects of all tether concepts, designs and mission architectures.
In 2017, NASA proposed the PHLOTE mission (Phobos L1 Operational Tether Experiment) to explore the surface of Phobos using a tether system ‘‘anchored’’ at the L1 libration point of the Mars–Phobos system [10]. The tether release point was proposed to be an orbiting spacecraft hovering in the vicinity of the L1 point. Once deployed, the small vehicle with a sensor package attached to the tether was expected to investigate Phobos. This mission concept is a synthesis of new technologies that would provide a unique platform for multiple sensors directed at Phobos as well as at Mars. The PHLOTE ConOps describes the PHLOTE mission and provides a key systems engineering document to support future mission development. However, such a complex innovative mission requires an additional theoretical justification and a variety of analytical models of the system motion. The work of [37] considers a mission similar to the PHLOTE mission, where a detailed study of the behavior of the tether system attached at the L1 collinear libration point was performed using the classical Nehvil equations.
The purpose of the present work is to find the stable equilibrium positions of the tether system attached at the L1 or L2 collinear libration points of the Mars–Phobos system and to study the features of the tether motion near these positions. The system consists of the tether and the end mass attached to its end. The mathematical model is based on the differential equations of the classical circular restricted three-body problem [38,39,40,41,42]. The equation of motion for the tether system of constant length under the action of two gravitational fields (Mars–Phobos) and the centrifugal force associated with the rotation of the frame of the Mars–Phobos system are obtained in polar coordinates. The first integral of this differential equation is found and used to determine the phase trajectories and the stable equilibrium positions. The approximate analytical solutions of the equation of motion for the tether system are obtained using Jacobi elliptic functions. Next, the dependence of the oscillation period on the length of the tether is found. Finally, analytical expressions for the tether tension are derived, and the influence of system parameters on this force is investigated for static and dynamic cases. The results of this work can be used for PHLOTE-like mission design. It is worth noting that the obtained solutions for small tether deflection angles are of interest for the creation of the space elevator at the L1 and L2 libration points of the Mars–Phobos system in the future.

2. Mathematical Model: Finding Sustainable Positions

In this section, the behavior of the tether system of constant length attached to the L 1 or L 2 libration point under the action of two massive attracting bodies, M 1 and M 2 (Mars and Phobos), is described using the differential equations of the classical circular restricted three-body problem [38,39,40,41,42]. It is assumed that the mass of the body M is much less than the mass of the bodies M 1 and M 2 . As a result, the body M has negligible influence on other bodies. In addition, it is assumed that the eccentricity of the two bodies of the primary orbit is e = 0 and the distance between them is
r = d ,
where d is the semilatus rectum. The orbiter is located at the L1 or L2 libration points, either of which can be the attachment point for the tether.
In the following subsections, we consider two cases characterized by different values of the tether deflection angle, namely, φ and ψ = φ + π .

2.1. Tether Deflection Angle φ

The equations of motion of the circular restricted three-body problem [42] in the polar reference frame l , φ (see Figure 1) for the constant length tether l = const can be written as
φ ¨ F i φ = 0 , i = 1 , 2
where
F i φ = n 2 sin φ a i l + G sin φ l m 1 d μ + a i r 1 3 + m 2 d μ 1 + a i r 2 3
n = G m 1 + m 2 d 3 is the mean motion, G is the Newtonian gravitational constant, m 1 and m 2 are masses of the bodies M 1 and M 2 , respectively, a 1 d [   1 μ 3 1 3   ] , a 2 d [   1 + μ 3 1 3   ] are the distances from the origin to the L 1 and L 2 libration points, respectively μ = m 2 m 1 + m 2 is the mass ratio,
r 1 is the distance between the primary 1 and the end mass,
r 1 = a i + d μ + l cos φ 2 + l sin φ 2 ,   i = 1 , 2 ,
r 2 is the distance between the primary 2 and the end mass,
r 2 = a i + d μ 1 + l cos φ 2 + l sin φ 2 ,   i = 1 , 2
Equation (2) has the following energy integral:
E = φ ˙ 2 2 + P φ = h = const
where E is the total energy. The potential energy can be written as
P φ = F i φ d φ = n 2 a i cos φ l G l 2 m 1 r 1 + m 2 r 2 ,   i = 1 , 2
It follows from the Equations (6) and (7) that the equation of phase trajectories has the form
φ ˙ = ± 2 h P φ .
Figure 2 shows the potential energy (7) and the corresponding phase portrait of the system (2) for the tether length l = 3000 m and the following parameters:
d = 9.4 10 6   m ,   a 1 = 9.38 10 6   m ,   a 2 = 9.42 10 6   m ,   μ = 1.67 10 8 .
The stationary positions for φ π , π can be found from the equation
F i φ * = 0 .
The stable equilibrium positions are φ s = π , 0 , π , and the unstable positions are φ u s = π 2 , π 2 .

2.2. Tether Deflection Angle ψ

To consider this case, let us represent the deflection angle of the tether as
ψ = φ + π .
The equations of motion in polar coordinates for the constant tether length are
ψ ¨ Q i ψ = 0 ,   i = 1 , 2
where
Q i ψ = n 2 sin ψ a i l G sin ψ l m 1 d μ + a i r 1 3 + m 2 d μ 1 + a i r 2 3
r 1 is the distance between the primary 1 and the end mass,
r 1 = a i l cos ψ + d μ 2 + l sin ψ 2 ,
r 2 is the distance between the primary 2 and the end mass,
r 2 = a i l cos ψ + d μ 1 2 + l sin ψ 2 .
Equation (11) has the following energy integral:
E = ψ ˙ 2 2 + P ψ = h = const
where E is the total energy. The potential energy is
P ψ = Q i ψ d ψ = n 2 a i cos ψ l G l 2 m 1 r 1 + m 2 r 2
It follows from the equations (15) and (16) that the equation of phase trajectories has the form
ψ ˙ = ± 2 h P ψ
Figure 3 depicts the potential energy (16) and the corresponding phase portrait of the system (11) for the tether length l = 3000 m and the following parameters:
d = 9.4 10 6   m ,   a 1 = 9.38 10 6   m ,   a 2 = 9.42 10 6   m ,   μ = 1.67 10 8
Equating the generalized force (12) to zero,
Q i ψ * = 0 ,
leads to two types of stationary positions for ψ π , π . The stable equilibrium positions are ψ s = π , 0 , π , and the unstable positions are ψ u s = π 2 , π 2 .

3. Approximate Analytical Solutions

In this section, the approximate analytical solutions of the equations of motion of the tether system in Jacobi elliptic functions [43] for small deflection angles are found and compared with the numerical solutions.

3.1. Tether Deflection Angle φ

The following notation for the small tether deflection angle is used:
φ α .
Let us expand the right hand side of Equation (2), which is an odd periodic function, into a Taylor series and keep the first two terms:
α ¨ = A α + B α 3 ,
where A and B are coefficients depending on the system parameters
A = n 2 a i l + G l m 1 d μ + a i l 1 3 + m 2 d μ 1 + a i l 2 3 ,   i = 1 , 2
B = n 2 a i 6 l + 3 G 2 l m 1 l d μ + a i 2 l 1 5 + m 2 l d μ 1 + a i 2 l 2 5
1 9 m 1 d μ + a i l 1 3 1 9 m 2 d μ 1 + a i l 2 3
l 1 = l + d μ + a i , l 2 = l + d μ 1 + a i
The phase trajectory equation for the Equation (20) in this case is
α ˙ = ± 2 E 1 + A α 2 + B α 4 2
where E 1 = A α 0 2 2 B α 0 4 4 = const and is determined from the initial conditions t = 0 , α ˙ = 0 , α = α 0 , α 0 being the initial tether deflection angle measured from the x-axis. Separating variables in the Equation (23) leads to
d t = ± d α 2 E 1 + A α 2 + B α 4 2
The polynomial under the root of the expression (24) [   2 E 1 + A α 2 + B α 4 2 ] can be factored as
B 2 α c 1 α c 2 α c 3 α c 4
The roots of the polynomial (25) have the form
c 1 , 4 = A N B , c 2 , 3 = A + N B
where N = A 2 4 B E 1 . The right part of expression (24) is an elliptic integral. To reduce it to the canonical form, it is necessary to calculate the modulus, which is determined by [43].
k = z z z + z , 0 < k 2 < 1
where z = c 13 c 24 , z = c 12 c 34 , c i j = c j c i   i = 1 , 2 , 3 , 4 ; j = 1 , 2 , 3 , 4 . Now the expression (24) can be reduced to the form
± z t = F α , k
where z = z + z 2 B 2 , F α , k is the elliptic integral of the first kind,
F α , k = 0 α d α 1 α 2 1 k 2 α 2
Converting the elliptic integral from the expression (28) and using an elliptic sine s n u , k , we obtain the approximate analytical solution:
α t = α 0 s n u , k
where u = z t .
Let us compare the obtained analytical solution Equation (30) with the results of numerical integration of the initial Equation (2). Figure 4 illustrates the simulation results for l = 3000 m and the following initial conditions for the tether system attached in the L 1 libration point:
α 0 = 0.25 rad, α ˙ = 0.00022732 rad/s
and
α 0 = 0.5 rad, α ˙ = 0.00043229 rad/s.
Figure 5 shows the simulation results for l = 3000 m and the following initial conditions for the tether system attached in the L 2 libration point:
α 0 = 0.25 rad, α ˙ = 0.00017327 rad/s
and
α 0 = 0.5 rad, α ˙ = 0.0003377 rad/s.

3.2. Tether Deflection Angle ψ

If the deflection angles are small, one can write, using Equations (10) and (19), that
β = α + π .
In this case, the approximate analytical solution can be written as
β t = β 0 s n u , k ,
where β 0 is the initial deflection angle of the tether measured from the x-axis,
u = z t , k = z z z + z ,   0 < k 2 < 1 , z = z + z 2 D 2 , z = c 13 c 24 , z = c 12 c 34 , c i j = c j c i , i = 1 , 2 , 3 , 4 ; j = 1 , 2 , 3 , 4 c 1 , 4 = C N D , c 2 , 3 = C + N D , N = C 2 4 D E 2 , E 2 = C β 0 2 2 D β 0 4 4 = const , C = n 2 a i l + G l m 1 d μ + a i l 3 3 m 2 d μ 1 + a i l 4 3 , D = n 2 a i 6 l + 3 G 2 l 1 9 m 1 d μ + a i l 3 3 + m 1 l d μ + a i 2 l 3 5 + 1 9 m 2 d μ 1 + a i l 4 3 + m 2 l d μ 1 + a i 2 l 4 5 , l 3 = l d μ a i , l 4 = l d μ 1 a i
Figure 6 shows the simulation results for l = 3000 m and the following initial conditions of the tether system attached in the L 1 libration point:
β 0 = 0.25 rad, β ˙ = 0.00017414 rad/s
and
β 0 = 0.5 rad, β ˙ = 0.00033936 rad/s.
Figure 7 shows the simulation result for l = 3000 m and the following initial conditions for the tether system attached in the L 2 libration point:
β 0 = 0.25 rad, β ˙ = 0.00022753 rad/s
and
β 0 = 0.5 rad, β ˙ = 0.00043273 rad/s.
The comparison of the numerical and analytical results shows that the approximate analytical solutions are able to accurately predict the amplitudes of tether angular oscillations, but not the frequencies. However, in the considered case, the amplitudes are much more important for the analysis of tether oscillations, so the approximate analytical solutions are quite consistent with the numerical solutions.

4. Oscillation Period of the Tether near the Stable Position

4.1. Tether Deflection Angle φ

According to Equation (28) and Ref. [43], the oscillation period of Equation (30) is determined by the formula
τ 1 = 4 K k z
where K π 2 , k = 0 π 2 d α 1 k 2 sin α 2 is the complete elliptic integral of the first kind,
k = z z z + z , 0 < k 2 < 1 , z = z + z 2 B 2 , z = c 13 c 24 , z = c 12 c 34 , c i j = c j c i , i = 1 , 2 , 3 , 4 ; j = 1 , 2 , 3 , 4 c 1 , 4 = A N B , c 2 , 3 = A + N B , A = n 2 a i l + G l m 1 d μ + a i l 1 3 + m 2 d μ 1 + a i l 2 3 , B = n 2 a i 6 l + 3 G 2 l m 1 l d μ + a i 2 l 1 5 + m 2 l   d μ 1 + a i 2 l l 2 5 1 9 m 1 d μ + a i l 1 3 1 9 m 2 d μ 1 + a i l 2 3 , N = A 2 4 B E 1 , E 1 = A α 0 2 2 B α 0 4 4 = const , l 1 = l + d μ + a i , l 2 = l + d μ 1 + a i
Figure 8 shows the dependence of the oscillation period on the length of the tether attached at the L 1 and L 2 libration points based on analytical and numerical calculations for the initial deflection angle of the tether α 0 = 0.25 rad.
The difference between the numerical and analytical calculations is approximately 300 s. As the tether length increases for the case when the tether system is fixed at the libration point L 1 , the period decreases. If the tether is attached at the L 2 libration point, the period of oscillation of the tether decreases at l < 200 m, but begins to increase at l > 200 m. This change in the character of the dependence of the oscillation period on the tether length requires further study. For the tether length l = 3000 m, the analytical calculation gives the oscillation period of about 7000 s ≈ 1.94 h for the tether attached at the L 1 point and 9081 s ≈ 2.52 h for the L 2 point, while the orbital period of Phobos around Mars is 27,540 s ≈ 7.65 h.

4.2. Tether Deflection Angle ψ

According to Ref. [43], the oscillation period of Equation (32) is determined by the formula
τ 2 = 4 K k z
where K π 2 , k = 0 π 2 d β 1 k 2 sin β 2 is the complete elliptic integral of the first kind,
k = z z z + z , 0 < k 2 < 1 , z = z + z 2 D 2 , z = c 13 c 24 , z = c 12 c 34 , c i j = c j c i , i = 1 , 2 , 3 , 4 ; j = 1 , 2 , 3 , 4 c 1 , 4 = C N D , c 2 , 3 = C + N D N = C 2 4 D E 2 , E 2 = C β 0 2 2 D β 0 4 4 = const , C = n 2 a i l + G l m 1 d μ + a i l 3 3 m 2 d μ 1 + a i l 4 3 , D = n 2 a i 6 l + 3 G 2 l 1 9 m 1 d μ + a i l 3 3 + m 1 l d μ + a i 2 l 3 5 + 1 9 m 2 d μ 1 + a i l 4 3 + m 2 l d μ 1 + a i 2 l 4 5 l 3 = l d μ a i , l 4 = l d μ 1 a i
Figure 9 shows the dependence of the oscillation period on the length of the tether attached in the L 1 and L 2 libration points based on analytical and numerical calculations for the initial tether deflection angle β 0 = 0.25 rad.
The difference between the numerical and analytical calculations is approximately 300 s. As the tether length increases for the case when the tether system is fixed at the libration point L 2 , the period rises. If the tether is attached at the L 1 libration point, the period of oscillation of the tether increases at l < 200 m, but begins to decrease at l > 200 m. At the tether length of l = 3000 m, attached in the L 1 or L 2 libration point, the oscillation period according to the analytical calculation is equal to, respectively, 6982 s ≈ 1.94 h и 9028 s ≈ 2.51 h.

5. Tether Tension Force

In this section, the equations of motion of the end mass in the polar reference frame are derived for the angle φ , in order to obtain the tension of the tether. Note that the equations will be the same for the angle ψ . Analytical expressions for the tether tension force are given and the influence of tether system parameters on this force is investigated in both dynamic and static cases.
Consider the equations of planar motion for the end mass in the Oxy coordinate frame within the scope of the classical restricted three-body problem [42]. In the presence of tether tension, these equations can be written as
x ¨ 2 n y ˙ n 2 x = U x 1 m T x ,
y ¨ + 2 n x ˙ n 2 y = U y 1 m T y ,
where T = T x , T y is the tether tension force acting on the end mass from the tether. Position of the body M relative to the origin of coordinates in the polar reference frame l , φ is defined by
x = a i + l cos φ ,   i = 1 , 2 y = l sin φ .
In the polar reference frame l , φ , the Equations (35) and (36) can be written as
φ ¨ + F φ = 0
l ¨ + F l = 1 m T
where T = T x 2 + T y 2 is the magnitude of the tether tension force acting on the end mass from the tether,
F φ = n 2 sin φ a i l G sin φ l m 1 d μ + a i r 1 3 + m 2 d μ 1 + a i r 2 3 + 2 l ˙ l n + φ ˙ ,   i = 1 , 2 ,
F l = n 2 cos φ a i + G m 1 l + cos φ d μ + a i r 1 3 + m 2 l + cos φ d μ 1 + a i r 2 3 l n + φ ˙ 2 ,
r 1 is the distance between the primary 1 and the end mass,
r 1 = a i + l cos φ + d μ 2 + l sin φ 2 ,
r 2 is the distance between the primary 2 and the end mass,
r 2 = a i + l cos φ + d μ 1 2 + l sin φ 2 .

5.1. Static Tension

The static tension of the constant length tether can be calculated using Equation (39) according to the following formula:
T s t = m n 2 l + a i m G m 1 l 1 2 + m 2 l 2 2
where l 1 = l + d μ + a i , l 2 = l + d μ 1 + a i .
Figure 10 shows that the tension force is almost proportional to the length of the tether.
For the 50 kg end mass and the tether length of 3000 m attached in the L 1 libration point, Equation (44) gives the tension force of 0.086 N. When considering the L 2 libration point, the tension force is 0.059 N.

5.2. Dynamic Tension

The dynamic tether tension force can be found using Equations (39) and (41) as
T = m F l = m n 2 cos φ a i + G m 1 l + cos φ d μ + a i r 1 3 + m 2 l + cos φ d μ 1 + a i r 2 3 l n + φ ˙ 2 ,
Figure 11 shows the tether tension force and the tether deflection angle from the gravitational vertical for the end mass of 50 kg and the tether length of 3000 m.
The graphs in Figure 11 allow us to make the following conclusions:
  • The tether is stretched ( T > 0 ) in all cases considered;
  • The greater the amplitude of oscillation of the tether, the greater the period of oscillation.

6. Conclusions

For the tether system attached at the L1 or L2 collinear libration points of the Mars–Phobos system, the equations of motion for the system for the case of massless and non-extensible tether with the end mass have been obtained. The first integrals of these differential equations have been found and used to determine the phase trajectories and the stable equilibrium positions. Simplified equations for small tether deflection angles in Jacobi elliptic functions have been obtained. The oscillation period of the system has been analytically found. It has been shown that the dependencies of the oscillation period on the tether length for L1 and L2 points are different. The obtained approximate analytical solutions and the results of the numerical integration of the original equations of motion for small angles of deflection of the tether are in good agreement. Analytical expressions have been obtained to determine the tether tension, and it has been shown that for the end mass of 50 kg, this force is small and does not exceed 1 N both for the static and dynamic states of the tether.
The results of this study confirm the possibility of a PHLOTE-like mission and give it some theoretical justification. The prospects of using similar tether systems in the future to create a space elevator anchored at the L1 or L2 libration point is a good stimulus for future research, which will also focus on the consideration of an elastic tether.

Author Contributions

Conceptualization: V.S.A.; investigation: V.S.A. and D.V.N.; methodology: V.S.A. and D.V.N.; software: D.V.N.; writing—original draft: D.V.N.; writing—review and editing: V.S.A. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Russian Science Foundation (Project No. 19-19-00085).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Olivieri, L.; Sansone, F.; Duzzi, M.; Francesconi, A. TED Project: Conjugating Technology Development and Educational Activities. Aerospace 2019, 6, 73. [Google Scholar] [CrossRef] [Green Version]
  2. Shi, G.; Zhu, Z.H. Cooperative game-based multi-objective optimization of cargo transportation with floating partial space elevator. Acta Astronaut. 2023, 205, 110–118. [Google Scholar] [CrossRef]
  3. Luo, S.; Cui, N.; Wang, X.; Fan, Y.; Shi, R. Model and Optimization of the Tether for a Segmented Space Elevator. Aerospace 2022, 9, 278. [Google Scholar] [CrossRef]
  4. Aslanov, V.S.; Ledkov, A.S.; Misra, A.K.; Guerman, A.D. Dynamics of Space Elevator After Tether Rupture. J. Guid. Control. Dyn. 2013, 36, 986–992. [Google Scholar] [CrossRef] [Green Version]
  5. Burov, A.A.; Guerman, A.D.; Kosenko, I.I. Tether orientation control for lunar elevator. Celest. Mech. Dyn. Astron. 2014, 120, 337–347. [Google Scholar] [CrossRef]
  6. Weiwei, W.; Zhigang, W.; Jiafu, L. Conceptual Design and Mechanical Analysis of a Lunar Anchored Cislunar Tether. Cosm. Res. 2023, 61, 80–89. [Google Scholar] [CrossRef]
  7. Ziegler, S.W.; Cartmell, M.P. Using Motorized Tethers for Payload Orbital Transfer. J. Spacecr. Rocket. 2001, 38, 904–913. [Google Scholar] [CrossRef]
  8. Aslanov, V.S.; Ledkov, A.S. Swing Principle in Tether-Assisted Return Mission from an Elliptical Orbit. Aerosp. Sci. Technol. 2017, 71, 156–162. [Google Scholar] [CrossRef]
  9. Luo, C.; Wen, H.; Jin, D.; Xu, S. Dynamics of a flexible multi-tethered satellite formation in a Halo orbit with uncertain parameters. Commun. Nonlinear Sci. Numer. Simul. 2021, 99, 105828. [Google Scholar] [CrossRef]
  10. Kempton, K.; Pearson, J.; Levine, E.; Carroll, J.; Amzajerdian, F. Phase 1 Study for the Phobos L1 Operational Tether Experiment (PHLOTE). End Report. NASA 2018, 1, 1–91. Available online: https://ntrs.nasa.gov/citations/20190000916 (accessed on 1 September 2021).
  11. Mashayekhi, M.J.; Misra, A.K. Optimization of tether-assisted asteroid deflection. J. Guid. Control. Dyn. 2014, 37, 898–906. [Google Scholar] [CrossRef]
  12. Baião, M.F.; Stuchi, T.J. Dynamics of tethered satellites in the vicinity of the Lagrangian point L2 of the Earth–Moon system. Astrophys. Space Sci. 2017, 362, 134. [Google Scholar] [CrossRef]
  13. Wong, B.; Patil, R.; Misra, A. Attitude dynamics of rigid bodies in the vicinity of the Lagrangian points. J. Guid. Control. Dyn. 2008, 31, 252–256. [Google Scholar] [CrossRef]
  14. Aslanov, V.S. Dynamics of a Phobos-anchored tether near the L1 libration point. Nonlinear Dyn. 2023, 111, 1269–1283. [Google Scholar] [CrossRef]
  15. Niccolai, L.; Caruso, A.; Quarta, A.A.; Mengali, G. Artificial Collinear Lagrangian Point Maintenance With Electric Solar Wind Sail. IEEE Trans. Aerosp. Electron. Syst. 2020, 56, 4467–4477. [Google Scholar] [CrossRef]
  16. Shan, M.; Shi, L. Comparison of Tethered Post-Capture System Models for Space Debris Removal. Aerospace 2022, 9, 33. [Google Scholar] [CrossRef]
  17. Lv, S.; Zhang, H.; Zhang, Y.; Ning, B.; Qi, R. Design of an Integrated Platform for Active Debris Removal. Aerospace 2022, 9, 339. [Google Scholar] [CrossRef]
  18. Huang, P.; Zhang, F.; Cai, J.; Wang, D.; Meng, Z.; Guo, J. Dexterous tethered space robot: Design, measurement, control, and experiment. IEEE Trans. Aerosp. Electron. Syst. 2017, 53, 1452–1468. [Google Scholar] [CrossRef]
  19. Feng, G.; Zhang, C.; Zhang, H.; Li, W. Theoretical and Experimental Investigation of Geomagnetic Energy Effect for LEO Debris Deorbiting. Aerospace 2022, 9, 511. [Google Scholar] [CrossRef]
  20. Bourabah, D.; Field, L.; Botta, E.M. Estimation of uncooperative space debris inertial parameters after tether capture. Acta Astronaut. 2023, 202, 909–926. [Google Scholar] [CrossRef]
  21. Zhang, F.; Sharf, I.; Misra, A.; Huang, P. On-line estimation of inertia parameters of space debris for its tether-assisted removal. Acta Astronaut. 2015, 107, 150–162. [Google Scholar] [CrossRef]
  22. Aslanov, V.S.; Yudintsev, V.V. Chaos in Tethered Tug–Debris System Induced by Attitude Oscillations of Debris. J. Guid. Control. Dyn. 2019, 42, 1630–1637. [Google Scholar] [CrossRef]
  23. Huang, P.; Zhang, F.; Chen, L.; Meng, Z.; Zhang, Y.; Liu, Z.; Hu, Y. A review of space tether in new applications. Nonlinear Dyn. 2018, 94, 1–19. [Google Scholar] [CrossRef]
  24. Cartmell, M.P.; McKenzie, D.J. A Review of Space Tether Research. Prog. Aerosp. Sci. 2008, 44, 1–21. [Google Scholar] [CrossRef]
  25. Shi, G.; Zhu, Z.H. Adaptive Anti-Saturation Prescribed-Time Control for Payload Retrieval of Tethered Space System. IEEE Trans. Aerosp. Electron. Syst. 2023, 99, 1–11. [Google Scholar] [CrossRef]
  26. Salazar, F.J.T.; Prado, A.B.A. Deployment and Retrieval Missions from Quasi-Periodic and Chaotic States under a Non-Linear Control Law. Symmetry 2022, 14, 1381. [Google Scholar] [CrossRef]
  27. Salazar, F.J.; Prado, A.F. Suppression of Chaotic Motion of Tethered Satellite Systems Using Tether Length Control. J. Guid. Control. Dyn. 2022, 45, 580–586. [Google Scholar] [CrossRef]
  28. Misra, A.K.; Modi, V.J. Deployment and Retrieval of Shuttle Supported Tethered Satellites. J. Guid. Control. Dyn. 1982, 5, 278–285. [Google Scholar] [CrossRef]
  29. Xu, S.; Sun, G.; Ma, Z.; Li, X. Fractional-order fuzzy sliding mode control for the deployment of tethered satellite system under input saturation. IEEE Trans. Aerosp. Electron. Syst. 2018, 55, 747–756. [Google Scholar] [CrossRef]
  30. Ma, Z.; Huang, P. Discrete-Time Sliding Mode Control for Deployment of Tethered Space Robot With Only Length and Angle Measurement. IEEE Trans. Aerosp. Electron. Syst. 2019, 56, 585–596. [Google Scholar] [CrossRef]
  31. Kang, J.; Zhu, Z.H.; Santaguida, L.F. Analytical and Experimental Investigation of Stabilizing Rotating Uncooperative Target by Tethered Space Tug. IEEE Trans. Aerosp. Electron. Syst. 2021, 57, 2426–2437. [Google Scholar] [CrossRef]
  32. Kang, J.; Zhu, Z.H. Passivity-Based Model Predictive Control for Tethered Despin of Massive Space Objects by Small Space Tug. IEEE Trans. Aerosp. Electron. Syst. 2022, 59, 1239–1248. [Google Scholar] [CrossRef]
  33. Ismail, N.A.; Cartmell, M.P. Three dimensional dynamics of a flexible Motorised Momentum Exchange Tether. Acta Astronaut. 2016, 120, 87–102. [Google Scholar] [CrossRef] [Green Version]
  34. Chen, S.; Chen, W.; Chen, T.; Kang, J. In-Plane Libration Suppression of a Two-Segment Tethered Towing System. Aerospace 2023, 10, 286. [Google Scholar] [CrossRef]
  35. Kim, M.; Hall, C.D. Control of a rotating variable-length tether system. J. Guid. Control. Dyn. 2004, 27, 849–858. [Google Scholar] [CrossRef] [Green Version]
  36. Singh, S.; Junkins, J.; Anderson, B.; Taheri, E. Eclipse-conscious transfer to lunar gateway using ephemeris-driven terminal coast arcs. J. Guid. Control. Dyn. 2021, 44, 1972–1988. [Google Scholar] [CrossRef]
  37. Aslanov, V.S. Prospects of a tether system deployed at the L1 libration point. Nonlinear Dyn. 2021, 106, 2021–2033. [Google Scholar] [CrossRef]
  38. Szebehely, V. The Restricted Problem of Three Bodies; Academic Press: New York, NY, USA, 1967; 668p. [Google Scholar]
  39. Koon, W.S.; Lo, M.W.; Marsden, J.E.; Ross, S.D. Dynamical Systems, the Three-Body Problem and Space Mission Design; Springer: New York, NY, USA, 2011; 327p. [Google Scholar]
  40. Schaub, H.; Junkins, J.L. Analytical Mechanics of Space Systems; Aiaa: San Diego, CA, USA, 2003; 819p. [Google Scholar]
  41. Roy, A.E. Orbital Motion; CRC Press: Boca Raton, FL, USA, 2020; 505p. [Google Scholar]
  42. Kluever, C.A. Space Flight Dynamics; John Wiley & Sons: Hoboken, NJ, USA, 2018; 583p. [Google Scholar]
  43. Janke, E.; Emde, F.; Losch, F. Tafeln Hoherer Funktionen; BG Teubner Verlagsgeselschaft: Leipzig, Germany, 1960; 322p. [Google Scholar]
Figure 1. Frame Oxy and the polar frame l , φ .
Figure 1. Frame Oxy and the polar frame l , φ .
Aerospace 10 00541 g001
Figure 2. (a) The potential energy P φ for the tether system attached at the L 1 libration point; (b) the potential energy P φ for the tether system attached at the L 2 libration point; (c) phase trajectories φ ˙ φ corresponding to different levels of the total energy E j j = 1 , 2 , 3 , 4 for the tether system attached at the L 1 libration point; (d) phase trajectories φ ˙ φ corresponding to different levels of the total energy E j j = 1 , 2 , 3 , 4 for the tether system attached at the L 2 libration point.
Figure 2. (a) The potential energy P φ for the tether system attached at the L 1 libration point; (b) the potential energy P φ for the tether system attached at the L 2 libration point; (c) phase trajectories φ ˙ φ corresponding to different levels of the total energy E j j = 1 , 2 , 3 , 4 for the tether system attached at the L 1 libration point; (d) phase trajectories φ ˙ φ corresponding to different levels of the total energy E j j = 1 , 2 , 3 , 4 for the tether system attached at the L 2 libration point.
Aerospace 10 00541 g002aAerospace 10 00541 g002b
Figure 3. (a) The potential energy P ψ for the tether system attached at the L 1 libration point; (b) the potential energy P ψ for the tether system attached at the L 2 libration point; (c) the separatrices ψ ˙ ψ in the phase space corresponding to different levels of the total energy E j j = 1 , 2 , 3 , 4 for the tether system attached at the L 1 libration point; (d) the separatrices ψ ˙ ψ in the phase space corresponding to different levels of the total energy E j j = 1 , 2 , 3 , 4 for the tether system attached at the L 2 libration point.
Figure 3. (a) The potential energy P ψ for the tether system attached at the L 1 libration point; (b) the potential energy P ψ for the tether system attached at the L 2 libration point; (c) the separatrices ψ ˙ ψ in the phase space corresponding to different levels of the total energy E j j = 1 , 2 , 3 , 4 for the tether system attached at the L 1 libration point; (d) the separatrices ψ ˙ ψ in the phase space corresponding to different levels of the total energy E j j = 1 , 2 , 3 , 4 for the tether system attached at the L 2 libration point.
Aerospace 10 00541 g003
Figure 4. Time history of the tether deflection angle for the tether system attached in the L 1 libration point.
Figure 4. Time history of the tether deflection angle for the tether system attached in the L 1 libration point.
Aerospace 10 00541 g004
Figure 5. Time history of the tether deflection angle for the tether system attached in the L 2 libration point.
Figure 5. Time history of the tether deflection angle for the tether system attached in the L 2 libration point.
Aerospace 10 00541 g005
Figure 6. Time history of the tether deflection angle for the tether system attached in the L 1 libration point.
Figure 6. Time history of the tether deflection angle for the tether system attached in the L 1 libration point.
Aerospace 10 00541 g006
Figure 7. Time history of the tether deflection angle for the tether system attached in the L 2 libration point.
Figure 7. Time history of the tether deflection angle for the tether system attached in the L 2 libration point.
Aerospace 10 00541 g007
Figure 8. The oscillation periods of the tether systems attached at the L 1 and L 2 libration points.
Figure 8. The oscillation periods of the tether systems attached at the L 1 and L 2 libration points.
Aerospace 10 00541 g008
Figure 9. The oscillation periods of the tether systems attached at the L 1 and L 2 libration points.
Figure 9. The oscillation periods of the tether systems attached at the L 1 and L 2 libration points.
Aerospace 10 00541 g009
Figure 10. Tension force of the tether, as a function of its length for the end mass of 50 kg: (a) tether attached at the L 1 libration point; (b) tether attached at the L 2 libration point.
Figure 10. Tension force of the tether, as a function of its length for the end mass of 50 kg: (a) tether attached at the L 1 libration point; (b) tether attached at the L 2 libration point.
Aerospace 10 00541 g010
Figure 11. (a) Tension force for the tether attached in the L 1 libration point; (b) tension force for the tether attached in the L 2 libration point; (c) tether deflection angle for the tether attached in the L 1 libration point; (d) tether deflection angle for the tether attached in the L 2 libration point.
Figure 11. (a) Tension force for the tether attached in the L 1 libration point; (b) tension force for the tether attached in the L 2 libration point; (c) tether deflection angle for the tether attached in the L 1 libration point; (d) tether deflection angle for the tether attached in the L 2 libration point.
Aerospace 10 00541 g011
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Aslanov, V.S.; Neryadovskaya, D.V. A Tether System at the L1, L2 Collinear Libration Points of the Mars–Phobos System: Analytical Solutions. Aerospace 2023, 10, 541. https://doi.org/10.3390/aerospace10060541

AMA Style

Aslanov VS, Neryadovskaya DV. A Tether System at the L1, L2 Collinear Libration Points of the Mars–Phobos System: Analytical Solutions. Aerospace. 2023; 10(6):541. https://doi.org/10.3390/aerospace10060541

Chicago/Turabian Style

Aslanov, Vladimir S., and Daria V. Neryadovskaya. 2023. "A Tether System at the L1, L2 Collinear Libration Points of the Mars–Phobos System: Analytical Solutions" Aerospace 10, no. 6: 541. https://doi.org/10.3390/aerospace10060541

APA Style

Aslanov, V. S., & Neryadovskaya, D. V. (2023). A Tether System at the L1, L2 Collinear Libration Points of the Mars–Phobos System: Analytical Solutions. Aerospace, 10(6), 541. https://doi.org/10.3390/aerospace10060541

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