Abstract

A calculation method of SCS wall which is used in the third generation of nuclear power plants to resist perforation from rigid projectile based on energy method is proposed in this paper. The energy is divided into four parts including the energy dissipated by front steel plate, concrete, back steel plate, and tie bars. The method accounts for the perforation of the concrete and steel plates separately and accounts for the interaction between them, and a practical antiperforation calculation formula of SCS wall with tie bars is given. The most formular results are close to the test results and the FEM results with a deviation less than 10%, which shows that the calculation formula given in this paper is reasonable and credible to effectively evaluate the perforation failure of the SCS wall and carry out a relevant design. The energy dissipated by the steel plate is much larger than that of the tie bars through a comparative analysis of dissipated energy. The effects of various factors on perforation velocity are analyzed according to finite element calculation results, which can be roughly divided into three categories: the influence of the thickness of steel plate and distance of tie bar is the largest effect, followed by that of yield strength of steel plate, yield strength of tie bar and diameter of tie bar, and that of compressive strength of concrete is the smallest effect.

1. Introduction

Steel-concrete-steel (SCS) sandwich wall is constructed with steel surface plates that act as concrete reinforcement. The tie bars and shear studs are welded to the steel faceplates to develop the composite behaviour of the steel faceplates and concrete. The SCS walls are used in the shield building of the third generation of nuclear power plants to keep the internal steel containment and reactor cooling system from damage by external events. It is important, therefore, that researchers need to study local damage caused by aircraft engine or missile impact, for example, perforation analysis is one of the important research contents. And engineers need an accurate and convenient calculation method to design SCS walls for resisting aircraft engine or missile impact.

Due to the demands of military and nuclear industries, reinforced concrete (RC) targets against perforation of missile have been widely studied with valuable test data and relevant empirical formulas by researchers, such as Petry formula, ACE formula, CEA-EDF formula, UKAEA formula, NDRC formula, and Kar formula. Kennedy [1] made a detailed review of empirical formulas of thick concrete targets against perforation of missile impact, which basically covered test data of European and American countries before the 1970s of the last century. Bangash [2] summarized empirical formulas of RC perforation in different periods, Ranjan et al. [3], Williams [4], and Wen [5] compared results of different empirical formulas for RC perforation with numerical analyses and test data, and Li et al. [6] discussed empirical formulas for RC perforation recently published in detail. These empirical formulas provide a direct and convenient suggestion for the perforation prevention design of RC structures, but the empirical formulas also have their own limitations, which cannot make perforation prediction results to be in good agreement with the test data under all conditions.

Song et al. [7] and Feng et al. [8] established a design and calculation method of antiperforation of RC containment and performed parameter analyses, which considered the effects of various reinforcing bars. Based on the cavity expansion theory, Li et al. [9, 10] and Chen et al. [11] gave the formulas for calculating the minimum concrete thickness and energy of a rigid projectile perforating a plain concrete wall. The calculated results of the formulas are in good agreement with the experimental results, and a perforation calculation method of RC wall against rigid projectile was given.

SCS wall is a new type of structure compared with RC wall, and the research and calculation method of perforation of SCS wall is not widely available. Grisaro and Dancygier [12] gave an antiperforation formula for calculating the minimum thickness of RC wall with a rear steel plate, which was obtained by combining the existing formulas for calculating the perforations of RC and steel plate against rigid projectile, respectively, and the result from the given formula was in good agreement with the experimental data. Siddiqui et al. [13] conducted a series of perforation tests of RC shielded steel plates against rigid projectiles; in such structures, the RC plate and steel plate were separated and RC plate surrounded the inner steel plate. The damage characteristics of the inner steel plate against the impact of ogive and biconical shaped projectiles were studied. Finally, quantitatively evaluate how safe the steel plate is under the impact of these projectiles by using the available perforation formulas of the RC plate and steel plate. Bruhl et al. [14] presented a calculation method for the perforation design of SCS wall subjected to rigid projectile impact without considering the effect of tie bars and the front steel plate.

