Abstract

High-speed penetrators carrying detection equipment impact planetary bodies at high speeds, and they are therefore buried at depths of up to several meters beneath the surface. During the friction and collision with the crust of the planet, the acceleration of the scientific instrumentation is significantly large. The vibration protection structure for scientific instrumentation is necessary for the reduction of the peak value of the acceleration response and the improvement of the survival rate. In this study, a penetrator with a multilayered energy absorbing structure was developed to improve the survival rate of the penetrator, of which the foam-filled thin-walled structure (FTS) was applied to the penetrating vibration-damping structure. The penetration process of the penetrator into the planetary medium was simulated using the LS-DYNA software platform. The results obtained using empirical formulas and theoretical derivations were compared with the results of the numerical analysis. The reliability of the penetrator limit element model was then verified by conducting an impulse response experiment and simulation. The results suggest that FTS has a positive influence on the isolation impact and energy absorption. Moreover, the vibration isolation effects of nine different FTSs were evaluated with respect to the following six factors: impact isolation efficiency, load efficiency, peak of acceleration, peak impact force, total energy absorption, and specific energy absorption. Furthermore, the design of the damping structure provides an indispensable solution for penetrator detection.

1. Introduction

Deep-space exploration involves the detection of Solar System, Milky Way, and outer space [13]. To better understand Earth, Solar System, and the universe, various research institutions [4, 5] have directed significant research attention toward the improvement of deep-space detection technology. Research on penetration detection was first carried out in the previous century. This form of detection involves more detailed and comprehensive scientific data with a lower detection cost. In particular, for several minor celestial bodies such as Moon and Europa, the detection equipment is difficult to land, which increases the complexity of the detection task and costs. Therefore, it is necessary to conduct in-depth research on the penetrator, for successful adaptation to future deep-space exploration tasks. There are various types of penetrators that carry different payloads, such as heat flow probes, accelerometers, magnetometers, and communication modules. Commonly used penetrators include Mars 96 (Russia), Deep Space 2 (US), Lunar-A (Japan), and MoonLITE (the United Kingdom).

The penetrator is accompanied by strong vibrations and impacts during operation, which result in the high peak acceleration of the onboard scientific instrumentation (modern electronics), which are very sensitive to acceleration. In such a harsh environment, modern electronics are easily damaged. Therefore, protective design is essential for scientific instrumentation, in addition to the optimization of the vibration-damping performance of the protective structure. The main research directions include penetration performance, material and structure of the penetrator, and vibration reduction methods.

The foam-filled thin-walled structure (FTS) is a widely used vibration-damping structure. Dirgantara et al. [6] described the collisional mechanical properties of foam-filled single-wall and double-wall columns. The interaction between the foam core and the column wall changes the deformation mode from one partial fold to multiple propagation folds. Ebrahimi et al. [7, 8] studied the crashworthiness of functionally graded foam-filled tubes and analyzed the foamed aluminum density with respect to total energy absorption and the influence of the peak collision force. Fang et al. [9] developed square and hexagonal cross-sectional foams. The thin-wall collision experiment was characterized by a progressive folding deformation mode, and the contribution of the number of pleats to the energy loss was studied. Yin et al. [10, 11] numerically simulated the energy absorption characteristics of foam-filled multiunit thin-walled structures. The optimal optimization structure was determined by the multiobjective optimization function of the specific energy absorption and the peak force of the members.

2. Materials and Methods

2.1. Kinematic Analysis
2.1.1. Working Principle of Penetrator

After the penetrator is separated from the orbiter in the descending orbit and cruises to the nearest point to the planet, the braking of the rocket ignition forces the penetrator to continuously decelerate. The brake rocket is then separated from the penetrator, and the penetrator continuously descends. Under the control system and gravity, the angle of flight attack should be small, which reduces the contact resistance between the penetrator and the subsurface of the planet, thus obtaining the ideal depth of dive. The penetrator is then separated, and it descends from the orbiter. The schematic diagram of the landing process is shown in Figure 1.

The scientific instrumentation inside the penetrator can detect the temperature, magnetic field, seismic wave, and other scientific objects of the celestial body. Then the detected information will be transmitted to Earth after being received by the satellite [12, 13].

2.1.2. Design of Penetrator

