Blast Loaded Columns—State of the Art Review

: The ever-present threat of terrorist attacks in recent decades gives way to research towards blast-resistant design of structures. Columns, as one of the main load-bearing elements in residential buildings and bridges, are becoming interesting targets in bombing attacks. Research of column blast load behavior leads toward increased safety by identifying shortcomings and problems of those elements and acting accordingly. Field tests and numerical simulations lead to the development of new blast load mitigation technics, either in the design process or as a retroﬁt and strengthening of existing elements. The article provides a state-of-the-art literature review of ﬁled blast load tests and numerical simulations of a bridge and building columns.


Introduction
In the last five decades, terrorist attacks have become more frequent. There are different types of terrorist attacks, but according to data provided by the National Consortium for the Study of Terrorism and Response to Terrorism [1], in the last two decades, explosive attacks exceed 50% of the total number of incidents, shown in Figure 1. The attacks on The Twin Towers of the World Trade Center on 11 September 2001 and bridges in California and New York have an impact on the design of structures in the United States and also in the rest of the world [2]. In every country, the transportation system is essential for performing everyday activities, so the Blue Ribbon Panel (BRP) indicates the transportation system as one of the viable targets for a bombing attack. Due to a large number of bridges worldwide, lots of potential casualties, high repair costs, and importance in everyday life, the bridges are increasingly in the focus of terrorist attacks. This is confirmed by the fact that in the last few attacks in Nigeria in 2020, seven bridges were destroyed. It is important to identify which bridges are vulnerable due to their easy accessibility to protect against attacks. Moreover, BRP states that the columns are one of the most critical components on all types of bridges [3]. As there are many types of bridges and many ways to attack the bridge, it is difficult to predict the construction's response to the blast loadings [4]. When detonation of an explosion is under the bridge, then columns are exposed to large lateral forces, depending on standoff distance, which can result in large deformations leading to flexural or shear failures. The contact explosion can breach the column to render it incapable of supporting the dead loads. For small standoff distance, blast waves can cause a serious reduction in a concrete cross-section in terms of spalling and cratering. Since the column failure depends on the position and amount of explosives, all examined attack scenarios were observed.
The main objective of this review article is an extensive literature overview of experimental and numerical research conducted on blast-loaded columns. Both building and bridge columns are considered due to differences in their static and blast behavior. A systematic summary is given of column behavior, possible damage and failure modes, and a review of software used for blast load simulation and analysis.

Experimental Testing
There are no experiments on real scale specimens of bridge columns due to the experimental setup complexity and high costs, while the building columns are mainly conducted in full scale due to their maximum height of 2.5 to 3 m. Moreover, these types of experiments require special testing ground (usually military field ranges) and trained personnel for handling explosives. Even these special testing grounds have limitations regarding the maximum amount of explosives that can be used in one detonation. This also limits the scale of specimens. Research conducted in recent decades has shown that scale tests provide reliable results and the necessary knowledge to analyze the effects of blast load on full-scale structures [5].
If considering building and bridge columns, except specimen dimensions, there is a difference in their behaviors due to different levels of axial load capacity. Therefore, it is recommended to analyze bridge and building columns separately [4]. A list of conducted experimental research is provided in Table 1.

