Next Article in Journal
Trajectory Optimization with Complex Obstacle Avoidance Constraints via Homotopy Network Sequential Convex Programming
Next Article in Special Issue
Prescribed Performance Adaptive Balance Control for Reaction Wheel-Based Inverted Pendulum-Type Cubli Rovers in Asteroid
Previous Article in Journal
Structural Efficiency Analysis of a Piston for Aviation Engines
Previous Article in Special Issue
Existence and Control of Special Orbits around Asteroid 4 Vesta
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability of a Flexible Asteroid Lander with Landing Control

School of Aerospace Engineering, Tsinghua University, Beijing 100084, China
*
Author to whom correspondence should be addressed.
Aerospace 2022, 9(11), 719; https://doi.org/10.3390/aerospace9110719
Submission received: 29 August 2022 / Revised: 31 October 2022 / Accepted: 10 November 2022 / Published: 16 November 2022
(This article belongs to the Special Issue Dynamics and Control Problems on Asteroid Explorations)

Abstract

:
Stable landing on asteroids is of considerable scientific and economic value but accompanied by huge difficulties. This paper proposes a novel flexible lander suitable for asteroids with microgravity and rugged surface. The gravity model with the artificial neural network and the surface model with the spherical harmonic method are introduced to establish the target asteroid’s dynamical environment. The flexible dynamics with the discrete shell model, the collision with the spring-damping model and viscous sliding friction, and the rigid coupling with the constraint violation stabilization method are elaborated for the lander. Combining the asteroid’s model with the lander’s dynamics, one successful landing scenario of the lander is presented. The lander’s landing stability of the final uncontrolled touching phase is studied through massive simulations. It is found that reasonable touching conditions can largely enhance the landing stability, and the lander can achieve a stable landing on the asteroid under a particular touching condition without control. The flexible lander’s comparison to the rigid lander is also discussed. It is concluded that the flexible lander does have higher adaptability and lower risk in asteroid landing. What is more, the attitude controller and position controller for the lander’s descent phase are also proposed and tested.

1. Introduction

Asteroid exploration is getting more and more attention due to its academic and economic value. It is widely recognized that some asteroids were born at the beginning of the solar system and therefore contained information about the solar system evolution [1,2], which is of great significance to planetary science. In addition, studies show that there exist rich rare metal resources on some asteroids [2]. With the development of deep space exploration, asteroids can be used as supply depots with abundant propulsion material.
With the development of aerospace technology and increasing interest in planetary protection, more and more asteroid exploration missions are being proposed and conducted. The exploration missions in the past can be roughly divided into circling exploration, short-range leap exploration, and landing exploration [3]. Landing exploration is a must to achieve comprehensive research about asteroids [4]. However, there exist a lot of difficulties in landing exploration. Different from planets with enough mass, asteroids’ gravity is so weak that the lander is likely to rebound or escape even at a slow velocity. Additionally, the asteroids’ terrain is usually rugged, which may cause the lander to overturn. The challenges above were first shown by the Philae lander, which tried but failed to stably land on Chury [5,6]. The Philae lander bounced for thousands of meters into a shaded crack after the failure of its anchoring harpoon system, and finally lost contact. In addition, the MINERVA-II and MASCOT landers of the Hayabusa2 mission also jumped for a long distance before their final stable landing [7,8]. Overall, the stable landing on asteroids is still tricky and worth further research.
Various feasible strategies for asteroid landing exploration have been proposed in past research. In the literature, there are several works that have focused on soft landing and in which asteroids’ orbit dynamics and the control scheme have been discussed [9,10,11,12]. There have also been some works about landing exploration with a CubeSat-based lander, as well as its simulation methods [13,14,15,16,17,18]. Some works have studied the improved anchoring method for asteroid landing based on Philae’s harpoon [19,20]. Recently, some researchers have noticed the unique advantages of flexible materials in asteroid landings. Some relevant works have focused on dissipating residual energy with flexible nets [21]; others have focused on the dynamical behavior of flexible materials in microgravity [4,22,23]. To expand on past research, this paper focuses on providing a feasible flexible lander of engineering significance for asteroid landing exploration, which should have both a flexible part to dissipate energy and a rigid part to carry out control and exploration. The asteroid flexible lander (shown in Figure 1) in this study is comprised of three rigid blocks connected with a thin layer of flexible materials, such as silicone rubber and vulcanized rubber. The necessary payloads, such as jets and the sampling instrument, are contained in the rigid blocks to conduct the control and sampling operations. When landing on the asteroid, the flexible lander will deform and dissipate energy via the deformation dissipation as well as the collision. Past research has indicated that the flexible structures can effectively avoid rebound and escape in weak gravity [4,21]. The symmetrical structure of the flexible lander can also avoid overturning on rugged terrain. Moreover, the symmetrical distribution and hexagonal prism shape of the blocks ensure that the blocks can be put together before the release phase to reduce space consumption, which is of great engineering significance. One asteroid landing mission can be divided into the launch of the exploration satellite, interplanetary transfer, orbital transfer to the target orbit, circling exploration, release of the lander, descent phase, and touching phase to the asteroid’s surface (see Figure 2). To avoid polluting the sampling area and damaging the lander during the landing, the jets are usually switched off after the descent phase, meaning the touching phase is uncontrolled. This paper mainly focuses on the touching phase, and the landing stability is discussed. In addition, a feasible controller for the descent phase is proposed and discussed briefly.
The main contributions of this paper can be summarized as follows. Firstly, this paper proposes an abstract physical model of a novel asteroid lander and develops the corresponding dynamical mode, including the asteroid’s dynamical model and the flexible landing dynamics, which are highly efficient and suitable for the flexible lander’s simulation. Secondly, this paper establishes the control scheme for the descent phase, in which reduced attitude equations are derived and the lander’s controllability is proved. Thirdly, this paper details the landing stability of the flexible lander for the touching phase, which gives accuracy requirements for the descent phase and shows the lander’s advantages in asteroid landing.
The rest of this paper is organized as follows. Section 2 establishes the asteroid’s dynamical environment, including the gravity and surface models. The following Section 3 details the flexible dynamics and the coupling dynamics of the flexible lander. Section 4 proposes a feasible controller for the descent phase for the flexible lander. Section 5 discusses one successful landing scenario, the landing stability of the proposed lander, the comparison to the rigid lander, and the simulation of the descent phase. Finally, the conclusion of this paper is given.

2. Asteroid’s Dynamical Environment

When exploring within the asteroid’s Hill radius for a short term, only the asteroid’s gravity and the collision need to be considered [24]. In this section, an efficient asteroid gravity model combing the artificial neural network (ANN) with the polyhedron method is proposed, and the asteroid surface model is also introduced.

2.1. Gravity Model