In this paper, a penetrator is designed as the research object, and the FTS is applied to the protective design of the penetrator. The detailed structure of the penetrator is shown in Figure 2.

The main components of the penetrator are an external protection structure and an FTS, as shown in Figure 3, in addition to scientific instrumentation, a damping layer, and a damping spring. The scientific instrumentation includes the detection device and its power module. The external protection structure is a critical protective structure for penetration protection, as it isolates the penetrating medium from the internal structure. It is the part under the most significant amount of stress. The FTS, which acts as a protective buffer during collision and penetration, is critical to the continued functionality of the onboard scientific instrumentation [14].

2.1.3. Young’s Formula

At present, the research methods of projectile penetration exhibit a diversified trend. The main methods are the empirical method, theoretical method, and the numerical analysis method. The empirical method (Young’s formula) and theory analysis method (Forrestal’s formula) can be employed for the prediction of the numerical analysis results.

Researchers in the United States of America summarized the experimental data into an empirical formula referred to as Young’s formula [15], which is a classic formula for projectile penetration into soil, rock, and concrete media, to calculate the penetration depth. In the engineering calculation of the soil penetration, it is still widely used. The specific calculation method can be expressed as follows:where m is the mass of the penetrator; A is the cross-sectional area of the penetrator; is the landing speed; N is the shape coefficient of the penetrator; and S is the index of penetrability, which is related to the hardness of the soil. The coefficient of elasticity is as follows:where h is the length of the projectile and 2a is the diameter of the projectile. In this study, N = 0.83. The penetration index can be expressed as follows:where Ea is the constrained modulus of the soil.

It can be seen from equations (1)–(3) that the penetration depth of the penetrator can be predicted using the shape parameters, in addition to the mass and soil parameters of the penetrator. However, Young’s formula has significant limitations, in which a detailed analysis of the physical quantities such as the velocity and acceleration of the penetrator cannot be conducted.

2.1.4. Forrestal’s Formula

Luk et al. [16, 17] applied the cavity expansion theory to the penetrator that penetrates the soil medium and then obtained the normal stress on the surface of the penetrator during the penetration process:where A and B are the material constants determined by the cavity expansion theory and and are the Tresca yield stress and initial density of the soil medium, respectively. Moreover, is the cavity expansion velocity, i.e., the normal velocity of the projectile surface.

The penetration depth of the limit can be expressed as follows:where P is the limit depth, m is the total mass, and are the penetration parameters, and is the landing speed.

The empirical formula and the theoretical derivation method can only be used to study the penetrating performance. The research on the impact isolation effect can be easily determined using a numerical simulation [1820].

2.2. Numerical Simulation
2.2.1. Penetration Model

The penetrator model was modeled in Hypermesh/LS-DYNA to obtain the finite element model shown in Figure 4. The contact between the outer protective shell and the soil medium was defined by the ERODING_SURFACE_TO_SURFACE contact type, the contact type between the internal scientific instrumentations was defined by the AUTOMATIC_ SURFACE _TO_SURFACE contact type, and the friction coefficient was set as 0.1 [2125]. The thin-walled material was aluminum alloy AA6060-T4, and the filling material was aluminum foam. A hexahedral element was used to divide the soil target into grids. The element size of the soil was set as 30 mm, thus totaling 1,188,000 elements. The whole model consists of 1,366,100 elements and 1,467,944 nodes.

2.2.2. Material Modeling

The material parameters of each part of the penetrator can be described as follows:

(1) Aluminum Foam. The hardening equation of aluminum foam [26] is a function of yield stress, plateau stress, and equivalent strain, which can be expressed as follows:where , , and are the material parameters; is the equivalent strain; is the yield stress; is the platform stress; and is the compact equivalent strain, which is related to the density of the foamed aluminum:where is the density of aluminum foam; is the density of the foamed aluminum-based material; is the plastic Poisson’s ratio; and is the material parameter, which is related to the compressive deformation volume and deformation form.

The equivalent strain is calculated as follows:where the strain path parameter . Moreover, the detailed material parameters of aluminum foam specimens with different densities are shown in Table 1.

(2) Aluminum Thin Wall. The tapered thin-walled structure is shown in Figure 5, in which the half of the taper angle is 5°. The material was AA6060-T4, the modulus of elasticity was 64.86 GPa, Poisson’s ratio was 0.28, the yield strength was 71 MPa, and the ultimate stress was 164 MPa. The engineering stress-strain curve of the material is shown in Figure 6.

