Improved SDOF Approach to Incorporate the Effects of Axial Loads on the Dynamic Responses of Steel Columns Subjected to Blast Loads

In this paper, a complicated single-degree-of-freedom (SDOF) approach was developed to determine the global response of steel columns under combined axial and blast-induced transverse loads. Nonlinear section and member analyses were incorporated into the suggested SDOF method to account for the complex features of the material behavior, the high strain rate effect, and the column geometry. )e SDOF technique was validated through comparisons with available finite element and experimental data, and a good consistency was obtained.)en, the validated SDOF approach was utilized to derive the pressure-impulse curves under various levels of axial loading. )e level of the axial load was shown to have a significant influence on the dynamic behavior of a steel column subjected to a blast load.


Introduction
Recently, increasingly frequent terrorist activities have made blast protection essential for engineering communities [1].
e failure of load-bearing columns can result in the progressive collapse of structures.erefore, it is important to improve the reliability of vulnerable columns under blasts to ensure the safety of structures [2,3].
One of the most widely used approaches to assess the effects of blasts on structures is the equivalent single-degreeof-freedom (SDOF) approach.In addition, the efficient SDOF method has proven to be a powerful tool for deriving isodamage pressure-impulse (P-I) diagrams, which can predict structural damage.For example, several technical manuals such as the unified facilities criteria (UFC) [4] based on linear resistance curves are used to design protective structures subjected to blast loads, and idealized load-deformation characteristics are given in the form of design charts.Dragos and Wu [5] proposed a simplified SDOF approach based on the new concept of a reduced resistance-deflection function to derive the P-I diagram for steel columns.Fallah and Louca [6] proposed analytical formulas for obtaining a normalized P-I diagram based on a simple SDOF method incorporating the linear resistance.Further, Nassr et al. [7] demonstrated the SDOF method's ability to study the response of steel beams subjected to blast loads, and Nassr et al. [8] later employed the SDOF approach to study the axial load effects on steel columns under blast loads.In addition, Wang et al. [9] studied the SDOF approach to evaluate the concrete structural members under close-in nonuniform blast loads, and a comparison with test results indicated that the SDOF method can give reasonable predictions.Bedon et al. [10] carried out the analytical calculations by means of SDOF formulations derived from structural dynamics theories and proposed a design approach for practical estimation of buckling strength of glass elements in the analyzed loading and boundary conditions.Liu et al. [11] performed the field tests to investigate the blast response of steel-reinforced concrete (RC) beams and validate the nonlinear SDOF approach to predict the maximum displacements of RC beams subjected to blast loading.However, the axial load effect was neglected.
Currently, structural elements are typically simplified as having perfectly elastic-plastic force-deformation relationships [7,8].However, this idealized member resistance has inherent limitations for capturing the gradual effect of plastic behavior on the dynamic properties of the system [9,12].erefore, employing this force-deformation characteristic causes some uncertainty in the calculations.
Further, in reality, when subjected to a lateral blast pressure, an axial load tends to induce additional lateral deformation and magnify bending within a column.us, a structural response to blast loads that ignores axial load effects may not very accurately reflect the true behavior of a structure [13].
In this paper, a more accurate and sophisticated flexural behavior (featuring both the moment-curvature relationship and the force-deformation relationship) that considers the complex features of the slab geometry and the material's behavior under high strain rate loading is incorporated into a nonlinear SDOF approach for the dynamic analysis of columns subjected to blast-induced shock waves under various levels of axial compressive loads.
is approach also has the capability to generate isocontour P-I diagrams of blast-loaded columns.To validate the suggested SDOF approach, the structural loaddeformation characteristics, predicted maximum displacements, and iso-displacement P-I diagrams in this paper were compared against experimental and finite element model (FEM) results from [5,14], respectively.

Theoretical Analyses of the Response of a Steel Column
2.1.Sectional Analysis.Assuming that plane sections remain planar under purely flexural behavior, a sectional analysis is performed to compute the moment-curvature characteristics for use in member analysis [15] as shown in Figure 1.
A solution procedure to generate the moment-curvature is outlined below, as shown in Figure 2.
(1) e cross section of the column is divided into a desirable number of slices n parallel to the column axis.(2) Strain values ε top are specified for the extreme steel compression fiber, and a corresponding depth x n for the neutral axis is selected.e strain profile within the section is then generated: where x i is the distance between the corresponding layer and the neutral axis.(3) After the strain profile within the section is generated, the internal forces on the cross section are generated by employing the material's stress-strain relationships: where b is the width of the section, σ i is the stress on the i-th layer of the steel, and h i is height of the i-th layer of the section.(4) e depth of the neutral axis is revised until the difference between the applied axial load N and the resultant internal force is within a convergence criterion: where ΔA i is the area of the i-th layer of the section.
(5) e resulting moments and curvatures are then calculated using these internal forces and strain profiles: (6) Another incremental extreme compression fiber strain is selected, and the procedure is repeated until the extreme compression of concrete up to its ultimate compressive strain is reached for a given axial load.

