Next Article in Journal
Flutter Optimization of Carbon/Epoxy Plates Based on a Fast Tree Algorithm
Next Article in Special Issue
Consensus SE(3)-Constrained Extended Kalman Filter for Close Proximity Orbital Relative Pose Estimation
Previous Article in Journal
The Bond Graph Modeling and Experimental Verification of a Hydraulic Inertial Vibration Isolator Including Nonlinear Effects
Previous Article in Special Issue
Modeling, Simulation and Control of a Spacecraft: Automated Rendezvous under Positional Constraints
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fuelless On-Orbit Assembly of a Large Space Truss Structure Using Repulsion of the Service Spacecraft by Robotic Manipulators

Keldysh Institute of Applied Mathematics, Russian Academy of Sciences, 125047 Moscow, Russia
*
Author to whom correspondence should be addressed.
Aerospace 2024, 11(8), 635; https://doi.org/10.3390/aerospace11080635
Submission received: 4 July 2024 / Revised: 29 July 2024 / Accepted: 31 July 2024 / Published: 2 August 2024

Abstract

:
A servicing spacecraft motion control approach for the problem of on-orbit truss structure assembly is developed in this paper. It is considered that a cargo container with a rod set and servicing spacecraft are in orbit initially. The assembly procedure is based on spacecraft free-flight motion between the structure’s specified points. The spacecraft is equipped with two robotic manipulators capable of attaching to the structure and holding rods. In addition, the spacecraft can repulse from the structure with a given relative velocity using a manipulator, so the spacecraft and the structure receive impulses. The repulsion velocity vector is calculated in order to reach the structure target point to deliver and install the rod into the truss structure, or to reach the cargo container and take a rod. The problem of searching the repulsion velocity is formulated as an optimization problem with constraints, taking into account the limited value of the repulsion velocity, collision avoidance with structure, restrictions on the angular velocity and translational motion of the structure in the orbital reference frame. This problem is solved numerically with an initial guess vector obtained analytically for simplified motion cases. The application of the proposed control scheme to the assembly of a truss-based antenna is demonstrated. It is shown that the servicing spacecraft is successfully transferred between the structure points by means of manipulator repulsion. Main features and limitations of the assembly problem using a spacecraft with two manipulators are discussed.

1. Introduction

An application of on-orbit assembly for new space systems is expected to become a trend. It allows building large structures in near-Earth orbit to improve the characteristics of scientific, communication, or observation missions. The size of the on-Earth assembled space structure is limited by the dimensions of the launcher payload volume, so many modern spacecrafts enlarge their size by deployable elements. However, any joints are greatly affected by launch vibrations and, after the launch, they can be out of order in orbit. These problems can be overcome by large structure on-orbit assembly; a current state-of-the-art in this field is analyzed in [1]. The space structure components can be delivered to the orbit by several launches and then assembled, increasing the overall size of the resulting space object. An additional advantage discussed in [2] is that the on-orbit assembled orbital structure is not exposed to on-the-ground conditions such as the on-Earth gravitational stress effect and air/oxygen influence on the structure materials.
Telescope or reflector mirrors are not limited in dimensions if their modules are transferred into orbit using several launch vehicles. This concept makes many project planning tasks feasible, as discussed in [3]. Furthermore, some new mission ideas such as telescope reflector elements acting as a solar sail [4] are being developed. As discussed in [1], deploying strategies should be classified based on economic advantage. Different strategies are presented in [5] and evaluated economically in [6]. One of the approaches of on-orbit assembly considers modules capable of docking that are placed by the launchers in close orbits [7,8]. For example, paper [9] proposes a space telescope that can be reconfigured, consisting of two 3U CubeSats capable of autonomous docking, undocking, and re-docking with a small central “15U” class microsatellite. Another approach of on-orbit assembly considers modules to be transferred from cargo container (CC) to the target position of the structure to be assembled [3]. In paper [10], a space solar power plant assembled by such an approach is proposed. The first approach implies onboard propulsion in order to control modules orbital motion, while the CC-based on-orbit assembly approach could rely on fuelless module motion control methods, such as using methods installed on the CC manipulators. The assembling space structure can be truss-based as for the OSAM-2 Mission [11], which involves the on-orbit 3D printing of a lattice truss to deploy solar cells on it. Another method to transfer modules or rods of the truss structure from the CC to the target position is to use a servicing spacecraft (SS). A review on the variety of SS applications is provided in [12].
The problem of the SS’s autonomous control is one of the most challenging in spaceflight dynamics, since it usually implies docking to space objects. In the case of active space debris removal missions, the SS’s relative motion control is aimed at docking to non-cooperative tumbling objects as in [13], which requires taking into account the object chaotic attitude motion. One of the approaches considered in the literature is to use robotic manipulators to perform the docking. The dynamics of spacecraft manipulators for the on-orbit assembly problem have been sufficiently investigated in many papers, particularly in [14,15,16,17,18,19]. So, in this paper, the problem of manipulator control is not considered; it is assumed that this separate task can be solved. A small SS for the on-orbit assembly is considered due to the trend of spacecraft miniaturization for modern space missions. Miniature manipulators [20,21], gas thrusters, and ion thrusters [22,23] are now available for CubeSat-type satellites, providing the necessary element basis for SS’s motion control. Several CubeSat-based in-orbit servicing missions are under development nowadays, for example, the SpEye mission for the demonstration of proximity operations capabilities of a nanosatellite [24].
This paper presents a new approach for SS’s motion control in the problem of truss structure (TS) assembly. The proposed approach is based on the use of manipulators for repulsion and catching, so the motion of SSs between TS points can be considered as free flight. The main advantage of this approach is that it does not require an onboard propulsion system to perform the SS’s transfers; the on-orbit assembly requires renewable energy with solar panels’ electrical power for repulsion and catching by manipulators. The proposed approach takes into account collision avoidance with the TS during the free flight. The SS transfer problem can be considered as boundary problem with restrictions, as in paper [25]. Two numerical search algorithms to solve this problem were applied. Since these algorithms are sensitive to the initial guess vector, simplified motion cases were considered, allowing to analytically calculate the repulsion velocity vector required for transfer.
The rest of the paper has the following structure. In Section 2, the problem statement and main assumptions are formulated. In Section 3, translational and attitude motion equations used in the paper are presented. In order to solve the SS transfer problem, the change in initial conditions due to repulsion is analytically derived in Section 4. The solution of the boundary problem for the planar motion case is provided in Section 5. For the spatial motion case, this problem is to be solved only numerically; two search algorithms are described in Section 6. Furthermore, finally, in Section 7, a demonstration example of antenna-like truss structure assembly is shown.

2. Problem Statement

The objects considered in this paper are the following: a CC with a set of rods, an installing platform placed on the CC, a TS to be assembled using rods and an SS with two manipulators. It is assumed that, using manipulators, the SS is capable of repelling from the CC or TS with the required velocity and to move relative to the structure during the free flight. Another function of the manipulator is to hold one rod and fix it to the TS. These objects are presented in Figure 1.
It is assumed that the common center of mass of the considered objects are in a circular Earth orbit at the initial time. During the assembly process, TS rods must be transferred from their initial position in the CC to their target position (rod installation point) in a predetermined sequence and fixed on the TS. The CC, installation platform, and TS are all considered as a single unified rigid body (UB) during the SS’s free motion. When the SS is fixed to the CC or TS, it is also a part of the UB.
During the free motion, the SS is considered as a point mass, even if the SS holds a rod. The attitude motion of the SS is not considered in this paper. It is assumed that, during the motion, the SS’s center of mass remains outside the UB convex hull, as shown in Figure 2. Additionally, it is assumed that the attitude control system is capable of maintaining the rod perpendicular to the normal position of the UB convex hull, which prevents the collision of the rod with the UB, as shown in Figure 2.
The considered objects’ motion is described using the following reference frames:
  • The structure reference frame (SRF) A x ˜ y ˜ z ˜ is fixed to the CC, point A is the defined point of the CC, and the axes are a right-handed triad. The coordinates of the TS rods’ edge points, of the CC’s center of mass, and of the installing platform’s position are defined in SRF and they are assumed to be constant. The SS’s start point is defined by the radius-vector r 0 S R F in SRF; the SS’s target point r 1 S R F for the SS’s free motion remains constant in SRF during the flight.
  • The body reference frame (BRF) C x ˜ y ˜ z ˜ axes are aligned with the SRF axes, although BRF’s origin point C is located at the UB center of mass.
  • The inertial reference frame (IRF) C E X Y Z origin C E is located in the Earth’s center of mass, the Z-axis is directed along the Earth’s rotation axis, the X-axis is directed towards the vernal equinox, and the Y-axis completes the right-handed triad.
  • The orbital reference frame (ORF) O x y z has its origin O at the system common center of mass at an initial time, and point O moves along the circular orbit with radius R o r b . During the motion, the the system’s common center of mass can deviate from point O. The z-axis of ORF is co-directed with the local vertical direction and the y-axis aligns with the orbital angular momentum; the x-axis completes the right-handed triad.