A calculation method of SCS wall to resist perforation from rigid projectile based on energy method is proposed in this paper, which accounts for the perforation of the concrete and steel plates separately and accounts for the interaction between them, and a practical antiperforation calculation formula of SCS wall with tie bars is given. The rationality and feasibility of the given calculation formula are verified by comparing its calculation results with existing test data and finite element analysis results, and the effect of various factors on perforation velocity is analyzed by the finite element method (FEM).

2. Establishment of Calculation Method of SCS Wall to Resist Perforation

This calculation method is based on the energy method. When a projectile normally impacts the SCS wall with tie bars under the condition that the translational constraint around the SCS wall is performed along the impact direction of the projectile, then the perforation energy dissipation can be divided into four parts: (I) energy dissipation of the front steel plate, (II) energy dissipation of the concrete (III) energy dissipation of the rear steel plate, and (IV) energy dissipation of the tie bars. The residual velocity is zero after perforation in case of the critical condition; then the following equation as is true according to the conservation of energy:where Ek is the initial kinetic energy of the projectile, Ws1 is the energy dissipation of the front steel plate, Wcon is the energy dissipation of the concrete, Ws2 is the energy dissipation of the rear steel plate, Wt is the energy dissipation of the tie bars. The local failure mechanism associated with projectile impact on SCS wall is that (I) the projectile is assumed to penetrate the front steel plate and dislodge a steel disk, (II) the original projectile and the steel disk become a projectile impacting the concrete and dislodge a circular truncated cone of concrete with cracked tie bars, (III) the original projectile, the steel disc and the concrete plug with cracked tie bars impact the rear steel plate. As shown in Figure 1, θ is the angle between the generatrix of the circular truncated cone and the impact direction.

2.1. The Energy Dissipation of the Front Steel Plate

Chen et al. [15] used cavity expansion theory to propose a formula of perforation velocity of the rigid projectile impacting metal plate. The calculation results of this formula are in good agreement with the relevant experimental results. The formula for calculating the energy dissipated by the front steel plate is shown in equation (2), and the perforation velocity is shown in equation (3) by Chen et al. [15]:where M1 is the mass of projectile, d is the diameter of projectile, Vbls1 is the perforation velocity, σ is the yield strength of steel plate, ρs is the density of steel plate, Hs is the depth of steel plate, A and B are the dimensionless material parameters of steel plate, A is defined as equation (4) by Forrestal and Warren [16], B is equal to 1.1, Es is the elastic modulus of steel plate, n is the strain hardening exponent, N1 and N2 are the geometric parameters of the projectile nose, and N is the geometric parameter of projectile defined by equation (5).

2.2. The Energy Dissipation of the Concrete

Kar [17] proposed a formula as shown in equation (6) for calculating θ of the circular truncated cone when plain concrete against a projectile. In this paper, θ of the circular truncated cone in SCS wall is approximately calculated by this formula, where Hc is the depth of concrete and d is the diameter of projectile. The perforation velocity of the plain concrete is given by Li et al. [6] as shown in equation (7), where Vblc is the perforation velocity of the concrete; fc is the compressive strength of concrete; k, S, and Y are the calculation parameters defined by equations (8)∼(10); and h is the length of projectile nose. The formula for calculating the energy dissipated by the concrete is shown in equation (11):

2.3. The Energy Dissipation of the Rear Steel Plate

The formula for calculating the energy dissipated by the rear steel plate is shown in equation (12), the perforation velocity Vbls2 is shown in equation (13), and N is defined by equation (14), where ρc is the density of concrete:

At this stage, the diameter of the formed projectile contact with the rear steel plate is taken as d to account for the fact that a portion of the impacting concrete surface will crush and thus the entire diameter of the circular truncated cone plug does not impact the steel plate as a rigid projectile. Compared with the total mass of concrete plug, projectile, and disk of the front steel plate, the mass of tie bars cracked in the concrete plug is so small that they can be omitted.

2.4. The Energy Dissipation of the Tie Bars

