Wave-based attitude control of in-orbit spacecraft with large-amplitude slosh

Abstract

This paper aims to reach high-quality attitude control and large-amplitude slosh suppression when a typical liquid-filled spacecraft executes three-axial large-angle maneuvers, by applying a wave-based attitude controller (WBAC). First, the sloshing dynamics is modeled by using a spherical pendulum, whose motion is expressed by splitting its coordinates. Thus, the large-amplitude lateral sloshing and the rotary sloshing as well as the over-all rigid motion of a liquid with respect to the tank in spacecraft can be approximately described. Second, the dynamics equations of system in terms of hybrid coordinate for the liquid-filled spacecraft are derived via the Lagrangian formulation. Third, an improved WBAC for three-axial large-angle maneuvers of in-orbit liquid-filled spacecraft is designed by adding a derivative control law and a gravity-gradient-torque compensation law to a wave-based control law. Simulations of three-axial large-angle maneuver demonstrate good performance of the WBAC in shortening completion time of attitude maneuvering, suppressing the jitter in the angular velocity of spacecraft, and accelerating the large-amplitude slosh suppression, simultaneously. The results indicate the potential applications of the proposed WBAC in the large-angle attitude maneuvering of liquid-filled spacecraft when there exists large-amplitude liquid slosh.

Keywords

liquid-filled spacecraft large-amplitude slosh composite spherical pendulum model wave-based attitude controller three-axial large-angle maneuver

1. Introduction

A new generation of spacecraft usually carries significant amounts of liquid propellant for long-life-span tasks, station keeping, large maneuvers, etc. Especially for three-axis-stabilized liquid-filled spacecraft, the attitude control should adapt to the requirements of high stability and high precision. However, the liquid in the propellant tanks tends to experience nonlinear sloshing motion (mainly including large-amplitude lateral sloshing, nonplanar sloshing, rotary sloshing, as well as the over-all rigid motion with respect to the tank) when the spacecraft executes large and rapid maneuvers (Dodge, 2000; Ibrahim, 2005; Liu et al., 2018, 2020). Then, the real-time evaluation of the generated slosh force and slosh torque acting on the rigid hub of spacecraft will be a very tricky problem in the engineering practice, and the disturbance of liquid slosh may cause the attitude motion of spacecraft accompanied by violent jitter (Liu et al., 2018, 2020). A famous example is the failure in the second flight of the SpaceX “Falcon-1” in 2007, and the cause was recognized as that the propellant sloshing in the liquid oxidizer tank strengthened the circular coning oscillation in the rolling motion of the rocket.

Abundant analytical, experimental, and numerical results suggest that the effect of the linear sloshing of a liquid within a specially shaped container can be accurately evaluated by an equivalent mechanical model (EMM) (Dodge, 2000; Ibrahim, 2005). The commonly used EMMs in engineering, including the mass–spring–damper model and the simple pendulum model, are mainly for describing the in-plane lateral sloshing. Moreover, a spherical pendulum model can be used to predict a liquid which experiences complicated nonplanar motion modes like out-plane lateral sloshing and/or rotary sloshing (Dodge, 2000; Kana, 1987, 1989; Liu et al., 2018; Shageer and Tao, 2008). In recent years, the large-amplitude slosh models were proposed and improved to meet with the demand of the development of spacecraft. Detailedly, a constraint-surface model has been validated to be suitable for the modeling of in-orbit large-amplitude slosh dynamics (Berry and Tegart, 1975; Zhou and Huang, 2015; Liu et al., 2020); a moving pulsating ball model was proposed by Vreeburg (1997), and its practicability has been validated by the in-orbit flight tests of the “Slohsat FLEVO” satellite; a 3DOF-rigid-pendulum model was proposed, on the base of the spherical pendulum model, to investigate more comprehensive nonlinear motion modes of a liquid in spherical propellant tanks (Liu et al., 2019). In recent years, some new numerical computation methods for large-amplitude liquid slosh, being important techniques for the validation of nonlinear EMMs, also have achieved great progress (Huang et al., 2018; Nguyen-Thanh et al., 2019; Tang and Yue, 2017).

The design of the real-time attitude control of liquid-filled spacecraft is almost based on the coupled dynamics equations of spacecraft, which incorporates the liquid sloshing dynamics by using an EMM. Shageer and Tao (2007) gave an overview of the pendulum analogy for fuel sloshing dynamics and developed an adaptive pole placement controller of liquid-filled spacecraft to asymptotically track the desired output. Yue and Zhu (2014) used dynamic inversion and input shaping techniques to design a hybrid controller for the reference tracking maneuvers of spacecraft with fuel slosh, and this hybrid control scheme was demonstrated to be a satisfactory solution of not only to implement the maneuver task but also to suppress the liquid fuel sloshing. Feddema et al. (1997) applied input shaping to suppress the liquid sloshing in a moving hemispherical container, and experiments found that a simple pendulum model is accurate for describing the motion of a liquid in hemispherical containers and that the rapid rotation of a container can induce an unexpected out-of-plane oscillation of liquid. Song and Lu (2019) derived the attitude dynamics of a rigid spacecraft with a partially fuel-filled elliptic tank by using an equivalent two-mode spring–mass model for small-amplitude lateral sloshing and designed an adaptive sliding mode control law with multiple mode input shapers to suppress the fuel sloshing. Thomas et al. (2010) derived a feedback control strategy via the Lyapunov theory for the reorientation of liquid-filled spacecraft, but the angular velocity jitter aroused by the fuel sloshing disturbance and the control efforts can’t be suppressed. However, the control system of three-axis-stabilized spacecraft usually requires attitude maneuvers with high imaging performance, resistance to disturbance, and low jitter. The wave-based control technique was first introduced for position control and vibration suppression of flexible robot arms, by O’Connor and Lang (1998) at the end of the last century. A practical advantage of wave-based control is that all measuring is done inside the control system, and the measurement of the liquid sloshing states is not necessary. Thompson and O’Connor (2016) designed a wave-based controller for in-plane small-angle attitude maneuver of an upper stage accelerating rocket with small-amplitude slosh in a single tank. Colagrossi and Lavagna (2021) improved a three-dimensional real-time MIMO wave-based controller for flexible spacecraft with liquid sloshing replaced by mass–spring–damper slosh model. The above works all consider the control of liquid-filled dynamics system by using small-amplitude slosh models; however, further research is demanded on the three-dimensional attitude control of liquid-filled spacecraft when the liquid is expected to experience more complicated motion such as large-amplitude nonplanar lateral sloshing and rotary sloshing.