The transition between these reference frames is described by the following rotation matrices:
A : I R F B R F , U : I R F O R F , S : O R F B R F .
If the SS is located at a certain point r 0 S R F on the UB, it is assumed that the UB and SS move as a single rigid body. It is assumed that a set of initial state parameters, denoted by X , is known before the SS’s repulsion time t 0 . The parameter set X is defined as follows:
X = [ ω U B I R F ( t 0 ) , r U B O R F ( t 0 ) , υ U B O R F ( t 0 ) , A ( t 0 ) , U ( t 0 ) ] T ,
where ω U B I R F is the UB’s angular velocity vector expressed in IRF, r U B O R F is the UB’s center of mass’s radius-vector in ORF, and υ U B O R F is the UB’s center of mass’s velocity vector in ORF. These defined initial conditions and repulsion velocity vector υ 0 S R F expressed in SRF determine the further angular and transnational motion of the UB and free motion trajectory of the SS’s center of mass. Hence, X and υ 0 S R F define the final SS trajectory end point r f S R F at given time t f . The mathematical problem considered in this paper is to find such a repulsion velocity vector υ 0 S R F that leads the SS’s trajectory end point r f S R F to the target point r 1 S R F with a defined acceptable error ε :
υ 0 S R F : Δ r = | r f S R F r 1 S R F | < ε ,
where ε can be defined as the length of the SS’s manipulator capable of catching the final TS point.
This problem statement is formulated for a single SS transfer between two UB points. On-orbit TS assembly can be implemented by a sequence of such SS transfers between the CC’s initial point and rods’ installation points. For each transfer problem to be solved, the UB’s initial conditions are different; during the assembly, the UB’s mass distribution is changing, and additional constraints on the SS’s relative trajectories are to be addressed in order to avoid collision.

3. Equations of Motion

In this paper, the translational motion of the UB’s center of mass and the SS’s center of mass are considered in the ORF. The attitude motion of the UB is also considered relative to ORF. The problem of the SS’s attitude motion can be separated from the translational motion and it is out of scope of this paper.
The motion of the UB or the SS’s center mass’s position vector r O R F = [ x , y , z ] T is described by the Clohessy–Wiltshire Equations [26]:
r ¨ O R F = x ¨ y ¨ z ¨ = 2 ω 0 z ˙ ω 0 2 y 3 ω 0 2 z + 2 ω 0 x ˙ ,
where ω 0 = μ R orb 3 is the orbital angular velocity, μ is the Earth gravitational parameter, and R orb 3 is the circular orbit radius. This system of differential equations has the following solution:
r O R F = C 4 3 C 1 ω 0 ( t t 0 ) + 2 C 2 cos ω 0 ( t t 0 ) 2 C 3 sin ω 0 ( t t 0 ) C 5 sin ω 0 ( t t 0 ) + C 6 cos ω 0 ( t t 0 ) 2 C 1 + C 2 sin ω 0 ( t t 0 ) + C 3 cos ω 0 ( t t 0 ) ,
where C 1 C 6 are the trajectory parameters, which depend on the initial conditions defined by r O R F ( t 0 ) and r ˙ O R F ( t 0 ) as follows:
C 1 = 2 z ( t 0 ) + x ˙ ( t 0 ) ω 0 , C 2 = z ˙ ( t 0 ) ω 0 , C 3 = 3 z ( t 0 ) 2 x ˙ ( t 0 ) ω 0 , C 4 = x ( t 0 ) 2 z ˙ ( t 0 ) ω 0 , C 5 = y ˙ ( t 0 ) ω 0 , C 6 = y ( t 0 ) .
These equations are valid for small distances relative to the ORF origin.
The attitude motion of the UB is described by Euler’s dynamical equations and kinematic relations expressed in attitude quaternions:
J ω ˙ B R F + [ ω B R F × J ω B R F ] = m e x B R F ,
Λ ˙ = 1 2 Λ ω B R F ,
where J is the UB’s inertia tensor calculated in BRF, ω B R F is the absolute angular velocity expressed in BRF, m e x B R F is the external torque, and Λ is the attitude quaternion. In this paper, external torque is a gravitational torque [27]:
m e x = m g r a v = 3 μ | r U B | 5 r U B × J r U B ,
where r U B is the radius-vector from Earth’s center to the UB’s center of mass written in BRF.
During the SS’s free-flight motion and when the SS is fixed at the UB, the inertia tensor J remains constant in SRF and BRF, though at time moments when the SS repels and catches the UB, the J changes instantly due to the changed mass distribution. The tensor J can be calculated as sum of the UB modules’ inertia tensors J k . Rods can be approximated as line segments with the following elements J i j k , of the inertia tensor:
J i j k = r 1 r 2 ( δ i j r 2 r i r j ) d m = m k 0 1 ( δ i j r 2 r i r j ) d τ , r = r 1 k τ + r 2 k ( 1 τ ) .
where r 1 k and r 2 k are the vectors of the kth rod ends in BRF, i , j = 1 , 2 , 3 , m k is the kth rod mass, and δ i j is the Kronecker delta. Since the UB’s center of mass changes its position during the assembly, the inertia tensor is calculated in SRF. According to the Huygens–Steiner theorem, the elements of J are translated into BRF as follows:
J i j B R F = J i j S R F m U B | r U B S R F | 2 δ i j r U B , i S R F r U B , j S R F ,
where m U B is the UB’s mass, and r U B S R F is the vector between the centers of SRF and BRF.

4. Initial Conditions after the SS’s Repulsion