The tie bars play a similar role as the stirrup in the RC wall for the circular truncated cone plug in the SCS wall. Under the impact load, the bond stress between tie bars and concrete is mainly concentrated in the crack edge area. Outside this area, the tie bars and concrete work together, and the strain of the tie bars in this part of the area is almost zero. The effective deformation length of the tie bar may be defined in equation (15), where le is the effective deformation length of the tie bar, εst is the fracture strain of tie bar, la is the transfer distance of bond stress, which can be defined in equation (16) from a Chinese code [18], where α and ζ are the conversion coefficients of anchorage length, fy1 is the design value of tensile strength of tie bar, ft1 is the design value of tensile strength of concrete, dt is the diameter of the tie bar. The number of tie bars through the side surface of the circular truncated cone plug is defined in equation (17), where s1 and s2 are the distances of the tie bar. The energy dissipation of the tie bars is defined in equation (18), fd is the dynamic factor of the tie bar, σt is the yield strength of the tie bar, (19) can be proposed by substituting equations (15)∼(17) into equation (18):

Under the condition of the SCS wall with tie bars against the rigid projectile, the perforation velocity Vp can be calculated by equation (20); the energy dissipation of each part has been defined in equations (2), (11), (12), and (19) through mechanism analysis:

3. Comparison Analysis between Formular Results and Test Results

Details of SCS walls subjected to projectiles were provided in a previous study [19]. The objectives of the actual impact test were to confirm the dynamic impact characteristics of SCS wall models. The projectile diameter was 150 mm with an overall length of 300 mm. The front end of the projectile had a slight ogive, and the projectile size is shown in Figure 2. The weight of the projectile was approximately 40 kg. The initial impact velocity was accelerated to an impact velocity of approximately 150 m/s for all cases. The overall size of the SCS wall model was 2000 mm × 2000 mm, the thickness of the wall was 175 mm and 250 mm, and the thickness of the steel plate was 6 mm. There was no tie bar in this test wall, and the studs having 13 mm diameter and 80 mm length were installed in 150 mm intervals.

There is no tie bar in the SCS wall in this test, so only the energy dissipation of the front and rear steel plates and concrete is considered. As the projectile nose length is 20 mm, which is much smaller than the projectile diameter of 150 mm, N1 and N2 can be approximately taken as 1.0 from Chen and Li [20]. In addition, the values of other relevant parameters in the above-mentioned calculation formulas are shown in Table 1, and a comparison between the calculation results and the test results is shown in Table 2.

The results calculated by the formula showed good agreement with the test results in terms of failure form and residual velocity. Therefore, the proposed calculation formula for perforation prevention of SCS wall based on the energy method is reasonable and feasible. However, there are only two examples in the above-mentioned comparison, the test examples need to be expanded in order to avoid contingency, and a study on SCS wall against projectile by finite element method is performed in the next section to verify the applicability of the formula.

4. Comparison Analysis between FEM Results and Test Results

4.1. Development of Finite Element Model

Liu and Zheng [21] used concrete constitutive models, MAT072R3, in FEM simulations of RC structures subjected to a GEJ79 aircraft engine. The analytical results showed good agreement with the test results. Liu et al. [22, 23] illustrated that the results correspond well with test [19, 24] and DEM results [25] in general for a projectile impacting a SCS structure in case of MAT072R3. So, MAT072R3 is used to model the concrete of the SCS wall in this paper. MAT072R3 is based on Karagozian & Case (K&C) Concrete Model-Release III, which is a three-invariant model, uses three shear failure surfaces, includes damage and strain-rate effects, and has origins based on the Pseudo-Tensor Model. MAT072R3 can generate model parameters based on inputting the unconfined compressive strength, density, and Poisson’s ratio of the concrete. The strength enhancement factor versus the effective strain rate is given a curve, which was proposed by Malvar and Ross [26]. The input parameters for MAT072R3 in this paper is shown in Table 1. A maximum strain for concrete is set as the element erosion criteria using MAT_ADD_EROSION in LS-DYNA [27].

