Open access
Technical Papers
Nov 10, 2023

Simulation Study of the Damage Effect of Tandem Warhead on Concrete Target

Publication: ASCE-ASME Journal of Risk and Uncertainty in Engineering Systems, Part A: Civil Engineering
Volume 10, Issue 1

Abstract

The process of a concrete target subjected to a tandem warhead was studied by means of field testing and numerical simulation. First, field testing of a shaped charge impacting the concrete target was carried out, in which the concrete targets included normal strength concrete and ultrahigh performance concrete. Then, the rigid projectile was used to penetrate the damaged concrete targets experimentally. Numerical simulation was adopted to simulate the two stages of a tandem warhead penetrating into the concrete targets. The combined effect of penetration and explosion in the open literature was also reproduced by simulation. By comparing the depth of test data and simulation results, the accuracy of the algorithm, material models, and corresponding parameters was proven. Finally, the numerical simulation method was adopted to simulate the process of a tandem warhead damaging a concrete target, including shaped charge impaction, penetration following the rigid projectile, and explosion inside the concrete. As seen from the comparative results, the parameters used in numerical simulation can appropriately describe the concrete damage, and the tandem warhead can effectively improve the damage degree of concrete targets. The greatest influence of concrete strength on the damage subjected to a tandem warhead is the stage following projectile penetration. The conclusion of this paper can provide important reference significance for the design of tandem warheads and protective engineering.

Introduction

Concrete is widely used in civil and military engineering due to its advantages of cheap materials, easy access, simplicity of manufacturing, and high strength. Compared to normal strength concrete (NSC), the impact resistance, toughness, stability, frost resistance, erosion resistance, permeability resistance, workability, and volumetric stability of ultrahigh performance concrete (UHPC) have been greatly improved. Therefore, UHPC has been widely used in military engineering such as in missile silos, aircraft shelters, weapon depots, military ports, airport runways, air defense buildings, underground command posts, headquarters, etc. In order to improve the weapon’s damage capability, a new type of earth penetration weapon consisting of a precursor shaped charge (SC) and a following kinetic energy (KE) projectile, namely a tandem warhead, has been developed in recent years. When attacking the target, the precursor SC is detonated by the intelligent fuse at an optimal standoff distance (SOD), and a high temperature, high pressure, ultrahigh speed metal is generated to predrill a hole. Subsequently, the following KE projectile advances along the formed hole and continues to penetrate, multiplying the total depth of penetration (DOP). Furthermore, the main charge of following KE projectile explodes subsequently inside the predamaged target to cause greater damage. The tandem warhead can be extensively used for the blockade and destruction of important and valuable targets such as airports, highways, ports, ships, etc.
Most of the existing research focuses on the damage of concrete targets, including the tests of SC into concrete-like targets (Murphy 1983; Wang et al. 2008, 2018; Huerta and Vigil 2006; Murphy et al. 2000, 2003; Esteban et al. 2015; Chen et al. 2023), the theories and simulations of KE projectile penetration into predrilled targets (Murphy 1984; Teland 2001; Bernard 1977; Guo and Wen 2012), as well as the penetration and explosion effects of the projectile in the concrete targets (Sun et al. 2021; Lai et al. 2015; Yang et al. 2018). However, little research has been performed on the full course of destruction of a tandem warhead affecting concrete targets (Heider and Hiermaier 2001; Sohn and Han 2021). For the assessment of concrete target damage, novel methods like acoustic frequency analysis to estimate the materials are worth being explored (Khoshouei et al. 2020; Yari et al. 2019; Yari and Bagherpour 2018a, b).
In this paper, the tests of SC penetrating into NSC and UHPC targets as well as the subsequent rigid projectile penetration into predamaged targets were conducted. The two stages of a tandem warhead attacking concrete targets were reproduced by adopting simulation methods. At the same time, the test in the published literature, including the projectile penetrating into the intact target and explosive exploding inside the predamaged concrete, is also simulated to verify the accuracy of the finite element (FE) model and material parameters and algorithms. On the basis of the verified simulation method, the damage effects of the explosive in the second stage of the tandem warhead exploding inside the predamaged concrete targets is further simulated. Thus the whole process of tandem warhead damage on concrete is studied numerically. Meanwhile, the effect of concrete strength against the tandem warhead is discussed and the damage effects of tandem and KE projectiles on the concrete targets is compared. The conclusions in this article can provide references for the design of protection structures and tandem warheads.

Damage of Concrete Tests

Impact Test of Tandem Warhead

For studying the damage effect of concrete targets subjected to the tandem warhead, an impact test of the tandem warhead was designed and conducted. First, a penetration test of explosively formed projectile (EFP) formed by SC into intact concrete targets was carried out in the field, and then a test of the following rigid projectile penetrating the predamaged concrete targets was carried out indoors.

EFP Impact Test (Hu et al. 2017, 2018)

The SC in the present test consists of the detonator, detonator base, booster, cover plate, explosive, casing, and metal liner (Fig. 1). The thickness of the liner is 1.8 mm. The SC can form a projectile-like metal fragment penetrator with a velocity between 1  km/s and 2  km/s.
Fig. 1. Shaped charge: (a) schematic diagram; and (b) dimensions.
The materials of detonator base, cover plate, and casing are 2A12 aluminum alloy, while the explosive in the casing and booster is Comp B. Fig. 2 is the SC with a liner material of T2 copper, of which the total mass is 62.6 g.
Fig. 2. Test photographs of shaped charge.
Fig. 3 illustrates the schematic diagram and dimensions of NSC and UHPC targets, and the corresponding compressive strengths are 10 MPa and 140 MPa, respectively. The test setup of the SC impact test is illustrated in Fig. 4 in which the standoff distance (SOD) is 70 cm. This test was also performed with a spaced target; see the literature (Hu et al. 2017, 2018) for a detailed test setup and analysis of the results.
Fig. 3. Schematic diagram of concrete targets.
Fig. 4. Test setup.

