Skip to main content

ORIGINAL RESEARCH article

Front. Robot. AI, 24 December 2021
Sec. Space Robotics
This article is part of the Research Topic Robotic In-Situ Servicing, Assembly and Manufacturing View all 8 articles

Modeling a Controlled-Floating Space Robot for In-Space Services: A Beginner’s Tutorial

  • 1Surrey Space Centre, Department of Electronic and Electrical Engineering, University of Surrey, Guildford, United Kingdom
  • 2Lincoln Centre for Autonomous Systems (L-CAS), University of Lincoln, Lincoln, United Kingdom

Ground-based applications of robotics and autonomous systems (RASs) are fast advancing, and there is a growing appetite for developing cost-effective RAS solutions for in situ servicing, debris removal, manufacturing, and assembly missions. An orbital space robot, that is, a spacecraft mounted with one or more robotic manipulators, is an inevitable system for a range of future in-orbit services. However, various practical challenges make controlling a space robot extremely difficult compared with its terrestrial counterpart. The state of the art of modeling the kinematics and dynamics of a space robot, operating in the free-flying and free-floating modes, has been well studied by researchers. However, these two modes of operation have various shortcomings, which can be overcome by operating the space robot in the controlled-floating mode. This tutorial article aims to address the knowledge gap in modeling complex space robots operating in the controlled-floating mode and under perturbed conditions. The novel research contribution of this article is the refined dynamic model of a chaser space robot, derived with respect to the moving target while accounting for the internal perturbations due to constantly changing the center of mass, the inertial matrix, Coriolis, and centrifugal terms of the coupled system; it also accounts for the external environmental disturbances. The nonlinear model presented accurately represents the multibody coupled dynamics of a space robot, which is pivotal for precise pose control. Simulation results presented demonstrate the accuracy of the model for closed-loop control. In addition to the theoretical contributions in mathematical modeling, this article also offers a commercially viable solution for a wide range of in-orbit missions.

1 Introduction

Innovative space system technologies have revolutionized the lives on Earth, which explains why the global space race continues. The space economy worldwide is booming, and there is now a paradigm shift to the in-orbit services and manufacturing (IOSM) market. There is a significant commercial boost in the in-orbit services market, which is predicted to be $1 billion by 2030 (Satellite Catapult Applications, Fair-Space and Astroscale, 2021). The candidate in-orbit missions include the following: servicing and repairing high-value space assets including operational spacecraft, life extension, refueling, orbit correction, in-space assembly of space telescopes for Earth observation and astronomical observations, space-based power generation, and active debris removal, to name a few (Oda et al., 1996; Whittaker et al., 2000; Yoshida, 2001; Friend, 2008; Shengwei, 2013; Flores-Abad et al., 2014; Jaekel et al., 2015; Lee et al., 2016; Reed et al., 2016; Medina et al., 2017; Taylor et al., 2018; Li et al., 2019; Wilde et al., 2019; Jackson et al., 2020; Nair et al., 2020; Romano, 2021; Xu, 2021). Building satellite servicing and debris removal capabilities will open up bigger and longer-term markets linked to assembly and manufacturing in space (Satellite Catapult Applications, Fair-Space and Astroscale, 2021). As the in-orbit economy evolves, robotics and autonomous systems will play a pivotal role in many future IOSM missions. However, controlling a space robot in an extreme environment is significantly more complex than its terrestrial counterpart; this is the foremost hurdle to the success of IOSM missions.

In the context of servicing a cooperative and noncooperative target spacecraft using a servicer space robot, safety is paramount during the close-range approach, target capture, and postcapture operations. Approaching a target spacecraft using a servicer space robot can be achieved by two different modes of operation: free-flying and free-floating (Dubowsky and Papadopoulos, 1993; Moosavian and Papadopoulos, 2007). The free-flying approach uses reaction jets to facilitate a stabilized and controlled base for the robot manipulator in motion. The stable platform is favorable for the manipulator’s motion, however, on the expense of excessive fuel consumption and limited workspace. On the other hand, the free-floating approach utilizes an uncontrolled base to limit fuel consumption on the expense of dynamic singularities and an undefined workspace for the robot. There is a well-established literature on these two modes of operation, and they have benefits depending on the nature of the mission (Longman et al., 1987; Papadopoulos and Dubowsky, 1991, 1993; Xu and Kanade, 1992; Xu, 1993; Yoshida and Abiko, 2002; Menon et al., 2007; Artigas et al., 2015; Guang et al., 2018; Wilde et al., 2018). For instance, when communication between the space robot and the ground station is paramount, the free-flying mode is preferred as its controlled base can keep the antennas pointing toward the Earth. On the other hand, when reducing fuel consumption is prioritized, the free-floating mode is more suitable because of its uncontrolled base. In addition, recent pertinent works include numerical simulations, hardware-in-the-loop experiments, and guidance algorithms to capture and detumble a space object (Virgili-Llop et al., 2017). Despite these current advancements, the free-flying and free-floating operation modes exhibit unavoidable disadvantages for practical in-orbit missions.

An inherent problem with free-flying and free-floating operation modes is the undesired dynamic coupling effect due to the manipulator’s motion. This dynamic coupling effect leads to changes in the pose of the spacecraft base that are not always corrected, depending on the selected approach mode. In the free-flying mode, reaction jets will control the spacecraft base to maintain a fixed pose or attain one alternate pose while the arm is in motion. On the other hand, in the free-floating mode, the spacecraft base is uncontrolled, and it is free to change its pose in reaction to the motion of the arm. The free-flying mode has a limited workspace, whereas the free-floating mode has an undefined workspace because of its uncontrolled base. As a result, it is highly challenging to precisely navigate and control the space robot while avoiding obstacles and singularities during the approach phase.

In the close vicinity of the target, well-defined trajectories for the position and attitude of the spacecraft base are needed along with the trajectories for the arm’s joints to perform collision-free navigation. Instead of maintaining a fixed pose for the base spacecraft or letting the base spacecraft float in an uncontrolled manner, a controlled motion of the spacecraft base (i.e., achieving desired translation and rotation with time) is highly desirable. A coordinated movement of this nature is distinct and challenging compared with both free-flying and free-floating modes. It corresponds to the “controlled-floating” mode, previously introduced by Seddaoui and Saaj (2019), which offers redundancy to operate the space robot in an unlimited but well-defined workspace. The space robot is referred to as the controlled-floating space robot (CFSR) when operated in this mode. Contrarily to the free-flying and free-floating space robots, the CFSR uses its base’s controlled translation and rotation to help the arm reach the target’s grasping point.

Literature review shows that space robots suffer from singularities, depending on the mode of operation. It is also known that dynamic singularities affect free-floating space robots, whereas kinematic singularities affect free-flying space robots. As the CFSR can control its base’s motion, only kinematic singularities occur during the arm’s movement. It is possible to avoid kinematic singularities using the extra degrees of freedom offered by the spacecraft base. In short, the controlled-floating mode offers the benefits of both free-flying and free-floating modes of operation, and an accurate model of a CFSR will be highly beneficial to the end-users.

Although modeling the kinematics and dynamics of a free-flying and free-floating space robots is well addressed in the literature, the same for CFSR is less extensively studied in literature. Recently, Virgili-Llop et al. (2019) and Virgili-Llop and Romano (2019) have obtained simulation and experimental results for a CFSR approach via convex programming for maneuvering and capturing a tumbling object. In order to help the beginners in the discipline and foster an increased use of the CFSR approach, this article presents a step-by-step tutorial on deriving the equation of motion of a CFSR during its final approach phase. The full-scale nonlinear model presented can be used for controlling the pose of a CFSR during any of the aforementioned in-orbit missions. This article’s benefits are twofold: (1) helps beginners gain a sound theoretical foundation in mathematical modeling of orbital robots and (2) creates a better awareness of the benefits of CFSR for practical space missions.