The previous work in the literature (Liu et al., 2019) proposed a 3DOF-rigid-pendulum model for describing more comprehensive nonlinear motion of a liquid in spherical propellant tanks, and the nonlinearities of liquid sloshing can be found by studying the model motion. Another previous work in the literature (Liu et al., 2018) has designed an output-feedback attitude controller (OFAC), based on the Lyapunov theory, for three-axial rest-to-rest maneuvers of spacecraft with a partially liquid-filled tank and two flexible solar panels. This kind of controller is capable of suppressing the jitter effect of flex and slosh on the attitude motion of spacecraft; however, the control–structure–slosh interaction problem is likely to be more prominent. Thus, more effective control methods for liquid-filled spacecraft with large-amplitude slosh need to be studied in future research.

In this paper, the large-amplitude lateral sloshing, the nonplanar sloshing, the rotary sloshing, as well as the over-all rigid motion of a liquid in spacecraft are considered by using the composite spherical pendulum model, and the coordinates of the spherical pendulum are split into two parts, one for describing the over-all rigid motion of liquid with respect to the tank, and the other one for describing the sloshing motion of liquid around the dynamic equilibrium position. Then, the rigid-liquid-control coupled dynamics equations are derived via the Lagrangian formulation. Finally, an improved wave-based attitude controller (WBAC) for three-axial large-angle maneuvers of in-orbit liquid-filled spacecraft is designed by adding a derivative control law and a gravity-gradient-torque compensation law to a wave-based control law. Simulations of large-angle rest-to-rest maneuver under the PD controller, the OFAC, and the WBAC are, respectively, given to illustrate the superiority of the WBAC.

2. Physical model of liquid-filled spacecraft

The liquid-filled spacecraft to be described, shown in Figure 1, consists of a rigid hub, a set of reaction wheels for attitude control, and a spherical tank which is partially filled with liquid fuel. Let $N : {n_{1}, n_{2}, n_{3}}$ be a Newtonian frame, and the direction of the local gravitational acceleration $g$ is opposite to $n_{3}$ , namely, $g = - g n_{3}$ . $B^{*}$ is the center of mass (c.m.) of the rigid hub, and the attitude motion of the spacecraft is described by a body-fixed frame $B : {b_{1}, b_{2}, b_{3}}$ , whose axis is the principal axis of inertial of the rigid hub, with regard to $N : {n_{1}, n_{2}, n_{3}}$ .

Figure 1.

Physical model of spacecraft, frames, and variables.

Here, the liquid fuel is approximately replaced by an EMM called the composite spherical pendulum model, which contains a lumped mass $Q$ for the motion description of non-sloshing liquid and a spherical pendulum $P$ to evaluate the mechanical effects of sloshing liquid. $O$ is the hinge point of the pendulum, which is located at the center of the spherical fuel tank, and the layout position of the tank can be given by $p^{B O} = X b_{1} + Y b_{2} + Z b_{3}$ .

The lumped mass $Q$ and the spherical pendulum $P$ can arbitrarily rotate about the hinge point when the spacecraft executes large-angle attitude maneuvers (Liu et al., 2018, 2019). The assignment of the equivalent parameters of the spherical pendulum $P$ , containing the mass parameter $m_{p}$ and the length parameter $l_{p}$ which depend on the radius of the tank $R_{tank}$ and the liquid fill ratio $δ_{liq}$ , can be referred to Dodge’s experimental results (Dodge, 2000). In addition, the mass parameter $m_{0}$ and the length parameter $l_{0}$ of the lumped mass $Q$ can be assigned according to the principles of mechanical equivalence (Ibrahim, 2005; Liu et al., 2019).

3. Attitude dynamics modeling

This work focuses on the coupling dynamics and control during the attitude maneuvering of a liquid-filled spacecraft, so the orbital translation of spacecraft is not considered. The attitude dynamics modeling of the controlled liquid-filled spacecraft system, based on the Lagrangian formulation, covers the motion description of system, the derivation of the system energy, and finally applying the Lagrangian equations.

3.1. Motion description of system

(1) Attitude motion of spacecraft

The attitude of spacecraft is described by the Euler quaternion $ε = {[ε_{0}, ε_{1}, ε_{2}, ε_{3}]}^{T}$ , and the kinematic evolution equation in terms of the quaternion is expressed as

\dot{ε} = \frac{1}{2} S^{T} (ε) {}^{N}{ω^{B}}