(3) Soil. The soil model was a cuboid with dimensions of 3 m × 4 m × 0.22 m to ensure that the penetrator stops. The deformation of soil was modeled using the MAT-147 MAT_FHWA_SOIL model on the LS-DYNA software platform. The model surface was a modified Mohr–Coulomb surface to determine the pressure-dependent peak shear strength. The specific parameters of the model are shown in Table 2.

(4) Others. The material of the external protection structure of the penetrator is alloy steel, and the material of the scientific instrument is aluminum alloy 7075. The material of the damping layer is polyurethane. The material properties are shown in Table 3.

3. Results and Discussion

3.1. Validation of Numerical Models
3.1.1. Shock Response Experiment

The experimental results validated the reliability of the numerical analysis results. Due to the complexity of the penetrator penetration experiment, significant preliminary preparation work was required. From the simulation results, the impact acceleration of the penetration process exceeded 10,000 g. However, the maximum impact acceleration of the experimental equipment was only 3000 g. Therefore, this experiment was only used as a verification of the reliability of the finite element model.

(1) Observation Points. The selection of two observation points was required during the test to facilitate the recording of the acceleration of the thin-walled structure and the FTS.

The selection of the observation points was carried out for a simple comparison between the data of the two observation points and to ensure that the sensor has a large attachment area (i.e., the sensor is firmly fixed during the impact test). The layout of the acceleration sensor is shown in Figure 7.

(2) Experimental Equipment and Impact Input. The penetrator and the resonance plate were fastened together using bolts, and the impact hammer was set up to strike the resonance plate after the pneumatic rod is pushed, as shown in Figure 8.

The input conditions of the impact test bench are shown in Figure 9. The peak value of the impact acceleration was set as 3000 g (in the frequency domain), and the rising slope was 6 db/oct. The impact response spectrum test was conducted on a penetrator with thin wall thicknesses of 0.2 mm, 0.3 mm, and 0.4 mm and a foam-filled aluminum thin wall with a thickness of 0.2 mm.

(3) Experimental Result. The deformation of the thin wall and FTS are shown in Figure 10. Given that the experimental input conditions of the four structural members were equal, the acceleration signals detected by using the acceleration sensor were similar in the time domain. The output result of one structural member is shown in Figure 11.

The acceleration was found to decay rapidly within 10 ms, and the peak acceleration reached 4000 g. This acceleration response was then used as an input to the numerical simulation analysis to obtain the acceleration response of the scientific instrumentation, i.e., the acceleration response of Point 2. The test response of Test Point 2 is presented with respect to the simulation results in Figure 12.

3.1.2. Numerical Simulation of Shock Response Spectrum Experiment

The finite element model is shown in Figure 13. The penetrator model in this model was the same as that in the penetration simulation. The acceleration signal of a point obtained by the test was inputted to the penetrator model along the direction of the penetrator axis. The simulated acceleration results were then compared with the test results to determine the accuracy of the response results [14].

4. Results

The deformation, which represents the length of the eroded above four structures along the axis direction, is shown in Figure 14. Table 4 presents a comparative data sheet for the simulation and experimental deformation, with a maximum error of 9.3%.

The acceleration of Point 2 is shown in Figure 12. Although the damping effect of the FTS was the least significant, the error between the numerical simulation and experimental results was small. Moreover, under the condition of 10,000 g, the numerical simulation results indicated that the impact absorption effect of the FTS improved and that the crushing of the thin wall would occur.

The results of the finite element analysis exhibited the same general trend as that of the test results. The time and peak size of the transient response were not significantly different, with an error of less than 9%, as shown in Table 5. The results revealed that the finite element analysis results contained a reference value, thus validating the reliability of the finite element model established in this study.

4.1. Simulation Results
4.1.1. Penetration Collision Process

The penetration depth of the penetrator is related to the mission requirements, the type of payload and the initial velocity, among others. The penetrator was designed to penetrate the surface of the planet by a depth greater than 4 m for the detection of the subsurface composition and changes in temperature. The penetration process of the soil medium at different initial velocities was simulated in this study. The trajectory of the penetrator with an initial velocity of 150 m/s in the soil is shown in Figure 15.