Test of the Following Rigid Projectile

As shown in Fig. 5(a), the length, diameter, and total mass of the projectiles used in the test are 152 mm, 25.3 mm and 344 g. The projectiles consist of D6A, which has high strength and can be seen as rigid during the penetration process. Fig. 5(b) shows the dimensions of the projectiles with a length/diameter ratio and caliber-radius-head (CRH) of 6 and 3, respectively. The inert polymer material is filled in the projectile for substituting the filled charge and adjusting the centroid of the projectiles. The projectiles are launched by a smooth-bore powder gun, and the striking velocity is designed between 300  m/s and 500  m/s. Fig. 6 illustrates the setup of the following projectile penetration test, and the projectile can accurately penetrate into the borehole.
Fig. 5. Projectile: (a) photograph; and (b) dimensions.
Fig. 6. Projectile impact test schematic diagram. (Image by Feng Hu.)
In total, the number of targets is 12, of which 8 are NSC and 4 are UHPC. As a comparison, the test of two rigid projectiles penetrating into intact NSC targets and the remaining ten into predamaged targets subjected to EFP is designed, and the detailed test conditions are provided in Table 1.
Table 1. Test setup
Concrete targetNo.Test typeStriking velocity
NSCI-8, I-11Projectile300  m/s, 500  m/s
NSCI-1, I-2, I-5, I-6, I-7, I-9Tandem warhead300  m/s, 417  m/s, 500  m/s
UHPCV-1–V-4Tandem warhead300  m/s, 500  m/s

Test Results

The sketch of damaged targets subjected to EFP and the following projectile (combined damage effect) is illustrated in Fig. 7, in which the ΔDOP is the increasing depth caused by the subsequent projectile individually, and Table 2 lists the corresponding test data (Hu et al. 2017, 2018).
Fig. 7. Diagrammatic sketch of damaged targets subjected to (a) EFP; and (b) following rigid projectile.
Table 2. Test data of tandem warhead
No.V0 (m/s)DOP (cm)ΔDOP (cm)
EFPProjectiles
I-129812.020.08.0
I-230511.023.012.0
I-541711.630.118.5
I-641812.533.721.2
I-750013.843.329.5
I-950414.143.829.7
I-830017.217.2
I-1149831.331.3
v-13057.69.92.3
V-23088.110.62.5
V-35008.015.77.7
V-45049.015.36.3

Penetration and Explosion Test (Sun et al. 2021)

In order to study the influence of penetration hole contour on the effect of inside explosion, Sun et al. (2021) carried out a test of explosion after penetration in plain concrete. The test consists of two stages, in which the concrete targets are normally penetrated by the projectile first, and then the explosive was placed at the bottom of the penetration tunnel area for the explosion test after removing the projectile.

Penetration Test

As shown in Fig. 8, the penetration test of projectile [Fig. 8(a)] and the target [Fig. 8(b)] are illustrated. The cylindrical targets with strengths of about 30 MPa and 40 MPa were used, in which the diameters were 0.8 m and thicknesses were 0.5 m (Sun et al. 2021). The shape of the projectile used is ogive-nosed, with a diameter, length, mass, CRH, and material composition of 20 mm, 100 mm, 168 g, 2.5, and 30CrMnSiNi2A steel, respectively. The strength of the projectile nose is about 1.5 GPa, and the range of test velocity is 479.2675.75  m/s.
Fig. 8. Penetration test: (a) projectile; and (b) target.

Internal Explosion Test

After the penetration test, the projectiles were removed from the penetration hole, and the damaged targets were moved outside. The predesigned explosive was loaded into the bottom of the penetration tunnel for the explosion test, see Fig. 9(a) (Sun et al. 2021). It should be noted that the explosive at the penetration bottom is filled with concrete in the fifth shot. The hexogen (RDX) explosive with detonation velocity of 8,750  m/s and detonation energy of 5,215  J/g was adopted. In order to adjust the charge, the explosive was made into a pillar with a 20-mm diameter, 20-mm thickness and 7.85 g, as shown in Fig. 9(b) (Sun et al. 2021). For each shot, the explosive consists of five superimposed pillars, and the total mass is 39.25 g.
Fig. 9. Explosion test: (a) setup; and (b) explosive.

Penetration and Explosion Test Results

The number of explosions after the penetration test is 10, in which the number of 30 MPa targets is six, and that of 40 MPa is four. The velocity in the test is designed between 479.2 and 675.75  m/s, and the total mass of the explosive is 39.25 g. Both the DOP of projectile penetrating and the subsequent explosion are listed in Table 3, where ΔDOP illustrates the net increase in DOP from the explosion test.
Table 3. Test data of penetration and explosion test
No.Target materialsVelocity (m/s)DOP (cm)ΔDOP (cm)
PenetrationExplosive
1C30479.2021.123.12
2C30488.7322.124.62.5
3C30512.3722.625.63
4C30525.2723.226.33.1
5C30551.1326.7
6C30567.4726.229.53.3
7C40548.8222.524.11.6
8C40566.2523.425.52.1
9C40612.5225.329.92.5
10C40675.7527.931.53.6

Numerical Simulations and Comparisons

Finite-Element Model and Materials Parameters