The translational motion Equation (2) and attitude motion Equations (4) and (5) define the SS and the UB’s trajectories if the initial conditions are specified. When the SS repulses from the UB or catches the TS target point, it sets new initial conditions for motion equations and it changes the UB’s mass distribution. The SS repulsion’s velocity vector υ 0 S R F is considered as an input parameter for achieving the required SS trajectory relative to the UB.
It is assumed that the SS’s repulsion and capture using manipulators occurs instantly. The control of the manipulator motion during the repulsion to achieve the desired velocity υ 0 S R F of the SS’s center of mass relative to the UB is not considered in this paper. It is required to find such a repulsion velocity vector υ S S S R F + in order to obtain the desired SS trajectory to get into the neighborhood of the target point r 1 S R F , i.e., to solve the problem (1).
The initial conditions for motion in Equation (2) are to be specified in ORF, though the position of the SS at the time of repulsion and repulsion velocity are defined in SRF. Thus, first, the SS and the UB’s centers of mass after repulsion should be expressed in ORF. Second, the SS’s initial velocity after the repulsion is to be calculated in ORF taking into account the UB’s angular velocity. Furthermore, the UB’s initial velocity taking into account the influence of the SS’s repulsion is to be corrected. Furthermore, finally, the UB’s angular velocity after the repulsion as an initial condition for the attitude motion of Equation (4) should be defined.
In this section, subscripts ( ) and ( ) + denote variables before and after the SS’s repulsion, respectively. Let r U B S R F be the radius-vector of the UB’s center of mass C U B with attached SS, and the radius-vector of the SS before the repulsion is denoted as r S S S R F . After the SS’s repulsion, the UB’s center of mass C U B + radius-vector changes to r U B S R F + . It is assumed that the SS’s radius-vector immediately after the repulsion r U B S R F + is the same as before the repulsion, i.e., r U B S R F = r U B S R F + . In this paper, the SS and UB’s center of mass’s trajectories are expressed in ORF. The radius-vectors of the UB and SS in ORF are defined by the following equations:
r U B O R F + = r U B O R F + S T ( r U B S R F + r U B S R F ) ,
r S S O R F + = r U B O R F + S T ( r S S S R F + r U B S R F ) .
where r U B O R F is the UB’s center of mass defined in ORF before the repulsion. Here, the SS’s position in SRF does not change instantly because of repulsion r S S S R F = r S S S R F + , but the UB’s center of mass is changed r U B S R F r U B S R F + due to the SS’s separation from the UB. The scheme of the system before and after the repulsion is presented in Figure 3.
To calculate the instant change in velocities of the centers of mass due to the SS’s repulsion, it is assumed that, during the short time of repulsion, the external forces do not change the momentum of the whole system, so the momentum conservation law is used. The momentum of the system expressed in ORF before the repulsion and after is as follows:
( m S S + m U B ) υ U B O R F = m S S υ S S O R F + + m U B υ U B O R F + ,
where m S S is the mass of the SS and m U B is the mass of UB without the SS, and υ U B O R F and υ U B O R F + are the UB’s center of mass’s velocity vectors defined in ORF right before the repulsion and immediately after the repulsion, respectively. Since the ORF is a non-inertial reference frame, the effect of the non-inertial forces is also neglected due to the short time period of the repulsion. The velocities of the UB and SS’s centers of mass in ORF is the sum of the instant velocity of the corresponding UB points before the repulsion and relative velocity due to the SS’s repulsion. Since the UB is rotating with an angular velocity vector ω U B O R F before the repulsion, and the UB’s center of mass changes its position, angular velocities υ U B O R F + and υ S S O R F + are defined as follows:
υ U B O R F + = υ U B O R F + [ ω U B O R F × ( r U B O R F + r U B O R F ) ] + Δ υ U B O R F + ,
υ S S O R F + = υ U B O R F + [ ω U B O R F × ( r S S O R F + r U B O R F ) ] + Δ υ S S O R F + ,
where Δ υ U B O R F + is the change in the UB’s center of mass’s velocity due to the SS’s repulsion, and Δ υ S S O R F + is the vector of SS’s repulsion velocity defined in ORF. Substituting (9) and (10) into (8) and taking into account the relation between the vectors of the UB’s center of mass before r U B O R F and after the repulsion r U B O R F + ,
r U B O R F = m S S r S S O R F + + m U B r U B O R F + m S S + m U B ,
the following relation between the change in the UB’s center of mass’s velocity Δ υ U B O R F + and the SS’s repulsion velocity Δ υ S S O R F + can be obtained:
Δ υ U B O R F + = m S S m U B Δ υ S S O R F + .
The problem considered in the paper (1) is formulated to obtain the SS’s repulsion velocity υ S S S R F + in SRF; it can be obtained by the transition of Δ υ S S O R F + into SFR using a transition matrix S as follows:
υ S S S R F + = S Δ υ S S O R F + .
The SS’s repulsion also results in the UB’s angular velocity change. To calculate the UB’s angular velocity after the repulsion ω U B O R F + , it is assumed that the external torques and torques due to non-inertial forces in ORF have no influence on the total angular momentum change of the system during the short repulsion time, so the angular momentum conservation law can be applied. Angular momentum in ORF is the sum of the UB’s angular momentum and the momentum of the center of mass’s impulses relative to the ORF. The following relation using the angular momentum conservation law written in ORF is obtained:
S T J S ω U B O R F + [ r U B O R F × ( m S S + m U B ) υ U B O R F ] = = S T J + S ω U B O R F + + [ r S S O R F + × m S S υ S S O R F + ] + [ r U B O R F + × m U B υ U B O R F + ] ,
where J is the UB’s (with attached SS) inertia tensor before the repulsion, and J + is the UB’s (without SS) inertia tensor after the repulsion. Since, in (13), all the values are written in ORF, the inertia tensor originally calculated in BRF is translated to ORF by S T J S . Thus, the angular velocity after the repulsion ω U B O R F + can be obtained from (13) as follows:
ω U B O R F + = ( S T J + S ) 1 S T J S ω U B O R F + [ r U B O R F × ( m S S + m U B ) υ U B O R F ] [ r S S O R F + × m S S υ S S O R F + ] [ r U B O R F + × m U B υ U B O R F + ] .
Thus, the UB and SS’s center of mass’s positions (6) and (7), the UB and SS’s velocities (9) and (10), and the UB’s angular velocity (14) all depend on the SS’s repulsion velocity υ S S S R F + . These initial conditions after the repulsion define the SS’s trajectory relative to the UB in SRF. So, the problem (1) of finding the SS’s repulsion vector υ S S S R F + to achieve the TS target point r 1 S R F can be formulated as boundary problem.
The change in the motion state parameters after the SS captures the TS target point is calculated using the same momentum and angular momentum conservation laws. The expressions for the UB’s center of mass’s position, velocity, and angular velocity are not provided in the paper for brevity. The initial conditions after the capture depend on the SS’s relative velocity before the capture; it defines the UB’s attitude and translational motion. Thus, UB’s motion parameters are also can be controlled by the choice of the SS’s repulsion velocity υ S S S R F + in the sequential order of the SS’s repulsion and capturing during the assembly process.

5. Planar Motion Case

First, we consider an approach to solve the problem (1) in a simplified motion case for a demonstration of the proposed on-orbit assembly scheme. Furthermore, the results of this section are later used as the initial guess for the numerical search algorithms in the spatial motion case.
The Clohessy–Wiltshire Equation (2) allows the separation the out-of-plane motion along the y-axis from the motion in the orbital plane O x z . In this section, several assumptions are made about the system’s motion:
  • The position vectors of the SS’s center of mass r S S S R F and the UB’s center of mass r U B S R F have zero values for their y-axis components;
  • The angular velocity ω U B O R F of the UB has zero values for its x- and z-axis components;
  • Gravitational torque is not considered;
  • The inertia tensor J is reduced to a single y-axis moment of inertia I y ;
  • The time duration of the SS’s free flight after the repulsion t f t 0 is fixed.