From Figure 16, it can be seen that the ultimate dive depth of the penetrator is related to the landing speed. When the landing speed is less than 280 m/s, the detection task is difficult to complete; thus, the landing speed of the selected penetrator was 300 m/s. In this study, the simulation results were compared with the results of the empirical method and the theory analysis method, as shown in Figure 17.

The results indicated that the trend of the changes in the calculation results of the theoretical method, empirical method, and simulation was the same, with an error of less than 5%, which verifies the accuracy of the penetration model established in this study, in addition to the reliability of the results.

4.1.2. Evaluation of Buffering Effect

Due to the influence of the numerical analysis calculation error, the initial velocity  = 300 m/s was set. In addition, 6060-T4 thin walls with thicknesses of t = 2.5 mm, 3 mm, and 3.5 mm, which had material densities of 0.17 g/cm3, 0.35 g/cm3, and 0.51 g/cm3, respectively, were selected for analysis. The deformation and energy absorption of the FTS with different density filling materials were studied. The acceleration produced by the landing resistance of a high-speed object allows for the object to withstand the impact of acceleration. The isolation system plays a role in reducing its maximum acceleration. Significant acceleration is the most critical indicator for the evaluation of the buffering performance of the high-speed impact penetrator.

(1) Acceleration Response. The vibration reduction effect of the FTS is directly reflected in the peak acceleration and the impact isolation efficiency of the internal load [29]. Moreover, it is possible to convert between the two parameters as follows:where is the impact isolation efficiency, is the peak of the impact acceleration response of scientific instrumentation, and is the peak of the impact acceleration response of the outer protective structure.

During the impact, the outer protective structure is first affected. The overload curve is shown in Figure 18(a). The maximum acceleration peak of the protective structure exceeded 10,000 g. With an increase in the penetration depth, the acceleration amplitude rapidly decayed.

The acceleration-time curve for scientific instrumentations under FTS protection is shown in Figures 18(b)18(d). It can be seen from the results that there were variations in the accelerated response curve of the scientific instrumentation. The magnitude of the variations reflects the isolation impact efficiency. With a decrease in the variation of the acceleration curve, there was an improvement in the buffer absorption effect. The instability of the acceleration is due to the fold within the thin wall and the compression deformation of the aluminum foam. The acceleration peaks of the S2.5A0.17 and S2.5A0.35 scientific instrumentations were excessively large, in addition to the variations, which is not a suitable design solution. For S2.5A0.51 and S3.0A0.35, the scientific instrumentation was found to have a stable acceleration and a minimum peak acceleration, which is suitable as a buffer structure.

The impact isolation efficiency of each structure is shown in Table 6. In addition, S2.5A0.51 was found to have the best barrier effect, with an impact isolation efficiency of 77.3%. S2.5A0.17 demonstrated the least significant barrier effect, with an impact isolation efficiency of 27.4%.

(2) Deformation of FTS. The deformation distance of the FTS represents the length of the eroded FTS along the axis direction. Aluminum foams of different densities and thin walls of different thicknesses were numbered in combination. The thin wall was labeled “S,” and the aluminum foam was labeled “A.” For example, S3.0A0.35 represents the FTS with a thin wall thickness of 3.0 mm and the aluminum foam with a density of 0.35 g/cm3. The deformation of each FTS is shown in Figure 19.

The upper end of the thin wall was in contact with the scientific instrumentation. Different types of deformation occur in the thin wall during impact. As shown in Figure 19, the thickness of the thin wall and the density of the aluminum foam have a significant influence on the impact deformation. With an increase in the density of the aluminum foam, the deformation distance (axial compression distance) of the FTS gradually decreased, for the same wall thickness. With an increase in the wall thickness, the deformation distance of the FTS decreased, for the same aluminum foam density. To ensure the buffering energy absorption effect of the FTS, the structure is required to ensure a sufficient deformation distance without crushing. Figures 19(a) and 19(b) reveal that the thin wall was crushed. Figures 19(c)19(f) reveal that the thin wall deformation was the most significant, by which buffering energy effect was realized. Figures 19(g)19(i) reveal that the thin wall deformation was not significant; thus, the expected energy absorption effect was not realized. The deformation distance data are shown in Table 7. The analysis result revealed that S2.5A0.51 and S3.0A0.17 demonstrated an improved deformation distance, and it can therefore be used as an alternative.

