# A Penetrating-Anchoring Mathematical Model for the Soft Asteroid Anchoring System

Authors: Zhijun Zhao; Shuang Wang; Jingdong Zhao
Journal: International Journal of Aerospace Engineering
 (2021)

---

## Abstract

The asteroid landing mechanism is necessary to be anchored to avoid flowing away. At present, the study on the anchoring system is mainly focused on the mechanical design, but there are few researches on the penetrating or anchoring mathematical model, and the researches on combining two models with each other are even more lacking. In the paper, based on the characteristics of Mohr-Coulomb material, a penetrating mathematical model of the anchoring system is established. This penetrating mathematical model can be used to calculate the penetrating depth of the anchor body according to the penetrating speed and the medium properties. Secondly, an anchoring mathematical model is established, which shows the relationships among the anchoring force, medium properties, and penetrating depth. Finally, a penetrating-anchoring mathematical model is built with the penetrating depth as the link. The model establishes a relationship between the anchoring force and the initial penetrating conditions.


---

## Body

## 1. Introduction

Asteroid exploration has high scientific and economic values [1–3], such as finding evolution history of the planetary system, finding the origin of the solar system and the life, and exploring the asteroid resources. Moreover, asteroid exploration opens up a new field of deep space exploring [4, 5].The landing mechanism must be fixed on the asteroid surface by the anchoring system to prevent flowing away due to the microgravity on the asteroid. The greater the anchoring force, the easier it is for the landing mechanism to be anchored on the asteroid surface. The anchoring force is closely related to the penetrating depth and the medium characteristics, and the penetrating depth is related to the mass of the anchor body, the shape of the anchor tip, the characteristics of the medium, and so on. Therefore, it is significant to analyze the penetrating-anchoring theory of the anchoring system on the asteroid surface. Penetrating analysis can predict penetrating depth of the anchoring system on the asteroid surface, and anchoring analysis can estimate the anchoring force that can be generated by the anchoring system at a certain penetration depth. Both of penetrating and anchoring analysis can guide the anchoring of the landing mechanism on the asteroid surface.At present, the anchoring theory research on the asteroid landing is still in the starting state, and the relevant theoretical model has not been reported. Currently, anchoring theory is mainly about the pile foundation engineering. Ilamparuthi and Dickin summarized some anchoring models established by the pioneers [6] (Majer et al. (1955) established a vertical slip surface model. Downs et al. (1966), Murray et al. (1987), and Clemence et al. (1977) established an inverted truncated cone model. Balla (1961) and Meyerhof et al. (1968) established a curved slip surface model.). Sutherland studied the anchoring force of pile-based anchors in clay and analyzed various models [7]. Dickin studied the anchoring of anchors in sand media from 1988 to 1992 [8]. In 1994, Rao and Kumar studied the anchoring capacity of pile-based anchors and analyzed various models [9]. Ilamparuthi and Dickin studied the anchorage effects of different ends in the sand medium in 2001 [6]. Merifield et al. studied the anchoring of pile anchors in clay [10]. Wu et al. studied the anchoring of pile-based anchors in plaster-cement mixed media [11]. Shanker et al. studied the anchoring of the light rod [12]. Niroumand et al. studied the anchoring of the late expanded pile foundation in the sand medium [13]. Liu et al. studied the deformation of the medium during the drawing process of the anchor [14]. Consoli et al. studied the anchoring force of the pile body in cement [15]. The anchorage of the pile foundation is larger (about several meters in length), and the ground buried depth is deep (tens of meters). It is distinctly different from the anchoring system in the volume and penetration depth. Pile anchorage is different from asteroid anchoring to some extent, but it can be taken as guidance and reference for the study of the asteroid anchoring.Currently, there are few researches on the anchoring model for the asteroid anchoring, let alone the penetrating-anchoring model. However, the mathematical model describing the relationship between the penetrating and the anchoring is necessary in asteroid anchoring. In the paper, the vertical penetration equation of the anchoring system on the asteroid surface is derived. This equation can describe the mathematical relationship among the penetration speed, penetration depth, and the medium properties. Thus, this equation can predict the penetration depth of the anchor body if the medium property is known. Then, the anchoring model of the anchoring system is established. The model can be used to estimate the anchoring force generated by the anchoring system. Finally, the penetrating-anchoring model is established, which displays the relationship between the penetration velocity, the medium characteristics, and the anchoring force.


## 2. Anchoring System

