Research Article  Open Access
Qingtao Wang, Dongping Jin, Xiaoting Rui, "Dynamic Simulation of Space Debris Cloud Capture Using the Tethered Net", Space: Science & Technology, vol. 2021, Article ID 9810375, 11 pages, 2021. https://doi.org/10.34133/2021/9810375
Dynamic Simulation of Space Debris Cloud Capture Using the Tethered Net
Abstract
Space debris, especially the space debris cloud, has threatened severely the safety of future space missions. In the framework of multibody system dynamics, a computational approach is proposed in this study to investigate the dynamics of net deployment and capture of space debris cloud using this net subject to large overall motions and large deformations. To obtain high simulation fidelity of capturing space debris cloud, the gradient deficient beam element of the absolute nodal coordinate formulation (ANCF) is employed to discretize threads which are woven into the net. The normal contact force between the net and the debris cloud and among debris particles is computed by using the penalty method. Some deployment examples are presented to investigate the influences of shooting velocity of bullets and microgravity as well as the angle between the traveling direction of the net and the microgravity direction on the deployment characteristics of the tethered net. Other capturing examples are given to clarify the effect of the deployment area of the net at the moment it starts to contact with the debris cloud on the capture rate and to demonstrate the effectiveness of the proposed approach for capturing space debris cloud using the net in microgravity conditions.
1. Introduction
Space debris in orbit around the Earth is produced by space activities and has been growing by nearly geometric progression [1]. These space debris may strike other space objects in use, significantly increase the number of small debris, and severely threaten the safety of future space missions [2–4]. Even more seriously, the chain reaction of collisions between space debris and other space objects possibly leads to a debris cloud composed of many tiny debris particles which may trigger the disastrous cascade effects of debris, named “Kessler syndrome” [5–7].
To mitigate this situation and keep safe space exploration, the space debris capturing and removal technique from overpopulated orbital regions has been one of the hot spots in space research for the entire international space community, and the past decades have witnessed significant theoretical and technical developments in this exciting research field [3, 8, 9]. Compared to the tentacle capturing method [10], single and multiple robotic arm capturing method [11, 12], tethergripper capturing method [13], and harpoon capturing method [14], the tethered net capturing method is regarded as one of the most promising capturing methods due to its features of flexibility, lightweight, efficient cost, and compatibility with various dimensions and shapes of space debris [8, 15].
Compared with the elastic continuum approach [16] and cubic Bspline approach [17], the massspring model [18] and the absolute nodal coordinate formulation (ANCF) [19] are the most commonly used modeling approaches of the net. Benvenuto et al. modeled the net via the massspring approach and studied different phases of the debris capture and disposal [20]. Botta et al. presented a scheme of contact detection between the net and the target based on the massspring approach and established a normal and frictional contact dynamic model [18]. Huang et al. and Zhang et al. proposed a massspring approach based on a maneuverable tethered space net used for noncooperative space target capture and removal, and the position and configuration of this net are controllable via four maneuvering units located at four corners of a square net [21, 22]. The massspring model exhibits linearity of dynamic ordinary differential equations and computational efficiency in net modeling. However, the approach is not an ideal modeling method for the net due to its low accuracy in net dynamic deployment and the lack of ability to describe the nonlinear and large deformation of the tether.
The finite elements of ANCF introduce the global position vector and position vector gradients as nodal coordinates, leading to a constant mass matrix and zero centrifugal and Coriolis forces [23]. Numerous important applications have been identified for the potential use of ANCF on the deformation of fibrous soft tissue [24], underwater vehicle [25], fluidstructure interaction [26], and topology optimization [27], just to name a few. Recently, ANCF has been employed to establish the dynamic model of a space net [28]. Shan et al. compared the ANCF model and massspring model of the space net by using a parabolic flight experiment and concluded that the ANCF model is superior to the massspring model in describing the dynamics of the inner knots of the net [29, 30]. Qingquan et al. investigated the dynamic modeling and ground test of the tethered net and showed that the ANCF model is more capable of describing the flexibility of the net with fewer nodes than the massspring model [31].
In addition to the dynamic modeling of the space net, researchers have focused their attention on the capture process of space debris, including net deployment [30], contact between the net and the target [32], and closing of the net [33]. Botta et al. established the models of the chaser, net, target, contact dynamics, and closing mechanism used for tethered net capture of space debris and simulated the process of capturing a small asteroid [34]. Si et al. established the selfcontact dynamic model of the tethered net based on its geometric characteristics [35]. They also proposed a new split closing mechanism model by using the threadring sliding joint and validated the effectiveness of the model by comparing the simulation results and the experimental ones [36]. Zhang and Huang improved the net capture performance by controlling the thrusters installed on the bullets [37]. To test the performance of active debris removal technologies on mock targets in low Earth orbit, Forshaw et al. and Aglietti et al. have developed techniques and devices, for example, navigation, hardware, and software, to capture noncooperative targets in their RemoveDEBRIS demonstration mission [38–40]. Furthermore, they captured a CubeSat in orbit successfully in 2018 using a net launched from another CubeSat.
The aforementioned literature focuses attention on the debris capture with only one target. Moreover, the target is assumed to be infinitely massive, and as a result, forces and moments transmitted to the target upon the impact of the net have a negligible effect on its motion [18]. Little attention has been paid to the capture of the space debris cloud, which could be produced when space debris impacts other debris or space targets, for example, the Whipple shield of spacecraft [41], as shown in Figure 1. The masses of these debris particles in the cloud are very small, and their dynamic responses upon the impact of the net cannot be neglected. The description of the space debris cloud is generally based on the finite element method and smoothed particle hydrodynamics by simplifying the tiny debris particles in the cloud into spheres with masses and radii [7]. Therefore, this paper is aimed at proposing a computational approach in the framework of multibody system dynamics to simulate the net deployment and capture of space debris cloud and to investigate the response of the debris cloud in the capture process.
The remaining part of the paper is organized as follows. In Section 2, the ANCF gradient deficient beam element is described to establish the dynamic model of the threads which are woven into the net. In Section 3, the multizone contact detection strategy and the penalty method are introduced to establish the normal contact model between the net and the debris cloud. Section 4 presents some numerical examples to investigate the influences of shooting velocity of bullets, microgravity, and the deployment area of the net at the contact moment on the deployment characteristics of the net and on the capture rate. In Section 5, some concluding remarks are made.
2. System Modeling
As threads woven into the net are usually very thin, it is very suitable to establish the dynamic model of the net undergoing large deformations and large overall motions in space based on the ANCF gradient deficient beam element [19]. Compared to the massspring model, the ANCF model exhibits a constant mass matrix of the dynamic system, geometric nonlinearity, zero centrifugal and Coriolis forces, and, consequently, good capability of describing the flexibility of the net [23].
Shan et al. applied the ANCF to model the net by introducing the local coordinate to obtain the position of an arbitrary point within only one element [30], which is not convenient for the multizone contact detection strategy used in this paper to detect the contact between the net and the debris cloud [19] as in this strategy the detecting point may skip from one element to its adjacent one. Moreover, in their work [30], only the elastic forces due to the longitudinal deformation of threads are taken into account. The bending deformation which is very important for contact detection is not considered. So in this paper, the ANCF gradient deficient beam element suitable for the contact detection is introduced, both the longitudinal and bending deformations are taken into consideration, and the elastic force and Jacobians due to these deformations are derived in detail.
For the convenience of further discussion on contact detection, as shown in Figure 2, a new bodyfixed arc coordinate is introduced to describe the global position of an arbitrary point on a thread meshed by the ANCF gradient deficient beam. Thus, the global position vector of an arbitrary point on the centerline of the entire thread yields where is the arc coordinate of the thread and is its length in the undeformed configuration. and are the integer part and decimal part of the dimensionless value , denoting the element the point belongs to and the position on the element the point locates, respectively. and are the total element number used to discretize the thread and element length in the undeformed configuration, respectively.
is the global nodal coordinate vector of the element , composed of the positions and gradient coordinates, as where both and are global position vectors of two nodes (red dots in Figure 2) of the element; and are gradient coordinates of the two ends.
is the shape function matrix of the element, expressed as where , , , and ; is the identity matrix of order 3.
The constant element mass matrix of the thread can be written as [30] where both and are the density of the thread material and the area of the thread crosssection, respectively.
The gravitational force acting on the thread element, , is deduced by using the virtual work principle, as where is the geocentric gravitational constant.
The Jacobian of , which is needed in solving the system dynamic equations using the implicit integration algorithm, yields
The strain energy of the thread element is deduced as [23] where both and represent the longitudinal strain energy and the bending strain energy of the thread element, respectively. and are Young’s modulus of the thread material and the second moment of the area of the thread crosssection, respectively. Besides, and are the longitudinal strain and the curvature of an arbitrary point in the element. is the gradient vector, and .
The nonlinear elastic force associated with the longitudinal deformation of the element can be expressed as where and .
The Jacobian of the elastic force is derived as
Accordingly, the elastic force and the Jacobian of the element due to the bending deformation can be cast as where , , , , , and .The elastic force and its Jacobian of the deformed net can be deduced as where and are elastic forces due to the longitudinal deformation and bending deformation of the net, respectively. and are their Jacobians. They can be obtained by the assembly of elastic forces and Jacobians of all threads composing the net.
The space debris cloud is generally described based on the finite element method and smoothed particle hydrodynamics by simplifying the tiny debris particles in the cloud into spheres with masses and radii, each of which has 3 degrees of freedom [7]. So for an arbitrary debris particle with the global position , its mass matrix can be expressed as where is the mass of this particle.
The gravitational force acting on the debris particle, , can be written as
Its Jacobian with respect to the is deduced as
3. Contact Dynamic Modeling
In the process of capturing space debris cloud using the tethered net, contact occurs between the net and the debris cloud and the contact force will influence their dynamic responses. So accurate contact detection is critical to the accuracy of contact force calculation and, as a consequence, to the simulation fidelity in the capture phase. The massspring model used to establish the net [18, 32] can detect the contact between mass points of the net and the target but is incapable of detecting the contact between the thread segments and the target, leading to the low fidelity of capture simulation. This study investigates the capture of space debris cloud using the net discretized by the ANCF gradient deficient beam. This type of element can describe the tension, bending, and large deformation of an arbitrary point of the net. Therefore, it can detect contacts between the net and the debris cloud accurately. The tiny debris particles in the cloud are commonly simplified into spheres with masses and radii [7]. So only the contact detection between the net and the spheres should be operated in this study. Wang et al. developed a multizone contact detection strategy between two thin beams based on ANCF [19]. The degenerated form of this model is applied to the contact detection between threads and debris particles in this study. The reader interested in the details of this model may refer to the work by Wang et al. [19].
In this study, the contact force model between the net and the debris cloud neglecting friction is based on the penalty method, which has been widely used to study the contact problems because of its simplicity, no introduction of additional equations, and easy interpretation from a viewpoint of physics [42]. According to the penalty method, the contact force acted at an arbitrary point in the thread by debris particles is defined as where is the penalty parameter.
in Equation (17) is the normal penetration between the thread and the particle, defined as where and are the global positions of the contact point in a thread and debris particle, respectively; and are the radii of the thread and the particle.
denotes the normal contact direction at contact points and yields
According to the principle of virtual work, the generalized contact force acted at the thread point and debris particle can be written as where . Accordingly, the Jacobian of the generalized contact force can be deduced as where is the assembly of the element nodal coordinate and the global position of the debris particle.
The derivation process of the contact force between debris particles is omitted in this study as it is similar to that between threads and debris particles.
4. Deployment and Capture Simulation
In this section, simulations of deployment of a square net and the capture of space debris cloud using this net in microgravity conditions are, respectively, performed based on the ANCF model and the penalty method. The software MATLAB is used to implement these simulations. Each thread is discretized into 8 ANCF gradient deficient beam elements to obtain high simulation fidelity. The four bullets which are shot by the shooting mechanism to deploy the net are assumed to be lumped masses attached at the four corners of the net. The debris cloud is right above the net in the initial configuration. In the work by He et al. [7], the debris cloud shape is represented by the distribution of the particles in SPH. So in this study, the debris cloud is composed of four layers with uniformly distributed particles for each layer. The relative distance between the net and the bottom layer of the debris cloud above the net is set to be 5 m. The distances of adjacent layers and adjacent particles are set to be 0.45 m and 0.35 m, respectively. The net has an initial orbit altitude of 43164 km with an orbit angular velocity of 0.704 rad/s to maintain the circular orbit. The traveling direction of the net is parallel with the microgravity direction. The folding scheme of the net is the same as that by Shan et al. [30]. Some other input parameters for the simulations are shown in Table 1. In this table, the shooting angle is defined as the angle between the direction of the shooting velocity and the traveling direction of the net relative to the debris cloud [30].