So, such a system’s motion is a plane-parallel motion in orbital plane O x z . The considered planar motion is not a common motion case, it is just a particular example, allowing to simplify the calculations of the repulsion velocity’s first guess for the spatial motion.
In this case, the final deviation of the trajectory point r f S R F = [ x f , 0 , z f ] T at time t f from the target point r 1 S R F = [ x 1 , 0 , z 1 ] T in the TS only has two components Δ x and Δ z :
Δ r S R F = [ Δ x , 0 , Δ z ] T = r 1 S R F r f S R F .
As the result of the integration of motion Equations (2), (4) and (5), the deviation components Δ x and Δ z depend on a set of initial conditions X , the mass of the SS m S S , the mass of the UB m U B , and the moment of inertia I y . This dependence is derived by analytically taking into account the change in the UB’s translational velocity (12) and change in the UB’s angular velocity (14). In the planar motion case, the SS’s repulsion velocity vector υ S S S R F + = [ υ 0 x , 0 , υ 0 z ] T also has only two nonzero components. The dependence of Δ x and Δ z on υ 0 x and υ 0 z has the following form:
Δ x = x 1 + A 1 cos B A 2 sin B , Δ z = z 1 + A 1 sin B + A 2 cos B ,
where
A 1 = a 1 υ 0 x + a 2 υ 0 z + a 3 , A 2 = a 4 υ 0 x + a 5 υ 0 z + a 6 , B = b 1 υ 0 x + b 2 υ 0 z + b 3 .
Here, a set of coefficients a i and b i are calculated using the parameter set X and mass-inertial parameters m S S , m U B , and I y . Analytical expressions of a i and b i are omitted from the text since they are too bulky. The value of the parameter B is a final increment of the UB’s angle after the SS’s free relative motion between the given time of repulsion t 0 and time of the SS attaching to the UB t f :
B = t 0 t f ω y d t ,
where ω y is an angular velocity of the UB after the SS’s repulsion, which is assumed to be constant since no external torques are taken into account in this section.
For the simplified case of motion considered, the problem (1) can be reformulated to find the exact solution that implies a zero-valued final trajectory deviation:
{ υ 0 x , υ 0 z } : Δ x = 0 , Δ z = 0 .
Even for this simplified case, the problem (16) cannot be solved in explicit form because Equation (15) is nonlinear with respect to the repulsion velocity components υ 0 x and υ 0 z . So, the solution of problem (16) is to be found numerically. However, using additional assumptions, the approximate solutions can be found analytically.
For example, if it is assumed that the UB’s rotation angle is small during the SS’s free motion, the value of parameter B can be considered as zero in (15). Equation (15) is simplified to the following equation:
Δ x A 1 + x 1 , Δ z A 2 + z 1 .
In this case, the approximate solution of the (16) is as follows:
υ 0 x a 2 ( a 6 + z 1 ) a 5 ( a 3 + x 1 ) a 1 a 5 a 2 a 4 , υ 0 z a 4 ( a 3 + x 1 ) a 1 ( a 6 + z 1 ) a 1 a 5 a 2 a 4 .
If coefficients a i are significantly larger than coordinates x 1 and z 1 , the approximate solution can be found from A 1 = 0 and A 2 = 0 and it can be further simplified:
υ 0 x a 2 a 6 a 3 a 5 a 1 a 5 a 2 a 4 , υ 0 z a 3 a 4 a 1 a 6 a 1 a 5 a 2 a 4 .
A particular example of a parameter set X of initial conditions is demonstrated in Figure 4. The value of the separation velocity components υ 0 x and υ 0 z is limited by the 5 cm/s in this example. This limitation of velocity components is based on the next simple evaluation. If it is supposed that the manipulator produces a force of F = 0.1 N along an acceleration trace of L = 50 cm, and pushes the SS loaded with a rod with a total weight of m = 35 kg, then the velocity of the repulsion can be calculated according to the following expression obtained under the assumption of constant acceleration:
υ = 2 L F m 5.3 cm / s .
The maximum length of the manipulator is defined as 1 m, so it is assumed that the acceptable final trajectory deviation must be less than this value. Using the dependence (15) of the trajectory deviation, components Δ x and Δ z on the separation velocity from the value of the final position deviation Δ r = | Δ r | can be calculated on a grid of υ 0 x , υ 0 z [−5; 5] cm/s for the specified parameters. Figure 4 shows a heatmap of Δ r and points of numerically obtained grid solution (with a resolution of 0.01 cm/s) and approximate solutions calculated using (17) and (18). For this particular case, the error of the repulsion velocity guesses are about 1 cm/s for (17) and about 1.5 cm/s for (18).
For a numerical study of the accuracy of approximate solutions (17) and (18), the random initial conditions and random initial and final points are considered; its standard deviation values are presented in the right side of Table 1. The resulting values of parameters a i , b i , x 1 , and z 1 for random initial conditions are presented in Figure 5. Particularly, it can be seen that the assumption used for solution (18) that the coefficients a i are significantly larger than coordinates x 1 and z 1 is correct. For each parameter values, the final deviation value Δ r is calculated. The results of the calculations are presented as boxplots in Figure 6. Here, the boxplot (a) shows the error of final position calculation and the boxplot (b) provides the deviation of the repulsion velocity calculated by the approximate solution from the grid search solution. Inside the box, there are 50% of the simulations results, below and under the box are 25%, and the red line is a mean value. The obtained accuracy of several meters is unacceptable to use it for the SS’s repulsion velocity calculation. However, this accuracy is considered acceptable as the initial guess for the numerical search of the solution (1) for the spatial motion case. As it can be concluded from Figure 4, the maximum value of the final position error can be up to 1200 m, and guess solutions (18) and (17) reduce this error by 10 2 times.
In the case where the SS’s repulsion velocity is limited by some feasible value, which can be implemented by the SS’s manipulator, no solution could be found to the problem (1). It can be demonstrated in the planar motion case. Distributions of the final position deviation Δ r in the repulsion velocity domain { υ 0 x , υ 0 z } [ 5 ; 5 ] cm/s for three values of UB’s initial angular velocities are presented in Figure 7. The left two cases in Figure 7 with UB’s initial angular velocities of 10 4 rad/s and 10 3 rad/s are solutions of the problem (16), while the right figure presents the case with an initial angular velocity of 7 × 10 3 rad/s, where the minimum of Δ r corresponds to deviation of 40 m, which is unacceptable and the problem has no solution in this case.

6. Assembly in Spatial Motion Case

In a common spatial motion case, it is impossible to obtain the final trajectory deviation Δ r in an analytical form as for the planar motion case (15), because the dynamical Euler Equation (4) has no analytical solution except for the three known attitude motion cases—the Euler, Lagrange, and Kovalevskaya cases. Since, in this study, the gravity-gradient torque is taken into account in the UB’s attitude motion equation, and the inertia tensor is arbitrary and it changes during the assembly, it does not meet requirements of any of these three special attitude motion cases. Thus, the Δ r dependence on the SS’s repulsion velocity υ S S S R F + can be obtained only by the numerical integration of the attitude and translational motion equations.

6.1. Constrained-Motion Parameters

In common motion case, the following set of additional requirements and constraints for the problem (1) is considered for taking into account the practical aspects of TS assembly:
  • The SS’s trajectory must not intersect the UB’s elements to avoid collision with them except for the target point r 1 S R F .
  • Since, at a high UB angular velocity, the problem (1) could have no solution in the limited domain of the repulsion velocity values as it has been demonstrated in the planar motion case, the value of the UB’s angular velocity ω U B O R F after the TS’s capture should not exceed an acceptable value.
  • To maintain the UB’s attitude in ORF in the vicinity of its initial attitude during the assembly, it can be additionally demanded that the UB’s attitude deviation after the TS’s capture should be inside the limited domain. It is reasonable to place initially the UB in the vicinity of the stable gravitational equilibrium attitude position and stay in this vicinity during the assembly.
  • During the assembly, the UB’s center of mass can significantly deviate from the origin of the ORF moving along the circular orbit. In order to prevent possibly undesirable deviation from the initial orbit, the UB’s center of mass’s position r U B O R F and velocity υ U B O R F must be kept within prescribed limits.
In order to address these limitations during the assembly, in this work, the following vector c of the UB’s motion-constrained parameters is taken into account in the numerical search of the problem (1)’s solution:
c = | r U B O R F | , | υ U B O R F | , | ω U B O R F | T .
To meet the requirements of the SS’s relative trajectory to prevent the collision with the TS, the following approaches are proposed:
  • Approach 1. During the numerical search of the repulsion velocity, the final position deviation Δ r is calculated using the integration of the attitude and translational motion equations. At each integration point of the SS’s center of mass’s position, the minimum distance with the UB’s points approximated by a grid is calculated. If this distance is less than the defined characteristic size of the SS, it is assumed that the trajectory intersects the UB. At this point r c S R F , the integration of the motion equations stops and the deviation is calculated as Δ r = r c S R F r 1 S R F . If | Δ r | > ε ; then, it is assumed that the trajectory collides with UB’s elements and the numerical search continues. If | Δ r | < ε , then it is assumed that solution of the problem (1) is found.
  • Approach 2. Another approach is to consider a continuous function g to describe the trajectory intersection with the UB’s elements. Function g is used to define the hyperquadratic potential field in the vicinity of the structure; this approach for structure assembly was proposed in [28]. Function g is positive if the SS is outside the UB’s elements, and negative if the SS is inside. Since the UB consists of a container and a set of rods, which can be considered as elongated bodies, the following function g is used to describe the potential field for a cylinder:
    g = d 2 ( ρ y 2 + ρ z 2 1 ) , 1 < ρ x < 1 , l 2 ( | ρ x | 1 ) , e l s e i f ρ y 2 + ρ z 2 < 1 , d 2 ( ρ y 2 + ρ z 2 1 ) + l 2 ( | ρ x | 1 ) , e l s e ,
    where d is a rod diameter, l is the rod length, and ρ x , ρ y , and ρ z , are dimensionless projections of the position vector r S S onto the rod-fixed reference frame. A schematic heatmap of the g function for a cylinder is shown in Figure 8. In this example, the length of the cylinder is 2 m and its diameter is 1 m; the boundaries of the rod are represented by blue lines in Figure 8. With defined UB element positions in SRF, a sum of the functions g i corresponding to each element forms a final potential field of the UB. If the SS’s trajectory does not intersect the UB, each function g i must remain strictly positive or be higher than some positive value to prevent the collision. These functions can also be included in the vector of constrained parameters c in (19) for a numerical search of the SS’s repulsion velocity.