The landing mechanism anchors on the asteroid via the anchoring system. The anchoring system contains penetrating unit, advancing unit, winding unit, damping unit, and so on. Its schematics and performances are shown separately in Figure1. The mass of the system is about 1.5 kg. The penetrating unit drafting a wire of about 2 m in length can penetrate the soft asteroid surface. Then, a wire connection is established between the landing mechanism and the asteroid.Figure 1
Schematics of anchoring system.Penetrating unit in the anchoring system is active and deformable, and its mechanical structure is shown in the top graph of Figure2. The penetrating unit is mainly composed of anchor tip, anchor claw, barb, handspike, piston, pyrotechnics, anchor body, etc. As the penetrating unit penetrated the asteroid surface, the control system ignites the pyrotechnics and produces high-pressure gas. The gas pushes the handspike via the piston, inducing the splay of the claw and the barbs. The deployed penetrating unit is shown in the lower graph of Figure 2. In soft media, the claws and barbs are fully opened, and the open barbs increase the contact area between the anchor tip and the media, increasing the anchoring force. In hard media, the claws and barbs will generate irregular and plastic deformation when opening due to the great thrust generated by the pyrotechnic. This phenomenon is similar to the effect of expansion bolt, and it can enhance the anchoring force.Figure 2
Schematic of the penetrating unit.The anchoring force is a key index of the anchoring system, directly related to its penetrating depth of the penetrating unit. However, this penetrating depth depends on the initial conditions such as the medium properties, the shape of the anchor tip, the mass, and penetrating velocity. Thus, it is necessary to establish the relationship between the anchoring force and the initial condition, which can provide theoretical foundation for reliable anchorage of the anchoring system.


## 3. Penetrating Model

### 3.1. Penetrating Resistance Force Analysis

The anchoring force is directly related to the penetrating depth. In order to establish the relationship between the anchoring force and the initial condition, it is firstly necessary to establish the relationship between the penetrating depth and the initial condition.The forces on arbitrary shape anchor tip are shown in Figure3 in penetrating.Figure 3
Force on arbitrary anchor tip.Based on the cavity expansion model [16–18] and citing the deduction process of the paper [19], we can get the axial force on the anchor tip [18].
(1)Fz1=2π∫0rAτ0+Bρ0Vz211+φ′r2rdr+2πμ∫0zAτ0+Bρ0Vz2f′z21+f′z2fzdz0≤z≤pFz2=2π∫0RAτ0+Bρ0Vz211+φ′r2rdr+2πμ∫0pAτ0+Bρ0Vz2f′z21+f′z2fzdzz>p,where
(2)A=1α1+τ0/2Eγ2α−1λ,B=3/1−η∗1−2α2−α+1/γ21+τ0/2E/γ2α3τ0/E+η∗1−3τ0/2E2−γ321−η∗2−α+3γ3/1−η∗1−2α2−α1+τ0/2E4.α=3λ/3+2λ, where λ=tanφ and φ is the internal friction angle.η∗=1−ρ0/ρ∗, where ρ0 is the volume density before distortion and ρ∗ is the volume density after distortion.
(3)γ=Vc=1+τ02E−1−η∗1/3,where E is the Young’s modulus and τ0 is the cohesion of the medium.


### 3.2. Penetrating Depth