To solve the problem using the Lagrange algorithm can result in severe distortions in large deformation calculations and the Euler algorithm cannot describe the material interface accurately; therefore, Ramaswamy and Kawahara (1987) and Hureta and Liu (1988) developed the multimaterial arbitrary Lagrange-Euler (ALE) algorithm. The ALE algorithm includes advantages of both Lagrange and Euler algorithms. Via the ALE algorithm, the shape and material interface of the continuous substance can be described and the serious distortions are avoided, thereby the high strain rate, high pressure, and nonlinear deformation characteristics can be simulated, accurately. The structured ALE (S-ALE) algorithm was embedded into the commercial FE software LS-DYNA in 2015. The S-ALE solver has exactly the same theoretical basis as the original ALE and uses the same transport and interface reconstruction algorithm. The S-ALE has advantages of simpler mesh generation, less memory required, 20%–40% less computation time, and higher parallel efficiency than the traditional ALE algorithm. In order to numerically reproduce the penetration of the following rigid projectile into the predamaged concrete targets subjected to the EFP, and the explosion in predamaged concrete target by projectiles, the restart input data method, which can divide the complete analysis process into multiple stages and carry over by a restart file, is used.
The damage effects of the tandem warhead to concrete targets includes three steps: precursor EFP forming and impacting concrete targets, following projectile penetrating into the predamaged concrete, and main explosive exploding inside the predrilled target. It includes explosive explosion, energy transfer, blast drive, hypervelocity impact, penetration of KE projectile, internal explosion, and structural damage, etc. In addition, it is a process that includes multiple substances, large deformation, high strain rate, and ultrahigh load.

EFP Impact Simulation

The test setup can be seen as axisymmetric and the 1/4 model is adopted to reduce the time of computer calculation. Fig. 10 shows the FE model and meshes of SC impacting on concrete target (units in cm-g-μs-K). The symmetric boundary condition is applied on the symmetric surface of the model. Meanwhile, the non-reflect boundary is applied around the target for the stress wave not reflecting when it reaches the target boundary, therefore the model excludes the steel culvert. In this simulation, the materials of air, explosives, liners, and casings are described by the multimaterial ALE algorithm, while the concrete materials are computed by the Lagrange algorithm. The interaction between the EFP and targets is calculated using the Fluid-Structure Interaction (FSI) method, which can efficiently solve problems relative to the interaction between the fluid and structures. As shown in Fig. 10, the penetration of EFP into the target is realized through coupling the Lagrange and ALE grids, thus the partial air and target grids are overlapped.
Fig. 10. FE model and grids of EFP impact test.

Simulation of Projectile Penetration

As a continuation of the section “EFP Impact Simulation,” by using the restart input data method which is embedded in LS-DYNA (2001), the penetration test of following projectiles into the predamaged concrete targets is carried out numerically. Similarly, considering that the test setup (Figs. 16) is axisymmetric and reduces the time of computer calculation, the 1/4 model is adopted. Fig. 11 illustrates the FE models and grids of the concrete targets in the tandem warhead impact experiment [Fig. 11(a)]. Fig. 11(b) illustrates the FE models mesh of projectiles into intact concrete in the penetration and explosion test (Sun et al. 2021). The projectiles and concrete targets are computed with the Lagrange algorithm. The keyword *CONTACT_ERODING_SURFACE_TO_SURFACE is used to define the contact behavior between rigid projectiles and concrete targets.
Fig. 11. FE models of projectile penetration into concrete target of (a) tandem warhead impact test; and (b) penetration and explosion test.

Internal Explosion

As a continuation of the section “Simulation of Projectile Penetration” [model in Fig. 11(b)], by adopting the restart input data method, the internal explosion inside concrete target after rigid projectile penetration is simulated in a similar way. Fig. 12(a) shows FE models of internal explosion inside concrete target, and it should be noted that the model consists of air and target, in which the explosive is calculated by S-ALE algorithm. Thus, Fig. 12(b) gives the generated ALE grids of explosive materials, in which the air is not displayed.
Fig. 12. Internal explosion inside concrete of (a) FE models; and (b) generated explosive grids.

Constitutive Models and Materials Parameters

