Next Article in Journal / Special Issue
Self-Organizing Control-Loop Recovery for Predictive Networked Formation Control of Fractionated Spacecraft
Previous Article in Journal
Minimisation of Failure Transients in a Fail-Safe Electro-Mechanical Actuator Employed for the Flap Movables of a High-Speed Helicopter-Plane
Previous Article in Special Issue
Implicit Extended Kalman Filter for Optical Terrain Relative Navigation Using Delayed Measurements
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Angular Trajectory Design for MR-SPS Using Bezier Shaping Approach

Department of Aerospace Engineering, Harbin Institute of Technology, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Aerospace 2022, 9(10), 528; https://doi.org/10.3390/aerospace9100528
Submission received: 1 July 2022 / Revised: 12 September 2022 / Accepted: 13 September 2022 / Published: 20 September 2022
(This article belongs to the Special Issue Recent Advances in Spacecraft Dynamics and Control)

Abstract

:
Solar power satellite (SPS) is a kind of large-scale on-orbit servicing spacecraft collecting solar energy in space and transmitting energy to the earth. The solar arrays of the SPS must point to the sun to collect enough solar energy and the antenna must point to the rectenna on the ground to transmit energy. But due to the limitation of the control effort, accurate solar and earth orientation may not be achieved. This paper focuses on the MR-SPS, and establishes the attitude kinematics and dynamics model of a MW-level MR-SPS. Angular trajectory based on Bezier shaping approach is generated at different time of a year to satisfy the control constraints. The simulation results demonstrate the effectiveness of the proposed method. Even if the control torque is limited to a small amount, the optimal angular trajectory can still ensure high average energy receiving efficiency.

1. Introduction

With the increasing demand for energy and the depletion of non-renewable energy, the world is facing a serious energy crisis. As a kind of abundant clean and renewable energy, the rational utilization of solar energy is of great significance to solve the energy crisis and environmental problems. Solar power satellite (SPS) also called space power station or space solar power satellite (SSPS) first proposed by Glaser in 1968 [1] is a class of large-scale on-orbit servicing spacecraft that can collect solar energy in space and transmit electrical energy in the form of microwave to the earth. A typical SPS consists of three main parts: the solar collector uses the solar arrays to receive the sunlight and convert the solar energy into electric energy, some SPSs are equipped with concentrators such as lenses or concave mirror to converge the sunlight; the conversion device converts the electric energy into microwave or laser; the antenna or laser transmitter transmits the microwave or laser to the ground, in which the rectenna of the receiving station receives the microwave or laser and converts it into electric energy which can be injected to power supply system.
Over the past couple of decades, the United States, Europe, Japan, and China have put forward a variety of construction schemes for SPSs. Since Glaser first proposed the idea of building a SPS in geostationary orbit (GEO), in 1970s, NASA continued to explore the feasibility and key technologies of SPS and proposed the 1979 SPS reference system [2]. In 1990s, NASA successively proposed a series of SPS configurations, including Sun Tower, Solar Disc, Abacus and Integrated Symmetrical Concentrated system (ISC) [3]. In 2011, John Mankins came up with the solar power satellite via arbitrarily large-phased array (SPS-ALPHA) [4]. ESA put forward the Sail Tower SPS, which effectively reduced the system mass but could not achieve sustainable energy supply [5]. JAXA began the study of SPS in 1990 and several configurations were given, for example, JAXA 2001, JAXA 2002, and Tethered SPS [6], etc. Tethered SPS has smaller mass and can supply energy continuously. China has carried out research on SPS since 2015 based on two types of SPS concept: Multirotary-joint Solar Power satellite (MR-SPS) by CAST [7] and Space Solar Power Station via Orb-shape Membrane Energy Gathering Array (SSPS-OMEGA) by Xidian University [8].
The core task of a SPS is to collect as much solar energy as possible and transmit it to the ground accurately. To accomplish the above goal, attitude of an SPS is supposed to be adjusted so that the solar collector tracks the sun and at the same time, the antenna or laser transmitter points to the rectenna in the ground accurately. Due to its large mass and inertia, large torque is required for attitude control of a SPS, which means more actuators and fuel need to be carried and further increases the mass of the system. In order to reduce fuel consumption and mass, some research has been done. [9] proposed a novel quasi-Sun-pointing attitude in Sun-frozen orbit. Although about 3% electricity must be given up, little control effort is required to deal with the solar radiation pressure and gravity-gradient torque. Refs. [10,11] developed an integrated orbit, attitude, and structural control system architecture for very large SSPS in GEO. A low-bandwidth attitude control system was proposed utilizing cyclic-disturbance accommodating control to provide precision pointing of the Abacus platform in the presence of dynamic modeling uncertainties and external disturbances. Moreover, the configuration of the actuators and the propellants required were also analyzed in detail. In [12,13], a detailed study of solar power satellites’ orbit dynamics was performed, and a direct comparison between the attitude dynamics of SPS in GLP orbit and in GEO was made. When both attitude and orbit control are considered, fuel consumption is less in GLP compared to in GEO.
MR-SPS is an SPS scheme located in GEO proposed by CAST, which utilizes a number of mutually independent rotating solar arrays continuously oriented to the sun to receive solar energy and converts it into electric energy that is transmitted to the power transmission bus through a number of independent low-power conductive rotating joints. The transmission medium of the energy is microwave.
There was some research on the attitude of the MR-SPS. In [14], the flexible multibody attitude dynamic model including external disturbances is derived, and a hybrid high/low bandwidth robust controller is proposed to achieve the attitude control of the SPS. Refs. [15,16] proposed a switched iterative learning controller (ILC) for a flexible SPS to periodically track the earth and the sun. The previous research established the dynamic model of MR-SPS based on analytical mechanics which need complicated derivatives with respect to generalized coordinates. The angular trajectory was not optimized which can cause quite large control torque.
In view of the dynamic characteristics of MR-SPS, if the microwave antenna is accurately oriented to the ground and the solar arrays are accurately oriented to the sun when the energy receiving efficiency is 100%, there will be a slew at some time around the vernal and autumnal equinox, which requires large control torques. However, subject to the total weight and degradation of the actuators, the torque generated by the actuators may fail to meet the orientation requirements. The previous attitude planning for SPS did not consider the constraints of dynamics. Aiming at this problem, an MR-SPS dynamic model was established. Using the method of Bezier shaping approach, the angular trajectory is optimized under the condition of limited control torque, which ensures the accurate earth orientation of the antenna and the quasi sun orientation of solar arrays as well as a high energy receiving efficiency. Simulation verifies the feasibility of the optimization method.
The rest of this paper is arranged as follows: In Section 2, the model of MR-SPS is presented; the angular trajectory planning based on Bezier shaping approach is described in Section 3; in Section 4, the simulation results are given. The results are summarized in Section 5.