Compared with the mathematical models of free-flying and free-floating space robots reported in the literature, the model of CFSR shown in this article includes a few new terms that represent different types of internal and external perturbations. More specifically, in this article, the motion of a CFSR is derived with respect to a reference frame attached to the moving target; thus, it complies with the close-proximity relative motion in-orbit. The center of mass (CoM) of the space robot changes because of both the arm and the base spacecraft’s motion. In addition, the overall inertia matrix at the CoM of the system is also not constant, but it changes during the motion of the CFSR. Corresponding changes in CoM, inertia matrix, Coriolis, and centrifugal terms are modeled in addition to accounting for the dynamic coupling effect. An accurate dynamic model together with a robust controller will enable the space robot, operating in the controlled-floating mode, to perform a safe and precise maneuver to approach the grasping point on the target. More details on the control architecture for a CFSR were published by Seddaoui and Saaj (2019).

In summary, the model of the CFSR presented in this article advances the state of the art of dynamic modeling a space robot. In addition to this article’s theoretical contributions, the highly accurate model presented is well suited for practical in-orbit missions. The rest of this article is organized as follows: Section 2 gives a recap of the CFSR and its mode of operation. The detailed tutorial that explains the mathematical modeling is covered in Section 3. The accuracy of the model presented is verified using simulations, and selected results are shown in Section 4. Finally, Section 5 summarizes the key inferences of the modeling methodology and added value of model presented for controlling the CFSR in orbit.

2 Background on CFSR

The concept of a CFSR originated with the idea of fusing the free-flying and free-floating operation modes of a space robot to utilize their individual pros. Both these modes of operation, individually, are undesirable for in-space services, which require precise operations; a limited or undefined workspace is unsatisfactory. Therefore, the identifiable challenges include collision-free navigation of the space robot with efficient fuel consumption.

A CFSR addresses these practical challenges by utilizing a closed-loop control for the motion of the space robot. Therefore, in addition to the arm’s joints, well-defined trajectories are needed to maintain the position and attitude of the spacecraft base. These desired trajectories are generated based on the linear and angular motion of the target spacecraft. In addition, the CFSR also takes into account the reaction forces and moments developed due to the actuation of the robotic arm. The dynamic coupling effect has to be controlled to withstand the undamped vibrations in the extremities of the space environment. The servicer spacecraft should match angular rates and keep zero relative attitude between the target to avoid misalignments.

Seddaoui and Saaj (2019) introduced the CFSR system first, but there is limited literature on its dynamic modeling. This tutorial aims at explaining the complex dynamics of CFSR, through a systematic mathematical formulation, at a high level of granularity. The details of path planning and robust control of CFSR are outside this article’s scope; details are available in Seddaoui (2020). A comparison between the existing modes of operation and the CFSR is required before any derivation. Figure 1 shows a comparison between the two main existing modes of operation, that is, free-flying and free-floating, as well as their subcategories introduced by Wilde et al. (2018) against the CFSR.

FIGURE 1
www.frontiersin.org

FIGURE 1. Comparison between existing modes of operation for space robot and the CFSR.

2.1 In-Orbit Relative Motion for the CFSR

The literature review reveals that modeling a space robot, which has a controlled spacecraft base, is always performed with respect to an inertial reference frame (Moosavian and Papadopoulos, 2007; Flores-Abad et al., 2014). This is to enable the utilization of Newton’s law of motion to derive the equation of motion of the space robot. Nevertheless, using a space robot for an in-orbit close-proximity approach is described as a three-body problem where it is preferable to reference the motion of the chasing body with respect to the moving and rotating frame of the target body, also known as local vertical local horizontal (LVLH) frame. Then it is necessary to transform the motion in LVLH back to the inertial frame to find the true actuating forces.

In this article, under the assumption that the space robot is operated at a very close proximity to the target (<10 m) and considering a short period (up to 150 s), three frames of reference are selected to derive the equation of motion of the CFSR as shown in Figure 2: the inertial frame I, the target LVLH frame T, and the base of the CFSR body frame B. The idea is to consider the motion of the CFSR in T as the controllable quantity in the closed-loop control and compute the fictitious forces, originating from the rotating frame T, as external forces affecting the system. This is to ensure that the true actuating forces and torques are identified in I.

FIGURE 2
www.frontiersin.org

FIGURE 2. Artistic illustration of the reference frames and vectors used in the mathematical model for the CFSR (Seddaoui, 2020).

The motion of the CFSR is described about its overall CoM, referred to here as ξ. Hence, the relative position is defined as the vector between the origin of T and the CoM ξ, referred to as rξrel. The basic equation of a relative position vector rξrel in I is (Curtis, 2013):

rξrel=rξr02,(1)

where vector rξR3 represents the inertial position vector of the CoM of the space robot ξ, r02R3 is the inertial position vector of the target, and rξrelR3 is the relative position vector between ξ and the origin of T, as shown in Figure 2. The corresponding relative velocity and acceleration vectors in I are defined as follows (Curtis, 2013):

vξv02=ṙξrel+Ω×rξrel(2a)
aξa02=r̈ξrel+Ω̇×rξrel+Ω×Ω×rξrel+2Ω×ṙξrel(2b)

where ΩR3 is the angular velocity of the rotating frame T with respect to I; vξR3 and v02R3 are the velocity vectors of the CoM ξ and the target, respectively; aξR3 and a02R3 are the acceleration vectors of the CoM ξ and the target, respectively. In the above equation and thereafter in this article, as it is often customary in the robotics literature, we use the word “vector” to indicate the column matrix of three scalar components of a physical vector in a specific Cartesian triad of axes. Furthermore, we use the symbol of vector cross product in order to indicate the three-by-three skew-symmetric matrix expressing in matricial form the cross product and built with the components of the algebraic vector preceding that symbol.

Eventually, the state vector that will be used in the closed-loop control of the CFSR includes the relative motion between the space robot and the target described by rξrel, ṙξrel and r̈ξrel. The fictitious velocity and acceleration involved in Eqs 2a, 2b are to be added as extra motion generated from the motion of T in order to find the true actuating forces and torques in I. The motion of the space robot’s base has both linear and angular components, as well as the generated interactions between the translational and rotational motion. This is computed through deriving the space robot’s motion about the overall CoM ξ.

3 Spacecraft Motion

The motion of the space robot’s base has both linear and angular components, as well as the generated interactions between the translational and rotational motion. This is computed through deriving the space robot’s motion about the overall CoM ξ.

3.1 The Linear Motion

During close-proximity maneuvers, the relative motion is of interest. Hence, the linear motion of the spacecraft base consists of its own relative motion, as well as that of the arm. It is translated by the relative linear momentum of the system as follows:

P=i=0nmiri̇P=Mtvξv02,(3)

where PR3 is the linear momentum of the system; mi and ri are, respectively, the mass and the position vector of the ith link of the space robot including the spacecraft base; and Mt is the total mass of the system. Differentiating Eq. 3and substituting Eq. 2b into Eq. 3 gives the equation describing the relative linear motion of the spacecraft base about the CoM (ξ) with respect to I:

FscI=Mtaξa02=Mtr̈ξrel+MtΩ̇×rξrel+MtΩ×Ω×rξrel+2MtΩ×ṙξrel,(4)

