Abstract
To better understand the penetration mechanism of the elliptical cross-section projectile (ECSP) into semi-infinite concrete target, penetration experiments using three types of ECSPs with different shape ratios (1, 1.25 and 1.61) and with striking velocities ranged from 550 m/s to 1050 m/s were conducted. Penetration depths, penetration trajectory and mass erosion rates of the projectile were obtained after the experiments. The experiment results show that the penetration performance and ballistic stability of the ECSP are equivalent to those of the circular cross-section projectile (CCSP). Based on the theory of complex variable function and conformal transformation, a semi-analytical model which can calculate the cavity boundary stress distribution of elliptical section cavity controlled by the displacement boundary condition was established and the model was validated by comparing the model degenerate solution with Kirsch problem results. Theoretical calculation results show that the radial stress of elliptical section cavity increases progressively from the minor axis to the major axis. In addition, a formula combining with the semi-analytical theoretical model and the local interaction theory was developed. The predicted penetration depths were compared with 30 groups of experiment data with different projectile parameters and striking velocities and coincide quite well with the corresponding experimental data. Finally, the influence of shape ratio and caliber-radius-head (CRH) on the penetration performance of projectile and the application prospect of ECSPs on hypersonic weapon platform were studied.
Keywords
Elliptical cross-section projectile; Resistance characteristics; Conformal transformation
1 Introduction
The penetration performance of Earth-penetrating weapons on rock or soil media has been a research hotspot for many years. When the inherent parameters of the projectile such as striking velocities, material and mass remain unchanged, changing the projectile shape becomes a significant way to improve penetration performance of the projectile. Compared with the circular cross-section projectile (CCSP), the non-circular cross-section projectile has the advantages of excellent platform adaptability, high bending strength and improving the space utilization rate of the internal cavity of the platform [11 Vedernikov Y.A., Khudyakov Y.S., (2005). Evolution and comparative analysis of group armor and aeroballistics of ancient and modern ruled and poly-wedge arrows. Math Model Geoph 10:93–116.
2 Bless S.J., (1996). Penetration mechanics of non-circular rods. American Institute of Physics 1996;370:1119-1122.-33 Du Z.H., Zhu J.S., Wang X.Z., et al (2009). Analytical Model on Non-circular Penetrator Impacting Semi-infinite Target.Acta Armamentarii 4:403-407. (in chinese)]. Therefore, it is crucial to study the penetration performance, resistance characteristic and penetration mechanism of the non-circular cross-section projectiles, especially the elliptical cross-section projectile (ECSP), into rock or soil media.
For CCSPs, there are a lot of experience and engineering models that have been widely used in penetration problems [44 Yankelevsky D.Z., Adin M.A., (1980). A simplified analytical method for soil penetration analysis. International journal for numerical and analytical methods in geomechanics 4(3):233-254.
5 Feldgun V.R., Yankelevsky D.Z., Karinski Y.S., (2017). A new simplified analytical model for soil penetration analysis of rigid projectiles using the Riemann problem solution. International Journal of Impact Engineering 101:49-65.
6 Luk V.K., Forrestal M.J., (1987). Penetration into semi-infinite reinforced-concrete targets with spherical and ogival nose projectiles. International Journal of Impact Engineering 6(4):291-301.
7 Forrestal M.J., Luk V.K., (1992). Penetration into soil targets. International Journal of Impact Engineering 12(3):427-444.-88 Bishop R.F., Hill R, Mott N.F., (1945). The theory of indentation and hardness experiments. Proceedings of the Physical Society 57(3):147-159.], and can be well predict the penetration depth of projectiles, such as soil-disc model [44 Yankelevsky D.Z., Adin M.A., (1980). A simplified analytical method for soil penetration analysis. International journal for numerical and analytical methods in geomechanics 4(3):233-254., 55 Feldgun V.R., Yankelevsky D.Z., Karinski Y.S., (2017). A new simplified analytical model for soil penetration analysis of rigid projectiles using the Riemann problem solution. International Journal of Impact Engineering 101:49-65.] and cavity expansion model [66 Luk V.K., Forrestal M.J., (1987). Penetration into semi-infinite reinforced-concrete targets with spherical and ogival nose projectiles. International Journal of Impact Engineering 6(4):291-301.
7 Forrestal M.J., Luk V.K., (1992). Penetration into soil targets. International Journal of Impact Engineering 12(3):427-444.-88 Bishop R.F., Hill R, Mott N.F., (1945). The theory of indentation and hardness experiments. Proceedings of the Physical Society 57(3):147-159.]. However, few studies are conducted to clarify the penetration mechanism of ECSPs. The research work of Ben-Dor et al. and Yakunina et al. [99 Ben-Dor G., Dubinsky A., (1997). Optimal 3D impactors penetrating into layered targets. Theor Appl Fract Mech 27(3):161–6.
10 Ben-Dor G., Dubinsky A., (2002). A model for predicting penetration and perforation of FRP laminates by 3-D impactors. Composite Structures 56(3):243-248.
11 Yakunina G.Y., (2020). The optimum non-conical and asymmetrical three-dimensional configurations. JOURNAL OF Mechanics and applied mathematics 64(4):583–91.-1212 Yakunina G.Y., (2000). The construction of optimum three-dimensional shapes within the framework of a model of local interaction. Journal of Applied Mathematics and Mechanics 64(2):289–98.] pointed out that when the projectile structure was formed by the surface of a circular cone and its tangential plane combination, the projectile has the smallest penetration resistance and the highest penetration efficiency. In 1978, Dhaliwal and Singh [1313 Dhaliwal R.S., Singh B.M., (1978). The axisymmetric boussinesq problem of an initially stressed neo-hookean half-space for punch of arbitrary profile. International journal of engineering science 16(6):379-385.] gave the derivation process of the penetration force for the rigid projectile with arbitrary cross-section shape to punch into an elastic object, and proposed a simple closed-form expression of the penetration force required for projectiles with six types of cross-section shapes. Furthermore, HYUNG [1414 HYUNG. Cavity expansion analysis of non-circular cross-sectional penetration problems. Austin: The University of Texas; 1997. p. 56-85.] proposed the cavity expansion problem with arbitrary cross section based on the least square method and concluded that the non-circular cross-section projectile structure can effectively reduce the penetration resistance of the target. However, due to the cavity boundary discretization method, this method cannot deal with continuous variable load problems and the accuracy of the calculation results is low. Bless [22 Bless S.J., (1996). Penetration mechanics of non-circular rods. American Institute of Physics 1996;370:1119-1122.] carried out high-speed penetration tests of rectangular cross-section rod and circular cross-section rod, compared and analyzed their penetration efficiency, and concluded that the penetration performance of asymmetric cross-section projectile was better below the critical velocity. In Petry and ACE (Army Corps of Engineers) [1515 Walter T.A., Wolde-Tinsae A.M., (1984). Turbine Missile Perforation of Reinforced Concrete. Journal of Structural Engineering 110(10):2439–55., 1616 Ben-Dor G., Dubinsky A., et al (2006). Applied high-speed plate penetration dynamics. Springer Netherlands 132(4):213-240.] empirical formulas, the residual velocity of non-circular projectiles is calculated using circular projectiles with the same cross-section area as non-circular projectiles. Wang et al. [1717 Wang W.J., Zhang X.F., Deng J.J., Zheng Y.M., Liu C., (2018). Analysis of projectile penetrating into mortar target with elliptical cross-section. Explosion and Shock Waves 38:164-73. (in chinese)] carried out the penetration experiment of ECSPs and CCSPs, compared and analyzed the penetration performance of the two projectiles, and proposed an analysis model to calculate the penetration depth of the projectile based on the cavity expansion model. Nevertheless, the solution of the projectile penetration is simplified and the resistance model was still the traditional circular cavity expansion model, the model cannot reasonably explain the mechanism of the ECSP penetrating into concrete target.
Wu et al. and Liu et al. [1818 Dong H., LIU Z.H., WU H.J., et al (2019). Study on penetration characteristics of high-speed elliptical cross-sectional projectiles into concrete. International Journal of Impact Engineering 132:1-12.
19 Dong H., WU H.J., et al (2020). Penetration characteristics of pyramidal projectile into concrete target. International Journal of Impact Engineering 143:1-14.-2020 Liu Z.H., Wu H.J., Gao X.D., Pi A.G., Huang F.L., (2019). Study on the resistance characteristics of elliptical cross-section ogive-nose projectile penetrating concrete. Transaction of Beijing Institute of Technology 39:135-6. (in chinese)] carried out a series of penetration experiments of ECSPs penetrating into semi-infinite concrete target in the striking velocity range of 700 m/s to 1000 m/s and analyzed the penetration performance of the ECSP and the damage characteristics of the target. They studied the radial stress distribution along the circumferential direction at the cross-section of projectile nose through numerical simulation, introduced the shape factor of the ECSP to modify the cavity expansion theory, analyzed the variation law of the resistance of the projectile nose, and established the model of the normal penetration depth of the ECSP. However, the shape factor of the ECSP is obtained by fitting the numerical simulation results and its rationality require further verified. And it is unable to make a reasonable explanation for the causes of the circumferential stress distribution difference of the ECSP. Dai et al. [2121 Dai X.H., Wang K.H., Li M.R., et al (2021). Rigid elliptical cross-section ogive-nose projectiles penetration into concrete targets. Defence Technology 7(3):1-12.] studied the penetration characteristics of the ECSP and analyzed the influence of CRH,axial stress, and resistance on penetration performance of the ECSP. Although the penetration mechanism of the ECSP was discussed in many aspects, the resistance expression in the theoretical model was still based on the dynamic cavity expansion theory, and the applicability of the resistance expression on ECSPs remains to be studied.
Practically, the penetration depth of the ECSP can be solved by converting the cross-sectional area of the ECSP into the area of an equivalent circle for engineering applications [1515 Walter T.A., Wolde-Tinsae A.M., (1984). Turbine Missile Perforation of Reinforced Concrete. Journal of Structural Engineering 110(10):2439–55., 1616 Ben-Dor G., Dubinsky A., et al (2006). Applied high-speed plate penetration dynamics. Springer Netherlands 132(4):213-240.]. However, due to the asymmetry of the ECSP, the radial stress distribution of the projectile nose and the response characteristics of the target are quite different from those of the CCSP [1818 Dong H., LIU Z.H., WU H.J., et al (2019). Study on penetration characteristics of high-speed elliptical cross-sectional projectiles into concrete. International Journal of Impact Engineering 132:1-12.
19 Dong H., WU H.J., et al (2020). Penetration characteristics of pyramidal projectile into concrete target. International Journal of Impact Engineering 143:1-14.-2020 Liu Z.H., Wu H.J., Gao X.D., Pi A.G., Huang F.L., (2019). Study on the resistance characteristics of elliptical cross-section ogive-nose projectile penetrating concrete. Transaction of Beijing Institute of Technology 39:135-6. (in chinese)]. The radial stress distribution characteristics of the ECSP nose are contributed to the structural optimization design of the projectile, the stability analysis of the charge and the oblique penetration trajectory of ECSPs. And the response characteristics of the target are crucial for understanding the asymmetric failure mode of the target. Unfortunately, the equivalent method focuses on macroscopic physical quantities such as penetration depth, and lacks discussion and research on the above problems.
As mentioned above, a wealth of research works on the penetration characteristics have been carried out, which promotes the study of the penetration mechanism of the ECSP. Nevertheless, there are still some limitations and more attention should be paid to the following aspects: (1) Experiments of ECSP penetration into semi-infinite concrete targets were conducted to have a better understanding of the penetration characteristics. (2) A semi-analytical model based on the theory of complex variable function and conformal transformation [2222 Zhou H., Kong G.Q., Liu H.L., (2016). Pressure-controlled elliptical cavity expansion under anisotropic initial stress:Elastic solution and its application. Science China-Earth Sciences 7:1-20., 2323 Muskhelishvili NI., (1995). Some Basic Problems of the Mathematical Theory of Elasticity, american mathematical monthly.] was developed, which can calculate the stress distribution characteristics of elliptical section cavity. (3) A formula for predicting penetration depth was developed by taking the striking velocity, shape ratio and CRH of projectile into account combining with the semi-analytical model and the force characteristics on the projectile nose. And the formula was evaluated by the corresponding experimental data. (4) The influences of shape ratio and CRH on the penetration performance of projectile and the application prospect of ECSP on hypersonic weapon platform were analyzed and discussed.
2 Experimental study on the ECSP penetrating concrete
In order to study the penetration performance of ECSPs on concrete targets, a series of penetration experiments were conducted using four kinds of projectiles with shape ratios of 1, 1.25 and 1.56 (C1, E1, E2-1, E2-2). The experiments of four projectiles penetrating semi-infinite concrete target with striking velocities ranged from 550 m/s to 1050 m/s were carried out. After the experiment, the ballistic trajectory was obtained by splitting the concrete target and all the projectiles used in the penetration experiment were recovered.
2.1 Experimental program
The main parameters of the four projectiles are shown in Table 1. Notably, as can be seen from Ref. [2121 Dai X.H., Wang K.H., Li M.R., et al (2021). Rigid elliptical cross-section ogive-nose projectiles penetration into concrete targets. Defence Technology 7(3):1-12.], the value of CRH of an ECSP varies with the azimuth. Namely, it increases progressively from the minor axis to the major axis along the circumference. Therefore, the values of CRH for each projectile at the major axis and the minor axis are noted in Table 1. In the following discussions, the value of CRH of the ECSP represents the CRH at the minor axis, that is, the maximum position on the projectile, and the value of CRH in the other directions can be calculated according to the size parameters of the projectile.
The diameter of C1 projectile is 24 mm and a caliber-radius-head (CRH) of 6.1. The semi-major and semi-minor of E2-1 and E2-2 projectiles are 15 mm and 9.6 mm. Particularly, the C1 projectile and the E2-1 projectile maintain the same cross-section area and mass, and the difference between them is that the shape ratio of cross-section is different, which is used to analyze the influence of the shape ratio of cross-section on the penetration performance.
The main difference between the two elliptical cross-section projectiles is that the mass of the E2-2 projectile is lower than that of the E2-1 projectile, which is realized by adjusting the depth and diameter of the inner hollows in their shanks, therefore the E2-2 projectile can have higher striking velocity in the experiments. The semi-major axis of the E1 projectile is 15 mm, and the semi-minor axis is 12 mm, which is used to study the influence of different shape ratios on the penetration performance. Four kinds of projectile head length are 58.1mm and the total length of the projectile are 144 mm, as shown in Fig. 1.
The experiment projectile is composed of projectile body and bottom-bolt. The projectile body material is 30CrMnSiNi2A, the density is 7.85 g/cm3, and the hardness is 42~45 HRC after quenching. The aluminum alloy is selected as the bottom-bolt material, which makes the projectile closely cooperate with the gun tube in the launch process to achieve excellent air closure effect and ensure the stability of the projectile launch.
The target used in the experiment is a concrete target of ϕ1000 × 1000 mm, and the average compressive strength is 40 MPa. The target is placed on the experiment target frame to ensure that the center height of the target is at the same height as that of the ballistic gun. The striking velocity was verified by the aluminum foil velocity measuring device and a high-speed photography. At the same time, the high-speed photography can observe the initial conditions of projectile entering the target. The experiment layout is shown in Fig. 2. In order to confirm the mechanical properties of concrete target, three cube samples with standard side length of 150 mm were poured and experimented. The compressive strength of the samples meets the design requirements.
2.2 Experiment results and analyses
A total of 16 projectiles were launched in the experiment with striking velocities ranged from 550 m/s to 1050 m/s. The E1 projectile in experiment 6 bounced due to the large incident angle, so a total of 15 effective data were obtained. As shown in Table 2, the projectile striking velocity, penetration depth and mass loss ratio are counted after the experiment. The flying posture and trajectory of the E1 projectile in Test NO.7 are shown in Fig. 3 and Fig. 4, which indicates that the penetration condition can be regarded as normal penetration.
From the damage picture of the target (Fig. 5), it can be seen that the penetration process of ECSPs on concrete targets is similar to that of ordinary CCSPs, which is also composed of a crater stage and a tunnel stage. In addition, under the penetration of the ECSP, the tunnel shapes of ECSPs were similar to elliptical rather than round, which is significantly different from the CCSP. This phenomenon shows that during the penetration process of ECSPs, the materials around the elliptical cross-section cavity have different expansion velocities, which is different from the materials around the circular cross-section cavity. Therefore, the radial stress of the ECSP must be asymmetric.
Fig. 6 is the penetration depth data of 15 projectiles under different striking velocities. Notably, the C1 projectile and E2-1 projectile have the same mass and cross-section area. It can be seen from the figure that the penetration depth of the C1 projectile is larger when the penetration velocity is about 750 m/s, and the average penetration depth is 7.3% higher than that of the E2-1 projectile. When the penetration velocity is about 600 m/s, the penetration depth of the E2-1 projectile is larger, and the average value is 13% higher than that of the C1 projectile. This shows that the ECSPs also have excellent penetration ability at high-speed striking velocity.
Compared with E1 projectiles, the values of CRH of E2-1 projectiles are larger, so the penetration ability of E2-1 projectile is significantly better than that of E1 projectile, which is consistent with the theory of Chen et al. [2424 Chen X.W., Li Q.M., et al (2018). Penetration trajectory of concrete targets by ogived steel projectiles experiments and simulations. International Journal of Impact Engineering 120:202-213., 2525 Chen X.W., (2009). Advances in the penetration/ perforation of rigid projectiles. Advances in Mechanics 39(3):316-351.].
The penetration depth of E2-2 projectiles are generally low than others, because the total mass of E2-2 projectile is 90 g lighter than that of other projectiles. The penetration depth data of E2-2 can be used to verify the subsequent theoretical calculation model in a wide range of striking velocity. These experimental results provide a strong experimental data support for the structural optimization of the projectile. Because the penetration performance of the ECSP is similar to that of the CCSP, it means that the ECSP can be a choice to improve the space utilization of the hypersonic weapon platform. Fig. 7 is a typical recovered projectile after experiment and the relationship between the mass loss rate of the recovered projectile and the striking velocity. It can be seen from the figure that the noses and shanks remained intact and with the increase of the striking velocity, the mass loss rate increases, and the mass erosion rates of all projectiles is within 5%. Therefore, it can be considered that the projectile is in a rigid body penetration within striking velocities ranged from 550 m/s to 1050 m/s.
3 Semi-analytical model of elliptical cavity expansion
3.1 Static elliptical cavity expansion model
3.1.1 Basic assumptions of the model
Considering a semi-infinite plane with an elliptical cavity (Fig. 8), the boundary of the elliptical cavity is controlled by displacement boundary conditions, and the stress state of the plane at infinity is known. Where 、、 and are structural parameters of the initial elliptical cavity and the elliptical cavity after expansion, respectively. is the principal stress at infinity and is the polar angle in the Cylindrical coordinate system. The basic assumptions of the theoretical model are as follows:
-
1
The material area can be approximated as an infinite single link and the deformation of the material is limited to a plane and satisfies the plane strain assumption;
-
2
The plastic zone stress of the material satisfies the Von-Mises yield criterion.
-
3
The shape ratio of elliptical cavity () remains constant during expansion.
-
4
The plastic stress field of the cavity is statically determined, and no elastic unloading of the material [2626 Galin L.A., (1946). Plane elasticplastic problem: plastic regions around circular holes in plates and beams. QUARTERLY JOURNAL OF MECHANICS AND APPLIED MATHEMATICS 10:777-783.
27 Detournay E., (1986). An approximate statical solution of the elastoplastic interface for the problem of Galin with a cohesive-frictional material. International journal of solids and structures 22(12):1435-1454.-2828 Zhuang P.Z., YU H.S., (2019). Two-dimensional elastoplastic analysis of cylindrical cavity problems in Tresca materials. International journal for numerical and analytical methods in geomechanics 2:1-20.].
3.1.2 Determination of the displacement boundary condition c assumptions of the model
The conformal mapping method developed by Muskhelishvili [2323 Muskhelishvili NI., (1995). Some Basic Problems of the Mathematical Theory of Elasticity, american mathematical monthly.] is used to solve the second kind boundary-value problem of elasticity.
As shown in Fig. 9, the conformal mapping technique can transform the complex cavity boundary shape in the physical plane (z plane) into the unit circle in the phase plane ( plane) and simplify the solving process of complex cavity boundary shape problem. r and θ are polar radius and polar angle in the physical plane, ρ and ϑ are polar radius and polar angle in the phase plane. The mapping transformation is
where R and m are real numbers, depending on the elliptic semi-majar axis a and semi-minor axis b.
In order to keep the shape ratio of elliptical cavity in expansion unchanged, the displacement at different positions should be a function of , as shown in Eq. (3)
where, is the displacement along the direction of major axis and .
The displacement boundary condition in the physical plane need to be converted to the phase plane. The variable in the physical plane is different from the variable in the phase plane. The conversion between the two parameters is as follows:
At the cavity boundary (),. Eq. (4) can be rearranged as
Therefore, the displacement boundary condition can be written as (in plane)
Conveniently, the Fourier transform of Eq. (6) were written as
where is the Fourier series expansion expression of the displacement boundary condition, and is the Fourier series coefficient.
3.1.3 Solving the stress function
Next, the stress function satisfying the boundary displacement condition needs to be solved. As no body force is applied, the plane strain problem of elasticity is attributed to solving the biharmonic equation of Airy stress function U.
By introducing complex variable and conjugate complex variable and substituting Eq. (8), the complex function representation of stress component can be obtained:
where, 、 and are the radial stress, hoop stress and shear stress of any point in the plane, respectively; Re is a real part of complex numbers, and are complex potential functions, which can be obtained by Eq. (10)
where, and are stress functions that can be solved by the displacement boundary condition.
According to Ref. [2323 Muskhelishvili NI., (1995). Some Basic Problems of the Mathematical Theory of Elasticity, american mathematical monthly.], we have
where, G is shear modulus of the material, u and v are displacements in the direction and , respectively; is Poisson’s ratio and .
The stress functions and are analytic and can be expressed as
The vectors and represent the pressure applied to the initial cavity wall. In this paper, the two vectors are equal to zero. B and are related to the stress at infinity, and their values in this paper are zero [2323 Muskhelishvili NI., (1995). Some Basic Problems of the Mathematical Theory of Elasticity, american mathematical monthly.]. Additionally, the functions and are analytical at infinity and can be written as Laurent series. From Eq. (12), we obtain
Substituting Eq. (7) into Eq. (11), we have
Multiplying the term of at both sides of Eq. (14) and taking integration of Eq. (14) along the initial cavity wall, we can get
By means of the Cauchy integral formula [2323 Muskhelishvili NI., (1995). Some Basic Problems of the Mathematical Theory of Elasticity, american mathematical monthly.], the integral terms of Eq. (15) can be written as follows
Substituting the Eq. (16) into Eq. (13) results in
Similarly, another stress function can also be obtained from Eq. (14)
Substituting Eqs. (17) and (18) into Eq. (9), we get
The elastic field stress of material can now be solved by Eqs. (19) and (20).
3.1.4 Solving the cavity boundary stress olving the stress function
The region outside the elliptical cavity satisfies the plastic yield criterion. Under the plane strain assumption, the equilibrium equation without body force is
Assuming that the Von-Mises yield criterion is satisfied in the plastic zone, the ela-plastic boundary c can be obtained by substituting Eqs. (19) and (20) into Eq. (22).
Integrating equation (21) in the plastic zone can obtain the plastic field stress distribution. As shown in Fig. 10, and are base vectors in cylindrical coordinates and we have and base vectors which are normal and tangent to the surface, respectively. The actual direction of the stress is along the outer normal direction of the contour, and Eq. (21) needs to be transformed. Assuming that l is the distance along the direction of , then . therefore:
where, and are directional cosines.
Eq. (21) can be rearranged as
where, the angle depends on the geometry of the contour.
On both sides of the elastic-plastic boundary, stress and displacement are continuous. As shown in Fig. 11, assuming that the elastic-plastic boundary contour is , the stress and displacement of the next node on the plastic zone contour can be obtained by iterative equilibrium equation. and can be obtained by
where, .
The plastic field stress of material can now be solved by Eqs. (25) and (26).
3.2 Validation of elliptical cavity expansion model
3.2.1 The displacement boundary condition expressed by the stress function
In the above solving process, the displacement boundary condition is expressed in three forms: the displacement boundary condition derived from assumption 3 (Eq. (6)), the Fourier series expansion of the displacement boundary condition (Eq. (7)) and the displacement boundary condition expressed by the stress function (Eq. (11)). Therefore, when the Fourier series expansion of the displacement boundary condition and the stress function solved by Eqs. (17) and (18) are correct, Eq. (6) 、Eq. (7) and Eq. (11) are three expressions of the same function, and the curves drawn according to three formulas should be consistent. When ,G=7.86GPa, , and , the curves are shown in Fig. 12. It can be seen from the figure that the three curves basically coincide, which proves the correctness of the stress function.
3.2.2 Comparison between model degenerate solution and Kirsch problem results
When a and b in Eq. (1) are equal, the model can map the outside of the circle with radius R to the outside of the unit circle. Using the conformal transformation formula, the above derived formula can be used to obtain the degenerate solution of elliptical cavity expansion-circular cavity expansion. Muskhelishvili [2323 Muskhelishvili NI., (1995). Some Basic Problems of the Mathematical Theory of Elasticity, american mathematical monthly.] gives the analytical solution of the elastic field stress of the material when the circular cavity expands to a certain distance u (Kirsch problem), so the reliability of the solution in this paper is verified by the degenerate solution. The stress components calculated by the method in this paper are
Fig. 13 shows the radial stress , hoop stress and radial displacement along the radial distribution of the calculation results in this paper compared with the calculation results of Muskhelishvili, this is completely consistent with the theoretical solution of Muskhelishvili.
Comparison between the calculation results and Kirsch problem results [2323 Muskhelishvili NI., (1995). Some Basic Problems of the Mathematical Theory of Elasticity, american mathematical monthly.].
3.3 Penetration depth analysis of the ECSP penetrating concrete target
3.3.1 Radial stress distribution characteristics of elliptical cavity boundary
Fig. 14 shows the calculation results of radial stress on the elliptical cavity boundary when the parameters are consistent with section 3.2.1. The radial stress is divided by the yield strength of the material so that the effect of the yield strength of the material on the calculation results can be eliminated, as shown in Fig. 14a.
Dong et al. [1818 Dong H., LIU Z.H., WU H.J., et al (2019). Study on penetration characteristics of high-speed elliptical cross-sectional projectiles into concrete. International Journal of Impact Engineering 132:1-12.] studied the radial stress distribution of ECSP with different striking velocities penetrating concrete target by numerical simulation. The results show that the radial stress at the projectile nose increases progressively from the minor axis to the major axis. In addition, in the process of projectile penetration, the nose of projectile always maintains the process of expanding from small elliptical cavity to large elliptical cavity with equal shape ratio. The displacement required in the major axis of the elliptical cavity is larger than that in the minor axis. Therefore, it is reasonable that the stress in the major axis of the elliptical cavity is larger than that in the minor axis. The calculation results of the semi-analytical model in this paper also have the same law.
Defined the stress distribution coefficient k (Fig. 14b) as the stress distribution of equal proportion expansion of the cavity of elliptic cross-section divided by the stress value required for the cavity of circular cross-section to expand the same distance and assuming that the stress distribution coefficient k of the projectile nose remains unchanged during the whole tunnel phase of penetration. Namely, k is related to the cross-section shape of the projectile and material properties of the target, but not to the striking velocity.
Specifically, with the parameters in section 3.2.1, the means of the stress distribution coefficient k is 3.84 and the stress value required for the cavity of circular cross-section to expand the same distance is 3.66. This means that due to the asymmetric stress distribution on the elliptical cavity boundary, the penetration resistance of the ECSP is larger than that of the CCSP with the same cross-section area
3.3.2 Prediction model of penetration-depth
When the projectile velocity is lower than the critical velocity, in order to simplify the analysis process, the theoretical model generally assumes that the projectile is rigid body. After exceeding the critical velocity, the penetration is dominated by a semi-hydrodynamic regime. From the relevant research results [2929 Liu C., Zhang X.F., Chen H.H., et al (2020). Experimental and theoretical study on steel long-rod projectile penetration into concrete targets with elevated impact velocities[J]. International Journal of Impact Engineering 138:103482.,3030 He X., Xu X.Y., Sun G.J., et al (2010). Experimental investigation on projectiles high-velocity penetration into concrete targets. Explosion Shock Waves 30(1): 1-6.] and the recovered projectile in this paper, it can be concluded that the projectile velocity does not exceed the critical velocity, and is still in the stage of rigid penetration.It is well known that the process of projectile penetrating into concrete target can be divided into two stages. Forrestal et al. [3131 Forrestal M.J., Frew D.J., et al (2003). Penetration of concrete targets with deceleration-time measurements. International Journal of Impact Engineering 28(5):479–97.] proposed that the depth of the crater area can be expressed by 2D when the projectile penetrates the concrete medium. The penetration resistance in the first stage is
where C is the parameter in Ref. [3131 Forrestal M.J., Frew D.J., et al (2003). Penetration of concrete targets with deceleration-time measurements. International Journal of Impact Engineering 28(5):479–97.], D is diameter of the projectile and is the axial position, as shown in Fig. 15.
The resistance calculation model of the ECSP [1818 Dong H., LIU Z.H., WU H.J., et al (2019). Study on penetration characteristics of high-speed elliptical cross-sectional projectiles into concrete. International Journal of Impact Engineering 132:1-12.].
At the second stage, the resistance of the projectile can be obtained by integrating over the projectile nose. is the angle between the axis and the normal direction of the microelement of projectile nose ds. It should be noted that for the CCSP, the value of is constant in the circumferential direction for a fixed penetration position , but for the ECSP, the value of is related to . The axial resistance of the projectile is determined by combining the normal stress and friction force acting on the nose of the projectile. As mentioned earlier, the surface stress of the nose can be expressed by the traditional sphere cavity expansion model (SCE) [3131 Forrestal M.J., Frew D.J., et al (2003). Penetration of concrete targets with deceleration-time measurements. International Journal of Impact Engineering 28(5):479–97.] combined with the stress distribution coefficient k.
where, , and are inertial, viscosity and strength coefficients, respectively [3131 Forrestal M.J., Frew D.J., et al (2003). Penetration of concrete targets with deceleration-time measurements. International Journal of Impact Engineering 28(5):479–97., 3232 Zhang S., et al (2017). High-velocity penetration of concrete targets with three types of projectiles: experiments and analysis. Latin American Journal of Solids and Structures 14(9):1614–28.]. is the density of concrete targets and is the normal velocity of the projectile.
The normal velocity and striking velocity of the projectile have the following relationship: . In the meantime, there is friction in the process of projectile penetration, the relationship between shear stress and friction coefficient is: .
The axial resistance can finally be expressed as
where S is the nose surface area of the projectile, is 0.01 [3333 Yinan W., (2009). The mechanism of high-speed kinetic energy projectile penetration into concrete. Transaction of Beijing Institute of Technology. (in chinese)], and can be obtained by integration
where, , and are integration constant.
The kinematic equation for the projectile with a mass is
Finally, the penetration depth P of the ECSP penetrating concrete target can be calculated by integral Eq. (33).
4 Penetration-depth predictions for circular and elliptical cross-section projectiles
In this section, the accuracy of penetration-depth prediction model in section 3 is verified by the experiment data in section 2 and Refs. [1717 Wang W.J., Zhang X.F., Deng J.J., Zheng Y.M., Liu C., (2018). Analysis of projectile penetrating into mortar target with elliptical cross-section. Explosion and Shock Waves 38:164-73. (in chinese), 1818 Dong H., LIU Z.H., WU H.J., et al (2019). Study on penetration characteristics of high-speed elliptical cross-sectional projectiles into concrete. International Journal of Impact Engineering 132:1-12.]. Table 3 shows the comparison between predictions and the experimental data, and represents the deviation of predictions from the experiment data. The predictions agree well with the experiment data, and the maximum deviation is 14.28%.
Fig. 16 shows the comparison between the predicted penetration depth in this paper and the experiment data with striking velocities ranged from 500 m/s to 1100 m/s. It can be seen from the prediction results that CCSPs (C1) have almost the same mass as ECSPs (E2-1), but in the striking velocity range of 500 m/s to 1100 m/s, the C1 projectile has a slight advantage in penetration depth, and the penetration depth of the CCSP is 6%~8% more than that of the ECSP at the same striking velocity, and with the increase of striking velocity, the difference of penetration depth between two types of projectiles increases progressively.
To verify the accuracy of the model, experimental data [1717 Wang W.J., Zhang X.F., Deng J.J., Zheng Y.M., Liu C., (2018). Analysis of projectile penetrating into mortar target with elliptical cross-section. Explosion and Shock Waves 38:164-73. (in chinese),1818 Dong H., LIU Z.H., WU H.J., et al (2019). Study on penetration characteristics of high-speed elliptical cross-sectional projectiles into concrete. International Journal of Impact Engineering 132:1-12.] are used to compare with the theoretical model and the results are shown in Table 4 and Fig. 17. Although there are differences in the projectile parameters and striking velocities, the penetration-depth predictions agree well with the experiment data and the maximum deviation is 14.56% and the calculation results show that the present model in this paper has good generalization.
5 Analysis and discussions
The shape ratio and CRH of projectile are the main factors that influence the penetration performance at the given projectile mass and target material. In this section, the influence of shape ratio and CRH of projectile on penetration performance are analyzed by taking some projectiles with different structures as computational model based on the validation in section 4. In addition, the ECSP is suitable for hypersonic weapon platform. Base on HIFiRE-4 boost-gliding hypersonic platform, the application prospect of the ECSP is studied.
5.1 Influence of shape ratio on penetration performance
The correctness of the theoretical model has been verified by abundant experimental data. The most representative feature of the ECSP is its shape ratio. In order to study the influence of shape ratio on the penetration performance of ECSP, five types of projectile with shape ratios from 1 to 2 are designed. When the shape ratio is 1, the projectile becomes an ordinary CCSP. Except for shape ratio, the mass, length and cross-section area of the five projectiles are consistent. Fig. 18 shows the penetration depth of five projectiles with striking velocities ranged from 500 m/s to 1200 m/s. At the same striking velocity, the projectile with larger shape ratio has lower penetration depth. In this striking velocity range, the penetration depth of the ECSP with shape ratio of 2 is about 15% lower than that of the CCSP.
The ECSP exhibits lower penetration performance than the CCSP. In order to analyze the reasons for the difference, the radial stress distribution of two typical cross-sections of the five projectiles are extracted, and the cross-section positions are 1/3 and 2/3 of the projectile nose length from the projectile tip, respectively. The striking velocity is set at 800 m/s and the remaining parameters are consistent with section 3. The dimensionless radial stress is shown in Fig. 19, it can be found the radial stress distribution of ECSP is significantly different from that of CCSP, showing a gradual increase from the minor axis direction to the major axis direction, which is similar to the sinusoidal distribution. However, the radial stress distribution is not completely symmetrical. The area surrounded by the curve in Fig. 19 represents the resistance of the section. If the radial stress distribution is completely symmetrical, the penetration performance of the ECSP and the CCSP should be equivalent, because their integral areas are equal.
The ratios of penetration resistance of five projectiles to penetration resistance of CCSP are 1, 1.018, 1.06, 1.12 and 1.19. It can be seen that the penetration resistance of ECSPs is higher than that of CCSPs, and the increase amplitude of penetration resistance increases with the increase of shape ratio. Therefore, under the same conditions, the CCSP has the optimal penetration performance.
5.2 Influence of the CRH on penetration performance
The penetration performance of the CCSP significantly depends on the CRH. As discussed in Refs. [3434 Frew D.J., Forrestal M.J., Cargile J.D., (2006). The effect of concrete target diameter on projectile deceleration and penetration depth. International Journal of Impact Engineering 32:1584-94., 3535 Chen X.W., Li Q.M., (2002). Deep penetration of a rigid projectile with different geometrical characteristics. International Journal of Impact Engineering 27:619-37.], a CCSP with a larger CRH has better penetration performance. In order to study the penetration performance of ECSPs with different CRH, six typical projectiles are designed, and their structural parameters are shown in Table 5. C1, S1 and Z1 represent oval-nosed projectile, sphere-shaped projectile and cone-shaped projectile, respectively. The suffix E represents that the cross-section shape of the projectile is an ellipse with a shape ratio of 2.
The penetration depths of six projectiles calculated by the theoretical model are shown in Fig. 20. For the oval-nosed projectile, the law is consistent with Section 5.1, and the penetration depth of the ECSP (E1-E) is lower than that of the corresponding CCSP (C1). Similarly, the cone-shaped projectile (Z1) is similar to the oval-nosed projectile, and the penetration performance decreases when the cross-section shape of the projectile becomes ellipse (Z1-E), which indicates that the elliptical cross-section affects the force distribution of the projectile, resulting in an increase in the penetration resistance.
However, sphere-shaped projectile (S1) appear to show different laws from oval-nosed projectile and cone-shaped projectile. When the cross-section shape of the projectile becomes ellipse (S1-E), the penetration performance has a slight improvement. It is noted that when the cross-section shape of the projectile becomes ellipse, the value of CRH of oval-nosed projectile and sphere-shaped projectile increases. According to Ref. [3232 Zhang S., et al (2017). High-velocity penetration of concrete targets with three types of projectiles: experiments and analysis. Latin American Journal of Solids and Structures 14(9):1614–28.], the increase of CRH will strengthen the penetration performance of the projectile, which may be the reason why the penetration performance of S1-E projectile is slightly better than that of S1 projectile.
Therefore, two projectiles (C2-E and S2-E) are added to compare and analyze the influence of CRH and shape ratio on penetration performance. C2-E projectile is similar to C1-E projectile, and the shape ratio is 2. By reducing the length of the projectile head, the value of CRH of C2-E projectile is consistent with C1 projectile, and the design method of S2-E projectile is the same. The penetration depths of the six projectiles are shown in Fig. 21. The penetration performance of C2-E projectile and S2-E projectile is the lowest, and their CRH are consistent with that of C1 projectile and S1 projectile, which indicates that after the projectile cross-section shape becomes ellipse, the penetration resistance increases due to the asymmetry of stress distribution. The penetration performance of C1-E and S1-E projectiles is better than that of C2-E and S2-E projectiles because of the influence of CRH. This also explains the reason why the penetration performance of S1-E projectile is slightly better than that of S1 projectile. The increase of CRH improves the penetration performance, and the increase of shape ratio reduces the penetration performance. The value of CRH and shape ratio collectively determine the penetration performance of projectile.
5.3 Applicative prospect
Hypersonic weapon platform is usually a flat waverider structure to meet its aerodynamic performance. If the projectile installed on the weapon platform is an ECSP instead of a CCSP, the installation space can be used more effectively. The mass of the projectile can be increased by increasing its size, and the initial kinetic energy of penetrating the target can be finally increased. Based on the HIFiRE-4 hypersonic weapon platform, a CCSP and an ECSP are designed to meet the size of the platform, which are used to study the application prospect of the ECSP in the hypersonic weapon platform.
The cross-section size of the cabin of HIFiRE-4 hypersonic weapon platform is about 512×192 mm, as shown in Fig.22. The minor axis of the two designed projectiles is tangent to the edge of the platform. The diameter of the CCSP (C3) is 192 mm, while the minor axis length of the ECSP (C3-E) is 192 mm, and the major axis length is 396 mm, and the shape ratio is 2.0625. Both the thickness of projectiles are 20 mm and filled with explosives. The calculation shows that the explosive area of the CCSP and the ECSP is 18145 mm2 and 42499 mm2, respectively, and the explosive area of the ECSP increases by 206%. From the perspective of charge area, due to the unique shape of ECSPs, it has an excellent application prospect in hypersonic weapon platform.
According to the area ratio of the projectile and explosive, the mass ratio of the two projectiles is estimated to be 1.8, assuming that the total mass of the CCSP is 100 kg and the value of CRH is 6.11, so the mass of the ECSP is 180 kg. Notably, due to the strong power of the hypersonic weapon platform, the hypersonic weapon platform is able to load the projectile with larger mass. Therefore, the initial velocity of the ECSP can be equivalent to than the CCSP. The target parameters are consistent with section 4 and the penetration depth of two projectiles with different striking velocities as shown in Fig. 23. It can be seen from the figure that the ECSP (C3-E) makes up for the lack of penetration performance caused by structure through total mass m0, and the penetration performance is better than that of the CCSP (C3), and the gap of penetration depth increases with the increase of striking velocity.
Although the above calculation has highly hypothetical, the calculation results still strongly prove that the ECSP is more suitable for hypersonic weapon platforms in terms of charge quantity and penetration performance.
6 Conclusions
The penetration characteristics of the ECSP into semi-infinite concrete target was studied. Experiments of CCSP and two types of ECSPs with shape ratios of 1.25 and 1.56 were conducted with striking velocities ranged from 550 m/s to 1050 m/s. A formula for predicting penetration depth was developed by combining with the force characteristics on the projectile nose and the semi-analytical model based on the theory of complex variable function and conformal transformation. The experimental and calculation results coincided fairly well. Further researches on the influences of shape ratio and CRH on penetration performance were conducted, and the application prospect of the ECSP on hypersonic weapon platform were studied as well.
-
1
The experimental results showed that ECSPs exhibit excellent penetration behavior and trajectory stability with striking velocities ranged from 550 m/s to 1050 m/s.
-
2
The semi-analytical model controlled by the displacement boundary condition can calculate the radial stress distribution characteristics of elliptical cross-section cavity and the calculation results show that the radial stress of elliptical cross-section cavity increases progressively from the minor axis to the major axis.
-
3
The proposed model was validated by comparing the predicted penetration depths with 30 groups of experiment data, and the maximum deviation was 14.56%.
-
4
With the increase of shape ratio and the decrease of the value of CRH, the penetration performance of the ECSP decreases gradually.
-
5
Based on the HIFiRE-4 hypersonic weapon platform, compared with the CCSP, the mass and charge quantity of the ECSP are increased by 1.8 times and 2.06 times respectively, which has better penetration performance and damage effect.
In conclusion, the research results can be used in the design of non-circular cross-section projectile and solve the ballistic deflection problem of the ECSP penetrating the target. The penetration characteristics of ECSPs, especially the dynamic analytical solution of circumferential stress distribution on the head surface of projectiles and the ballistic trajectory deflection problem, remains to be further explored.
Acknowledgements
This work was supported by the Fundamental Research Funds for the Central Universities [No. 30919011401].
References
-
1Vedernikov Y.A., Khudyakov Y.S., (2005). Evolution and comparative analysis of group armor and aeroballistics of ancient and modern ruled and poly-wedge arrows. Math Model Geoph 10:93–116.
-
2Bless S.J., (1996). Penetration mechanics of non-circular rods. American Institute of Physics 1996;370:1119-1122.
-
3Du Z.H., Zhu J.S., Wang X.Z., et al (2009). Analytical Model on Non-circular Penetrator Impacting Semi-infinite Target.Acta Armamentarii 4:403-407. (in chinese)
-
4Yankelevsky D.Z., Adin M.A., (1980). A simplified analytical method for soil penetration analysis. International journal for numerical and analytical methods in geomechanics 4(3):233-254.
-
5Feldgun V.R., Yankelevsky D.Z., Karinski Y.S., (2017). A new simplified analytical model for soil penetration analysis of rigid projectiles using the Riemann problem solution. International Journal of Impact Engineering 101:49-65.
-
6Luk V.K., Forrestal M.J., (1987). Penetration into semi-infinite reinforced-concrete targets with spherical and ogival nose projectiles. International Journal of Impact Engineering 6(4):291-301.
-
7Forrestal M.J., Luk V.K., (1992). Penetration into soil targets. International Journal of Impact Engineering 12(3):427-444.
-
8Bishop R.F., Hill R, Mott N.F., (1945). The theory of indentation and hardness experiments. Proceedings of the Physical Society 57(3):147-159.
-
9Ben-Dor G., Dubinsky A., (1997). Optimal 3D impactors penetrating into layered targets. Theor Appl Fract Mech 27(3):161–6.
-
10Ben-Dor G., Dubinsky A., (2002). A model for predicting penetration and perforation of FRP laminates by 3-D impactors. Composite Structures 56(3):243-248.
-
11Yakunina G.Y., (2020). The optimum non-conical and asymmetrical three-dimensional configurations. JOURNAL OF Mechanics and applied mathematics 64(4):583–91.
-
12Yakunina G.Y., (2000). The construction of optimum three-dimensional shapes within the framework of a model of local interaction. Journal of Applied Mathematics and Mechanics 64(2):289–98.
-
13Dhaliwal R.S., Singh B.M., (1978). The axisymmetric boussinesq problem of an initially stressed neo-hookean half-space for punch of arbitrary profile. International journal of engineering science 16(6):379-385.
-
14HYUNG. Cavity expansion analysis of non-circular cross-sectional penetration problems. Austin: The University of Texas; 1997. p. 56-85.
-
15Walter T.A., Wolde-Tinsae A.M., (1984). Turbine Missile Perforation of Reinforced Concrete. Journal of Structural Engineering 110(10):2439–55.
-
16Ben-Dor G., Dubinsky A., et al (2006). Applied high-speed plate penetration dynamics. Springer Netherlands 132(4):213-240.
-
17Wang W.J., Zhang X.F., Deng J.J., Zheng Y.M., Liu C., (2018). Analysis of projectile penetrating into mortar target with elliptical cross-section. Explosion and Shock Waves 38:164-73. (in chinese)
-
18Dong H., LIU Z.H., WU H.J., et al (2019). Study on penetration characteristics of high-speed elliptical cross-sectional projectiles into concrete. International Journal of Impact Engineering 132:1-12.
-
19Dong H., WU H.J., et al (2020). Penetration characteristics of pyramidal projectile into concrete target. International Journal of Impact Engineering 143:1-14.
-
20Liu Z.H., Wu H.J., Gao X.D., Pi A.G., Huang F.L., (2019). Study on the resistance characteristics of elliptical cross-section ogive-nose projectile penetrating concrete. Transaction of Beijing Institute of Technology 39:135-6. (in chinese)
-
21Dai X.H., Wang K.H., Li M.R., et al (2021). Rigid elliptical cross-section ogive-nose projectiles penetration into concrete targets. Defence Technology 7(3):1-12.
-
22Zhou H., Kong G.Q., Liu H.L., (2016). Pressure-controlled elliptical cavity expansion under anisotropic initial stress:Elastic solution and its application. Science China-Earth Sciences 7:1-20.
-
23Muskhelishvili NI., (1995). Some Basic Problems of the Mathematical Theory of Elasticity, american mathematical monthly.
-
24Chen X.W., Li Q.M., et al (2018). Penetration trajectory of concrete targets by ogived steel projectiles experiments and simulations. International Journal of Impact Engineering 120:202-213.
-
25Chen X.W., (2009). Advances in the penetration/ perforation of rigid projectiles. Advances in Mechanics 39(3):316-351.
-
26Galin L.A., (1946). Plane elasticplastic problem: plastic regions around circular holes in plates and beams. QUARTERLY JOURNAL OF MECHANICS AND APPLIED MATHEMATICS 10:777-783.
-
27Detournay E., (1986). An approximate statical solution of the elastoplastic interface for the problem of Galin with a cohesive-frictional material. International journal of solids and structures 22(12):1435-1454.
-
28Zhuang P.Z., YU H.S., (2019). Two-dimensional elastoplastic analysis of cylindrical cavity problems in Tresca materials. International journal for numerical and analytical methods in geomechanics 2:1-20.
-
29Liu C., Zhang X.F., Chen H.H., et al (2020). Experimental and theoretical study on steel long-rod projectile penetration into concrete targets with elevated impact velocities[J]. International Journal of Impact Engineering 138:103482.
-
30He X., Xu X.Y., Sun G.J., et al (2010). Experimental investigation on projectiles high-velocity penetration into concrete targets. Explosion Shock Waves 30(1): 1-6.
-
31Forrestal M.J., Frew D.J., et al (2003). Penetration of concrete targets with deceleration-time measurements. International Journal of Impact Engineering 28(5):479–97.
-
32Zhang S., et al (2017). High-velocity penetration of concrete targets with three types of projectiles: experiments and analysis. Latin American Journal of Solids and Structures 14(9):1614–28.
-
33Yinan W., (2009). The mechanism of high-speed kinetic energy projectile penetration into concrete. Transaction of Beijing Institute of Technology. (in chinese)
-
34Frew D.J., Forrestal M.J., Cargile J.D., (2006). The effect of concrete target diameter on projectile deceleration and penetration depth. International Journal of Impact Engineering 32:1584-94.
-
35Chen X.W., Li Q.M., (2002). Deep penetration of a rigid projectile with different geometrical characteristics. International Journal of Impact Engineering 27:619-37.
Edited by
Publication Dates
-
Publication in this collection
15 Apr 2022 -
Date of issue
2022
History
-
Received
05 Jan 2022 -
Accepted
21 Mar 2022