6.2. Target Function for the Numerical Search

For a numerical search of the SS’s repulsion velocity υ 0 = υ S S S R F + taking into account constrained-motion parameters, the original problem (1) is reformulated using the following target function Φ ( υ 0 ) :
υ 0 : Φ ( υ 0 ) = Δ r 2 m i n , c i m a x c i 0 , i = 1 N ,
where υ 0 D , D R 3 is a feasible region of velocity vector υ 0 , Δ r is final trajectory deviation vector Δ r = r 1 S R F r f S R F , c = [ c 1 , c N ] T is vector of constrained variables, and c m a x = [ c 1 m a x , c N m a x ] T is the vector of maximal acceptable values of the constrained parameters. The problem (21) can then be transformed into an unconditional problem using the barrier or penalty function introduced in [29]:
Φ b a r r i e r ( υ 0 , c ) = Φ ( υ 0 ) μ i = 1 N l o g 1 c i c i m a x ,
Φ p e n a l t y ( υ 0 , c ) = Φ ( υ 0 ) μ i = 1 N 1 c i c i m a x + ,
where μ > 0 is a small positive parameter and ( ) + represents the rectified linear unit (ReLU) function. However, the Φ b a r r i e r function is undefined if the constraints are not fulfilled, which may cause the non-convergence during the minimization of the function. The Φ p e n a l t y function highly depends on parameter μ and it only counts constraints in the constraint-violation area. These features may lead to a local minimum of the function that is not the solution of the problem (1). Thus, a new function is required that is characterized by the given functions’ advantages, but does not have their shortcomings. That is why, in this paper, the following barrier function is proposed:
Φ ( υ 0 , c ) = Φ ( υ 0 ) μ i = 1 m l o g ( ε i ) , ε i = 1 c i c i m a x + δ c i c i m a x , ( 1 c i c i m a x ) 0 , δ e 1 c i c i m a x , ( 1 c i c i m a x ) < 0 ,
where δ is the hyperparameter that connects the function Φ ( υ 0 ) defined in the constraint-fulfilled domain with the function Φ ( υ 0 ) defined in the constraint-violated domain. In the case where δ 1 , the function Φ ( υ 0 ) represents a penalty function in the constraint-violated domain. An example of the logarithm argument function ε ( c i c i m a x ) and the function Φ ( c i c i m a x ) dependence on the ratio c i c i m a x using different values of δ is shown in Figure 9. From the figure, it can be concluded that the function Φ ( c i c i m a x ) is smooth and the gradient minimization methods can be applied for the numerical solution search.

6.3. Repulsion Velocity’s Initial Guess

The performance of many numerical minimization methods strongly depends on the initial guess. Gradient-based methods are prone to converge to a local minimum when searching the parameter area where the initial guess vector is defined. So, this aspect of numerical search is of high importance in this work. This paper proposes two approaches to calculate an initial guess for the SS’s repulsion velocity υ 0 using the following simplified motion cases:
(1)
The planar motion case described in Section 5 is used for the initial guess υ 0 = [ υ 0 x , 0 , υ 0 z ] T ’s calculation according to (18).
(2)
In the spatial motion case, under the assumption that the UB’s position and attitude are fixed in ORF, the following initial guess υ 0 = [ υ 0 x , υ 0 y , υ 0 z ] T using the solution (3) of the Clohessy–Wiltshire Equation (2) is calculated according to [25]:
υ 0 x = ω 0 ( 6 Δ t ω 0 z 0 + x 1 x 0 ) sin ( Δ t ω 0 ) + ( 14 z 0 2 z 1 ) ( cos ( Δ t ω 0 ) 1 ) 3 Δ t ω 0 sin ( ω 0 Δ t ) 8 cos ( ω 0 Δ t ) + 8 , υ 0 y = ω 0 ( y 0 cos ( ω 0 Δ t ) + y 1 ) sin ( ω 0 Δ t ) , υ 0 z = ω 0 ( 3 Δ t ω 0 z 0 + 2 x 1 2 x 0 ) cos ( Δ t ω 0 ) 3 Δ t ω 0 z 1 + 2 x 0 2 x 1 + 4 ( z 1 z 0 ) sin ( Δ t ω 0 ) 3 Δ t ω 0 sin ( ω 0 Δ t ) 8 cos ( ω 0 Δ t ) + 8 ,
where r 0 O R F = [ x 0 , y 0 , z 0 ] T is the SS’s start point, r 1 O R F = [ x 1 , y 1 , z 1 ] T is the SS’s target point, and Δ t = t t 0 is defined time of the SS’s free flight. It can be concluded from (23) that the selection of the flight time Δ t strongly affects the SS’s repulsion velocity and, as a sequence, it affects the SS’s trajectory. Furthermore, it follows from (23) that the flight time multiplication to half of the orbital period Δ t = π k / ω 0 , k = 1 , 2 , 3 , , n is unacceptable; it leads to infinite repulsion velocity and it means that there is no solution to the boundary problem. In this paper, for an acceptable initial guess search, a set of SS’s flight times is considered Δ t = [ Δ t 1 , , Δ t k ] T ; for each Δ t i , the repulsion velocity Υ = [ υ 0 1 , υ 0 k ] is calculated using (23) and the trajectory collision avoidance requirement is roughly checked. It is assumed that the forbidden zone for the trajectory is a sphere with a diameter with endpoints r 0 O R F and r 1 O R F . The SS’s trajectory with repulsion velocity υ 0 i can be roughly considered as a set of SS’s positions r j i with m points. The repulsion velocity υ 0 i is specified as being invalid if the following inequality is true:
j : r j i 1 2 ( r 0 O R F + r 1 O R F ) < 1 2 r 0 O R F r 1 O R F , j = 1 , , m .
This condition corresponds to the case when the SS’s radius-vector from the sphere center is less than the sphere radius. Using (24), the repulsion velocity for the acceptable initial guess in terms of preliminary rough collision avoidance requirements and corresponding flight time Δ t is selected.

6.4. Numerical Search Algorithms

Two numerical search algorithms were used in this paper to solve problem (21) depending on the proposed above approaches to address the SS’s collision avoidance aspect:
  • Algorithm 1. For the case of using Approach 1 from Section 6.1, when the numerical integration stops when the SS’s trajectory intersects the UB’s element during the search for the solution, and when the vector of constrained parameters c from (19) is short (it does not include the g functions), then Newton’s method was used in this paper to solve the problem (21). Newton’s method iteratively calculates the repulsion velocity using the following formula [30]:
    υ 0 i + 1 = υ 0 i Φ υ 0 | υ 0 = υ 0 i 1 Φ | υ 0 = υ 0 i ,
    where i is the iteration number of the method, and Φ is the gradient vector of the target function (22):
    Φ = Φ υ 0 x , Φ υ 0 y , Φ υ 0 z T .
    The iterative procedure stops if the norm of the difference of two sequential values of υ 0 i + 1 is less than the defined small value ε υ 0 : υ 0 i + 1 υ 0 i < ε υ 0 .
  • Algorithm 2. For the case of using the Approach 2 from Section 6.1, when the functions g (20) are added to the vector of constrained parameters c to address the collision avoidance, the trust region [30] method from the Python library “scipy” was used. This method is based on a subset of the objective function region that is approximated using a model quadratic function. If an appropriate approximation of the objective function is found within the trust region, then the searching region is expanded; conversely, if the approximation is poor, then the region is contracted. The criterion for appropriate approximation is the improvement of the objective function value.
Both algorithms were used for the numerical search for the problem’s solution. Algorithm 1 is faster than Algorithm 2, though Algorithm 1 is more sensitive to the initial velocity guess. The main reason that Algorithm 2 runs significantly slower is that the solution for the problem Φ ( υ 0 ) m i n (21) is searched on the boundary of fulfilling constraints c i m a x c i 0 . So, the more computationally complex Algorithm 2 is used for the repulsion velocity search if Algorithm 1 fails to find a valid solution.