(3) Crushing Force. The crushing force between the scientific instrumentation and the FTS can be evaluated using peak crushing force, average crushing force, and load efficiency. Under an axial load, the peak crushing force of the filled structure occurs at two moments of thin wall fold and foam aluminum compaction [30].

The total energy absorption is equal to the work of the deformation due to the crushing force:where EA is the total energy absorption; F is the crushing force; is the deformation distance, which is a function of time; and d is the final deformation distance, which is the steady state value of .

The average crushing force can be expressed as follows:where is the average crushing force, which reflects the magnitude of the inertial force during the crushing, and it is related to the instance of the crushing.

The loading efficiency can be expressed as follows:where is the maximum crushing force.

The crushing force under different combinations is shown in Figure 20. The results indicated that there were significant variations in the impact forces of S2.5A0.17, S2.5A0.35, and S3.0A0.17 with large collision force peaks; and the load efficiencies were 40.3%, 40.6%, and 47.1% with low energy absorption efficiencies. The collision force peaks of S2.5A0.51, S3.0A0.35, S3.0A0.51, S3.5A0.35, and S3.5A0.51 were approximately 200 kN, and the load efficiency exceeded 74%, in which the energy absorption rate is better and is a suitable and efficient combination of the energy absorption structure. The load efficiency of the FTS is shown in Table 8.

(4) Specific Energy Absorption. In the optimization of the structural performance, it is preferable to achieve the desired goal with a smaller volume or quality. In this study, the specific energy absorption was used as an objective function of optimization. The specific energy absorption [31] is defined as follows:where SEA is the specific energy absorption, EA is the total energy absorption, and M is the mass of the energy absorption structure.

In the range of the wall thickness and density of the foam-filled aluminum, the energy absorption and specific energy absorption of the structure exhibited a decreasing trend in accordance with an increase in the thickness and density. The results are shown in Figures 21 and 22. The absorption energy of the comparison of S2.5A0.17 and S2.5A0.35 with the crush were obviously reduced, so the absorption energy is not much different. However, the other two sets of simulation data can be concluded that for the same thickness wall of the structural, with an increase in the density of the foam-filled aluminum, the total energy absorption decreased.

Overall, the performance of S2.5A0.51 was superior with respect to the barrier efficiency, load efficiency, peak acceleration, peak collision energy, total energy absorption, and specific energy absorption.

5. Conclusions

In this study, a penetrator with a multilayer damping structure was designed, which includes rubber, polyurethane, spring and aluminum foam, a thin wall, and other damping materials. The FTS was the main vibration-damping structure. In addition, the effects of the wall thickness and packing density on the vibration isolation and energy absorption effects of the FTS were also quantitatively studied. The conclusions of this paper are as follows.

According to the experimental results of the impact response, the penetrator finite element model was accurate. In addition, based on a comparison between the simulation results and the results of the empirical method, in addition to the theoretical analysis methods during the analysis of the penetration process of the penetrator, the analysis results of the scientific instrumentations during the penetration process were validated. The ultimate penetration depth at different landing speeds was obtained. The landing speed of the penetrator was determined to be 300 m/s, according to the detection requirements.

The vibration isolation effects of nine different FTSs were evaluated using the following six factors: the impact isolation efficiency, load efficiency, peak acceleration, peak impact force, total energy absorption, and specific energy absorption. Although the nine structures can alleviate the negative effects of the impact on the scientific instrumentation, the isolation effects were different. With respect to the evaluation criteria, the combined structure with a thin wall thickness of 2.5 mm and aluminum foam with a density of 0.51 g/cm3 demonstrated the best performance.

The FTS can effectively reduce the overload shock of the scientific instrumentation and improve the survival probability. The shock response test cannot reach the working environment of the penetrator, and the ground test is needed. The penetration process is a thermal-mechanical coupling process, and the design of the penetrator needs to consider this part in future.

Data Availability

The data from this manuscript may be made available upon request to the authors.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant nos. 51975567 and 51505470), State Key Laboratory of Robotics (Grant no. Y7A1207301), Youth Innovation Promotion Association, CAS (Grant no. 2018237), and “JXS” innovation fund, SIA (Grant no. 20180504).