Research Article  Open Access
A New Simple Method for the Strength of HighStrength Steel ThinWalled Box Columns Subjected to Axial Force and Biaxial End Moments
Abstract
The strength of Q460 steel welded thinwalled box columns under biaxial bending is investigated by an FE model. The numerical results, together with available experimental data, are compared with the American Specification ANSI/AISC 36010. It shows that ANSI/AISC 36010 provides a good estimate over a wide range of column slenderness ratios; but it severely underestimates the strengths of beamcolumns with large slenderness ratios and slightly underestimates those of beamcolumns made from HSS with f_{y} = 741 MPa, and the calculation is more complicated because of using an effective width concept. Therefore, a simple and effective calculation method is proposed. Also, a comparison is made between the proposed formulas and the available experimental results. It shows that the proposed formulas can precisely evaluate the localoverall interactive buckling strength of HSS beamcolumns and, in most cases, are also suitable for beamcolumns fabricated from mild steels. The proposed formulas are simple and, meanwhile, achieve the same level of accuracy as ANSI/AISC 36010.
1. Introduction
Structural steel with a nominal yield stress (denoted by f_{y}) greater than or equal to 420 MPa is often defined as highstrength steel (HSS). HSSs, such as Q460 and Q690, the nominal yield stress of which is 460 MPa and 690 MPa, respectively, have been successfully applied in largespan structures, highrise buildings, and transmission towers in China owing to having high plasticity, good toughness, and excellent weldability. Economical design of HSS structures leads to thinwalled members. These thinwalled members will, very likely, fail in a mode of interaction between local plate buckling and overall member buckling, also called localoverall interactive buckling, in which local buckling of plate elements occurs prior to overall buckling of members and the local buckling weakens the overall stiffness and reduces the overall buckling strength.
A large amount of research work has been done on the localoverall interaction buckling strength of welded I and boxsection columns [1–14], and beamcolumns [2, 6, 11, 12, 15–23]. Nonetheless, at present, there are a few research studies about HSS beamcolumns [11, 12, 20–23]. In the early 1980s, an experimental study was performed by Usami and Fukumoto [11, 12] on the localoverall interaction buckling strength of HSS welded box columns, and an empirical formula for the buckling strength was suggested. A total of 23 specimens with f_{y} = 460 MPa and 27 specimens with f_{y} = 690 MPa were tested to failure under concentric load and uniaxial eccentric load; material properties, residual stresses, and column initial deflections were also measured. The column slenderness ratios ranged from 30 to 65 and the plate widthtothickness ratios from 22 to 58. Until recent years, a series of numerical simulations of inplane and outofplane ultimate carrying capacities of Q460 steel welded I and boxsection columns loaded with an eccentricity in a plane of symmetry were conducted by Shen and Yang [20], Shen and Liu [21], Shen [22], and Shen and Zhao [23], respectively. The beamcolumns were selected to cover a wide range of plate widthtothickness ratios ranging from 40 to 70 and a wide range of member slenderness ratios from 20 to 120. The effects of plate widthtothickness ratio, column slenderness ratio, and eccentricity on the ultimate carrying capacities were investigated, and the calculation formulas for the ultimate strengths were proposed. However, all these studies focused on the HSS beamcolumns under uniaxial bending, and more attention should be paid to those under biaxial bending.
In addition, no strength formulas are provided specially for HSS welded beamcolumns in some existing standards GB 500172017 [24], ANSI/AISC 36010 [25], EN 199311 [26], and EN 199315 [27]. These standards are all based on the research results of mild steels and are also used for HSSs, but it is unknown whether they are applicable for HSS members. Furthermore, when the strength of thinwalled box beamcolumns under biaxial bending is calculated in accordance with the above standards, a cumbersome calculation has to be made because of the effective width method, which is used in these standards. Hence, a simple method needs to be provided.
In this study, the localoverall interaction buckling behavior of Q460 steel welded thinwalled box beamcolumns under biaxial bending is analyzed by the commercial software ANSYS 8.0 [28, 29]. The numerical results, as well as available experimental results of different steel grades, are compared with ANSI/AISC 36010 to check whether ANSI/AISC 36010 is suitable for such types of beamcolumns. Simple formulas are presented for predicting the localoverall interactive buckling strength of Q460 steel beamcolumns subjected to axial compression and biaxial bending moments. The proposed formulas are verified by comparing them with the available experimental results.
2. Finite Element Analysis
2.1. An FE Model and Its Verification
A pinended beamcolumn, as shown in Figure 1, with a welded thinwalled box cross section, fabricated from Q460, subjected to combined axial compression force, , and biaxial bending moments, M_{x} and M_{y}, is investigated in this paper. Here, the bending moments around the x and yaxes, M_{x} and M_{y}, are taken as M_{x} = Pe_{y} and M_{y} = Pe_{x}, where e_{x} and e_{y} are the eccentricities in the x and y directions, respectively. Also given in Figure 1(b) are the dimensions of the welded box section, in which b and h are the widths of plate elements and t is the thickness of the plate elements.
(a)
(b)
A geometrical and material nonlinear FE model was established by Shen [9, 22] for HSS box columns loaded with axial compression and uniaxial eccentric compression, respectively, by using the ANSYS program. In this study, the FE model is modified for HSS box columns loaded with biaxial eccentric compression.
As described by Shen [9, 22], the effects of initial geometrical imperfections and residual stresses are taken into account. The initial geometrical imperfections include both overall member imperfections and local plate imperfections. The member imperfections are assumed to bend in a half sine wave about the major and minor axes, respectively, as follows:where l is the length of the member.
In equations (1a) and (1b), the maximum value of outofstraightness of a member is taken equal to l/1000 as specified by GB 500172003 [30] and ANSI/AISC 36010 [25].
The plate imperfections are assumed aswhere , , and m is the buckling halfwave number along the longitudinal axis of the member.
The magnitudes of plate imperfections were taken as b/100 and b/1000 (b is the plate width of a square box section) [8], and it showed that the magnitude of plate imperfections had a significant effect on the compressive strength of welded box columns. No uniform value was specified for the maximum of plate imperfections, and it was taken as 0.01 [4] and b/1000 [3], but it was demonstrated that the value of b/1000 seemed to be reasonable [3] when residual stresses were taken into consideration. Hence, for a square box section, . But for a rectangular one, according to the deformation compatibility condition, the angle between the two adjacent plates remained unchanged before and after the plate buckling, so and herein.
Residual stress measurements of welded square and rectangular box sections were made by Usami and Fukumoto [11] and Pavlovčič et al. [6], respectively. According to their measurement results, a simplified pattern of residual stress distribution, as shown in Figure 2, is adopted, where tensile residual stresses are positive values and compressive stresses are negative values. An SM58 steel plate was used [11]; its nominal yield strength was 460 MPa, and its measured yield strength was 568 MPa. The maximum value of tensile stresses was about 80% of the measured yield strength, i.e., 0.8 × 568 = 454 MPa, which is very close to the nominal yield strength of 460 MPa. The value of compressive residual stresses decreased along with the increase of the widthtothickness ratio of the plate [6, 11]. The distribution width of tensile residual stresses was about 3t [9]. As a consequence, the maximum value of tensile residual stresses is taken as the nominal yield stress of steel and those of compressive residual stresses and , along the flange and web, respectively, may be obtained according to equilibrium conditions.
The Shell181 element is employed to take account of plate imperfections. Shell181 is suitable for analyzing thin to moderately thick shell structures. It is a 4node element with six degrees of freedom at each node: translations in the x, y, and z directions and rotations about the x, y, and zaxes. It is well suited for linear, largerotation, and/or largestrain nonlinear applications. Also, the effects of plasticity of steel, member imperfections, and residual stresses can be considered by using Shell181. As a result, the Shell181 element satisfies the requirements of this study.
The experiment results [31] showed that stressstrain curves for Q460 steel, as shown in Figure 3, were similar to those for mild steel and also had a yield plateau. The material models with and without strain hardening had almost no influence on the compressive strength of welded box columns [8]. The similar finding was also noted by Pircher et al. [7] in the simulation of the buckling of thinwalled steel box tubes. So, an ideal elasticplastic material model is assumed to represent the stressstrain relationship for Q460 steel. The nominal material properties: yield stress f_{y} = 460 N/mm^{2}, Young’s modulus E = 206000 N/mm^{2}, Poisson’s ratio = 0.3, and tangent modulus E_{t} = 0, are used.
The effect of geometrical imperfections is taken into account by direct modeling according to the following steps: Firstly, the number of nodes required in the FE model and the order in which the nodes should be generated are determined. Secondly, all nodes are generated according to the x, y, and z coordinates considering geometrical imperfections. Thirdly, the element attributes, including the element type, material model, and real constant, have been set, and shell elements are automatically generated within each area defined by four nodes. For any plate element in a box column, element sizes are b/8 or h/8 along its width and h/4 along its length, and thus, a finite element grid is b/8 × h/4 or h/8 × h/4. In this paper, the minimum grid is approximately 20 mm × 40 mm and the maximum one 50 mm × 100 mm, which satisfies the accuracy requirement of the numerical simulation. Figure 4 shows the initial geometrical imperfections generated by direct modeling, in which the initial deflections are magnified 5,000 times only to be observed clearly.
(a)
(b)
Residual stresses are treated as initial stresses that are loads and must be applied at the integration points of Shell181 elements. For a Shell181 element, there are four integration points in its plane, together with 5 integration points in the thickness direction to take account of nonlinear bending properties, and thus, the total number of integration points is 20. In order to input initial stresses in the FE model, two steps must be performed: Firstly, an initial stress file needs to be written based on the integration point locations. Secondly, the initial stress file must be read by using the ISFILE command only at the first substep of the first load step, as specified in the ANSYS program. In this study, the same residual stress distributions are assumed, along the thickness direction, because of thin plate elements.
The loading and boundary conditions of the beamcolumns studied in this paper are different greatly from those of the beamcolumns [22]. Figure 5 displays the loads and boundary conditions in the FE model. As shown in Figure 5, a 30 mm thick end plate is attached to both ends of a beamcolumn [32], and an eccentric load at every end of the member is equivalent to combined axial compressive force and biaxial bending moments. These forces, as well as the boundary conditions, are applied at the center points of both end plates [32]. At one end, z = 0, besides the bending moments M_{x} and M_{y}, the three translations U_{x}, U_{y}, and U_{z} in the x, y, and z directions and the rotation Rot_{z} about the zaxis are prevented, i.e., U_{x} = U_{y} = U_{z} = Rot_{z} = 0. At the other end, z = l, besides , M_{x}, and M_{y}, the boundary conditions satisfy U_{x} = U_{y} = Rot_{z} = 0. In addition, both ends of the member are restrained against the rotation about the longitudinal axis by applying U_{x} = U_{y} = 0 at the four corner points of the end plate.
(a)
(b)
(c)
The available experimental data are used to verify the FE model. 35 experimental specimens fabricated from two types of HSSs having f_{y} = 460 MPa and 690 MPa, of which 21 were centrally loaded and 14 eccentrically loaded in a symmetry plane by Usami and Fukumoto [11, 12], were simulated by Shen [22], and the ratios of the numerical and experimental ultimate strength change from 0.901 to 1.084 with an average value of 1.012 and a standard deviation of 4.53%. Here, the FE model [22] was modified for biaxial bending, but no published experimental information could be found on the beamcolumns made from HSSs and subjected to biaxial bending. Hence, in this section, 28 experimental specimens tested by Richard Liew et al. [16], fabricated from mild steels with the nominal yield strength of 353 MPa, 268 MPa, and 293 MPa, were simulated to verify the effectiveness of the modified FE model. Among these specimens, 10 were loaded with uniaxial eccentricity and labeled with a single letter X or Y, as shown in Table 1, and 18 were loaded with biaxial eccentricity and marked with two letters XY. The second field S of the specimen number indicates a square box section, while R indicates a rectangular box section. The third field refers to the web depthtothickness ratio, and the last field denotes the column slenderness ratio about the minor axis. Also listed in Table 1 are the main parameters of the specimens, including the widthtothickness ratio b/t; depthtothickness ratio h/t; maximum initial deflections δ_{0x}/L and δ_{0y}/L about x and yaxes, respectively; eccentricities e_{x} and e_{y} in x and y directions, respectively; and measured yield stress f_{y}. All these parameters, as well as the material’s Young’s modulus, Poisson’s ratio, and residual stress patterns measured by Richard Liew et al. [16], are used.