where FscIR3 is the force acting on the spacecraft base about ξ in I. The CoM of the full system ξ is different from the CoM of the base spacecraft. Also, the actuation forces and torques are applied on the spacecraft base rather than ξ. Therefore, in order to derive the equation of the relative linear motion of the spacecraft base, one needs to introduce the motion of B in T. Hence, using the equation for the relative velocity and acceleration, as seen in Eqs 2a, 2b, gives the equation for the relative velocity and acceleration ṙξrelR3 and r̈ξrelR3 when considering the motion of B. In other words, the following represents the expression of the first term of Eqs 2a, 2b when introducing the motion of B:

rξrel=rT+rξB(5a)
ṙξrel=ṙT+ṙξB+ωsc×rξB(5b)
r̈ξrel=r̈T+r̈ξB+ω̇sc×rξB+ωsc×ωsc×rξB+2ωsc×ṙξB(5c)

where rTR3 is the relative position vector of the spacecraft base defined as rT=r01r02, vector rξBR3 represents the relative position vector between ξ and B, and vector ωscR3 is the angular velocity of the spacecraft base with respect to T. Multiplying Eq. 5c by the total mass of the system gives the following expression for the linear motion of the spacecraft base in T, which is function of the spacecraft base’s position vector rT instead of the overall CoM’s position vector rξrel, as follows:

Fsc=Mtr̈T+Mtr̈ξBMtrξB×ω̇scMtωsc×rξB×ωsc+2Mtωsc×ṙξB,(6)

where FscR3 represents the force acting on the spacecraft base about ξ in T.

Substituting Eqs 5a, 5b into the last three remaining terms of Eq. 4 gives the following expression for the fictitious forces resulting from the rotation of frame T:

Fscfict=MtΩ̇×rT+rξB+MtΩ×Ω×ṙT+ṙξB+2MtΩ×ṙT+ṙξB+ω̇sc×rξB,(7)

where FscfictR3 represents the virtual forces generated from the rotation of the frame T.

Substituting Eqs 6, 7 into Eq. 4 gives the relative linear motion of the spacecraft base in I when considering the motion of frame T. This is expressed as follows:

FscI=Fsc+Fscfict.(8)

The first term of Eq. 8 represents the controlled linear motion of the space robot with respect to T. The second term is related to the rotation of T, which must be included in the equation of motion to account for the virtual forces that exist only in the rotating frame of reference. Expressing the linear motion of the space robot, represented by Eq. 6, in matrix form gives:

Fsc=MtEMt[rξB]×r̈Tω̇sc+0Mtωsc×rξB×ṙTωsc++Mtr̈ξB+2Mtωsc×ṙξB,(9)

where ER3×3 is the identity matrix. Eq. 9 describing the relative motion of the spacecraft base is one of the components constituting the equation of the relative motion of the CFSR.

3.2 The Angular Motion

The angular momentum of the multibody space robot about ξ, in terms of the inertia tensor and the angular velocity, is defined, with respect to a frame attached to the CoM of the spacecraft base and that does not rotate with the B, as follows (see Supplementary Material):

Lξ=Iξωξ+i=0nIiωi+i=0nIimi[riB]×[riB]×ωξB(10)

where IξR3×3 and IiR3×3 are, respectively, the inertia tensors of the full system expressed at the CoM ξ and the ith link of the arm; ωξ is the angular velocity about ξ with respect to a frame attached at the origin of B; and ωξB is the angular velocity about ξ with respect to B.

Let IiB=i=1nIimi[riB]×[riB]×, then Eq. 10 becomes:

Lξ=Iξωξ+i=1nIiωi+IiBωξB(11)

Given ωξ = ωsc + ωξB, Eq. 10 becomes:

Lξ=Iξωsc+IξωξB+i=1nIiωi+IiBωξB.(12)

The second and last terms of Eq. 12 are related to the rotation about ξ, generated by the motion of the arm. Differentiating Eq. 12 gives:

L̇ξ=ddtIξωsc+i=1nddtIiωi+ddtIξωξB+ddtIiBωξB.(13)

As the CFSR is one redundant system, finding the total inertia tensor Iξ involves two steps: use the parallel axis theorem to express the total inertia matrix of both the spacecraft base and the arm at the origin of B and then express the result at ξ. The total inertia matrix at the origin of B is as follows:

IξB=Isc+i=1nIimi[riB]×[riB]×.(14)

Considering the space robot with its overall CoM ξ, which is different from the CoM of the spacecraft base alone, the total inertia tensor of the system at ξ, using the parallel axis theorem, is:

Iξ=IξB+Mt[rξB]×[rξB]×.(15)

Substituting Eq. 14 into Eq. 15 gives:

Iξ=Isc+Mt[rξB]×[rξB]×+i=1nIimi[riB]×[riB]×=Isc+Mt[rξB]×[rξB]×+IiB.(16)

It is clear from Eq. 16 that the inertia tensor of the space robot is not constant as it depends on the position of the moving CoM ξ and the moving arm’s links. Hence, the derivatives of the first, second, and third terms of Eq. 13 are given as follows:

ddtIξωsc=Iξω̇sc+İξωsc,(17a)
i=1nddtIiωi=i=1nIiω̇i(17b)
ddtIξωξB=Iξω̇ξB+İξωξB(17c)
ddtIiBωξB=IiBωξḂ+İiBωξB(17d)

The derivatives of the inertia tensors Iξ and IiB (see Supplementary Material) still appear in the equation because they are function of position vectors that vary during the motion of the space robot, as seen in Eq. 16.

Substituting 17c and 17d into Eq. 13 gives the equation for the rotational motion of the space robot about ξ as follows:

τsc=L̇ξ=Iξω̇sc+İξωsc+i=0nIiω̇i+Iξω̇ξB+İξωξB++IiBω̇ξB+İiBωξB(18)

Eq. 18 describes the rotational dynamics of the space robot about the overall CoM ξ with respect to ξ. Expressing this vector of torques relative to the CoM of the spacecraft base (origin of B), with respect to B, requires the following transformation (Curtis, 2013):

τscI=τsc+rξB×FscI.(19)

Substituting Eqs 818 into Eq. 19 and simplifying give:

τscI=Mt[rξB]×r̈T+Isc+IiBω̇scMtrξB×ωsc×rξB×ωsc++İξωsc+i=0nIiω̇i+İξωξB+Iξω̇ξB+İiBωξB+IiBω̇ξB++MtrξB×r̈ξB+2MtrξB×ωsc×ṙξB+rξB×Fscfict.(20)

From Eq. 20, the angular motion of the space robot in T can be represented in matrix form as follows:

τsc=Mt[rξB]×Isc+IiBr̈Tω̇sc+0İξωscMt[rξB]×ωsc×rξB×ṙTωsc+i=0nIiω̇i+İξωξB+Iξω̇ξB+İiBωξB+IiBω̇ξB+MtrξB×r̈ξB+2MtrξB×ωsc×ṙξB.(21)

3.3 The Overall Motion of the Spacecraft Base

From Eqs 921, the overall equation describing the linear and angular motion of the spacecraft base in T is:

fsc=Fscτsc=MtEMt[rξB]×Mt[rξB]×Isc+IiBr̈Tω̇sc++0Mtωsc×rξB×0İξωscMt[rξB]×ωsc×rξB×ṙTωsc++Mtr̈ξBMtrξB×r̈ξB+i=0nIiω̇i+Iξω̇ξB+IiBω̇ξB++2Mtωsc×ṙξBİξωξB+İiBωξB+2MtrξB×ωsc×ṙξB,(22)

which can be written in a compact form as follows:

fsc=Fscτsc=DvDvωDωvDωr̈Tω̇sc+0Cv0CωṙTωsc+DξB+CξB=DscẌ+CscẊ+DξB+CξB,(23)

where DvR3×3, DvωR3×3, DωvR3×3 and DωR3×3 are mass submatrices related the linear and angular motion of the spacecraft base, as well as the interaction between the linear and angular motion. Vectors CvR3×3 and CvωR3×3 involve the Coriolis and centrifugal terms originating from the motion of the spacecraft base, and vectors DξB and CξB are related to the motion of the CoM ξ. The compacted matrix DscR6×6 is related to the linear and angular motion of the spacecraft base, and the compacted vector CscR6×6 involves the Coriolis and centrifugal forces.

4 Manipulator Dynamics

The mathematical model for the dynamics of a terrestrial robotic arm usually involves both the kinetic and the potential energy (Spong and Vidyasagar, 2008). However, space robots operate in a microgravity environment. For this reason, the terms related to the potential energy in the equation of motion are ignored. Hence, only the kinetic energy was used in the Lagrange–Euler method to derive the equation of motion of the robotic arm of the CFSR.

The Lagrange–Euler equation for an n DoF robotic arm is (Spong and Vidyasagar, 2008):

τi=ddtLθ̇iLθiwithi=1,,n,(24)

where τi is the torque applied the ith joint; L is the Lagrangian function of system L = KE, where KE is the kinetic energy of the arm; θi represents the ith joint angle; and θ̇i is the velocity of the ith joint.

4.1 Kinetic Energy of an n DoF Manipulator

All moving objects have a kinetic energy, and this energy varies according to the mass of the object and the rate of change of the motion (velocity). For a multilink robotic arm, the velocity of each link is represented by the Jacobian matrix and the joint velocity θ̇. The expression for the linear and angular velocities for the ith joint is expressed as follows (Spong and Vidyasagar, 2008):

vi=Jvmiθ̇,ωi=Jωmiθ̇,,(25)

where JvmiR3×n and JωmiR3×n are, respectively, the linear and angular terms of the 6 × n Jacobian matrix of the ith CoM of the arm, and RLiR3×3 is the rotation matrix of the ith link with respect to T.

The general equation for the overall kinetic energy for an n DoF manipulator is:

KE=12i=1nmivivi+ωiIiωi,(26)

where Ii is the inertia tensor of the ith link about a frame attached at the CoM of the ith link. In order to express the inertia tensor in T, one has to perform a transformation using the rotation matrix RLi: IiT=RLiIiRLi. Hence, by substituting Eq. 25 into Eq. 26, the expression for the arm’s total kinetic energy is:

KE=12θ̇i=1nmiJviJvi+JωiRLiIiRLiJωiθ̇.(27)

Eq. 27 can be written in matrix form:

KE=12θ̇Dmθ̇,(28)

where DmRn×n is the symmetric positive definite mass matrix of the robotic manipulator, and it is expressed as:

Dm=i=1nmiJviJvi+JωiRLiIiRLiJωi.(29)

4.2 The Coriolis and Centrifugal Forces

The joints’ rotational motion results in extra forces known as the Coriolis and centrifugal forces. These forces are computed as follows (Spong and Vidyasagar, 2008):

ckj=i=1ncijkθi̇=i=1n12dkjθj+dkiθj+dijθkθi̇,(30)

where ckj are terms constituting the matrix of Coriolis and centrifugal forces denoted as CmRn×n.

From Eqs 29, 30, the equation describing the linear and angular motion of the arm, the target frame T is:

τm=Dmθ̈+Cmθ̇.(31)

5 Jacobian Matrix for the Kinematics of the CFSR

The Jacobian matrix relates the velocity of the end-effector, in the Cartesian space, with the velocity of the space robot, in the configuration space. It describes the kinematics of multibody chain, such as the space robot, at the velocity level. It is expressed as:

ṙeωe=JscJmẊθ̇.(32)

The submatrices constituting the Jacobian matrix are defined as follows:

Jsc=JvscJωsc,Jm=JvmJωm,

where JvscR3×6 and JωscR3×6 are, respectively, the linear and rotational components of the Jacobian matrix for the spacecraft base; and JvmR3×n and JωmR3×n are, respectively, the linear and rotational components of the Jacobian matrix for the arm.

The kinematics of a spacecraft can be expressed using either Euler angles or quaternions. Depending on which of these two methods is selected to describe the motion of the space robot, the derivation of the Jacobian matrix varies accordingly. Moreover, when designing a path for the end-effector, the Jacobian matrix related to the end-effector’s position is required, whereas the Jacobian involved in the dynamics is related to the position of the ith CoM of the chain. In the following section, the Jacobian matrix of the CFSR is derived for a potential use with both Euler and quaternions. Also, the end-effector’s velocity, as well as the ith CoM velocity, is described in T.

5.1 Jacobian Matrix of the End-Effector in T Using Euler Rate

Deriving the Jacobian matrix of the space robot for use with Euler rate involves several partial derivatives of the rotation matrices as well the spacecraft base’s transformation matrix Rω defined as follows (Wie, 2008):

Rω=10sin(β)0cos(α)sin(α)cos(β)0sin(α)cos(α)cos(β)(33)

5.1.1 The Linear Jacobian of the End-Effector

The linear velocity of the end-effector, in the Cartesian space with respect to T, is derived using Eq. 2a as follows:

ṙe=ṙT+ṙeB+ωsc×reB,(34)

where reBR3 is the vector from the origin of B to the end-effector, as seen in Figure 3 and described as:

reB=Rscs0+i=1nRLi1li1,(35a)
ṙeB=Ṙscs0+i=1nṘLi1li1,(35b)

where s0R3 is the vector from the origin of B to the first joint of the arm. and liR3 is the vector from the ith joint to the (i + 1)th joint.

FIGURE 3
www.frontiersin.org

FIGURE 3. Vector representation for the space robot’s physical parameters (Seddaoui, 2020).

The steps of the differentiation of a rotation matrix were presented by Umetani and Yoshida (1989) as follows:

Ṙsc=Rscαα̇+Rscββ̇+Rscγγ̇=Rscϕϕ̇.(36a)
ṘLi1=k=0i1RLi1qkq̇k=RLi1ϕϕ̇+k=1i1RLi1θkθ̇k.(36b)

where alpha is the roll angle, beta is pitch angle, and gamma is the yaw angle of the spacecraft base. In addition, phi is the vector of spacecraft base orientation in the target frame.

Using Eq. 36, the linear velocity of the end-effector described by Eq. 34 becomes:

ṙe=ṙT+Rscϕs0+i=1nRLi1ϕli1ϕ̇+i=1nk=1i1RLi1θkli1θ̇[reB]×ωsc.(37)

The last term of Eq. 37 involves the body rate ωsc which can be a function of Euler rate and Rω using Eq. 33. The linear velocity described by Eq. 37 then becomes:

ṙe=ṙT+Rscϕs0+i=1nRLi1ϕli1ϕ̇+i=1nk=1i1RLi1θkli1θ̇[reB]×Rωϕ̇.(38)

From Eq. 38, the linear Jacobian submatrices for both the spacecraft base and the arm are:

Jvsc=EJvsc(1)Jvsc(2)Jvsc(3)Jvsc(1)=Rscαs0+i=1nRLi1αli1[reB]×Rωe1,Jvsc(2)=Rscβs0+i=1nRLi1βli1[reB]×Rωe2,Jvsc(3)=Rscγs0+i=1nRLi1γli1[reB]×Rωe3,Jvm=Jvm(1)Jvm(2)Jvm(n)Jvm(j)=i=1jk=1i1RLi1θkli1,j=1,2,,n,(39)