2. Model Description

2.1. MR-SPS Structure

In terms of structure, MR-SPS can be divided into three parts, as shown in Figure 1. For a MW-level MR-SPS, there are two solar arrays in each side in the north and the south direction, each solar array is composed of 10 solar subarrays, the size of each subarray is 30 m × 100 m, and the distance between two subarrays is 2 m. The main structure is a truss structure, which consists of the upper and lower north-south main truss structure and the longitudinally connected main truss structure, with a total length of 800 m and width of 51 m. The upper north-south main truss structure is used to support the microwave transmitting antenna that is fixed to the main structure. The longitudinally connected main truss structure connects the north-south main structures together. The main structure and the antenna form the central body structure. The diameter of the microwave transmitting antenna is 140 m.
It is assumed that the structural flexibility is not considered. Since the attitude of the 10 subarrays in each side is consistent, they can be considered to be fixed to each and form a big array. Torque actuators mounted on the central body generate three-dimensional control torque, and motors between the solar arrays and the main structure generate control torque only along the ya-axis.

2.2. Sun and Earth Orientation Analysis

Regardless of the dynamic constraints, the real-time and precise orientation to the sun for the arrays and to the rectenna for the antenna can be achieved by adjusting the attitude of the MR-SPS and in this case, the energy receiving efficiency reaches the peak. The energy receiving efficiency is defined as the area weighted average of the cosine of the angle between the normal direction of the arrays and the solar direction.
As shown in Figure 2, five kinds of reference frame are introduced to describe the attitude of the MR-SPS.
  • The origin oi of the geocentric inertial reference frame oi-xiyizi is located at the center of mass of the earth, the xi-axis is defined along the direction of vernal equinox, the zi-axis is normal to the equatorial plane and points to the north pole, and the yi-axis completes a right-handed coordinate system;
  • The origin oo of the orbital reference frame oo-xoyozo is at the center of mass of the main structure and antenna, the zo-axis points to the center of the earth, yo-axis is consistent with the negative normal direction of the orbital plane, and the xo-axis completes a right-handed system;
  • The origin oa of the main structure body-fixed reference frame oa -xayaza coincides with oo, the three axes are aligned with the principal axes of inertia, the ya-axis is along the longest side of the central truss, and the za-axis is perpendicular to the antenna surface.
  • The origin ocj of the jth intermediate reference frame ocj-xcjycjzcj is the hinged point of the jth solar array and the three axes are parallel to those of the main structure frame;
  • The origin osj of the jth array body-fixed reference frame osj-xsjysjzsj is at the center of mass of the jth solar array, the three axes are aligned with the principal axes of inertia, and ysj-axis is parallel to ya-axis.
φ is the right ascension of the SPS, α is the angle that the sun runs on the ecliptic from the vernal equinox, and β is the obliquity of the ecliptic. Limited by the area of the rectenna and the distance between the SPS and the rectenna, the antenna must point to the rectenna precisely, and in other words, the za-axis must be parallel to the zo-axis at any time. Because the direction of the za-axis is determined, only one-axis motion is considered for the main structure. After rotating γ around zo-axis, oo-xoyozo coincides with oa-xayaza. θj is the rotation angle of the jth array around the main truss from the initial position where osj-xsjysjzsj coincides with ocj-xcjycjzcj. x ^ denotes the unit vector of the vector x or the x-axis. The position unit vector of the sun relative to the earth center is r ^ s = [ cos α sin α cos β sin α sin β ] T , and the unit vector from the earth center to the centroid of the main structure and antenna is r ^ a = [ cos φ sin φ 0 ] T . The two unit vectors are both expressed in the inertial frame. The zo-axis is parallel to the za-axis, so z ^ a = r ^ a . For a spacecraft in GEO, it is obvious that y ^ o = [ 0 0 1 ] T and x ^ o = [ sin φ cos φ 0 ] T expressed in the inertial frame. In fact, as long as γ and θj are confirmed, the attitude of the MR-SPS is determined, and according to the geometric relationship, γ and θj can be calculated by
γ = arctan ( x ^ o y ^ a y ^ 0 y ^ a )
θ j = arctan ( x ^ a z ^ s j z ^ a z ^ s j )
For the purpose of precise solar orientation of the arrays and earth orientation of the antenna, the ya-axis needs to be vertical to r ^ s and r ^ a simultaneously and the zsj-axis needs to be parallel to r ^ s which means z ^ s j = r ^ s . In view of the symmetrical structure of the SPS, y ^ a and x ^ a can be obtained as follows:
y ^ a = ± r ^ s × r ^ a r ^ s × r ^ a
x ^ a = y ^ a × z ^ a
Substitute Equations (3) and (4) into Equations (1) and (2) and in this case, the ideal angular trajectory is obtained and the energy receiving efficiency reaches to 100%.

