Abstract
This study uses the unified strength theory to analyse the elastoplastic stage and plastic stage of a linear strain-hardening target material while considering the effects of the intermediate principal stress and the free lateral boundaries of the target. In this investigation, analytical solutions of the radial stress in the cavity wall are obtained, and a unified penetration model of the target material is built. On this basis, penetration resistance formulas and penetration depth formulas for rigid projectiles with various nose shapes penetrating into thick, finite-radius, metallic targets are deduced, the solutions of which are obtained by utilizing the Simpson method. Accordingly, the proposed method offers a broader scope of application and higher accuracy than previous methods. Through this method, a series of analytical solutions based on different criteria can be obtained, and the penetration depth ranges of targets under different striking velocities can be effectively predicted. Moreover, penetration processes under different conditions are numerically simulated using the software ANSYS/LS-DYNA to study the motion law of the projectiles and the dynamic response of the targets. From the theoretical and numerical approaches, a list of influencing factors for terminal ballistic effects are analysed, including the strength criterion differences, the strength parameter b, the striking velocity , the projectile nose shape, and the target radius-to-projectile radius ratio rt/a. The results indicate that, as b changes from 1 to 0, the penetration depth Dmax increases by 22.45%. Additionally, Dmax increases by 40.76% when rt/a changes from 16 to 4; hence, it cannot be calculated as an unlimited-size target anymore when rt/a ≤ 16. In weapons field tests, the radius of the metallic target can be conservatively designed to be greater than 28 times the projectile radius to ignore the effect from the free lateral boundaries of the target.
1. Introduction
The cavity expansion theory which has been widely used in the field of terminal ballistics can be sorted into two theories according to the different hypothesized cavity expansion styles: spherical cavity expansion theory and cylindrical cavity expansion theory. By applying these two theories respectively, investigators have studied penetration phenomena in several materials, such as metals, rocks, concrete, and soil, and many results have been obtained [1–10]. However, most studies on penetration have focused on semi-infinite targets or the finite-thickness targets (considering the free-surface on the back of the target), and systematic studies on penetration into thick, finite-radius, targets have been limited to experiments [11, 12]. The few theoretical studies on penetration mechanics in such targets have generally ignored the effects of the lateral boundaries, and the experiments conducted by Littlefield [11] have shown that this assumption makes the results obviously deviate from reality when the in-plane dimensions of the target are small; that is, the ballistic performance of certain targets is known to be sensitive to confinement. On the other hand, the armoured targets encountered in practical applications and the real targets in experiments often have finite dimensions; hence, edge effects are bound to exist. Accordingly, the finite cylindrical cavity expansion theory was proposed by Jiang et al. [13] in 2011; moreover, they conducted the first systematic investigation on the penetration of long-rod projectiles into thick metal targets with finite in-plane dimensions. Then, an engineering model for the penetration of rigid sharp-nosed projectiles into thick, finite-radius, metal targets was presented by Song et al. [14]. However, this model adopting the von Mises yield criterion is imperfect and can only be used for projectiles with a single nose shape and for materials with an ultimate shear strength that is 0.577 times its yield strength in tension or compression. Even if adopting the Tresca yield criterion or Twin shear stress yield criterion, this model is only applied for materials with an ultimate shear strength that is 0.5 or 0.667 times its yield strength in tension or compression, respectively. Hence, the application scope of this model is limited.
Since the unified strength theory, which considers all the stress components acting on a twin-shear element and their different effects on material failure, can comprehensively reflect the basic strength characteristics of different target materials and the influence of different strength criteria [15, 16], it is widely used as a more reasonable strength criterion for solving penetration problems with complex stress states [8–10]. Note that the establishment and selection of material strength criteria are crucial when studying the antipenetration performance of targets. Hence, to expand the application scope of materials and fully utilize their potential, this paper uses the unified strength theory to investigate the elastic-plastic stage and plastic stage of linear strain-hardening target materials while considering the effects of the intermediate principal stress and the free lateral boundaries of the target based on the finite cylindrical cavity expansion theory which neglects the effect of the inertia terms and the strain rate of the target material. Through this study, analytical solutions of the radial stress in the cavity wall are obtained, and a unified penetration model for linear strain-hardening target materials is built. On this basis, penetration resistance formulas and penetration depth formulas are presented for rigid projectiles penetrating into thick, finite-radius, metallic targets. Thus, the model presented in this work, which is also suitable for the penetration of rigid projectiles into semi-infinite metallic targets, has a broader scope of application and higher accuracy than previous models. Through simulations of the penetration process under different conditions with the finite element software ANSYS/LS-DYNA, the motion law of projectiles and the dynamic response of targets are studied. The analytical results in this paper are compared with the experimental and analytical results in relevant documents, which reveals that the solution in [14] is only a special case of the solution in this paper. Moreover, a series of analytical solutions based on different strength criteria are obtained to effectively predict the penetration depth ranges of all kinds of targets with equal tensile and compressive strength under different striking velocities. As the model proposed herein can be used for sharp-nosed projectiles, conical-nosed projectiles, and spherical-nosed projectiles, the effect of nose shape on the penetration depth is analysed. Moreover, the effects of the strength parameter, strength criterion, striking velocity, and target radius on the penetration depth limit are discussed and quantified. The suggestion that the condition of ignoring the effect of the free lateral boundaries of the target in target dimension design for weapons field tests is put forward through quantifying this effect on penetration performance in various working conditions. The conclusions of this study can be used as a reference of designing for metallic armor protective structures such as tanks and ships.
2. Twin Shear Unified Strength Theory
The twin shear unified strength theory is a new structural theory that considers the effects of the intermediate principal stress and can be applied to a variety of different materials. This theory can be mathematically expressed as [17]where , , and are the maximum principal stress, intermediate principal stress, and minimum principal stress in the unit, respectively; is the tensile strength-to-compressive strength ratio; is a failure criterion parameter that reflects the influence of the intermediate principal stress, wherein ; and , , and are the tensile yield strength, compressive yield strength and shear yield strength of the material, respectively.
When taking , equations (1a) and (1b) can be simplified as
Since this strength theory includes infinitely many yield criteria, it is also known as the twin shear unified yield criterion [17].
3. Finite Cylindrical Cavity Expansion Theory
3.1. Analytical Model Based on the Unified Strength Theory
Figure 1 shows the finite cylindrical cavity expansion model [13]. The target is assumed to have a radius of , a cavity radius of at time , and a final cavity radius of . Note that (c is the elastic wave velocity) is the radius of the wavefront for elastic waves, is the radius of the elastoplastic interface, and is a constant. For incompressible materials, we can take the elastic wave velocity as and Poisson’s ratio as . When the cavity radius increases from to , the entire expansion can be decomposed into two parts: the elastoplastic stage and the plastic stage . It is cracks beginning to form along the surfaces of the cylinder that can be seen as the end of the plastic stage.
Since the cylindrical cavity expansion model belongs to axisymmetrical plane strain problems, is the intermediate main stress, where is the intermediate main stress parameter (0 ≤ m ≤ 1). From generalized Hooke’s law, one obtains in the elastic zone and m ≈ 1 in the plastic zone (i.e., [18]). For incompressible materials, it is suggested that because the volumetric strain is 0. If the pressure is positive, then , , and according to σ1 ≥ σ2 ≥ σ3. Because σ2 ≤ (σ1 + σ3)/2, the equivalent stress and strain of metal materials derived from (2a) can be written aswhere and are the equivalent stress and strain based on the unified strength theory, respectively, and are the radial stress and strain, respectively, and and are the hoop stress and strain, respectively.
From the constitutive relations of linear hardening materials, one obtains
Combining equations (3a) and (3b) and (4a) and (4b), the relationship of effective stress-strain for linear hardening materials can be written as
Geometric relations in the elastic zone (small strain) and the plastic zone (large strain) are expressed as [13]where and are the spatial coordinates and displacement at time , respectively.
If the material density is assumed to be constant, from mass conservation, one will obtain
Integrating the above equation and using the cavity boundary condition gives the displacement field for the particle, which can be expressed as
In the elastic zone, combining (6) and (9), one will obtain
Substituting (10) into (5a) and (5b), the relationship of effective stress-strain for linear hardening materials is given as
The velocity field of the cylindrical cavity expansion model is expressed as [13]
The momentum conservation equation of the target in the cylindrical coordinate system can be expressed aswhere is the radius of the cylinder, is the time, and is the density of the target material.
Equations (6), (7), and (9) and ((11a), (11b))∼(13) are the fundamental equations of the finite cylindrical cavity expansion model for incompressible linear hardening materials based on unified strength theory.
3.2. Cavity Expansion Stress Calculation
3.2.1. Elastoplastic Stage
Because of the continuity of elastoplastic boundary stresses, substituting into (11a) and (11b), one obtains
Setting , the radius of cavity expansion at the end of the first stage is expressed as
In the plastic zone , substituting (11a), (11b), and (12) into (13), we have
Integrating (16) yieldswhere is a constant.
Considering the boundary condition of the cavity wall , can be written as
Then, substituting (18) into (17), the radial stress is given as
Considering the elastoplastic boundary condition , the radial stress in the cavity wall in the elastoplastic stage derived from (19) can be written aswhere is the radial stress in the plastic zone side of the elastoplastic boundary.
In the elastic zone , substituting (11a), (11b), and (12) into (13), we have
Integrating (21), one obtainswhere is a constant.
Considering the free lateral boundaries of the cylinder and (22), one can obtain
Substituting (23) into (22), we obtain the radial stress in the elastic zone as
Considering the elastoplastic boundary condition and (19), one obtainswhere is the radial stress in the elastic zone side of the elastoplastic boundary.
Because of the continuity conditions , substituting (25) into (20), one obtains
Substituting (14) into (26), the radial stress in the cavity wall in the elastoplastic stage can be written as
3.2.2. Plastic Stage
Substituting the boundary condition into (19), the radial stress in the cavity wall in the plastic stage can be given aswhere is the cavity radius at the end of the second stage.
According to when the boundary condition , we obtain from (1a) and (1b), and the equivalent stress based on unified strength theory can be computed as
Based on the hypothesis of plastic work [19, 20] and the stress-strain state of the finite cylindrical cavity expansion problem, one obtainswhere is the equivalent strain based on unified strength theory. Combining (29) and (30), this equivalent strain can be given as
Substituting (9) into (7), the strain components can be expressed as
Combining (31) and (32) and considering the fracture criterion in [21], the cavity radius at the end of the second stage when can be expressed aswhere is the fracture strain of the target material under uniaxial tension.
3.3. Penetration of Rigid Projectiles into a Thick Cylindrical Metallic Target
3.3.1. Penetration Model Analysis for Rigid Projectiles
Suppose that the generatrix equation for the rotating warhead of a rigid sharp-nosed projectile is expressed as , where and are the radius and nose length of the projectile, respectively, as shown in Figure 2. Accordingly, a simulation of the penetration process of the projectile can be regarded as a series of finite cylindrical cavity expansion models with a cavity radius of and a target radius of [14]. Assuming the initial velocity of the projectile is and the initial penetration depth of the target is , the velocity of the projectile is and the penetration depth of the target is at time t; the radius of the target cavity and the velocity of cavity expansion at time t can be expressed as [14]
3.3.2. Penetration Resistance Calculation
Ignoring the variation in resistance when D < l, when D ≥ l, the penetration resistance can be expressed as [14]where is the angle between the tangent of a point on the projectile arc and the projectile axis and is the radial stress in the surface of the projectile calculated by the analytical model in this paper.
Rigid projectiles can be divided into three types according to their nose shape: conical-nosed projectiles, sharp-nosed projectiles, and spherical-nosed projectiles, as shown in Figure 3. For conical-nosed projectiles, sharp-nosed projectiles with CRH = 3 (CRH is the ratio between the arc radius of the projectile nose and the diameter of the projectile), and spherical-nosed projectiles, the cavity radius can be expressed as , , and , respectively.
(a)
(b)
(c)
Substituting (27) and (28) into (35), one obtainswhere and .
Supposing at the end of the elastoplastic stage, a ≤ rc1; hence, the penetration resistance can be calculated by (36a) because the finite cylindrical cavity expansion model only undergoes the elastoplastic stage when x1 ≥ l. In contrast, when a > rc1, the penetration resistance can be calculated in sections by (36a) and (36b) because the finite cylindrical cavity expansion model undergoes the elastoplastic stage and the plastic stage fully when x1 < l.
Combining (34a) and (34b) and (36a) and (36b), formulas with the same form as those in [14] can be given aswhere and are the static strength resistance and flow resistance, respectively; , , , and are the static strength coefficients related to the performance of the projectile and target materials in the elastoplastic stage; , , , , , and are the static strength coefficients related to the performance of the projectile and target materials in the plastic stage; and is the coefficient of the flow resistance. Note that projectiles with different nose shapes have different values of .
The coefficients in (37) and (38) can be calculated aswhere can be determined with the method of digital integration using the Simpson formula.
3.4. Penetration Depth Calculation
Supposing the projectile mass is , according to [14], the penetration depth of a projectile penetrating into a finite-radius metallic target at time is derived from Newton’s second law and the initial conditions, which can be expressed as
Considering the penetration terminating condition , the above formula can be written as
3.5. Experimental Verification
Substituting the corresponding test data in [12] into the formulas in this paper, for the target, rt = 127 mm (i.e., the target radius of 6061-T6511 aluminium alloy in the test), E = 68.9 GPa (elastic modulus), Ep = 46 MPa (tangent modulus), σoy = 365 MPa (initial yield stress), and ρ = 2710 kg/m3 (density). Comparisons between the analytical results and the test results in the document are summarized in Table 1.
Figure 4 compares the ultimate penetration depth results calculated with the formulas in this paper when , the test results in [12] and the results calculated with the Forrestal formula [14] for different values of the projectile striking velocity . The mean relative deviation between the results from the theoretical formula derived in this paper and the test results is 0.0094, and the corresponding mean square deviation is 0.0427. In contrast, the mean relative deviation between the Forrestal formula results [14] and the test results is 0.1469, and the corresponding mean square deviation is 0.0669. Accordingly, the theoretical values in this paper have a 13.75% smaller deviation from the test results than the Forrestal formula results. Hence, the results calculated with the formulas in this paper match the test results well and have higher precision than the results calculated with the Forrestal formula.
3.6. Parameter Discussion
3.6.1. Strength Parameter b and Striking Velocity
Since different materials have different values of b, we can obtain a series of different penetration depth and penetration resistance solutions for which the unified strength theory can be applied to various materials. Figure 5 shows the relation curves between the ultimate penetration depth and the striking velocity for different values of the strength parameter b, which reveals that the strength parameter b has obvious influences on the ultimate penetration depth: the higher the value of b is, the more obvious the intermediate principal stress effect is and the smaller the value of is. Hence, the value of b can objectively indicate the strength potential of the material and ensure the penetration resistance of the material is fully utilized due to the consideration of the intermediate principal stress effect. Thus, this parameter should not be neglected in penetration analyses and associated calculations; otherwise, the accuracy of precision guided weapon designs will be affected.
Because the unified strength theory parameter b represents the selection of different strength criteria, when b takes different values, the unified strength theory can degenerate into different strength criteria (e.g., it can degenerate into the Tresca yield criterion when b = 0, the von Mises yield criterion when b = 0.366, and the Twin shear stress yield criterion when b = 1), for which there are large differences among the computed results. Taking = 569 m/s as an example, the penetration depth calculated based on the Tresca yield criterion is 22.45% higher than that calculated based on the Twin shear yield criterion . This finding shows that the selection of the strength criteria plays an important role in predicting the penetration effects. Thus, a suitable strength criterion should be selected for practical engineering designs.
For the penetration of thick, finite-radius, metallic targets under special conditions, unlike other calculation methods with a unique solution, the calculation method in this paper can obtain a series of analytical solutions. The results in [14] (von Mises yield criterion results), which are only applicable for materials with , are only a special case of the results in this paper when . Moreover, nearly all of the experimental values fall within the range of the abovementioned series of analytical solutions. Therefore, under special conditions, the lower limit and upper limit of the penetration depth according to the calculation method in this paper can be used to effectively predict the range of penetration depth. For instance, the penetration depth ranges of all kinds of targets with equal tensile strength and compressive strength under different striking velocities are determined, and the results are summarized in Table 2.
In addition, from Figure 5, we can also see that when ≤ 1147 m/s (meeting the rigid penetration conditions in this paper), the higher the striking velocity of the projectile is, the greater the penetration depth is.
Figure 6 shows the penetration depth versus the projectile penetration velocity for different projectile striking velocities , from which we can see that as the penetration process elapses, the projectile velocity continuously decreases and the penetration depth continuously increases; however, the rate of increase in the penetration depth decreases over time.
3.6.2. Ratio of the Target Radius to the Projectile Radius
Hereafter, for comparison with the ballistic tests results, the variation in the penetration depth is studied under different projectile striking velocities and target radii (i.e., different ) while keeping the other conditions the same as those in the tests, as shown in Figure 7. During the calculation, when rt/a ≤ 10, a ≤ rc1, and the finite cylindrical cavity expansion model undergoes the elastoplastic stage and the plastic stage fully; hence, the penetration resistance can be calculated in sections by (36a) and (36b). In contrast, when rt/a > 10 and a > rc1, and the finite cylindrical cavity expansion model only undergoes the elastoplastic stage; hence, the penetration resistance can be calculated by (36a). Figure 7 shows that as the target radius-to-projectile radius ratio increases, the penetration depth-to-projectile radius ratio decreases, but the rate of decrease slows when rt/a ≥ 16. Compared to the penetration depth when rt/a = 16, the penetration depth when rt/a = 70 is reduced by only 3.15%, whereas when rt/a = 4, the penetration depth is increased by 40.76%. These findings indicate that the influence of the target boundary size cannot be neglected (i.e., the calculations cannot still assume that the target has an unlimited-size) when rt/a < 16, which is in good agreement with the experimental results in [12].
3.6.3. Projectile Nose Shape
Figure 8 shows a comparison of the ultimate penetration depth Dmax of a sharp-nosed projectile, a conical-nosed projectile, and a spherical-nosed projectile (the diameter of each column body is the same) under different projectile striking velocities . It is quite evident that the projectile nose shape has an obvious influence on the antipenetration performance of thick, finite-radius, metallic targets. The penetration depth of the conical-nosed projectile is greater than that of the sharp-nosed projectile and is substantially greater than that of the spherical-nosed projectile. When the striking velocity is less than approximately 570 m/s, the penetration depth of the sharp-nosed projectile and that of the conical-nosed projectile are approximately equal. Although the column body diameters and lengths are equal, the spherical-nosed projectile, which has a flat nose and a relatively large surface area, experiences a greater amount of resistance during the penetration process than the other two projectile types. This indicates that the penetration capacities of the sharp-nosed projectile and conical-nosed projectile are obviously much stronger than that of the spherical-nosed projectile under the same conditions, and the sharp-nosed projectile has the strongest penetration capacity. There is little difference between the penetration capacity of the sharp-nosed projectile and that of the conical-nosed projectile when the striking velocity is less than approximately 570 m/s.
Figure 8 also shows that all penetration depths increase with increasing projectile striking velocity. However, the penetration depths of the sharp-nosed projectile and conical-nosed projectile are more substantially affected by the striking velocity than that of the spherical-nosed projectile.
4. Numerical Simulation
To further analyse the effect of the free lateral boundaries of the target on penetration into thick, finite-radius, metallic targets and investigate the projectile/target morphologies, the penetration process, and the instantaneous response, the software LS-DYNA is applied to carry out three-dimensional numerical calculations, which can greatly reduce experimental expenses and allow data and phenomena to be observed that cannot be obtained through other methods [22].
4.1. Model Verification and Results Analysis
The target centreline is the initial contact point between the projectile and the target, leaving a gap between them to eliminate the initial permeability. Considering the amount and efficiency of calculation, quarter-models of the entities are established by using the Lagrange method for the numerical simulations based on the ballistic test parameters in [12], as shown in Figure 9 (take member 1 for example). SOLID164 elements are adopted to mesh the projectile and the target; these elements have eight nodes, wherein every node has nine degrees of freedom. For the projectile, there are 11657 nodes and 9728 elements in total. For the target, there are 96128 nodes and 88900 elements in total, which divide the target evenly along the thickness in 127 layers and divide the target unevenly along the cross section in 40 regions. The meshing is fine near the centre of the target, whereas it is coarse near the boundaries: the elements near the boundaries are 2.175 times the size of those in the centre. In addition, a symmetric boundary constraint is applied on the symmetric plane.
(a)
(b)
Since the experimental results showed that the projectile did not sustain any obvious deformation or mass loss when the striking velocity was less than 1240 m/s, it can be considered to be rigid and defined by the rigid material model MAT_RIGID, whereas the Johnson–Cook model is adopted to define the target. Moreover, using EOS_GRUNEISEN and equivalent plastic strain as the state equation and failure criterion, respectively, when the equivalent plastic strain of the element reaches the critical strain, this element is considered to be damaged and can be deleted. The parameters of the material model are shown in Table 3.
The velocity of the projectile showed an extremely small negative value caused by its few rebound at the end of the penetration. In order to increase efficiency, we set the total time to 0.5 ms because the maximum penetration depth of models all appeared, the velocity of the projectile was 0 and the projectile appeared rebound at this moment, and this velocity (very small negative value) was far less than 2% of the initial velocity.
Figure 10 shows a comparison between the numerical simulation results (the ultimate penetration depth Dmax) in this paper and the test results in the document for different projectile striking velocities . From this comparison, we can see that the mean relative deviation between the calculated results and test results is 0.0265, and the corresponding mean square deviation is 0.0537. Accordingly, this error is less than that for the Forrestal formula calculations. This greater agreement with the test results indicates that the numerical model and the material parameters adopted in this paper are more reasonable and effective.
The results of the numerical simulations based on the material parameters of member 1 are plotted in Figures 11 and 12. From Figure 12(b), we can intuitively see that the stress wave that spreads from the contacting point between the rod and target will arrive at the lateral boundaries of the target in a short time (if we set the total time to 0.5 ms, this time will be t = 34 µs for a target with rt/a = 10). However, then, the stress wave will reflect and spread back and forth constantly between the projectile and the target. In the process of the projectile contacting the target, entering the target, and finally stopping, the highest stress occurs at the point of contact, which is the reason that the material between the projectile and target is easily destroyed.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
Figure 13 shows the time-history curves of the penetration depth D and projectile velocity . It is quite evident that the projectile velocity decreases continuously over time, during which the penetration depth increases correspondingly. This increase in penetration depth is faster in the early stage before stabilizing.
Hereafter, to get the approximate influence curve of the penetration caused by the target size while taking computer’s computational efficiency into consideration, the size of the target is changed (, 5, 10, 15, 20, 35.7, and 60) while keeping the other parameters consistent with those of member 2 to investigate the influence of different target radius-to-projectile radius ratios on the ultimate penetration depth Dmax, as shown in Figure 14, and the relevant detailed influence extents are analysed in part 4.2. From the results in this figure, we can intuitively observe that as the ratio decreases, the penetration depth increases, whereas the penetration depth increases only slightly when is high (the penetration depth when is only 5.01% higher than that when ). At this time, the impact of the target size on the ultimate penetration depth is small, even negligible. However, when is low, especially when , the penetration depth suddenly increases greatly as decreases (the penetration depth when is 27.54% higher than that when ), and the impact of the target size on the ultimate penetration depth is increased, which is basically consistent with the theoretical calculation results.
4.2. Analysis of the Lateral Boundary Influence and Its Eliminating Condition
The theoretical and numerical results both show that when is small, there are substantial differences in the penetration performance between a rigid projectile penetrating into a thick, finite-radius, metallic target and that of the same projectile penetrating into an unlimited-size target. To analyse this difference further and quantify it, a large number of simulations of penetration under different working conditions when and penetration of a semi-infinite target with the corresponding impact velocities are performed. Adding initial velocities in large scale to each model with different target radius-to-projectile radius ratios and comparing the penetration depth (D) results of the finite-radius aluminium-alloy targets with those of the semi-infinite aluminium-alloy targets, regular results based on sufficient simulation data are obtained, as shown in Figures 15–24. In these figures, the deviation value of the penetration depth is defined as the percentage of the difference between the penetration depth of the semi-infinite target and the penetration depth of the finite-radius target.
The penetration depth comparisons in Figures 15 and 16 indicate that the free lateral boundaries of the target have a substantial influence on the penetration damage when and that the deviation in the penetration under these two boundary conditions can still reach 16.8% (far exceeding the given threshold of 5%), although the projectile penetrates the target with a low initial velocity ( m/s). Thus, it can be considered that the effect of the free lateral boundaries of the target should not be ignored even if the striking velocity is low. Figures 17 and 18 show that, for , the deviation in the penetration under these two boundary conditions is small when m/s; hence, the effect of the free lateral boundaries of the target on the penetration damage is slight (≤5%). However, this effect increased significantly when m/s, for which the deviation in the penetration is 9.2%. Therefore, it can be considered that, for , the effect of the free lateral boundaries of the target can be ignored when m/s. Similarly, for , as shown in Figures 19 and 20, although the effect of the free lateral boundaries of the target on the penetration damage is slight (≤5%) when m/s, the deviation in the penetration can reach 20.5% when m/s; thus, for , the effect of the free lateral boundaries of the target can be ignored when m/s. Similarly, the plots in Figure 21 show that, for , the effect of the target lateral free boundary on the penetration can be ignored when m/s (the deviation in the penetration under these two boundary conditions is 5.2%, which begins to slightly exceed the given threshold of 5%). However, it can be seen from Figures 22–24 that, for , the effect of the free lateral boundaries of the target on the penetration damage is slight because the deviation in the penetration is small when m/s; this deviation is 2.1% when and 2.6% when (). Consequently, the effect of the free lateral boundaries of the target on the penetration damage is considered to become insensitive when .
As the striking velocity is generally less than two Mach numbers in weapon site tests, combining the theoretical results and the numerical simulation results, for thick, finite-radius, metallic targets, the effect of the free lateral boundaries of the target on the penetration performance can be ignored by conservatively designing the target radius to be greater than 28 times the projectile radius.
5. Conclusions
(1)Considering the effect of the free lateral boundaries of the target, the intermediate principal stress, and the different strength criteria, the finite cylindrical cavity expansion model based on the unified strength theory was built in this paper. Then, penetration resistance formulas and penetration depth formulas for rigid projectiles penetrating into thick, finite-radius, metallic targets were deduced, which can be applied to resolve the complicated stress problems of all kinds of targets with equal tensile strength and compressive strength and different nose shapes. Additionally, these formulas are suitable for penetration problems of semi-infinite metallic targets. According to a comparative analysis, the calculation results in this paper are in better agreement with the experimental results (i.e., have higher precision) than the calculation results in relevant documents.(2)The calculation method in this paper can obtain a series of analytical solutions based on different strength criteria to predict penetration depth ranges of all kinds of targets with equal tensile strength and compressive strength. The results in [14] are only a special case of the results deduced in this paper when . The strength parameter b had an obvious influence on the penetration results. It is considering the intermediate principal stress effect that can objectively show the strength potential of the material and ensure that the antipenetration performance of the member is fully utilized.(3)When , the finite cylindrical cavity expansion model fully undergoes the elastoplastic stage and the plastic stage, whereas when , it only undergoes the elastoplastic stage. As decreases, increases, and the extent of this increase is obvious when . The results showed that the penetration depth was obviously affected by the target boundary size under these conditions, and the calculations cannot be performed under the assumption of an unlimited-size target.(4)According to the combined theoretical analysis and numerical simulation results, if the target radius is designed to be greater than 28 times the projectile radius, the effect of the free lateral boundaries of the target on the penetration performance can be ignored in weapons field tests. Furthermore, the results showed that the striking velocity of the projectile and the shape of the projectile nose also greatly influenced the penetration depth.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This study was financially supported by the National Natural Science Foundation of China (Grant no. 51878056), the Social Development Foundation for Science and Technology Planning Project of Shaanxi Province (Grant no. 2019SF-256), and the Fundamental Research Funds for the Central Universities, CHD (Grant nos. 300102289105 and 300102280203).