Numerical Study on Dynamic Response of RC Slabs Subjected to Cyclic Blast Loading

Finite element (FE) model which can be used to investigate accumulated damage of the reinforced concrete (RC) slabs under cyclic blast loading is proposed. Firstly, accuracy of the FE model is veri�ed through an experiment. Secondly, effects of the explosive inclination angle, standoff distance, and the explosive mass on the dynamic response of RC slabs are studied. At last, accumulated damage of the reinforced concrete (RC) slabs under cyclic blast loading is presented using the FE model. It is indicated that parameters of the explosive inclination angle, the scale distance, and the explosive mass have strong in�uences on the dynamic behavior of RC slabs. Under cyclic blast loading with small explosive charge, damage of the RC slab mainly appears in the �rst blast cycle. Displacement of the RC slab, axial force and strain of the rebars increase remarkably, which could be the mainly reason that damage evolution of the RC slab is not obvious, in the second blast cycle.


Introduction
Reinforced concrete is widely used in civil and military buildings because of its good blast resistance and integrity, and its safety under blast load is always one of the most concerned issues.In recent years, there has seen a signi cant rise in explosions caused by terrorist attacks and accidental explosions globally 1- 3 .These explosions have a huge impact on RC structures, resulting in tremendous casualties and economic losses 4,5 .However, due to the randomness of explosion events, it is hard to predict the location and number of the explosions and the blast wave intensity.Thus, it is of great urgency and importance to study dynamic response of RC structures or members under single and cyclic blast loads.
At present, the following methods are primarily adopted to investigate dynamic behaviours of RC structures under blast load: theoretical analysis, eld test and numerical simulation.The more commonly used theoretical analysis mainly includes single-degree-of-freedom (SDOF) methods [6][7][8] , modal approximation 9 and lumped mass model 10 .However, owing to complexity of the structure and its actual boundary conditions, the simpli ed theoretical method cannot fully re ect the actual state of the structure, nor it can be easily used to study the damage characteristics of concrete structural members under blast load.
Field blast test is the most effective method to understand the dynamic behaviour of concrete members subjected to near-eld blast load 11,12 presented details of a large-scale test program on composite structural panels under the action of explosive loading.Luccioni et al. 13 presented behaviour of concrete slab lying on the ground due to blast loads.Ngo et al. 14 investigated the blast-resistance of prestressed super high strength concrete panels subjected to blast load.Ohkubo et al. 15 evaluated the effect of ber reinforcement on improving the performance of explosive-resistant of concrete slabs, experimentally.Lu and Silva 16,17 presented a procedure which could estimate the damage level of RC slabs, by considering the effect of the charge weights and standoff distances.Wang et al. 18 made experimental studies on the dynamic behaviour of small-size RC slabs under close-in blast loads.Yao et al. 19 analyzed through blast experiments the explosion-proof performance and damage behavior of reinforced concrete slabs at various reinforcement ratios.In order to better simulate real in-situ conditions, eld tests of real scale bridge decks were performed by Foglar et al. 20 .Shi et al. 11 examined fragmentation characteristic and local damage of the RC slabs under close-in explosions through full-scale eld tests.Wu et al. 21nvestigated the fragments size distributions of two RC slabs due to airblast loads in a blast chamber.Contact explosion tests were carried out by Li et al. 22 to observe crater and spall damage of ultra-high performance concrete slabs.Failure patterns and dynamic response of the small-scale slabs under contact detonation were obtained by Zhao et al. 23 , experimentally.Ichino et al. 24 investigated effect of parameters such as TSC strength, type and diameter of coarse aggregate on the local failure of two-stage concrete plate under contact explosion.
As a useful supplement to the blast impact test, numerical simulation has been used by many scholars to simulate the blast resistance and dynamic response of structural components under blast load.Among the great deal of numerical analysis software [25][26][27][28][29] , LS-DYNA is the most widely used.Tai et al. 30 studied the dynamic response of RC slabs under the action of the blast wave.Using the ConWep code proposed by Kingery and Bulmash 31 , Yao et al. 19 investigated the damage characteristics and blast resistance of RC slabs with different reinforcement ratios.Strain-rate effects on the compressive and tensile strength of concrete are investigated numerically by Lin et al. 32 .Li and Hao 33 studied spall damage of RC columns subjected to blast loads.Using the ConWep air blast function, Zhao et al. 34,35 made a numerical investigation on the dynamic behaviours of concrete slabs under close-in blast load.Kong et al. 36 studied effect of AFRP on the blast resistance of RC slab.Using Arbitrary Lagrangian-Eulerian (ALE) coupling method, Yang et al. 37 and Feng et al. 38 investigated dynamic response of rubber concrete slabs under close-in blast loading.Wu et al. 39 studied dynamic failure process of RC slab under internal blast load.
In this paper, a numerical analysis is proposed to investigate dynamic behaviours and damage of RC slabs subjected to close-in blast loading, by using the ALE method.FE model is veri ed rstly by simulation of a previous test reported by Zhao et al. 23 .Then, parameters of explosive inclination angle, scale distance, and explosive mass on the dynamic behaviour of RC slabs are studied systematically.Finally, cumulative damage of the RC slab under the cyclic blast load is investigated, and damage evolution of the RC slab during cyclic blast load is quantitatively evaluated.