where ER3×3, Jvsc(1),Jvsc(2),Jvsc(3)R3, Jvm(j)R3 and e1, e2, e3 are unit vectors.

5.1.2 The Rotational Jacobian of the End-Effector

The angular velocity of the end-effector can be expressed as:

ωe=ωsc+ωeBωeB=i=1nRLi00θi̇,(40)

where ωeR3 is the angular velocity of the end-effector in T, ωeBR3 is the angular velocity of the end-effector in B, and ωscR3 is the angular velocity of the spacecraft base in T.

Using Eqs 3340, ωe becomes:

ωe=Rωα̇β̇γ̇+i=1nRLi00θi̇.(41)

From Eq. 41, the rotational Jacobian submatrices for both the spacecraft base and the arm are:

Jωsc=0Jωsc(1)Jωsc(2)Jωsc(3)Jωsc(1)=Rωe1Jωsc(2)=Rωe2Jωsc(3)=Rωe3Jωm=Jωm(1)Jωm(2)Jωm(n)Jωm(j)=i=0jRLie3j=1,2,,n,(42)

where 0R3×3Jωsc(j)R3 and Jωm(j)R3 and e1, e2 and e3 are unit vectors.

5.1.3 The Full Jacobian Translating the Velocities of the End-Effector

From Eqs 3942, the overall Jacobian matrix for the space robot is as follows:

JscJm=EJvsc|Jvm0Jωsc|Jωm.(43)

5.2 Jacobian Matrix of the ith CoM in T Using Euler Rate

The Jacobian described by Eq. 43 is used to find the position of the end-effector in the Cartesian space, given a configuration or the inverse procedure for path planning. However, the matrices Jvmi and Jωmi of Eq. 29 relate to the velocity of the CoM of the ith link with the velocities in the configuration space. For this reason, the equations for Jvmi and Jωmi in Eq. 29 are different from those for Jvm and Jωm in Eqs 3942.

5.2.1 The Linear Jacobian of the ith CoM

Using Eq. 2a, the linear velocity of the ith CoM, with respect to T, is determined as follows:

ṙi=ṙT+ṙiB+ωsc×riB,(44)

where riB is the vector from the origin of B to the ith CoM, as seen in Figure 3. It is represented by the following equation:

riB=j=1iRLjbj+RLj1sj1,(45)

where biR3 is the vector from the ith joint to the CoM of the ith link, and si1R3 is the vector from the CoM of the (i − 1)th link to the ith joint as shown in Figure 3.

Differentiating Eq. 45 and using the derivatives of the rotation matrix in Eq. 36 give:

ṙiB=j=1iRLjϕbj+RLj1ϕsj1ϕ̇+j=1ik=1iRLjθkbj+k=1j1RLj1θksj1θ̇.(46)

By substituting Eq. 46 into Eq. 44, the equation for the linear velocity of the ith CoM becomes:

ṙi=ṙT+j=1iRLjϕbj+RLj1ϕsj1ϕ̇+j=1ik=1jRLjθkbj+k=1j1RLj1θksj1θ̇[riB]×Rωϕ̇.(47)

From Eq. 47, the linear Jacobian matrices, derived for the dynamics, for both the base spacecraft and the arm are:

Jvsci=EJvsci(1)Jvsc(2)Jvsc(3)Jvsci(1)=j=1iRLjαbj+RLj1αsj1[riB]×Rωe1Jvsci(2)=j=1iRLjβbj+RLj1βsj1[riB]×Rωe2Jvsci(3)=j=1iRLjγbj+RLj1γsj1[riB]×Rωe3Jvmi(j)=j=0ik=1jRLjθkbj+k=1j1RLj1θksj1,

where Jvsci is the ith linear Jacobian submatrix of the spacecraft base resulting from the motion of the base itself and the arm, and Jvmi(j) is the jth column of the linear Jacobian matrix of the ith CoM of the arm.

5.2.2 The Rotational Jacobian of the ith CoM

The angular velocity of the ith CoM in T is expressed as follows:

ωe=ωsc+ωiscωisc=j=1iRLj00θi̇,(48)

where ωiscR3 is the angular velocity of the ith CoM in B. Based on Eq. 42, the rotational Jacobian for the ith CoM is:

Jωsc=0Jωsc(1)Jωsc(2)Jωsc(3)Jωsc(1)=Rωe1Jωsc(2)=Rωe2Jωsc(3)=Rωe3Jωmi(j)=j=1iRLje3,(49)

where Jωmi(j) is the jth column of the rotational Jacobian matrix of the ith CoM of the manipulator.

5.3 Jacobian Matrix of the End-Effector in T Using Body Rate

It is known that singularities can occur when using the Euler rate ϕ̇ to express the attitude of the spacecraft base. Hence, quaternions are preferred to avoid singularities. For this reason, the Jacobian matrix has to be expressed with respect to the body rate ωsc to allow an easy transformation to quaternions. In this case, there are no partial derivatives of the rotation matrices involved.

5.3.1 The Linear Jacobian of the End-Effector Using Body Rate

The linear velocity of the end-effector, described by Eq. 34, is here derived in a different manner to maintain the body rate term in the equation. This is performed using another method for the derivative of the rotation matrices Rsc and RLi, hereafter presented as (Spong and Vidyasagar, 2008):

Ṙsc=[ωsc]×Rsc=[e1]×Rscωscx+[e2]×Rscωscy+[e3]×Rscωscz,(50)

and

ṘLi1=j=1i1RLj1[ωj]×RLi1=[ωsc]×RLi1li1+j=1i1RLj1[ωj]×RLi1=[e1]×RLi1li1ωscx+[e2]×RLi1li1ωscy+[e3]×RLi1li1ωscz++j=1i1RLj1[e3]×RLi1θ̇j.(51)

Substituting Eqs 50, 51 into Eq. 34 gives:

ṙe=ṙT+[ωsc]×Rscs0+i=1n[ωsc]×RLi1li1[reB]×ωsc,=ṙT+[e1]×Rscs0+i=1n[e1]×RLi1li1[reB]×e1ωscx+[e2]×Rscs0+i=1n[e2]×RLi1li1[reB]×e2ωscy+[e3]×Rscs0+i=1n[e3]×RLi1li1[reB]×e3ωscz+i=1nj=1i1RLj1[e3]×RLi1li1θj̇.(52)

Expressing Eq. 52 in matrix form gives the following Jacobian submatrices for both the spacecraft base and the arm:

Jvsc=EJvsc(1)Jvsc(2)Jvsc(3)Jvsc(1)=[e1]×Rscs0+i=1n[e1]×RLi1li1[reB]×e1,Jvsc(2)=[e2]×Rscs0+i=1n[e2]×RLi1li1[reB]×e2,Jvsc(3)=[e3]×Rscs0+i=1n[e3]×RLi1li1[reB]×e3,Jvm=Jvm(1)Jvm(2)Jvm(n)Jvm(j)=k=1ji=1i1RLk1[e3]×RLi1li1,j=1,2,,n.(53)

5.3.2 The Rotational Jacobian of the End-Effector Using Body Rate

Using Eq. 48, the rotational Jacobian submatrix, used to find the angular velocity of the end-effector, with respect to body rate expressed in matrix form is:

Jωsc=0e1e2e3Jωm=Jωm(1)Jωm(2)Jωm(n)Jωm(j)=i=1jRLie3j=1,2,,n,(54)

where 0R3×3 and e1, e2, e3 are unit vectors.

5.4 Jacobian Matrix of the ith CoM in T Using Body Rate

