Abstract
Cross-section behaviour for cold-formed high strength steel (HSS) irregular hexagonal hollow section (IHexHS) stub columns under combined compression and bending is studied and presented in this paper. Finite element models were developed and validated using the existing experimental data collated from the previous research. Upon the validated finite element models, extensive parametric studies were subsequently carried out to generate more numerical data covering a wider range of cross-section dimensions, steel grades and load combinations from pure compression to pure bending. The obtained numerical results were utilised to assess the accuracy and the applicability of the current design codes, such as Eurocode of EN 1993-1-12 (EC3) and the North American code of ANSI/AISC 360-16 (AISC) for cold-formed HSS IHexHS stub columns under combined loading. It was demonstrated that the existing design codes can be safely applied and can be extended for cold-formed HSS IHexHS stub columns design under combined loading. In cross-sectional resistance predictions, conservative results were provided from the existing design codes. The over-predictions were primarily due to the neglect of the strain hardening and plate element interaction. The end points used in the interaction curves of EC3 and AISC adopt an idealised elastic-plastic material model to derive the corresponding resistance in cross-sectional level. The employment of Continuous Strength Method (CSM) leads to improved accuracy in cross-sectional resistance prediction with updated end points in the interaction curve. More consistent and reliable predictions were revealed by carrying out reliability analysis in accordance with EN 1990.
Keywords
Introduction
Steel tubular structural members are widely adopted in the construction industry including buildings, bridges and offshore structures because of the appealing appearance and prominent performance in resisting the buckling loads (Chan and Gardner, 2008; Liu et al., 2022a; Ma et al., 2015). The family of tubular members typically comprise conventional profiles of rectangular- (RHS), square- (SHS), circular- (CHS) and elliptical hollow sections (EHS). With the advancement in production technology and welding techniques, tubular sections in different cross-section configurations with combined structural efficiency and architectural aesthetics have attracted increasing interests from the researchers. Cold-formed steel oval hollow sections have been investigated by Zhu and Young (2011) and cold-formed semi-oval hollow sections have been studied by Chen and Young (2019, 2020). In recent years, increasing interests have also been made focusing on the applications of polygonal sections such as high strength steel (HSS) hexagonal-(HexHS), irregular hexagonal- (IHexHS), octagonal- (OctHS) and irregular octagonal hollow sections (IOctHS). The advantages of using HexHS/IHexHS and OctHS/IOctHS have been illustrated that the width of the flat portion is smaller than counterparts of square and rectangular hollow sections (SHS/RHS) if they are delivered with similar perimeters, therefore exhibiting stronger local buckling resistance than that of SHS and RHS (Liu et al., 2022g, 2022h). In addition, the flat surfaces allow for easier construction by providing an operating platform for welding or bolt-connection with end plate or gusset plate.
Experimental investigations on material properties and residual stresses for high strength steel (HSS) HexHS/IHexHS have been conducted by Liu et al. (2022a, 2022b) with sections made of Q690 steel grade (nominal yield strength fy = 690 MPa). Material properties and residual stresses investigations of OctHS/IOctHS have been experimentally studied by Zhu et al. (2019) with sections made of conventional mild steel with nominal yield strength of 355 MPa, and by Fang et al. (2018) and Liu et al. (2022c) with sections comprising Q690 steel plates with nominal yield strength of 690 MPa, and by Chen and Young (2020) with specimens made of Q460 steel plates with nominal yield strength of 460 MPa. In addition, local buckling behaviour of the HSS HexHS/IHexHS and OctHS/IOctHS stub columns under concentric compression have also been studied. The applicability of the existing structural steel design codes of EN 1993-1-1 (CEN, 2005), EN 1993-1-5 (CEN, 2006), EN 1993-1-12 (CEN, 2007), AISI/AISC 360-16 (AISC, 2016) and AS 4100 (AS, 2020) have been evaluated and compared. The previous studies demonstrated that the cross-section classification slenderness limit value cannot simply be extended to classify the cross-section of HexHS/IHexHS and OctHS/IOctHS under pure compression. The design equations in terms cross-sectional resistance under pure compression also exhibit either over-conservatism for stocky sections or unsafe predictions for sections subject to local buckling prior to the attainment of cross-section yielding. Despite the illustrative prominent advantages of HexHS/IHexHS and OctHS/IOctHS, limited studies have been carried out which further need to be systematically explored. Existing studies for cold-formed HSS IHexHS have been performed to study material properties, residual stresses and its structural behaviour at cross-sectional level under concentric compression. To better understand the structural behaviour and promote the application of HSS IHexHS, investigations of local buckling performance under combined compression and bending should be carried out.
Conventional tubular profiles of RHS/SHS and CHS with high strength steel materials are codified in current design codes of EN 1993-1-12, AISI/AISC 360-16. Note that provisions of high strength steel sections with nominal yield strength 690 MPa have been provided in AISI/AISC 360-16 (AISC, 2016), and EN 1993-1-12 (CEN, 2007) which extend the design framework for mild steels in EN 1993-1-1 (CEN, 2005) to high strength steel with strength grade up to steel grade of S700. Nevertheless, stringent requirements are proposed to tubular sections in AISI/AISC 360-16. AISI/AISC 360-16 specifies that hollow sections with nominal yield strength not greater than 485 MPa are only allowed to be designed. Whether the current structural steel design codes can be extended to cover the design of HSS IHexHS stub columns under combined compression and bending is unknown. The applicability of the existing design codes to the HSS HexHS stub columns needs to be assessed and evaluated.
Moreover, the design approach of Continuous Strength Method (CSM), which can rationally exploit the effects of strain hardening and redistribution of plasticity, was also evaluated and assessed for extending the application scope covering cold-formed HSS IHexHS stub columns under combined compression and bending. This paper therefore presents numerical investigation on cross-section behaviour of cold-formed HSS IHexHS stub columns under combined loading. The finite element models were firstly developed and validated against existing test results collated from cold-formed HSS SHS/RHS sections under combined compression and bending (Ma et al., 2019) and cold-formed HSS OctHS under combined loading (Fang et al., 2021). The validated finite element models were subsequently employed to conduct extensive parametric studies, covering a wider range of loading combinations and steel grades of Q460 and Q690. The obtained ultimate resistances from parametric studies were used to assess the applicability and accuracy of the current design codes and design approach of CSM.
Numerical investigation
General
A comprehensive numerical modelling programme was conducted including the development of the finite element models (FEM), validation and parametric studies. The FEM was developed using the software package of ABAQUS (ABAQUS, 2014). The adopted notations for the geometries of the cold-formed HSS IHexHS stub column are illustrated in Figure 1, where H is the height of the cross section, B is the outer width of the IHexHS, BL is the longer edge length of the IHexHS, bL is the clear width of the vertical flat portion with welding bead excluding the corner regions, Bs is the shorter edge length, bs is the clear width of an inclined short flat side excluding the corner portions, t is the nominal thickness of the cross section, ro and ri are the outer and inner corner radius of the cross section respectively, e0,y and e0,z are the initial eccentricity magnitudes around the major and minor axis respectively. The FEM was firstly validated against the existing test results of HSS cold-formed SHS and RHS in Ma et al. (2019) and Fang et al. (2021). Subsequently, extensive parametric studies were carried out to generate additional structural performance data on cold-formed HSS IHexHS stub columns under combined loadings. Definition of symbols for HSS irregular hexagonal hollow section.
Development of finite element model
The four-noded shell element with reduced integration (S4R) was employed in modelling cross-sections, which have successfully applied in modelling shell elements (Chan et al., 2015; Liu et al., 2022e, 2022f, 2023). Based on the preliminary mesh sensitivity studies, mesh size with a uniform value of bs/20 was adopted for the flat portions of the modelled specimens and a finer mesh of bs/30 was assigned to the corner portions to explicitly account for the curved geometry. To accurately represent the structural behaviours of the modelled specimens, measured geometries and experimentally tested engineering stress-strain relationships were used. The measured stress-strain curves were converted to the true stress-plastic strain curves in accordance with equations (1)–(2). Plastic material model with isotropic hardening was used for modelling steel materials.
To explicitly account for the effects of the cold-forming fabrication process, the enhanced material strengths were incorporated into the corner portions. Moreover, it was found that the enhancement of the strength was not limited to the corner portion only, observations have been made that strength increase can be extended to a certain distance away from the corner portion.
It was found that extension of the corner strength to adjacent flat portions with a magnitude of 2t achieve the best agreement between the numerical models and the experimental results for HSS cold-formed SHS/RHS in Ma et al. (2016). As for the cold-formed HSS IHexHS specimens examined in this study, no corner extension was needed basing on the preliminary sensitivity study of extension of cold-forming effect.
Linear elastic buckling analysis was performed to obtain the lowest elastic buckling eigen mode under the combined compression and bending, which was taken as the representative distributed profile of initial local geometric imperfection. The measured imperfection magnitudes were used and scaled to the corresponding profile from linear perturbation. For the modelled cold-formed HSS IHexHS stub columns, the imperfection measurement results were extensively analysed in Liu et al. (2022d) with predictive model of the initial local geometric imperfection (ωpred) proposed, as given in equation (3).
To simulate the pin-ended boundary conditions for eccentrically loaded stub column specimens, the displacement-controlled loads were applied to the reference points locating at the centroid of the modelled specimens, offsetting from the corresponding end sections at a distance equal to the sum of the thickness of the knife edge, end plate and wedge plate in longitudinal direction. The initial eccentricities were incorporated into the FE models by subsequently offsetting the reference points from the centroid of the end section to the measured values in the tests in the direction perpendicular to the bending axis. The reference points were constrained to the nodes at the edges of the end sections through kinematic coupling constraints. The reference points were restrained against all other degree of freedoms, except for the longitudinal displacements at the loading point and the rotation around the bending axis to simulate the pin-ended boundary conditions. The compressive load was applied by specifying the axial displacement of the reference point at one end with nonlinear geometric parameter activated to allow for large displacement analysis. Typical loading and boundary conditions of the cold-formed HSS IHexHS stub columns under combined compression and bending are shown in Figure 2. Typical loading and boundary conditions of the HSS IHexHS stub columns under combined compression and bending.
The fabrication processes of cold-forming and welding generate inter-locked existing residual stresses in an unloaded state. The residual stresses can be decomposed into membrane residual stresses and bending residual stresses. It is known that the bending residual stresses are considered to be inherently reintroduced to the stress-strain curves during tensile coupon tests in straightening the curved coupon specimen (Liu et al., 2022g). The effect of the residual stress on the structural performance were therefore considered by incorporating the membrane residual stresses with magnitudes determined from predictive model of residual stress distribution provided by Ma et al. (2015) for cold-formed HSS SHS/RHS and by Fang et al. (2018) for HSS cod-formed OctHS and by Liu et al. (2022b) for cold-formed HSS IHexHS.
Validation
Summary of the tensile coupon test results.

Comparisons between experimental and numerical load-end rotation curves for cold-formed SHS stub columns under combined loading (Ma et al., 2019). (a) H120×120×4-e5 (b) H120×120×4-e60.

Comparisons between experimental and numerical load-end rotation curves for cold-formed octagonal stub columns under combined loading (Fang et al., 2021). (a) OHS-70-6-e50 (b) OHS-70-6-e100.

Comparison of typical failure mode from experimental test and numerical simulation for cold-formed SHS stub column of H120 × 120 × 4-e30 (Ma et al., 2019).

Comparison of typical failure modes from experimental test and numerical simulation for cold-formed octagonal hollow section stub column of OHS-70 × 6-e75 (Fang et al., 2019).
Comparison of the experimental and numerical results.
Parametric study
Key material properties adopted in parametric study.

Stress-strain curves of cold-formed HSS IHexHS stub columns employed in parametric study (a) Stress-strain curves of flat portions (b) Stress-strain curves of corner portions.
The initial local geometric imperfection was incorporated in the FE model by scaling the predicted imperfection magnitude to the lowest elastic buckling mode. The membrane residual stresses were included in the FE models using the simplified predictive model of residual stress distribution for cold-formed HSS IHexHS provided in Liu et al. (2022b), as shown in Figure 8. Predictive model of membrane and bending residual stresses distribution and amplitudes in FE modelling (Liu et al., 2022b).
The length of all the specimens were designed as 2.5 times the average cross-section dimensions, which are short enough to avoid interaction effect from global buckling but sufficiently long to contain the representative initial local geometric imperfection and residual stresses. The thickness of the modelled IHexHS was kept constant as 6 mm, and the inner corner radius of the cold-bent corner of each specimen ri was taken equal to three times of its thickness. The effect of the cross-section slenderness was investigated by varying the lengths of the clear width of the vertical flat portion bL and the clear width of an inclined short flat side excluding the corner portions bs. The cold-formed HSS IHexHS is a novel cross-section which is not covered in existing design codes for traditional profile. Thus, cross-sectional slenderness parameter of λp = (fcr/fy)0.5 was used, where fcr is the elastic buckling stress of the cross-section under the applied loading conditions, which can be derived using finite strip software of CUFSM or ABAQUS, and fy is the measured yield strength from the flat portion. The cross-section slenderness of the numerical models for parameter studies were within the range of λp ≤ 0.64, based on the experimental investigation on cross-section classification of cold-formed HSS IHexHS in Liu et al. (2022d). Furthermore, the slenderness limit value of λp = 0.68 is specified as the boundary between slender and non-slender sections for CSM, the numerical results from parametric studies were thereof within the application range of CSM for non-slender sections. The load eccentricities were designed to yield a wider spectrum of compression and bending combinations. The ratios of the compression load to bending was represented using a dimensionless parameter, defined as radial angle θ, as illustrated in Figure 9. This dimensionless parameter can be derived as per equation (4), where PRd and MRd are the design resistances for stub column specimens under pure compression and pure bending respectively, and PEd and MEd are the design compression load and bending moment respectively. Based on the definition of the radial angle θ, the pure compression case can be expressed as θ = 0° whereas the pure bending scenario can be expressed as θ = 90°. Definition of on the load and moment interaction curve.

The obtained numerical results were used to evaluate and assess the applicability of the existing codifed provisions in Eurocode EC3, North American code of ANSI/AISC 360-16 as well as the newly developed design approach of CSM to cold-formed HSS IHexHS stub column specimens under combined compression and bending.
Assessment of design methods
General
Comparison of FE results with predicted design resistances for HSS IHex stub columns under combined compression and bending.
European code EN 1993-1-12 (EC3)
The Eurocode EN 1993-1-1 specifies design provisons for conventional mild steel with grades up to S460. EN 1993-1-12 extend the scope of applicability of the design rules for conventional mild steels to high strength steels with grade above S460 up to S690, though the design provisions in EN 1993-1-12 extensively refer to the design framework and provisions in EN 1993-1-1. Hence, for cross-sections made up of conventional strength steel and high strength steel, same design rules are provided for cross-sections under combined compression and bending. EN 1993-1-1 codified linear interaction relationship betwee the compressive load and the bending moment for cross-sectional design under combined loading, as given in equation (5).
In the linear interaction equation, PEd and MEd are the design compressive load and design bending moment respectively. PRd and MRd are the design resistance according to the provisions codified in design codes. As per EN 1993-1-1, the design resistance of cross-sections classified as Class 1–3 are able to attain the yield load Py, and the sections of Class 1–2 are able to achieve the plastic bending moment Mpl = fyWpl, whereas elastic bending moment Mel = fyWel is used for Class 3 cross-sections, where Wpl and Wel are the plastic and elastic section modulus respectively. As the investigated specimens were carefully designed to be short enough to avoid any effect from global buckling, the interaction factor due to the second order effect was approximately equal to the unity and the amplified second order moment is negligible. Thus, the design bending moment was taken as the initial eccentricity value multiplied by the design compressive load. At present, no specifications for determination of cross-section classification of HSS IHexHS under bending are developed, thus the Mpl was used for the HSS IHexHS with maximum bending moment larger than the corresponding Mpl when the e0/H and e0/B are eight and six respectively. It was found that the compression resistance of the investigated specimens under such eccentricities were far lower than 2% of its yield load, leading to negligible effect to the bending moment capacity. The predicted ultimate resistances were compared with the FE results generated from parametric studies. The compressive resistance was normalised to the predicted compressive resistance and is plotted against the normalised ratio between design bending moment capacity and the plastic bending moment, as shown in Figures 10 and 11. Figures 10 and 11 present the relationship between the normalised compression and normalised bending moment capacity about major axis and minor axis for cold-formed HSS IHexHS respectively with inclusion of EC3 interactive equation in Figures 10(b)−11(b). It is shown that EC3 design strengths are generally conservative, particularly for stocky sections, principally due to the conservative predictions of the end points of the design interaction curve. The end points of the EC3 interaction curve were derived without consideration of strain hardening and element interaction effect. In comparison with Q690 steel materials, specimens made from Q460 exhibit pronounced degree of conservatism for stockier sections as the increased level of strain hardening is associated with the stress-strain relationship of Q460 steel materials. As for Class 3 cross-sections, predictions from EC3 provide overly-conservative results indicating the current interaction design curve for IHexHS under combined compression and bending may be too optimistic. In particular, the increased level of conservatism is proportional to the strength grades of the steel materials. The cold-formed HSS IHexHS made of Q460 is overly-predicted by EC3 in comparison with the counterparts made of Q690 steel as shown in Figures 10 and 11 for combined compresion and both major axis and minor axis bending. Compared with the predictions for the specimens under combined compression and bending about major axis, less scattered predictions are obtained for sections under combined compression and bending about minor axis, though over-predicted resistances are still observed. The mean value of Pu/PEC3 are 1.21 and 1.28 for HSS IHexHS stub columns under combined compression and bending around major axis made of Q690 and Q460 materials with corresponding coefficients of variations (CoV) of 0.11 and 0.13 respectively. For HSS IHexHS stub columns under combined compression and bending around minor axis, relatively lower mean values of Pu/PEC3 were obtained which are 1.15 and 1.22 with the CoV of 0.03 and 0.07 for specimens comprising Q690 and Q460 materials. Comparison of the predicted cross-sectional resistance from EC3 for cold-formed HSS IHexHS under combined compression and major axis bending (a) Relationship between normalised ratio of ultimate load from FE to predictions from design code against radial angle (b) Interaction curve of normalised value of ultimate load to yield load against normalised bending moment. Comparison of the predicted cross-sectional resistance from EC3 for cold-formed HSS IHexHS under combined compression and minor axis bending (a) Relationship between normalised ratio of ultimate load from FE to predictions from design code against radial angle (b) Interaction curve of normalised value of ultimate load to yield load against normalised bending moment.

American specification ANSI/AISC 360-16
The North American design specification of ANSI/AISC 360-16 employ a pair of equations to design doubly symmetric cross-section structural components subject to combined compression and bending, as given in equations (6) and (7).
For compact sections (equivalent to Class 1–2 sections in EC3), plastic bending moment Mpl can be formed and elastic bending moment Mel are designed for non-compact sections in AISC design provisions. The design moment is determied based on the end moment of the cross-sections multiplied by amplificaiton factor to account for the second order effect, which can be dervied in accordance with the equation (8), where Pcr is the elastic buckling load of the column. In present study, this amplication factor is approximately euqal to the unity.
The numerical results were also used to assess and evaluate the applicability of the design code of ANSI/ISC 360-16 to HSS IHexHS. As the investigated HSS IHexHS are not covered by ANSI/AISC 360-16, the estimation of the design bending moment follow the same approach employed in EN 1993-1-12. Comparison of the FE results with the predictions from ANSI/ISC 360-16 (strength-to-predicted strength) is plotted in Figure 12 against the radial angle θ for HSS IHexHS under compression and major axis bending and the comparison of HSS IHexHS under compression and minor axis bending is shown in Figure 13. In addition, the compression resistance is normalised to the yield load and design moment is normalised to the plastic moment capacity with interaction design curve added in Figures 12(b)–13(b). As can be seen from the figures and comparison results, the ultimate resistances were significantly under-estimated by AISC design curves. The mean value of Pu/PAISC for Q690 and Q460 specimens are 1.16 and 1.22 with corresponding CoV of 0.10 and 0.13 respectively. In terms of the sections under compression and minor axis bending, the mean value of Pu/PAISC for Q690 and Q460 specimens are 1.11 and 1.17 with the CoV of 0.05 and 0.07 respectively. Compared with the predictions from EC3, the ANSI/AISC 360-16 provides more consistent and accurate predictions. Comparison of the predicted cross-sectional resistance from AISC for cold-formed HSS IHexHS under combined compression and major axis bending (a) Relationship between normalised ratio of ultimate load from FE to predictions from design code against radial angle (b) Interaction curve of normalised value of ultimate load to yield load against normalised bending moment. Comparison of the predicted cross-sectional resistance from AISC for cold-formed HSS IHexHS under combined compression and minor axis bending (a) Relationship between normalised ratio of ultimate load from FE to predictions from design code against radial angle (b) Interaction curve of normalised value of ultimate load to yield load against normalised bending moment.

Continuous Strength Method
The Continuous Strength Method (CSM) is a deformation-based design method that can account for strain hardening and element interaction (Afshan and Gardner, 2013; Gardner, 2008). The employment of the CSM, relating the cross-sectional resistance to the deformation capacity, replace the traditional methodology of cross-section classification. The CSM is a systematic design approach which originally grounded up focusing on stainless steel (Afshan and Gardner, 2013; Zhao et al., 2017) and is extended to the design of carbon steel sections for both hot-rolled sections and cold-formed sections (Yun et al., 2018; Yun and Gardner, 2018). The observations of disparity and over-conservatism in cross-sectional resistance predictions from design codes underpin the necessities to improve the accuracy of the predictions using design approach of CSM. The design approach of CSM comprise two key components, (i) Base curve, which build a continuous relationship between the cross-sectional slenderness λp and deformation capacity (maximum attainable strain εcsm) and (ii) CSM material model, which can explicitly capture the strain hardening of stress-strain curve. The design approach of CSM allows for the spread of plasticity and strain hardening, enabling the cross-section resistance particular for stocky sections to be considered explicitly. Hence, the CSM was used herein to generate more accurate end points anchored to the interaction curves, and the design equation of CSM for cold-formed HSS IHexHS under combined loading is proposed in a form analogous to EC3, given in equation (9).
The CSM have been developed and extended to be applicable for HSS tubular sections under pure compression (Lan et al., 2018) and pure bending (Lan et al., 2019). The base curve propsoed in Lan et al. (2018) for HSS tubular members was provided as follows in equation (10).
Two upper bounds of 15 and C1εu/εy are proposed on the normalised cross-section deformation capacity (εcsm/εy). The first limit of 15 is to avoid excessive plastic strain and subject to the material ductility requirement in EN 1993-1-1 (CEN, 2005) and the second limit of C1εu/εy prevents over-prediction of material strengths. The material model, applicable for design approach of CSM, comprises quad-linear and bi-linear material model. The key characteristics of the stress-strain curves of Q460 steel materials comprise yield plateau and moderate level of strain hardening, which conform to the typical mode of quad-linear material model demonstrated in Figure 14(a) and bi-linear model as shown in Figure 14(b) was used to represent the relatively rounded response of Q690 steel material with no sharply defined yield point and plateau. To determine the CSM design stress of fcsm, strain hardening slope Esh and the predicted strain εu corresponding to the ultimate stress are needed, as shown in equation (13). For quad-linear material model, the expressions given in equation (11) originally developed for hot-rolled steel by Yun and Gardner (2017) and bi-linear model employed in Yun and Gardner (2018) was used in this paper as shown in equation (12) for Q690 steel. Material model for design approach of CSM (a) Quad-linear material model for stress-strain curve with yield plateau and sharply defined yield point (b) Bi-linear material model for stress-strain with typically rounded response without yeild plateau. CSM material model for HSS materials.

The determination of the bending resistance of a cross-section from CSM largely depends on whether the CSM strain has entered the range of strain hardening (i.e., εsh < εcsm) or not. For steel materials with distinctive yield plateau, the CSM moment capacity can be calculated from equations (15) and (16).
The numerical results were used to assess the accuracy and the applicability of the proposed CSM method with updated end points anchored to the interaction curve. Comparisons of the FE results with the predictions from the proposed CSM (strength-to-predicted strength) are plotted in Figure 15(a) against the radial angle θ for HSS IHexHS under compression and major axis bending amd in Figure 16(a) for HSS IHexHS under compression and minor axis bending. The comparisons illustrate that the employment of the CSM leads to more consistent and less scatted resistance predictions compared with the design codes of EC3 and AISC. The mean Pu/PCSM for cold-formed HSS IHexHS under major axis of Q690 and Q460 are 1.06 and 1.10 with corresponding CoVs of 0.04 and 0.04 respectively. The mean values of Pu/PCSM for cold-formed HSS IHexHS under minor axis are also lower than the existing design codes with 1.08 and 1.12 for Q690 and Q460 steel materials with CoVs of 0.03 and 0.05 respectively. In addition, the interaction diagrams of normalised compression and normalised moment capacity are also plotted in Figures 15(b)–16(b) with notably improved accuracy in cross-section predictions. Nevertheless, Figures 15–16 generally reveal pronounced degree of conservatism as the applied loading varies from pure compression to pure bending for both design codes and proposed CSM interaction curve. Though more accurate and consistent predictions are obtained from CSM, predictions still exhibit conservatism for load combinations with lower axial ratios (larger radial angle). This can be principally attributed to two main issues, (i) the maximum deformation capacity was restricted to the strain of 15εy which impede the further exploitation of strain hardening, this restrictions can be relaxed based on certain case and (ii) the equations to derive the CSM bending moment was developed originally for SHS/RHS, the degree of plasticity formation and extent of strain hardening to which is different for IHexHS due to the effect of constituent plate element interaction from less onerous stress gradients. Comparison of the predicted cross-sectional resistance from design approach of CSM for cold-formed HSS IHexHS under combined compression and major axis bending (a) Relationship between normalised ratio of ultimate load from FE to predictions from design approach against radial angle (b) Interaction curve of normalised value of ultimate load to CSM predictions against normalised bending moment. Comparison of the predicted cross-sectional resistance from design approach of CSM for cold-formed HSS IHexHS under combined compression and minor axis bending (a) Relationship between normalised ratio of ultimate load from FE to predictions from design approach against radial angle (b) Interaction curve of normalised value of ultimate load to CSM predictions against normalised bending moment.

Reliability analysis
Reliability analysis results for HSS IHexHS stub columns under combined compression and major axis bending.
Reliability analysis results for HSS IHexHS stub columns under combined compression and minor axis bending.
Vδ is the CoV of the FE results relative to the resistance design model, which can be determined as follows, in equations (19)–(21).
The effect of the variability of the geometric and material in resistance function grt(X) are accounted for through their coefficient of variation parameter Vrt. The calculation of this parameter is shown in equation (22)
Vr is the combined CoV incorporating both model and basic variable uncertainties, and γM0 is the partial safety factor for cross section resistance. To determine the partial safety factor, design resistance value needs to be derived as follows,
Based on the statistical results reported in Tables 6 and 7, the design method of CSM provide partial factor γM0 of 1.15, 1.16 for cold-formed HSS IHexHS stub columns made of Q690 and Q460 under combined compression and major bending, which were lower than the partial factor needed for design codes of EC3 and AISC. Similarly, the resultant partial factors for cold-formed HSS IHexHS stub columns made of Q690 and Q460 under combined compression and minor bending are 1.09 and 1.11 respectively, all lower than the required partial factor derived for EC3 and AISC. The average values of the ratio between the numerical resistances to predictions from CSM are significantly lower than the design codes of EC3 and AISC, resulting in improved accuracy in CSM resistance predictions. Furthermore, lower values of the Vδ were determined for CSM leading to notable improvement in cross-section predictions and increased consistency.
Conclusions
A comprehensive numerical investigation into the cross-section behaviour of the cold-formed HSS IHexHS stub columns under combined compression and bending has been conducted and presented in this paper. The finite element models were developed firstly and validated against the experimental results collated from previous investigations for cold-formed HSS SHS/RHS as well as cold-formed HSS OctHS. The validated FE models were subsequently used to conduct parametric studies covering a wider range of load combinations with a total of 1800 numerical models. The obtained numerical results were then used to assess the accuracy and applicability of the current design codes of EN 1993-1-1, EN 1993-1-12, and ANSI/AISC 360-16. It was shown that the existing design codes yield consistent but unduly conservative predictions, primarily due to the conservative predictions of the end points of the pure compression and pure bending in isolated loading conditions. Compared with the predictions obtained from the design codes, the design method of CSM provides improved accuracy in cross-section predictions in terms of both compression and bending, and end points with increased accuracy were anchored in the interaction curves. Based on the reliability analysis, design method of CSM provides more accurate, consistent and reliable predictions than the existing design codes with lower partial factors and statistical parameters. The design codes and the proposed design approach of CSM can be extended to cover the design of cold-formed HSS IHexHS under combined loading, despite the unduly conservative predictions from design codes, particular for stocky sections experiencing strain hardening.
Footnotes
Acknowledgements
The research work presented in this paper was fully supported by a grant from the Research Grants Council of the Hong Kong Special Administrative Region, China (Project No. PolyU 15218918). The authors would also like to thank the support from Mr Shuai Li and Mr Jiachen Guo.
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the Research Grants Council of the Hong Kong Special Administrative Region, China (PolyU 15218918).
Data availability
Some or all data that support the findings of this study are available from the corresponding author upon reasonable request.