Previous Tests
Zhao et al. 23 conducted three groups of blasting tests on C30 RC slabs, and studied the effects of reinforcement ratio and explosive mass on the blast resistance of RC slabs.One of the RC slabs measured 100cm×100cm×4cm (length×width×height).A single-mesh layer HRB335 reinforcement bars was applied, with a diameter of 8mm and a perpendicular spacing of 100mm.Dimensions of the specimen are shown in Fig. 1.
The specimen was xed at two opposite sides on the steel frame by U-shape loop clamps.The cylindrical emulsion explosive was hung 40cm over the concrete slab (Fig. 2).The blast-induced vertical acceleration and vertical displacements were measured by 3 acceleration sensors (ASs) and 3 linear voltage differential transducers (LVDTs), respectively.

Numerical Model Calibration
The three-dimensional numerical model of the RC slab in literature Zhao et al. 23 was established using the software of ANSYS/LS-DYNA, and dynamic response of the RC slab subjected to blast load was studied.The results were compared with the test results in the literature to verify rationality of the numerical model.Where Δσ = 3J 2 and J 2 is the second invariant of deviatoric stress; is the hydrostatic pressure; and a y , a y1 , a y2 , a m , a m1 , a m2 , a r1 and a r2 are strength parameters of the material.
With the linear interpolation method, the yield surface of concrete under the current stress state can be obtained and expressed as: Where the value of η is determined by the damage variable λ.When λ increases from 0 to λ m , η increases from 0 to 1, representing the strengthening stage; when λ increases from λ m to in nity, η decreases from 1 to 0, representing the softening stage.
In the LS-DYNA nite element software, this model is coupled with the state equation *EOS_TABULATED_COMPACTION to de ne the functional relationship between the concrete pressure and volumetric strain.The pressure is de ned as: Where the volumetric strain ε v is de ned as the natural logarithm of the relative volume, and are the functions of the volumetric strain, E 1 is the initial internal energy and γ is the heat capacity ratio.
The relationship of concrete pressure and volumetric strain is shown in Fig. 4. It is indicated that while the tensile stress is greater than the threshold value of the hydrostatic pressure, tensile failure occurs.Once the volumetric strain goes over the elastic limit, the concrete is compacted and then becomes a granular material.The material input parameters for the concrete in this paper are shown in Table 1.The dynamic tensile strength of RC changes signi cantly under the action of blast loading, indicating its high strain rate sensitivity.For this reason, the strain rate effect should be considered in the selection of material model for blast-resistance analysis.This is achieved in this paper by adopting the dynamic increase factor (DIF) recommended in the European concrete speci cation CEB 41 .
For concrete in compression: Where f cd is the dynamic compressive strength of concrete under the dynamic strain rate ε; and f cs is the compressive strength of concrete under the static strain rate εs .In this study, εs =30×10 −6 s −1 , logγ 1 =6.156α-2 and α=(5+9f cs /10) −1 .

Reinforcement bar
The reinforcement bar was modelled using the *MAT_PLASTIC_KINEMATIC.This model can simulate isotropic and kinematic hardening plasticity.The isotropic and kinematic hardening can be determined by seamlessly changing the hardening parameter between 0 and 1.The constitutive model of the reinforcement bar is shown in Fig. 5.
The Cowper-Symonds model 32 is adopted to incorporate the strain-rate effect under blast loading.The yield stress after hardening is de ned by Eq. ( 9).
Where σ 0 is the initial yield stress and ε is the strain rate, C s and P s are the strain rate parameters for Cowper-Symonds model, ε eff p is the effective plastic strain, χ is the hardening parameter, E p is the plastic hardening modulus, E t is the tangent modulus, and is the modulus of elasticity.
Parameters for reinforcement material are shown in Table 2.