Asteroids usually have irregular shapes and mass distribution, which result in irregular gravity. To numerically model the asteroid’s irregular gravity, Werner et al. [25] established the polyhedron method, which has high accuracy but low efficiency. Meanwhile, ANN can fit and predict complex functions with high efficiency and accuracy [26]. This paper adopts the polyhedron method to calculate the sample data for ANN and then predicts the asteroid’s gravity.
The asteroid coordinate system is utilized, with the centroid as the origin, the minimum inertial axis as the x axis, the middle inertial axis as the y axis, and the maximum inertial axis as the z axis. In this coordinate system, the gravity is expressed as follows:
g ( r ) = G ρ f F S F f r f r ω f G ρ e E S E e r e L e
where g ( r ) is the asteroid’s gravity at position r , ρ is the asteroid’s average density, F S represents all polyhedral units of the polyhedral model, E S stands for all edges of F S , r f r is one vector from any point on f to the target point r , r e is one vector from any point on e to the target point, F f and E f are the tensors of f and e , respectively, and L e and ω f are two dimensionless coefficients of the polyhedral units, as shown in ref [21,25].
This paper focuses on the gravity near the asteroid’s surface, from which the sample data for the ANN are taken. Take Bennu as the target asteroid [27], whose size is 565   m   ×     535   m   ×     508   m . Accordingly, a ball with Bennu’s centroid as the center and a radius of 500     m was chosen as the sample space. The sample point and its gravity calculated with the polyhedron method were set as the ANN’s input and output, respectively. The number of sample data was 637,000, of which 90% were training data and 10% were validation data for the ANN; see Figure 3 for details.
Factoring in accuracy and efficiency, the authors chose a network with 4 hidden layers and 256 nodes in each layer. The activation function for the input and output layers was the hyperbolic tangent (tanh). The activation function for the hidden layers was the rectified linear unit (ReLU). In the training, the learning rate was 0.0001, the batch case was 200, the training epoch was 60, and the training strategy was stochastic gradient descent with momentum (SGDM).
f a c c u r a c y = 1 1 N i = 1 N N e t g ( r i ) g ( r i ) / g ( r i ) was adopted as the fitting accuracy to evaluate the ANN, where g ( r i ) is the gravity at r i , N means the number of sample data, and N e t g ( r i ) stands for the gravity calculated with ANN. After training, the fitting accuracy was 98.63%, which means a highly accurate prediction of the irregular gravity. Calculating 10,000 points’ gravity in the same computer, the ANN took 0.16 s while the polyhedron method took 360.43 s, proving the ANN’s great advantage in efficiency. What is more, the result of the ANN is an explicit continuous function, and therefore it is easy to be transplanted to other computing environments, and the gravitational gradient can also be estimated easily. Bennu’s gravity calculated by the ANN is shown in Figure 4.

2.2. Surface Model