Bridge Columns
The experiments were carried out on standard RC columns, additionally retrofitted columns, and improved composite concrete columns. The columns are exposed to various scenarios of explosive attacks. In addition to the type of column, the scenarios also differ in the position, type, and amount of explosives.
Williamson et al. [32] provided the list of possible terrorist courses of action and indicated that the hand placed explosives on the column and large truck-bomb below the bridge superstructure can destroy columns and cause bridge collapse.
Due to similarities between the effects of the explosions and earthquakes, Bruneau et al. [6] developed a multi-hazard pier concept that they expect to provide a satisfying level of protection against failure under both loadings. All specimens were concrete-filled circular steel columns (CFCSC) with three different diameters (10.16 cm (4"), 12.7 cm (5"), 15.24 cm (6")) and a minimum steel thickness of 3.2 mm. Specimens were made in 1:4 scale of the prototype bridge columns. Due to security reasons, the actual values of charge weights and standoff distances are not provided. Experimental results showed that even a minimal increase in standoff distance and column diameter significantly reduces column deformation. The CFCS columns showed ductile behavior and high resistance to the effects of the explosion [6]. The same scale and scenario when the explosive was located under the bridge in a car placed near to the column were examined in [7]. They assumed the charge weight similar to the blast weights predicted in FEMA (2003) [33] and FHWA (2003) [34]. Charges are set at heights of 0.25 m and 0.75 m, which correspond to the actual height of 1 m (car bomb) and the half column height, respectively. They concluded that only steel jacketing is not enough to provide adequate resistance to the large shear forces influencing the bottom of the column. Therefore, they found that a better solution is using a fully concrete-filled steel tube (CFST) continuously embedded into the footing. CFST columns provided ductile behavior and sufficient resistance to the lateral forces from earthquakes and explosions. Moreover, the advantage of CFST columns is that they do not have a breach and a spall of concrete, i.e., they do not produce flying debris [7]. Figure 2 shows the connection concept between the foundation beam and the CFST column, which provides the full moment capacity of the column. At a rotation of 3.8 • of the bottom of the column, plastic deformation is visible but without cracking of the concrete. The first cracks occur at a rotation of 8.3 • , while the fracture of the steel tube occurs at 17 • . At the height of the explosive charge, pits and notches appeared on the steel tube, while concrete cracks occurred on the tension side at the bottom and top of the column due to the rigid boundary conditions [7]. They assumed the same blast scenario as Fujikura et al. [7] in their work at the same scale of 1:4, but there were four columns in the test specimen while the bridge prototype has three. Figure 3 shows the test setup for the same blast load scenario but on different types of columns. The RC column exhibited shear failure at the base and cracking of concrete along the column [9], RC SJ shear failure [9], and CFST column flexural failure and buckling [7]. Fujikura and Bruneau [9], in their work, presented experimental and analytical investigations of seismically ductile RC columns and non-ductile RC columns retrofitted with steel jackets. The charge was set to a height of 0.25 m which corresponds to the actual height of 1 m (car bomb), and at this height, the column experienced the maximum deflection. All columns failed in direct shear at the base, but RC columns with steel jackets did not experience any structural damage, and the RC column experienced spalling of concrete at the bottom. Compared to the CFST columns, these columns did not exhibit a ductile behavior. Davis et al. [8] conducted an experiment in two phases. In the first phase, they tested eight small-scale columns where they first changed the standoff distance and then the amount of explosive charge while the scaled distance was kept fixed. In the second phase, they tested 16 columns in half scale (1:2), where 10 of them were set at a small standoff to observe the mode of failure (flexure or shear) like in [4], and the remaining six were to sustain local damage (spall and breach patterns). In all samples, concrete strength, clear cover, concrete class, and reinforcement grade were unchanged. Boundary conditions for tested columns were assumed to be pinned at the top and fixed at the bottom. The test setup did not include axial load because low levels of axial loads provide greater capacity to the column, and without axial load, tests are on the conservative side. Five experimental observations and guidelines for the design of blast-loaded columns are provided in [4,8,35]: 1.
Using protective fences and barriers for vehicles to increase the standoff distance; 2.
The circular cross-section can maintain a lower pressure of up to 1/3 concerning a square cross-section of the same dimensions, so the second guideline is to use circular columns, and also, the pressure reducing factors on the circular column were proposed by Winget et al. [2], Marchand et al. [36], and Fujikura et al. [7], respectively, as 0.80, 0.75, and 0.45; 3.
Increase in the reinforcement in the column, as this increases the shear capacity, ductility, and confinement of the concrete; 4.
Use of continuous reinforcement because discrete hoops can be extracted during a blast load; 5.
Placing longitudinal splices away from the charge if they cannot be completely avoided.
They also proposed three design categories (A, B, C) that depend on the scaled distance and require a different approach to designing, i.e., gravity, seismic, blast. In the C category (Z ≤ 0.6 m/kg 1/3 ), columns are exposed to the higher loads than columns in A (Z > 1.2 m/kg 1/3 ) and B (0.6 < Z ≤ 1.2 m/kg 1/3 ). Williamson et al., in Part II [35], provided a review of experiments represented in Part I [4]. Square columns experienced greater net resultant impulse than circular columns under the same blast loads and also had a larger cross-section area, so less shear occurred at the base. To increase the shear capacity of the column at the base and the ductility, it is necessary to increase the amount of transverse reinforcement. Moreover, columns with continuous spirals had a better performance than columns with discrete ties, which confirms the recommendation given in [8]. Due to changes in column design, the weight of charge, and standoff distance, several levels of damage were obtained [37,38]. The test setup is shown in Figure 4. Superficial damage means that the column performed well and has only surface damage and cracks, while minor damage means spalling of the concrete cover and cracks along the column. Deformations, flexural cracking, and spalling of concrete all together are extensive damage, while the failure of the column means that a shear occurred at the base [4]. Echevarria et al. [18] tested RC and Concrete-Filled Fiber-Reinforced Polymer (FRP) tube (CFFT) bridge columns at a scale of 1:5. CFFT columns are reinforced only with longitudinal reinforcement, while RC columns have spiral hoops in addition to longitudinal reinforcement, shown in Figure 5. The quantity and distance of the explosives were not provided for safety reasons. The columns experienced minimal visual damage, but measurements showed that both concrete and steel experienced large loads and strains. In the residual test, CFFT columns showed greater strength and ductility than RC columns. Fouché et al. [5] made the same prototype of the bridge columns as [9] and exposed them to the same blast load scenario. To improve the behavior of the steel jacked columns to the explosion, they added structural steel collars at the top and bottom of the column. The steel collar at the base of the steel jacked RC column is shown in Figure 2 in [5]. They concluded that the modified column was effective in preventing direct shear and had increased ductility. The obtained maximum base rotation ranges from 8.6 to 10.3 • , and all the columns showed satisfactory behavior.
Yuan et al. [23] experimentally tested two RC bridge columns in scale 1:3 under 1kg TNT contact explosion placed at the height of 33 cm from the ground. As a retaining structure that prevents the rotation and displacement of the column at the top, they used a wall with the opening that was placed at a distance of 1.4 m from the detonation point. Both columns experienced spallation and crushing of concrete cover, but the stirrup fracture was observed only in the square column. The damages to the front and back sides of the circular and square columns are shown in Figure 13. Therefore, the test results showed that the square column had more severe damage than the circular column.
Wang et al. [29] investigated the impact of contact explosions of 1 kg, 2 kg, and 3 kg of TNT on the mode of failure and original and residual axial capacity of Ultra-High Performance Cementitious Composite Filled Steel Tube (UHPCC-FST) bridge column. The columns were made in 1:4 scale and are tested horizontally, and the views of the test setup are shown in Figure 6. The top of the column is pinned, and the bottom of the column is fixed. The cylindrical explosive is placed at a distance of 25 cm from the lower support, which represents the actual position of the explosive in the vehicle at the height of 1 m. Quantities of explosives of 1 and 2 kg made only a crater in the column, while 3 kg fractured the tube and crushed the core, as shown in Figure 7. In the axial compression test, all columns experienced diagonal shear failure.  Based on the reviewed literature of bridge columns, maximum support rotations in field blast tests are shown in Table 2. The information can be used as a preliminary, fast damage assessment of blast-loaded columns based on the measured post-blast column rotations.

Building Columns
In addition to bridges, interesting targets of terrorist attacks are buildings. By damaging the ground floor columns, the whole building loses stability, so in most of the new buildings, the ground floor columns are designed also considering the blast loads. Building columns differ from bridge columns in the magnitude of the axial loads. Moreover, the dimensions of building columns are significantly smaller, so most of the experiments conducted on building columns were in full or half-scale.
Burrell et al. [11] tested two half-scale reinforced concrete columns and six Steel Fiber-Reinforced Concrete columns (SFRC) with steel fiber content from 0 to 1.5% by volume of concrete (non-seismic and seismic detailing) at shock tube. In the experiment, the axial load equal to 30% load capacity was applied using a hydraulic jack. According to their experiments, columns designed seismically (38 mm distance between transverse reinforcement) have smaller maximum displacements and can withstand larger blast loads. Moreover, SFRC columns with non-seismic detailing (75 mm distance between transverse reinforcement) showed smaller maximum displacements and no secondary blast fragments.
Aoude et al. [13] experimentally tested nine Ultra-High-Performance Fiber Reinforced Concrete (UHPFRC) columns designed with Compact Reinforced Composite (CRC). Tests were performed in a shock tube, and a hydraulic jack was used to input the axial load. The applied blast pressures varied from 69 kPa to 689 kPa. Dimensions of the cross-section were 152 mm × 152 mm, and the free span of the column was 1980 mm. The results showed that the UHPFRC columns reduce secondary blast fragments. Increasing the proportion of fibers from 2% to 4% had a positive effect on the decrease in displacement, but a higher proportion of fibers did not result in improvements in blast behavior.
Zhang et al. [12] experimentally blast tested three square columns and one circular column. The columns were made from steel tubes filled with concrete (CFST). Specimens were placed horizontally with a simple boundary condition at both ends. A pneumatic jack was used to input the axial load, shown in Figure 8. The entire length of the specimen was 2.5 m. In the experiment, an emulsion explosive, which has a TNT-equivalent of 0.7, was used. The maximum and residual column deformations are provided for blast loads utilizing 17.5 to 35 kg TNT equivalence at a standoff distance of 1.5 m. The concrete inside the steel tube reduced local deformations, and the energy was dissipated through the global response of the element. Zhang et al. [20] tested two types of CFDST columns, circular and square, with inner and outer tubes, as shown in Figure 9. They concluded that columns filled with normal strength concrete experienced greater crushing of concrete and higher steel buckling than columns filled with UHPC. Moreover, UHPC proved to be very resistant to spalling or crushing. The tested CFDST and CFST samples have similar oscillation periods and displacements, so it is concluded that they behave similarly. Zhang et al. [21] tested six ultra-high-performance concrete-filled double skin tube columns with square hollow sections. At a standoff distance of 1.5 m, the specimen exposed to 35 kg of TNT did not experience any localized damage or steel buckling. Moreover, the axial load (25% of the maximum load) contributed to the reduction in maximum deflection in the middle of the column. They concluded that the ratio of the cavity and the section influence the overall column deflection and period of oscillation; therefore, it is recommended not to go above 0.5. Wang et al. [19] exposed four circular Reactive Powder Concrete Filled Steel Tubular columns to explosion and fire durations of 0, 60, and 105 min. Fixed supports were simulated, and on one side, an axial load was introduced. Steel tube protects columns against cracks and spalling of concrete. After the detonation of the 17.5 kg explosive, the column experienced bending and after 35 kg bending-shearing deformations. With an increase in the number of explosives and with longer exposure to fire, maximum and residual displacements increased. Wang et al. [24] studied four square and four circular CFST columns under close-range blasts. The column length was 2.5 m, and the thickness of the steel tube was 2.8 mm or 3.8 mm. The standoff distance of the emulsion explosive was 1.5 m, and the charge weight expressed through the TNT equivalence was from 25 to 50 kg. Only a 10% increase in the amount of explosive in square columns increases mid-span deflection by 200%, which is assumed to be caused by a large surface exposed to the blast load. Moreover, increasing the thickness of the steel tube by 1 mm (from 2.8 to 3.8 mm) significantly reduced the displacement (by over 50%). Global failure mode was a flexural failure, and after removing the steel, the square columns sustained spalling and crushing of concrete. Circular columns were broken into several parts. Figure 10 shows the damage to the concrete after removing the steel tube after the blast load.  [24] and test pit [21]. Copyright permission obtained from authors. Figure 9. Cross-sections of two CFDST specimens [20]. Copyright permission obtained from authors. Figure 10. Columns after blast test [24]. Copyright permission obtained from authors.
Codina et al. [14] investigated the effects of close-in blast loading on the full-scale reinforced concrete column. They conducted experimental tests and numerical simulations to calibrate the numerical model. The observed column had a square cross-section of 230 mm × 230 mm and a free height of 2.44 m. The column was tested in a horizontal position, and the standoff distance from the center of the charge (8 kg of equivalent TNT-the used explosive is Gelamon VF65, which is equivalent in a mass to 65% TNT) to column 1 was 100 cm and to column 2 was 60 cm. Both columns experienced flexural damage, spallation of the concrete on the bottom side, and crushing of concrete on the exposed side, shown in Figure 11. Figure 11. Damage of the RC columns after detonation at 100 cm and 60 cm standoff [14]. Copyright permission obtained from authors.
Codina et al. [15,16] investigated the behavior of RC columns and RC columns with two types of protection (steel jacket and reinforced polyurethane bricks) exposed to a near field explosion. At a standoff distance of 60 cm, 8 kg of TNT equivalent shaped into a cylinder was placed, resulting in a scaled distance Z = 0.30 m/kg 1/3 . Comparing the test results of three types of columns, the steel-jacked column had the best results in residual capacity and in reducing final deflection. Polyurethane bricks are lighter and have cheaper protection for columns but give three times worse results than steel jacketing. It is recommended to set bricks of higher density to improve the effect in the area of blast load. In [16], they also examined reinforced resin panels with insulation layers as a possible improvement of the column. This protection system gave the best results, minimizing column damage and the greatest deflection reduction. The protection system in the field test is shown in Figure 12. Codina et al. [22], after their research on plain RC columns [14], tested RC columns covered with reinforced resin panels with an insulation layer and steel jacketing. From the obtained results, it can be concluded that a significant reduction in damage and displacement was achieved with the cladding system, but the spalling and burst of concrete cannot be prevented. Figure 14 shows a schematic representation of all specimens and their damage covered by the studies of Codina et al.  [15,16] and (b) resin panels and insulation layer [22]. Copyright permission obtained from authors.  Xu et al. [17] subjected four Ultra-High-Performance Fiber Reinforced Concrete (UH-PFRC) and four High Strength Reinforced Concrete (HSRC) square columns to the effects of blast loads. Columns were tested under varied charge weights from 1.4 kg to 48 kg of emulsion explosive (TNT equivalence factor is 1.4). The standoff distance was fixed at 1.5 m in all tests. Specimens were placed horizontally, and the axial load was applied using a pneumatic jack. The results showed that UHPFRC columns could better withstand overpressure and shock waves, reducing the maximum displacements.
Li et al. [25] tested 10 Ultra-High-Performance Concrete (UHPC; six reinforced with twisted fiber and four reinforced with micro fiber) and 5 High Strength Reinforced Concrete (HSRC) columns. The length of the specimens was 2.5 m, and the cross-section was square with dimensions 0.2 m × 0.2 m. In experiments, the standoff distances of the explosive from the columns were constantly 1.5 m in all tests, but the explosive weight was changed. For the UHPC columns, weights of 17.5, 25, and 35 kg, and for the HSRC 8, 17.5 and 25 kg were used. Residual load capacity tests showed that UHPC columns after blast loads did not lose much on the axial load capacity. UHCP columns also showed much better load capacity after 35 kg TNT detonation than HSRC columns after 8 kg TNT detonation.
Fouché et al. [26] experimentally tested 12 columns at a scale of 1:4 under the blast load. They varied the void ratio, diameters, and thicknesses of the outer and inner steel tubes. The specimens' cross-section generally experienced denting, and that deformation helped to energy absorption from the overpressure from the near-contact explosions. The inner steel tube played the role of a dowel preventing direct shear failure, and this is the advantage of the CFDST columns over CFST and RC columns. On the tensile side of the column, the concrete was crushed, or horizontal flexure cracks appeared.
Dua et al. [27] tested three RC columns (30 cm × 30 cm × 375 cm) in full scale with the same material and geometrical properties on contact explosion at the bottom of the column. They used 0.1 kg of plastic explosive (PEK-TNT equivalent 1.15) and 0.5 kg and 1 kg of TNT. The TNT charge of 1 kg made a hole in the column; 0.5 kg destroyed the concrete cover, and the remaining core has no residual capacity; 115 g TNT equivalent caused the spalling of the concrete cover. Column damage profiles are shown in Figure 15. A contact explosion causes significant local damage on at least three sides of the column, while a far-field explosion causes the worst damage on the front, exposed side. Dua et al. [28] investigated the same blast load scenario on the column concerning the increase in the cross-sectional width of the column. They experimentally tested columns of dimensions 50 cm × 30 cm, 70 cm × 30 cm, and 90 cm × 30 cm. Rectangular columns showed better behavior under contact explosion than squares. They examined the residual load-bearing capacity of the column and determined the column damage index. When the width dimension of the column subjected to the blast load is greater two and more times from depth, the damage index is lower. Kadhom et al. [30] examined five half-scale RC columns. Three columns are strengthened with unidirectional and woven CFRP laminates while the other two remained unprotected. They were first tested in the shock tube on the induced blast load, and thereafter, their residual axial capacity was examined. The column strengthened with CFRP laminates with ±45 • woven fibers showed the best blast behavior and the most ductile response.
Fujikake and Aemlaor [43] investigated how longitudinal and shear reinforcement ratios, concrete strength, and the number of explosives affect the demolition of RC building columns. Their primary research goal is not terrorist attacks but the demolition of dilapidated concrete buildings. They used a Composite 4 (C4) explosive because of its stability and ease of shaping and placed it in the core of the column. They found that shear reinforcement plays a significant role in the residual bearing capacity after blasting; the strength of the concrete also affects the increase in residual compressive and flexural resistance capacities. The quantities of explosives and reinforcement cannot be applied to the external action of the explosion, but certainly, conclusions about the role of reinforcement and the strength of concrete are useful.
Roller et al. [44] observed two scenarios: first is a contact explosion with the amount of PETN that fits in the suitcase, and the second is the close-in scenario when the explosive is in the vehicle. They investigated the impact on RC columns and strengthened columns. Polymer concrete, SIFCON (Slurry Infiltrated Fiber Concrete), DUCON (Ductile Concrete), and Ultra-High-Performance Concrete (UHPC) were used to improve the resistance of the bridge and building columns. The appearances of column damages after contact explosion and the residual load capacities are shown in Table 3. The results showed an increase in residual bearing capacity by up to 70%. Xu et al. [45] tested five columns in an explosion containment vessel (ECV). The columns were exposed to an explosion of 40 g charge mass, and the distance of the explosive was changed in each test from a contact explosion to a standoff distance of 50 cm. They installed four smart aggregates (SAs) in each specimen for internal damage detection, shown in Figure 16. The propagation of the stress wave energy decreases with the formation of cracks under the blast load, and hence the amplitude of the time-domain signal recorded by piezometric smart aggregate sensors decreases with the appearance of cracks. This method of detecting internal damage has proven to be useful for completing the picture of the condition of the structure because internal cracks have a greater impact on the damage index of the structure than surface cracks. Vapper and Lasn [31] examined columns measuring 100 mm × 100 mm × 1000 mm on the action of a different amount of explosion placed at a vertical standoff distance of 300 mm. They tested four types of columns: reinforced concrete columns, columns reinforced with steel fibers, and both types strengthened with Glass Fiber Reinforced Polymer (GFRP).
Plain concrete showed higher compressive strength compared to steel fiber reinforced concrete. GFRP in reinforced concrete columns did not contribute to the increase in the residual strength, while in steel fiber reinforced columns, the contribution was significant. A reduction in surface damage to GFRP-wrapped columns was also observed.

Numerical Modeling
Due to the increasing number of terrorist attacks, new challenges were posed to engineers. It is necessary to have a good understanding of computer programs, their capabilities, and their limitations to predict individual attack scenarios. The use of numerical simulations gives a clearer insight into the blast effects on the entire bridge and individual components. It provides the possibility of determining the most critical parts and problems that cannot be numerically simulated and need to be examined experimentally. Thus, one of the most important and difficult parts of this analysis is to properly define air blast loadings. Table 4 summarizes the software used to predict the effects of explosions on different types of columns.  [52] 2011 RC bridge columns LS-Dyna Crawford [10] 2013 RC building columns + FRP + SJ LS-Dyna Magali et al. [53] 2013 RC building columns Abaqus Eisa [54] 2014 RC building columns Abaqus Abladey and Braimah [55] 2014 RC building columns Autodyn Li and Hao [56] 2014 RC column LS-Dyna Shi and Stewart [57] 2015 RC building column LS-Dyna Liu et al. [58] 2015 RC bridge pier-bent model LS-Dyna, ConWep Cui et al. [59] 2015 RC column LS-Dyna Zhang et al. [12] 2015 CFST building columns LS-Dyna Codina et al. [14] 2016 RC building column Autodyn Zhang et al. [21] 2016 CFDST columns LS-Dyna Arowojolu et al. [60] 2017 RC bridge column LS-Dyna Eamon and Aslendi [61] 2017 RC bridge columns + SFRP LS-Dyna Kravchenko et al. [62] 2017 RC building columns LS-Dyna Kyei and Braimah [63] 2017 RC building columns LS-Dyna Yuan et al. [23] 2017 RC bridge columns LS-Dyna Abedini et al. [64] 2018 RC building columns LS-Dyna Li et al. [65] 2018 CFDST bridge columns LS-Dyna Liu et al. [66] 2018 RC bridge piers LS-Dyna, ConWep Li et al. [67] 2019 CFDST bridge columns LS-Dyna Liu et al. [68] 2019 RC building columns Autodyn, LS-Dyna Liu et al. [69] 2019 RC bridge column + CFRP LS-Dyna Thai et al. [70] 2019 RC column + SJ LS-Dyna Abedini et al. [71] 2019 RC column LS-Dyna Dua et al. [72] 2019 RC columns LS-Dyna Dua et al. [28] 2020 RC columns LS-Dyna Li et al. [73] 2020 CFDST columns LS-Dyna Rajkumar et al. [74] 2020 RC columns LS-Dyna Vavilala et al. [75] 2020 RC building columns + polymeric foam Abaqus Zhang et al. [76] 2020 Segmental CFST column LS-Dyna Yuan et al. [77] 2020 RC column LS-Dyna Yan et al. [78] 2020 RC columns + CFRP LS-Dyna Hu et al. [79] 2021 RC column + CFRP LS-Dyna Conventional Weapon Effects Predictions (ConWep) [80] and BlastX [81] are programs used to calculate the effect of a blast wave from different types of detonation. ConWep is more used for air-blast calculations, including free-field and reflected blast pressure histories from the free-air, surface, and hemispherical burst explosions, and BlastX calculates internal blast pressure histories. BlastX is based on semi-empirical methods, including nonlinear addition laws for blast pressures from multiple reflecting surfaces based on computational fluid dynamics. Second-order Hydrodynamic Automatic Mesh Refinement Code (SHAMRC) [82] is also used to investigate high explosive and blast effects based on finite-difference computational fluid dynamics (CFD) code. Nonlin [83] does not make the empirical adjustments just for blast loads because it is initially designed for earthquake loads. It has similarities with SPAn32 [84] because it performs a nonlinear dynamic response history analysis taking bilinear material properties. Both programs are based on the analysis of a system with a single degree of freedom (SDOF). MSC.Dytran [85] is an explicit finite element analysis (FEA) solution for simulating blast load effects and analyzing the complex nonlinear behavior that structures undergo during detonation. Ansys Autodyn [86], LS-Dyna [87], and Abaqus [88] are programs that provide the ability to simulate detonation, wave propagation from an explosion, interaction with a structure, and nonlinear material behavior, which are known as hydrocode programs specialized for simulations in fluid dynamics.
Ray et al. [46], in their research, compared three methods with three different resolutions for air blast prediction. In the scenario of below-deck detonation, ConWep has a low resolution of air blast prediction, and the charge can be observed as hemispherical or spherical, while BlastX has a medium resolution and considers the shape of the charge and reflections of the blast pressure. The 3D bridge and blast load can be modeled in the SHAMRC because it is an advanced Eulerian-based finite difference code that has high resolution. Research shows that the highest resolution is not always necessary, as it is a mostly low resolution that provides a conservative design. The authors stated that additional analysis is needed to determine the most economical and sufficiently precise tool for a particular problem [46]. The shape of the explosive also drastically affects the resulting pressure and impulse, so it is necessary to use a program that allows the input of charge geometry.
Marchand et al. [36] determined concrete breaching using ConWep and calculated flexural response and support rotation on a reduced diameter column in SPAn32. They concluded that the strength of concrete does not significantly affect the maximum rotation of the support, but it does affect breaching, i.e., the lower strength causes greater breaching.
Winget et al. [2] use SPAn32 to calculate the flexural response of the columns and to define the equivalent SDOF stiffness and mass parameters based on the column properties. For the calculation of the blast load pressure history, the BlastX program was used. Other useful programs are AT Blast [89] and Nonlin. For the calculation of the pressure-impulse history using AT Blast, it is necessary to know the charge weight, angle of incidence, and standoff distance. However, AT Blast does not consider the effects of multiple reflections under the bridge explosions. They list four categories of bridge design concerning their importance, where category 1 represents very important bridges, and category 4, unimportant. Winget et al. did not take the real conditions of the ground and energy absorption by creating craters but the ideal reflecting surface. Footing instability, however, could also result from large ground deformations, and this aspect of behavior must also be addressed.
Wu et al. [48] numerically simulated RC and composite columns in LS-Dyna for the contact-placed TNT charges from 2.5 to 25 kg. In the simulations, they obtained a higher residual bearing capacity of the column when the explosive was placed at the height of 1.5 m from the bottom than when it was placed at the bottom.
Fujikura et al. [7], for calculation of impulse variations per unit length along the height of the column, were using the Bridge Explosive Loading (BEL) [90] program. BEL also considers the reflected pressure of the blast wave on the surface of the superstructure and on the ground.
Rutner et al. [47] studied the behavior of four types of column cross-sections in MSC.Dytran software on blast load: single-cell hollow steel section, multi-cell hollow steel section, single-cell hollow composite column, and multi-cell composite column. Compared to steel columns, the composites showed negligible deflection on blast load. The best stress distribution in the element was achieved with multi-cell composite columns, which is also visible by the displacements shown in Figure 17. Hao et al. [49] analyzed three reinforced concrete columns with the same dimensions, material strengths, and reinforcement ratios but subjected to blast loads of different scaled distances. They wanted to find the failure probability using the computer code CARLER, which is verified with Monte Carlo simulations. They defined four damage levels (D) that depend on the ratio of residual axial load carrying capacity (N residual ) when the column is damaged and the axial load of the undamaged column (N 0 ), shown in Table 5. Through numerical simulations, they concluded that neglecting some of the material properties of the column has minimal impact on the probability of failure, while the random changes in the blast loading have a much greater role.  [51] emphasized the spalling of side-cover concrete because, in previous works, only the spalling of concrete off the back of reinforced concrete columns was mentioned. The aim of their research was to make and validate a numerical model with respect to the experiment explained in [38] and justified the cross-sectional response mechanisms that cause loss of side-cover concrete. For numerical simulations, they used the LS-Dyna program and the Karagozian and Case concrete (KCC) material model.
Numerical simulations in LS-Dyna showed that the shape of the column cross-section has a large influence on the resulting impulse. The authors developed expressions for calculating column shape factors for circular and square columns. The expressions are used when the R / D ratio is less than 4.5 because they provide sufficiently conservative loads but less than those experimentally determined on the walls [52].
Crawford [10] performed numerical simulations in LS-Dyna of RC columns and columns retrofitted with fiber-reinforced plastic (FRP). FRP increases the resistance of the RC columns to the blast load. For numerical modeling, the choice of concrete material model is very important, and the analysis was performed with four different concrete models (KC, Winfrith, Continuous Smooth Cap, and RHT model). The best results were obtained with the KC model.
Magali et al. [53] performed a parametric numerical analysis in Abaqus to see which of the six varied parameters (section ratio, compressive strength of concrete, column height and thickness, charge radius, and ratio between standoff distance to the charge radius) had the greatest impact on column damage. It was shown that the column thickness, charge radius, and the ratio of standoff distance to charge radius play a significant role in the column response. They give an empirical formula based on the conducted simulations to predict the damage index of the column. Comparing the results obtained by formula and numerical simulations, the deviations are up to a maximum of 15% what is in an acceptable range.
Eisa [54] modeled RC columns in Abaqus. The position of the spherical charge remained unchanged, but the charge weight, column height, longitudinal steel reinforcement, columns aspect ratio, and transverse steel ratio were varied. Fixed boundary conditions are provided at the top and bottom of the column. In addition to column damage, they measured acceleration and displacement concerning varied parameters. Four quantities of explosives were used (45.36, 226.80, 453.59, and 1016.05 kg) and placed at a distance of 4.87 m. The increase in the lateral reinforcement in the column had the effect of reducing the displacement in the middle of the column. It is recommended to examine the influence of the axial force during the blast load and to include additional parameters such as the variation in the standoff distance of the charge.
In numerical simulations using Autodyn, Abladey and Braimah [55] tested three columns designed for different loads in accordance with the Canadian reinforced concrete design code [91]. The first type of column is designed only for gravity loads, and the distance between the transverse reinforcement is 300 mm, and the other two types of columns are designed for seismic loads, but in the second type, the distance between transverse reinforcement is 150 mm, and in the third 75 mm. Column damage is significantly less with denser reinforcement, especially at small-scaled distances. They have proven in research that regardless of the same scaled distance, in a situation where a larger amount of explosive is detonated, the column has a higher deflection.
Li and Hao [56] calibrated the numerical model for RC slab in LS-Dyna according to a previously performed experiment and then used that numerical model for RC column simulation. For good simulation of concrete spallation, the erosion criterion using principle tensile strain of 0.01 was defined. Through simulations, they concluded that denser reinforcement and greater column depth reduced spall damage, i.e., increase the confinement of concrete. The boundary conditions and column height do not play a significant role in the level of spall damage in close-in cases.
Shi and Stewart [57] analyzed a spatial and non-spatial simulation of the blast load on RC columns in LS-Dyna. They used three quantities of ANFO explosives, 50, 100, and 1000 kg, at distances from 0 to 30 m. The analysis is based on axial load-carrying capacity and concluded that the variability of the results in the spatial model is lower, and the probability of damage is significantly higher. They consider the spatial model more reliable and recommended it for future research with an additional assessment of the scale of fluctuation.
Elsanadedy et al. [50] used LS-Dyna for analyzing the behavior of the exterior building RC circular column and strengthened column with Carbon Fiber Reinforced Polymer (CFRP) sheets under blast load. Four different charge weights (100, 200, 500, and 1000 kg) of TNT at three different standoff distances (1, 4, and 15 m) and at 1 m height from the ground were analyzed. They modeled columns with different boundary conditions, first with both fixed ends and a second type with both hinged ends. For calculation of blast load parameters in all assumed scenarios, they used the software ConWep. The use of CFRP is increasing the shear capacity of the column and the strength of the column, which results in less lateral displacement, and more layers of CFRP can undergo more intense blast loads [50]. Moreover, at the scaled distance from 0.50 to 0.68 m/kg 1/3 , columns with CFRP showed better behavior than RC columns.
Ashalekshmi and Subha [92] modeled a bridge column in the Ansys Autodyn software to analyze the impact of concrete grade and spacing of ties under the blast load. They observed total deformation and principal stress for concrete grades M40 and M50 and ties spacings of 10 and 20 cm in the near and far-field. In the near field, the explosive is placed at the same standoff distance of 2 m, and the charge weight varies from 250 to 1500 kg. The influence of the concrete grade on the maximum deformation is visible only at the weight of explosives greater than 750 kg. With concrete M50, the principal stress is higher, but there is no big difference in it when increasing the weight of explosives. In the far-field, the explosive was placed at a standoff distance of 10 m, and the weight also varied from 250 to 1500 kg. There is a slight difference in total deformation concerning the grade of concrete. The principal stress exceeds the strength of concrete, and the difference in stress is visible for the grade of concrete. The effect of tie spacing is visible in the near field only in the increase in the maximum principal stress, while in the far-field, there is no effect on either deformation or principal stress.
Liu et al. [58] modeled the bridge column and bent it in the LS-Dyna for three design categories provided in [35]. They determined six damage mechanisms in the models, four of the column and two of bent. In all three design categories (A, B, C), spalling of concrete and crushing of the bent concrete were observed. Plastic joints in the column and shear of bent occur only in B and C categories. The shear or flexure failure of the column is most probable in category C, i.e., at the highest blast load. They found that increasing the transverse reinforcement reduces damage. In LS-Dyna, they received underestimated blast load, and therefore they used the ConWep program to calculate pressure-time diagrams.
Cui et al. [59] concluded through numerical simulations of columns in LS-Dyna that a larger cross-section and reinforcement ratio, smaller spacing between the stirrups, and a thinner concrete cover for columns exposed to close-in explosions give less damage.
Zhang et al. [12] analyzed circular and square CFST (Concrete Filled-Steel Tube) columns with tube thicknesses of 2.8 and 3.8 mm. The Emulsion explosive was used, with TNT equivalences of (0.7) 17.5, 25, and 35 kg. A numerical simulation of the columns was performed in LS-Dyna, see Figure 18, but to reduce the computation time and model congestion, an air blast model was made in the ConWep program and then imported into LS-Dyna. By comparing the obtained periods of oscillation and maximum displacements, a good match between the numerical model and the experiment was obtained. Differences are found only in residual deflections but are not considered crucial to the accuracy of the model. Zhang et al. [21] numerically modeled CFDST columns with steel-fiber reinforced concrete using LS-Dyna and compared them with experimentally obtained results. For concrete, the KC model was used, and for the steel tube, the Material model 24 was used. The parameters for the concrete model were modified (f t -tensile strength, B 1 -parameter for residual strength, w lz -strain softening, ω-confinement effect, λ, and η-damage parameters) because UHPC was used. Column erosion occurs when the maximum shear strain value reaches 0.045. Numerical research concluded that the axial load up to a certain limit has a favorable effect on the deflection in the middle of the column. It is not recommended to use columns with a hollow section ratio greater than 0.5. Increasing the thickness of the outer steel tube affects the decrease in deflection in the middle of the column. Columns filled with UHPC have less plastic deformation than columns filled with normal strength concrete, shown in Figure 19. Codina et al. [14] compared the numerical simulation with an already conducted experiment. The used explosive was Gelamon VF65 (65% TNT equivalence), 8 kg TNT equivalence at distances of 1 m and 0.60 m from the column, which is classified as a near-field range. Overpressure, impulse, deflection, acceleration, and visual damage to the column were measured and compared. For numerical modeling, the Ansys Autodyn program was used in which air and explosives were modeled through an Euler processor and the column through Lagrange. The optimal mesh size of concrete, steel, and the air was 10 mm. The authors made models with default RHT values and with the values proposed in [93] but concluded that the parameters that affect the strength degradation (damage factors D 1 and D 2 , and e f ail min ) and the residual strength, Y* fric (parameters B and M) should be changed in the model. For good prediction of spallation, the instantaneous geometric strain was used for erosion type with a value of 0.5. The obtained column damage with different RHT parameters is shown in Table 6. The parameters of the RHT material model (shown in Table 3) are validated for scaled distances (Z) from 0.5 to 0.3 m/kg 1/3 . Table 6. Column damage obtained: (a) experimentally and numerically using (b) default RHT model parameters; (c) parameters provided in [93]; (d) modified parameters by [14]. Arowojolu et al. [60] studied, using LS-Dyna numerical models, the influence of axial and blast load on the RC column of the bridge. For the concrete model, they used CSCM (Continuous Surface Cap Model) and for reinforcement MAT 24. Exact quantities of explosives and distances are not given, but scaled distances from 1.77 to 0.45 m/kg 1/3 . They concluded that when an axial load ratio of 0.25 is applied, the displacement in the middle of the column decreases but the damage of the column increases.

Parameters
Eamon and Aslendi [61] made a numerical model of the column, experimentally tested in [4] using LS-Dyna software. The influence of concrete strength, reinforcement ratio, axial load, and the column wrapping with SFRP (Steel Fiber Reinforced Polymer) was observed. The Johnson-Holmquist-Cook (JHC) model was used for concrete modeling and the elastic-plastic kinematic model for steel reinforcement. The impact of the blast load was determined in the ConWep software, and the detonation point was placed 5 cm above the ground, 40 cm from the column. SFRP proved to be an inexpensive and ductile retrofit. One layer has the largest contribution in blast capacity, while all additional layers have a small effect on increasing the capacity. They obtained a linear relationship between the concrete strength and the increase in the blast load capacity.
Kravchenko et al. [62] performed numerical simulations of the RC column in LS-Dyna. The concrete was modeled using the CSCM concrete model (type 159 material) and reinforced using the plastic-kinematic model (type 3 material). They observed the influence of detonation of 10 kg of TNT at a distance of 1.2 m from the ground and 1 m from the column. They also concluded that the reinforcement ratio has a significant impact on the behavior of columns under the blast load and that the columns in the ground floor and bases need to be better reinforced due to their easy accessibility.
Kyei and Braimah [63] modeled three RC columns in LS-Dyna, which differ in the distance between the transverse reinforcement. They designed the columns according to the instructions for the level of seismicity in the Canadian concrete design code [91]. Concrete was modeled using the Continuous Surface Cap Model (MAT_CSCM_159), and for reinforcement, they used Material Piecewise Linear Plasticity (MAT_024) model, while blast load was calculated in ConWep and then imported with Load Blast Enhanced (LBE) in LS-Dyna. They performed simulations with mesh sizes from 5 to 100 mm, and with 15 mm, they obtained a good ratio of the time spent for the calculation and the accuracy of the results compared with the experiment in [94]. The used explosive was ANFO (100, 250, 500, and 1000 kg) at scaled distances of 0.8, 1.0, and 1.5 m/kg 1/3 . In near-field explosions, the distance between the transverse reinforcement has a significant effect on the reduction in displacement, while in far-field explosions, this effect is negligible. At a high axial load ratio (0.35), the seismically designed columns showed better behavior at scaled distances than the standard ones.
Yuan et al. [23] performed numerical simulations in LS-Dyna of the circular and square columns of the bridge, exposed to the contact explosion of 1 kg of TNT. To reduce the computation time, at the height of 1 m, in the area of the contact explosion, they placed a denser mesh (8 mm), while on the rest of the column, the mesh size was courser (20 mm), shown in Figure 20. The principal strain of 0.5 was used as the erosion criterion. Numerical simulations well described the damage on the front sides of the column, while on the back, there are differences. They concluded that the damage was greater on the square column due to the flat surface and the higher stress concentration than on the circular column. The damage of the columns is shown in Figure 21.  Yuan et al. [77] investigated the effect of axial load on RC bridge columns subjected to far-field, close-in, and contact explosion. Columns in the far-field have mainly a flexure response, and the axial load affects the reduction in the maximum displacement in the middle of the column. In the case of a close-in detonation, a shear failure is expected, and the axial load affects the increase in the damage of the column and should not be neglected. In contact detonation, the concrete covers at the front and backside of the column spall off. The axial load reduces the damage of the concrete but increases the stress in the reinforcement and must be considered.
Li et al. [65] conducted numerical investigations on CFDST columns under contact explosion in LS-Dyna. CFDST columns have proven to be good for two reasons: the first is that the confinement of concrete by steel tube allows better energy absorption, and the second is the prevention of the spallation of the concrete cover. In work [73], the researchers concluded that increasing the cross-sectional area and the ratio of reinforcement plays a significant role in the post-blast residual capacity of CFDST columns under contact explosion. In [67], the behavior of CFDST columns subjected to close-in blast loading was studied. They concluded that the influence of the charge shape significantly affects the response and behavior of the column at scaled distances from 0.079 m/kg 1/3 to 0.175 m/kg 1/3 .
Liu et al. [66] performed a dynamic and static analysis of the bridge columns in LS-Dyna. Dynamic analysis is based on the comparison of accelerations and static analysis on the determination of damage through the ratio of the residual to the ultimate axial bearing capacity. They examined the damage to the column concerning the position of the explosives, at the bottom, in the middle, and at the top, shown in Figure 22. In all three cases, bending deformation occurs at the column, but the position of the crack formation differs. Liu et al. [68] performed a parametric analysis in LS-Dyna to determine the influence of longitudinal and transverse reinforcement ratios, longitudinal force, and boundary conditions. They concluded that at smaller-scaled distances, the increase in transverse and longitudinal reinforcement reduces the displacement in the middle of the column. It is recommended that the percentage of longitudinal reinforcement does not exceed 6% of the cross-sectional area of the column because too much reinforcement can lead to brittle failure. Analyzing the influence of the axial compressive load, they found that in an amount of up to 40%, it reduces the maximum displacement in the mid-span of the column due to the increase in moment capacity. The conclusions are based on close-in blast loading; a recommendation for future research is to conduct a parametric analysis for near and far filed scenarios.
Liu et al. [69] made numerical models of RC columns strengthened with CFRP in LS-Dyna. The size of the mesh elements for all materials was 10 mm, and for air, 20 mm. The 1 kg and 2 kg TNT charges were placed in contact with the column at the height of 30 cm. The numerical results show that the dragging force of the blast load separates the CFRP from the concrete, but despite this, CFRP protects the column from contact explosion. Column damage is shown in Figure 23. Thai et al. [70] modeled in LS-Dyna RC steel jacked columns (25 cm × 25 cm × 360 cm) and observed the influence of steel thickness, scaled distance, and longitudinal compressive force on the behavior of the columns under the blast load. The columns are designed according to Eurocode 2. A 10 mm mesh was used for the column, while 5 mm elements were used for the explosive. Placing 8 kg of TNT in the middle of the column causes less global damage to the column, while placing the same amount at a distance of 32 cm from the ground causes significant local damage. The scaled distance was varied from 0.10 to 0.40 m/kg 1/3 . Increasing the steel thickness from 3 mm to 6 mm did not result in less damage.
Abedini et al. [33] investigated numerically in LS-Dyna the influence of charge and scaled distance on the level of column damage and influence of column width, concrete strength, and reinforcement ratio on the residual axial load capacity.
Dua et al. [72] performed a parametric analysis in LS-Dyna on the RC columns exposed to the contact explosion. They used from 115 g to 1000 g TNT, varied the compressive strength of concrete from 20 to 55 MPa, and reduced the distance of the transverse reinforcement from 200 to 50 mm. Increasing the transverse reinforcement reduces the damage to the concrete core, and higher compressive strength of concrete contributes to the reduction in cracks and peeling. The columns under contact explosion have local damage, while the global damage is negligible. Dua et al. [28] investigated in LS-Dyna the influence of the column cross-sectional width under the contact explosion. They concluded that a larger cross-sectional width (greater by two or more times from the depth) has a favorable effect on the behavior and damage of the column.
Rajkumar et al. [74] examined 45 numerical models of reinforced concrete columns in scale 1:4 (85 mm × 85 mm × 900 mm) in LS-Dyna. In the models, the behavior of four different cross-sections (circular, octagonal, hexagonal and square) on the blast load was examined. The circular column retains the lowest peak pressures and has the smallest deflection in the middle, while the square has the highest pressures and the largest deflection. The edges in cross-sections at small-scaled distances play a significant role in pressure retention, while with increasing scaled distance, the shape influence decreases. An increase in reinforcement in cross-section affects the improvement of the behavior of the columns during the blast load, especially in shapes that retain higher pressures.
Vavilala et al. [75] numerically simulated polymeric foam retrofitted RC columns (35.5 cm × 35.5 cm × 348 cm) in Abaqus. They used a 10 mm mesh for reinforcement and 20 mm for concrete. Columns coated with 5 mm, 8 mm, and 10 mm thick foam were exposed to 10, 25, and 50 kg of TNT. The greatest reduction in displacement in the middle of the column was obtained under 10 kg of TNT when the column was coated with 10 mm thick foam.
Zhang et al. [76] compare in LS-Dyna the behavior of a segmental CFST column with monolithic and prestressed monolithic columns. The columns were exposed to 20 and 50 kg of TNT at a standoff distance of 1.5 m. The segmental column showed a smaller residual displacement, and numerical analysis proved that a larger number of segments had a more favorable effect on the behavior of the column. Moreover, increasing the steel thickness had a beneficial effect on reducing concrete damage.
Yan et al. [78] and Hu et al. [79] used LS-Dyna for numerical simulations of RC columns retrofitted with CFRP subjected to the close-in explosion. CFRP sheets reduce the deformation and spalling of concrete. The CFRP thickness, wrapping, and dimensions ratio of the charge also had a large impact on the damage of the column and peak pressure. Debonding of CFRP is the most common form of failure, but despite this, CFRP has a role in reducing column damage during the blast load, and it may even have a role in changing the failure mode from shear to flexural deformation. The setting of the CFRP on both sides of the column needs to be further investigated because direct shear is possible due to over-reinforcement.

Experimental Testing
Most of the experiments were conducted for the scenario of an attack by an auto-bomb located near a bridge or building column. Experimental tests on building columns are mostly full-scale, while tests on bridge columns predominate on scales of 1:3 and 1:4, due to the high cost of performing such experiments, the need for trained personnel to handle explosives, large quantities of explosives, and the field where such tests can be carried out.
Experiments showed that even a minimal increase in the cross-sectional dimensions of a column could favorably affect the behavior of the column under the blast load. Moreover, a minimal increase in the standoff distance reduces the impact and intensity of the blast load, and therefore it is necessary whenever possible to fence the column, increase visibility around the column, and reduce its accessibility. In addition to the dimensions of the column, the shape of the column plays a significant role. Circular columns retain less impulse from the blast load than square ones of the same dimensions. Squares have a larger cross-sectional area and, therefore, can better withstand shear. On the action of the contact explosion, the circular column suffered less damage than the square. However, comparing a square and a rectangular column, when the width dimension of the rectangular column is greater two and more times from depth, damage of the column is lower. The shape of the column on the blast load needs to be further investigated. In RC columns, position, quantity, anchoring, and reinforcement overlap have a great impact. Therefore, seismically designed columns have better blast load behavior than standard designed columns.
Comparing steel jacked columns, CFST and CFDST with RC columns, all showed better ductility, less cracking of the concrete, and the absence of flying debris. However, CFDST showed the best behavior, as the inner steel tube contributes to the prevention of direct shear. With this type of column, it is important that the ratio of the cavity to the column cross-section is not greater than 0.5 and that the thickness of the steel tube is well determined (a thickness of 3.8 mm gave satisfactory results).
Columns with Ultra-High-Performance Concrete and with various Fiber-Reinforced Polymers showed better load-bearing capacity than High-Strength Reinforced Concrete.
The influence of axial load up to 30% of the total load capacity has shown a favorable effect on the reduction in the maximum displacement, but this percentage needs to be further investigated.
Recent research applies smart aggregates to measure internal cracks and internal damage because the column can be damaged and reduced load-bearing capacity without being visible on the outside. Therefore, post-blast tests are performed to determine the residual strength and ductility of the column.

Numerical Modeling
Numerical simulations make it possible to study the effect of large amounts of explosions (more than 1000 kg) on columns in full-scale. The most widely used software for analyzing the nonlinear behavior of elements on blast load is LS-Dyna. The most accurate simulation of the blast pressure requires as finer mesh as possible, which leads to a long duration of simulations and congestion of the computer processor. Due to the use of a larger mesh size than recommended, the LS-Dyna underestimates the pressures. Most researchers use the ConWep program to calculate pressures and import the resulting pressures into LS-Dyna.
Defining a model of concrete material is the most demanding because many parameters affect its behavior. Karagozian and Case (KC) concrete is mostly used in LS-Dyna. The proper definition of erosion criteria has proven to be very important in modeling column damage.
The position of the explosive plays a significant role in the behavior of the column. Placing the charge in the far-field causes a uniformly distributed load per column and global response, while a charge placed in the near-field, close-in, and contact creates local damage.
Moreover, the columns exposed to the charge placed in the lower half showed greater damage and lowered residual capacity than the columns where the charge was placed in the middle.
In all fields (far, near, close-in), changes in the quantity, shape, and position of the explosives showed a great impact on column behavior. The shape of the column, the ratio of reinforcement, and the concrete grade showed an influence only at small-scaled distances. The concrete grade does not affect the rotation of the column but does affect the reduction in concrete breach and spallation.
Until recently, the axial load on bridge columns was neglected in the calculation because it was considered to be on the safe side. However, numerical simulations showed that axial load has a large impact on increasing damage when the charge is placed near the column or in contact with the column and should not be neglected. Data Availability Statement: Data available on request due to restrictions, e.g., privacy or ethics. The data presented in this study are available on request from the corresponding author. The data are not publicly available due to limited access to certain literature used in the article.