(1)

in which

\dot{ε}

is the time derivative of the quaternion,

{}^{N}{ω^{B}}

is the spacecraft angular velocity and

{}^{N}{ω^{B}} = ω_{1} b_{1} + ω_{2} b_{2} + ω_{3} b_{3}

, and

S (ε) = [\begin{array}{c} - ε_{1} & ε_{0} & ε_{3} & - ε_{2} \\ - ε_{2} & - ε_{3} & ε_{0} & ε_{1} \\ - ε_{3} & ε_{2} & - ε_{1} & ε_{0} \end{array}]

Additionally, the quaternion relates to the direction cosine matrix through the following formula

{}^{N}{C^{B}} = (2 {ε_{0}}^{2} - 1) E_{3} + 2 ε_{v} {ε_{v}}^{T} - 2 ε_{0} {\tilde{ε}}_{v}

(2)

in which

E_{3}

is the third-order identity matrix,

ε_{v} = {[ε_{1}, ε_{2}, ε_{3}]}^{T}

is the vector part of the quaternion, and

{\tilde{ε}}_{v}

represents a skew-symmetric tilde matrix operator which is given as

{\tilde{ε}}_{v} = [\begin{array}{c} 0 & - ε_{3} & ε_{2} \\ ε_{3} & 0 & - ε_{1} \\ - ε_{2} & ε_{1} & 0 \end{array}] .

(2) Motion of the slosh model

First, the position vector of the lumped mass $Q$ in frame $N : {n_{1}, n_{2}, n_{3}}$ is $p^{N Q} = p^{N B} + p^{B O} + p^{O Q}$ ; as the orbital translation of spacecraft is not considered here, it can be assumed that the origin of frame $N : {n_{1}, n_{2}, n_{3}}$ locates at $B$ and then $p^{N B} = 0$ . Additionally, as $p^{O Q}$ always points the direction of the local gravitational acceleration, thus $p^{O Q} = - l_{0} n_{3}$ . Then, the velocity of $Q$ in $N : {n_{1}, n_{2}, n_{3}}$ can be found by differentiating $p^{N Q}$ with respect to time in $N : {n_{1}, n_{2}, n_{3}}$

{}^{N}{v^{Q}} = \frac{{}^{N}d}{d t} p^{{N Q}^{p^{N B = 0}}} = \frac{{}^{N}d}{d t} p^{B O} + \frac{{}^{N}d}{d t} (- l_{0} n_{3}) = - p^{B O} \times {}^{N}{ω^{B}}

(3)

Second, the position vector of the lumped mass of the spherical pendulum

P

in frame

N : {n_{1}, n_{2}, n_{3}}

p^{N P} = p^{N B} + p^{B O} + p^{{O P}^{p^{N B = 0}}} = p^{B O} + p^{O P}

(4)

Here, the motion of the spherical pendulum with respect to the tank is described by a set of angle coordinates

ρ = {[φ, ψ]}^{T}

, referring to the previous work in the literature (Liu et al., 2018), and the position vector

p^{O P}

can be given as

p^{O P} = l_{p} (- b_{1} \sin φ \sin ψ + b_{2} \sin φ \cos ψ - b_{3} \cos φ)

(5)

Then, the velocity of $P$ in $N : {n_{1}, n_{2}, n_{3}}$ can be found by differentiating $p^{N P}$ with respect to time in $N : {n_{1}, n_{2}, n_{3}}$

{}^{N}{v^{P}} = \frac{{}^{N}d}{d t} p^{B P} = - (p^{B O} + p^{O P}) \times {}^{N}{ω^{B}} + \frac{{}^{B}d}{d t} p^{O P}

(6)

In order to separately describe the essence of the liquid motion with respect to the tank, including the large-amplitude lateral sloshing, the rotary sloshing, and the over-all rigid motion,

ρ

is splitted into

ρ_{r} = {[φ_{r}, ψ_{r}]}^{T}

for representing the over-all rigid motion of a liquid with respect to the tank (namely, the dynamic equilibrium position) and

ρ_{s} = {[φ_{s}, ψ_{s}]}^{T}

for representing the sloshing motion of the liquid around the dynamic equilibrium position. Then

ρ = ρ_{r} + ρ_{s} = {[φ_{r}, ψ_{r}]}^{T} + {[φ_{s}, ψ_{s}]}^{T}

(7)

in which

ρ_{r}

is determined by the gravitational equilibrium point associated to the attitude of spacecraft, and the solution to

ρ_{r}

follows (Liu et al., 2018)

[\begin{array}{c} \sin φ_{r} \sin ψ_{r} \\ - \sin φ_{r} \cos ψ_{r} \\ \cos φ_{r} \end{array}] = [\begin{array}{c} 2 (ε_{1} ε_{3} - ε_{0} ε_{2}) \\ 2 (ε_{2} ε_{3} + ε_{0} ε_{1}) \\ 2 {ε_{0}}^{2} + 2 {ε_{3}}^{2} - 1 \end{array}]

(8)

3.2. Derivation of system energy

First, the kinetic energy of system, containing that of the attitude motion of the rigid hub, the spin of the reaction control wheels, and the liquid slosh, is defined by