Material Modeling.
e stress-strain relationships of the member's material are incorporated into a member analysis to relate the member's deformation to the developed internal forces.
e stress-strain curve is identical to the material properties described in [15], which have been thoroughly validated.

Steel.
Compressive stress-strain relationships are assumed to be the same as tensile stress-strain relationships.
e stress-strain relationship for reinforcing steel in tension consists of two linear segments, as shown in Figure 3.

High Strain Rate Effect on Steel.
Accounting for the effects of high strain rates is critical in the analysis and design of members subjected to blast loads, as materials experience an increase in strength under high strain rates.
e dynamic increase factor, DIF, is defined as the ratio of the dynamic material strength to the static material strength.Based on the relationship proposed by Malvar and Crawford [16], the dynamic material increase factors for both the yield strength and the ultimate strength are given by the following: where the DIF for the yield strength of steel may be determined by substituting α � α f y into equation ( 5) as follows: where f y is the static yield strength of steel (in MPa).Similarly, the DIF for the ultimate strength of steel may be determined by substituting α � α f u in equation ( 5) as follows: 2 Advances in Civil Engineering

Member Analysis.
To capture the full time history response of a steel-reinforced concrete column subjected to dynamic loading, the complete force-displacement relationship rather than the ultimate resistance is needed.e moment-curvature curve generated by sectional analysis is used to perform the member analysis of the steel-reinforced columns.e following discussion outlines a procedure for obtaining the complete nonlinear resistance functions for steel concrete columns, which is also shown in Figure 4: (1) Establish the moment-curvature relationship for the section.A column with a length l is divided equally into a desirable number of segments (equal to m), each of which has a length Δl, as shown in Figure 5. Assume a curvature on the midspan of the column φ m ; the initial φ m is zero, and at every step, φ m � φ m + Δφ. e moment M m on the midspan is then computed based on the moment-curvature relationships.m � 500 was used for the analysis.Convergence study shows that further decrease of mesh size only has little effect on the numerical results but leads to a much longer calculation time.(2) Assume an initial incremental distributed deflection δ i,k (i � 0, 1, . .., m).(3) e assumed horizontal load P is calculated by the axial load P a , M m , and δ i,0 : (4) e incremental distributed moment M i is calculated by P, P a , and δ i,k .en, the distributed curvature φ i is determined based on the moment-curvature relationship: (5) In the conjugate column method, each segment divided above is considered a small imaginary column having the same length as the original segment.e midspan displacement of the segment can be calculated by the following expression: Advances in Civil Engineering 3 (6) Assume the new incremental distributed deflection is δ i,k+1 , repeat steps ( 4)-( 6), and then determine δ i,2 , δ i,3 , . ..δ i,k , δ i,k+1 .If δ i,(k+1) and δ i,k are close enough (where the difference is less than an allowable value), then enter the next cycle.Finally, the whole loaddeflection curve is determined [1].

Dynamic SDOF Analysis.
e dynamic response of the columns was predicted using dynamic inelastic singledegree-of-freedom (SDOF) analysis after defining the nonlinear resistance curve for the structure [11] e SDOF system is shown in Figure 6. e idealized SDOF system is described by the following equation of motion: where x(t) and € x(t) represent the displacement and acceleration, respectively, of the column at the midspan, M is the mass of the column, R(x) is the resistance as a function of the displacement, and F c (t) is the loading as a function of time.
K M is the mass factor required to transfer a system with a distributed mass into an equivalent SDOF system and is computed by the following equation: where m is the mass per unit length of the member and Φ(x) is the shape function of the member.K L is the load factor used to transform the distributed force and stiffness into a single point load and is defined as follows: It is convenient to define a load mass factor [2], K LM , as a ratio of the mass factor to the load factor: Moreover, the response of the system may be written as follows by incorporating the axial load effect: where P a is the axial load.e load mass transformation factor, K LM , is defined based on the Biggs analysis [17].Newmark proposed the following numerical solution technique for this equation [18]: (1) Assume that y , _ y, and € y are known at the time t � t i and assign them the subscript i.
(2) t i+1 � t i + Δt and assume a value € y i+1 for the acceleration at time t i+1 , _ Check equilibrium: Calculate and store displacement and load  Advances in Civil Engineering solution is achieved.If they are not close, repeat steps 2 through 6 until convergence is attained.Since one needs to assume a new value for € y i+1 in step 2, it is quite efficient to select the value obtained in step 5 for this purpose.e β factor determines how the acceleration changes within a time step.In step 2, it is assumed that € y varies linearly over Δt. us, the default value for β is 1/6 [19].