Table 2
The material parameters for the reinforcement bar 23 Where is the pressure, A 1 , B 1 , R 1 , R 2 and ω are constants, and E e0 are the relative volume and the internal energy per unit volume, respectively.Parameters of the model are shown in Table 3.
Where μ=ρ/ρ 0 − 1 (ρ is the air density and ρ 0 is the initial air density); E a0 is the internal energy per unit volume, and

Erosion criteria
For simulate the failure mode of RC slab after explosion accurately, the erosion algorithm is used to de ne the failure criteria of material through the *MAT_ADD_EROSION keywords 43 .Once the dynamic response of a material element reaches a threshold value, that element will be removed immediately.
Here, the maximum principal strain is selected to de ne the failure criteria of material.The threshold of the maximum principal strain determined by simulation trial is 0.015, which means that the maximum principal strain of the element will be deleted immediately if it reaches or exceeds this value.

Finite element model
The FE model is composed of concrete, reinforcement bars, air and explosive, as illustrated in Fig. 6.In this model, concrete, air and explosive were simulated by SOLID164 element, and reinforcement bars were simulated by BEAM161 element.Reinforcement bars and concrete were modeled separately, and it was assumed that the concrete and reinforcement bars bond well within a very short time without any slip under high blast pressure.The interaction between the reinforcement bars and concrete was simulated by constraining and coupling the reinforcement bars into the concrete using the keyword *CONSTRAINED_BEAM_IN_SOLID.Fixed boundary condition was applied on two sides of the concrete slab and non-re ection boundary condition (NRBC) was applied on the six surfaces of the nite element model to simulate the in nite air domain.

Comparison of numerical and experimental results
The experimental data from specimen P1-2 in literature 23 was selected to verify the applicability of the numerical model established in this paper.In the literature, the standoff distance was 40cm, the mass of the emulsion explosive was 0.2kg, and the scale distance was 0.684 m/kg 1/3 .In this paper, the angle was de ned as θ, which changed clockwise at an interval of 15° between 0° and 90°, as shown in Fig. 7(a)..It can be seen that a longitudinal crack parallel to the xed boundary appears in the middle of the slab, concrete spalling is observed at the four corners of the slab, and no penetration damage across the slab is detected.By comparison, it is found that the calculation results when θ = 90° are consistent with the damage situation of the RC slab tested.
The displacement comparison of the measuring points is shown in Fig. 9.When θ = 0°, the displacement values of points D 1 , D 2 and D 3 obtained by numerical simulation at 4ms are -1.28cm,-0.80cm, -0.95cm, respectively, which are well consistent with the experimental results of -1.15cm, -0.89cm and -1.0cm, and thus lead to reasonable evaluation of the RC slabs response subjected to blast loads.

Results Analysis
With the above numerical model veri ed, the effect of different factors such as standoff distance, explosive mass and cyclic blast load on the dynamic behaviour and cumulative damage of the RC slab is studied, when θ = 60°.

The effect of standoff distance
The dynamic response of the RC slab subjected to blast load was studied when the standoff distance was 40cm, 60cm, 80cm and 100cm, respectively.The explosive mass was 0.2kg.
Figure 10 shows the damage contours of the slab at different standoff distances.With the change of the standoff distance, the slab shows no penetration or spalling damage, and xed boundaries of the top surface are found with serious damage.As the standoff distance increases, the damage zone of the slab gradually decreases and the distribution range of cracks on the bottom surface gets smaller and shrinks toward the center of the slab.The slab shows obvious bending deformation, but the bending degree decreases with the increase of the standoff distance.
Figure 11 shows the displacement vs time history of the RC slab at different standoff distances.As demonstrated in the gure, for the same measuring point, the peak displacement decreases non-linearly with the increase of the standoff distance.Compared with measuring points D 2 and D 3 , the peak vertical displacement of D 1 is the largest at the same standoff distance.The reason might be that D 1 is less affected by the boundary conditions.
Figure 12 shows the axial stress contours of the rebars in the RC slab at t=5ms.Under the action of blast load, stress of the rebars in the middle area of the slab parallel to the xed boundary is larger, because this area is less affected by the boundary in the blasting process.With increase of the standoff distance, bending deformation of the slab decreases, and the number of rebars under tension also decreases.
When the standoff distance grows from 40cm to 100cm, the peak axial stress of the rebars decreases gradually.Especially when the standoff distance increases from 80cm to 100cm, the peak axial stress of the rebars decreases most signi cantly, from 422.3MPa to 234.8MPa.During this period, the detonation wave weakens dramatically.