\begin{array}{l} T = \frac{1}{2} \int_{m_{hub}} ({}^{N}{ω^{B}} \times p_{hub}) \cdot ({}^{N}{ω^{B}} \times p_{hub}) d m_{hub} + \frac{1}{2} \sum_{i = 1}^{3} {[I_{W_{i}} \cdot ({}^{N}{ω^{B}} + {}^{B}{ω^{W_{i}}})] \cdot ({}^{N}{ω^{B}} + {}^{B}{ω^{W_{i}}})} \end{array} + \frac{1}{2} m_{0} {}^{N}{v^{Q}} \cdot {}^{N}{v^{Q}} + \frac{1}{2} m_{p} {}^{N}{v^{P}} \cdot {}^{N}{v^{P}}

(9)

in which

d m_{hub}

represents a mass element of the hub, and

p_{hub}

represents the position vector of

d m_{hub}

in frame

N : {n_{1}, n_{2}, n_{3}}

;

I_{W_{i}} = I_{W_{i}} b_{i} b_{i}

is the dyadic of the moment of inertia of the

i_{t h}

reaction wheel; and

{}^{B}{ω^{W_{i}}}

is the spinning velocity of the

i_{t h}

reaction wheel with respect to the body-fixed frame

B : {b_{1}, b_{2}, b_{3}}

, and it can be given as

{}^{B}{ω^{W_{i}}} = Ω_{i} b_{i}

From here on, for the sake of matrix manipulation, vectors and dyadic tensors can be represented by coordinate matrices in the body-fixed frame $B : {b_{1}, b_{2}, b_{3}}$ , and the kinetic energy of system can be written in the following matrix form

T = \frac{1}{2} ω^{T} I_{mb} ω + \frac{1}{2} {(Ω + ω)}^{T} I_{w} (Ω + ω) + \frac{1}{2} {\dot{ρ}}^{T} M_{p} \dot{ρ} + ω^{T} P \dot{ρ}

(10)

in which

ω = {[ω_{1}, ω_{2}, ω_{3}]}^{T}

Ω = {[Ω_{1}, Ω_{2}, Ω_{3}]}^{T}

, and

\dot{ρ} = {[\dot{φ}, \dot{ψ}]}^{T}

;

I_{mb} = I_{hub} + I_{q} + I_{p}

, where

I_{hub} = diag {I_{hub 11}, I_{hub 22}, I_{hub 33}}

is the moment of inertia matrix of the rigid hub about the body-fixed frame

B : {b_{1}, b_{2}, b_{3}}

I_{q} = m_{0} [{\tilde{p}}^{B O}] {[{\tilde{p}}^{B O}]}^{T}

and

I_{p} = m_{p} [{\tilde{p}}^{B P}] {[{\tilde{p}}^{B P}]}^{T}

, and here,

[{\tilde{p}}^{B O}]

and

[{\tilde{p}}^{B P}]

are the skew-symmetric tilde matrices of vectors

p^{B O}

and

p^{B P}

, respectively;

I_{w} = diag (I_{W_{1}}, I_{W_{2}}, I_{W_{3}})

;

M_{p} = diag (m_{p} {l_{p}}^{2}, m_{p} {l_{p}}^{2} \sin^{2} φ)

; and

P

is a

3 \times 2

matrix which denotes the coupling coefficient matrix between the attitude motion and the liquid slosh, and its elements are given by

\begin{array}{c} P_{11} = m_{p} l_{p} (l_{p} \cos ψ - Z \cos φ \cos ψ + Y \sin φ); \\ P_{12} = m_{p} l_{p} (- l_{p} \sin φ \cos φ \sin ψ + Z \sin φ \sin ψ); \\ P_{21} = m_{p} l_{p} (l_{p} \sin ψ - Z \cos φ \sin ψ - X \sin φ); \\ P_{22} = m_{p} l_{p} (l_{p} \sin φ \cos φ \cos ψ - Z \sin φ \cos ψ); \\ P_{31} = m_{p} l_{p} (Y \cos φ \sin ψ + X \cos φ \cos ψ); \\ P_{32} = m_{p} l_{p} (l_{p} \sin^{2} φ + Y \sin φ \cos ψ - X \sin φ \sin ψ) . \end{array}

Second, the potential energy of system is the gravitational potential energy caused by the liquid sloshing motion, and it is given by

V = m_{p} l_{p} g (1 - \cos φ_{s})

(11)

3.3. Lagrangian formula

The Lagrangian equations in terms of quasi-coordinates and generalized coordinates are used to derive the attitude dynamics equation and the sloshing dynamics equation (Meirovitch, 1991). The Lagrangian function of the liquid-filled spacecraft is given by

L = T - V

(12)

First, the Lagrangian formulation for the attitude motion of spacecraft is given as

\frac{d}{d t} (\frac{\partial L}{\partial ω}) + [\tilde{ω}] \frac{\partial L}{\partial ω} = T_{e}

(13)

in which

T_{e}

represents the coordinate matrix in the body-fixed frame

B : {b_{1}, b_{2}, b_{3}}

of total external disturbance torques, and only the gravity gradient torque is considered here, namely,

T_{e} = T_{g}

Second, the Lagrangian formulation for the liquid sloshing motion is given as

\frac{d}{d t} (\frac{\partial L}{\partial \dot{ρ}}) - \frac{\partial L}{\partial ρ} = T_{p}

(14)

in which

T_{p}

represents the coordinate matrix in the body-fixed frame

B : {b_{1}, b_{2}, b_{3}}

of the damping torque acting on the spherical pendulum about the pivot point, and here, it is replaced by a linear viscous torque model

T_{p} = {[- β_{φ} \dot{φ}, - β_{ψ} \dot{ψ}]}^{T}

(15)

in which

β_{φ}

and

β_{ψ}