The generatrix equation of the conic anchor tip is as follows:(4)rz=Rp=tanθ,where θ is the semivertical angle of the anchor tip and θ=9.46.Obtaining(5)r=fz=ztanθ,z=φr=rcotθ,f′z=tanθ,φ′r=cotθ.Integrating the above equation with equation (4) will obtain
(6)Fz1=πr2Aτ0+Bρ0Vz2⋅sin2θ+2πμ∫0zAτ0+Bρ0Vz2⋅sin2θztanθdz0≤z≤pFz2=πR2Aτ0+Bρ0Vz2⋅sin2θ+2πμ∫0pAτ0+Bρ0Vz2⋅sin2θztanθdzz>p.(1)
When penetrating depth0≤z≤pTo the conic anchor tip, according to equation (6),
(7)Fz1=πr2Aτ0+Bρ0Vz2⋅sin2θ1+μcotθ0≤z≤p.It can be written as follows:(8)Fz1=πr2Aτ01+μcotθ+πr2Bρ0sin2θ1+μcotθVz2.Substitutingr=ztanθ into the equation above, it is obtained that
(9)Fz1=α1z2+β1z2Vz2,where
(10)α1=πtan2θAτ01+μcotθ,β1=πtan2θBρ0sin2θ1+μcotθ.According to Newton second law, we can get(11)mVzdVzdz=−α1+β1Vz2z2,m2dVz2α1+β1Vz2=−z2dz.Integrating(12)−m2β1lnα1+β1Vz2+C1=13z3.Whent=0, Vz=V0, and z=0, we find
(13)C1=m2β1lnα1+β1V02.SubstitutingC1 into equation (12), we can get
(14)z=3m2β1lnα1+β1V02α1+β1Vz21/3.Whenz=p,
(15)Vzp=α1β1+V02e−2β1p3/3m−α1β1.When maximum penetrating depthh1<p, setting Vz=0, we can get
(16)h1=3m2β1ln1+β1α1V021/3.By the way, according to the Newton second law, we can get(17)md2zdt2=−α1+β1Vz2z2.AsVz=dz/dt, the above equation can be written as
(18)md2zdt2=−α1+β1dzdt2z2.Integrating the above equation, we can get the relation between the depthz and the time t. And then differentiating, we can get the relation between the acceleration and the time. In the paper, we are only interested in the relationship between the initial penetration velocity and the final penetration depth. Thus, the solution to the equation is not discussed in detail.
(2)
When the penetrating depthz>pWhenz>p, we can get the antiforce of the anchor tip by equation (6).
(19)Fz2=πR2Aτ01+μcotθ+πR2Bρ0sin2θ1+μcotθVz2z>p.It can be written as follows:(20)Fz2=α2+β2Vz2where
(21)α2=πR2Aτ01+μcotθ,β2=πR2Bρ0sin2θ1+μcotθ.According to the Newton second law, we can get(22)mVzdVzdz=−α2+β2Vz2.Obtaining(23)m2dVz2α2+β2Vz2=−dz.Integrating(24)−m2β2lnα2+β2Vz2+C2=z.Whenz=p, Vz=Vzp, getting
(25)C2=m2β2lnα2+β2Vzp2+p.SubstitutingC2 into equation (26), we can get
(26)z=m2β2lnα2+β2Vzp2α2+β2Vz2+p.SettingVz=0, we can get the maximum penetrating depth h2h2>p.
(27)h2=m2β2ln1+β2α2Vzp2+p.By the way, according to the Newton second law, we can get(28)mdVzdt=−α2+β2Vz2.Integrating the above equation, we can get(29)Vz=α2β2tan−α2β2mt+C3.Taking the initial conditionVztp=Vzp (tp is the time of z=p) into the above equation, we can obtain the relationship between Vz and the time t. Then, the relationship between the acceleration and the time can be known.


### 3.3. Penetrating Test

In order to verify the abovementioned penetration equation, the relationship between the penetrating depth and the penetrating speed of the anchor body in hard clay and fine sand media is tested experimentally. Testing platform of the anchoring system is shown in Figure4. The anchoring element is fired by the powder, which is ignited via the trigger, and the trigger is controlled by the electric cable. When the anchor body is fired by the anchoring element, the firing speed of the anchor body is measured by the speedometer, and the penetrating depth is measured by the thread on the rear of the anchor body. Then, the anchor body is pulled out of the medium by the rocker-arm via the pulley. There is a force gauge between the rocker-arm and the pulley; thus, the anchoring force of the anchoring element will be measured. The diameter of anchor body (penetrating unit) used in the test is about 20 mm, and the mass is about 160 g.Figure 4
Testing platform of the anchoring system.The parameters of the media are shown in Table1. Hard clay and fine sand are used in penetrating tests, and soft clay is used in anchoring tests in Section 4.2. Figure 5 shows the prediction of the penetration equation and the corresponding experimental results. It can be seen that the penetration equation can logically describe the correspondence between the penetrating depth and the penetrating speed. As the depth of the experimental media is about 0.5 m, the experiment cannot penetrate the depth that exceeds 0.5 m. Experimental result errors are mainly from the media parameters error and equipment measurement error. Simultaneously, it can be found that the depth increase in fine sand is less than that in hard clay as the penetrating speed increases. This phenomenon is related to the fluidity (characterized by the cohesion) of the fine sand. The smaller the cohesion, the easier it is to flow. Thus, the cavity formed in fine sand as the anchor body penetrating is difficult to maintain. Compared with hard clay, flowing fine sand is easier to consume the penetration energy of the anchor body.Table 1
Parameters of the clay and fine sand.
ParameterDensityρ0 (kg/m3)Internal frictional angleφ (°)Cohesionτ0(kPa)Elasticity moduleE (MPa)Poisson ratioWater contentVoid ratioLocked volumetric strainη∗MediumHard clay1.98×1031416023.80.319.91%0.7090.15Fine sand1.71×10335.50.245.00.3500.5850.1Soft clay1.85×103145.0—————Figure 5
Results of penetration equation and experiments of anchor penetration.


## 3.1. Penetrating Resistance Force Analysis