A comparison between the experimental and numerical results is also presented in Table 1, where P_{exp} and P_{fem} are the ultimate strengths obtained from the experimental result and numerical simulation, respectively. The ratios of the experimental and numerical ultimate strength, P_{exp}/P_{fem}, vary within a range of 0.86 to 1.12 with an average value of 0.97 and a standard deviation of 6.8%. It can be seen that the numerical results are in good agreement with the experimental ones, and the FE model established in this study can precisely predict the localoverall interactive buckling strength of welded box beamcolumns under both uniaxial bending and biaxial bending.
2.2. Buckling Behavior of Q460 Steel Welded ThinWalled Box BeamColumns
By using the FE model mentioned above, a great amount of numerical simulation work with different parameters has been performed on the localoverall interactive buckling of Q460 steel thinwalled box columns under the biaxial eccentric loading condition. The parameters include the crosssectional shape, plate widthtothickness ratio, column slenderness ratio, and applied eccentricities e_{x} and e_{y}.
In the calculation, the cross sections include square and rectangular box sections and the depthtowidth ratio of the rectangular boxes is 1.25. The thickness of the plate elements, t, is taken as 4 mm and kept constant.
Plate widthtothickness ratios, together with column slenderness ratios, are selected elaborately to ensure that the localoverall interaction buckling occurs. In accordance with the Chinese Code GB 500172003 [30], the limiting widthtothickness ratio is for compression elements of members subjected to axial compression and bending. This limiting value is in ANSI/AISC 36010 [25]. Both of them are basically the same. When the widthtothickness ratio is greater than the limiting value, the localoverall interaction buckling will occur. Therefore, the widthtothickness ratio, b/t, is taken as 40, 50, 60, 70, and 80, respectively, whether a square or a rectangular box cross section. If it is a rectangular box, the corresponding depthtothickness ratio, h/t, will be 50, 62.5, 75, 87.5, and 100, respectively.
To obtain the predicated failure mode of localoverall interaction buckling, the beamcolumns with medium slenderness ratios should be selected. Therefore, the member slenderness ratio about the major axis, λ_{x} (λ_{x} = l/i_{x}, where l is the effective length of a member and i_{x} is the gyration radius about the major axis, i.e., xaxis), is taken equal to 40, 60, 80, and 100, respectively.
Eccentricities e_{x} and e_{y} are taken as two of 10 mm, 15 mm, 20 mm, 25 mm, 30 mm, 35 mm, 40 mm, and 50 mm, respectively, satisfying e_{x} = e_{y} or e_{y} − e_{x} = 5 mm in order to save materials.
The beamcolumns are marked with numbers, such as S4040 (10, 10), in which the first S indicates a square box section, the second is the widthtothickness ratio, the third is the column slenderness ratio about the xaxis, and the figures in brackets are the eccentricities e_{x} and e_{y}, respectively. If it is a rectangular box section, S is replaced by R, and another figure following the second part is added representing the depthtothickness ratio.
2.2.1. LoadDisplacement or LoadRotation Curves
Figure 6 shows the flexural deformation of a beamcolumn, where μ_{max} and are the maximum deflections in the x and y directions, respectively, and the section 11 is the cross section of the member at midlength. Points A and B are the midpoints of the upper flange and left web of the cross section 11, respectively. Owing to the interaction between local buckling and biaxial bending, the maximum deflections in the x and y directions of the points on the crosssection 11 are different from each other. In this study, μ_{max} is taken as the deflection of the point B and as that of the point A.
(a)
(b)
(c)
There is a large difference in values of the displacement between a long and a short column, so their loaddisplacement and loadrotation curves are plotted separately. Figures 7 and 8 show the typical ones for the member S4040 (40, 40) that represents a short column and the ones for the member S50100 (30, 35) that represents a long column, respectively. In Figures 7 and 8, the curves of the axial compressive force, , versus the axial compressive displacement, Δ, the maximum deflection in the x direction, μ_{max}, the maximum deflection in the y direction, , and the maximum rotation about the zaxis, θ_{max}, which is taken as the rotation of the point A as above, are shown, respectively. It can be seen from Figures 7 and 8 that when reaching the ultimate strengths (440.36 kN and 338.45 kN), for the member S4040 (40, 40), Δ = 7.21 mm, μ_{max} = 12.24 mm, = 12.86 mm, and θ_{max} = 0.005 rad, while for the member S50100 (30, 35), Δ = 12.71 mm, μ_{max} = 64.92 mm, = 84.93 mm, and θ_{max} = −0.001 rad. From these values, a conclusion can be drawn that the bending deflections about the x and yaxes and axial compressive displacement are large and the rotation about the zaxis is very small for both the short and long boxsection beamcolumns subjected to biaxial bending. Such a small twisting effect is due to the fact that a box section has large torsional rigidity.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
2.2.2. Von Mises Stress Distribution of the Cross Section at Midlength
A biaxial bending beamcolumn of this type might be in the elastic region when it reaches the maximum strength. It can be determined by viewing von Mises stress distribution. A typical member R5062.560 with b/t = 50, h/t = 62.5, and λ_{x} = 60 is selected to exhibit the von Mises stress distribution of the crosssection 11, as shown in Figure 6, where the maximum bending moments occur.
To make it easier to describe the von Mises stress distribution, the coordinate system, as shown in Figure 6, is defined again. The topleft corner point of the crosssection 11 is taken as the coordinate origin, the xaxis is horizontal with values increasing from left to right, and the yaxis is vertical with values increasing from top to bottom. The von Mises stress distributions at the ultimate loads are plotted in Figure 9 under different eccentricities. In Figure 9, all eccentricities e_{x} and e_{y} are in mm. A little difference in von Mises stress distributions exists in the bottom flange and right web under eight different eccentricities, but obvious difference is observed in the top flange and left web. The maximum von Mises stress of 452.77 MPa occurs at the corner point between the top flange and the left web when e_{x} = 50 mm and e_{y} = 50 mm. All these stresses are smaller than the steel yield stress of 460 MPa, and the beamcolumn is still in the elastic range and does not enter the elasticplastic range under any eccentricity condition.
(a)
(b)
(c)
(d)
2.2.3. Interaction between Bending Moment and Axial Force
It is very important for a beamcolumn subjected to biaxial bending to find the relationship equation between the axial load and end bending moments. The scatter diagrams between P_{u}/P_{y} and M_{ux}/M_{xy} (or M_{uy}/M_{yy}) of Q460 steel thinwalled boxsection beamcolumns with five different widthtothickness ratios are essentially the same, so Figures 10 and 11 only show the ones for beamcolumns with square box sections and b/t = 40, 60, and 80 and for those with rectangular box sections and b/t = 60 and 80, respectively, where , , , M_{xy} = W_{x}f_{y}, and M_{yy} = W_{y}f_{y} (A is the gross area and W_{x} and W_{y} are the gross section moduli about x and yaxes, respectively). In Figures 10 and 11, the scatter plots are the finite element results, and also, the straight lines that predict the changing trends of all the coordinate points are added. It is seen in Figures 10 and 11 that the straight lines can fit very well the overall trend of the scatter plots between and M_{ux}/M_{xy} (or M_{uy}/M_{yy}).
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
(c)
(d)
3. Evaluation of the American National Standard
3.1. Comparison between EN 199311 and ANSI/AISC 36010
In the Eurocode EN 199311 [26], for a member under combined axial compression and biaxial bending, the following conditions should be satisfied:where N, M_{x}, and M_{y} are the axial compression force and the maximum moments about the x and yaxes, respectively; ΔM_{x} and ΔM_{y} are the moments due to the shift of the centroidal axis for class 4 sections; N_{R} is the resistance to axial compression force; M_{x,R} and M_{y,R} are the resistances to bending moments about the x and yaxes, respectively; χ_{x} and χ_{y} are the reduction factors due to flexural buckling; χ_{LT} is the reduction factor due to lateral torsional buckling; and k_{xx}, k_{xy}, k_{yx}, and k_{yy} are the interaction factors.
This type of cross section studied in this paper is called class 4 section in EN 199311 [26]. For class 4 cross sections, these conditions are expressed aswhere , , and are the effective crosssectional area and effective section moduli about the x and yaxes, respectively, and and are the eccentricities with respect to the neutral axes.
According to the American National Standard ANSI/AISC 36010 [25], for members subjected to axial force and biaxial end moments, the following conditions should be met:where is the required compressive strength; is the nominal compressive strength; M_{rx} and M_{ry} are the required flexural strengths about the x and yaxes, respectively; and M_{nx} and M_{ny} are the nominal flexural strengths about the x and yaxes, respectively.
In ANSI/AISC 36010 [25], the cross section studied herein is called the slender section. For slender sections, , M_{nx}, and M_{ny} are determined on the basis of the effective width concept (see Section 3.2).
Comparing the two standards, the following can been found: (1) Linear interaction equations between the compression force and biaxial bending moments are used in both the standards, but two formulas in EN 199311 and only one formula in ANSI/AISC 36010. (2) The secondorder effect of δ is taken into account. It is considered in the interaction factors k_{xx}, k_{xy}, k_{yx}, and k_{yy} in EN 199311, while it is included in the required flexural strengths M_{rx} and M_{ry} in ANSI/AISC 36010, but the calculations of k_{xx}, k_{xy}, k_{yx}, and k_{yy} are very complicated, while those of M_{rx} and M_{ry} are simple. (3) For a beamcolumn with a thinwalled section, both of the specifications adopt the effective width concept. In short, compared to EN 199311, the calculation according to ANSI/AISC 36010 becomes much simpler. Hence, a comparison between ANSI/AISC 36010 and the available experimental results and that between ANSI/AISC 36010 and the numerical results will be made in the following sections to evaluate the applicability of ANSI/AISC 36010 to HSS thinwalled beamcolumns.
3.2. Comparison between ANSI/AISC 36010 and the Available Experimental Results
When making a calculation by employing equations (5a) and (5b), the values of , M_{nx}, M_{ny}, M_{rx}, and M_{ry} must be determined first.
For slender sections, the nominal compressive strength should be taken aswhere is the flexural buckling stress and determined bywhere K is the effective length factor; L is the laterally unbraced length of the member; r is the governing radius of gyration; is the steel yield stress; is the elastic critical buckling stress; and Q is the reduction factor, for slender stiffened elements, Q = Q_{a}.
Q_{a} is defined aswhere A_{e} is the effective crosssectional area based on effective width b_{e}.
The effective width b_{e} is expressed as follows:where
It requires iteration to calculate the strength of a column according to equations (6)–(10). In order to simplify the calculation, f_{1} is taken equal to as specified in ANSI/AISC 36010 [25]. This will result in a slightly conservative estimate of available column strengths.
For slender sections, the nominal flexural strengths M_{nx} and M_{ny} shall be determined as follows:where S_{ex} and S_{ey} are the effective section moduli about the x and yaxes and determined with the effective width b_{e}, respectively.
In order to consider the secondorder effect of δ, M_{rx} and M_{ry} are calculated using the following equations:where C_{mx} and C_{my} are the equivalent moment factors and P_{Ex} and P_{Ey} are Euler’s critical loads.
Equations (5a) and (5b) are verified by the available experimental results (for slender box sections), as shown in Table 2. These experimental results are taken from different literatures with the measured steel yield stress of 268∼741 MPa. By substituting these experimental values into the lefthand side of equations (5a) and (5b), the calculation results marked with “(1)” in Table 2 are obtained and are closer to 1 except for the three specimens with λ_{x} = 85 [16] and the two specimens with f_{y} = 741 MPa [12]. The average values of “(1)” in Table 2 are about 1.58 and 1.26 for the specimens with λ_{x} = 85 (f_{y} = 353 MPa) and those with f_{y} = 741 MPa, respectively. Comparisons show that the standard ANSI/AISC 36010 agrees, in most instances, very well with the experimental results. The ANSI/AISC 36010 standard provides a better evaluation of the strengths of beamcolumns with medium slenderness ratios and f_{y} = 248 MPa∼373 MPa. It underestimates the strengths of beamcolumns with f_{y} = 741 MPa but overestimates the capacities of members with f_{y} = 568 MPa.
 