P-I Diagram.
A P-I diagram as an iso-response curve is typically employed to assess the dynamic response of a structure subjected to a blast load.P-I curves consist of combinations of the peak reflected overpressure and the impulse that cause the same damage level based on a predefined failure mode.
For this study, the midspan displacement of the column is considered the damage criterion.One point on the P-I diagram is generated by specifying an initial positive load duration and iterating a reflected overpressure until the maximum midspan displacement predicted by the SDOF analysis reaches the predetermined damage value.

Validation of the SDOF Method Using FEM Results
Available in [1].Dragos and Wu [5] performed an FE analysis of two steel columns under blast loads.One of the two steel columns presented in [5] was also utilized in this study to validate the present SDOF approach.e span of the abovementioned steel column was 4 m, and the width and thickness of the flanges were 102 mm and 10.3 mm, respectively.e columns were simply supported.e total depth of the section was 160 mm, and the thickness of the web was 6.6 mm. e yield strength and Young modulus of the steel specimens were 470 MPa and 180 GPa, respectively.A 359 kPa reflected pressure and 23.88 ms positive blast duration for the blast conditions were chosen for the calculation.e maximum displacements of the steel column specimens under an axial load of 350 kN were calculated by employing the proposed SDOF analysis.
Figure 7 shows a comparison of the deflection time history between the calculated and FEM results [5].e calculated and FEM midspan maximum displacements of the steel column are 208.3mm and 207.7 mm, respectively.However, a significant discrepancy occurs in the initial portions of the midspan deflection histories by both analyses.
e negligence of membrane action and damping effects is a possible reason for the significant discrepancies in the initial portions of the midspan deflection histories by both analyses.
P-I diagrams for the maximum midspan displacement of 208.8 mm calculated above using the proposed SDOF approach were generated under an axial load of 350 kN. Figure 8 compares the P-I diagram of a column obtained by the proposed method with the FEM results available in [5].As the figure illustrates, the P-I combination in the present SDOF analysis is very close to the FEM results from [5]. is result indicates that the present SDOF model gives a reliable prediction of the P-I diagram of a steel column subjected to blast loads.

Validation of the SDOF Method Using the Experimental
Results in [14].To further validate the accuracy of the abovementioned SDOF approach, maximum midspan displacements generated from the proposed SDOF analysis were compared with the results of the test reported by Nassr et al. [14].
irteen typical wide-flange steel columns with two different section sizes, W150 × 24 and W200 × 71, were chosen by Nassr et al. [14] to investigate the dynamic responses of the steel columns to blast loads.Only the steel columns with a section size of W150 × 24 and a height of 2413 mm under the blast conditions in shot 1 and shot 3 are analyzed herein.All of the steel columns were pinned at the top and supported by rollers at the bottom.erefore, the load mass factor, K LM , was taken as being equal to 0.78 and 0.66 before and after yielding, respectively.e axial load in the experiment was 270 kN for the two specimens.
Figure 9 shows the layouts of the steel column's sections [14].
e static nominal yield strength and ultimate strength of the test specimens were 393 and 537 MPa, respectively.In shot 1, a charge size of 50 kg at a stand-off distance of 10.3 m was used to generate the blast environment.
e average of the measured reflected pressures was 307 kPa, and the positive phase duration had an average value of 7.3 ms.In shot 3, a pressure time history was obtained by a charge size of 150 kg at a standoff distance of 9 m. e peak pressure and positive phase Advances in Civil Engineering duration were 1560 kPa and 6.2 ms, respectively.A typical reflected pressure time history obtained in shot 1 is shown in Figure 10.
Comparisons of the calculated and test deflection time histories at the middle of the column are shown in Figure 11.
e differences in the maximum deflections between the present analysis and experiments are 2.5% and 2.6% for shot 1 and shot 3, respectively.is result indicates that the SDOF model also gives reliable predictions of the responses of steel columns to blast loads.However, significant discrepancies occur in the initial portions of the midspan deflection histories by both analyses.e difference in the load-time histories under the real blast load and the SDOF analysis represents one possible reason for the reduced accuracy.e reflected pressure as a function of time was modeled using an equivalent triangular blast load having the same peak reflected pressure and impulse found in the experiments.us, more accurate load-time histories should be employed.In addition, the shear damage mode is not considered in the present study; accordingly, a more accurate SDOF model considering other damage modes should be analyzed in further research.