The anchoring force is directly related to the penetrating depth. In order to establish the relationship between the anchoring force and the initial condition, it is firstly necessary to establish the relationship between the penetrating depth and the initial condition.The forces on arbitrary shape anchor tip are shown in Figure3 in penetrating.Figure 3
Force on arbitrary anchor tip.Based on the cavity expansion model [16–18] and citing the deduction process of the paper [19], we can get the axial force on the anchor tip [18].
(1)Fz1=2π∫0rAτ0+Bρ0Vz211+φ′r2rdr+2πμ∫0zAτ0+Bρ0Vz2f′z21+f′z2fzdz0≤z≤pFz2=2π∫0RAτ0+Bρ0Vz211+φ′r2rdr+2πμ∫0pAτ0+Bρ0Vz2f′z21+f′z2fzdzz>p,where
(2)A=1α1+τ0/2Eγ2α−1λ,B=3/1−η∗1−2α2−α+1/γ21+τ0/2E/γ2α3τ0/E+η∗1−3τ0/2E2−γ321−η∗2−α+3γ3/1−η∗1−2α2−α1+τ0/2E4.α=3λ/3+2λ, where λ=tanφ and φ is the internal friction angle.η∗=1−ρ0/ρ∗, where ρ0 is the volume density before distortion and ρ∗ is the volume density after distortion.
(3)γ=Vc=1+τ02E−1−η∗1/3,where E is the Young’s modulus and τ0 is the cohesion of the medium.


## 3.2. Penetrating Depth

The generatrix equation of the conic anchor tip is as follows:(4)rz=Rp=tanθ,where θ is the semivertical angle of the anchor tip and θ=9.46.Obtaining(5)r=fz=ztanθ,z=φr=rcotθ,f′z=tanθ,φ′r=cotθ.Integrating the above equation with equation (4) will obtain
(6)Fz1=πr2Aτ0+Bρ0Vz2⋅sin2θ+2πμ∫0zAτ0+Bρ0Vz2⋅sin2θztanθdz0≤z≤pFz2=πR2Aτ0+Bρ0Vz2⋅sin2θ+2πμ∫0pAτ0+Bρ0Vz2⋅sin2θztanθdzz>p.(1)
When penetrating depth0≤z≤pTo the conic anchor tip, according to equation (6),
(7)Fz1=πr2Aτ0+Bρ0Vz2⋅sin2θ1+μcotθ0≤z≤p.It can be written as follows:(8)Fz1=πr2Aτ01+μcotθ+πr2Bρ0sin2θ1+μcotθVz2.Substitutingr=ztanθ into the equation above, it is obtained that
(9)Fz1=α1z2+β1z2Vz2,where
(10)α1=πtan2θAτ01+μcotθ,β1=πtan2θBρ0sin2θ1+μcotθ.According to Newton second law, we can get(11)mVzdVzdz=−α1+β1Vz2z2,m2dVz2α1+β1Vz2=−z2dz.Integrating(12)−m2β1lnα1+β1Vz2+C1=13z3.Whent=0, Vz=V0, and z=0, we find
(13)C1=m2β1lnα1+β1V02.SubstitutingC1 into equation (12), we can get
(14)z=3m2β1lnα1+β1V02α1+β1Vz21/3.Whenz=p,
(15)Vzp=α1β1+V02e−2β1p3/3m−α1β1.When maximum penetrating depthh1<p, setting Vz=0, we can get
(16)h1=3m2β1ln1+β1α1V021/3.By the way, according to the Newton second law, we can get(17)md2zdt2=−α1+β1Vz2z2.AsVz=dz/dt, the above equation can be written as
(18)md2zdt2=−α1+β1dzdt2z2.Integrating the above equation, we can get the relation between the depthz and the time t. And then differentiating, we can get the relation between the acceleration and the time. In the paper, we are only interested in the relationship between the initial penetration velocity and the final penetration depth. Thus, the solution to the equation is not discussed in detail.
(2)
When the penetrating depthz>pWhenz>p, we can get the antiforce of the anchor tip by equation (6).
(19)Fz2=πR2Aτ01+μcotθ+πR2Bρ0sin2θ1+μcotθVz2z>p.It can be written as follows:(20)Fz2=α2+β2Vz2where
(21)α2=πR2Aτ01+μcotθ,β2=πR2Bρ0sin2θ1+μcotθ.According to the Newton second law, we can get(22)mVzdVzdz=−α2+β2Vz2.Obtaining(23)m2dVz2α2+β2Vz2=−dz.Integrating(24)−m2β2lnα2+β2Vz2+C2=z.Whenz=p, Vz=Vzp, getting
(25)C2=m2β2lnα2+β2Vzp2+p.SubstitutingC2 into equation (26), we can get
(26)z=m2β2lnα2+β2Vzp2α2+β2Vz2+p.SettingVz=0, we can get the maximum penetrating depth h2h2>p.
(27)h2=m2β2ln1+β2α2Vzp2+p.By the way, according to the Newton second law, we can get(28)mdVzdt=−α2+β2Vz2.Integrating the above equation, we can get(29)Vz=α2β2tan−α2β2mt+C3.Taking the initial conditionVztp=Vzp (tp is the time of z=p) into the above equation, we can obtain the relationship between Vz and the time t. Then, the relationship between the acceleration and the time can be known.