Similar to the Jacobian matrices of the links’ CoMs using Euler rate, deriving these matrices using body rate is necessary when quaternions are utilized.

5.4.1 The Linear Jacobian of the ith CoM Using Body Rate

Differentiating Eq. 45 using the derivatives of the rotation matrices in Eqs 50, 51 gives:

ṙiB=j=1i[ωsc]×RLjbj+[ωsc]×RLj1sj1+j=1ik=1jRLk1[ωk]×RLjbj+k=1j1RLk1[ωk]×RLj1sj1.(55)

Substituting Eq. 55 into Eq. 44 gives the expression for the linear velocity of the ith link using the spacecraft body rate as follows:

ṙi=ṙT+j=1i[e1]×RLjbj+[e1]×RLj1sj1[riB]×e1ωx+j=1i[e2]×RLjbj+[e2]×RLj1sj1[riB]×e2ωy+j=1i[e3]×RLjbj+[e3]×RLj1sj1[riB]×e3ωz+j=1ik=1jRLk1[e3]×RLjbj+k=1j1RLk1[e3]×RLj1sj1θ̇.(56)

Expressing Eq. 56 in matrix form gives the following linear Jacobian submatrices to compute the linear velocity of the ith link:

Jvsci=EJvsci(1)Jvsci(2)Jvsci(3)Jvsci(1)=j=1i[e1]×RLjbj+[e1]×RLj1sj1[riB]×e1Jvsci(2)=j=1i[e2]×RLjbj+[e2]×RLj1sj1[riB]×e2Jvsci(3)=j=1i[e3]×RLjbj+[e3]×RLj1sj1[riB]×e3Jvmi(j)=j=1ik=1jRLk1[e3]×RLjbj+k=1j1RLk1[e3]×RLj1sj1j=1,2,,n.(57)

5.4.2 The Angular Jacobian of the ith CoM Using Body Rate

It is similar to the rotational Jacobian in Eq. 49, and it is expressed as follows:

Jωsc=0e1e2e3Jωmi(j)=j=1iRLje3j=1,2,,n(58)

5.4.3 Simulation Validation of the Kinematic Equations

It is important to validate the kinematic model of the space robot before integrating it with the equations of the dynamics. As presented in this article, the equations for the motion of the end-effector and the ith CoM are different. The inverse kinematics computation uses Eq. 35a as it represents the end-effector’s motion. The motion of the links’ CoM, represented by Eq. 44, is used to compute the robot’s dynamics. The simulation-based validation uses these two equations to define the motion of the manipulator. Testing these two equations for different space robot configurations will help validate the kinematics equations by verifying if, for each configuration, the position of the ith CoM corresponds to the ith link.

For the simulations, the joint angles were computed separately from the positions of the end-effector and the links’ CoM and were plotted in the same plot to check if they coincide. The check was performed not only visually but also by comparing the positions of each CoM with respect to its link, in six arbitrary sets of space robot configurations, including the arm’s joints and spacecraft’s attitude. This process was carried out to validate the accuracy of the two sets of equations before integrating them with the kinematics and dynamics.

The results are shown in Figure 4, where the first three configurations have a fixed spacecraft base and joints, respectively, equal to: θ = [11° 11° 11° 0° 0° 0°]′, θ = [0° 0° 0° 0° 0° 0°]′, and θ = [0° 0° 30° 0° 0° 0°]′. The last three configurations have a fixed joint angle equal to θ = [0°–30° 0° 30° 30° 0°]′ and three different sets of base rotations defined, respectively, as rotation about z axis, rotations about z and y axes, and rotations about z, y and x axes. This simulation-based validation proves the correctness of the equations used in the derivation of the Jacobian matrices that describe the kinematics of the CFSR.

FIGURE 4
www.frontiersin.org

FIGURE 4. Links CoMs’ positions using different joint angles for arm and base attitude (Seddaoui, 2020).

6 The Dynamic Coupling Between the Arm and the Spacecraft Base

The dynamic coupling describes the interactions between the arm and its base during the motion of the space robot. The impact of the motion of the base on the arm is introduced through the Jacobian matrix. In other words, desired changes in the pose of the spacecraft base will systematically change the end-effector’s position through the Jacobian matrix. This feature is used to aid the arm to reach the target. However, the motion of the arm produces undesired changes in the pose of the spacecraft base, which in return affect the end-effector’s position. This results from the motion of the system’s CoM, which corresponds to the last two terms of Eq. 22. These two terms are:

DξB+CξB=Mtr̈ξBMtrξB×r̈ξB+i=1nIiω̇i+Iξω̇ξB+IiBω̇ξB++2Mtωsc×ṙξBİξωξB+İiBωξB+2MtrξB×ωsc×ṙξB.(59)

The angular velocity ωξB is expressed as follows (Curtis, 2013):

ωξB=rξB×ṙξBrξB2.(60)

The derivative of ωξB, described by Eg (60), is defined as follows:

ω̇ξB=rξB×r̈ξBrξB2+2rξB×ṙξBrξB3.(61)

Vector rξB and its first derivative are defined as follows:

rξB=1Mti=1nmiriB,(62a)
ṙξB=1Mti=1nmiṙiB.(62b)

In order to derive the first term of Eq. 59, the second derivative of vector riB has to be computed. Its first derivative is defined by Eq. 55, which was used to find the linear term of the Jacobian matrix of the ith link in T. Contrarily, vector riB is here defined in B. This changes the origin of the rotation matrices from T to B. Hence, the first derivative of vector riB is computed as follows:

ṙiB=j=1ik=1jRLk1[ωk]×RLjbj+k=1j1RLk1[ωk]×RLj1sj1.(63)

An inspection of Eq. 63 shows that it is similar to the linear part of the ith link’s Jacobian matrix as described in Eq. 57. Therefore, the first derivative of vector riB is given as:

ṙiB=i=1nJvmiθ̇.(64)

The second derivative of vector riB is then defined as:

r̈iB=i=1nJvmiθ̈+J̇vmiθ̇.(65)

The expression for the first term of Eq. 59 is defined as follows:

r̈ξB=1Mti=1nmiJvmiθ̈+1Mti=1nmiJ̇vmiθ̇.(66)

The termi=1nIiω̇i is to be further derived to make it function of the joint velocity θ̇. This is achieved through finding the equation for ωi as follows:

ωi=j=1iRLjθj.(67)

Eq. 67 is similar to the rotational part of the ith link’s Jacobian matrix. Hence:

i=1nIiωi=i=1nIiJωmiθ̇(68a)
i=1nIiω̇i=i=1nIiJωmiθ̈+i=1nIiJωmi̇θ̇.(68b)

Substituting Eqs 62a, 62b, (Eq. 66), and (Eq. 68b) into (Eq. 59) and rearranging give:

DξB+CξB=i=1nmiJvmi1+Iξ+IiBMtrξB2[rξB]×i=1nmiJvmi+i=1nIiJωmiθ̈+Csc.mvCsc.mωθ̇,(69)

where matrices Csc.mv and Csc.mω are defined as follows:

Csc.mv=i=1nmiJ̇vmi+2ωsc×Jvmi,Csc.mω=i=1nIiJ̇ωmi+2rξB×ωsc×Jvmi+İξ+İiBMtrξB2[rξB]×i=1nmiJvmi+2Iξ+IiBMtrξB3[rξB]×i=1nmiJvmi.

Writing Eq. 69 in a compact matrix form gives:

DξB+CξB=Dsc.mvDsc.mωθ̈+Csc.mvCsc.mωθ̇=Dsc.mθ̈+Csc.mθ̇,(70)