are constants to assess the damping effect of liquid slosh.

By substituting equations (10) and (11) into equation (12), and then into equations (13) and (14), respectively, the dynamics equations of system can be formulated as

I_{mb} \dot{ω} + P \ddot{ρ} + I_{w} (\dot{Ω} + \dot{ω}) = - [\tilde{ω}] [I_{mb} ω + I_{w} (Ω + ω) + P \dot{ρ}] - {\dot{I}}_{mb} ω - \dot{P} \dot{ρ} + T_{g}

(16a)

M_{p} \ddot{ρ} + P^{T} \dot{ω} = - C_{p} \dot{ρ} - K_{p} (ρ - ρ_{r}) + \frac{1}{2} ω^{T} \frac{\partial I_{mb}}{\partial ρ} ω - \frac{1}{2} {\dot{ρ}}^{T} \frac{\partial M_{p}}{\partial ρ} \dot{ρ}

(16b)

I_{w} (\dot{ω} + \dot{Ω}) = - T_{w}

(16c)

in which

C_{p} = diag {β_{φ}, β_{ψ}}

;

K_{p}

is the generalized stiffness matrix of the slosh pendulum model, and here, the linear sloshing stiffness is adopted, and

K_{p}

can be approximately given as

K_{p} = diag {m_{p} l_{p} g, 0}

;

T_{w} = {[T_{W 1}, T_{W 2}, T_{W 3}]}^{T}

represents the output of the reaction control wheels.

4. Control design

This paper mainly focuses on investigating the attitude control for the three-axial large-angle maneuvers of a liquid-filled spacecraft with large-amplitude slosh, aiming to improve the response performance of the attitude motion of spacecraft, as well as to enhance the active suppression of liquid slosh. Here, the wave-based control method in the literature (Thompson and O’Connor, 2016), for in-plane small-angle attitude maneuver of an upper stage accelerating rocket with small-amplitude slosh, is introduced, and an improved WBAC, which includes a wave-based control law, a derivative control law, and a gravity-gradient-torque compensation law, has been designed for three-axial large-angle maneuvers of in-orbit spacecraft with large-amplitude slosh.

In order to demonstrate the effectiveness of the WBAC, the results under a PD controller and an OFAC will be given for comparison in the following section.

4.1. Wave-based attitude controller (WBAC)

Referring to the literature (Thompson and O’Connor, 2016), a WBAC for single channel attitude control of a spacecraft with a single liquid-filled tank can be designed based on the wave model shown in Figure 2, in which the notional mass $m_{0}$ represents the attitude control system of the spacecraft, $m_{1}$ represents the rigid hub of spacecraft, $m_{2}$ represents the sloshing liquid inside the spacecraft, $B_{0}$ is the returning wave component at the actuator (the attitude control system of the spacecraft) in the Laplace domain, and $G_{0}$ and $H_{0}$ are the transfer functions associated with the control output wave and the returning wave, respectively, and the second-order uniform system approximation are adopted here (O’Connor, 2011; Thompson and O’Connor, 2016)

G_{0} = \frac{{ω_{G}}^{2}}{s^{2} + ω_{G} s + {ω_{G}}^{2}}, ω_{G} = \sqrt{2 k_{0} / m_{1}}

(17a)

H_{0} = \frac{{ω_{H}}^{2}}{s^{2} + ω_{H} s + {ω_{H}}^{2}}, ω_{H} = \sqrt{2 k_{0} / m_{0}}

(17b)

Figure 2.

Wave model for single channel attitude control of the spacecraft.

First, according to the inherent properties of the wave model, the returning wave component $B_{0 i}$ in the $i_{t h} (i = 1,2,3)$ attitude channel can be calculated in the Laplace domain by using the measured values of $ε_{i}$ and $T_{W i}$ (Colagrossi and Lavagna, 2021)

B_{0 i} = H_{0 i} (\frac{ε_{i} - P_{0 i} T_{W i}}{1 - P_{0 i} Q_{0 i}})

(18)

in which

H_{0 i}

P_{0 i}

, and

Q_{0 i}

are the transfer functions for the

i_{t h}

attitude channel, and they are given as

P_{0 i} = \frac{G_{0 i}}{k_{0 i} (1 - G_{0 i})}

(19a)

Q_{0 i} = k_{0 i} (H_{0 i} - 1)

(19b)

G_{0 i} = \frac{{ω_{G i}}^{2}}{s^{2} + ω_{G i} s + {ω_{G i}}^{2}}, ω_{G i} = \sqrt{2 k_{0 i} / m_{1 i}}

(19c)

H_{0 i} = \frac{{ω_{H i}}^{2}}{s^{2} + ω_{H i} s + {ω_{H i}}^{2}}, ω_{H i} = \sqrt{2 k_{0 i} / m_{0 i}}

(19d)

in which the associated physical parameters

m_{1 i}

m_{0 i}

, and

k_{0 i}

should be tuned to reach the desired response of the system under control. For the liquid-filled spacecraft, these control parameters are chosen by

[m_{11}, m_{12}, m_{13}] = [I_{hub 11}, I_{hub 22}, I_{hub 33}]

(20a)

m_{0 i} = δ_{i} m_{1 i}

(20b)

ω_{H i} = \sqrt{2 k_{0 i} / m_{0 i}} = η_{i} ω_{slosh}

(20c)

in which

ω_{slosh} = \sqrt{g / l_{p}}

is the natural angular frequency of the first-order mode of lateral sloshing,

δ_{i}