## 3.3. Penetrating Test

In order to verify the abovementioned penetration equation, the relationship between the penetrating depth and the penetrating speed of the anchor body in hard clay and fine sand media is tested experimentally. Testing platform of the anchoring system is shown in Figure4. The anchoring element is fired by the powder, which is ignited via the trigger, and the trigger is controlled by the electric cable. When the anchor body is fired by the anchoring element, the firing speed of the anchor body is measured by the speedometer, and the penetrating depth is measured by the thread on the rear of the anchor body. Then, the anchor body is pulled out of the medium by the rocker-arm via the pulley. There is a force gauge between the rocker-arm and the pulley; thus, the anchoring force of the anchoring element will be measured. The diameter of anchor body (penetrating unit) used in the test is about 20 mm, and the mass is about 160 g.Figure 4
Testing platform of the anchoring system.The parameters of the media are shown in Table1. Hard clay and fine sand are used in penetrating tests, and soft clay is used in anchoring tests in Section 4.2. Figure 5 shows the prediction of the penetration equation and the corresponding experimental results. It can be seen that the penetration equation can logically describe the correspondence between the penetrating depth and the penetrating speed. As the depth of the experimental media is about 0.5 m, the experiment cannot penetrate the depth that exceeds 0.5 m. Experimental result errors are mainly from the media parameters error and equipment measurement error. Simultaneously, it can be found that the depth increase in fine sand is less than that in hard clay as the penetrating speed increases. This phenomenon is related to the fluidity (characterized by the cohesion) of the fine sand. The smaller the cohesion, the easier it is to flow. Thus, the cavity formed in fine sand as the anchor body penetrating is difficult to maintain. Compared with hard clay, flowing fine sand is easier to consume the penetration energy of the anchor body.Table 1
Parameters of the clay and fine sand.
ParameterDensityρ0 (kg/m3)Internal frictional angleφ (°)Cohesionτ0(kPa)Elasticity moduleE (MPa)Poisson ratioWater contentVoid ratioLocked volumetric strainη∗MediumHard clay1.98×1031416023.80.319.91%0.7090.15Fine sand1.71×10335.50.245.00.3500.5850.1Soft clay1.85×103145.0—————Figure 5
Results of penetration equation and experiments of anchor penetration.


## 4. Anchoring Model

### 4.1. Anchoring Force Analysis

After the anchor body has penetrated, the barbs on the anchor tip are opened to produce a greater anchoring force. The shape of the open barbs formed can be approximately regarded to be round and like a disc. There will be shear damage between the round disc anchor tip and the medium. It is found that a circular truncated cone concave was formed on the surface of the medium after the anchor body is pulled out. Thus, an anchoring model shown in Figure6 is established. The maximum principal stress of the element in medium is σ1, which is the passive stress of the Rankine. The minimum principal stress σ3 is the compressive stress caused by the gravity of the upper layer. The depth of the round disc on the anchor tip is H. The medium is sheared layedly when the anchor is pulled out. When the round disc of the anchor tip is pulled out by the distance of H1 and reaches A position, the anchoring force F reaches the maximum value. Then, under the effect of this maximum force, the surface of the medium is subjected to the overall shear failure and forms the circular truncated cone concave. The height of the circular truncated cone concave is H-H1. Its bottom radius is R, and its top radius is equal to the round disc radius on the anchor tip, which is r. The angle between medium surface normal and the concave generatrix is δ=π/4+φ/2 (based on soil mechanics, the angle between the shear failure surface and the maximum principal stress surface is δ=π/4+φ/2). Finally, the anchor tip is separated from the medium by holding the circular truncated cone.Figure 6
Anchoring model of the anchoring system after barbs stretching.The maximum principal stressσ1 and the minimum principal stress σ3 are
(30)σ1=ρgH−ztan2δ+2τ0tanδ,σ3=ρgH−z,where ρ is the density of the medium and g is the gravity coefficient.From the Mohr-Coulomb yield criterion, the shear strength of the medium is(31)τ=τ0+σtanφ.The positive pressureσ of the medium in column shear failure stage (0≤z≤H1) and circular truncated cone shear failure stage (H1<z≤H), respectively, are
(32)σ=σ10≤z≤H1,σ1+σ32+σ1−σ32cos2δH1<z≤H.Thus, the shear strength is as follows:(33)τ1=τ0+σ1tanφ0≤z≤H1,τ2=τ0+σ1+σ32−σ1−σ32sinφtanφH1<z≤H.The anchoring forceF1 generated in column shear failure stage (0≤z≤H1) can be expressed as follows:
(34)F1=∫0z2πrτ1dz.The anchoring forceF2 generated in circular truncated cone shear failure stage (H1<z≤H) can be expressed as follows:
(35)F2=∫H1H2πr+z−H1tanδτ2dzcosδ.According to the assumption of the anchoring model, whenz=H1, the anchoring force reaches the maximum value, and this force causes the shear failure of part (H-H1), resulting in a circular truncated cone concave. Thus,
(36)F1z=H1=F2cosδ.Jointing equations (30) and (33)~(36), we can firstly obtain H1. Then, the value of the anchoring force in 0≤z≤H1 interval can be obtained.When the circular truncated cone is generated after shear failure of the medium (z>H1), the anchoring force F3 is the gravity of the circular truncated cone pulled out by the anchor tip. It can be expressed as follows:
(37)F3=∫H1Hρgπr+z−H1tanδ2dz.Thus, the anchoring forceF in the whole pulling out process of the anchor body can be expressed as follows:
(38)F=F10≤z≤H1,F3z>H1.