7. On-Orbit Assembly Example

We consider an example of the large truss-based structure on-orbit assembly using the proposed approach for fuelless control of the servicing spacecraft motion. As described in Section 2, it is assumed that, initially, in orbit, there is a CC with a set of rods, installing a platform placed on the CC; the TS is to be assembled using rods and an SS with two manipulators. The structure to be assembled is the antenna-like construction presented in Figure 10.
The main parameters of the considered system and parameters of the algorithms are provided in Table 2. The initial circular orbit height of the UB is 800 km, which corresponds to an orbital period of about 5700 s. It is assumed that the SS is initially attached to the CC. During the assembly process, the motion of the system is influenced by the central gravity field and by the repulsion and capturing of the SS. The effects of any orbital disturbances such as the second harmonic of Earth’s gravity field, aerodynamic forces, solar pressure, etc., are neglected. Initially, the UB is stabilized in ORF; the axes of BRF coincide with ORF’s axes.

7.1. Acceptable Repulsion Velocity Restriction

For the proposed assembly method, it is the most important to estimate the required maximum value for the SS’s repulsion velocity to perform the rod transfers. This value can be estimated using the simplified motion cases described above. In Equation (23), for the calculation of the initial guess for the repulsion velocity, the x - and z-components of repulsion velocity have same denominator. Furthermore, the same denominator is applied to the x - and z-components of repulsion velocity in Equation (18) for the planar motion case. Values of the denominator (23) for the x - and z-components depending on transfer time are represented in Figure 11. Transfer time regions’ near-zero value of the denominator requires an extremely large repulsion velocity and should be avoided by choosing an appropriate transfer time. Additionally, the y-component denominator (23) is zero during one-half of the orbital period. Thus, transfer time regions with high denominator values should be chosen.
Examples of the required repulsion velocity values (23) and (18) according to simplified cases for randomly generated one hundred transfers between structure points are shown in Figure 12 and Figure 13. Blue boxplots correspond to formulas (23), and yellow boxplots correspond to formulas (18). Half of the results are inside the box and the horizontal line is the median value of the required repulsion velocity. Figure 12 demonstrates the repulsion velocities’ boxplots for transfers between points with a structure of an overall size of L = 20 m, which corresponds to the structure’s overall size of the antenna-like structure considered in this section. For comparison, Figure 13 presents the repulsion velocities’ boxplots for the transfers with a structure of an overall size of L = 100 m. It can be noticed that the short transfer time when the denominator is small results in higher values of the repulsion velocities. In the case of a time transfer of 3000 s, the values of the repulsion velocities are also high because this is close to the half of the orbital period and close to the zero-value of the denominator for the y-component’s velocity in (23). For the other cases of time transfer, the required repulsion velocity is less than the value of 5 cm/s for the structure with a characteristic size of L = 20 m, and 20 cm/s for L = 100 m.
Thus, in this example of space structure assembly, the repulsion velocity is restricted by the reasonable value of 5 cm/s, though, in practice, this restriction is the result of the choice of a particular SS’s manipulators and its parameters.

7.2. Simulation of the Assembly Process

The assembly process consists of a sequence of predetermined transfers of the SS between the UB’s points. Initially, the SS is located at the CC, it takes the rod, and then repulses from the UB using the manipulator. When the SS’s trajectory is in the vicinity of the target point, the SS’s manipulators capture the UB; then, the SS installs the rod to the structure and repulses by manipulator to transfer to the initial point at the CC to take the next rod. Then, the cycle is repeated until the structure is fully assembled.
Figure 14 presents an example of the SS’s position deviation from the target point during the SS’s motion (top plot) and the values of the constrained-motion parameters (19) divided by their maximum values c i / c i m a x (bottom plot). The vector of the constrained-motion parameters includes the value of the UB’s angular velocity ω U B O R F , and the UB’s center of mass’s position r U B O R F and velocity υ U B O R F ; their maximum values are provided in Table 2. In Figure 14, red vertical dashed lines indicates time points of the SS’s repulsion and catching by the UB. The catching is performed when the SS’s position deviation is less than the length of the manipulator capture radius of 0.5 m. The SS and the UB’s free-flight motion takes place between the repulsion and catching points. As can be seen from Figure 14, at the time moments of the SS’s repulsion from the UB and catching by the UB, the UB’s angular and translational velocities change simultaneously according to (12) and (14). The change in the UB’s angular velocity is caused by oscillations of the UB in the vicinity of the gravity-stable orientation.
The repulsion velocities for the presented transfers are found using the two algorithms described above of numerical search. Figure 15 presents an example of the convergence process of reducing the SS’s final position deviation vs. the number of required iterations in the case where Algorithm 1 is used to solve the problem using an initial velocity guess calculated for simplified cases. Algorithm 1 stops when the Δ r ( t ) is less than the manipulator’s capture radius of 0.5 m, or when the difference between the two values of the repulsion velocities υ 0 i + 1 υ 0 i < ε υ 0 is less than ε υ 0 = 10 4 m/s. In most cases presented in Figure 15, Algorithm 1 successfully finds the repulsion velocity with Δ r ( t ) < 0.5 m taking into account the collision avoidance aspect according to Approach 1 described in Section 6.1. However, in some cases, Algorithm 1 fails to find a solution as is shown by the red line in Figure 15: the Δ r ( t ) is more than the manipulator’s capture radius and the repulsion velocity is converged. In that case, Algorithm 2 is used to solve the problem (21). An example of the convergence process of Algorithm 2 is presented in Figure 16. It requires many more iterations, though, since it uses the more complicated Approach 2 to take into account the collision avoidance with the structure, it finds solution to the transfer problem more steadily. It can be noted from Figure 16 that the position error decreases insignificantly from 10 to 45 iterations. It can be explained by the slow convergence rate of the trust region method in the area close to the violation of condition g i 0 (20) related to structure collision avoidance.
Figure 17 presents the SS’s position deviation from the target points during several transfers. Furthermore, in Figure 18, the relative constrained-motion parameters during multiple transfers are shown. From Figure 18, it can be seen that the relative constrained values do not exceed the value of 1, so the relative constrained values are not violated. However, the SS’s repulsion impulses and the change of the UB’s mass distribution may lead to the significant deviation of the UB from the gravity-stabilized attitude and result in an increase of the angular velocity.
Intermediate assembly stages are shown in Figure 19 and Figure 20. In these figures, the CC, the SS, trajectories of the SS, and a set of the installed rods are presented. In Figure 19, most of the rods are still inside the CC, and the symmetry axis of the cylindrical CC with a minimal inertia moment is almost aligned with the local vertical axis z of ORF. So, the UB’s attitude is in the vicinity of the gravity-gradient stabilization equilibrium. After a while, the mass distribution is changed, and the attitude of the UB due to SS’s repulsion is also changed significantly according to Figure 20.
Figure 21 shows the values of the calculated SS’s repulsion velocities during the assembly process. It can be concluded that, due to slightly increasing UB’s angular velocity, the maximal values of the SS’s repulsion velocity also steadily increase, though they do not exceed the limit of 5 cm/s. It means that the estimations of the maximum repulsion velocity obtained using the simplified motion cases (see Figure 12) are valid for the considered example.

8. Conclusions

This paper proposed a new approach to the problem of the space structure construction using the servicing spacecraft with two manipulators. With these manipulators, the servicing spacecraft is able to move between the structure points using repulsion and catching. In the case of planar motion, the boundary problem was solved analytically using simplifying assumptions. In the spatial motion case, the repulsive velocity was searched for using numerical minimization techniques, taking into account a set of restrictions on the motion parameters and considering collision avoidance of the servicing spacecraft with the space structure. A numerical example demonstrated the successful construction of the particular antenna-like truss structure. However, it should be noted that, in the case of a complex truss-based structure or at a high angular velocity of the structure, the acceptable repulsion velocity for the rod transfers could not be found in a feasible region of velocity vector. In this case, the servicing spacecraft would need to be equipped with onboard propulsion and the repulsion by manipulator can be considered as a means for reducing the fuel consumption. This case will be considered in further studies by the authors.

