Research of Steel-Concrete Composite Bridge under Blasting Loads

*is paper presents a study to simulate the performance of steel-plate composite bridge under various blasting loads. *e multiEuler domain method based on the fully coupled Lagrange and Euler models is adopted for the structural analysis of explosion process with the commercial software Autodyn. Due to the difference of material characteristics and space distribution between the concrete and steel part, the most adverse position is estimated to be above and below detonation. A remarkable difference between these two explosive denotations for steel-concrete composite bridge is noted, and the failuremode above denotation is the damage of local concrete deck, with the compressionmode near the denotation point showing a standard trigonometric curve.*e failure mode below denotation includes damage of steel girders and concrete failure near junction.


Introduction
Blasting incidents are becoming an increasing threat for urban infrastructure in the past years.Terrorist attacks and transport vehicles with explosive and flammable materials cause potential blasting loads structures subjected to.Terrorism and accidental explosion have become a major threat to the security of international community.As crucial thoroughfare for urban traffic system, disastrous bridge failures can bring enormous economic and related loss especially in populated regions.Bridge engineers are focused on bridge antiexplosion and blasting damage assessment especially after the incident of 911.However, consideration of prominent explosive resistances in both official and nongovernmental design codes for nonmilitary bridges is rarely found.Steel-concrete composite bridges have been widely adopted in urban construction for its advantages like highly construction speed and cross ability in recent years.However, the research of blast on steel-concrete composite bridge is insufficient.erefore, it is extremely necessary to figure out the stress and failure mechanism of the steelconcrete composite bridge under blasting load.
Many researchers [1][2][3] focused on the structural performance of concrete piers or slabs under blasting loads.Wu et al. [4] investigated the residual capability of localized blast-damaged RC columns by LS-DYNA and the common failure mode of structures subjected to blast load and summarized the residual capacity index by a service axial load index.Yao et al. [5] investigated the damage features and dynamic response of RC beams under blasting loads.Pan et al. [6] investigated three modern types of reinforced concrete bridges under various blasting loads, including a slab-on-girder bridge, a box-girder bridge, and a long-span cable-stayed bridge and discussed the protection performance of carbon fiber-reinforced polymer on bridges.Fujikura and Bruneau [7] studied the blast resistance of steel jacketing bridge piers by experiments.Li and Hao [8] analyzed steel wire mesh-reinforced concrete slab under contact explosion.
Although the researches about bridge or other constructions under explosive listed above have obtained some results and suggestions, the research about the performance of concrete-steel composite bridges under blasting loads is still very limited.In this study, a steel-concrete composite bridge is analyzed thoroughly to obtain the damage patterns of bridge and predict propagation rule of explosive wave accurately.To do this, a three-dimensional finite element model for concrete-steel composite bridge is established by ANSYS Autodyn.Two analysis modules, the traditional preprocessor and AUTODYN module which is professional in explosive area, are utilized in this study..Explosive process is a high-rate chemical phenomenon in which energy is released rapidly within a limited scope, which is obviously different from general structural analysis.Carta and Stochino [9] investigated the flexural failure of reinforced concrete beams under blasting loads through theoretical models, which indicated that the material constitutive relation and conservation equation have a decisive effect on the result of explosive calculation.e materials considered in this study include concrete, steel plate, ideal air, and high-explosive TNT.In order to obtain an accurate and reliable performance of the steel-concrete composite bridge under blasting loads, it is obviously necessary to simulate the behavior of explosion in air with an appropriate method.Generically, during the explosion process (usually in several milliseconds), the nearby air expands rapidly with high energy and temperature, forming a shock wave which evolves into high pressure on structure timely.e spread of shock wave in air can be described by nonviscous flow, deciding by Euler equation as follows:

Conservation Equation
where q is the state vector about time t and f(q), g(q), and h(q) are flux of conservative state variables.
where ρ is the density; u, v, and ω are the velocities in the x, y, and z direction, respectively; p is the pressure; and Q is the total energy.In order to achieve the solution, an ideal state equation is introduced by Zhou et al. [10] as follows: where e is the internal energy, value of which is set to 2.068 × 10 8 to ensure that the structure is at a standard atmosphere pressure at the beginning of the analysis, c is the adiabatic exponent of ideal air, and ρ is the material density.Material parameters of air are shown in Table 1.

