1. Introduction
The traditional launch tube layout is difficult to meet the diversified use needs when launching underwater vehicles. In the future, aerial delivery or rocket-assisted delivery will become more common. The underwater vehicle will cross the gas–liquid interface at high speed and be subjected to huge impact load when using these two methods, which will damage the shell and internal instruments. Therefore, it is urgent and necessary to carry out research on the impact and load reduction of vehicles entering water.
The research methods of water entry impact mainly include theoretical derivation, numerical calculation and experimental testing. The research on water inflow theory is mostly based on potential flow theory, and the approximate plate theory was first proposed by von Karman [
1] in 1929, which became the theoretical basis for many scholars to study related issues later. Wagner [
2] improved the theory and proposed a self-similar solution, which concluded that the peak value of the water inflow impact pressure is proportional to the square of the velocity. Yu et al. [
3,
4,
5] have further improved and perfected Wagner’s theory. Von Karman and Wagner’s water entry theory is only applicable to the initial stage of water entry of structures with small angle of rise, but not the water entry impact problem of flat plate structures. The boundary element method developed since then is applicable to structures with large angle of rise. Park [
6] studied the water entry impact load of a cylinder based on the inviscid potential flow theory. Wu et al. [
7] developed the similarity theory of expansive water body proposed by Longuet Higgins and studied the water entry problem of structures by using the self-similarity method. Jalalisendi et al. [
8] used Donnel’s thin shell theory to establish the model of the cylindrical shell and calculated the water entry load of the cylindrical shell based on Wagner’s theory, and its theoretical calculation results were basically consistent with the test results. Mirzaei [
9], based on the principle of independent expansion of cavitation and considering the memory effect of cavitation, established a theoretical model to solve the cavitation shape of structures entering water. Wang et al. [
10] proposed a theoretical analytical model of elastic wedge impact load and realized high-precision and rapid prediction of structural response based on modal superposition method.
The theoretical research, which is mostly based on the assumption of the actual situation, has certain limitations. With the development of computational fluid dynamics theory, researchers began to widely use numerical simulation methods to carry out water entry research on various structures. At present, the developed numerical methods for water inflow mainly include the grid method, meshless method and hybrid method. In 1974, Hirt [
11] proposed the arbitrary Lagrange–Euler (ALE) method to solve the fluid structure coupling problem with structural deformation and liquid surface movement. Sun et al. carried out a series of studies on the stress load and deformation of small balls during water entry [
12,
13,
14,
15]. Xie et al. carried out a series of studies on water impact problems of wedges and hulls [
16,
17,
18,
19,
20]. Li et al. [
21,
22,
23,
24,
25,
26,
27] carried out diversified studies on the water entry of cylinders, torpedoes and other revolving structures. Among them, Li [
21] calculated the scale model of MK46 torpedo and obtained its impact load laws under different shell materials, different head types and different water entry velocities. Fang et al. [
22] conducted a numerical study on the vertical water entry process of high-speed projectiles with different head types and obtained the conclusion that the peak value of water entry resistance coefficient decreases in the order of flat head, intercepting head, intercepting cone head, cone head and ball head. Song et al. [
23] studied the low-speed inclined water entry process of the revolving body with different head cone angles based on the finite volume method and found that the minimum pressure value in the flow field raises with the increase of the head cone angle. Shi et al. [
26] used the ALE method to conduct numerical simulation on the water entry process of flat-headed and tapered head cylinders with a diameter of 40 mm, and obtained the change laws of head deformation, pressure distribution, cavitation morphology and water entry motion state.
Although numerical calculation is widely used in water entry research, experiments are still needed to verify its correctness and accuracy. Worthington [
28], a British scholar, used flash photography technology to observe the splash and cavitation effects of a small ball entering the water for the first time. The definitions are still used today. Yang et al. [
29] conducted a low-speed water entry test on 10 mm diameter projectiles with different head types. The results showed that round-nosed projectiles were not easy to form cavitation. For the case of a 120° cone nose projectile entering the water vertically, the higher the velocity is, the more likely the surface closure phenomenon and complete detached bubbles will occur. Lu et al. [
30] conducted an experimental study on the cavitation and ballistic characteristics of projectiles with different heads inclined into water at low speed and found that the cavitation diameter raises with the increase of the head cone angle. Lu et al. [
31] carried out a low-speed vertical water entry experiment on streamlined and flat-headed revolving bodies with a diameter of 50 mm, and obtained the variation laws of bubble closure, collapse, jet and splash with the water entry speed.
When vehicles enter water at high speed, it is necessary to take buffering and load shedding measures to ensure the safety of the structure. The researchers took various measures to reduce the water entry load, such as optimizing the shape design and installing elastic expansion parts, air bags and buffer caps on the head of the vehicle [
32]. Liu et al. [
33] adopted the method of head air injection to reduce the impact load of a cylinder with a diameter of 7 mm when entering the water at low speed. Wei et al. [
34] studied the impact of different jet structures on high-speed water entry impact of vehicles and found that the load reduction effect of nozzles with equal diameters was the best. Compared with other schemes, the buffer head cap has higher practical value. Hinckley et al. [
35,
36,
37,
38,
39,
40,
41,
42] carried out a lot of research on the head cap with different buffer materials. Among them, Xuan et al. [
37] carried out an experimental study on the disassembly of the head cap when it entered water at a speed of 26 m/s, which explained the cause of the head cap’s fracture and proved that the buffer head cap can significantly reduce the impact load. Qian et al. [
38] studied the reducing effect of polyurethane foam with different densities on water entry stress through experiments and simulations. Li et al. [
39] designed a composite structure buffer to reduce the load during the high-speed water entry (100 m/s) of a 250 mm diameter vehicle. They respectively used rigid polyurethane foam, polystyrene and aluminum foam as the buffer materials; studied the buffer performance of different materials; and found that aluminum foam had better buffer performance when entering the water at high speed. Sun et al. [
40] carried out a study on the cushioning performance of aluminum foam for the water entry load reduction of a 324 mm diameter vehicle and found that the load reduction performance first increased and then decreased with the raise of the relative density of aluminum foam. Wei et al. [
41] carried out a study on the delamination of the buffer foam in the process of high-speed water entry (150 m/s) of a 324 mm diameter vehicle and found that the layered buffer foam was not easy to be completely destroyed; additionally, there was a secondary buffering phenomenon, and the load reduction performance was better than that of the single-layer foam. Shi et al. [
42] carried out gradient density head cap design and load reduction performance research to solve the problem of load reduction of vehicles with a diameter of 534 mm entering water at high speed. It is found that the energy absorption effect raised with the increase of the number of layers, and the buffer performance of negative density gradient arrangement buffers was stronger than that of positive density gradient.
At present, although there are a lot of research studies on different head types entering the water, the research object is basically the model with smaller diameter and low entering speed. The water entry impact load is strongly related to the head shape, water impact area and water entry velocity. It will rise significantly when the inlet velocity increase or the impact area increase while the model’s diameter remains the same. Therefore, this paper carried out the numerical simulation of vehicles (diameter 200 mm) with flat and round heads entering the water at high speed (100 m/s–175 m/s) to obtain the effects of entering speed and angle on the impact load. In terms of load shedding research, although a large number of scholars have carried out systematic research on various buffer caps, most of them are theoretical and numerical studies. Limited by the experimental conditions, the existing experimental studies mostly consider low-speed water entry, and the experimental objects are mostly small models that do not have enough space to install large test systems. Large test systems have continuous test capability and sensors with larger range, so the model needs larger space to accommodate them. Furthermore, when the diameter of the model increases, the mass of the model raises dramatically, which requires higher launch capacity. At the same time, the improvement of the water entry speed of the model also requires a higher launching capacity. Therefore, this paper uses a high-speed water entry launcher to accelerate the 200 mm diameter model with a large test system inside, and carries out an experimental study on the model with polyurethane foam at the head entering the water at high speed (60 m/s–100 m/s) to obtain the load reduction performance of the polyurethane head cap. The paper is expected to be a reference for the research on the load reduction of vehicles entering the water at high speed.
3. Analysis of Load Reduction Characteristics of Head Cap
According to the simulation results, the load reduction ability brought by the head type is strongly related to the water entry angle. The load reduction effect is obvious only when entering water vertically, making the load reduction effect not stable. However, it is difficult to ensure the vertical entry attitude of underwater vehicles during the launching process. Therefore, in engineering practice, more attention is paid to the reliable load reduction scheme when there is oblique entry. The buffer head cap is a kind of load reduction device installed on the head of the vehicle, which can absorb the impact energy through the buffer material in the process of water entry. A lot of research has been carried out on the buffer head cap, and it has been applied to the engineering practice of low-speed (50 m/s) water entry. However, limited by the experimental conditions, the experimental studies on the high-speed water entry of the head cap are limited, and the energy absorption and breaking process of the head cap still remain unclear. Therefore, this paper carried out an experimental study on the load reduction performance of a 200 mm diameter vehicle with buffer caps when entering water at high speed (100 m/s). The damage mechanism of buffer head cap can be intuitively obtained through analyzing the damage degree of polyurethane buffer material under different working conditions.
3.1. Test Conditions
In order to test the load reduction and energy absorption performance of the buffer head cap, the processed head cap is shown in
Figure 10. The material properties of the three head caps are the same, all of which are polyurethane foam with a density of 70 kg/m
3. They are numbered as No. 1, No. 2 and No. 3, respectively. The length of head cap No. 1 is 200 mm, and the length of head cap No. 2 and No. 3 is 300 mm. As shown in
Figure 11, use No. 88 glue to fix the head cap on the head of the 200 mm diameter model. During this period, use clamp and tooling to keep them stable. After the glue is dried, an experimental model with a head cap is obtained.
As shown in
Table 2, a total of seven tests were carried out, numbered T01–T07 according to the sequence of test conditions, in which T01–T03 are variable pressure tests, T04 is the repeatability test of the launcher and T05–T07 is the load reduction performance test with buffer head caps.
3.2. Analysis of Experimental Results
Taking T01 as an example, the acceleration data of axial (
Ox axis), normal (
Oy axis) and lateral (
Oz axis) are recorded using a three-axis accelerometer during the experiment. The results obtained after Butterworth low-pass filtering (cut-off frequency 4.5kHz) are shown in
Figure 12a. Within 0–0.2 s, the model is still in the launch tube. Within 0.2–0.4 s, the model starts to accelerate under the action of high-pressure gas until entering the water. After the completion of the water entry process, the model will automatically rise to the surface, and then the test is completed. Among them, the water entry process is the key data process to be analyzed. Therefore, this moment is taken as the initial time 0, and the acceleration change processes in the previous and subsequent 50 ms are analyzed, respectively.
Figure 12b shows the test results of axial acceleration. The peak value of axial load is −95.97 g, and the peak pulse width is 0.82 ms.
Figure 12c shows the test results of normal acceleration. The peak value of normal load is −14.3 g, and the peak pulse width is 0.74 ms.
Figure 12d shows the test results of lateral acceleration. The peak value of lateral load is −56 g, and the peak pulse width is 0.96 ms. The impact acceleration of the test model in three directions at the moment of entering the water has the characteristics of large peak value and short pulse width, which is a typical impact signal. It also proves that the process of entering the water is a transient process with dramatic changes in the mechanical environment. The axial acceleration is obviously greater than the normal and lateral acceleration. Theoretically, the lateral acceleration should be zero, but the peak value of the lateral acceleration obtained is slightly greater than the normal value. It is analyzed that the model has a certain roll motion when accelerating in the launch tube. At the moment of entering the water, an included angle is generated between the normal and lateral acceleration sensors on the model and the defined normal and lateral directions. The data processing of other working conditions is similar to T01 and will not be described repeatedly. The specific load peak value and pulse width results are shown in
Table 3.
Due to rolling phenomenon, normal and lateral data have no obvious regularity; therefore, the results of axial acceleration in
Table 3 are mainly analyzed. The following conclusions can be drawn:
(1) Comparing the results of T02 and T04, in which launch speed and angle are basically the same, the relative error of axial test results is about 0.05%. It can be concluded that the test system has good consistency and repeatability.
(2) Comparing the results of T01–T03, it can be seen that the peak value of axial load rises significantly with the increase of water inlet velocity, which is the same as the conclusion obtained by simulation. The ratio of the peak axial load of T03 to T01 is 2.58, and the ratio of the square velocity of the two working conditions is 2.59, which is consistent with Wagner’s water entry theory; that is, the water entry impact is proportional to the square velocity. The peak load pulse width of axial is basically the same, which is about 0.8 ms.
(3) Compare the results of T02 and T05 (with No. 1 head cap). The peak axial load of T02 is −185.9 g, and the peak pulse width is 0.81 ms. The peak axial load of T05 is −177.1 g, and the peak pulse width is 1.29 ms. After using No. 1 buffer head cap, the peak value of axial load decreases by 4.7%, and the corresponding peak pulse width increases by 0.48 ms. We can learn from the relevant knowledge of impact dynamics that the key parameters to describe a typical half sine pulse energy are the impact peak value and its pulse width. The existence of the buffer head cap reduces the peak value, increases the peak pulse width and absorbs part of the impact energy acting on the vehicle. This proves that the head cap can effectively absorb and reduce the impact load acting on the vehicle.
Figure 13 shows the cavitation morphology diagram of T05 (with No. 1 head cap) in the water entry process, and the time interval of each diagram is 5 ms. Within 15 ms (
Figure 13b–d) after entering the water, the head cap deflects towards the upstream side of the vehicle under the impact of water. Then the head cap moves to the lower left under the combined action of buoyancy and initial inertia and separates from the vehicle (
Figure 13e). After that, the head cap and its debris generate irregular cavitation near the water surface. The cavitation generated by the vehicle, and the head cap debris continues to expand. The cavitation is spindle shaped as a whole and is wrapped by numerous small air bubbles (
Figure 13f–l).
Figure 14 compares the changes of No. 1 head cap before and after entering the water. It can be seen from the figure that under the impact load, the tip of the head cap is squeezed and slightly damaged. The contact area between the bottom surface of the head cap and the front-end surface of the vehicle appears obviously peeling after being squeezed.
(4) Compare the results of T02 and T06 (with No. 2 head cap). The peak axial load of T02 is −185.9 g, and the peak pulse width is 0.81 ms. The peak axial load of T06 is −151.5 g, and the peak pulse width is 1.31 ms. After using No. 2 buffer head cap, the peak value of axial load decreases by 18.5%, and the corresponding peak pulse width increases by 0.5 ms.
Figure 15 compares the changes of No. 2 head cap before and after entering the water. It can be seen from the figure that under the impact load, the tip of the head cap is damaged. The contact area between the bottom of the head cap and the front-end face of the vehicle is significantly compressed. The foam material is compacted, and large cracks appear with extrusion.
Compare the results of two groups using No. 1 and No. 2 head caps. The water entry speed of the two groups is basically the same, but the length of the buffer head cap is different. When the length of the buffer head cap is increased from 200 mm to 300 mm, the load reduction rate increases from 4.7% to 18.5%, indicating that increasing the length of the buffer head cap can raise the load reduction effect.
(5) Compare the results of T03 and T07 (with No. 3 head cap). The peak axial load of T03 is −247.3 g, and the peak pulse width is 0.75 ms, while the peak axial load of T07 is −179.9 g, and the peak pulse width is 1.27 ms. After using No. 3 buffer head cap, the peak value of axial load decreased by 27.3%, and the corresponding peak pulse width is increased by 0.52 ms.
Figure 16 shows the change of No. 3 head cap before and after entering the water. It can be seen from the figure that the head cap is split into upper, middle and lower parts under the impact load. Obvious erosion can be observed on the surface of the upper and middle parts. Most of the back of the middle part is stripped from the body. Half of the lower part is broken into small pieces due to the impact load; thus, only half can be collected.
Compare the results of two groups of tests using No. 2 and No. 3 head caps. The buffer head caps used in the two groups of tests are identical, but the water entry speed is different. The raise of the water entry speed increases the impact load, leading to the sharply increase of damage degree of the head cap, and the load reduction rate also increases to a certain extent.
4. Conclusions
Based on ALE method, the numerical simulation of high-speed water entry of flat-nosed and round-nosed vehicles is carried out, the variation rules of cavitation and axial acceleration are analyzed and the influence of head shape on water entry impact load is analyzed. The following conclusions are drawn:
(1) When the vehicle enters the water at high speed, the cavitation continues to expand in a spindle shape within 15 ms. The maximum diameter of the cavitation is always near the liquid surface, and the change rate shows a trend of first increasing and then decreasing, while the velocity attenuation is not obvious.
(2) The peak value of the axial acceleration of the vehicle is greatly affected by the water entry velocity and angle, while the stable value of the axial acceleration is greatly affected by the water entry velocity.
(3) Flat-nosed and round-nosed vehicles have different impact areas at different water entry angles; therefore, when the water entry angle is greater than 80°, the impact load on round-nosed vehicles is significantly smaller than that on flat-nosed vehicles, and the peak axial load can be reduced by 22% when entering the water vertically at 100 m/s.
The way of changing the head shape of the vehicle to reduce the water entry load depends on the water entry speed and angle, which makes the load reduction effect unstable. Therefore, further experimental research on the load reduction performance of the vehicle with buffer head cap entering the water at high speed is carried out, and the following conclusions are obtained:
(1) The peak value of the axial load of the vehicle increases with the water entry velocity, which is consistent with Wagner’s water entry theory. However, the peak pulse width is basically the same, both being about 0.8 ms.
(2) The front end of the buffer head cap is damaged after being impacted by water in the process of entering the water. The contact area between the bottom surface of the head cap and the front-end surface of the vehicle shows significant compression deformation. The foam material is compacted, and large damages and cracks occur along with the squeezing action. In this process, part of the water impact energy can be absorbed, so that the axial load peak value is reduced and the peak pulse width is increased, which proves that the buffer head cap has an obvious load reduction effect.
(3) The load reduction rate increases from 4.7% to 18.5% when the length of the buffer head cap is increased from 200 mm to 300 mm, while the water inlet speed is the same, which indicates that increasing the length of the buffer head cap can improve the load reduction effect.
(4) The damage extent of the buffer head cap increases sharply when the water entry speed is raised while the length of the buffer head cap is the same. Due to the increase of the water entry impact load, the damage ranges from cracks and damage at low speed to splitting into three parts, which means that most of the middle and lower parts are broken into small pieces at high speed. The load reduction rate rises to a certain extent when the water entry speed increases.