Author Contributions

Software, validation, and writing—original draft, V.O.; formal analysis, U.M.; methodology and project administration, M.O.; supervision, conceptualization and writing—review and editing, D.I. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Moscow Center of Fundamental and Applied Mathematics, Agreement with the Ministry of Science and Higher Education of the Russian Federation, No. 075-15-2022-283.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
IRFInertial reference frame
ORFOrbital reference frame
BRFBody reference frame
SRFStructure reference frame
SSServicing spacecraft
TSTruss structure
CCCargo container
UBUnited body

References

  1. Piskorz, D.; Jones, L.K. On-orbit assembly of space assets: A path to affordable and adaptable space infrastructure. Aerosp. Corp. 2018, 1–10. Available online: https://aerospace.org/sites/default/files/2018-05/OnOrbitAssembly_0.pdf (accessed on 30 July 2024).
  2. Davis, J.P.; Mayberry, J.P.; Penn, J.P. On-orbit servicing: Inspection repair refuel upgrade and assembly of satellites in space. Aerosp. Corp. 2019, 1–14. Available online: https://aerospace.org/sites/default/files/2019-05/Davis-Mayberry-Penn_OOS_04242019.pdf (accessed on 30 July 2024).
  3. Oegerle, W.R.; Purves, L.R.; Budinoff, J.G.; Moe, R.V.; Carnahan, T.M.; Evans, D.C.; Kim, C.K. Concept for a large scalable space telescope: In-space assembly. Space Telesc. Instrum. I Opt. Infrared Millim. 2006, 6265, 755–766. [Google Scholar]
  4. Izzo, D.; Pettazzi, L.; Ayre, M. Mission Concept for Autonomous on Orbit Assembly of a Large Reflector in Space. In Proceedings of the 56th International Astronautical Congress, Fukuoka, Japan, 17–21 October 2005. IAC-05-D1.4.03. [Google Scholar]
  5. Gralla, E.L.; De Weck, O.L. Strategies for on-orbit assembly of modular spacecraft. J. Br. Interplanet. Soc. 2007, 60, 219–227. [Google Scholar]
  6. Sternberg, D.; Chodas, M.; Jewison, C.; Jones, M.; De Weck, O. Multidisciplinary system design optimization of on orbit satellite assembly architectures. In Proceedings of the 2015 IEEE Aerospace Conference, Big Sky, MT, USA, 7–14 March 2015. [Google Scholar]
  7. McInnes, C.R. Distributed control of maneuvering vehicles for on-orbit assembly. J. Guid. Control Dyn. 1995, 18, 1204–1206. [Google Scholar] [CrossRef]
  8. Izzo, D.; Pettazzi, L. Equilibrium shaping: Distributed motion planning for satellite swarm. In Proceedings of the 8th International Symposium on Artificial Intelligence, Robotics and Automation in Space, Munich, Germany, 5–8 September 2005. [Google Scholar]
  9. Underwood, C.; Pellegrino, S.; Lappas, V.J.; Bridges, C.P.; Baker, J. Using CubeSat/micro-satellite technology to demonstrate the Autonomous Assembly of a Reconfigurable Space Telescope (AAReST). Acta Astronaut. 2015, 114, 112–122. [Google Scholar] [CrossRef]
  10. Cheng, Z.; Hou, X.; Zhang, X.; Zhou, L.; Guo, J.; Song, C. On-orbit assembly mission for the space solar power station. Acta Astronaut. 2016, 129, 299–308. [Google Scholar] [CrossRef]
  11. On-Orbit Servicing, Assembly, and Manufacturing 2 (OSAM-2). Available online: https://www.nasa.gov/mission/on-orbit-servicing-assembly-and-manufacturing-2-osam-2/ (accessed on 2 November 2023).
  12. Li, W.J.; Cheng, D.Y.; Liu, X.G.; Wang, Y.B.; Shi, W.H.; Tang, Z.X.; Gao, F.; Zeng, F.M.; Chai, H.Y.; Luo, W.B.; et al. On-orbit service (OOS) of spacecraft: A review of engineering developments. Prog. Aerosp. Sci. 2019, 108, 32–120. [Google Scholar] [CrossRef]
  13. Akhloumadi, M.; Ivanov, D. Satellite Motion Control System Parameters on Performance of Space Debris Capturing. Aerospace 2020, 7, 160. [Google Scholar] [CrossRef]
  14. Rutkovskii, V.Y.; Sukhanov, V.M.; Glumov, V.M. Motion equations and control of the free-flying space manipulator in the reconfiguration mode. Autom. Remote Control 2010, 71, 70–86. [Google Scholar] [CrossRef]
  15. Glumov, V.M.; Krutova, I.N.; Sukhanov, V.M. Some features of powered gyrostabilization of a large space structure assembled in orbit. Autom. Remote Control 2018, 79, 524–534. [Google Scholar] [CrossRef]
  16. Glumov, V.M.; Rutkovskii, V.Y.; Sukhanov, V.M. Formation of the control strategy for the large space structure assembled in orbit. Autom. Remote Control 2007, 68, 2113–2127. [Google Scholar] [CrossRef]
  17. Uzo-Okoro, E.; Erkel, D.; Manandhar, P.; Dahl, M.; Kiley, E.; Cahoy, K.; De Weck, O.L. Optimization of On-Orbit Robotic Assembly of Small Satellites. In Proceedings of the Accelerating Space Commerce, Exploration, and New Discovery Conference, Las Vegas, NV, USA, 16–18 November 2020; American Institute of Aeronautics and Astronautics: Reston, VA, USA, 2020. [Google Scholar]
  18. Pisculli, A.; Felicetti, L.; Gasbarri, P.; Palmerini, G.B.; Sabatini, M. A reaction-null/Jacobian transpose control strategy with gravity gradient compensation for on-orbit space manipulators. Aerosp. Sci. Technol. 2014, 38, 30–40. [Google Scholar] [CrossRef]
  19. Palma, P.; Seweryn, K.; Rybus, T. Impedance control using selected compliant prismatic joint in a free-floating space manipulator. Aerospace 2022, 9, 406. [Google Scholar] [CrossRef]
  20. Liu, J.; Zhao, P.; Chen, K.; Zhang, X.; Zhang, X. 1U-Sized Deployable Space Manipulator for Future On-Orbit Servicing, Assembly, and Manufacturing. Space Sci. Technol. 2022, 2022, 1–14. [Google Scholar] [CrossRef]
  21. Hubert Delisle, M.; Christidi-Loumpasefski, O.-O.; Yalçın, B.C.; Li, X.; Olivares-Mendez, M.; Martinez, C. Hybrid-Compliant System for Soft Capture of Uncooperative Space Debris. Appl. Sci. 2023, 13, 7968. [Google Scholar] [CrossRef]
  22. Magarotto, M.; Manente, M.; Pavarin, D. Electric Propulsion for CubeSats: A Review. In Proceedings of the 72nd International Astronautical Congress (IAC), Dubai, United Arab Emirates, 25 October 2021. IAC-21,C4,8-B4.5A,1,x64088. [Google Scholar]
  23. Schäfer, F.; Montag, C.; Herdrich, G.; Laufer, R. Flight results of the PETRUS pulsed plasma thruster on the 3u CubeSat GreenCube. In Proceedings of the 73rd International Astronautical Congress (IAC), Paris, France, 18–22 September 2022. IAC-22,C4,8-B4.5A. [Google Scholar]
  24. Borelli, G.; Gaias, G.; Nakajima, Y.; Colombo, C.; Capuano, V.; Saggiomo, F.; Leccese, G.; Natalucci, S. Mission analysis and guidance and control for the SpEye inspection CubeSat. Acta Astronaut. 2024, 220, 75–87. [Google Scholar] [CrossRef]
  25. Ivanov, D.; Shestakov, S.; Ovchinnikov, M. Satellite formation flying control by mass exchange. Acta Astronaut. 2014, 102, 392–401. [Google Scholar]
  26. Clohessy, W.H.; Wiltshire, R.S. Terminal Guidance System for Satellite Rendezvous. J. Aerosp. Sci. 1960, 27, 653–658. [Google Scholar] [CrossRef]
  27. Sazonov, V.V.; Barbashova, T.F. Gravitational Orientation of Artificial Earth Satellites: A Special Course; CURS: Moscow, Russia, 2023; pp. 1–144. [Google Scholar]
  28. Badawy, A.; McInnes, C.R. On-orbit assembly using superquadric potential fields. J. Guid. Control Dyn. 2008, 31, 30–43. [Google Scholar] [CrossRef]
  29. Nesterov, Y.E. Convex Optimiation Methods (Metody Vypukloj Optimizacii); ICNMO Publishing House: Moscow, Russia, 2010; pp. 1–281. [Google Scholar]
  30. Fletcher, R. Practical Methods of Optimization; John Wiley & Sons: New York, NY, USA, 1987; pp. 1–455. [Google Scholar]