### 4.2. Anchoring Test

The anchoring model is verified by testing the anchoring forces at different depths in fine sand and soft clay media. After several times of penetratings, the state of the hard clay is destroyed and its parameters are changed. We remeasure the parameters of the clay and call the clay in this state as soft clay. Moreover, soft clay is convenient for operating in test. The key parameters of the fine sand and the soft clay are shown in Table1. Discs with 90 mm diameter and 40 mm diameter are separately mounted on the top of the anchor body, which is shown in Figure 7. Then, the anchoring forces of the anchor body in two different penetrating depths are tested. There are two penetrating depth values: one is 285 mm which represents full penetrating of the anchor body, and the other is 165 mm representing partial penetrating.Figure 7
40 mm and 90 mm diameter disc mounted on the anchor body.The medium on the anchor body after it is pulled out from the soft soil in 285 mm penetrating depth is shown in Figure8, where the anchoring hole shape is also displayed. As the medium on the anchor body would collapse under the effect of the gravity, it is difficult to find the integral shear shape of the medium on the anchor body at the time of pulling out. Nonetheless, it is obvious that the medium shear shape on the left side of the anchor body in Figure 8 is very similar to that shown in Figure 6.Figure 8
Medium (soil) and its shape when the anchor body is pulled out.When the anchor body is pulled out from the soft clay and the fine sand, the relationship between the anchoring force and the pullout displacement is shown in Figures9 and 10, and the maximum anchoring force can be found in two figures. We can conclude from the testing that the anchoring force increases gradually during the pullout process, and the surface of the media is ruptured suddenly when the anchoring force is up to the maximum (due to strong liquidity of the sand, the phenomenon is not obvious). After the anchoring force decreases to a constant value, this value equals to the total weight of the anchor body, the disc, and the medium on its surface. At the same time, it can be seen that the cohesion of the medium, the internal friction angle, the penetration depth of the anchor body, and the diameter of the disc have great influence on the anchoring force.Figure 9
Anchoring experiments in soft clay.
(a)(b)(c)(d)Figure 10
Anchoring experiments in fine sand (also called soft sand).
(a)(b)(c)(d)The calculated maximum anchoring force can be obtained by substituting the parameters of the experimental medium and the penetrating depth of the anchor body into the anchoring model. The calculated and the testing maximum anchoring force are shown in Table2. It can be found that the calculated maximum anchoring force is close to the maximum anchoring force obtained by the experiment, so the anchoring model is considered as reasonable.Table 2
Comparison of anchoring model with experiments in soft clay and sand.
H=285 mm2r=90 mmH=165 mm2r=90 mmH=285 mm2r=40 mmH=165 mm2r=40 mmSoft clayMaximum anchoring force in the anchoring model348.8 N131.5 N176.1 N67.0 NMaximum anchoring force in the tests355.0 N133.6 N195.4 N70.1 NFine sandMaximum anchoring force in the anchoring model84.8 N20.9 N61.0 N13.7 NMaximum anchoring force in the tests77.0 N27.1 N59.3 N17.5 N


## 4.1. Anchoring Force Analysis