Effects of Axial Loads.
When a column is subjected to blast loads, the static compressive axial load induced by the gravity load affects the dynamic response of the column.e specimen analyzed above (W150 × 24) and the blast conditions (shot 1 and shot 3) available in [14] were employed here to investigate the effects of axial loads on the dynamic responses of steel columns.e effects of axial loads were investigated by applying three different axial loads, namely, 0 kN, 270 kN, and 400 kN.

Effects of Axial Loads on the Maximum Displacements.
Figure 12 illustrates the effects of different axial loads on the displacement responses of columns under various load conditions.When the impulsive loading and corresponding deformation are small, the midspan displacement decreases with an increase in the axial load because the axial load amplifies the moment capacity of the column.However, when the column experiences large deformation, an increasing axial load results in a greater midspan displacement.
is is caused by the P-Δ effect, which is that when a column is subjected to both lateral and axial loading, the axial load creates an additional bending moment and thus enhances the lateral deflection of the column.e blast load induces lateral deformation in the column, causing the applied and resisting axial forces to form a couple, as they are no longer collinear.is couple causes additional lateral deformation and bending in the column, thereby inducing a greater midspan displacement.6 Advances in Civil Engineering

Effects of Axial Loads on P-I Diagrams.
To investigate the effects of axial loads on the P-I diagram of a steel column, P-I diagrams were generated under the abovementioned axial loads.e maximum midspan displacements of 10 mm and 90 mm were determined as damage criteria.Both P-I diagrams are shown in Figure 13.As seen from Figure 13(a), with a decrease in the axial load applied to the columns, the P-I diagram shifts to the left. is result indicates that an increase in the axial load can increase the ability of a column to resist blast loads.However, Figure 13(b) indicates that a larger axial load can decrease the ability of a steel column to resist blast loads when the maximum midspan displacement is large because of the P-Δ effect.

Conclusion
In this paper, an SDOF model was developed to determine the global response of blast-loaded steel columns.
e model also incorporates nonlinear momentcurvature and force-deformation relationships instead of a perfectly linear elastic-plastic force-deformation relationship.
e proposed member analysis accounted for the material constitutive model, the strain rate effects, the gradual effect of the formation and progression of plastic behavior, and the effects of various axial loads.To validate the SDOF method, the theoretical SDOF maximum displacements were compared with the FEM and experimental results in [5,14], respectively.e results Advances in Civil Engineering agreed well, indicating that the proposed SDOF model can accurately predict the deflection and P-I diagram of a blast-loaded column.e maximum midspan displacements and P-I diagrams for steel columns under various axial loads were generated based on the proposed SDOF approach.e level of the axial compressive load was found to have a significant effect on the response of a steel column.When the level of the blast load is small, an increase in the axial load can reduce the maximum displacement and enhance the ability of the column to resist the blast load, which is caused by the enhanced moment capacity induced by the axial load.However, when the level of the blast load is sufficiently large, an increase in the axial load can increase the maximum displacement and reduce the ability to resist the blast load, which is caused by the P-D effect.

Data Availability
e origin data used to support the findings of this study are included within the article.e Matlab code data used to support the findings of this study were supplied by Junbo Yan under license and thus cannot be made freely available.Requests for access to these data should be made to Junbo Yan (e-mail address: 3120140060@bit.edu.cn).

Figure 2 :Figure 1 :Figure 3 :
Figure 2: Flow chart of the numerical procedure for section analysis.

Figure 4 :
Figure 4: Flow chart of the numerical procedure for calculating the resistance function.

Figure 10 :Figure 11 :
Figure 10: A typical reflected pressure time history obtained in shot 1.

Figure 13 :Figure 12 :
Figure 13: P-I curves for varying axial loads.(a) P-I curves for a 10 mm midspan displacement.(b) P-I curves for a 90 mm midspan displacement.