where Dsc.mvR3×n and Dsc.mωR3×n are matrices related to the dynamic coupling between the arm and its spacecraft base, and matrices Csc.mvR3×n and Csc.mωR3×n involve the Coriolis and centrifugal terms originating from the interaction between the arm and the spacecraft base.

The model presented in this article has been extensively used in simulations to test the controllers presented in Seddaoui et al. (2019) and Seddaoui and Saaj (2019). The controllers are model-based controllers designed as a combined H controller (Seddaoui and Saaj, 2019) and its extension as an adaptive controller to minimize the control forces and torques (Seddaoui et al., 2019). Moreover, the CFSR model and the adaptive H controller were used to test the trajectory planning algorithm introduced by Seddaoui and Saaj (2021). Simulation results showcasing the desired and actual spacecraft position and orientation trajectories, the desired and actual trajectories of the arm joints, and the final Cartesian motion of the CFSR when using the model presented in this article are presented by Seddaoui and Saaj (2021).

7 Conclusion

An accurate mathematical model that captures the complex dynamics of a space robot is helpful for effective pose control of space robot in extreme conditions. The controlled-floating mode of operation of a space robot is distinct from the free-flying and free-floating modes in terms of the nature of motion and corresponding system model. This tutorial guides the reader through a systematic process to derive the nonlinear dynamic model of CFSR. The equation of motion derived with respect to a frame attached to the target spacecraft complies with the close-proximity relative motion in orbit. The changes in the CoM and subsequent changes to the overall inertia matrix of the coupled system were considered. The refined model of the CFSR presented includes a few additional mathematical terms that account for these perturbations. In short, this tutorial helps users to compute the CFSR model as accurately and efficiently as possible.

It is important to note that the model presented here can be used for space robots of any size and mass. This is possible because the physical parameters of the space robot are the inputs to the generated matrices, which are the building blocks of the final model. One exception is the array of rotation matrices specific to each space robot as it needs to be defined beforehand. Unlike free-flying and free-floating space robots, a CFSR will require an advanced path-planning algorithm to generate optimal trajectories for both the arm and its base. This path planner must consume as little energy as possible and produce a minimal dynamic coupling effect. Moreover, a controller capable of executing the desired motion with a small control effort is also vital for the mission’s success.

The equation of motion of the controlled-floating space robot cannot be reduced to that of the freefloating space robot because they are fundamentally different. The latter is based on the generalized Jacobian matrix that is specific to free-floating robots with an assumed fixed CoM. In this method, the dynamics and kinematics of the space robot are tightly related and part of the same equation that is the generalized Jacobian matrix. Alternatively, the free-flying mode could be a subset of the CFSR when only controlled translation of the base spacecraft is required without the need to change its attitude. This state is different from free-floating mode where both translation and rotation of the base spacecraft are uncontrolled. This scenario is possible theoretically, but more research is needed to verify this statement, which is outside the scope of this article.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author Contributions

AS - Lead researcher CS - PhD Supervisor MN - Simulation and editorial support.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors wish to acknowledge the United Kingdom. Space Agency, the Surrey Satellite Technology Ltd., and the Algerian Space Agency for funding this research.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/frobt.2021.725333/full#supplementary-material

References

Artigas, J., De Stefano, M., Rackl, W., Lampariello, R., Brunner, B., Bertleff, W., et al. (2015). “The Oos-Sim: An On-Ground Simulation Facility for On-Orbit Servicing Robotic Operations,” in IEEE International Conference on Robotics and Automation (ICRA). Seattle, 26–30. doi:10.1109/icra.2015.7139588

CrossRef Full Text | Google Scholar

[Dataset] Satellite Catapult Applications, Fair-Space and Astroscale (2021). A $1bn uk Opportunity for In-Orbit Services. Available at: https://sa.catapult.org.uk/news/in-orbit-servicing-capability (Last Accessed June 13, 2021).

Google Scholar

Curtis, H. D. (2013). Orbital Mechanics for Engineering Students. Florida: Butterworth-Heinemann.

Google Scholar

Dubowsky, S., and Papadopoulos, E. (1993). The Kinematics, Dynamics, and Control of Free-Flying and Free-Floating Space Robotic Systems. IEEE Trans. Robot. Automat. 9, 531–543. doi:10.1109/70.258046

CrossRef Full Text | Google Scholar

Flores-Abad, A., Ma, O., Pham, K., and Ulrich, S. (2014). A Review of Space Robotics Technologies for On-Orbit Servicing. Prog. Aerospace Sci. 68, 1–26. doi:10.1016/j.paerosci.2014.03.002

CrossRef Full Text | Google Scholar

Friend, R. B. (2008). “Orbital Express Program Summary and mission Overview,” in SPIE Defense and Security Symposium (Orlando, FL: International Society for Optics and Photonics), 695–803. doi:10.1117/12.783792

CrossRef Full Text | Google Scholar

Guang, Z., Heming, Z., and Liang, B. (2018). Attitude Dynamics of Spacecraft with Time-Varying Inertia during On-Orbit Refueling. J. Guidance, Control Dyn., 1–11. doi:10.2514/1.g003474

CrossRef Full Text | Google Scholar

Jackson, L., Saaj, C. M., Seddaoui, A., Whiting, C., Eckersley, S., and Hadfield, S. (2020). Downsizing an Orbital Space Robot: A Dynamic System Based Evaluation. Adv. Space Res. 65, 2247–2262. doi:10.1016/j.asr.2020.03.004

CrossRef Full Text | Google Scholar

Jaekel, S., Lampariello, R., Panin, G., Sagardia, M., Brunner, B., Porges, O., et al. (2015). “Robotic Capture and De-orbit of a Heavy, Uncooperative and Tumbling Target in Low Earth Orbit,” in 13th Symposium on Advanced Space Technologies in Robotics and Automation (Noordwijk, Netherlands: ASTRA), 11–13.

Google Scholar

Lee, N., Backes, P., Burdick, J., Pellegrino, S., Fuller, C., Hogstrom, K., et al. (2016). Architecture for In-Space Robotic Assembly of a Modular Space Telescope. J. Astronomical Telescopes, Instr. Syst. 2, 041207. doi:10.1117/1.jatis.2.4.041207

CrossRef Full Text | Google Scholar

Li, W.-J., Cheng, D.-Y., Liu, X.-G., Wang, Y.-B., Shi, W.-H., Tang, Z.-X., et al. (2019). On-orbit Service (Oos) of Spacecraft: A Review of Engineering Developments. Prog. Aerospace Sci. 108, 32–120. doi:10.1016/j.paerosci.2019.01.004

CrossRef Full Text | Google Scholar

Longman, R. W., Lindbergt, R. E., and Zedd, M. F. (1987). Satellite-Mounted Robot Manipulators - New Kinematics and Reaction Moment Compensation. Int. J. Robotics Res. 6, 87–103. doi:10.1177/027836498700600306

CrossRef Full Text | Google Scholar

Medina, A., Tomassini, A., Suatoni, M., Avilés, M., Solway, N., Coxhill, I., et al. (2017). Towards a Standardized Grasping and Refuelling On-Orbit Servicing for Geo Spacecraft. Acta Astronautica 134, 1–10. doi:10.1016/j.actaastro.2017.01.022

CrossRef Full Text | Google Scholar

Menon, C., Busolo, S., Cocuzza, S., Aboudan, A., Bulgarelli, A., Bettanini, C., et al. (2007). Issues and Solutions for Testing Free-Flying Robots. Acta Astronautica 60, 957–965. doi:10.1016/j.actaastro.2006.11.014

CrossRef Full Text | Google Scholar