Materials Constitutive.
With an explosion taking place nearby, bridges accept a dynamic transient shock wave.erefore, the damage of concrete occurs in a flash.To analyze and simulate the destruction of concrete under the sudden impact during explosion, the concrete failure model of Riedel-Hiermaier-oma is considered for concrete constitutive properties, which include elastic limit surface, failure surface, and residual strength surface.
e constitutive model of RHT is shown in Figure 1.
e equation is as follows: where Y * TXC (p) is the compressing meridian, F rate (_ ε) is the augmentation factor of strain rate, and R 3 (θ) is the corner function.Researchers have shown that there are some defects in the original RHT model in spite of wide application on explosion analysis.Tu and Lu [11] modified the residual strength surface and tension meridian.Nyström and Gylltoft [12] introduced the bilinear principal tensile stress failure criterion in the RHT model.Leppänen [13] modified the tensile stress-strain relationship and tensile strength strain rate of the RHT model.In this study, the principal stress tensile failure was considered instead of original tensile failure model of RHT, value of which was set at 5.0e3.Crack softening option was set, and the value of fracture energy was 100.In explosive analysis, the RHT constitutive relation is usually used with the p-alpha state equation together.e p-alpha state equation is shown as follows: where A 1 , A 2 , A 3 , B 0 , and B 1 are calculation parameters, ρ 0 is the initial density, and e is the internal energy.e parameters are mostly obtained through experiments.In this research, Adobe 1.8-RHT is used in this analysis for its adaptability and widespread applicability from the literature.Material constitutive of steel is much more simple than concrete.Bearing steel 2.7YS is used in explosive analysis.Values of the parameters and other failure mechanism of these two constitutive relations are found in the literature [14].

Advances in Civil Engineering
Jones-Wilkins-Lee (JWL) state equation is used to model to explosive material and ambient air, which can indicate the interaction between explosive pressure and air energy clearly.
e air range surfaces adopt the flow-out boundary condition.e JWL equation is as follows: where p is the pressure during explosion, V is the volume ratio of detonation products to initial explosive, E 0 is the internal energy of unit volume, and A, B, R 1 , R 2 , and ω are material parameters, values of which are shown in Table 2.

Finite Element Modeling.
Figure 2 shows the cross section of a 40 m single span steel-concrete composite bridge whose deck width is 1250 cm.ere are two I-shaped steel girders placed 670 cm apart in the bridge.Concrete bridge decks are adopted for its high pressure resistance and driving comfort.Shear studs are adopted for the connection between I-shaped steel girders and concrete deck, which is considered as ideal connection in this study.e concrete deck is fully restrained with the steel girders.e height of the box girder is 220 cm, while top plate width of I-shaped is 90 cm and bottom width is 110 cm.Transverse beams are set at every 8 meters in longitudinal directions to enhance the structural integrity.In this study, the concrete and steel girders are modeled with 8-node solid Lagrange elements and 4-node shell elements.Because of the complex nonlinearity of the macroscopic concrete, the improved Riedel-Hiermaieroma (RHT) concrete failure model and the porous equation of state (EOS) model are utilized to consider the large-scale heterogeneity and porosity.Due to the complexity of structure, the three-dimensional finite element model was first established by using preprocessor module of traditional ANSYS predisposition.Explosion execution was completed in the analysis module of AUTODYN through "K file." e accuracy of the analysis is highly affected by the cell size, especially for nonlinearity multi-Euler analysis.eoretically, a smaller mesh size of the finite element model will obtain more precise structural performance under blasting loads.However, it is infeasible to model components with unreasonable small mesh size.Pan et al. [6,15] pointed out that a 10 cm mesh size can have an adequate assurance rate to ensure a reliable pressure-time history and dynamic structural response.e range of air is also a key factor for computational efficiency, and it is impossible to simulate an infinite space of air.A more economical approach is used to define the air flow interface.Flow-out boundary in this study is introduced for the air's free edge.Explosive analysis begins with a standard atmosphere which is achieved by defining the initial energy of air. Figure 3 shows the AUTODYN analysis model.

Bridge Performance and Result Discussion
Terrorist attacks and accidental explosion are always unpredictable, for which reason it is necessary to make a large number of detonation scenarios.In this study, considering the particularity of composite structure, detonation locations are the key factor of damage mode.100 kg weight of TNT equivalent is considered for a car bomb.Above-deck and below-deck detonation are simulated in this study.Because of the instantaneity and tremendousness of explosion energy, the effect of gravity and other live load are neglected during such short period.