; ; f_{y} is the measured yield stress. 
3.3. Comparison between ANSI/AISC 36010 and the Numerical Results
The standard ANSI/AISC 36010 is also compared with the numerical results, as shown in Figure 12. S40 represents a series of square box beamcolumns with b/t = 40, and R4050 represents another series of rectangular ones with b/t = 40 and h/t = 50. The meaning of other symbols (S50, S60, …) is the same as that of the two members. The maximum value of the lefthand side of equation (5) is 1.79, the minimum one is 0.83, the average one is 1.09, and the standard deviation is 14.8%. For Q460 steel, the same phenomenon, as observed in the comparison between ANSI/AISC 36010 and the available experimental results, is found in Figure 12, and ANSI/AISC 36010 severely underestimates the strengths of beamcolumns with large slenderness ratios and small widthtothickness ratios, such as S40 with λ_{x} = 80 and 100, S60 with λ_{x} = 100, R4050 with λ_{x} = 80 and 100, and R5062.5 with λ_{x} = 100, especially for rectangular box sections. Except for these members, ANSI/AISC 36010 gives good estimation of the strength of Q460 steel thinwalled box beamcolumns subjected to axial force and biaxial end moments. It seems that ANSI/AISC 36010 is not suitable for the rectangular box sections with a large depthtowidth ratio owing to the symmetry of the calculation formulas for the x and yaxes.
4. Proposal Strength Formulas for Q460 Steel ThinWalled Box BeamColumns Subjected to Axial Force and Biaxial End Moments
4.1. Proposed Strength Formulas
The ANSI/AISC 36010 standard is also suitable for Q460 steel thinwalled box beamcolumns under biaxial bending over a wide range of column slenderness ratios. However, a complex calculation must be done to acquire the effective section properties to (1) calculate the effective width, (2) determine the centroidal axis of the effective section, and (3) calculate the effective section properties. Moreover, the available experimental results, together with the finite element results, show that ANSI/AISC 36010 severely underestimates the strengths of beamcolumns with large slenderness ratios (f_{y} = 353 MPa and 460 MPa) and slightly underestimates those of beamcolumns with a yield stress of 741 MPa. As a result, a simple method, which uses the gross section properties, is proposed.
For a beamcolumn, the interaction formula between the axial compression load and the bending moment gives a theoretical basis to develop the strength formulas. As mentioned before, a linear relationship between P_{u}/P_{y} and M_{ux}/M_{xy} or between P_{u}/P_{y} and M_{uy}/M_{yy} is found, in most cases, for Q460 steel box crosssectional biaxial bending beamcolumns with slender plates. In addition, Salem et al. [17] investigated the ultimate strength of biaxially loaded slender Isection beamcolumns made from mild steel and concluded that the linear interaction relationship could predict satisfactorily the strength of beamcolumns with medium slenderness ratios. Therefore, the following linear interaction equation that is used in EN 199311 [26] and NAS AISI 2007 [33] is adopted:where is the axial force; is the allowable pure axial force; M_{xc} and M_{yc} are the allowable pure bending moments about the x and yaxes, respectively; and α_{x} and α_{y} are the amplification factors for bending moments.
When equation (13) is used for class 4 cross sections in EN 199315 [27], it is changed into equations (4a) and (4b).
When equation (13) is used for coldformed sections in NAS AISI 2007 [33], it is modified as follows:where and are the equivalent moment factors about the x and yaxes, respectively; is the smaller one of and ; and are the effective section moduli about the x and yaxes, respectively; and and are Euler’s critical loads.
Also, the effective section properties are adopted in equation (14), thus leading to a complex calculation. Hence, a simple method is proposed in the following section.
Being in the elastic range is found when the Q460 steel beamcolumns reach the ultimate carrying capacities from the von Mises stress distributions in Figure 9. The same phenomenon was also discovered by Shen [22]. A simple formula was proposed by Shen [22] to evaluate the strength of HSS rectangular crosssectional beamcolumns with slender webs under uniaxial bending. The formula is expressed as follows:where is the strength reduction factor of a centrally loaded member about the xaxis; is a parameter, (where 1.1 is a partial factor for resistance); is the modification factor for the steel yield strength; and is the steel yield strength.
In addition, for boxsection beamcolumns under biaxial bending, as mentioned in Section 2.2.1, the bending deflections about the x and yaxes and axial compressive displacement are two major deformations and the rotation about the zaxis is very small. For this reason, using equation (14) as the reference, equation (15) can be modified for the ultimate strength of Q460 steel thinwalled boxsection columns subjected to axial compressive force and major and minoraxis bending moments. The modified equation is expressed as follows:where is the strength reduction factor of a centrally loaded member about the yaxis and is a parameter, .
The ultimate carrying capacity of columns with nonslender sections and made from Q345, Q390, Q420, and SM58 steels (the nominal yield stresses are 345 MPa, 390 MPa, 420 MPa, and 460 MPa, respectively) was simulated to determine which column curve should be used when calculating the values of and by Shen [9]. The results show that the curve a in GB 500172003 [30] should be adopted. Degée et al. [3] also suggested that the curve a in EN 199311 [26] should be employed when the local and global interaction buckling of welded boxsection compression members fabricated from S355, S460, and S690 steels (the nominal yield stresses are 355 MPa, 460 MPa, and 690 MPa, respectively) was investigated. There is a little difference in the value of the equivalent imperfection between GB 500172003 and EN 199311 when determining the stability reduction factor, but the values of three design column curves a, b, and c in GB 500172003 are very close to those of three column curves a, b, and c in EN 199311. As a result, in equation (16), and are obtained according to and , respectively, and the column curve a recommended in the Chinese Code GB 500172003 [30].
Besides, in equation (16), for a beamcolumn with equal end flexural moments; N is taken as ; M_{x} as M_{ux}, i.e., ; and M_{y} as M_{uy}, i.e., ; and then the numerical result, , is substituted into equation (16). After finishing all these, the modified factor, , can be obtained by dividing the lefthand side of equation (16) by the steel yield strength, .
The ultimate strength of Q460 steel welded thinwalled box crosssectional beamcolumns is related to the slenderness ratios, λ_{x} and λ_{y}, widthtothickness ratio, b/t, and depthtothickness ratio, h/t. For rectangular box sections, in this paper, the depthtowidth ratio is a constant of 1.25; hence, h/t is related to b/t and, similarly, λ_{y} to λ_{x}. Therefore, the modification factor, , may be treated as a function of only two variables b/t and λ_{x}. The study shows that, in most cases, the nondimensional ultimate carrying capacity, ( is the ultimate carrying capacity), is nearly linear to each of λ_{x} and b/t, respectively. By using a curvefitting technique, the following equation can be obtained:
Substituting the numerical results and equation (17) into equation (16) and introducing the symbols , , and , the ratios of the left side to the right side of equation (16) are obtained and plotted in Figure 13. In Figure 13, the ratio values change over a range of 0.818 to 1.239, with an average value of 1.015 and a standard deviation of 7.92%. Hence, equation (16) coincides quite well with the numerical results.
4.2. Verification of the Proposed Formulas
The available experimental results, as listed in Table 2, are also used to verify the proposed formulas. In order to compare them easily, equation (16) is changed into the following equation:
Similarly, substituting the experimental results and equation (17) into the lefthand side of equation (18), the calculation results labeled with “(2)” in Table 2 are obtained. These values are very close to 1 except for four beamcolumns with larger widthtothickness ratios (b/t = 80 and 85) and fabricated from the mild steels with f_{y} = 248.6 MPa and f_{y} = 268 MPa. For the four members, the average value of “(2)” in Table 2 is about 1.40 and equation (18) gives a low estimate of their strengths. In general, the simple formulas proposed herein give a reasonably accurate estimate of the strength of the beamcolumns over a wide range of the steel yield stress, and for the most part, they are applicable to the thinwalled box beamcolumns made from both mild and highstrength steels. Meanwhile, a standard deviation of 15.7% listed in Table 2, compared to 17.2%, shows that the simple formulas achieve the same level of accuracy as the American Standard ANSI/AISC 36010. For the HSS specimens with f_{y} = 741 MPa, the simple formulas, based on the numerical results, provide a good estimate of their strengths; for those with f_{y} = 568 MPa, they give a slight overestimate; however, all these specimens were loaded with uniaxial eccentricity, so more experiments, on the HSS beamcolumns under biaxial bending, are needed to be conducted to verify the proposed formulas.
5. Conclusions
The behavior and strength of Q460 steel welded thinwalled boxsection beamcolumns under combined compression force and biaxial bending are analyzed by an established FE model. The American Specification ANSI/AISC 36010 is verified with the numerical results and the available experimental results. Simple and effective formulas are proposed for predicting the maximum strength of the beamcolumns. Also, the proposed formulas are compared with the available experimental results. Some conclusions can be drawn from this paper:(1)The established FE model can simulate the localoverall interaction buckling behaviors of the beamcolumns subjected to combined axial compression force and biaxial bending.(2)It is observed that, in most cases, the von Mises stresses are much smaller than the steel yield strength when the beamcolumns reach the ultimate loads.(3)The correlation curves between the axial force and the flexural moment M_{ux}/M_{xy} (or M_{uy}/M_{yy}) of Q460 steel thinwalled boxsection beamcolumns under compression and biaxial bending moments are nearly linear.(4)The code ANSI/AISC 36010 agrees very well with both the numerical and experimental results except for two types of beamcolumns, one with large slenderness ratios and the other made from the HSS with a yield stress of 741 MPa. For the former, ANSI/AISC 36010 seriously underestimates the strengths, and for the latter, it slightly underestimates.(5)After having been modified, the formula proposed for the localoverall interactive buckling strength of HSS box beamcolumns under uniaxial bending can be used for Q460 steel box beamcolumns under biaxial bending. The proposed modified formulas are in good agreement with the numerical and experimental results over a wide range of the steel yield stress. The maximum of nominal yield stress in the available tests is 690 MPa, and hence, f_{y} ≤ 690 MPa has been adopted as a limit of application of the proposed formulas. h/b = 1.0 or 1.25 in this study, and h/b or b/h = 1.0∼1.33 in the available experiments [2, 6, 11, 12, 16]. In practical projects, for economic consideration, h/b or b/h of box sections should be close to or slightly greater than 1 to avoid that the carrying capacity in one direction is much higher than that in the other direction. Therefore, the modification factor , which should be a function of b/t, h/t, λ_{x}, and λ_{y}, can be simplified as a function of b/t and λ_{x}, as given in equation (17). But the available experimental data on HSSs are very limited, and more experiments are needed to be conducted to verify these proposed formulas.
Data Availability
The data used to support the findings of this study are available from the author upon request.
Conflicts of Interest
The author declares that there are no conflicts of interest.
Acknowledgments
The research reported in this paper was supported by National Science Basic Research Plan in Shaanxi Province, China (Grant no. 2018JM5079), and the Steel Structure Innovation Team of Xi’an University of Architecture and Technology, China. Their financial support is highly appreciated.
References
 L. Bai and M. A. Wadee, “Slenderness effects in thinwalled Isection struts susceptible to local–global mode interaction,” Engineering Structures, vol. 124, pp. 128–141, 2016. View at: Publisher Site  Google Scholar
 S. P. Chiew, S. L. Lee, and N. E. Shanmugam, “Experimental study of thinwalled steel box columns,” Journal of Structural Engineering, vol. 113, no. 10, pp. 2208–2220, 1987. View at: Publisher Site  Google Scholar
 H. Degée, A. Detzel, and U. Kuhlmann, “Interaction of global and local buckling in welded RHS compression members,” Journal of Constructional Steel Research, vol. 64, no. 78, pp. 755–765, 2008. View at: Publisher Site  Google Scholar
 Y. B. Kwon, N. G. Kim, and G. J. Hancock, “Compression tests of welded section columns undergoing buckling interaction,” Journal of Constructional Steel Research, vol. 63, no. 12, pp. 1590–1602, 2007. View at: Publisher Site  Google Scholar
 E. L. Liu and M. A. Wadee, “Geometric factors affecting Isection struts experiencing local and strongaxis global buckling mode interaction,” ThinWalled Structures, vol. 109, pp. 319–331, 2016. View at: Publisher Site  Google Scholar
 L. Pavlovčič, B. Froschmeier, U. Kuhlmann, and D. Beg, “Finite element simulation of slender thinwalled box columns by implementing real initial conditions,” Advances in Engineering Software, vol. 44, no. 1, pp. 63–74, 2012. View at: Google Scholar
 M. Pircher, M. D. O’Shea, and R. Q. Bridge, “The influence of the fabrication process on the buckling of thinwalled steel box sections,” ThinWalled Structures, vol. 40, no. 2, pp. 109–123, 2002. View at: Publisher Site  Google Scholar
 H. X. Shen, “Ultimate capacity of welded box section columns with slender plate elements,” Steel & Composite Structures, vol. 13, no. 1, pp. 15–33, 2012. View at: Publisher Site  Google Scholar
 H. X. Shen, “On the direct strength and effective yield strength method design of medium and high strength steel welded square section columns with slender plate elements,” Steel and Composite Structures, vol. 17, no. 4, pp. 497–516, 2014. View at: Publisher Site  Google Scholar
 J. J. Shen, M. A. Wadee, and A. J. Sadowski, “Interactive buckling in long thinwalled rectangular hollow section struts,” International Journal of NonLinear Mechanics, vol. 89, pp. 43–58, 2017. View at: Publisher Site  Google Scholar
 T. Usami and Y. Fukumoto, “Welded box compression members,” Journal of Structural Engineering, vol. 110, no. 10, pp. 2457–2470, 1984. View at: Publisher Site  Google Scholar
 T. Usami and Y. Fukumoto, “Local and overall buckling of welded box columns,” Journal of Structural Division, vol. 108, no. 3, pp. 525–542, 1982. View at: Google Scholar
 J. J. Shen and M. A. Wadee, “Length effects on interactive buckling in thinwalled rectangular hollow section struts,” ThinWalled Structures, vol. 128, pp. 152–170, 2018. View at: Publisher Site  Google Scholar
 J. J. Shen and M. A. Wadee, “Imperfection sensitivity of thinwalled rectangular hollow section struts susceptible to interactive buckling,” International Journal of NonLinear Mechanics, vol. 99, pp. 112–130, 2018. View at: Publisher Site  Google Scholar
 S. L. Lee, N. E. Shanmugam, and S. P. Chiew, “Thinwalled box columns under arbitrary end loads,” Journal of Structural Engineering, vol. 114, no. 6, pp. 1390–1402, 1988. View at: Publisher Site  Google Scholar
 J. Y. Richard Liew, N. E. Shanmugam, and S. L. Lee, “Behavior of thinwalled steel box columns under biaxial loading,” Journal of Structural Engineering, vol. 115, no. 12, pp. 3076–3094, 1989. View at: Publisher Site  Google Scholar
 A. H. Salem, M. El Aghoury, F. F. El Dib, and M. T. Hanna, “Strength of biaxially loaded slender Isection beamcolumns,” Canadian Journal of Civil Engineering, vol. 34, no. 2, pp. 219–227, 2007. View at: Publisher Site  Google Scholar
 N. E. Shanmugam, S. P. Chiew, and S. L. Lee, “Strength of thinwalled square steel box columns,” Journal of Structural Engineering, vol. 113, no. 4, pp. 818–831, 1987. View at: Publisher Site  Google Scholar
 N. E. Shanmugam, J. Y. Richard Liew, and S. L. Lee, “Thinwalled steel box columns under biaxial loading,” Journal of Structural Engineering, vol. 115, no. 11, pp. 2706–2726, 1989. View at: Publisher Site  Google Scholar
 H. X. Shen and C. H. Yang, “Analysis of localoverall interaction buckling of highstrength steel welded Isection compressionbending members,” Building Structure, vol. 43, no. 22, pp. 33–38, 2013, in Chinese. View at: Google Scholar
 H. X. Shen and X. Liu, “Analysis on the ultimate carrying capacities of localoverall interaction of highstrength steel welded squarebox section columns loaded with eccentricity,” Building Structure, vol. 44, no. 4, pp. 35–38, 2014, in Chinese. View at: Google Scholar
 H. X. Shen, “Behavior of highstrength steel welded rectangular section beam–columns with slender webs,” ThinWalled Structures, vol. 88, pp. 16–27, 2015. View at: Publisher Site  Google Scholar
 H. X. Shen and K. X. Zhao, “Finite element analysis of local and overall flexuraltorsional interactive buckling of Q460 high strength steel welded Isection beamcolumns,” Progress in Steel Building Structures, vol. 17, no. 4, pp. 1–9, 2015, in Chinese. View at: Google Scholar
 GB 500172017, Standard for Design of Steel Structures, Guobiao Standards, Beijing, China, 2017.
 ANSI/AISC 36010, Specification for Structural Steel Buildings, AISC, Chicago, IL, USA, 2010.
 EN 199311, Eurocode 3: Design of Steel StructuresPart 11: General Rules and Rules for Buildings, European Committee for Standardization, Brussels, Belgium, 2005.
 EN 199315, Eurocode 3Design of Steel StructuresPart 15: Plated Structural Elements, European Committee for Standardization, Brussels, Belgium, 2007.
 ANSYS Inc, ANSYS Theory Reference, Electronic Release 8.0, SAS IP, Inc., Cary, NC, USA, 1998.
 ANSYS Inc, ANSYS Element Reference, Electronic Release 8.0, SAS IP, Inc., Cary, NC, USA, 1998.
 GB 500172003, Code for Design of Steel Structures, Guobiao Standards, Beijing, China, 2003.
 H. Y. Ban, G. Shi, Y. J. Shi, and Y. Q. Wang, “Overall buckling behavior of 460 MPa high strength steel columns: experimental investigation and design method,” Journal of Constructional Steel Research, vol. 74, pp. 140–150, 2012. View at: Publisher Site  Google Scholar
 Y. Liu and L. B. Hui, “Finite element study of steel single angle beam–columns,” Engineering Structures, vol. 32, no. 8, pp. 2087–2095, 2010. View at: Publisher Site  Google Scholar
 NAS AISI, North American Specifications for the Design of ColdFormed Steel Structural Members, AISI, Washington, DC, USA, 2007.
Copyright
Copyright © 2019 HongXia Shen. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.