2.3. Attitude Dynamic Model

The Newton–Euler approach based on the balance of all the forces and torques acting on each part of the structure is widely used to solve inverse dynamics of robots. The method avoids complicated derivatives with respect to generalized coordinates, which is common with Lagrangian-based approaches [17]. It can be adopted to model the complex dynamics of the MR-SPS.
The force balance equations of the central body and the jth (j = 1, 2) solar array are listed as follows:
F s o j + F s i j = m s j C s j C a a s j = m s j C s j C a ( r ¨ a + C a T ω ˙ a × r c j + C a T ω a × ω a × r c j )
F a o + k = 1 2 F a i k = m a C a a a = m a C a r ¨ a
F s i j = C s j F a i j
where ma and msj are the mass of the central body and jth solar array respectively. Csj is the rotation matrix from the central body frame to the jth solar array frame, and Ca denotes the orientation of the central body frame respect to the inertial frame. Fao and Fsoj are external forces acting on the central body and the jth solar array frame. Faij and Fsij are the forces exerted on the central body by the jth solar array and on the jth solar array by the central body. rais the position vector of the central body and rsj is the vector from the centroid of the central body to the centroid of the jth solar array expressed in the inertial frame. ωa denotes the angular velocity of central body respect to inertial frame expressed in the central body frame. The notation x × for x = [ x 1 x 2 x 3 ] T represents the skew symmetric matrix:
x × = 0 x 3 x 2 x 3 0 x 1 x 2 x 1 0
According to Equations (5)–(7), forces exerted on the central body by the jth solar array are given as follows:
F a i j = A o j + A ω ˙ a j ω ˙ a + A o t h e r j
where
A o j = 1 m a + k = 1 2 m s k ( ( m a + k = 1 n m s k ) C s j T F s o j m s j k = 1 2 C s k T F s o k m s j F a o ) A ω ˙ a j = m s j m a + k = 1 2 m s k ( k = 1 2 ( m s k r c k ) + ( m a + k = 1 n m s k ) r c j ) × A o t h e r j = m s j ω a × ω a × m a + k = 1 2 m s k ( k = 1 2 ( m s k r c k ) ( m a + k = 1 n m s k ) r c j )
The external forces acting on the MR-SPS include gravity, solar pressure, and microwave radiation force. The magnitude of solar pressure and microwave radiation pressure is far less than that of gravity, so they can be ignored. The external can be expressed as follows:
F a o = μ m a r a 3 C a r a
F s o j = μ m s j r a + C a T r c i 3 ( C s j C a r a + C s j r c j )
where μ denotes the gravitational coefficient of the earth.
The torque balance equations of the central body and the jth (j = 1,2) solar array are listed as follows:
T c + T a o + k = 1 2 ( T a i k + r c k × F a i k ) = I a ω ˙ a + ω a × I a ω a
T s o j + T s i j = I s j ( ω s j × C s j ω a + C s j ω ˙ a + ω ˙ s j ) + ( C s j ω a + ω s j ) × I s j ( C s j ω a + ω s j )
T s i j = C s j T a i j
τ s i j = D j T T s i j
where I a = diag ( I a 1 , I a 2 , I a 3 ) denotes the moment of inertia of the central body relative to the center of mass, and I s j = diag ( I s j x , I s j y , I s j z ) denotes the moment of inertia of the jth solar array relative to the center of mass. Tc and τscj are the control torques of the central body and the jth solar array. Taoj and Tsoj are external disturbance torques acting on the central body and the jth solar array. Taij and Tsij are the torques exerted on the central body by the jth solar array and on the jth solar array by the central body. ω s j = D j θ ˙ j denotes the angular velocity of the jth solar array respect to the central body frame expressed in jth solar array frame, where D j = [ 0 1 0 ] T represents the unit vector of the direction of the rotation axis of the jth solar array and θj represents the angle that the jth solar array rotates. Csj can be written as:
C s j = cos θ j 0 sin θ j 0 1 0 sin θ j 0 cos θ j
The attitude dynamics are expressed in a form convenient for the subsequent research on attitude control. Substituting Equations (8), (13) and (14) into Equations (11) and (12), the attitude dynamics can be expressed as follows:
B ω ˙ a 0 ω ˙ a + B θ ¨ 1 0 θ ¨ 1 + B θ ¨ 2 0 θ ¨ 2 + B o t h e r 0 = T c + B o 0
B ω ˙ a j ω ˙ a + B θ ¨ j j θ ¨ j + B o t h e r j = τ s c j + B o j
where
B o 0 = T a o + k = 1 n ( r c k × A o k + C s k T T s o k ) B o j = D j T T s o j B ω ˙ a 0 = I a + k = 1 n ( C s k T I s k C s k r c k × A ω ˙ a k ) B ω ˙ a j = D j T I s j C s j B θ ¨ j 0 = C s j T I s j D j B θ ¨ j j = D j T C s j D j B other 0 = ω a × I a ω a + k = 1 n C s k T ( I s k ω s k × C s k ω a + ( C s k ω a + ω s k ) × I s k ( C s k ω a + ω s k ) r c k × A other k ) B other j = D j T C s j ( I s j ω s j × C s j ω a + ( C s j ω a + ω s j ) × I s j ( C s j ω a + ω s j ) )
The external torques acting on the MR-SPS include gravity gradient torque, solar pressure torque, and microwave radiation torque caused by the deviation of centroid and pressure center. Due to the symmetrical structure, the deviation is small, so the magnitude of solar pressure torque and microwave radiation torque is far less than that of gravity gradient torque and cannot be known in advance, so they can be ignored. The external torques can be expressed as follows:
T a o = 3 μ r a 5 ( C a r a ) × I a ( C a r a )
T s o j = 3 μ C a r a + r c j 5 ( C s j C a r a + C s j r c j ) × I s j ( C s j C a r a + C s j r c j )