After the anchor body has penetrated, the barbs on the anchor tip are opened to produce a greater anchoring force. The shape of the open barbs formed can be approximately regarded to be round and like a disc. There will be shear damage between the round disc anchor tip and the medium. It is found that a circular truncated cone concave was formed on the surface of the medium after the anchor body is pulled out. Thus, an anchoring model shown in Figure6 is established. The maximum principal stress of the element in medium is σ1, which is the passive stress of the Rankine. The minimum principal stress σ3 is the compressive stress caused by the gravity of the upper layer. The depth of the round disc on the anchor tip is H. The medium is sheared layedly when the anchor is pulled out. When the round disc of the anchor tip is pulled out by the distance of H1 and reaches A position, the anchoring force F reaches the maximum value. Then, under the effect of this maximum force, the surface of the medium is subjected to the overall shear failure and forms the circular truncated cone concave. The height of the circular truncated cone concave is H-H1. Its bottom radius is R, and its top radius is equal to the round disc radius on the anchor tip, which is r. The angle between medium surface normal and the concave generatrix is δ=π/4+φ/2 (based on soil mechanics, the angle between the shear failure surface and the maximum principal stress surface is δ=π/4+φ/2). Finally, the anchor tip is separated from the medium by holding the circular truncated cone.Figure 6
Anchoring model of the anchoring system after barbs stretching.The maximum principal stressσ1 and the minimum principal stress σ3 are
(30)σ1=ρgH−ztan2δ+2τ0tanδ,σ3=ρgH−z,where ρ is the density of the medium and g is the gravity coefficient.From the Mohr-Coulomb yield criterion, the shear strength of the medium is(31)τ=τ0+σtanφ.The positive pressureσ of the medium in column shear failure stage (0≤z≤H1) and circular truncated cone shear failure stage (H1<z≤H), respectively, are
(32)σ=σ10≤z≤H1,σ1+σ32+σ1−σ32cos2δH1<z≤H.Thus, the shear strength is as follows:(33)τ1=τ0+σ1tanφ0≤z≤H1,τ2=τ0+σ1+σ32−σ1−σ32sinφtanφH1<z≤H.The anchoring forceF1 generated in column shear failure stage (0≤z≤H1) can be expressed as follows:
(34)F1=∫0z2πrτ1dz.The anchoring forceF2 generated in circular truncated cone shear failure stage (H1<z≤H) can be expressed as follows:
(35)F2=∫H1H2πr+z−H1tanδτ2dzcosδ.According to the assumption of the anchoring model, whenz=H1, the anchoring force reaches the maximum value, and this force causes the shear failure of part (H-H1), resulting in a circular truncated cone concave. Thus,
(36)F1z=H1=F2cosδ.Jointing equations (30) and (33)~(36), we can firstly obtain H1. Then, the value of the anchoring force in 0≤z≤H1 interval can be obtained.When the circular truncated cone is generated after shear failure of the medium (z>H1), the anchoring force F3 is the gravity of the circular truncated cone pulled out by the anchor tip. It can be expressed as follows:
(37)F3=∫H1Hρgπr+z−H1tanδ2dz.Thus, the anchoring forceF in the whole pulling out process of the anchor body can be expressed as follows:
(38)F=F10≤z≤H1,F3z>H1.


## 4.2. Anchoring Test

The anchoring model is verified by testing the anchoring forces at different depths in fine sand and soft clay media. After several times of penetratings, the state of the hard clay is destroyed and its parameters are changed. We remeasure the parameters of the clay and call the clay in this state as soft clay. Moreover, soft clay is convenient for operating in test. The key parameters of the fine sand and the soft clay are shown in Table1. Discs with 90 mm diameter and 40 mm diameter are separately mounted on the top of the anchor body, which is shown in Figure 7. Then, the anchoring forces of the anchor body in two different penetrating depths are tested. There are two penetrating depth values: one is 285 mm which represents full penetrating of the anchor body, and the other is 165 mm representing partial penetrating.Figure 7
40 mm and 90 mm diameter disc mounted on the anchor body.The medium on the anchor body after it is pulled out from the soft soil in 285 mm penetrating depth is shown in Figure8, where the anchoring hole shape is also displayed. As the medium on the anchor body would collapse under the effect of the gravity, it is difficult to find the integral shear shape of the medium on the anchor body at the time of pulling out. Nonetheless, it is obvious that the medium shear shape on the left side of the anchor body in Figure 8 is very similar to that shown in Figure 6.Figure 8
Medium (soil) and its shape when the anchor body is pulled out.When the anchor body is pulled out from the soft clay and the fine sand, the relationship between the anchoring force and the pullout displacement is shown in Figures9 and 10, and the maximum anchoring force can be found in two figures. We can conclude from the testing that the anchoring force increases gradually during the pullout process, and the surface of the media is ruptured suddenly when the anchoring force is up to the maximum (due to strong liquidity of the sand, the phenomenon is not obvious). After the anchoring force decreases to a constant value, this value equals to the total weight of the anchor body, the disc, and the medium on its surface. At the same time, it can be seen that the cohesion of the medium, the internal friction angle, the penetration depth of the anchor body, and the diameter of the disc have great influence on the anchoring force.Figure 9
Anchoring experiments in soft clay.
(a)(b)(c)(d)Figure 10
Anchoring experiments in fine sand (also called soft sand).
(a)(b)(c)(d)The calculated maximum anchoring force can be obtained by substituting the parameters of the experimental medium and the penetrating depth of the anchor body into the anchoring model. The calculated and the testing maximum anchoring force are shown in Table2. It can be found that the calculated maximum anchoring force is close to the maximum anchoring force obtained by the experiment, so the anchoring model is considered as reasonable.Table 2
Comparison of anchoring model with experiments in soft clay and sand.
H=285 mm2r=90 mmH=165 mm2r=90 mmH=285 mm2r=40 mmH=165 mm2r=40 mmSoft clayMaximum anchoring force in the anchoring model348.8 N131.5 N176.1 N67.0 NMaximum anchoring force in the tests355.0 N133.6 N195.4 N70.1 NFine sandMaximum anchoring force in the anchoring model84.8 N20.9 N61.0 N13.7 NMaximum anchoring force in the tests77.0 N27.1 N59.3 N17.5 N