represents the proportional coefficient of the notional moment of inertia of the control system to that of the rigid hub about the axis

b_{i}

, and

η_{i}

represents the proportional coefficient of the angular frequency of the returning wave in the

i_{t h}

attitude channel to

ω_{slosh}

Second, the wave-based control law is designed as (Thompson and O’Connor, 2016; Colagrossi and Lavagna, 2021)

τ = [\begin{array}{c} k_{01} (0.5 ε_{1}^{tar} + b_{01} - ε_{1}) \\ k_{02} (0.5 ε_{2}^{tar} + b_{02} - ε_{2}) \\ k_{03} (0.5 ε_{3}^{tar} + b_{03} - ε_{3}) \end{array}]

(21)

in which

ε_{i}^{tar}

is the quaternion of the target attitude and

b_{0 i}

is the value in the time domain of

B_{0 i}

Finally, an improved WBAC for the three-axial large-angle maneuvers of in-orbit liquid-filled spacecraft is designed by combining a derivative control law and a gravity-gradient-torque compensation law with the wave-based control law. Thus, the WBAC is given as (the output torque of the reaction wheels)

T_{w 1} = τ - k_{d 1} I_{hub} ω - T_{g}

(22)

4.2. Output-feedback attitude controller (OFAC)

Previous work in the literature (Liu et al., 2018) has designed an OFAC for rigid-flexible-liquid spacecraft, based on the Lyapunov stability theorem. For the liquid-filled spacecraft studied here, this kind of controller (the output torque of the reaction wheels) with a state estimator is given as

\begin{array}{l} [\begin{array}{c} \dot{\hat{ρ}} \\ \dot{\hat{η}} \end{array}] = [\begin{array}{c} 0_{2 \times 2} & E_{2} \\ - M_{p}^{- 1} K_{p} & - M_{p}^{- 1} C_{p} \end{array}] [\begin{array}{c} \hat{ρ} \\ \hat{η} \end{array}] + [\begin{array}{c} - M_{p}^{- 1} P^{T} \\ M_{p}^{- 1} C_{p} M_{p}^{- 1} P^{T} \end{array}] ω + B_{2}^{- 1} ([\begin{array}{c} K_{p} M_{p}^{- 1} P^{T} \\ C_{p} M_{p}^{- 1} P^{T} \end{array}] + B_{1} [\begin{array}{c} M_{p}^{- 1} P^{T} \\ M_{p}^{- 1} C_{p} M_{p}^{- 1} P^{T} \end{array}]) ω + [\begin{array}{c} 0_{2 \times 2} \\ M_{p}^{- 1} K_{p} \end{array}] ρ_{r} \end{array}

(23a)

T_{w 2} = - k_{p 2} I_{hub} ε_{v} - k_{d 2} I_{hub} ω - F [\begin{array}{c} \hat{ρ} \\ \hat{η} \end{array}] + P M_{p}^{- 1} K_{p} ρ_{r} - T_{g}

(23b)

with

F = [- P M_{p}^{- 1}, P M_{p}^{- 1} C_{p} M_{p}^{- 1}] B_{1} + [P M_{p}^{- 1} K_{p}, P M_{p}^{- 1} C_{p}]

in which

η = \dot{ρ} + {M_{p}}^{- 1} P^{T} ω

represents the generalized total velocity of the pendulum;

\hat{ρ}

and

\hat{η}

are the estimations of

ρ

and

η

, respectively;

E_{2}

is the second-order identity matrix; and

B_{1}

and

B_{2}

are the symmetric positive definite (or semi-definite) matrices to be determined, which can be calculated by solving the following matrix equation

B_{i} [\begin{array}{c} 0_{2 \times 2} & E_{2} \\ - M_{p}^{- 1} K_{p} & - M_{p}^{- 1} C_{p} \end{array}] + {[\begin{array}{c} 0_{2 \times 2} & E_{2} \\ - M_{p}^{- 1} K_{p} & - M_{p}^{- 1} C_{p} \end{array}]}^{T} B_{i} = - 2 D_{i}

in which

D_{i} (i = 1,2)

are any fixed symmetric positive definite (or semi-definite) matrices. Here,

D_{1}

is chosen as

D_{1} = [\begin{array}{c} E_{2} & 0_{2 \times 2} \\ 0_{2 \times 2} & K_{p}^{- 1} M_{p} \end{array}]

for determining

B_{1}

, and

B_{2}

is given by

B_{2} = λ B_{1}

, where

λ

is the proportional constant representing the accuracy of the state estimator. For the simulations in Section 5,

λ = 1.0 \times 10^{7}

4.3. PD controller

In addition, the PD controller for this spacecraft can be designed as (the output torque of the reaction wheels)

T_{w 3} = - k_{p 3} I_{hub} ε_{v} - k_{d 3} I_{hub} ω - T_{g}

(24)

5. Numerical simulations and results

In this section, the dynamic responses during the spacecraft undergoing a large-angle rest-to-rest maneuver under different attitude controllers, including the PD controller, the OFAC, and the WBAC, are simulated.

In detail, the spacecraft is maneuvered from the initial attitude described by

ε_{0} (t = 0) = 0.1736, ε_{v} (t = 0) = {[- 0.5264, - 0.2632,0.7896]}^{T}

to the target attitude

ε_{0}^{tar} = 1, ε_{v}^{tar} = 0

, and the attitude of spacecraft is still at the beginning of the maneuver, namely,

ω (t = 0) = 0

. In addition, the initial equilibrium position of the slosh pendulum corresponding to the initial attitude is represented by