The effect of explosive mass
The standoff distance is set to 100cm.In this section, dynamic behaviour of the RC slab is studied when the explosive mass is 0.2kg, 0.1kg, 0.05kg and 0.025kg, respectively.
Damage of the RC slab with different explosive masses is shown in Fig. 13.As shown in Fig. 13, due to the large distance of the detonation wave from of the slab surface (the standoff distance-span ratio is 1), damage pattern of the RC slab is characterized by serious damage at the xed boundary of the top surface, and damage of the bottom surface concentrated in the middle area of the slab parallel to the xed boundary direction, accompanied by a small number of cracks.With decrease of the explosive mass, damage zone of the bottom surface decreases, and the bending deformation of the RC slab reduces as well.When the explosive mass is 0.025kg, the damage zone of the bottom surface is approximately bilaterally symmetrical.
Figure 14 shows the relationship between displacement and time history of the measuring points under different explosive mass.As can be seen from the gure, with the explosive mass decreasing, the scale distance increases, energy of the detonation wave and the peak displacement of the same measuring point decrease.Under the same explosive mass, the peak displacement of D 1 is the largest, and those of D 2 and D 3 are close to each other.The smaller the explosive mass, the longer the at section of the displacement-time history curves.
Figure 15 shows the axial stress contours of the rebars at t=5ms.With the decrease of explosive mass, the peak axial stress of the rebars in both the tensile and compressive zones of the slab decreases gradually.Particularly, when the explosive mass decreases from 0.2kg to 0.1kg, the peak value of the axial tensile stress of the rebars changes most signi cantly, from 234.8MPa to 95.4MPa.Under the same explosive mass, the peak axial tensile stress is greater than the peak axial compressive stress of the rebars.Affected by the xed boundary, the maximum axial tensile stress appears in the middle of the rebars perpendicular to the xed boundary direction on the top and bottom sides of the slab.

Dynamic response of RC slab under cyclic blast load
The standoff distance is set to 100cm and the explosive mass 0.025kg, so as to investigate the dynamic response of the slab subjected to cyclic blast load.The cumulative damage of the slab under cyclic blast load is studied by adopting the complete restart technology of the nite element analysis software ANSYS/LS-DYNA.After trial calculations, the calculation time of a single cycle is set to 45ms, which ensures that the RC slab reaches a stable state before the next blast load is applied.

Damage evolution
In the concrete material models selected in this paper, the damage variables de ned are different, depending on whether the concrete is under compression or tension, and the expressions are 44 : , for p≥0 ( 13) , for p<0 (14)   Where is the increment of effective plastic strain, and ε p ij is the plastic strain of the concrete material; f t is the quasi-static tensile strength of the concrete material; r f is the strain enhancement factor; b 1 and b 2 are the parameters that control the softening section of the compressive and tensile stress-strain curves, respectively; and is the hydrostatic pressure.This model describes the damage of concrete through scaled damage factor (SDF), which is expressed as follows:

15
SDF is a positive value of monotonic increase between 0 and 2. 0<SDF<1 shows that the concrete material enters the strengthening section without damage, SDF≥1 shows that the concrete material enters the softening section and begins to damage, and SDF=2 indicates that the material is completely destroyed.
Through the SDF that comes with the model above, the damage and failure of the RC slab can be analyzed through the damage contours.Fig. 16 shows the damage evolution contours of the slab after the rst explosion.After the explosion, obvious damage rst appears on the xed boundaries and bottom surface of the slab.When t=7ms, part of the concrete on the bottom surface undergoes tensile failure and cracks begin to propagate.When t=9ms, longitudinal cracks parallel to the xed boundary direction appear on the slab top surface, accompanied by symmetrically distributed damage zones.As the blast wave continues to release energy, more cracks appear and damage of the RC slab gradually intensi es.When t=15ms, visible cracks appear in the middle of the top surface of the RC slab along the xed boundary direction.When t=45ms, the RC slab tends to stay stable, and the cracks become visible on the bottom surface of the slab.Figure 18 compares the failure elements of the RC slab after two explosions.After the calculation for the rst explosion cycle was completed, some elements at the xed boundaries of the RC slab were deleted due to failure, but the concrete slab was still xed on the support below it.After the calculation for the second cycle was completed, all the elements at the xed boundaries of the RC slab were deleted, which led to failure of the constraint.The RC slab fell as a whole, making the calculation for the third cycle impossible.
In order to dig deeper into the RC slab damage subjected to cyclic blast load, the elements with the same plane coordinates on the surfaces of the slab were selected and their damage evolution process was calculated.The position and number of the selected elements are shown in Fig. 19.
Figure 20(a) shows the damage evolution of the elements on path 1 at the top surface of the slab.During the rst explosion, the damage value of elements H303901, H303921 and H303941 had two sharp increases.The second increase was particularly obvious (at around t=9ms).The damage value of elements H303961 and H303981 had only one increase at around t=9.5ms.After the calculation for the rst cycle was completed, the damage value of element H303901 at the center of the slab was the largest, and that of element H303981 near the boundary was the smallest.During the second explosion, no change was found in the damage value of the elements on path 1 at the top surface of the slab.During the rst explosion, the damage value of elements H311901, H307901, H319901 and H315901 had two sharp rises.The second increase was particularly obvious (at around t=9ms).Element H307901 failed at around t=9.5ms.During the second explosion, the damage value of element H319901 showed a third rise at t=46.5ms, and reached the maximum (SDF=1.48)until t=49ms.After the calculation for the second cycle was completed, the damage values of elements H319901 and H315901 are close, and the damage value of element H311901 is the smallest.
In general, the damage value of the elements at the center of the top surface is larger, while that near the boundary is smaller.The damage value of each element at the top surface increases sharply at around t=9ms.This is because the RC slab rebounds at this time, causing the top surface to be under tension and the damage value of each unit at the top surface to increase.
Figure 21(a) shows the damage evolution of the elements on path 3 at the bottom surface.During the rst explosion, the damage value of element H23901 had two sharp rises.The rst increase was particularly obvious (at around t=2ms).The damage value of elements H23921 and H23941 had one sharp rise and elements H23961 and H23981 showed no damage.After the calculation for the rst cycle was completed, the damage value of element H23901 at the center of the slab was the largest, with its SDF being 1.88.During the second explosion, the damage value of elements H23901 and H23921 demonstrated no change, that of elements H23941 and H23981 had one slight increase and that of element H23961 had two slight increases.Compared with the top surface, the bottom surface of the slab is subject to tension earlier during the explosions, and the elements on the bottom surface are also damaged earlier.
Through the analysis of the damage evolution pattern of the selected elements, it can be seen that the damage of the elements on the concrete slab shows an irreversible increase.The damage value of the elements near the center of the slab is usually larger.Under cyclic blast load, the damage of the bottom surface elements of the concrete slab appears earlier, and the damage evolution pattern of the top surface elements is greatly affected by the rebound of the slab.The damage of the elements in the concrete slab mostly occurs during the rst explosion.The damage value of some elements increases during the second explosion, but with a relatively small increment.

Displacement
The vertical displacement-time history curves for the measuring points of the RC slab during the two explosions are shown in Fig. 22.During the rst explosion, the concrete slab shows an obvious rebound (t=8.5ms-24ms), and then the displacement of the measuring points increases, but at a low rate.After the calculation for the rst explosion cycle is completed, the vertical displacement of the measuring point D1 is 0.79cm, and the vertical displacements of the measuring points D2 and D3 are close to each other, about 0.38cm.The rst blast load leads to local cracking of the slab, which destroys the integrity of the slab and subjects the rebars to yield deformation.After the second explosion, the rebars are subjected to a greater force and gradually enter the strengthening stage.The vertical displacements of the three measuring points on the RC slab increase rapidly.After the calculation for the second explosion cycle is completed, the vertical displacements of the measuring points D1, D2 and D3 are 4.62cm, 4.74cm and 4.70cm, respectively.