MAT003 is suited to model isotropic and kinematic hardening plasticity with the option of including rate effects. MAT003 is employed for material models of projectile and steel plates of SCS wall, in this model the strain rate is accounted for using the Cowper-Symonds model [27, 28]. The hardening parameter is equal to zero in this paper; the dynamic yield strength is defined by equation (21). CONTACT_AUTOMATIC_SINGLE_SURFACE is used in the contact between projectile, steel plates, and concrete. Concrete is modelled with solid elements, whereas the steel plates and shear studs are modelled with shell elements and beam elements. CONTACT AUTOMATIC NODES TO SURFACE is used in the contact between projectile and shear studs and used in the contact between steel plates and shear studs. The boundary conditions of SCS wall in the impact direction are fixed:whereyd is the dynamic yield strength, y is the static yield strength, and Q and P are material constants. The parameters of MAT003 input in this paper include density, yield strength, elastic modulus, tangent modulus, Poisson’s ratio, Q, P, and failure strain. The values of the parameters are shown in Table 1.

4.2. Validation of Finite Element Model

The impact test for SCS-175-6T is the test for the SCS wall which has concrete wall thickness of 175 mm and steel plate thickness of 6 mm. The test and FEM simulation results of the SCS wall are shown in Figure 3.

The measured impact velocity of the projectile is 152.4 m/s and the failure mode is a perforation that is consistent with simulation. The steel plate of the front side is clearly cut and the back plate is torn out as shown in Figure 3. The impact test for SCS-250-6T is the test for the SCS wall which has concrete wall thickness of 250 mm and steel plate thickness of 6 mm. The test and numerical simulation results of the SCS wall are shown in Figure 4. The measured impact velocity of the projectile is 147.7 m/s and the failure mode is bulging that is consistent with simulation. The steel plate of the front side is clearly cut, but the rear plate is not torn out as shown in Figure 4. The projectile is completely inserted into the SCS wall and the SCS wall is bulged significantly at the center.

A comparison of the FEM analysis results and the test results in terms of residual velocity and failure mode is shown in Table 2. The measured residual velocity of the projectile after peroration of SCS-175-6T is 33.9 m/s, which corresponds closely with 31.1 m/s in the simulation. For the two cases, the failure mode and residual velocity simulated by FEM correspond well with those from the test as shown in Figures 3 and 4 and Table 2.

4.3. SCS Wall without Tie Bars Subjected to Blunt Projectile

SCS wall without tie bars in this section is divided into four cases, only the thickness of steel plate and SCS wall is changed, and the projectile material becomes rigid material instead of the elastic-plastic material mentioned above, and other parameters and contact types of the finite element model are consistent with those of Section 4.1. These settings have been tested and verified in Section 4.2.

It requires tedious trial work to get perforation velocity by the finite element calculation and analysis. Because the perforation velocity is that the SCS can be perforated, the residual velocity of the projectile is exactly zero. As shown in Figure 5, when the impact velocity is 134 m/s, the perforation of the SCS wall does not occur, but when the impact velocity is 135 m/s, it occurs. It can be seen that the velocity of the projectile at 134 m/s becomes zero and has a small rebound velocity after impact, while the final residual velocity of the projectile is about 9 m/s in the case of impact velocity at 135 m/s. In order to get the perforation velocity quickly and conveniently without losing the rationality of the analysis, the method of determining the perforation velocity by FEM in this paper is that the perforation does not occur at a certain impact velocity V, the impact velocity will increase by 1 m/s, and if perforation can occur, the penetration velocity will be V + 1 m/s. So, the penetration velocity can be accurate to units. As shown in Table 3, a comparison of calculation results of four cases is performed. The calculation results of the formula are close to those obtained by finite element analysis with a deviation of less than 10%. From the comparison results, it can be seen that the calculation formula of antiperforation of SCS wall based on energy method is reasonable and the calculation results are credible.

4.4. SCS Wall with Tie Bars Subjected to Blunt Projectile

The perforation velocity of SCS wall with tie bars instead of shear studs subjected to blunt projectile is proposed in this section. The calculation parameters α is 0.16, ζ is 0.6, the fracture strain εst is 0.3, and fd is taken equal to 1.2, respectively. The Design value of the tensile strength of the tie bar is 300 MPa, and other parameters are consistent with the shear stud. The distance of tie bars is changed in order to discuss and analyze the effect of tie bars on perforation velocity. A comparison of the perforation velocity of SCS-175-4T subjected to blunt projectile is shown in Table 4.