Figure 1. Considered objects: (1) cargo container, (2) installing platform, (3) truss structure, and (4) servicing spacecraft.
Figure 1. Considered objects: (1) cargo container, (2) installing platform, (3) truss structure, and (4) servicing spacecraft.
Aerospace 11 00635 g001
Figure 2. Angular motion of SS with rod: (1,2) two positions of servicing spacecraft with rod, (3) united body, and (4) united body convex hull.
Figure 2. Angular motion of SS with rod: (1,2) two positions of servicing spacecraft with rod, (3) united body, and (4) united body convex hull.
Aerospace 11 00635 g002
Figure 3. Scheme of system before and after servicing spacecraft’s repulsion.
Figure 3. Scheme of system before and after servicing spacecraft’s repulsion.
Aerospace 11 00635 g003
Figure 4. Example of distribution of final position deviation value Δ r .
Figure 4. Example of distribution of final position deviation value Δ r .
Aerospace 11 00635 g004
Figure 5. Values of motion parameters a i , b i , x 1 , and z 1 .
Figure 5. Values of motion parameters a i , b i , x 1 , and z 1 .
Aerospace 11 00635 g005
Figure 6. Boxplots of variation of the final position deviation (a) and of repulsion velocity error (b) ( N s a m p l e s = 100 ).
Figure 6. Boxplots of variation of the final position deviation (a) and of repulsion velocity error (b) ( N s a m p l e s = 100 ).
Aerospace 11 00635 g006
Figure 7. Solution existence depending on the united body’s angular velocity ω .
Figure 7. Solution existence depending on the united body’s angular velocity ω .
Aerospace 11 00635 g007
Figure 8. Heatmap of function g; the boundaries of the rod are represented by blue lines.
Figure 8. Heatmap of function g; the boundaries of the rod are represented by blue lines.
Aerospace 11 00635 g008
Figure 9. Example of dependence ε ( c i c i m a x ) and Φ ( c i c i m a x ) on ratio c i c i m a x for several δ values.
Figure 9. Example of dependence ε ( c i c i m a x ) and Φ ( c i c i m a x ) on ratio c i c i m a x for several δ values.
Aerospace 11 00635 g009
Figure 10. Large space truss structure to be assembled.
Figure 10. Large space truss structure to be assembled.
Aerospace 11 00635 g010
Figure 11. Repulsion velocity of x - and z-components’ denominator in (23) depending on transfer time.
Figure 11. Repulsion velocity of x - and z-components’ denominator in (23) depending on transfer time.
Aerospace 11 00635 g011
Figure 12. Repulsion velocity values at fixed motion time for randomly generated transfers between points of the structure with characteristic size L = 20 m.
Figure 12. Repulsion velocity values at fixed motion time for randomly generated transfers between points of the structure with characteristic size L = 20 m.
Aerospace 11 00635 g012
Figure 13. Repulsion velocity values at fixed motion time for randomly generated transfers between points of the structure with characteristic size L = 100 m.
Figure 13. Repulsion velocity values at fixed motion time for randomly generated transfers between points of the structure with characteristic size L = 100 m.
Aerospace 11 00635 g013
Figure 14. The SS’s deviation from the target point Δ r ( t ) (top) and the relative values c i / c i m a x of the constrained-motion parameters (19) (bottom) during the SS’s motion.
Figure 14. The SS’s deviation from the target point Δ r ( t ) (top) and the relative values c i / c i m a x of the constrained-motion parameters (19) (bottom) during the SS’s motion.
Aerospace 11 00635 g014
Figure 15. A set of examples of final position deviation during the numerical search by Algorithm 1. Each colored line corresponds to different example.
Figure 15. A set of examples of final position deviation during the numerical search by Algorithm 1. Each colored line corresponds to different example.
Aerospace 11 00635 g015
Figure 16. An example of final position deviation during the numerical search by Algorithm 2.
Figure 16. An example of final position deviation during the numerical search by Algorithm 2.
Aerospace 11 00635 g016
Figure 17. The SS’s deviation from the target point Δ r ( t ) during multiple transfers.
Figure 17. The SS’s deviation from the target point Δ r ( t ) during multiple transfers.
Aerospace 11 00635 g017
Figure 18. Relative values c i / c i m a x of the constrained-motion parameters (19) during multiple transfers.
Figure 18. Relative values c i / c i m a x of the constrained-motion parameters (19) during multiple transfers.
Aerospace 11 00635 g018
Figure 19. Intermediate stage of assembly: cargo container, servicing satellite, and its trajectories during the transfers.
Figure 19. Intermediate stage of assembly: cargo container, servicing satellite, and its trajectories during the transfers.
Aerospace 11 00635 g019
Figure 20. Demonstration of attitude change of the assembling structure.
Figure 20. Demonstration of attitude change of the assembling structure.
Aerospace 11 00635 g020
Figure 21. Repulsion velocity values during the space structure assembly.
Figure 21. Repulsion velocity values during the space structure assembly.
Aerospace 11 00635 g021
Table 1. Numerical simulation parameters.
Table 1. Numerical simulation parameters.
ParameterValueValue of VectorsStandard Deviation
m S S 20 kg υ U B O R F 10 2 m/s
m U B 650 kg r U B O R F 10 m
I y 5600 kg·m2 r 0 S R F , r 1 S R F 10 m
t f t 0 4000 s r U B S R F 5 m
R o r b 7000 km ω U B O R F 10 4 rad/s
Table 2. Main system parameters.
Table 2. Main system parameters.
VariableValue
The SS’s mass20 kg
The UB’s initial mass650 kg
Number of rods78
Length of rods5–10 m
Linear mass density of rods1 kg/m
Maximum SS repulsion velocity υ 0 B R F 0.05 m/s
m a x | ω U B O R F | 0.001 rad/s
m a x | r U B O R F | 500 m
m a x | υ U B O R F | 0.1 m/s
Radius of collision avoidance zone around rods30 cm
Manipulator capture radius50 cm
Maximum motion time t f m a x 10 4 s
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

Orlov, V.; Monakhova, U.; Ovchinnikov, M.; Ivanov, D. Fuelless On-Orbit Assembly of a Large Space Truss Structure Using Repulsion of the Service Spacecraft by Robotic Manipulators. Aerospace 2024, 11, 635. https://doi.org/10.3390/aerospace11080635

AMA Style

Orlov V, Monakhova U, Ovchinnikov M, Ivanov D. Fuelless On-Orbit Assembly of a Large Space Truss Structure Using Repulsion of the Service Spacecraft by Robotic Manipulators. Aerospace. 2024; 11(8):635. https://doi.org/10.3390/aerospace11080635

Chicago/Turabian Style

Orlov, Vladislav, Uliana Monakhova, Mikhail Ovchinnikov, and Danil Ivanov. 2024. "Fuelless On-Orbit Assembly of a Large Space Truss Structure Using Repulsion of the Service Spacecraft by Robotic Manipulators" Aerospace 11, no. 8: 635. https://doi.org/10.3390/aerospace11080635

APA Style

Orlov, V., Monakhova, U., Ovchinnikov, M., & Ivanov, D. (2024). Fuelless On-Orbit Assembly of a Large Space Truss Structure Using Repulsion of the Service Spacecraft by Robotic Manipulators. Aerospace, 11(8), 635. https://doi.org/10.3390/aerospace11080635

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