Moosavian, S. A. A., and Papadopoulos, E. (2007). Free-flying Robots in Space: an Overview of Dynamics Modeling, Planning and Control. Robotica 25. doi:10.1017/s0263574707003438

CrossRef Full Text | Google Scholar

Nair, M. H., Saaj, C. M., Esfahani, A. G., Nanjangud, A., Eckersley, S., and Bianco, P. (2020). “In-space Robotic Assembly and Servicing of High-Value Infrastructure,” in Proc. International Astronautical Congress. Cyberspace Edition (Dubai: IAC2020).

Google Scholar

Oda, M., Kibe, K., and Yamagata, F. (1996). “Ets-vii, Space Robot In-Orbit experiment Satellite,” in Proc. 1996 IEEE International Conference on Robotics and Automation, Minnesota, April 22-28, 739–744.

Google Scholar

Papadopoulos, E., and Dubowsky, S. (1991). “Coordinated Manipulator/spacecraft Motion Control for Space Robotic Systems,” in IEEE International Conference on Robotics and Automation, Sacramento, April, 9–11, 1696–1701.

Google Scholar

Papadopoulos, E., and Dubowsky, S. (1993). Dynamic Singularities in Free-Floating Space Manipulators. J. Dynamic Syst. Meas. Control 115, 44–52. doi:10.1115/1.2897406

CrossRef Full Text | Google Scholar

Reed, B. B., Smith, R. C., Naasz, B., Pellegrino, J., and Bacon, C. (2016). “The Restore-L Servicing mission,” in AIAA Space Conference, Long Beach California, September, 13–16. doi:10.2514/6.2016-5478

CrossRef Full Text | Google Scholar

Romano, M. (2021). Space Robotics. Cham: Springer International Publishing, 2090–2097. doi:10.1007/978-3-030-44184-5_100035

CrossRef Full Text | Google Scholar

Seddaoui, A. (2020). Precise Motion Control of a Space Robot for In-Orbit Close Proximity Manoeuvres. Ph.D. thesis (Guildford, UK: University of Surrey).

Google Scholar

Seddaoui, A., and Saaj, C. M. (2021). Collision-free Optimal Trajectory Generation for a Space Robot Using Genetic Algorithm. Acta Astronautica 179, 311–321. doi:10.1016/j.actaastro.2020.11.001

CrossRef Full Text | Google Scholar

Seddaoui, A., and Saaj, C. M. (2019). Combined Nonlinear H∞ Controller for a Controlled-Floating Space Robot. J. Guidance, Control Dyn. 42, 1878–1885. doi:10.2514/1.g003811

CrossRef Full Text | Google Scholar

Seddaoui, A., Saaj, C. M., and Eckersley, S. (2019). “Adaptive H Infinity Controller for Precise Manoeuvring of a Space Robot,” in Proc. IEEE International Conference on Robotics and Automation, Montreal, May, 20–24 (Reston, VA: ICRA).

Google Scholar

Shengwei, L. X. P. (2013). On-orbit Capture Technology of Spacecraft. Spacecraft Eng. 1, 242.

Google Scholar

Spong, M. W., and Vidyasagar, M. (2008). Robot Dynamics and Control. United States: John Wiley & Sons.

Google Scholar

Taylor, B., Aglietti, G., Fellowes, S., Ainley, S., Salmon, T., Retat, I., et al. (2018). “Remove Debris mission, from Concept to Orbit,” in 32nd Annual AIAA/USU Conference on Small Satellites, Utah, August, 4–9.

Google Scholar

Umetani, Y., and Yoshida, K. (1989). Resolved Motion Rate Control of Space Manipulators with Generalized Jacobian Matrix. IEEE Trans. Robot. Automat. 5, 303–314. doi:10.1109/70.34766

CrossRef Full Text | Google Scholar

Virgili-Llop, J., Drew, J. V., Zappulla, R., and Romano, M. (2017). Laboratory Experiments of Resident Space Object Capture by a Spacecraft-Manipulator System. Aerospace Sci. Technology 71, 530–545. doi:10.1016/j.ast.2017.09.043

CrossRef Full Text | Google Scholar

Virgili-Llop, J., and Romano, M. (2019). Simultaneous Capture and Detumble of a Resident Space Object by a Free-Flying Spacecraft-Manipulator System. Front. Robot. AI 6, 14. doi:10.3389/frobt.2019.00014

PubMed Abstract | CrossRef Full Text | Google Scholar

Virgili-Llop, J., Zagaris, C., Zappulla, R., Bradstreet, A., and Romano, M. (2019). A Convex-Programming-Based Guidance Algorithm to Capture a Tumbling Object on Orbit Using a Spacecraft Equipped with a Robotic Manipulator. Int. J. Robotics Res. 38, 40–72. doi:10.1177/0278364918804660

CrossRef Full Text | Google Scholar

Whittaker, W., Urmson, C., Staritz, P., Kennedy, B., and Ambrose, R. (2000). “Robotics for Assembly, Inspection, and Maintenance of Space Macrofacilities,” in Proc. AIAA Space Conference and Exposition.

Google Scholar

Wie, B. (2008). Space Vehicle Dynamics and Control. American Institute of Aeronautics and Astronautics.

Google Scholar

Wilde, M., Kwok Choon, S., Grompone, A., and Romano, M. (2018). Equations of Motion of Free-Floating Spacecraft-Manipulator Systems: An Engineer's Tutorial. Front. Robot AI 5, 41. doi:10.3389/frobt.2018.00041

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilde, M., Harder, J., and Stoll, E. (2019). Editorial: On-Orbit Servicing and Active Debris Removal: Enabling a Paradigm Shift in Spaceflight. Front. Robot. AI 6, 136. doi:10.3389/frobt.2019.00136

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, Y. (2021). Adaptive Attitude-Tracking Control of Spacecraft Considering On-Orbit Refuelling. Trans. Inst. Meas. Control. 43, 1298–1309. doi:10.1177/0142331220973132

CrossRef Full Text | Google Scholar

Xu, Y., and Kanade, T. (1992). Space Robotics: Dynamics and Control, 188. Springer Science & Business Media.

Google Scholar

Xu, Y. (1993). “The Measure of Dynamic Coupling of Space Robot Systems,” in IEEE International Conference on Robotics and Automation, Atlanta, May, 2–6.

Google Scholar

Yoshida, K., and Abiko, S. (2002). “Inertia Parameter Identification for a Free-Flying Space Robot,” in AIAA Guidance, Navigation, and Control Conference and Exhibit, California, August, 5–8. doi:10.2514/6.2002-4568

CrossRef Full Text | Google Scholar

Yoshida, K. (2001). “Ets-vii Flight Experiments for Space Robot Dynamics and Control,” in Experimental Robotics VII (Springer), 209–218.

Google Scholar

Keywords: space robot, free-flying, free-floating, controlled-floating, pose control, dynamic modeling

Citation: Seddaoui A, Saaj CM and Nair MH (2021) Modeling a Controlled-Floating Space Robot for In-Space Services: A Beginner’s Tutorial. Front. Robot. AI 8:725333. doi: 10.3389/frobt.2021.725333

Received: 15 June 2021; Accepted: 02 November 2021;
Published: 24 December 2021.

Edited by:

Marcello Romano, Naval Postgraduate School, United States

Reviewed by:

Simone Chesi, Maxar Technologies, United States
Jennifer Hudson, Naval Postgraduate School, United States
Fabio Curti, Sapienza University of Rome, Italy

Copyright © 2021 Seddaoui, Saaj and Nair. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Chakravarthini Mini Saaj , MSaaj@lincoln.ac.uk

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.