Air
The air material is described by the NULL material model (*MAT_NULL (009#)) and linear polynomial equation of state (EOS), which relates the pressure P with E0 as
P=c0+c1μ+c2μ2+c3μ3+(c4+c5μ+c6μ2)E0
(1)
where c0c6 = equation constants; and E0 = energy per unit volume. The parameters used for air are given in Table 4.
Table 4. Parameters of material model and EOS of air
Density, ρ0 (kg/m3)CoefficientsEnergy per unit volume, E0 (J/m3)
c0c1c2c3c4c5c6
1.29300000.40.402.5×105

Source: Data from Alia and Souli (2006).

Explosives
The material of explosive is described by material model [*MAT_HIGH_EXPLOSIVE_BURN (008#)] and Jones-Wilkens-Lee (JWL) EOS. The JWL EOS can describe the relationship of the pressure, volume, and energy of detonation products as
P=A(1ωR1V)eR1V+B(1ωR2V)eR2V+ωE0V
(2)
where A, B, R1, R2, and ω = materials constants; E0 = detonation energy per unit volume; and V = relative volume. The materials parameters of Comp B, TNT, and RDX are shown in Table 5.
Table 5. Parameters of material model and EOS of explosives
ParametersComp BTNTRDX
ρ0 (kg/m3)1,7171,6301,850
D (m/s)7,9806,9308,300
PCJ (GPa)29.52134.7
A (GPa)52371376.5
B (MPa)7.68×1033.21×1030
R14.24.154.062
R21.10.951.546
ω0.340.30.5
E0 (J/m3)8.5×1094.3×1099×109
V111

Sources: Data from Sun et al. (2021); ANSYS/Autodyn-2D and 3D (2007).

Metallic Materials
The Johnson-Cook (JC) material model (Johnson and Cook 1983) [*MAT_JOHNSON_COOK (015#)] is used to describe the liner and casing materials. The JC model can also describe the strength behavior of the metallic materials subjected to large strains, high temperatures, and high strain rates, which is expressed as
σy=(A+Bεpn)(1+clnε˙*)(1T0*m)
(3)
where A, B, n, c, and m = materials constants; and σy and εp = flow stress and effective plastic strain, respectively. ε˙*=ε˙p/ε˙0 illustrates the normalized effective plastic strain rate and ε˙0=1  s1 illustrates reference plastic strain rate. T0*=(TTroom)/(TmeltTroom) is homologous temperature, and the Tmelt and Troom denote the melting and room temperatures, respectively.
The EOS of Gruneisen is adopted to describe the metallic material under high strain rate (>105  s1), in which the pressure P of materials under compressed conditions is expressed as (LS-DYNA 2001)
P=ρ0C2μ(1+(1γ02)μa2μ2)1(s11)μs2μ2μ+1s3μ3(μ+1)2+(γ0+aμ)E
(4)
While for the expanded condition, the pressure P is given as (LS-DYNA 2001)
P=ρ0C2μ+(γ0+aμ)E
(5)
where C denote the curve intercept of the vs-vp, units in m/s, while the vs and vp indicate the cubic shock velocity and particle velocity, respectively. The s1, s2, and s3 variables denote the dimensionless curve slope coefficients of vs-vp. μ=ρ/ρ01 denotes the volumetric strain, while ρ and ρ0 represent the current and initial densities, respectively. E = internal energy; γ0 = dimensionless Gruneisen parameter; and a = corresponding first-order volume correction. As shown in Table 6, the parameters of SC liner and casing are listed.
Table 6. Parameters of the liners and casing
ParametersT2 copper2A12 aluminum
Density, ρ0 (kg/m3)8,9602,785
Yield strength, A (MPa)90265
Hardening modulus, B (MPa)292426
Hardening coefficient, n0.310.34
Strain rate sensitivity coefficient, c0.0250.015
Thermal softening coefficient, M1.091.0
Melting temperature, Tmelt (K)1,356775
Room temperature, Troom (K)293293
Intercept, C (m/s)3,9405,328
Coefficient, s11.491.34
Coefficient, s20.60
Coefficient, s300
Gruneisen parameter, γ01.992.0
Constant, a0.470
Internal energy, E (J)00

Source: Data from ANSYS/Autodyn-2D and 3D (2007).

Concrete Targets
In 1999, Riedel et al. (1999) proposed the RHT (Riedel-Hiermaiser-Thoma) constitutive model for concrete, which introduced pressure-related elastic limit surfaces, failure surfaces, and residual strength surfaces, as shown in Fig. 13. The variation patterns of the initial yield strength, failure strength, and residual strength of the materials are described. The RHT concrete model (Riedel et al. 1999) can describe the irreversible deformation of certain zones of the concrete materials under impact load, thus the RHT model is adopted to calculate the concrete targets in the tandem warhead impact test and the penetration and explosion test in this paper.
Fig. 13. Schematic diagram of elastic limit surface, residual strength, and failure surface.
The RHT model is divided into three phases, i.e., elastic phase, linear hardening phase, and damage weakening phase, as shown in Fig. 14. The relevant equations involved in the three phases are as follows:
1.
Failure surface equation
The failure surface stress strength of the material is usually determined experimentally in the quasi-static loading mode. In the RHT constitutive model, the normalized pressure is introduced to eliminate the effect of the concrete strength grade, and the failure stress strength σfail is calculated as follows:
σfail=fcσTXC*(ps)R3(θ)Frate(ε˙)
(6)
where fc = quasi-static uniaxial compressive strength; σTXC*(ps) = quasi-static failure surface compression meridian equivalent force strength; R3(θ) = lode angle parameter; Frate(ε˙) = dynamic amplification factor; and ps=p/Frate(ε·) is the quasi-static pressure.
2.
Elastic limit surface equation
The elastic limit surface equivalent force of the material is deduced from the equivalent stress of the failure surface, and the elastic limit stress strength σelastic is calculated as follows:
σelastic=fcσTXC*(ps,el)R3(θ)Frate(ε˙)FelasticFcap
(7)
where Felastic and Fcap represent elastic compression function and the “cap” function, respectively; and ps,el=ps/Felastic is the quasi-static elastic limit pressure.
3.
Damage weakening function
The advantage of the RHT dynamic constitutive model is that it gives the softening equation after damage while defining the material strain hardening effect, and is able to simulate the softening effect of material damage in concrete materials after the equivalent stress strength exceeds the failure stress strength. The damage variable D is defined as the ratio of the cumulative equivalent plastic strain increment to the final failure equivalent plastic strain
0D=Δεpεpfailure1
(8)
where Δεp = increment of equivalent plastic strain; and εpfailure denotes the plastic strain when the concrete fails. The concrete parameters of the RHT model used in this paper are listed in Table 7.
Fig. 14. The three phases sketch map of the constitutive RHT model.
Table 7. Parameters used in the RHT model for concrete
ParametersNSCNSCUHPC
Shear modulus (GPa)4.44726.816.6
Compressive strength (fc)1030/40140
Tensile strength (ft/fc)0.10.10.1
Shear strength (fs/fc)0.180.180.18
Intact failure surface constant A1.61.41.6
Intact failure surface exponent N0.610.610.61
Tens./Comp. meridian ratio (Q)0.68050.68050.6805
Brittle to ductile transition0.01050.01050.0105
G (elastic)/(elastic–plastic)222
Elastic strength/ft0.70.70.7
Elastic strength/fc0.530.530.53
Fractured strength constant B0.70.71.6
Fractured strength exponent M0.80.80.61
Compressive strain rate exponent α0.0320.0320.00909
Tensile strain rate exponent δ0.0360.0360.0125
Max. fracture strength ratio1.0 × 10201.0 × 10201.0 × 1020
Erosion strain111
Reference density (kg/m3)2,7502,7502,750
Porous density (g/cm3)2,3002,4002,500
Porous sound speed (m/s)2,9202,9203,242
Initial compaction pressure (MPa)23.323.39.33
Solid compaction pressure (GPa)5.9995.9996
Compaction exponent333
Bulk modulus, A1 (GPa)35.2735.2735.27
Parameter, A2 (GPa)39.5839.5839.58
Parameter, A3 (GPa)9.049.049.04
Parameter B01.221.221.22
Parameter B11.221.221.22
Parameter, T1 (GPa)35.2710035.27
Parameter, T2 (GPa)000
Reference temperature (K)300300300
Specific heat (J/kg·K)654654654
Thermal conductivity (W/mk)000

Sources: Data from Sun et al. (2021); Hu et al. (2017, 2018); Jeon et al. (2015); Wang and Zhang (2014).

Note: Tens./Comp. = Tensile/Compressive.

Comparison of the Test Data and Simulation Result

EFP Impacts on the Concrete Targets

Fig. 15 shows the progress of T2 copper EFP forming and flying. When the shock wave and gaseous product produced by the detonation of the explosive arrive, the metallic liner is wrapped and accelerated with the high-speed EFP penetrator forming. Taking the NSC target as an example, Fig. 16 further illustrates the process of EFP impacting on the NSC target with compressive strength of 10 MPa. When the penetration depth does not increase further, the penetration process terminates by generating a restart file named “d3dump”, which saves the results of the previous calculation, and can be used for subsequent rigid projectile penetration simulation.
Fig. 15. Progress of T2 copper EFP forming and flying.
Fig. 16. Typical impact process of EFP on concrete.
Fig. 17 gives the simulation results of the EFP impacting the NSC and UHPC targets, in which the concrete is described by RHT model. The dimensions of DOP and spalling crater are illustrated graphically. Furthermore, Fig. 18 shows comparisons of the terminal damage contour map of NSC and UHPC targets, in which the black solid lines represent the average test data and the red dashed lines represent the simulation results of RHT models. This indicates that the simulated results agree well with the test data. Thus, it can be verified that the adopted simulation methods, the constitutive models, and corresponding materials parameters are appropriate.
Fig. 17. Simulated results with targets type of (a) NSC; and (b) UHPC.
Fig. 18. Damage contour map of EFP impact on targets of (a) NSC; and (b) UHPC.

Following Projectile Penetration into Predamaged Concrete Targets

Similarly, taking the NSC target with an impact velocity of 300  m/s for example, Fig. 19 gives the typical process of following projectile penetration into the predamaged concrete target subjected to the EFP (Fig. 16). The projectile penetrates into the predamaged target vertically, while the projectile velocity decreases to zero and a restart file generates the penetration process ends. Fig. 20 further gives damage contours of the target before and after the projectile impact, which represents the target damage subjected to EFP and tandem warhead, respectively. It is seen that the damage of the concrete target induced by EFP [Fig. 20(a)] is recorded and further aggravated by the following rigid projectile impacting [Fig. 20(b)].
Fig. 19. Typical process of following projectile penetration into the predamaged NSC concrete targets: (a) impact instant; (b) penetration process; and (c) penetration termination.
Fig. 20. Damage contours of the NSC targets impacted by (a) EFP; and (b) following rigid projectile.
The simulation results of EFP penetrating into the concrete target were shown and evaluated in the section “EFP Impacts on the Concrete Targets.” Furthermore, Table 8 compares the impact test data and simulation results of the tandem warhead (DOP) and individual following projectile (ΔDOP), in which the average test data is from Table 2. In addition, the simulation results of projectile penetration into the intact NSC targets are listed, and it should be noted that ΔDOP refers to the depth caused by the individual projectile. Considering the largest deviation is just 9.9%, good agreement is derived for the final DOP. Thus, the restart input data method, the adopted constitutive models, and materials parameters are proven to be appropriate.
Table 8. Comparisons of test data and simulation results
Target typeV0 (m/s)DOP (cm)ΔDOP (cm)
ExperimentSimulationDeviation (%)ExperimentSimulationDeviation (%)
I30021.521.5010.09.28.0
I40031.932.82.819.920.53.0
I50043.639.39.929.6278.8
I30017.218.04.7
I50031.330.33.2
V30010.39.93.92.42.68.3
V50015.514.56.57.07.22.9

Rigid Projectile Penetration into Concrete Targets in Penetration and Explosion Test

Fig. 21 illustrates the simulation penetrating process of projectile into the intact concrete target. Similarly, the penetration process terminates with penetration depth not increasing and a restart file generated. Table 9 lists the tolerance of simulation results and test data. It can be seen that the simulation results and the test data are in good agreement for the largest difference of 7.5%, thus the constitutive model and materials parameters, especially C30 and C40 concrete, are verified.
Fig. 21. Typical process of following projectile penetrating into the C30 intact concrete targets with impact velocity of 488  m/s: (a) impact instant; (b) penetration process; and (c) penetration termination.
Table 9. Comparisons of test data and simulation results
TargetV0 (m/s)DOP (cm)
ExperimentSimulationDeviation (%)
C30479.2021.121.41.4
C30488.7322.121.62.3
C30512.3722.623.12.2
C30525.2723.224.34.7
C30551.1326.725.35.2
C30567.4726.226.20
C40548.8222.521.64.0
C40566.2523.423.92.1
C40612.5225.326.33.9
C40675.7527.9307.5

Internal Explosion inside Concrete Targets

Based on the section “Rigid Projectile Penetration into Concrete Targets in Penetration and Explosion Test,” the explosions inside predamaged concrete targets are further simulated by using the restart algorithm. Taking the C30 target with a striking velocity of 488  m/s as an example, Fig. 22 gives the concrete damage caused by the rigid projectile impacting and combined effect of projectile and explosive. Table 10 further compares the simulation results and test data of the combined effect of penetration and explosion (DOP) and individual internal explosion (ΔDOP), in which the test data is from Table 3. It can be seen that the experimental data and numerical simulation results are in very good agreement, because the maximum difference of DOP is only 12%. While the difference of ΔDOP between the experiment and simulation is great, the reason is that the absolute value of the increase is too small. The constitutive model and materials parameters of C30 and C40 are further verified.
Fig. 22. Damage of concrete targets: (a) projectile penetration; and (b) internal explosion.
Table 10. Comparisons of test data and simulation results
Target typeV0 (m/s)DOP (cm)ΔDOP (cm)
ExperimentSimulationDeviationExperimentSimulationDeviation
C30479.2023.124.24.7%22.840%
C30488.7324.624.602.5320%
C30512.3725.627.98.9%34.860%
C30525.2726.328.16.8%3.13.822.5%
C30551.13
C30567.4729.530.01.6%3.33.815.2%
C40548.8224.124.51.7%1.62.981.3%
C40566.2525.527.16.2%2.13.252.4%
C40612.5229.930.93.3%2.54.60
C40675.7531.535.312%3.65.347.2%

Damage Effect of Tandem Warhead on Concrete

In the above sections, the penetration process of tandem warhead into concrete, which consists of the SC impaction on intact concrete target and the following projectile penetration along the hole induced by EFP, was simulated. Meanwhile the concrete damage caused by previous rigid projectile and subsequent internal explosion is simulated. The simulation work was carried out to verify the adopted FE models, constitutive models, algorithms, and corresponding material parameters. In this section, the whole process of tandem warhead damaging the concrete targets, which consists of SC, following, and mean explosive, is further studied by simulation method, i.e., a third stage is added.

Simulation of the Process of a Tandem Warhead Damaging a Concrete Target

Based on the simulation of the section “Numerical Simulations and Comparisons,” the third stage of exploding inside the damaged concrete targets caused by a rigid projectile is added, where the concrete target with compressive strength of 40 MPa is added. Fig. 23 gives concrete target damage of the three different stages of tandem warhead impact, in which the compressive strength is 10 MPa and the impacting velocity is 300  m/s. As a comparison, the projectile with subsequent internal explosion inside concrete is given simultaneously, shown in Fig. 24. In order to be consistent with the test setup in the section “Impact Test of Tandem Warhead,” the explosive used in the third stage is TNT with a diameter of 16 mm and a length of 10 cm, thus the total quality of the main explosive is 32.8  g. Similarly, the explosive is modeled using the S-ALE algorithm.
Fig. 23. Damage of C40 targets subjected to (a) EFP; (b) following projectile; and (c) internal explosion.
Fig. 24. Damage of 10 MPa concrete targets subjected to (a) projectile; and (b) internal explosion.

Comparison of Individual KE Projectile and Tandem Warhead

Fig. 25 shows the total damage effect (depth) of the tandem warhead and individual KE projectile on the concrete target of 10 MPa at different velocities via simulation. As stated in the above sections, the masses of the SC, following projectile, and main explosive are 62.6 g, 344 g, and 32.8 g, thus the weight of SC is just 16.6% of the following KE projectile, while the total DOP increases by 27.6% to 37.2% with a striking velocity between 300  m/s and 500  m/s.
Fig. 25. Damage of tandem warhead and individual KE projectile.

Influence of Concrete Strength against the Tandem Warhead

For studying the effect of concrete strength against the tandem warhead, Fig. 26 gives the simulated concrete target damage in the three different stages, in which the compressive strengths of concrete targets are 10 MPa, 40 MPa, and 140 MPa, respectively. It can be seen that with increasing concrete strength, the damage degree of the tandem warhead on concrete targets is less. The damage depth of tandem warhead on the UHPC target is 34.2%–43.1% and 50.2%–56.4% of those of NSC targets with 10 MPa and 40 MPa; thus, the UHPC plays an obvious role in resisting the tandem warhead.
Fig. 26. Influence of concrete strength against the tandem warhead.

Influence of Concrete Strength on Different Stages of Tandem Warhead Impact

Fig. 27 illustrates the influence of concrete strength on different stages of the tandem warhead impacting on concrete targets. It can be seen that with increasing concrete strength, the damage degree in all three stages decreases. It can be also concluded that the concrete strength influence is greatest in the stage of following projectile penetration, which is the same conclusion as KE projectile penetration.
Fig. 27. Influence of concrete strength on different stages of damage by tandem warhead.

Conclusion

Based on the impact experiment of tandem warhead (a precursor EFP and a subsequent projectile) on concrete targets, as well as the penetration and explosion test, the whole process of a concrete target subjected to a tandem warhead was studied by experiment and simulation methods. The tandem warhead impact test with 12 shots was conducted, in which two were the individual projectiles. Both the tandem warhead impact test as well as the penetration and explosion test are simulated by using the restart input data method implemented in the FE program LS-DYNA. The influences of concrete strength and the following projectile striking velocity on the damage effect by tandem warhead are compared. For the condition discussed in this article, the main conclusions can be derived as follows:
1.
The tandem warhead has prominent impact performance, since the damage depth of the tandem warhead into NSC targets increases, on average, by 27.6% to 37.2% in comparison to the individual KE projectile;
2.
UHPC targets shows much higher impact resistance, considering that the damage depths of the tandem warhead into UHPC targets are 34.2%–43.1% and 50.2%–56.4% of those into NSC targets with strengths of 10 MPa and 40 MPa, respectively;
3.
The greatest influence of concrete strength on the stage of tandem warhead impacting on the concrete targets is the penetration of the following projectile.
The conclusion of this article has important reference significance for the design of tandem warhead and protective engineering.

Data Availability Statement

The experimental data that support the simulation works and the findings of this study are openly available in literature Hu et al. (2017, 2018) and Sun et al. (2021). Some or all data, models, or code generated or used during the study are available in a repository or online in accordance with funder data retention policies.

Acknowledgments

The project was supported by the National Natural Science Foundation of China (12102476).

References

Alia, A., and M. Souli. 2006. “High explosive simulation using multi-material formulations.” Appl. Therm. Eng. 26 (10): 1032–1042. https://doi.org/10.1016/j.applthermaleng.2005.10.018.
ANSYS/Autodyn-2D and 3D. 2007. ANSYS/Autodyn-2D and 3D, version 6.1. User documentation. Canonsburg, PA: ANSYS.
Bernard, R. S. 1977. Empirical analysis of projectile penetration in rock. Vicksburg, MS: US Army Waterways Experiment Station.
Chen, J.-Y., D.-L. Feng, Q.-Y. Sun, C. Peng, Y.-H. Zhu, and S.-Y. Yu. 2023. “Numerical modeling of shaped charge jet penetration into ceramic–metal double-layered medium using smoothed particle hydrodynamics.” Int. J. Impact Eng. 175 (May): 104526. https://doi.org/10.1016/j.ijimpeng.2023.104526.
Esteban, B., L. M. Lenhart, L. Rüdiger, and N. Gebbeken. 2015. “An evaluation of shaped charge experiments using concrete components.” Int. J. Prot. Struct. 6 (3): 439–455. https://doi.org/10.1260/2041-4196.6.3.439.
Guo, X. J., and H. M. Wen. 2012. “Performance analysis and optimization of a dual warhead system.” Int. J. Nonlinear Sci. Numer. Simul. 13 (1): 49–54. https://doi.org/10.1515/ijnsns-2011-100.
Heider, N., and S. Hiermaier. 2001. “Numerical simulation of the performance of tandem warheads.” In Proc., 19th Int. Symp. on Ballistics, 1493–1499. San Antonio: International Ballistics Society.
Hu, F., H. Wu, Q. Fang, and J. C. Liu. 2018. “Impact resistance of concrete targets predamaged by explosively formed projectile (EFP) against rigid projectile.” Int. J. Impact Eng. 122 (Dec): 251–264. https://doi.org/10.1016/j.ijimpeng.2018.08.014.
Hu, F., H. Wu, Q. Fang, J. C. Liu, B. Liang, and X. Z. Kong. 2017. “Impact performance of explosively formed projectile (EFP) into concrete targets.” Int. J. Impact Eng. 109 (Nov): 150–166. https://doi.org/10.1016/j.ijimpeng.2017.06.010.
Huerta, M., and M. G. Vigil. 2006. “Design, analyses, and field test of a 0.7 m conical shaped charge.” Int. J. Impact Eng. 32 (8): 1201–1213. https://doi.org/10.1016/j.ijimpeng.2004.10.002.
Hureta, A., and W. K. Liu. 1988. “Viscous flow with large free surface motion.” Comput. Methods Appl. Mech. Eng. 69 (3): 277–324. https://doi.org/10.1016/0045-7825(88)90044-8.
Jeon, S., T.-H. Kim, and K.-H. You. 2015. “Characteristics of crater formation due to explosives blasting in rock mass.” Geomech. Eng. 9 (3): 329–344. https://doi.org/10.12989/gae.2015.9.3.329.
Johnson, G. R., and W. H. Cook. 1983. “A constitutive model and data for metals subjected to large strains, high strain rates and high temperatures.” In Vol. 21 of Proc., 7th Int. Symp. on Ballistics, 541–547. San Antonio: International Ballistics Society.
Khoshouei, M., R. Bagherpour, M. H. Jalalian, and M. Yari. 2020. “Investigating the acoustic signs of different rock types based on the values of acoustic signal RMS.” Rudarsko-Geološko-Naftni Zbornik 35 (3): 29–38. https://doi.org/10.17794/rgn.2020.3.3.
Lai, J. Z., X. J. Guo, and Y. Y. Zhu. 2015. “Repeated penetration and different depth explosion of ultra-high performance concrete.” Int. J. Impact Eng. 84 (Oct): 1–12. https://doi.org/10.1016/j.ijimpeng.2015.05.006.
LS-DYNA. 2001. LS-DYNA keywords user’s manual. Livermore, CA: Livermore Software Technology Corporation.
Murphy, M. J. 1983. Shaped charge penetration in concrete: A unified approach. Livermore, CA: Lawrence Livermore National Laboratory.
Murphy, M. J. 1984. “Performance analysis of two-stage munitions.” In Proc., 8th Int. Symp. on Ballistics. San Antonio: International Ballistics Society.
Murphy, M. J., D. W. Baum, D. B. Clark, E. M. Mcguire, and S. C. Simonson. 2000. “Numerical simulation of damage and fracture in concrete from shaped charge jets.”. Livermore, CA: Lawrence Livermore National Laboratory.
Murphy, M. J., G. Randers-Pehrson, R. M. Kuklo, T. A. Rambur, L. L. Switzer, and M. A. Summes. 2003. “Experiments and simulations of penetration into granite by an aluminum shaped charge.” J. Phys. IV 110 (Sep): 603–608. https://doi.org/10.1051/jp4:20020759.
Ramaswamy, B., and M. Kawahara. 1987. “Arbitrary Lagrangian–Eulerianc finite element method for unsteady, convective, incompressible viscous free surface fluid flow.” Int. J. Numer. Methods Fluids 7 (10): 1053–1075. https://doi.org/10.1002/fld.1650071005.
Riedel, W., K. Thoma, and S. Hiermaiser. 1999. “Penetration of reinforced concrete by BETA-B-500 numerical analysis using a macroscopic concrete model for hydrocodes.” In Proc., 9th Int. Symp. Interaction of the Effect of Munitions with Structures. Wuhan, China: Scientific Research Publishing.
Sohn, D., and J. Han. 2021. “An empirical approach for penetration of tandem warheads into concrete targets.” Eng. Fail. Anal. 120 (Feb): 105043. https://doi.org/10.1016/j.engfailanal.2020.105043.
Sun, S. Z., H. Lu, S. L. Yue, H. Geng, and Z. Z. Jiang. 2021. “The composite damage effects of explosion after penetration in plain concrete targets.” Int. J. Impact Eng. 153 (4): 103862. https://doi.org/10.1016/j.ijimpeng.2021.103862.
Teland, J. A. 2001. “Cavity expansion theory applied to penetration of targets with pre-drilled cavities.” In Proc., 19th Int. Symp. on Ballistics. San Antonio: International Ballistics Society.
Wang, C., T. B. Ma, and J. G. Ning. 2008. “Experimental investigation of penetration performance of shaped charge into concrete targets.” Acta Mech. Sin. 24 (3): 345–349. https://doi.org/10.1007/s10409-008-0160-3.
Wang, C., W. Xu, and S. C. K. Yuen. 2018. “Penetration of shaped charge into layered and spaced concrete targets.” Int. J. Impact Eng. 112 (Feb): 193–206. https://doi.org/10.1016/j.ijimpeng.2017.10.013.
Wang, G., and S. Zhang. 2014. “Damage prediction of concrete gravity dams subjected to underwater explosion shock loading.” Eng. Fail. Anal. 39 (Apr): 72–91. https://doi.org/10.1016/j.engfailanal.2014.01.018.
Yang, G. D., G. H. Wang, W. B. Lu, P. Yan, M. Chen, and X. X. Wu. 2018. “A SPH-Lagrangian-Eulerian approach for the simulation of concrete gravity dams under combined effects of penetration and explosion.” KSCE J. Civ. Eng. 22 (8): 3085–3101. https://doi.org/10.1007/s12205-017-0610-1.
Yari, M., and R. Bagherpour. 2018a. “Implementing acoustic frequency analysis for development the novel model of determining geomechanical features of igneous rocks using rotary drilling device.” Geotech. Geol. Eng. 36 (Jun): 1805–1816. https://doi.org/10.1007/s10706-017-0433-3.
Yari, M., and R. Bagherpour. 2018b. “Investigating an innovative model for dimensional sedimentary rocks characterization using acoustic frequencies analysis during drilling.” Rudarsko-Geološko-Naftni Zbornik 33 (2): 25. https://doi.org/10.17794/rgn.2018.2.2.
Yari, M., R. Bagherpour, and M. Khoshouei. 2019. “Developing a novel model for predicting geomechanical features of carbonate rocks based on acoustic frequency processing during drilling.” Bull. Eng. Geol. Environ. 78 (Apr): 1747–1759. https://doi.org/10.1007/s10064-017-1197-y.

Information & Authors

Information

Published In

Go to ASCE-ASME Journal of Risk and Uncertainty in Engineering Systems, Part A: Civil Engineering
ASCE-ASME Journal of Risk and Uncertainty in Engineering Systems, Part A: Civil Engineering
Volume 10Issue 1March 2024

History

Received: Nov 7, 2022
Accepted: May 12, 2023
Published online: Nov 10, 2023
Published in print: Mar 1, 2024
Discussion open until: Apr 10, 2024

Authors

Affiliations

Feng Hu, Ph.D. [email protected]
Engineer, Institute of Defense Engineering, Academy of Military Sciences, People’s Liberation Army of China, Beijing 100850, China. Email: [email protected]
Yi Zhang, Ph.D. [email protected]
Senior Engineer, Institute of Defense Engineering, Academy of Military Sciences, People’s Liberation Army of China, Beijing 100850, China. Email: [email protected]
Engineer, Institute of Defense Engineering, Academy of Military Sciences, People’s Liberation Army of China, Beijing 100850, China (corresponding author). Email: [email protected]
Shuang Liu, Ph.D. [email protected]
Engineer, Army 92859 of People’s Liberation Army of China, No. 4 Xiyuan Rd., Hexi District, Tianjin 300060, China. Email: [email protected]
Bukui Zhou, Ph.D. [email protected]
Senior Engineer, Institute of Defense Engineering, Academy of Military Sciences, People’s Liberation Army of China, Beijing 100850, China. Email: [email protected]
Xiao Yu, Ph.D. [email protected]
Engineer, Institute of Defense Engineering, Academy of Military Sciences, People’s Liberation Army of China, Beijing 100850, China. Email: [email protected]

Metrics & Citations

Metrics

Citations

Download citation

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.

View Options

Media

Figures

Other

Tables

Share

Share

Copy the content Link

Share with email

Email a colleague

Share