4.1. Deployment Phase
Shan et al. studied in detail the influence of initial parameters, such as shooting velocity, shooting angle, and mass of bullets, on the deployment characteristics of the tethered net, namely, maximum area, deployment time, traveling distance, and effective period [30]. In this study, a new parameter to describe the dynamic characteristics of the net deployment phase, effective distance, is introduced. Similar to the definition of the effective period [30], the effective distance is defined as the distance in which the deployment area of a net is beyond 80% of its designed maximum area, as shown in Figure 3. To increase the probability of success of capturing debris cloud, it is expected that the debris cloud is within the effective distance which is as larger as possible.
The influence of the shooting velocity of bullets relative to the debris cloud on the deployment area of the net at the deployment phase is investigated in this study. The dimensionless deployment areas, calculated by dividing the deployment area by the designed maximum area (16 m^{2} in this study), versus both time and relative displacement with different shooting velocities, are shown in Figures 4 and 5, respectively. It can be seen from the two figures that with the increase of the shooting velocity, the maximum area is increasing, the deployment time is decreasing, and, accordingly, the capture is more efficient, which coincides with the conclusions in Shan et al. [30]. While the effective distance is increasing with the increase of the shooting velocity, as shown in Figure 5, the effective period for capturing is decreasing, as depicted in Figure 4. Shan et al. [30] pointed out that a shorter effective period leads the capturing to a higher risk of failing and is less reliable. Therefore, a reasonable shooting velocity should be adopted so as to achieve an optimal combination of efficient capture and high reliability. So in the next study, the shooting velocity is chosen as 21 m/s.
The influences of both the microgravity acted on the net by the Earth and the angle between the traveling direction of the net and the microgravity direction on the deployment area are also investigated in this study. The dimensionless deployment areas versus both time and relative displacement with and without microgravity are shown in Figures 6 and 7, respectively. In the two figures, the red curves depict the deployment area of the net considering microgravity, with whose direction the net traveling direction is parallel. The green curves, on the other hand, depict that with the microgravity direction perpendicular to the net traveling direction.
It can be concluded that compared to the results without consideration of microgravity, the maximum deployment area of the net with the traveling direction parallel with the microgravity direction is smaller, while that with the traveling direction perpendicular to the microgravity direction is larger. Both the deployment time and the traveling distance become longer and larger for the two cases. The effective period and effective distance for the first case almost do not change, while those for the second case become longer and larger. It can also be seen from the two figures that all the abovementioned differences are very slight.
The dynamic configurations of the net at six specific moments at the deployment phase are depicted in Figure 8. It can be seen from the figure that at the initial moment, the four bullets with a shooting velocity of 21 m/s drive threads at four corners of the net to deploy and then the rest of the threads of the net deploy successively. The relative velocities of bullets to the debris cloud decrease with time, and when their components in the plane perpendicular to the relative traveling direction decrease to zero (0.341 s, as shown in Figure 7), the deployment area of the net reaches its maximum value and the net starts to shrink due to the tensions of the threads. The figure also indicates that the ANCF model can describe the large flexibility and large deformation of the net.
4.2. Capture Phase
After the deployment of the net with a shooting velocity of 21 m/s and a traveling distance of 5 m, as shown in Figure 5, the net reaches its maximum area and starts to contact the debris cloud. The system dynamic configurations at eight specific moments at the capture phase are depicted in Figure 9. It can be seen from the figure that at 0.3 s, the net has not yet contacted the debris cloud, while at 0.35 s, the net contacted the first layer of the debris cloud. Later debris particles impact each other, causing the debris cloud a state of disorder, as shown in Figure 9 at the moment of 0.4 s. As the four bullets move with one velocity component towards the center of the net and the other one in the traveling direction, the deployment area of the net is decreased and the envelope surface made up of the net wraps some debris particles from the moment 0.5 s to the moment 0.7 s. In this process, more and more debris particles escape from the envelope surface due to the impact between particles and would not be captured successfully. With the bullets moving away from the center of the net after the moment the deployment area of the net reaches its maximum value, some other debris particles escape from the envelope surface successively while the rest of the particles are wrapped by the net more tightly and captured successfully, as shown in Figure 9 from the moment 0.8 s to the moment 0.9 s. From this viewpoint, it is difficult to capture the whole debris cloud due to the impact of particles. Increasing the designed maximum area of the net may be an effective way to capture more debris particles because the net can make a larger envelope space and can wrap more particles.
A debris particle that may be captured should be located within an envelope space made up of the net and a surface containing the four bullets. It should be noted that some of the debris particles within the envelope space can also escape from the space as this space is not closed, as shown in Figure 9. Moreover, the deployment area of the net at the moment it starts to contact with the debris cloud may affect the number of debris particles captured by the net. So the initial relative distance between the net and the bottom layer of the debris cloud is changed from 5 m to 4 m to investigate this influence. It can be seen from Figure 5 that the initial relative distance of 5 m is corresponding to the maximum deployment area of the net, about 92% of the designed maximum area, while the 4 m is corresponding to 80% of the designed one.
The comparison of the time histories of the number of debris particles that may be captured by the net with 92% (maximum deployment area) of the designed maximum area of the net and with 80% of that at the moment it starts to contact with the debris cloud is shown in Figure 10. It can be seen from the figure that after the moment the deployment area of the net is 92% (80%) of the designed maximum area, the number of debris particles that may be captured is increased sharply to 49 because all the particles of the first layer of the debris cloud are located in the interior of the envelope space. The rest of the layers of the debris cloud are captured successively, but the increasing rate of the number of captured particles becomes slow as some debris particles escape from the envelope space. At the moment of 0.47 s (0.36 s), the number reaches its maximum value and starts to decrease due to the impact among particles within the envelope space. From the moment about 0.95 s (1.05 s), the number of debris particles captured remains unchanged, about 72 (32), as particles are wrapped by the net tightly, as shown in Figure 9. These debris particles are the particles that are captured successfully by the net in the final configuration. Furthermore, it can be seen clearly that the net with 92% (maximum deployment area) of the designed maximum area at the contact moment can capture more debris particles than the net with 80% of that.
5. Concluding Remarks
A computational approach is proposed in this study to simulate the dynamics of net deployment and capture of space debris cloud using this net subject to large overall motions and large deformations. To obtain high simulation fidelity of capturing space debris cloud, the gradient deficient beam element of the absolute nodal coordinate formulation (ANCF) is employed to discretize threads which are woven into the net. The multizone contact detection strategy proposed by the authors in their previous work is used to detect contacts between threads and debris particles. The normal contact force between the net and the debris cloud and among debris particles is computed by using the penalty method.
A numerical example is made to investigate the influences of shooting velocity of bullets and microgravity on the deployment characteristics of the tethered net. The numerical results indicate that with the increase of the shooting velocity, the maximum area and the effective distance are increasing while the deployment time and the effective period are decreasing. Another example shows that the microgravity and the angle between the traveling direction of the net and the microgravity direction have slight influences on the net dynamic characteristics in the deployment phase. The numerical results also validate the good performance of the proposed approach on the description of large flexibility and large deformation of the net.
The last numerical example is presented to verify the effectiveness and applicability of the proposed approach on the dynamics of capturing the space debris cloud using the net in microgravity conditions. The results show that the debris cloud exhibits a state of disorder due to the impact of debris particles in the capture process and some of the debris particles are wrapped by the net while the others escape from the envelope surface. In the final configuration, some of the debris particles are wrapped tightly and captured successfully by the net. Furthermore, the net with the maximum deployment area of the designed maximum area at the contact moment can capture more debris particles than the net with 80% of that.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
Acknowledgments
This work was supported in part by the National Natural Science Foundation of China (Grant no. 11802130). It was also supported in part by the China Postdoctoral Science Foundation (Grant no. 2019M661849).
References
 Z. Huang, Y. Lu, H. Wen, and D. Jin, “Groundbased experiment of capturing space debris based on artificial potential field,” Acta Astronautica, vol. 152, pp. 235–241, 2018. View at: Publisher Site  Google Scholar
 J. C. Liou, N. L. Johnson, and N. M. Hill, “Controlling the growth of future LEO debris populations with active debris removal,” Acta Astronautica, vol. 66, no. 56, pp. 648–653, 2010. View at: Publisher Site  Google Scholar
 C. Bonnal, J. M. Ruault, and M. C. Desjean, “Active debris removal: recent progress and current trends,” Acta Astronautica, vol. 85, pp. 51–60, 2013. View at: Publisher Site  Google Scholar
 H. Wen, Z. H. Zhu, D. Jin, and H. Hu, “Constrained tension control of a tethered spacetug system with only length measurement,” Acta Astronautica, vol. 119, pp. 110–117, 2016. View at: Publisher Site  Google Scholar
 E. Levin, J. Pearson, and J. Carroll, “Wholesale debris removal from LEO,” Acta Astronautica, vol. 73, pp. 100–108, 2012. View at: Publisher Site  Google Scholar
 J. N. Pelton and W. H. Ailor, Space debris and other threats from outer space, Springer, 2013.
 Q. He, X. Chen, and J. F. Chen, “Finite elementsmoothed particle hydrodynamics adaptive method in simulating debris cloud,” Acta Astronautica, vol. 175, pp. 99–117, 2020. View at: Publisher Site  Google Scholar
 M. Shan, J. Guo, and E. Gill, “Review and comparison of active space debris capturing and removal methods,” Progress in Aerospace Sciences, vol. 80, pp. 18–32, 2016. View at: Publisher Site  Google Scholar
 C. P. Mark and S. Kamath, “Review of active space debris removal methods,” Space Policy, vol. 47, pp. 194–206, 2019. View at: Publisher Site  Google Scholar
 A. Chiesa, F. Fossati, G. Gambacciani, and E. Pensavalle, Enabling technologies for active space debris removal: the CADET project, Springer, 2015.
 T. Debus and S. Dougherty, “Overview and performance of the frontend robotics enabling nearterm demonstration (FREND) robotic arm,” in AIAA Infotech@Aerospace Conference, Reston, VA, USA, April 2009. View at: Publisher Site  Google Scholar
 K. Yoshida, “Achievements in space robotics,” IEEE Robotics & Automation Magazine, vol. 16, no. 4, pp. 20–28, 2009. View at: Publisher Site  Google Scholar
 P. Huang, J. Cai, Z. Meng, Z. Hu, and D. Wang, “Novel method of monocular realtime feature point tracking for tethered space robots,” Journal of Aerospace Engineering, vol. 27, no. 6, p. 04014039, 2014. View at: Publisher Site  Google Scholar
 C. Billot, S. Ferraris, R. Rembala, F. Cacciatore, A. Tomassini, and R. Biesbroek, “e.Deorbit: feasibility study for an active debris removal,” in 3rd European Workshop on Space Debris Modeling and Remediation, Paris, France, 2014. View at: Google Scholar
 P. Huang, F. Zhang, L. Chen et al., “A review of space tether in new applications,” Nonlinear Dynamics, vol. 94, no. 1, pp. 1–19, 2018. View at: Publisher Site  Google Scholar
 K. K. Mankala and S. K. Agrawal, “Dynamic modeling and simulation of satellite tethered systems,” Journal of Vibration and Acoustics, vol. 127, no. 2, pp. 144–156, 2005. View at: Publisher Site  Google Scholar
 S. Gao, Y. Yin, X. Sun, and Y. Sun, “Dynamic simulation of fishing net based on cubic Bspline surface,” Communications in Computer and Information Science, vol. 325, pp. 141–148, 2012. View at: Publisher Site  Google Scholar
 E. M. Botta, I. Sharf, and A. K. Misra, “Contact dynamics modeling and simulation of tether nets for spacedebris capture,” Journal of Guidance, Control, and Dynamics, vol. 40, no. 1, pp. 110–123, 2017. View at: Publisher Site  Google Scholar
 Q. Wang, Q. Tian, and H. Hu, “Dynamic simulation of frictional multizone contacts of thin beams,” Nonlinear Dynamics, vol. 83, no. 4, pp. 1919–1937, 2016. View at: Publisher Site  Google Scholar
 R. Benvenuto, M. Lavagna, and S. Salvi, “Multibody dynamics driving GNC and system design in tethered nets for active debris removal,” Advances in Space Research, vol. 58, no. 1, pp. 45–63, 2016. View at: Publisher Site  Google Scholar
 P. Huang, Z. Hu, and F. Zhang, “Dynamic modelling and coordinated controller designing for the manoeuvrable tethernet space robot system,” Multibody System Dynamics, vol. 36, no. 2, pp. 115–141, 2016. View at: Publisher Site  Google Scholar
 F. Zhang, P. Huang, Z. Meng, Y. Zhang, and Z. Liu, “Dynamics analysis and controller design for maneuverable tethered space net robot,” Journal of Guidance, Control, and Dynamics, vol. 40, no. 11, pp. 2828–2843, 2017. View at: Publisher Site  Google Scholar
 Q. Wang, Q. Tian, and H. Hu, “Dynamic simulation of frictional contacts of thin beams during large overall motions via absolute nodal coordinate formulation,” Nonlinear Dynamics, vol. 77, no. 4, pp. 1411–1425, 2014. View at: Publisher Site  Google Scholar
 L. P. Obrezkov, M. K. Matikainen, and A. B. Harish, “A finite element for soft tissue deformation based on the absolute nodal coordinate formulation,” Acta Mechanica, vol. 231, no. 4, pp. 1519–1538, 2020. View at: Publisher Site  Google Scholar
 T. Z. Htun, H. Suzuki, and D. GarcíaVallejo, “Dynamic modeling of a radially multilayered tether cable for a remotelyoperated underwater vehicle (ROV) based on the absolute nodal coordinate formulation (ANCF),” Mechanism & Machine Theory, vol. 153, p. 103961, 2020. View at: Publisher Site  Google Scholar
 F. Liu, G. Liu, and C. Shu, “Fluid–structure interaction simulation based on immersed boundarylattice Boltzmann flux solver and absolute nodal coordinate formula,” Physics of Fluids, vol. 32, no. 4, p. 047109, 2020. View at: Publisher Site  Google Scholar
 J. Sun, Q. Tian, and H. Hu, “Topology optimization based on level set for a flexible multibody system modeled via ANCF,” Structural and Multidisciplinary Optimization, vol. 55, no. 4, pp. 1159–1177, 2017. View at: Publisher Site  Google Scholar
 L. Liu, J. Shan, Y. Ren, and Z. Zhou, “Deployment dynamics of thrownet for active debris removal,” in 65th International Astronautical Congress, Toronto, Canada, 2014. View at: Google Scholar
 M. Shan, J. Guo, E. Gill, and W. Gołębiowski, “Validation of space net deployment modeling methods using parabolic flight experiment,” Journal of Guidance, Control, and Dynamics, vol. 40, no. 12, pp. 3319–3327, 2017. View at: Publisher Site  Google Scholar
 M. Shan, J. Guo, and E. Gill, “Deployment dynamics of tetherednet for space debris removal,” Acta Astronautica, vol. 132, pp. 293–302, 2017. View at: Publisher Site  Google Scholar
 C. Qingquan, Z. Qingbin, and T. Qiangang, “Dynamic modeling and ground test of tetherednet,” MATEC Web of Conferences, vol. 220, no. 4, p. 07003, 2018. View at: Publisher Site  Google Scholar
 M. Shan, J. Guo, and E. Gill, “Contact dynamic models of space debris capturing using a net,” Acta Astronautica, vol. 158, pp. 198–205, 2019. View at: Publisher Site  Google Scholar
 I. Sharf, B. Thomsen, E. M. Botta, and A. K. Misra, “Experiments and simulation of a net closing mechanism for tethernet capture of space debris,” Acta Astronautica, vol. 139, pp. 332–343, 2017. View at: Publisher Site  Google Scholar
 E. M. Botta, I. Sharf, and A. K. Misra, “Simulation of tethernets for capture of space debris and small asteroids,” Acta Astronautica, vol. 155, pp. 448–461, 2019. View at: Publisher Site  Google Scholar
 J. Si, Z. Pang, Z. du, and C. Cheng, “Dynamics modeling and simulation of selfcollision of tethernet for space debris removal,” Advances in Space Research, vol. 64, no. 9, pp. 1675–1687, 2019. View at: Publisher Site  Google Scholar
 J. Si, Z. Pang, Z. du, and J. Fu, “Dynamics modeling and simulation of a net closing mechanism for tethernet capture,” International Journal of Aerospace Engineering, vol. 2021, Article ID 8827141, 16 pages, 2021. View at: Publisher Site  Google Scholar
 F. Zhang and P. Huang, “Releasing dynamics and stability control of maneuverable tethered space net,” IEEE/ASME Transactions on Mechatronics, vol. 22, no. 2, pp. 983–993, 2016. View at: Publisher Site  Google Scholar
 J. L. Forshaw, G. S. Aglietti, N. Navarathinam et al., “RemoveDEBRIS: an inorbit active debris removal demonstration mission,” Acta Astronautica, vol. 127, pp. 448–463, 2016. View at: Publisher Site  Google Scholar
 J. L. Forshaw, G. S. Aglietti, S. Fellowes et al., “The active space debris removal mission RemoveDEBRIS. Part 1: from concept to launch,” Acta Astronautica, vol. 168, pp. 293–309, 2019. View at: Publisher Site  Google Scholar
 G. S. Aglietti, B. Taylor, S. Fellowes et al., “The active space debris removal mission RemoveDEBRIS. Part 2: in orbit operations,” Acta Astronautica, vol. 168, pp. 310–322, 2019. View at: Publisher Site  Google Scholar
 X. T. Zhang, G. H. Jia, and H. Huang, “A fast numerical approach for Whipple shield ballistic limit analysis,” Acta Astronautica, vol. 93, pp. 112–120, 2014. View at: Publisher Site  Google Scholar
 L. Przemysław, “Finite element analysis of beamtobeam contact,” Lecture Notes in Applied & Computational Mechanics, vol. 53, pp. 1–159, 2010. View at: Google Scholar
Copyright
Copyright © 2021 Qingtao Wang et al. Exclusive Licensee Beijing Institute of Technology Press. Distributed under a Creative Commons Attribution License (CC BY 4.0).