The asteroid’s surface model is required when detecting the collision between the flexible lander and the asteroid. Usually, a polyhedron model is used to capture the asteroid’s irregular shape, and reconstruction is used to obtain a continuous surface expression [13]. The spherical harmonic method is an efficient and mature reconstruction method that transforms the irregular surface into a series of spherical harmonics [28]. Here, this paper briefly introduces the spherical harmonic method.
The asteroid coordinate system mentioned above is utilized in this section. In this coordinate system, the formula for the spherical harmonic method is:
R ( θ , φ ) = n = 0 m = n n C n m Y n m ( θ , φ )       θ [ 0 , π ] , φ [ 0 , 2 π ]
where R ( θ , φ ) , θ , and φ are the surface point’s three components in the spherical coordinates, and C n m stands for the spherical harmonic complex coefficient, which can be calculated by the least square method with the surface data. Y n m ( θ , φ ) represents the spherical harmonic basis function, and its formula is:
Y n m ( θ , φ ) = 2 n + 1 4 π · ( n m ) ! ( n + m ) ! P n m ( cos θ ) e i m φ
where P n m ( · ) is the Legendre polynomial; see Formula (4) for details:
P n m ( cos θ ) = ( 1 ) m 2 n n ! ( 1 cos 2 θ ) m / 2 d n + m d x n + m ( cos 2 θ 1 ) n
With Formulas (2)–(4), an explicit and continuous function for points on the asteroid’s surface was obtained. The formula for the surface’s local tangential and normal directions is:
{ τ 1 = R θ r ^ + R r ^ θ τ 2 = R φ r ^ + R r ^ φ λ = τ 1 × τ 2
where r ^ is the unit vector from the origin to the surface point, τ 1 and τ 2 are two tangent vectors, and λ stands for the normal vector of the surface point; see [28] for details. Figure 5 shows the surface of Bennu reconstructed with the spherical harmonic method.

3. Flexible Landing Dynamics

There are flexible and rigid parts in the flexible lander. A shell characteristic exists in the flexible part, resulting in low computational efficiency when simulating with the finite element method or absolute node coordinate method. At the final landing on the asteroid, the collision between the lander and the asteroid’s surface also needs to be taken into consideration. In this section, the discrete shell model suitable for the simulation of the flexible shell under large deformation is adopted. The tangential and normal collisions between the lander and the asteroid’s surface are modeled, and the coupling dynamics for coupling the flexible and rigid parts are also introduced.

3.1. Dynamics of the Flexible Part

The discrete shell model divides the flexible part along the middle plane into triangular units. The mass is assigned to the triangular units’ vertex, called the mass point; see Figure 6 for details. The lander’s movement is equivalent to the movement of the mass points during the simulation. It is worth noting that in this paper, the surface of the whole lander, including the rigid parts, is divided into units. The mass points corresponding to the rigid parts are fixed relative to the centroid through constraints. To describe the shell’s deformation, the discrete shell model calculates the membrane energy and bending energy with the relative variation among units, and then obtains the inner forces with the two deformation energies [29,30,31,32].
When the flexible part deforms, one triangular unit’s membrane energy [31] is calculated by the following formula:
W m = E h 2 ( 1 μ 2 ) [ μ ( tr ε ) 2 + ( 1 μ ) tr ε 2 ] d A
where W m is the membrane energy for a deformed unit, E stands for Young’s modulus of the flexible material, h represents the thickness of the unit, μ means the Poisson’s ratio of the flexible material, tr ( · ) is the trace function, and ε stands for the strain. According to ref [31], tr ε and tr ε 2 can be rewritten as:
tr ε = ( l 1 2 L 1 2 ) cot α 1 + ( l 2 2 L 2 2 ) cot α 2 + ( l 3 2 L 3 2 ) cot α 3 2 A tr ε 2 = i j 2 Δ 2 l i Δ 2 l j i = 1 3 ( Δ 2 l i ) 2 64 A 2
where l 1 , l 2 , and l 3 are the three side lengths of the deformed unit, L 1 , L 2 , and L 3 represent the three side lengths of the undeformed unit, α 1 , α 2 , α 3 stand for the three angles of the undeformed unit, A is the area of the undeformed unit, and Δ 2 l i is equal to ( l i 2 L i 2 ) . More details about the symbols above are shown in Figure 7.
According to the Rayleigh-Ritz analysis, the inner force resulting from the membrane energy is the negative derivative of that membrane energy with respect to the positions of the unit’s mass points. The formula is:
F m i = j i [ K i j ( l k 2 L k 2 ) + K i ( l j 2 L j 2 ) + K j ( l i 2 L i 2 ) ] ( x j x i )
where F m i is the membrane inner force at point Q i   ( i = 1 , 2 , 3 ) , ( i , j , k ) represents the non-repeated arrangement of ( 1 , 2 , 3 ) , and x i stands for the position vector of Q i . K i j and K i are constants related to the undeformed unit:
K i j = E h ( 2 cot α k + 1 μ ) 16 ( 1 μ 2 ) A K i = E h ( 2 cot α j cot α k + μ 1 ) 16 ( 1 μ 2 ) A
With Formulas (6)–(9), the membrane inner force caused by one deformed unit is derived. The total membrane inner force is obtained by summing the forces resulting from all deformed units. It is worth noting that only the units’ side lengths and position vectors are required to update during the simulation, which ensures high computational efficiency.
When the flexible part deforms, the bending energy between two connected units [29] is calculated by the following formula:
W b = E h 3 24 ( 1 μ 2 ) e ¯ 0 h ¯ e ( θ ¯ e θ e ) 2
where W b is the bending energy for the two connected units, θ e is the angle between the two connected units with deformation, θ ¯ e is the angle between the two connected units without deformation, and e ¯ 0 is the length of the undeformed common edge. h ¯ e is equal to ( h ¯ f 0 + h ¯ f 1 ) / 6 , where h ¯ f 0 and h ¯ f 1 represent the undeformed heights of the two units on the common edge. respectively. See Figure 8 for more details about the symbols above.
Like the membrane energy, the inner force resulting from the bending energy is the negative derivative of that bending energy with respect to the position of the unit’s mass points. Noting that only θ e varies during the simulation, the formula is written as follows:
F b i = E h 3 12 ( 1 μ 2 ) e ¯ 0 h ¯ e ( θ ¯ e θ e ) θ e x i i = 0 , 1 , 2 , 3
where F b i stands for the bending force on mass point x i . According to ref [30], θ e / x i is written as:
θ e x 0 = cot α 03 n 0 + cot α 04 n 1 e 0 θ e x 1 = cot α 01 n 0 + cot α 02 n 1 e 0 θ e x 2 = e 0 2 A f 0 n 0 θ e x 3 = e 0 2 A f 1 n 1
where α 01 represents the angle between x 0 x 1 and x 0 x 2 , α 02 stands for the angle between x 0 x 1 and x 0 x 3 , α 03 is the angle between x 1 x 0 and x 1 x 2 , α 04 represents the angle between x 1 x 0 and x 1 x 3 , A f 0 stands for the area of unit f 0 , and A f 1 means the area of unit f 1 . The other symbols are the same as mentioned above.
Shells dissipate energy via bending oscillations [32]. Here. the dissipation force is given directly as a supplement to the bending force:
F c i = E h 3 γ 12 ( 1 μ 2 ) e ¯ 0 h ¯ e θ ˙ e θ e x i
where F c i is the dissipation force on the mass point x i , and γ is the dissipation coefficient, which is set as 5   ×   10 4 in this paper.
By combining Formulas (10)–(13), the force resulting from bending is obtained. The total bending force is then obtained by calculating all two connected units. In each simulation step, only the θ e / x i and θ e are recalculated, and therefore, high efficiency is ensured.

3.2. Collision Model

The landing exploration must consider the collision between the flexible model and the asteroid’s surface. A precise collision model will contribute a lot to the choices of the flexible material, landing site, and so on. In this study, the collision was equivalent to the normal and tangential interactions between the mass point and the surface. The spring-damping model based on the Hunt-Crossley theory and the Hertz model was adopted for the normal collision [4,33,34,35]. Such a collision model has been widely used in past research, so only the main equations are presented here. The formula for the normal collision is as follows:
F n = ( K n d + C n d ˙ ) n · δ
where F n is the normal collision force of the mass point, and d stands for the embedding depth of the mass point into the surface and is equal to R ( θ , φ ) r , where r is the mass point’s position vector. R ( θ , φ ) means the asteroid surface distance to the origin with the same zenith angle and azimuth angle to the mass point; see Section 2.2 for details. d ˙ is the component of the mass point’s velocity in the surface normal direction, n is the unit normal vector of the asteroid surface, which is introduced in Section 2.2. δ is the collision detection function; if d < 0 , δ = 1 , otherwise δ = 0 . K n and C n are the stiffness coefficient and damping coefficient from the Hertz model and Hunt-Crossley theory, respectively:
K n = 4 3 κ E C n = 2 K n m π 2 + ( ln e ) 2 ln e
where m is the mass point’s mass, e means the coefficient of restitution, κ stands for the equivalent curvature radius between the mass point and the asteroid surface, and E represents the equivalent Young’s modulus. The formulas for κ and E are:
κ = ( 1 κ 1 + 1 κ 2 ) 0.5 E = ( ( 1 μ 1 2 ) E 1 + ( 1 μ 2 ) E ) 1
where κ 1 is the asteroid’s curvature radius at the collision point, κ 2 means the mass point’s radius, which can be equivalent to a constant [21]. E 1 and μ 1 are the estimated Young’s modulus and Poisson’s ratio of material on the asteroid, and E and μ represent Young’s modulus and Poisson’s ratio of the flexible material, respectively.
The viscous sliding friction was adopted for the tangential collision of the mass point and fully considers the asteroid surface’s characteristics. The formula is as follows:
F t = { ( ζ F n · v τ / v τ ) δ , K t s < μ F n ( K t s + C t v τ ) δ , K t s μ F n
where F t stands for the force resulting from the tangential collision, ζ means the coulomb friction coefficient, v τ is the mass point’s velocity with respect to the collision point’s tangential direction. δ represents the collision detection function, as mentioned above. s stands for the mass point’s slippage on the surface, which is equal to t 0 t v τ d t , where t 0 represents the start time of the collision. K t is the tangential friction coefficient, C t means the damping coefficient of the tangential frictions, which are considered to be equal to K n /   10 , and C n /   100 , respectively.

3.3. Coupling Dynamics

This paper presents the completed asteroid gravity model, the flexible dynamics for the flexible part, and the collision model. The coupling between the flexible and rigid parts was then solved with the constraint violation stabilization method [36]. In the asteroid coordinate system, the dynamical equation for the mass point i without connection to the rigid parts is:
m i r ¨ i = m i [ 2 ω × r ˙ i + ω × ( ω × r i ) ] + m i g ( r i ) + ( F n ( i ) + F t ( i ) ) + F i i n
where m i is the mass of the mass point i , r i represents the position vector of the mass point, ω stands for the asteroid’s angular velocity, g ( r i ) means the asteroid gravity, F n ( i ) and F t ( i ) are the normal and tangential collision forces on the mass point respectively, and F i i n represents the inner force from the flexible part. The dynamical equations for the mass point j connected to the rigid parts are:
m j r ¨ j = m j [ 2 ω × r ˙ j + ω × ( ω × r j ) ] + m j g ( r j ) + ( F n ( j ) + F t ( j ) ) + F j i n ϕ j = R c ( r j r c ) + C = 0
where ϕ j is the constraint on the mass point j , λ j is the Lagrange multiplier, r c is the rigid parts’ centroids, R c represents the transformation matrix from the asteroid’s ontology coordinate system to the rigid parts’ ontology coordinate system, and C is a constant decided by the lander’s undeformed state. The dynamical equations for the rigid part k are:
m k r ¨ k = m k [ 2 ω × r ˙ k + ω × ( ω × r k ) ] + m k g ( r k ) Λ ¨ k = ( R T J R ) 1 ( 1 2 R T M r 2 R T R R ˙ T J R Λ ˙ k ) ϕ k = Λ k T Λ k 1 = 0
where Λ k is the quaternion of the rigid part, J stands for the inertia matrix, M r is the moment on the rigid part, ϕ k is the constraint for Λ k , and R represents a 3 × 4 matrix consisting of Λ k ; see ref [37] for details.
The generalized coordinates q   =   [ r 1 r N p r c 1 r c 2 r c 3 Λ 1 Λ 2 Λ 3 ] T were adopted to rewrite Formulas (18)–(20), where N p is the number of the mass point, r 1 r N p stands for the position vectors of all mass points, and [ r c 1 r c 2 r c 3 Λ 1 Λ 2 Λ 3 ] T represent the position and attitude vectors of the three rigid parts. The rewritten formula is:
M q ¨ + Φ q T λ = F Φ = 0
where M is the generalized mass matrix, F is the generalized force, Φ is the vector consisting of all constraints, and Φ q T is the transposed Jacobian matrix of Φ . The constraint violation stabilization method solves Equation (21) as follows (see ref [38] for more details):
[ M Φ q T Φ q 0 ] [ q ¨ λ ] = [ F ( Φ q q ˙ ) q q ˙ 100 Φ ˙ q 2500 Φ q ]
By solving the linear Equation (22), the generalized acceleration q ¨ is solved. Combined with the initial condition, the lander’s movement over time is obtained via numerical integration.

4. Control Scheme for the Descent Phase

In this section, the control scheme for the descent phase is proposed, which guides the flexible lander from the mother spacecraft to an appropriate state for the final touching phase. The corresponding simulation and test are set in Section 6.

4.1. Attitude Control

The reduced attitude equations governing the flexible lander’s attitude and its control scheme are derived in this part. The complete attitude formula for the flexible lander is:
I ˜ ( t ) ω ˜ ˙ + ω ˜ × ( I ˜ ( t ) ω ˜ ) = M ˜
where I ˜ ( t ) is the inertia matrix, ω ˜ is the angular velocity of the lander, and M ˜ is the total moment from the control thruster and the environment. It is worth noting that I ˜ ( t ) is not a constant matrix in the lander’s ontology coordinate system since the lander will deform during the motion. To reduce the difficulty in designing the control scheme, the deformation was not considered and the simplified attitude formula is as follows:
I ω ˙ + ω × ( I ω ) = M
where I represents the inertia matrix, ω is the angular velocity of the lander, and M stands for the moment caused by the control thruster and the environment under the premise that the flexible lander will not deform, respectively. In the lander’s ontology coordinate system, I and Formula (24) can be written as:
I = [ I x x 0 0 0 I y y 0 0 0 I z z ]
I x x ω ˙ x = ( I z z I x x ) ω y ω z + M c x + M e x I y y ω ˙ y = ( I z z I y y ) ω x ω z + M c y + M e y I z z ω ˙ z = M c z + M e z
where I x x =     134.58   kg / m 2 , I y y =     134.58   kg / m 2 , and I x x =     73.81   kg / m 2 in this research, ω x , ω y , and ω z are the three elements of ω , M c = [ M c x , M c y , M c z ] T stands for the moment from the control thruster, and M e = [ M e x , M e y , M e z ] T means the moment from the environment. I x x , I y y , and I z z are the elements of I , and I x x is equal to I y y , so ω z is the decoupling, as shown in Formula (26). As mentioned above, M c z is equal to 0 if there is no deformation, and only the ω x , ω y are controllable. According to ref [39], this study adopted the reduced attitude equations for the flexible lander, which means controlling a single unit vector of the attitude instead of controlling the whole attitude.
The equation governing the unit vector n = [ n x , n y , n z ] T , which is expressed in the lander’s ontology coordinate system but fixed in the inertial system, is as follows:
n ˙ = ω × n
where n = 1 . Combining Formula (26) with (27), and disregarding the environmental disturbance, the reduced attitude equations are proposed:
I x x ω ˙ x = ( I z z I x x ) ω y ω z + M c x + M e x I y y ω ˙ y = ( I z z I y y ) ω x ω z + M c y + M e y n ˙ x = ω y n z + ω z n y n ˙ y = ω x n z ω z n x
x = [ ω x , ω y , n x , n y ] T [ ω ¯ x , ω ¯ y , n ¯ x , n ¯ y ] T was chosen as the state vector, where [ ω ¯ x , ω ¯ y , n ¯ x , n ¯ y ] T represents the targeted state. u = [ M x , M y ] was chosen as the input. Under the condition of one order of simplification, Formula (28) is derived into a linear time-invariant system:
x ˙ = A x + B u A = [ 0 I x x I z z I x x ω ¯ z 0 0 I z z I y y I y y ω ¯ z 0 0 0 0 n ¯ z 0 ω ¯ z n ¯ z 0 ω ¯ z 0 ] , B = [ 0 1 I x x 1 I y y 0 0 0 0 0 ]
By checking the controllability matrix’s rank, it is easy to find that when n ¯ z 0 , the system is controllable. Considering the time invariance of the attitude equation and that n ¯ z is designable, it can be concluded that the flexible lander’s reduced attitude is controllable. With Equation (29), the controller can be easily designed by using the linear-quadratic regulator (LQR) method. LQR is a traditional and mature method for the linear time-invariant system, and therefore, this paper does not discuss it deeply [40]. The descent simulation with the LQR controller is discussed in the following section. With the attitude controller, n is controllable to the target, which means that the lander’s attitude is controllable.

4.2. Position Control

The position controller was designed under the premise that the flexible lander will not deform, and the lander’s normal direction is controllable since the thrust of the proposed lander is along the lander’s normal direction. A feasible control is that the response time of the attitude controller is much shorter than the position controller to ensure an appropriate thrust direction.
This study adopted the PID controller for the position control, which makes the lander’s ideal acceleration d ¨ d e s as a second-order system, as shown in Formula (30).
d ¨ d e s + 2 η ω n d ˙ + ω n 2 d = 0
where d stands for the vector from the lander to the target point, η is the damping ratio, and ω n represents the natural frequency of the PID controller. The total thrust f is derived as follows:
f = m n d e s · d ¨ d e s
where m is the total mass of the lander, and n d e s stands for the instantaneous normal direction of the lander.
Combining the input u = [ M x , M y ] of the attitude controller and the total thrust f of the position controller, the thrust f 1 , f 2 , f 3 that the three thrusters should produce during the descent phase is as follows:
[ M x M y f ] = [ 0 l / 2 l / 2 l / 3 l / 2 3 l / 2 3 1 1 1 ] [ f 1 f 2 f 3 ]
where l means the distance between each two rigid parts’ centroid when the lander is not deformed.

5. Landing Simulation and Discussions

This paper presents the completed method for simulating asteroids’ dynamic environments and the motion of the flexible lander, which contributes to the foundation for the landing simulation of the asteroid. In this section, a successful landing scenario on the asteroid, the landing stability of the flexible lander, and the comparison between the flexible lander and the rigid lander are presented and discussed. The simulation of the descent phase is also shown in this section.

5.1. Stable Landing Scenario

This section presents one successful landing of the flexible lander, demonstrating its feasibility. The asteroid Bennu was chosen as the landing target, which was also the target of NASA’s OSIRIS-Rex. The material with a Young’s modulus of 0.001     GPa , a Poisson’s ratio of 0.4, and a density of 1000     kg / m 3 was chosen as the flexible part’s material, which corresponds to the silicone rubber. The thickness of the flexible part was 0.005     m and the total mass of the flexible material was 53.19     kg . The mass and the moment of inertia for three rigid parts were all set as 150     kg and diag ( [ 10 10 10 ] )     kg / m 3 , respectively. The initial center of the flexible lander was located at [ 13 13 255 ]   T   m , the initial attitude quaternion between the asteroid’s ontology coordinate system and the lander’s ontology coordinate system was [ 0 . 9582 0 0 . 2860 0 ] T , and the initial speed was [ 0 . 0199 0 - 0 . 0980 ] T     m / s . Under the above conditions, the first collision happened at [ 13.3 13.0 253.6 ] T   m , which is a relatively flat area on the Bennu, and therefore was chosen as an ideal landing point in this study. The whole simulation time was 2000     s .
It is difficult to use the lander’s centroid to describe the movement due to its deformation during the collision, so this study adopted the centroid of the three rigid parts instead. Figure 9a,b shows the trajectory of the flexible lander on the asteroid. Figure 9c,d shows the flexible lander at 17 s and 2000 s, respectively. Figure 10 shows the variation in the distance between the center of the lander and the targeted landing site.
Figure 9 and Figure 10 show that with the conditions above, the first collision happened at 17 s, during which the flexible lander deformed obviously, as shown in Figure 9c, and its energy dissipated through the collision and deformation. The asteroid’s gravity was so weak that the lander could not achieve a stable landing through one collision. After the first collision, the lander moved in the form of a parabola under the microgravity and collided with the asteroid at 596 s again. The lander returned to the asteroid’s surface and then bounced several times after the second collision. The lander finally achieved a stable landing, as shown in Figure 9d. During a series of collisions, the lander did not leave the target point too far during the whole landing, and the final distance between the lander and the target point was just 0.8506   m , which shows the potential of the flexible lander in avoiding rebound and escape under microgravity. See Figure 9 and Figure 10 for more details.
All in all, the flexible lander completed a stable landing on the asteroid through passive deformation in the uncontrolled touching phase.

5.2. Landing Stability of the Flexible Lander

This section focuses on the flexible lander’s landing stability through massive simulations, guiding the touching condition of the lander. Too many factors affect the lander’s landing, such as the lander’s attitude, the uneven surface, etc. Taking so many factors into account is impossible, and therefore this study focused on three main elements, which are the flexible lander’s initial velocity v 0 , the angle β between the ground and the velocity, and the angle α between the ground and the lander. The three elements are also the main control variables in the practical touching of the asteroid lander, as shown in Figure 11. In this section, the rugged terrain is not considered as a main research factor since a flat area is normally chosen for landing in practice. In microgravity, the influence of rugged terrain and the error in observation of the asteroid’s surface can be equivalent to the errors of α and β . The irregular asteroid’s gravity is also assumed to be constant since the primary component of the gravity changes little on the asteroid’s surface, especially on the asteroid Bennu, which is close to spherical.
The rebound distance between the targeted landing site and the next touching point was considered to be the index for the landing stability. This study took 5     m as the criterion for a stable landing. In the simulation, v 0 was set as 0.1   0.3   m / s , α was set as 45 °     45 ° , β was set as 45 ° 45 ° , the parameters of the flexible lander were set as mentioned in Section 5.1, and the gravity was set as [ 0 0 8   ×   10 5 ]     m / s 2 , which is an estimate of the mean gravity on Bennu’s surface. Figure 12a–c illustrates the lander’s tangential velocity, normal velocity after the first collision, and its rebound distance in the condition of v 0 =   0.1   m / s , respectively. Local amplification of Figure 12c is shown in Figure 12d,e. Figure 12f illustrates the lander’s rebound distance in the condition of v 0 =   0.2   m / s . Figure 12g illustrates the lander’s rebound distance in the condition of v 0 = 0.3   m / s . The red contour in Figure 12c,f,g marks the feasible scope of α and β for a stable landing using the mentioned criterion. It is worth noting that the lander’s rebound distance in this study was solved by assuming the lander’s movement as a parabola with its tangential and normal velocity so as to reduce the simulation time.
Figure 12 shows that v 0 , α , and β significantly impacted the landing, proving the necessity of this research. Here, this paper presents a detailed discussion of Figure 12a–e, which corresponds to v 0 = 0.1   m / s . From Figure 12a, it can be concluded that the smaller α and β are, the larger the normal velocity will be, which brings more escaping risk to a landing. In Figure 12b, there is an opposite trend in the tangential velocity, so α and β will be neither too small nor too large. In Figure 12c, it is easy to find that a stable landing only happens in some areas of α and β . Reasonable α and β values can enhance the flexible lander’s adaptability on the asteroid through higher energy dissipation, while a wrong touching condition of the lander will cause a long-rebound distance. When landing under the conditions of v 0 = 0.1   m / s , α = 0 , and β = 0 , the allowable errors in α , β , and terrain slopes are ( 0.03 , 0.01 ) , ( 0.12 , 0.14 ) , and ( 0.03 , 0.01 ) , respectively, as illustrated in Figure 4d. When landing under the state of v 0 =   0.1   m / s , α = 0.58 , and β = 0.2 , corresponding to the initial condition in Section 5.1, the allowable errors in α , β , and terrain slopes are ( 0.12 , 0.11 ) , ( 0.23 , 0.17 ) , and ( 0.12 , 0.15 ) , respectively, which means much better adaptability on the asteroid, as seen in Figure 12e. α = 0 and β = 0 are usually chosen when landing on plants with larger gravity. However, this is not an appropriate choice under the microgravity of the flexible lander. Under the condition of v 0 = 0.2   m / s , the feasible scopes of α and β are much smaller than that of v 0 = 0.1   m / s , as shown in Figure 12f. When v 0 = 0.3   m / s , a stable landing is nearly impossible, as seen in Figure 12g. It is evident that the higher the v 0 is, the lower the possibility of a successful landing is, which is also reflected in these results.

5.3. Comparison with the Rigid Lander

This section discusses the flexible lander’s comparison with the rigid lander, showing the flexible lander’s potential advantage in asteroid exploration. The simulation was almost the same as that in Section 5.2, and only the material of the lander’s flexible part was substituted with steel with a Young’s modulus of 201   GPa and a Poisson’s ratio of 0.4. The rigid lander’s deformation was negligible with this stiffness. Figure 13a–c shows the rigid lander’s tangential velocity, the normal velocity after the first collision, and its rebound distance under the condition of v 0 = 0.1   m / s , respectively. Figure 13d,e shows the local amplification of the origin and the best α and β , respectively.
Compared to Figure 12, it was found that the rigid lander’s feasible scopes for α and β at v 0 = 0.1   m / s were much smaller than the flexible lander, meaning more risk in landing. The longest rebound distance was nearly double the flexible lander’s, indicating more divergence to the targeted landing site. The sensitivity to α and β are also remarkable for the rigid lander, showing the significance of this research. When landing under the conditions of v 0 = 0.1   m / s , α = 0 , and β = 0 , the allowable errors in α , β , and the terrain slopes are ( 0.02 , 0.02 ) , ( 0.38 , 0.38 ) , and ( 0.02 , 0.02 ) , respectively, as seen in Figure 13d. When landing under the conditions of v 0 = 0.1   m / s , α = 0.22 , and β = 0.66 , which are better touching conditions for the rigid lander, the allowable errors in α , β , and the terrain slopes are ( 0.23 , 0.08 ) , ( 0.07 , 0.06 ) , and ( 0.08 , 0.10 ) , respectively, meaning a higher possibility of a stable landing than landing under the conditions of v 0 = 0.1   m / s , α = 0 , and β = 0 , as seen in Figure 13e. Comparing the aforementioned results for the flexible and rigid lander, it is clear that the flexible lander had better adaptability on the uneven terrain and lower requirements for control.
In conclusion, the flexible lander does have an edge in asteroid landing.

5.4. Simulation of the Descent Phase

A feasible controller for the descent phase has been proposed above in theory. This section applies the controller to the flexible lander to prove its feasibility.
In the simulation, this study assumed that the three rigid parts’ position and velocity were measured in real-time, and the state vector needed by the controller was estimated accordingly. The weighted matrices Q and R of the LQR controller were both set as the unital matrix. The PID controller’s damping ratio and natural frequency were set to 2 and 0.01, respectively. In the beginning, the flexible lander was located at [ 0 0 280 ] T   m , the initial attitude quaternion between the asteroid’s ontology coordinate system and the lander’s ontology coordinate system was [ 1 0 0 0 ] T , and the initial speed was [ 0 0 0 ]   T   m / s . The target point was [ 2 2 288 ] T   m , and [ 0 . 2357 0 . 2357 0 . 9428 ] T   was chosen as target n . The parameters of the flexible lander were set as mentioned in Section 5.1. The simulation lasted for 1000 s, and the attitude controller began working at t =   0   s while the position controller began working at t =   100   s .
Figure 14 shows the distance between the lander and the target point, as well as the angle between the lander’s normal direction and the target n . The lander’s normal direction was calculated according to the three rigid parts’ position in case the lander deformed during the control. Figure 14 shows that the controller could control the flexible lander to the target point with centimeter position accuracy, and its normal direction to the target direction had nearly 0 attitude error above the asteroid.
In conclusion, the controller is feasible for the descent phase.

6. Conclusions

In this work, a novel flexible asteroid lander is proposed, and the corresponding dynamics, including the asteroid’s dynamical environment, the flexible dynamics, and the coupling dynamics, are established. Taking Bennu as the target asteroid, simulations of the final touching phase were performed and discussed, of which one successful stable landing scenario is shown, the landing stability is detailed, and the comparison to the rigid lander is presented. It is concluded that reasonable touching conditions can enhance the flexible lander’s adaptability, and the flexible lander has better stability as well as lower requirements for control than a rigid lander. Furthermore, the controller for the descent phase was designed, and a corresponding simulation was conducted. The results of the simulation show that the controller is feasible and the touching condition is controllable with the controller. To expand on this research, it will be valuable to consider an optimum design of the flexible lander, including the thickness of the flexible material and the distribution of the rigid blocks. Designing an experimental prototype and test beds to simulate the asteroid’s dynamical environment is also significant for future research. In addition, an optimal controller for the descent phase is worth further study.

Author Contributions

Investigation, methodology, software, writing—original draft, W.Y.; methodology, software, writing—review and editing, R.F.; methodology, project administration, writing—review and editing, H.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key R&D Program of China, grant number 2019YFA0706500.

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. Castillo-Rogez, J.C.; Pavone, M.; Nesnas, I.; Hoffman, J.A. Expected science return of spatially-extended in-situ exploration at small Solar system bodies. In Proceedings of the Aerospace Conference, Big Sky, MT, USA, 3–10 March 2012. [Google Scholar]
  2. Peloni, A.; Ceriotti, M.; Dachwald, B. Solar-Sail Trajectory Design for a Multiple Near-Earth-Asteroid Rendezvous Mission. J. Guid. Control Dyn. A Publ. Am. Inst. Aeronaut. Astronaut. Devoted Technol. Dyn. Control 2016, 39, 2712–2724. [Google Scholar] [CrossRef] [Green Version]
  3. Luo, Y.; Jiang, X.; Zhong, S.; Ji, Y.; Sun, G. Multisatellite Task Allocation and Orbit Planning for Asteroid Terminal Defence. Aerospace 2022, 9, 364. [Google Scholar] [CrossRef]
  4. Yan, W.; Baoyin, H. A Numerical Research on a New Type of Asteroid Flexible Probe. IEEE Access 2021, 9, 129863–129873. [Google Scholar] [CrossRef]
  5. Dorofeeva, V.A. Chemical and Isotope Composition of Comet 67P/ChuryumovGerasimenko: The RosettaPhilae Mission Results Reviewed in the Context of Cosmogony and Cosmochemistry. Sol. Syst. Res. 2020, 54, 96–120. [Google Scholar] [CrossRef]
  6. Hand, E. Philae probe makes bumpy touchdown on a comet. Science 2014, 346, 900–901. [Google Scholar] [CrossRef]
  7. Yoshikawa, M.; Kawaguchi, J.; Kuninaka, H. The Results of Asteroid Exploration Mission Hayabusa. Ieice Tech. Rep. 2010, 110, 219–222. [Google Scholar]
  8. Fujiwara, A.; Kawaguchi, J.; Yeomans, D.; Abe, M.; Mukai, T.; Okada, T.; Saito, J.; Yano, H.; Yoshikawa, M.; Scheeres, D. The rubble-pile asteroid Itokawa as observed by Hayabusa. Science 2006, 312, 1330–1334. [Google Scholar] [CrossRef]
  9. Li, X.; Qiao, D.; Huang, J.; Han, H.; Meng, L. Dynamics and control of proximity operations for asteroid exploration mission. Sci. China Phys. Mech. Astron. 2019, 49, 084508. (In Chinese) [Google Scholar] [CrossRef]
  10. Tang, G.; Jiang, F.; Li, J. Low-thrust trajectory optimization of asteroid sample return mission with multiple revolutions and moon gravity assists. Sci. China Phys. Mech. Astron. 2015, 58, 114501. [Google Scholar] [CrossRef]
  11. Chen, Y.; Baoyin, H.; Li, J. Design and optimization of a trajectory for Moon departure Near Earth Asteroid exploration. Sci. China Phys. Mech. Astron. 2011, 54, 748–755. [Google Scholar] [CrossRef]
  12. Huang, J.; Li, X.; Qiao, D.; Jia, F.; Meng, L. Mission design for multi-target and multi-mode rendezvous missions to small bodies. Sci. China Phys. Mech. Astron. 2019, 49, 084512. (In Chinese) [Google Scholar] [CrossRef]
  13. Wal, S.V.; Reid, R.G.; Scheeres, D.J. Simulation of Nonspherical Asteroid Landers: Contact Modeling and Shape Effects on Bouncing. J. Spacecr. Rocket. 2019, 57, 1–22. [Google Scholar]
  14. Bazzocchi, M.C.; Hakima, H. A CubeSat-based Robotic Asteroid Sampling Mission. In Proceedings of the AIAA Scitech 2020 Forum, Orlando, FL, USA, 6–10 January 2020. [Google Scholar]
  15. Kim, Y.-B.; Jeong, H.-J.; Park, S.-M.; Lim, J.H.; Lee, H.-H. Prediction and Validation of Landing Stability of a Lunar Lander by a Classification Map Based on Touchdown Landing Dynamics’ Simulation Considering Soft Ground. Aerospace 2021, 8, 380. [Google Scholar] [CrossRef]
  16. Zhang, Y.; Li, J.; Zeng, X.; Wen, T.; Li, Z. High-fidelity landing simulation of small body landers: Modeling and mass distribution effects on bouncing motion. Aerosp. Sci. Technol. 2021, 119, 107149. [Google Scholar] [CrossRef]
  17. Wen, T.; Zeng, X.; Circi, C.; Gao, Y. Hop reachable domain on irregularly shaped asteroids. J. Guid. Control Dyn. 2020, 43, 1269–1283. [Google Scholar] [CrossRef]
  18. Zen, X.; Li, Z.; Gan, Q.; Circi, C. Numerical study on the low-velocity impact of the asteroid lander onto deformable regolith surfaces. J. Guid. Control Dyn. 2022, 45, 1644–1660. [Google Scholar]
  19. Zhao, Z.; Zhao, J.; Hong, L. An asteroid landing mechanism and its landing simulation. In Proceedings of the 2012 IEEE International Conference on Robotics and Biomimetics (ROBIO), Guangzhou, China, 11–14 December 2012. [Google Scholar]
  20. Wang, Y.; Jiang, W.; Long, L.; Zhu, Q.; Feng, R.; Wang, L.; Antonio, C. Design and Experimental Research of a New Type of Asteroid Anchoring System. Int. J. Aerosp. Eng. 2021, 2021, 1–7. [Google Scholar]
  21. Zhang, Y.; Yu, Y.; Baoyin, H. Dynamical behavior of flexible net spacecraft for landing on asteroid. Astrodynamics 2021, 5, 249–261. [Google Scholar] [CrossRef]
  22. Bartlett, N.W.; Tolley, M.T.; Overvelde, J.T.B.; Weaver, J.C.; Mosadegh, B.; Bertoldi, K.; Whitesides, G.M.; Wood, R.J. A 3D-printed, functionally graded soft robot powered by combustion. Science 2015, 349, 161–165. [Google Scholar] [CrossRef] [Green Version]
  23. Feng, R.; Zhang, Y.; Liu, J.; Zhang, Y.; Li, J.; Baoyin, H. Soft Robotic Perspective and Concept for Planetary Small Body Exploration. Soft Robot. 2022, 9, 889–899. [Google Scholar] [CrossRef]
  24. Hamilton, D.P.; Burns, J.A. Orbital stability zones about asteroids. Icarus 1991, 92, 118–131. [Google Scholar] [CrossRef] [Green Version]
  25. Werner, R.A.; Scheeres, D.J. Exterior gravitation of a polyhedron derived and compared with harmonic and mascon gravitation representations of asteroid 4769 Castalia. Celest. Mech. Dyn. Astron. 1996, 65, 313–344. [Google Scholar] [CrossRef]
  26. Song, Y.; Cheng, L.; Gong, S. Fast estimation of gravitational field of irregular asteroids based on deep learning and its applications. In Proceedings of the 29th AAS/AIAA Space Flight Mechanics Meeting, Maui, HI, USA, 13–17 January 2019. [Google Scholar]
  27. Chesley, S.R.; Farnocchia, D.; Nolan, M.C.; Vokrouhlický, D.; Chodas, P.W.; Milani, A.; Spoto, F.; Rozitis, B.; Benner, L.A.; Bottke, W.F.; et al. Orbit and bulk density of the OSIRIS-REx target Asteroid (101955) Bennu. Icarus 2014, 235, 5–22. [Google Scholar] [CrossRef] [Green Version]
  28. Yu, Y.; Michel, P.; Hirabayashi, M.; Schwartz, S.R.; Zhang, Y.; Richardson, D.C.; Liu, X. The Dynamical Complexity of Surface Mass Shedding from a Top-shaped Asteroid Near the Critical Spin Limit. Astron. J. 2018, 156, 59. [Google Scholar] [CrossRef]
  29. Grinspun, E. A Discrete Model of Thin Shells; Birkhäuser: Basel, Switzerland, 2005. [Google Scholar]
  30. Wardetzky, M.; Bergou, M.; Harmon, D.; Zorin, D.; Grinspun, E. Discrete quadratic curvature energies. Comput. Aided Geom. Des. 2007, 24, 499–518. [Google Scholar] [CrossRef]
  31. Delingette, H. Triangular Springs for Modeling Nonlinear Membranes. IEEE Trans. Vis. Comput. Graph. 2008, 14, 329–341. [Google Scholar] [CrossRef]
  32. Mohamed, A.-N.A.; Shabana, A.A. A nonlinear visco-elastic constitutive model for large rotation finite element formulations. Multibody Syst. Dyn. 2011, 26, 57–79. [Google Scholar] [CrossRef]
  33. Shan, M.; Guo, J.; Gill, E. Contact dynamic models of space debris capturing using a net. Acta Astronaut. 2017, 158, 198–205. [Google Scholar] [CrossRef]
  34. Zhang, F.; Huang, P.; Meng, Z.; Zhang, Y.; Liu, Z. Dynamics Analysis and Controller Design for Maneuverable Tethered Space Net Robot. J. Guid. Control Dyn. 2017, 40, 1–16. [Google Scholar] [CrossRef]
  35. Si, J. Dynamics Modeling and Simulation of a Net Closing Mechanism for Tether-Net Capture. Int. J. Aerosp. Eng. 2021, 2021, 1–16. [Google Scholar] [CrossRef]
  36. Baumgarte, J.W. A New Method of Stabilization for Holonomic Constraints. J. Appl. Mech. 1983, 50, 869. [Google Scholar] [CrossRef]
  37. Liu, Y. Advanced Dynamics; Academic: Beijing, China, 2003; p. 113. [Google Scholar]
  38. Chiou, J.C.; Wu, S.D.; Wit, S.D. Constraint violation stabilization using input-output feedback linearization in multibody dynamic analysis. In Proceedings of the Guidance, Navigation, and Control Conference, San Diego, CA, USA, 29–31 July 1996. [Google Scholar]
  39. Mueller, M.W.; D’Andrea, R. Stability and control of a quadrocopter despite the complete loss of one, two, or three propellers. In Proceedings of the 2014 IEEE International Conference on Robotics and Automation (ICRA), Hong Kong, China, 31 May–7 June 2014. [Google Scholar]
  40. Tuna, S.E. LQR-Based Coupling Gain for Synchronization of Linear Systems. arXiv 2008, arXiv:0801.3390. [Google Scholar]
Figure 1. (a) The flexible asteroid landing lander; (b) the lander’s profile.
Figure 1. (a) The flexible asteroid landing lander; (b) the lander’s profile.
Aerospace 09 00719 g001
Figure 2. The concept of the asteroid exploration mission using the flexible lander.
Figure 2. The concept of the asteroid exploration mission using the flexible lander.
Aerospace 09 00719 g002
Figure 3. The distribution of sample points.
Figure 3. The distribution of sample points.
Aerospace 09 00719 g003
Figure 4. The learned gravity on Bennu’s surface.
Figure 4. The learned gravity on Bennu’s surface.
Aerospace 09 00719 g004
Figure 5. Bennu’s reconstructed asteroid surface.
Figure 5. Bennu’s reconstructed asteroid surface.
Aerospace 09 00719 g005
Figure 6. The units of the flexible lander.
Figure 6. The units of the flexible lander.
Aerospace 09 00719 g006
Figure 7. (a) The undeformed unit; (b) the deformed unit.
Figure 7. (a) The undeformed unit; (b) the deformed unit.
Aerospace 09 00719 g007
Figure 8. Illustration of the two connected units.
Figure 8. Illustration of the two connected units.
Aerospace 09 00719 g008
Figure 9. (a) The flexible lander on the asteroid; (b) the movement of the flexible lander on the asteroid; (c) the first collision at 17 s; (d) the flexible lander completing a stable landing.
Figure 9. (a) The flexible lander on the asteroid; (b) the movement of the flexible lander on the asteroid; (c) the first collision at 17 s; (d) the flexible lander completing a stable landing.
Aerospace 09 00719 g009
Figure 10. The distance between the flexible lander and the targeted landing site.
Figure 10. The distance between the flexible lander and the targeted landing site.
Aerospace 09 00719 g010
Figure 11. Illustration of the simulation conditions.
Figure 11. Illustration of the simulation conditions.
Aerospace 09 00719 g011
Figure 12. (a) The flexible lander’s tangential velocity after the first collision at v 0 = 0.1   m / s ; (b) the flexible lander’s normal velocity after the first collision at v 0 = 0.1   m / s ; (c) the flexible lander’s rebound distance at v 0 = 0.1   m / s ; (d) local amplification of the origin in Figure 4c; (e) local amplification of α = 0.58 , β = 0.2 in Figure 4c; (f) the flexible lander’s rebound distance at v 0 = 0.2   m / s ; (g) the flexible lander’s rebound distance at v 0 = 0.3   m / s .
Figure 12. (a) The flexible lander’s tangential velocity after the first collision at v 0 = 0.1   m / s ; (b) the flexible lander’s normal velocity after the first collision at v 0 = 0.1   m / s ; (c) the flexible lander’s rebound distance at v 0 = 0.1   m / s ; (d) local amplification of the origin in Figure 4c; (e) local amplification of α = 0.58 , β = 0.2 in Figure 4c; (f) the flexible lander’s rebound distance at v 0 = 0.2   m / s ; (g) the flexible lander’s rebound distance at v 0 = 0.3   m / s .
Aerospace 09 00719 g012
Figure 13. (a) The rigid lander’s tangential velocity after the first collision at v 0 = 0.1   m / s ; (b) the rigid lander’s normal velocity after the first collision at v 0 = 0.1   m / s ; (c) the rigid lander’s rebound distance at v 0 = 0.1   m / s ; (d) local amplification of the origin in (c); (e) local amplification of α = 0.22 , β = 0.66 in (c).
Figure 13. (a) The rigid lander’s tangential velocity after the first collision at v 0 = 0.1   m / s ; (b) the rigid lander’s normal velocity after the first collision at v 0 = 0.1   m / s ; (c) the rigid lander’s rebound distance at v 0 = 0.1   m / s ; (d) local amplification of the origin in (c); (e) local amplification of α = 0.22 , β = 0.66 in (c).
Aerospace 09 00719 g013
Figure 14. (a) The distance between the flexible lander and the target; (b) the angle between the lander’s normal direction and the target n .
Figure 14. (a) The distance between the flexible lander and the target; (b) the angle between the lander’s normal direction and the target n .
Aerospace 09 00719 g014
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yan, W.; Feng, R.; Baoyin, H. Stability of a Flexible Asteroid Lander with Landing Control. Aerospace 2022, 9, 719. https://doi.org/10.3390/aerospace9110719

AMA Style

Yan W, Feng R, Baoyin H. Stability of a Flexible Asteroid Lander with Landing Control. Aerospace. 2022; 9(11):719. https://doi.org/10.3390/aerospace9110719

Chicago/Turabian Style

Yan, Weifeng, Ruoyu Feng, and Hexi Baoyin. 2022. "Stability of a Flexible Asteroid Lander with Landing Control" Aerospace 9, no. 11: 719. https://doi.org/10.3390/aerospace9110719

APA Style

Yan, W., Feng, R., & Baoyin, H. (2022). Stability of a Flexible Asteroid Lander with Landing Control. Aerospace, 9(11), 719. https://doi.org/10.3390/aerospace9110719

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