As shown in Table 4, the perforation velocity increases with reducing the distance of tie bars. The results of formula calculation are close to those of finite element analysis with a deviation less than 10%. The energy dissipation of the tie bars corresponding to the distance of 150 mm, 75 mm, and 50 mm are 15232 J, 60927 J, and 137087 J, accounting for 4.51%, 16.23%, and 30.87% of the total dissipated energy of SCS wall perforation, respectively. The relationship between shear stud distance and steel plate thickness is defined in equation (22) from Yang [29]; the distance of shear studs is about 92 mm for SCS-175-4T. The anchorage capacity of tie bars is better than that of studs, and the distance between tie bars is usually larger than that of studs under the same condition. If the distance between tie bars is 75 mm, the energy dissipated by tie bars accounts for 16.23% of the total energy dissipation, while the energy dissipated by steel plates accounts for 62.04% of the total energy dissipation. It indicates that the energy dissipated by the steel plate during perforation is much larger than that dissipated by tie bars in normal design:where Dd is the distance of stud, Hs is the thickness of steel plate, Es is the elastic modulus of steel plate, σ is the yield strength of steel plate.

Some comparisons of perforation velocity for different cases are shown in Tables 57. As shown in the tables, perforation velocity increases with the decrease of tie bar distance, and the perforation velocity calculated by the formula is close to the result obtained by finite element analysis with a deviation less than 10%. After calculation, the energy dissipated by the tie bars also accounts for a small proportion of the total energy dissipation.

4.5. SCS Wall with Tie Bars Subjected to Flat Projectile

The working conditions expanded by the test for SCS wall against blunt projectile in the above-mentioned section because of a blunt projectile in the test. The projectile may also be flat-headed in the design of calculation and analysis, so this section mainly focuses on the flat projectile impacting SCS wall, which changes the distance, diameter, and yield strength of tie bar; the thickness and yield strength of steel plate; and the compressive strength of concrete, to perform a comparison of perforation velocity.

4.5.1. The Distance, Diameter, and Yield Strength of Tie Bar

A study on SCS-250-6T subjected to a flat rigid projectile with a mass of 47.57 kg is performed in this section. The compressive strength of concrete is 50 MPa; the distance, diameter, design value of tensile strength, and yield strength of tie bar are 75 mm, 13 mm, 300 Mpa, and 345 MPa, respectively. Other calculation parameters are the same as Section 4.4. Perforation velocity comparison of SCS-250-6T with different conditions of the tie bar is shown in Table 8.

As shown in Table 8, the formular results are close to the FEM results with a small deviation. The results of formular calculation corresponding to 150 mm, 75 mm, and 50 mm of the tie bar distance are 152 m/s, 160 m/s, and 173 m/s, and the results of finite element analysis are 147 m/s, 167 m/s, and 178 m/s, respectively. Compared with 150 mm, 75 mm, and 50 mm of the tie bar distance is reduced by 50% and about 67%, respectively, and the corresponding perforation velocity calculated by FEM is increased by about 13.61% and 21.09%, respectively. The diameters of 13 mm, 16 mm, and 20 mm of tie bar are increased from 10 mm by 30%, 60%, and 100%, respectively, and the corresponding perforation velocity calculated by FEM is increased by 9.15%, 13.73%, and 15.03%, respectively. The yield strength of 345 MPa and 400 MPa of the tie bar is increased from 307 MPa by 12.38% and 30.29%, and the corresponding perforation velocity calculated by FEM is increased by 2.45% and 4.91%, respectively.

4.5.2. The Thickness and Yield Strength of Steel Plate

The masses of the flat rigid projectile are 47.57 kg and 88.35 kg. The compressive strength of concrete is 50 MPa; the distance, diameter, design value of tensile strength, and yield strength of tie bar are 75 mm, 13 mm, 300 Mpa, and 345 MPa, respectively. Perforation velocity comparison for different conditions of steel plate is shown in Table 9.

