1. Introduction
In recent years, the utilisation of aluminium columns in both civil and aerospace engineering has experienced significant growth, propelled by technological innovations and a steadfast commitment to sustainability, efficiency, and performance. This surge in usage underscores the inherent versatility and advantages offered by aluminium alloys, which have become indispensable in a wide array of structural applications. Aluminium, one of the most abundant elements on the earth, renowned for its corrosion resistance, favourable strength-to-weight ratio and similarities in design with stainless steel, has emerged as a preferred material choice in various industries (Georgantzia et al. 2021) [
1]. Builders and developers are incorporating aluminium columns into green building projects to meet sustainability goals and reduce carbon footprints due to their light weight compared to steel (almost less than 2.5 times). Moreover, innovations in fabrication technologies, such as extrusion and 3D printing, have expanded the design possibilities for aluminium members in civil engineering applications compared to steel counterparts due to lower energy requirement for manufacturing. Complex geometries and custom profiles can be produced, allowing for more efficient structural designs and optimised material usage. Furthermore, in regions prone to seismic activity, aluminium columns are preferred for their ductility and ability to withstand dynamic loads. However, despite its numerous benefits, the inherent characteristics of aluminium alloys, such as lower stiffness compared to steel alloys, necessitate meticulous analysis and design considerations, particularly in structural applications where load-bearing capacity is paramount. One example would be the instability causing to lower the load-bearing capacity of aluminium columns under axial compressive loading.
Recently, numerous studies have explored compression members constructed from cold-formed steel (CFS), particularly those with perforations, generating a substantial amount of data to guide standards in the design of steel columns. Therefore, before delving into the available literature on aluminium columns, we present a brief overview of the literature concerning CFS columns. Shanmugam and Dhanalakshmi (2001) [
2] examined the influence of web plate slenderness ratio and opening area ratio on the ultimate compressive strengths of perforated CFS channel stub columns. Moen and Schafer (2007) [
3] found that the presence of slotted holes in CFS columns affected the post-peak response and ductility, influenced by the column's cross-section type and length. Yao and Rasmussen (2012) [
4] explored the effects of perforations on inelastic stress distributions, load transfers, and failure modes of perforated simply supported plates and C-section columns. They observed distinct failure modes and stress distribution changes due to perforations. Kulatunga and Macdonald (2013) [
5] investigated the influence of perforation positions, while Kulatunga et al. (2014) [
6] studied the effect of perforation shapes on the ultimate compressive strengths of CFS columns with lipped channel cross-sections. Singh et al. experimentally studied the axial compressive capacity of CFS square hollow section (SHS) and rectangular hollow section (RHS) columns containing circular perforations, noting conservative but generally scattered predictions by existing design equations. Additionally, various studies investigated the compression behaviour of columns under geometrical imperfections using design methodologies such as the direct strength method (DSM) [
8], continuous strength method (CSM) [
9], and design codes such as Eurocode 3 (EC3) and Eurocode 9 (EC9).
Limited research has been conducted on the compressive capacity of aluminium columns, with most existing studies concentrating on SHS, RHS, and circular hollow sections (CHS). Zhu and Young [
13,
14,
15] conducted finite element (FE) investigations on various aluminium alloy hollow sections (SHS, RHS, and CHS), both welded and non-welded, focusing on their axial compressive capacity. They developed tailored design equations for aluminium alloy tubes with transverse welds at column ends. Zhou and Young [
16] also examined the effect of circular holes on the web crippling strength of aluminium alloy SHSs through experimental and numerical analyses. Mohandas et al. [
17] explored the axial compressive capacity of SHS stub columns made from aluminium alloy, finding that those constructed from the 6061-T6 grade aluminium exhibit promising structural behaviour as a potential substitute for steel stub columns. Su et al. [
18] conducted an experimental study on aluminium alloy SHS and RHS columns, to examine their cross-section capacity and to explore the potential leverage of strain hardening in design, with and without internal cross stiffeners. Later, Su et al [
19] investigated the influence of strain hardening and moment redistribution on the compressive behaviour and design of aluminium alloy structures. Using the CSM, they analysed approximately 900 experimental and numerical results. Their findings indicate that CSM yields more accurate mean resistance predictions and reduces variability for both determinate and indeterminate aluminium alloy structures, in comparison to the Aluminum Design Manual [
20], the Australian/New Zealand Standard [
21], and Eurocode 9 [
22]. Additionally, Feng et al. [
23] and Feng and Liu [
24] studied the flexural buckling performance of 6061-T6 and 6063-T5 normal-strength SHS and RHS columns with circular openings, comparing experimental and numerical findings with current standards. Furthermore, various column shapes, including RHS [
25], SHS [
26], CHS [
27,
28], angle-sections [
29,
30], I-sections [31-33] and others [
34,
35,
36], were investigated. The analysis revealed that existing standards generally offered inaccurate predictions of resistance for columns failing due to flexural buckling, local buckling, torsional buckling, or combined buckling.
This study draws upon the pioneering work by Tsavdaridis et al. [
37,
38] that resulted a series of highly optimised cross-sections tailored for aluminium columns, with a prime objective of augmenting their performance. Given aluminium lower stiffness relative to steel, conventional structural elements may necessitate optimisation of their cross-sections to enhance structural efficacy. While certain aluminium alloys, like AL 6061-T6, boast yield and ultimate strength akin to or surpassing common structural steels, their inferior stiffness underscores the need for cross-sectional optimisation strategies.
Expanding upon the groundwork laid by Tsavdaridis et al. [
37], who introduced 16 innovative aluminium section beam and column profiles, this study advances the investigation. Marinopoulou et al. [
38] utilised finite element analysis (FEA) with ABAQUS software to determine the ultimate compressive resistance values for stub columns, aligning FEA outcomes with those derived from CSM, DSM, and EC9 methodologies. With a focus on pin-ended columns featuring topology-optimised cross sections, this research builds upon the modeling approach outlined in Georgantzia et al.'s [
39] study, which explored the compressive behaviour of C-section aluminium columns. This study extends its scope to encompass stub and pin-end columns, with the FEA model of pin-ended aluminium columns validated against Georgantzia et al.'s findings. Following experiment validation, a parametric inquiry delves into both standard and novel cross sections, assessing the efficacy of modern design techniques. This extensive investigation spans 16 diverse cross sections across two column lengths, subjecting each to varying global and local imperfections in FEM simulations, totaling 288 distinct FEA models to determine the maximum load capacity of these columns.
The structure of the paper is as follows:
Section 2 outlines the validation of the FEA model against the experiments conducted by [
39].
Section 3 details the FEA models and the subsequent parametric investigation. The findings and their discussion are presented in
Section 4, followed by the conclusions in
Section 5.
Figure 1.
Aluminium columns of different structures from Greece [
40] (a) and Colombia [
41] (b).
Figure 1.
Aluminium columns of different structures from Greece [
40] (a) and Colombia [
41] (b).
2. Validation of the FEA Model
Georgantzia et al. [
39] conducted both experimental and numerical investigations on C-section columns made from 6082-T6 heat-treated aluminium alloy. In this paper, the FE models were validated against the experimental results as shown in
Figure 2.
Columns inherently possess geometric imperfections due to manufacturing processes and other external factors. aluminium columns tend to be more imperfect than the steel columns due to its nature. Aluminium has smaller elastic modulus, and it deforms easily. It triggers higher imperfection values for aluminium columns after its production, especially for the thin sections. Moreover, aluminium alloys are more ductile than the steel and alloys, the plastic deformations occurred during the manufacturing along with the thin sections. To incorporate these imperfections into the FEA, a two-step approach is adopted, comprising eigenvalue buckling analysis in advance of the nonlinear static analysis. The eigenvalue buckling analysis is employed to determine the imperfect shapes of the columns, which are subsequently factored and integrated before utilised as the initial conditions for the nonlinear static analysis.
The FEA model incorporates two types of geometric imperfections: global and local imperfections. Initial conditions for the static analysis are determined using the first global and local eigen buckling modes, which are then adjusted using imperfection factors. These factors are dependent on column dimensions such as effective length, wall thickness, and flange thickness. Specifically, the global imperfection factor is set at
, where
represents the column length, while the local imperfection factor is defined as
, with
representing the flange thickness of the column cross-section. Similar parameter values as those used in [
39] were employed here as they encompass typical configurations of aluminium columns.
The aluminium C-section column is represented using shell elements, depicted in
Figure 3. Utilising shell elements is a widely adopted practice for thin-walled columns, as demonstrated in the literature [
8,
10,
12,
39]. S4R element from the ABAQUS element library was employed. To simulate the pin-ended column test accurately, one of the rotational degrees of freedom is left unconstrained. This is achieved using KINEMATIC COUPLING constraints via the RBE2 element, as illustrated in
Figure 4, since both ends of the column are pin-ended. Specifically, the UX and ROTY degrees of freedom (DOF) are set free based on the column's orientation during the column test.
Mesh convergence studies were conducted to optimise computational efficiency. A column with a square hollow section measuring 100 mm × 100 mm × 5 mm and extending 1 m in length was modeled using four different mesh sizes ranging from 2 mm to 5 mm. After comparing the results, it was found that the buckling load value varied within approximately ±1%. Consequently, a mesh size of 5 mm was selected for all simulations to improve computational efficiency. Considering the substantial number of simulation cases—totalling 288—that must be addressed, computational efficiency is a crucial aspect of this study.
In the validation studies, a multilinear hardening material model was adopted. However, for the subsequent parametric studies, a bilinear hardening model was employed due to the limited tangent modulus of AL6063-T6. The material properties required to establish the bilinear hardening model were sourced from EC9 [
22].
The comparison between the maximum load and mid-height deflection results of the validation study is depicted in
Figure 5, showing an exceptional alignment between the outcomes. Specifically,
Table 1 provides a detailed comparison of the ultimate load between the findings of [
39] and those derived from the present investigation.
4. Parametric Study
In the parametric investigation, the novel cross-sections examined by Marinopoulou et al. [
15] were utilised with slight adjustments. Square hollow sections (SHS) featuring both inner and outer fillets have been included, as they offer better manufacturability compared to profiles without radii or with radii only on the inner or outer edges. The incorporation of corner radii enhances the effectiveness of extrusion manufacturing processes. Imperfection parameters are sourced from Georgantzia et al. [
39], and the various types of cross-sections are illustrated in
Figure 6 and
Figure 7. Detailed cross-sectional dimensions are provided in
Table 3.
Table 4 presents the comprehensive conditions encompassing cross-section types, imperfection specifications, and column lengths considered in this study. Each cross-section variant is tested across both column lengths. Also, for each length, three of the global imperfections have been applied to columns. For each global imperfection value, three of the local imperfections have been combined for the study. Subsequently, all specified local imperfection values are individually applied to each setup. Global imperfection values are determined based on the effective length of the column, while local imperfection values are calculated using the flange thickness.
For considering imperfections, the same approach as in the validation study earlier, is followed. To ensure mesh convergence, column A1 is simulated with four different mesh sizes ranging from 2 mm to 5 mm. Upon comparison, the buckling load values exhibit a variation of approximately ±1%. Consequently, to optimise computational efficiency, a mesh size of 5 mm is selected for all simulations. Given the substantial number of simulation cases computational efficiency is a paramount concern in this study.
A bilinear plasticity material model was adopted, as it is well-suited for this analysis due to the material's minimal tangent modulus. The model employs a 5 mm free mesh utilising S4R elements in ABAQUS. Boundary conditions are applied from both ends of the column, with a center node coupled using KINEMATIC COUPLING.
5. Results and Discussion
A set of 288 finite element models were carried out to assess the impact of imperfection values. Across all cases, there exists an average deviation of 3% between the lowest and highest imperfection values, with a maximum deviation of 5.4%. Given that the first global mode has been applied as the global imperfection factor, resulting in columns exhibiting a half-sine wave shape upon buckling. Notably, the alteration of local factors yields a minimal impact on the results, with variations less than 1% observed across all models.
Figure 9 portrays the behaviour exhibited by the columns, demonstrating a global buckling pattern with observable deformations on the walls and flanges induced by local imperfections. This underscores the predominant influence of global imperfections on the failure behaviour, attributed to the manufacturing process. In instances where columns have shorter lengths, the behaviour aligns with the first mode of global imperfection, leading to failure primarily at the walls and/or flanges. Conversely, longer columns exhibit a half-sine deflection along their length, indicative of Euler failure mode dependence.
In the comparative analysis, the finite element results for each column are averaged across all imperfection scenarios.
Table 5 provides a comparison between the results obtained from FE modeling and calculations based on EC9. It is observed that the disparity between FEM and EC9 calculations increase with the slenderness value, with EC9 values consistently surpassing FEM results.
Table 6 presents the results of Direct Strength Method (DSM) alongside their comparison with FEM results. Across all cases, FEM results are higher than those obtained from DSM, albeit with less sensitivity to slenderness compared to EC9.
Table 7 offers a comparison between FEM results and those derived from Continuous Strength Method (CSM). Notably, this comparison yields lower error values, with FEM results generally exceeding those of CSM. However, this trend shifts within certain slenderness ratios (0.32 - 0.43).
Figure 11 and
Figure 12 present the comparison results based on slenderness values outlined in the preceding section. The DSM exhibits precision for slenderness values below 0.4, while EC9 yields satisfactory results for slenderness values exceeding 0.22. The CSM generally outperforms other methods, consistent with findings in the literature. Notably, the CSM results exhibit the lowest coefficient of variation (COV), indicating superior precision. By directly calculating cross-sectional strength with plasticity correction, CSM offers adequate results without requiring significant modifications.
Figure 12 provides the precision of methods with respect to slenderness values. As it is mentioned before, CSM methods works perfectly within all slenderness values. EC9 works relatively better between the slenderness values 0.4-0.43, however it may need some improvement especially the other slenderness values. On the other hand, DSM works relatively better between the slenderness values 0.2-0.215 but it also needs some improvement. At the other slenderness values, it has divergent characteristics.
The nonlinear relationship between maximum load and slenderness values primarily stems from variations in cross-sectional inertia across different slenderness values. Since slenderness values encompass three distinct geometries, each with varying inertia-to-cross-sectional area ratios, this variability contributes to the observed nonlinearity in results. Additionally, differences in local imperfection modes across various cross-sectional types, influence both maximum load and column failure mode.