Axial force and strain of rebars
In order to analyze the force and deformation of the rebars in the concrete slab during the two explosions, the rebar elements near the center of the concrete slab, the xed boundary and the free boundary were selected for analysis.The numbers and locations of the elements are shown in Fig. 23.
Figure 24 shows the relationships between axial stress and time history of the rebar elements in two explosions.During the two explosions, the axial stress of rebar element B1911 increases gradually, but at a very low rate.The axial stresses of rebar elements B2991 and B3711 follow basically the same change trend.During the rst explosion, the peak axial stress of rebar element B2991 is 267.25MPa, which occurs at t=15.5ms, and that of rebar element B3711 is 294.23MPa, which occurs at t=12.5ms.After the second explosion, the axial stresses of rebar elements B2991 and B3711 increase sharply at around t=47.5ms, and keep going up until they reach the maximum value of 491.26MPa and 474.45Mpa.Compared to the rst blast load, after the calculation for the second blast load is completed, the axial stress of rebar elements B2991 and B3711 increases by 28.47Mpa and 212.78Mpa, respectively.In view of the cumulative damage evolution of the concrete slab, the increase in the axial force of the rebars is one of the important factors that prevent signi cant increase in the cumulative damage of the RC slab after the second explosion.In other words, rebars play a crucial part in reducing the structural damage and failure of RC slabs under cyclic blast load.
Figure 25 shows the axial strain-time history curves of the rebar elements in two explosions.As can be seen from the gure, axial strain of the rebar near the xed boundary of the RC slab shows no evident change, while that of the rebar near the center of the slab and the free boundary increases greatly during the second explosion.Therefore, in case of cracking or overall damage of the RC simply supported slab under the action of cyclic blast load, re-application of the blast load leads to signi cant rise in axial deformation of the rebars, especially in the middle area and free boundaries of the slab.

Conclusions
Dynamic behaviors of RC slabs under blast load were numerically investigated by LS-DYNA in this paper.Validity of the numerical model is calibrated through a blasting test.The effect of the explosive inclination angle, scale distance, and explosive mass on the response of reinforced concrete slabs is explored by intensive numerical simulations.Damage evolution of the slabs under cyclic blast loading is studied using the numerical model.The conclusions are summarized as follows: (1) The parameters such as the explosive inclination angle, scale distance and explosive mass can signi cantly affect the dynamic behavior of RC slabs.When θ takes different values, obvious damage can be found at the xed boundaries of the slab and there is also a strip damage zone perpendicular to the length direction of the strip explosive on the top surface of the slab.
(2) With increase of the standoff distance, the damage zone of the slab gradually decreases and the distribution range of cracks on the bottom surface gets smaller and shrinks toward the center of the slab.
The slab shows obvious bending deformation, and the bending deformation of the RC slab decreases with the standoff distance increases.
(3) Damage zone of both surface of the RC slab decreases with decrease of the explosive mass, and the bending deformation of the RC slab reduces as well.When the explosive mass is 0.025kg, damage zone of the RC slab is approximately bilaterally symmetrical.

3. 1 .
Material models 3.1.1.Concrete *MAT_CONCRETE_DAMAGE_REL3 is used to model concrete40 .This model de nes the plastic properties of concrete through three yield failure surfaces, as shown in Fig.3.The initial yield surface, ultimate strength surface and softened strength surface are expressed as follows:

Figures 8
Figures 8(a)-(g) show the damage of RC slab on the top surface with different values of θ.Under the action of blast load, obvious damage can be found at the xed boundaries of the RC slab when θ takes different values, and there is also a strip damage zone perpendicular to the length direction of the strip

Figure 16 .
Figure 16.Damage evolution of RC slab during the rst explosion

Figure 17 .
Figure 17.Damage evolution of RC slab during the second explosion

Figure 20 (
Figure 20(b) shows the damage evolution of the elements on path 2 at the top surface of the RC slab.During the rst explosion, the damage value of elements H311901, H307901, H319901 and H315901 had

Figure 21 (
Figure21(b) shows the damage evolution of the elements on path 4 at the bottom surface.The damage value of elements H27901, H35901 and H39901 had a sharp increase at around t=2ms, and then the elements failed.The damage value of element H31901 had a sharp increase at around t=1.5ms, with its SDF reaching 1.57, and remained unchanged during the second explosion.

Figures Figure 1 Figure 2 Figure 3 Strength model for concrete material Figure 4
Figures

Figure 9 Comparison
Figure 9

Figure 10 Damage
Figure 10

Figure 13 Damage
Figure 13

Figure 16 Damage
Figure 16

Figure 17 Damage
Figure 17

Figure 23 Selection of rebar elements Figure 24
Figure 23

Table 1
The material parameters for the concrete model. .
42e emulsion explosive is described by the *MAT_HIGH_EXPLOSIVE_BURN material model and JWL state equation42, which is

Table 3
The material parameters for the explosive C 5 , C 6 and C 7 are the state equation parameters.Parameters of the model are shown in Table4.