Above-Deck Detonation.
Dynamite is located in the middle of the bridge, 80 centimeters above the deck.Figure 4 shows the wave propagation of a 100 kg weight TNT by using the three-dimensional multi-Euler domain method.
e blasting wave has spread to whole lateral of the bridge system within only 6.0 milliseconds.e section at the top of the explosive shows the flow-out circumstance.
In order to illustrate the transmission rules of blast pressure on structure clearly, gauges in the air were defined to record the high blast pressure and impulse at different locations during calculation, as shown in Figure 5. Figure 6 shows the pressure history of Gauge #28, where the curve is similar to a triangle.e peak reaches 1.546E6 kPa, with the time corresponding to the peak value is only 1.512E − 1 millisecond.Generally, the area of this triangle can be utilized to describe the degree of explosion.Figure 7 shows the velocity history of Gauge #28, where the air impacts down rapidly towards the bridge deck in a very short time and then is reflected.An opposite velocity can be got in the diagram.Moments later, velocity curve vibrates around the zero line.
rough the calculation results of monitoring points, the overall trends are that a smaller the pressure will be obtained with a greater distance.Figure 8 shows the pressure-time history of some gauges.Meanwhile, Gauge #16 seams have  Advances in Civil Engineering a larger pressure peak than Gauge #17, and it is due to the negative pressure phenomenon which occurs in a small distance from detonation. Figure 9 summarizes the ultimate damage states in the case of above-deck detonation.It can be seen clearly that due to the full con nement e ect and the highly re ected pressure, a 4.6 m × 4.6 m damage area was generated.e surface of the deck is predominantly damaged by direct crushing of RC elements, and at the same time, there is no damage observed in the I-shaped steel girder.Figures 10 and  11 present stress distribution.It can be seen that the blue area in the middle of the bridge re ects the quitting of concrete under blasting loads.

Below-Deck Detonation.
Another situation for the urban bridge is carried out.Overpass bridges can su er a belowdeck detonation because of heavy volume of tra c below them.In this section, explosive detonation is arranged right down the composite bridge.Figure 12 presents a below-deck explosion process: explosive wave forms within a rectangle space between the I-shaped steel girders in the process of transmission.
Figure 13 shows the structure deformation under blasting loads.A large local deformation occurs in the steel girder, value of which reaches 158 mm.Because of the transverse beam, the max deformation is not in the midspan.Figure 14 shows the stress distribution on damage state.Structural failure occurs mainly in two positions: the rst one is the local failure of steel girder itself and the other is the juncture position of steel girder and concrete deck.
In a similar way, Gauges are de ned to gure out the stress characteristics and damage mechanism on structure.Gauges 113, 114, and 50 to 54 are arranged on the longitudinal direction and Gauges 54, 75, and 96 are arranged on the transverse direction, which is shown in Figure 15. Figure 16 shows velocity-time history of Gauges on z direction.e curve of #53 is even higher than Gauge #54 on local scale.It is because that the transverse beam right near the point limits the local deformation.Other curves of Gauge points decreased with distance.Figure 17 shows velocity-time history of Gauge points of the horizontal setting.e velocity of Gauge #75 presents some irregularity, because of the I-shaped steel beam and the local re ection e ect of the explosion wave.e pressuretime history of Gauges #113 and #114 (horizontal near the detonation point) illustrates the same phenomenon.eir pressure curves cannot be simply concluded as typical triangle, and there are two peaks which are shown in Figures 18 and 19.Apparently, both failure characteristics and blasting loads below the detonation are very di erent from those above the detonation.Advances in Civil Engineering bending moments in the instant of the explosion.When explosive happens below the deck, the overpressure characteristic of the structure is more complex.Shock wave re ections are produced in local space.e overpressure time curve near the blast point showed multiple peaks, and the I-shaped girder deformed considerably.

Conclusions
is paper presents the blast analysis results of concrete-steel composite bridge by an accurate simulation of explosion process.e conclusions are as follows.
In the case of a steel-concrete composite bridge, a 100 kg TNT above-deck explosion right on the midspan of the bridge is identi ed as the most critical case.e force of explosion lead to a 4.6 × 4.6 m damage area signi cantly.e pressure history curve of the point right below the explosive can be described as a typical triangle because of the hollowness above the bridge deck.It means that a simple triangle time load can be utilized in explosive analysis instead of complex simulation of the whole blasting system.Arranging the detonation point at the center, it can be shown that a smaller pressure peak will be obtained with a farther distance.
ere is not a remarkable disruption on steel I-shaped beam.
For below-deck detonation of composite bridge, the blast propagation is more complex than that of the above-deck detonation.Blasting loads lead to more deformations and stresses for steel girder.e damage pattern mainly occurs e pressure-time history curve is not a simple triangle, which means a simpli ed time load cannot be utilized to replace the actual explosive process.e transverse beam plays an important role in displacement limitation, but the overall destruction of the structure still lays on the local steel yield strength and the joint failure.
Data Availability e data used in the manuscript can be replicated through calculation method as described in the manuscript.e data supporting the conclusions of the study can be obtained in the manuscript.

Figure 1 :
Figure 1: e RHT constitutive model used for concrete.

Table 1 :
Material parameters of air.
x y z Figure 13: Deformation su ering from explosive load.