ρ_{r} (t = 0) = {[1.2585, - 0.8907]}^{T}

Previous works in the literatures (Liu et al., 2018, 2020) indicate that just small-amplitude liquid slosh will be aroused during the attitude maneuvering when there is no liquid slosh initially, and Tang and Yue (2017) have proved that the damping effect can be ignored for small-amplitude slosh but not for large-amplitude slosh case. Thus, the following two cases are, respectively, simulated, and the results under the three different controllers are compared here. The parameters of the spacecraft system for numerical simulations are contained in Table 1.

(1) Case A: there is no initial liquid slosh, namely, $ρ (t = 0) = ρ_{r} (t = 0) = {[1.2585, - 0.8907]}^{T}$ , and the damping effect of liquid slosh is ignored by setting $β_{φ} = β_{ψ} = 0$ .

(2) Case B: there is initial large-amplitude liquid slosh by setting $ρ (t = 0) = {[1, - 1]}^{T}$ , and the damping effect of liquid slosh is considered by setting $β_{φ} = β_{ψ} = 5 kg \cdot m^{2} \cdot s^{- 1}$ .

Table 1.

Parameters of the system.

Symbol	Value
$I_{hub}$	$diag {1500,3000,3000} (kg \cdot m^{2})$
$I_{w}$	$diag {20,20,20} (kg \cdot m^{2})$
$(X, Y, Z)$	$(0,0, - 0.5) (m)$
$R_{tank}$	$0.6 (m)$
$δ_{liq}$	50%
$ρ_{liq}$	$1000 (kg \cdot m^{- 3})$
$(m_{p}, m_{0})$	$(273.7,178.7) (kg)$
$(l_{p}, l_{0})$	$(0.393, - 0.0323) (m)$
$g$	$7.689 (m \cdot s^{- 2})$

5.1. Results of Case A

In this case, the damping effect of liquid slosh is ignored as the amplitude of liquid slosh is small. The liquid-filled spacecraft is actually a non-damping system, and the WBAC tends to achieve the ideal control effect when the angular frequency of the returning wave is tuned to be equal to the natural angular frequency of liquid slosh, namely,

ω_{H i} = ω_{slosh}

. The control parameters of the WBAC and the PD controller are contained in Table 2.

Table 2.

Parameters of the controllers for Case A.

Controller	Symbol	Value
WBAC	$(δ_{1}, δ_{2}, δ_{3}; η_{1}, η_{2}, η_{3}; k_{d 1})$	$(0.002,0.002,0.002; 1,1,1; 0.2)$
PD controller	$(k_{p 3}, k_{d 3})$	$(0.1,1)$

By comparing with the PD controller, as shown in Figures 3 and 4, the WBAC significantly lowers the requirement of the maximum output capacity of control system (the maximum control torque of the WBAC is just about 20% of that of the PD controller), simultaneously ensuring almost the same completion time of large-angle attitude maneuvering (the completion time of the two controllers is 100–120 s). In this case, the induced liquid slosh is very slight and also the jitter in the attitude motion of the spacecraft (Figure 3).

Figure 3.

Angular velocity of spacecraft of Case A.

Figure 4.

Attitude control torque of Case A.

5.2. Results of Case B

In this large-amplitude slosh case, the damping effect of liquid slosh is included by setting

β_{φ} = β_{ψ} = 5 kg \cdot m^{2} \cdot s^{- 1}

, and the WBAC brings the spacecraft system to a desired response when the returning wave frequency is tuned to

ω_{H i} = 1.75 ω_{slosh}

. The related parameters of the WBAC, the OFAC, and the PD controller are contained in Table 3.

Table 3.

Parameters of the controllers for Case B.

Controller	Symbol	Value
WBAC	$(δ_{1}, δ_{2}, δ_{3}; η_{1}, η_{2}, η_{3}; k_{d 1})$	$(0.005,0.005,0.003; 1.75,1.75,1.75; 0.2)$
OFAC	$(k_{p 2}, k_{d 2}; λ)$	$(0.1,1; 1.0 \times 10^{7})$
PD controller	$(k_{p 3}, k_{d 3})$	$(0.1,1)$

Figure 5 shows the time history responses of the angular velocity of spacecraft under the three different controllers with the same level of maximum output torque (Figure 6). Comparing with the result of the PD controller, the OFAC performs good capacity to reduce the jitter effect of liquid slosh on the attitude motion at the expense of greater control efforts (Figure 6). However, the WBAC shows potential advantages in shortening completion time of attitude maneuvering, suppressing the jitter in the angular velocity of spacecraft (Figure 5), simultaneously accelerating the large-amplitude slosh suppression (Figure 7). Figure 7 shows the liquid has experienced large-amplitude lateral sloshing and finite-angle rotary sloshing during the large-angle attitude maneuvering of spacecraft, and the liquid slosh will be completely dissipated due to the damping effect. The results indicate the potential applications of the proposed WBAC in the large-angle attitude maneuvering of liquid-filled spacecraft when there exists large-amplitude liquid slosh.

Figure 5.

Angular velocity of spacecraft of Case B.

Figure 6.

Attitude control torque of Case B.

Figure 7.

Position coordinates of the slosh pendulum of Case B.

6. Conclusions