## 5. Penetrating-Anchoring Mathematical Model

At the time of the asteroid landing mechanism landing, the anchor body of the anchoring system is projected into the asteroid surface at a certain speed, which produces the anchoring force. The anchoring force is the key performance index for the asteroid landing mechanism. Establishing the relationship between the initial penetrating speed of the anchor body and the anchoring force is of great importance.The anchoring model mainly studies the relationship between the penetrating depthH and the anchoring force F, and the penetration equation mainly describes the relationship between the penetration depth H and the initial penetrating velocity V0. Therefore, on the base of the anchoring model and the penetration equation, the mathematical relationship between the initial penetrating velocity V0 and the anchoring force F can be established with the penetration depth H as the link. Through this mathematical relationship, we can easily calculate the minimum initial penetrating speed needed to obtain the required anchoring force.The penetrating depth of the anchor body can be obtained by the penetrating equations (16) and (27), which can be expressed as follows:
(39)H=h10≤z≤p,h2z>p,written as(40)H=f1m,V0,θ,μ,τ0,E,φ,ρ,ρ∗.When the anchor body is anchored at the depthH, the maximum anchoring force occurrence point H1 can be obtained by combining equations (30) and (33)~(36). It is written as follows:
(41)H1=f2H,τ0,φ,ρ,g.And then substituting thez=H1 into equation (34), we can obtain the maximum anchoring force Fmax.
(42)Fmax=f3H,τ0,φ,ρ,g.Substituting equation (40) into equation (42), we obtain
(43)Fmax=f4m,V0,θ,μ,τ0,E,φ,ρ,ρ∗,g.Equation (43) is the penetrating-anchoring model. This model incorporates the penetration equation and the anchoring model with the penetrating depth as a link, which combines two independent processes (penetrating and anchoring) as a whole. This model describes the mathematical relationship among the maximum anchoring force Fmax, the initial penetrating velocity V0, and the medium properties. The penetrating-anchoring model can easily predict the required minimum initial penetrating speed that is adopted for generating the expected maximum anchoring force.


## 6. Conclusions

In the paper, based on the cavity expansion model and the Mohr-Coulomb yield criterion, the force on the anchor body during the penetration is analyzed. A penetrating equation is introduced for calculating the resistance force and penetrating depth of the anchor body when it is penetrating. The penetrating depth estimated by the model is approximate to the testing results in clay and sand media. Utilizing soil mechanics, the force on the deployed anchor body in the process of pulling out is analyzed, and the anchoring model is built. This model describes the relationship between the penetrating depth and the anchoring force. The anchoring force on the anchor body equipped with 40 mm disc and the 90 mm disc which represent deployed anchor tip are tested in the clay and fine sand. The maximum anchoring force obtained by the tests is close to the calculated maximum anchoring force. The maximum anchoring force in the test appears at the time of media surface failure, and the residual medium’s shape on the anchor body after it is extracted is similar to the medium’s shear shape in the anchoring model. The penetrating-anchoring model is built by uniting the anchoring model and the penetration equation with the penetration depth as the link. This model establishes the relationship among the anchoring force, the initial penetrating velocity, and the media properties, and it provides a theoretical foundation for the design and performance estimation of the anchoring system in the asteroid or the comet landing exploration. Tests from different angles or close to the actual situation are the follow-up plan, and we are planning to carry out similar tests in the future.



---

*Source: 8649172-2021-02-20.xml*