As shown in Table 9, the formular results are close to the FEM results with a deviation less than 10% except that the deviation of SCS-250-20T is 10.55%. Compared with 4 mm, 6 mm, and 8 mm of the steel plate thickness is increased by 50% and 100%, respectively, and the corresponding perforation velocity calculated by FEM is increased by about 21.01% and 33.33%, respectively. The thickness of 15 mm and 20 mm of steel plate is increased from 10 mm by 50% and 100%, respectively, and the corresponding perforation velocity calculated by FEM is increased by 12.00% and 32.67%, respectively. The yield strength of 307 MPa, 345 MPa, and 400 MPa of steel plate is increased from 250 MPa by 22.8%, 38%, and 60%, and the corresponding perforation velocity calculated by FEM is increased by 6.37%, 7.64%, and 12.74%, respectively.

4.5.3. The Compressive Strength of Concrete

The masses of the flat rigid projectile is 47.57 kg. The distance, diameter, design value of tensile strength, and yield strength of tie bar is 75 mm, 13 mm, 300 MPa, and 345 MPa, respectively. The thickness and yield strength of the steel plate are 6 mm and 307 MPa. The compressive strength of concrete is from 30 MPa to 60 MPa in 5 MPa increments.

As shown in Table 10, the formular results are close to the FEM results with a deviation less than 10% for different compressive strengths of concrete. The compressive strength of 35 MPa, 40 MPa, 45 MPa, 50 MPa, 55 MPa, 60 MPa of concrete is increased from 30 MPa by 16.67%, 33.33%, 50%, 66.67%, 83.33%, and 100%, respectively, and the corresponding perforation velocity calculated by the formula is increased by 3.27%, 3.92%, 4.58%, 9.15%, 9.80%, and 11.11%, respectively.

4.5.4. Analysis of Effect of Different Factors on Perforation Velocity

In order to simply analyze the effect of different factors on perforation velocity, a net increased ratio index can be defined to evaluate this influence, and a calculation formula of net increased ratio is shown in equation (23), where is the change rate of perforation velocity and rf is the change rate of influence factor:

Generally, the bigger the net increased ratio is, the greater influence of this factor on perforation velocity is. According to equation (23), the calculation results of the net increased ratio of various factors are shown in Table 11. As shown, the influence of factors on perforation velocity can be roughly divided into three categories: the influence of the thickness of steel plate and distance of tie bar is the largest effect, followed by that of yield strength of steel plate, yield strength of tie bar and diameter of tie bar, and that of compressive strength of concrete is the smallest effect.

5. Conclusions and Suggestions

(1)In order to study local damage of SCS walls which is used in the third generation of nuclear power plants caused by aircraft engine or missile impact, a calculation method of SCS wall to resist perforation from rigid projectile based on energy method is proposed in this paper, which accounts for the perforation of the concrete and steel plates separately and accounts for the interaction between them; the energy is divided into four parts, namely, the energy dissipated by front steel plate, concrete, rear steel plate, and tie bars. A practical antiperforation calculation formula of SCS wall with tie bars is given.(2)The rationality and feasibility of the given calculation formula are verified by comparing its calculation results with existing test data and dozens of finite element analysis results. The formular results are close to the FEM results with a deviation less than 10% except for a working condition with a deviation of 10.55%, which shows that the calculation formula given in this paper is reasonable and credible to effectively evaluate perforation failure of SCS wall and carry out a relevant design.(3)In the comparative analysis of dissipated energy, the energy dissipated by the steel plate is much larger than that of the tie bars. The effect of various factors on perforation velocity is analyzed according to finite element calculation results, which can be roughly divided into three categories: the influence of the thickness of steel plate and distance of tie bar is the largest effect, followed by that of yield strength of steel plate, yield strength of tie bar, and diameter of the tie bar, and that of compressive strength of concrete is the smallest effect.(4)Additional research is recommended to further evaluate the influence of concrete thickness on perforation velocity and the adaptability of the antiperforation calculation formula to different projectile noses and higher impact velocities.

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 they have no conflicts of interest.

Acknowledgments

This study was supported by the National Science and Technology Major Project of China (Grant no. 2018ZX06902016).