The attitude control for three-axial large-angle maneuvers of in-orbit spacecraft with large-amplitude slosh has been studied in this paper. First, a systemized attitude dynamics modeling process of liquid-filled spacecraft was implemented by using the Lagrangian formulation, where the large-amplitude liquid slosh dynamics are equivalently modeled by adopting the composite spherical pendulum model. Then, an improved WBAC, including a wave-based control law, a derivative control law, and a gravity-gradient-torque compensation law, was designed for high-quality control of three-axial large-angle maneuvers of in-orbit liquid-filled spacecraft. In practice, the measurements of the feedback information of the WBAC, including the attitude quaternion and angular velocity of spacecraft and the actual achieved torque, are all done at the control system. Comparing with the PD controller and the OFAC, simulation results show that the WBAC has advantages in shortening completion time of attitude maneuvering, suppressing the jitter in the angular velocity of spacecraft, and accelerating the large-amplitude slosh suppression, simultaneously.

Footnotes

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Natural Science Foundation of China [grant numbers 11772049, 12132002, and 12202044] and the 13th Five-Year Civil Space Technology Advance Research Project of the Administration of Science, Technology and Industry [grant number D020201].

ORCID iD

Feng Liu

Appendix

References

Berry

Tegart

(1975) Experimental Study of Transient Liquid Motion in Orbiting Spacecraft. NASA Report NAS8-30690. Huntsville, Alabama: Marshall Space Flight Center.

Colagrossi

Lavagna

(2021) Integrated vibration suppression attitude control for flexible spacecraft with internal liquid sloshing. Multibody System Dynamics 51: 123–157.

Dodge

(2000) The New Dynamic Behavior of Liquids in Moving Containers. San Antonio TX: Southwest Research Institute.

Feddema

Dohrmann

Parker

, et al. (1997) Control for slosh-free motion of an open container. IEEE Control Systems 17(1): 29–36.

Huang

Nguyen-Thanh

Zhou

(2018) An isogeometric-meshfree coupling approach for contact problems by using the third medium method. International Journal of Mechanical Sciences 148: 327–336.

Ibrahim

(2005) Liquid Sloshing Dynamics: Theory and Applications. Cambridge England, UK: Cambridge University Press.

Kana

(1987) Validated spherical pendulum model for rotary liquid slosh. Journal of Spacecraft and Rockets 26(3): 169–177.

Kana

(1989) A model for nonlinear rotary slosh in propellant tanks. Journal of Spacecraft and Rockets 24(2): 188–196.

Liu

Yue

Zhao

(2018) Attitude dynamics and control of spacecraft with a partially filled liquid tank and flexible panels. Acta Astronautica 143: 327–336.

10.

Liu

Yue

Tang

, et al. (2019) 3DOF-rigid-pendulum analogy for nonlinear liquid slosh in spherical propellant tanks. Journal of Sound and Vibration 460: 114907.

11.

Liu

Yue

Banerjee

, et al. (2020) Large motion dynamics of in-orbit flexible spacecraft with large amplitude propellant slosh. Journal of Guidance, Control, and Dynamics 43(3): 438–450.

12.

Meirovitch

(1991) Hybrid state equations of motion for flexible bodies in terms of quasi-coordinates. Journal of Guidance, Control, and Dynamics 14(5): 1008–1013.

13.

Nguyen-Thanh

Huang

, et al. (2019) An adaptive isogeometric analysis meshfree collocation method for elasticity and frictional contact problems. International Journal for Numerical Methods in Engineering 120(2): 209–230.

14.

O’Connor

Lang

(1998) Position control of flexible robot arms using mechanical waves. Journal of Dynamic Systems, Measurement, and Control 120(3): 334–339.

15.

O’Connor

(2011) Wave-like modelling of cascaded, lumped, flexible systems with an arbitrarily moving boundary. Journal of Sound and Vibration 330(13): 3070–3083.

16.

Shageer

Tao

(2007) Modeling and adaptive control of spacecraft with fuel slosh: overview and case studies. In: AIAA Guidance, Navigation and Control Conference and Exihibit, 2007, Hilton Head, South Carolina, 20–23 August 2007. AIAA Paper 2007-6434.

17.

Shageer

Tao

(2008) Zero dynamics analysis for spacecraft with fuel slosh. In: AIAA Guidance, Navigation and Control Conference and Exihibit, 2008, Honolulu, Hawaii, 18–21 August 2008. AIAA Paper 2008-6455.

18.

Song

(2019) Attitude maneuver control of liquid-filled spacecraft with unknown inertia and disturbances. Journal of Vibration and Control 25(8): 1460–1469.

19.

Tang

Yue

(2017) Simulation of large-amplitude three-dimensional liquid sloshing in spherical tanks. AIAA Journal 55(6): 2052–2059.

20.

Thomas

Kang

Victoria

(2010) Analytical control law for spacecraft reorientation via lyapunov theory. In: AIAA Guidance, Navigation and Control Conference, 2010, Toronto: Canada, 02–05 August 2010. AIAA Paper 2010-7894.

21.

Thompson

O’Connor

(2016) Wave-based attitude control of spacecraft with fuel sloshing dynamics. Archive of Mechanical Engineering 63(2): 263–275.

22.

Vreeburg

JPB

(1997) Dynamics and control of a spacecraft with a moving pulsating ball in a spherical cavity. Acta Astronautica 40(2): 257–274.

23.

Yue

Zhu

(2014) Hybrid control of liquid-filled spacecraft maneuvers by dynamic inversion and input shaping. AIAA Journal 52(3): 618–626.

24.

Zhou

Huang

(2015) Constraint surface model for large amplitude sloshing of the spacecraft with multiple tanks. Acta Astronautica 111: 222–229.