3. Trajectory Planning for the MR-SPS Considering Dynamic Constraints

3.1. State Approximation Using Bezier Curve

In Bezier shaping approach, it is assumed in advance that the angular trajectory of the MR-SPS follows Bezier curve functions. The approximation of γ and θj can be expanded by applying Bezier curve functions. To avoid repetition, only the processing approach for γ is given, and θj can be processed in the same way. The approximation of γ is expanded as follows:
γ ( τ ) = k = 0 n γ B γ , k ( τ ) P γ , k
where τ = t / T   ( 0 τ 1 ) is the scaled time, and T is the flight time; nγ is the order of the Bezier curve function. P γ , k   ( k [ 0 , n γ ] ) are the Bezier coefficients; B γ , k τ are the Bezier basis functions given by:
B γ , k ( τ ) = n γ ! k ! ( n γ k ) ! τ k ( 1 τ ) n γ k
The first and second τ-derivative of γ can be derived as follows:
γ ( τ ) = k = 0 n γ B γ , k ( τ ) P γ , k
γ ( τ ) = k = 0 n γ B γ , k ( τ ) P γ , k
where B γ , k ( τ ) and B γ , k ( τ ) are first and second τ-derivative of the Bezier basis functions that can be derived according to Equation (20) as:
B γ , k ( τ ) = n γ ( 1 τ ) n γ 1 k = 0 n γ ! ( k 1 ) ! ( n γ k ) ! τ k 1 ( 1 τ ) n γ k n γ ! k ! ( n γ k 1 ) ! τ k ( 1 τ ) n γ k 1 k [ 1 , n γ 1 ] n γ τ n γ 1 k = n r
B γ , k ( τ ) = { n γ ( n γ 1 ) ( 1 τ ) n γ 2 k = 0 n γ ( n γ 1 ) ( n γ 2 ) τ ( 1 τ ) n γ 3 2 n γ ( n γ 1 ) ( 1 τ ) n γ 2 k = 1 n γ ! ( k 2 ) ! ( n γ k ) ! τ k 2 ( 1 τ ) n γ k 2 n γ ! ( k 1 ) ! ( n γ k 1 ) ! τ k 1 ( 1 τ ) n γ k 1 + n γ ! k ! ( n γ k 2 ) ! τ k ( 1 τ ) n γ k 2 k [ 2 , n γ 2 ] n γ ( n γ 1 ) ( n γ 2 ) τ n γ 3 ( 1 τ ) 2 n γ ( n γ 1 ) τ n γ 2 k = n γ 1 n γ ( n γ 1 ) τ n γ 2 k = n γ
The start time is selected as when the right ascension between the sun and the SPS is about 90° because the value of γ is equal to the inclination of the sun at equatorial plane and γ ˙ is equal to 0. After a solar day (86400s) instead of the orbital period considering the movement of the sun, γ and γ ˙ are almost equal to the ones at start time. So, T is chosen as a solar day. To ensure the convergence, the following boundary conditions are required to be satisfied:
γ ( τ = 0 ) = γ 0 γ ( τ = 1 ) = γ f γ ( τ = 0 ) = T γ ˙ 0 γ ( τ = 1 ) = T γ ˙ f
where x ˙ denotes the derivative of x with respect to time t; x denotes the derivative of x with respect to scaled time τ. γ0, γf, γ ˙ 0 , and γ ˙ f can be calculated according to Equations (1)–(4) so that the approximation of γ  γ ˙ is continuous after a solar day.
Substituting τ = 0 and τ = 1 into Equations (20)–(22) and (26), four Bezier coefficients P γ , 0 , P γ , 1 , P γ , n γ 1 , and P γ , n γ can be determined as follows:
P γ , 0 = γ 0 P γ , 1 = γ 0 + T γ ˙ 0 / n γ P γ , n γ = γ f T γ ˙ f / n γ P γ , n γ = γ f
In order to meet the dynamic constraints at every moment, some discrete points can be adopted to calculate the dynamic constraints. The discrete points can be obtained by applying Gaussian Legendre distribution, which is defined as the roots of m-order Legendre polynomial to be scaled to [ 0 , 1 ] and to be arranged in ascending order as follows: [ τ ] m × 1 = [ τ 1 τ 2 τ m 1 τ m ] T . The discrete points of γ and its first and second τ –derivatives can be given in compact matrix notation form as follows:
[ γ ] m × 1 = [ B γ ] m × ( n γ + 1 ) [ P γ ] ( n γ + 1 ) × 1
[ γ ] m × 1 = [ B γ ] m × ( n γ + 1 ) [ P γ ] ( n γ + 1 ) × 1
[ γ ] m × 1 = [ B γ ] m × ( n γ + 1 ) [ P γ ] ( n γ + 1 ) × 1
where [ P γ ] ( n γ + 1 ) × 1 = [ P γ , 0 P γ , 1 [ X γ ] ( n γ 3 ) × 1 T P γ , n γ 1 P γ , n γ ] T is the Bezier coefficient column matrix and [ X γ ] = [ P γ , 2 P γ , n r 1 ] T is the unknown coefficient column matrix to be optimized. [ B γ ] m × ( n γ + 1 ) , [ B γ ] m × ( n γ + 1 ) , and [ B γ ] m × ( n γ + 1 ) are Bezier basis function matrix and the first and second τ-derivative matrices of the Bezier basis function which can be calculated by substituting [ τ ] m × 1 into Equations (21), (24) and (25).
Once the order of Bezier curve and the number of discrete points are determined, [ B γ ] m × ( n γ + 1 ) , [ B γ ] m × ( n γ + 1 ) , and [ B γ ] m × ( n γ + 1 ) are determined, only calculated once and stored. In the process of optimization, they can be used in every iteration without repeated calculation, which reduce the calculation burden of the computer.

3.2. Description of the Optimization Problem

As a large-scale spacecraft, the MR-SPS requires large torque in the attitude tracking process. During a short period of time around the time of the equinox, if the arrays strictly point to the sun at any moment, there will be a large angle of ration in a short time for the main structure when the right ascension between the sun and the SPS is about 0° or 180°, resulting in large control effort. From the perspective of economic benefits, the large-scale spacecraft needs to limit the mass. The larger the maximum control torque is, the heavier the actuators are which leads to an increase in the mass of actuators as well as the launch cost. Moreover, supposing that the actuators fail, the maximum output torque decreased, which causes degradation of the output of the actuator. Therefore, there is a need to plan the angular trajectory when the control torque is limited.
From the geometric relationship shown in Figure 2, after rotating −π/2 around x-axis and −φπ/2 around y-axis, the inertial frame is parallel to the orbital frame. The orbit frame is parallel to the main structure frame after it rotates γ around z-axis. The jth intermediate frame is parallel to the jth array frame after rotating θj around y-axis. The rotation matrix from the inertial frame to the jth array frame can be written as:
C j = C s j C a = C y ( θ j ) C z ( γ ) C y ( φ π / 2 ) C x ( π / 2 ) = cos φ sin θ j sin φ cos γ cos θ j sin φ sin θ j + cos φ cos γ cos θ j   sin γ cos θ j sin φ sin γ   cos φ sin γ   cos γ cos φ cos θ j sin φ cos γ sin θ j sin φ cos θ j + cos φ cos γ sin θ j sin γ sin θ j
One rotation can be expressed by a unit quaternion Q = [ q 0 q T ] T , where q 0 R and q R 3 are the scalar and vector part of the unit quaternion respectively. The rotation of the main structure frame relative to the inertial frame can be expressed as:
Q a = [ q a 0 q a T ] T = 2 / 2 cos ( φ / 2 + γ / 2 + π / 4 ) 2 / 2 cos ( φ / 2 γ / 2 + π / 4 ) 2 / 2 cos ( φ / 2 γ / 2 π / 4 ) 2 / 2 cos ( φ / 2 + γ / 2 π / 4 )
The first and second derivatives can be expressed as:
Q ˙ a = 2 / 4 sin ( φ / 2 + γ / 2 + π / 4 ) ( φ ˙ + γ ˙ ) 2 / 4 sin ( φ / 2 γ / 2 + π / 4 ) ( φ ˙ γ ˙ ) 2 / 4 sin ( φ / 2 γ / 2 π / 4 ) ( φ ˙ γ ˙ ) 2 / 4 sin ( φ / 2 + γ / 2 π / 4 ) ( φ ˙ + γ ˙ )
Q ¨ a = 2 / 4 sin ( φ / 2 + γ / 2 + π / 4 ) ( φ ¨ + γ ¨ ) 2 / 8 cos ( φ / 2 + γ / 2 + π / 4 ) ( φ ˙ + γ ˙ ) 2 2 / 4 sin ( φ / 2 γ / 2 + π / 4 ) ( φ ¨ γ ¨ ) + 2 / 8 cos ( φ / 2 γ / 2 + π / 4 ) ( φ ˙ γ ˙ ) 2 2 / 4 sin ( φ / 2 γ / 2 π / 4 ) ( φ ¨ γ ¨ ) + 2 / 8 cos ( φ / 2 γ / 2 π / 4 ) ( φ ˙ γ ˙ ) 2 2 / 4 sin ( φ / 2 + γ / 2 π / 4 ) ( φ ¨ + γ ¨ ) 2 / 8 cos ( φ / 2 + γ / 2 π / 4 ) ( φ ˙ + γ ˙ ) 2
The attitude kinematics equation of the main structure can be expressed as follows:
Q ˙ a = 1 2 G a ω a
where G a = [ q a ( q a 0 E 3 + q a × ) T ] T , and it is easy to find that G a T G a = E 3 where E3 is the identity matrix. So ωa and ω ˙ a can be expressed as:
ω a = 2 G a T Q ¨ a
ω ˙ a = 2 G ˙ a T Q ˙ a + 2 G a T Q ¨ a
Substituting Equations (28)–(34), (36) and (37) into Equations (16) and (17), the compact matrix notation form of Tc and τscj can be written as [ T c ] 3 m × 1 and [ τ s c j ] m × 1 . Compared with Tc, the value of τscj is much smaller, so only the value of Tc is constrained. Considering that the maximum torque of Tc is T c max = [ T c 1 max T c 2 max T c 3 max ] T and the minimum torque of Tc is T c min = [ T c 1 min T c 2 min T c 3 min ] T , 6 m inequality constraints are obtained as follows:
1 m × 1 T c min [ T c ] 3 m × 1 1 m × 1 T c max
where 1 m × 1 is a matrix with all the elements equal to 1. denotes the Kronecker product operator.
In theory, if the integrated attitude energy control system [18] is employed for the MR-SPS, the process of generating the control torque does not consume any energy. Continuous and efficient energy supply is the main task, so the energy receiving efficiency is selected as the performance index. The cosine of the angle between the normal vector of the array and the vector of the sun can measure the energy receiving efficiency. The following equation can calculate the average energy receiving efficiency in a solar day for one array:
η j = 0 T z ^ s j T C j r ^ s d t = 0 1 z ^ s j T C j r ^ s d τ
By calculation at the discretization points, the continuous angular trajectory optimization problem for the MR-SPS can be converted into the following small-scale nonlinear programming problem:
min J = k = 1 2 η k [ X γ ] ( n γ 3 ) × 1 [ X θ j ] ( n θ 3 ) × 1 ( j = 1 , 2 ) s . t . 1 m × 1 T c min [ T c ] 3 m × 1 1 m × 1 T c max
[ X γ ] ( n γ 3 ) × 1 and [ X θ j ] ( n θ 3 ) × 1   ( j = 1 , 2 ) are the unknown Bezier coefficients, a total of n r + 2 n θ 9 variables need to be optimized.

3.3. Initial Value Selection

The selection of the initial value of the variables has a great impact on the result of the optimization. The initial value close to the optimal solution can avoid local optimization and improve the operation speed. Section 2.2 has given the ideal angular trajectory, namely the analytical angular trajectory without dynamic constraints which can be used as a reference for the selection of initial value. At most time, by calculating ideal γ and θj at discrete time, the compact matrix notation forms [ γ i n i t i a l ] m × 1 and [ θ i n i t i a l j ] m × 1 are obtained. The initial values [ X γ i n i t i a l ] ( n γ 3 ) × 1 and [ X θ j i n i t i a l ] ( n θ 3 ) × 1   ( j = 1 , 2 ) can be calculated by solving the following equation to get the least square solution:
[ γ i n i t i a l ] m × 1 = [ B γ ] m × ( n γ + 1 ) [ P γ , 0 P γ , 1 [ X γ i n i t i a l ] ( n γ 3 ) × 1 T P γ , n γ 1 P γ , n γ ] T
[ θ i n i t i a l j ] m × 1 = [ B θ ] m × ( n θ + 1 ) [ P θ j , 0 P θ j , 1 [ X θ j i n i t i a l ] ( n θ 3 ) × 1 T P θ j , n θ 1 P θ j , n θ ] T
However, around equinox, the rate of change of the ideal angular trajectory is large when the right ascension between the sun and the SPS is about 0° or 180°, so the inversely calculated initial value violates the dynamic constraints a lot. Nevertheless, at this moment, the altitude angle of the sun is small, and quite high energy receiving efficiency is guaranteed as long as the ya-axis is perpendicular to the orbital plane. [ γ i n i t i a l ] m × 1 and [ θ i n i t i a l j ] m × 1 are given by: γ = 0 and θ j = π / 2 + 2 π τ . By solving Equations (41) and (42), the initial value [ X γ i n i t i a l ] ( n γ 3 ) × 1 and [ X θ j i n i t i a l ] ( n θ 3 ) × 1   ( j = 1 , 2 ) can be acquired.
By solving the optimization using two sets of initial value, the better result is selected as the optimal solution. The following criterion can help choose the initial value:
σ = 0 1 γ ¨ 2 d τ
where σ represents the smoothness of the curve.

4. Simulation Results

To verify the effectiveness of the proposed Bezier curve-based shaping approach, the trajectory planning for the MR-SPS in different time of the year is conducted. The order of the Bezier curve function is selected as n γ = n θ = 20 and the number of LG discrete points is selected as m = 50 . The nonlinear programming problem is solved using sequential quadratic programming. The parameters of the MR-SPS are listed in Table 1.

4.1. Angular Trajectory around Solstice

Assume that the MR-SPS is on orbit at 18 June 2025 and the maximum torque and minimum torque are set as T c max = [ 80 80 80 ] T   N m and T c min = [ 80 80 80 ] T   N m . The optimal angular trajectory of γ and θj generated by the proposed Bezier approach and the ideal angular trajectory of γ and θj are shown in Figure 3. It can be obtained that the optimal angular trajectory generated by the Bezier approach is very close to one ideal trajectory because the rate of change of the rotation angle is always small and the control effort is enough to track the ideal trajectory. The average energy receiving efficiency is 99.97% which is very close to 100%.
The control torque history is shown in Figure 4. It can be observed that all the components of Tc satisfy the constraints, and the peak value of Tc of the ideal trajectory is 150 N·m and is larger than that of the optimal trajectory which is 80 N·m. Although the peak value of control torque is less for the optimal trajectory, the loss of the efficiency is negligible. The value of the component Tcy is much less than the other components. The histories of τsc1 and τsc2 are almost the same due to the symmetrical structure and the history of τscj of the optimal trajectory and ideal trajectory is similar. If the constraints are less strict, the ideal trajectory can also be applied.

4.2. Angular Trajectory around Equinox

Assume that the MR-SPS is on orbit at 15 September 2025 and the maximum control torque and minimum control torque are set as T c max = [ 100 100 100 ] T   N m and T c min = [ 100 100 100 ] T   N m . The optimal angular trajectory of γ and θj generated by the proposed Bezier approach and the ideal angular trajectory are shown in Figure 5. It can be obtained that the optimal angular trajectory generated by the Bezier approach is very close to one ideal trajectory at the beginning and the end but in the middle it switches to another ideal trajectory, because the rate of change of the rotation angle is quite large at about 21,600 s and 64,800 s and the control effort is not enough to track the ideal trajectory. The average energy receiving efficiency is 99.97% which is very close to 100%.
The control torque history is shown in Figure 6. It can be observed that all the components of Tc satisfy the constraints of the optimal trajectory, and the peak value of Tc of the ideal trajectory is over 104 N·m and is far larger than that of the optimal trajectory which is 100 N·m. The value of the component Tcy is much less than the other components. The history of τsc1 and τsc2 is almost the same due to the symmetrical structure but the history of τscj of the optimal trajectory and ideal trajectory is different during the period when the trajectory is switched.

5. Conclusions

The generation of angular trajectory for a MW-level MR-SPS using Bezier shaping approach is demonstrated in this paper. The precise solar and earth orientation for the MR-SPS is analyzed and the ideal angular trajectory without dynamic constraints is given which is helpful for the initial value selection of the following optimization. The attitude dynamics and kinematics are modeled for calculating the inequality constraints of the optimization. The angular trajectory of the MR-SPS is assumed to be described by Bezier curve functions. The angular trajectory design around equinox and solstice is performed. The numerical simulation verifies that the Bezier shaping approach can generate an angular trajectory that meet the dynamic constraints. The energy receiving efficiency can be near 100% and the demand for torque actuators is reduced.

Author Contributions

Conceptualization, S.X.; methodology, S.X. and M.H.; software, S.X. and M.H.; validation, S.X. and W.F.; formal analysis, S.X. and M.H.; investigation, S.X. and W.F.; resources, S.X.; data curation, T.L.; writing—original draft preparation, S.X.; writing—review and editing, Z.L. and T.L.; visualization Z.L. and S.X.; supervision, N.Q. and M.H.; project administration, N.Q.; funding acquisition, N.Q. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by the National Science Foundation of China, grant number 11702072; in part by the China Postdoctoral Science Foundation, grant number 2017M611372; in part by the Heilongjiang Postdoctoral Fund, grant number LBH-Z16082; and in part by the National Natural Science Foundation of China, grant number 61903245.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Glaser, P.E. Power from the sun: Its future. Science 1968, 162, 857–861. [Google Scholar] [CrossRef] [PubMed]
  2. Koomanoff, F.A. Satellite power system concept development and evaluation program. Doe Nasa Rep. 1981, 1–6. [Google Scholar]
  3. Carrington, C.; Fikes, J.; Gerry, M.; Perkinson, D.; Feingold, H.; Olds, J. The Abacus/Reflector and integrated symmetrical concentrator—Concepts for space solar power collection and transmission. In Proceedings of the 35th Intersociety Energy Conversion Engineering Conference and Exhibit, Las Vegas, NV, USA, 24–28 July 2000. [Google Scholar]
  4. Mankins, J.; Kaya, N.; Vasile, M. SPS-ALPHA: The First Practical Solar Power Satellite via Arbitrarily Large Phased Array (A 2011–2012 NIAC Project). In Proceedings of the 10th International Energy Conversion Engineering Conference, Atlanta, GA, USA, 30 July–1 August 2012. [Google Scholar]
  5. Seboldt, W.; Klimke, M.; Leipold, M.; Hanowski, N. European Sail Tower SPS concept. Acta Astronaut. 2001, 48, 785–792. [Google Scholar] [CrossRef]
  6. Sasaki, S.; Tanaka, K.; Higuchi, K.; Okuizumi, N.; Kawasaki, S.; Shinohara, N.; Senda, K.; Ishimura, K. A new concept of solar power satellite: Tethered-SPS. Acta Astronaut. 2007, 60, 153–165. [Google Scholar] [CrossRef]
  7. Hou, X.B.; Wang, L.; Zhang, X.H.; Zhou, L. Concept design on multi-rotary joints SPS. J. Astron. 2015, 36, 1332–1338. [Google Scholar]
  8. Yang, Y.; Zhang, Y.; Duan, B.; Wang, D.; Li, X. A novel design project for space solar power station (SSPS-OMEGA). Acta Astronaut. 2016, 121, 51–58. [Google Scholar] [CrossRef]
  9. Li, Q.; Wei, Y.; Wu, Z.; Jiang, J. Novel orbit-attitude combination mode for solar power satellites to reduce mass and fuel. Chin. J. Aeronaut. 2022, 35, 132–142. [Google Scholar] [CrossRef]
  10. Wie, B.; Roithmayr, C.M. Integrated Orbit, Attitude, and Structural Control System Design for Space Solar Power Satellites. In Proceedings of the AIAA Guidance, Navigation, and Control Conference and Exhibit, Montreal, QC, Canada, 6–9 August 2001. [Google Scholar]
  11. Wie, B.; Roithmayr, C.M. Attitude and Orbit Control of a Very Large Geostationary Solar Power Satellite. J. Guid. Contr. Dyn. 2005, 28, 439–451. [Google Scholar] [CrossRef]
  12. McNally, I.J.; Scheeres, D.J.; Radice, G. Attitude Dynamics of Large Geosynchronous Solar Power Satellites. In Proceedings of the AIAA/AAS Astrodynamics Specialist Conference, San Diego, CA, USA, 4–7 August 2014. [Google Scholar]
  13. McNally, I.; Scheeres, D.; Radice, G. Locating Large Solar Power Satellites in the Geosynchronous Laplace Plane. J. Guid. Contr. Dyn. 2015, 38, 489–505. [Google Scholar] [CrossRef]
  14. Zhang, K.M.; Wu, S.N.; Wu, Z.G. Multibody Dynamics and Robust Attitude Control of a MW-level Solar Power Satellite. Aerosp. Sci. Technol. 2021, 111, 106575. [Google Scholar] [CrossRef]
  15. Li, Q.; Deng, Z.; Zhang, K.; Wang, B. Precise Attitude Control of Multirotary-Joint Solar-Power Satellite. J. Guid. Control. Dynam. 2018, 41, 1435–1442. [Google Scholar] [CrossRef]
  16. Li, Q.; Wei, Y.; Deng, Z.; Wu, Z.; Jiang, J. Switched iterative learning attitude and structural control for solar power satellites. Acta Astronaut. 2021, 182, 100–109. [Google Scholar] [CrossRef]
  17. Crassidis, J.L.; Junkins, J.L. Optimal Estimation of Dynamic Systems; Chapman and Hall/CRC: Boca Raton, FL, USA, 2004. [Google Scholar]
  18. Roes, J.B. An Electro-Mechanical Energy Storage System for Space Application. In Energy Conversion for Space Power; Snyder, N.W., Ed.; Academic Press: Cambridge, MA, USA, 1961; pp. 613–622. [Google Scholar]
Figure 1. Structure of the MW-level MR-SPS.
Figure 1. Structure of the MW-level MR-SPS.
Aerospace 09 00528 g001
Figure 2. The definitions of the five reference frames.
Figure 2. The definitions of the five reference frames.
Aerospace 09 00528 g002
Figure 3. Angular trajectory around solstice. (a) Angular trajectory of γ; (b) angular trajectory of θj.
Figure 3. Angular trajectory around solstice. (a) Angular trajectory of γ; (b) angular trajectory of θj.
Aerospace 09 00528 g003
Figure 4. Control torque history around solstice. (a) Control torque history of the components of Tc; (b) Control torque history of τscj.
Figure 4. Control torque history around solstice. (a) Control torque history of the components of Tc; (b) Control torque history of τscj.
Aerospace 09 00528 g004
Figure 5. Angular trajectory around equinox. (a) Angular trajectory of γ; (b) angular trajectory of θj.
Figure 5. Angular trajectory around equinox. (a) Angular trajectory of γ; (b) angular trajectory of θj.
Aerospace 09 00528 g005
Figure 6. Control torque history around equinox. (a) Control torque history of the components of Tc; (b) control torque history of τscj.
Figure 6. Control torque history around equinox. (a) Control torque history of the components of Tc; (b) control torque history of τscj.
Aerospace 09 00528 g006
Table 1. Parameters of the MR-SPS.
Table 1. Parameters of the MR-SPS.
ParametersValue
Ia diag   ( 4.4 × 10 9 ,   1 × 10 8 ,   4.4 × 10 9 )   kg m 2
Isj diag   ( 1.88 × 10 8 ,   2 . 08 × 10 7 ,   2.08 × 10 8 )   kg m 2
ma 8 × 10 4   kg
msj 2.5 × 10 4   kg
rc1 [ 0 230 30 ] T   m
rc2 [ 0 230 30 ] T   m
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xu, S.; Huo, M.; Qi, N.; Feng, W.; Lin, T.; Li, Z. Angular Trajectory Design for MR-SPS Using Bezier Shaping Approach. Aerospace 2022, 9, 528. https://doi.org/10.3390/aerospace9100528

AMA Style

Xu S, Huo M, Qi N, Feng W, Lin T, Li Z. Angular Trajectory Design for MR-SPS Using Bezier Shaping Approach. Aerospace. 2022; 9(10):528. https://doi.org/10.3390/aerospace9100528

Chicago/Turabian Style

Xu, Song, Mingying Huo, Naiming Qi, Wenyu Feng, Tong Lin, and Zheng Li. 2022. "Angular Trajectory Design for MR-SPS Using Bezier Shaping Approach" Aerospace 9, no. 10: 528. https://doi.org/10.3390/aerospace9100528

APA Style

Xu, S., Huo, M., Qi, N., Feng, W., Lin, T., & Li, Z. (2022). Angular Trajectory Design for MR-SPS Using Bezier Shaping Approach. Aerospace, 9(10), 528. https://doi.org/10.3390/aerospace9100528

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop