1. Introduction
Concrete-filled steel tubes (CFT) exhibit a satisfactory load capacity, sufficient ductility, and adequate energy absorption capability. Additionally, the steel tube functions as a mold for pouring the concrete, leading to cost savings in construction. With the tube serving as both longitudinal and lateral reinforcement for the concrete core, no additional reinforcement is required. The enhanced structural characteristics of CFT members are primarily attributed to the combined effect of the steel hollow section and the core concrete. This composite action results in the core concrete being subjected to triaxial stress, while also preventing the steel hollow section wall from buckling inward [
1]. Hatzigeorgiou [
2] proposed analytical models for predicting the strength of circular CFT beam–columns, assuming that confinement enhances both the strength of the concrete core and its ductility. Furthermore, Zhang and Shahrooz [
3], and Hajjar and Gourley [
4] have proposed analytical models for predicting the strength of rectangular CFT beam–columns where these models assume that confinement only enhances the ductility of the concrete core, not the strength. Sakino et al. [
5] proposed design formulas to predict the axial compressive load capacity and ultimate moment of CFT columns with circular or square cross-sections. Mukai and Nishiyama [
6] suggested that the hoop stress on steel tubes for CFT columns is about 20% of the tensile yield stress of steel. Tang et al. [
7] developed a model for circular tubes that considers the effect of geometry and material properties on strength enhancement and post-peak behavior. Susantha et al. [
8] proposed appropriate stress–strain relationships for concrete confined by steel tubes of various shapes, adopting Tang et al.’s [
7] model for circular columns. Watanabe et al. [
9] conducted model tests to determine a stress–strain relationship for confined concrete and proposed a method to analyze the ultimate behavior of concrete-filled box columns, considering the local buckling of component plates and initial imperfections. However, there is limited research on the pure flexural behavior of concrete-filled hollow structural section beams. Furlong [
10] proposed that the flexural strength derived solely from the capacity of steel tubes was approximately 50% lower than that of the tested specimens. Subsequently, Elchalakani et al. [
11], Lu and Kennedy [
12], and Prion and Boehme [
13] conducted experimental investigations on numerous tests, revealing that the beam specimens exhibited adequate behavior and failed in a highly ductile manner.
This research offers a thorough and detailed examination aimed at forecasting the complete behavioral range of both rectangular and circular concrete-filled tube (CFT) members when subjected to a combination of axial forces and bending moments. The study is systematically organized into two distinct yet complementary sections, each significantly enhancing the overall comprehension and practical application of CFT structural components. The initial section focuses on the creation and application of various advanced finite element models. These models are intended to replicate a diverse array of loading conditions, including both concentric and eccentric axial loads on columns, as well as beams experiencing pure flexural loading. The main goal of this section is to establish a comprehensive and reliable response database that encapsulates the intricate behavior of CFT members under different loading scenarios. Several critical factors are meticulously considered and integrated into these advanced numerical models: the confinement effect, which is crucial for CFT behavior, is addressed by accounting for how it enhances the strength and ductility of the concrete core, while also considering the reduction in yield stress of the steel tube due to the biaxial stress state induced by confinement; the potential for slippage between the steel tube and the concrete core is also taken into account, which is essential for accurately predicting the composite action and load transfer mechanisms within the CFT member; finally, the models incorporate the possibility of local or global buckling of the steel tube, which is particularly significant for understanding the behavior of CFT members under high axial loads or slender configurations. By addressing these intricate interactions and behaviors, the finite element models deliver a comprehensive representation of CFT performance across a wide range of conditions. Therefore, this study advances the understanding of composite members beyond ultimate carrying capacity, offering a comprehensive behavioral prediction model for rectangular and circular CFTs under combined loads. This innovative approach integrates advanced finite element modeling with often-overlooked factors, such as confinement effects, material slippage, and steel tube buckling. By simulating diverse loading conditions, a comprehensive response database has been created that captures intricate CFT behavior. This holistic methodology, combining multiple critical factors in a single model, enhances the understanding of CFT behavior and provides valuable insights for practical design and analysis applications.
2. Materials and Methods
In this Section, the finite element analysis of CFTs is examined.
Figure 1 depicts the composite members sections examined here, which are the most applied in engineering practice, i.e., rectangular and circular sections.
For the rectangular section, B, H, t
f, and t
w are the section width, height, flange thickness, and web thickness, respectively, as shown in
Figure 1. Furthermore, for the circular section, D and t denote the diameter and thickness of the steel tube, respectively.
To develop an extended and reliable database for the behavior of CFT members with rectangular and circular sections, a 3D nonlinear finite element model was created, utilizing the ATENA [
14] finite element analysis. To ensure precision in the simulation, factors such as the nonlinear behavior of confined concrete, the local buckling of steel tubes, and the interface action between concrete and steel tubes were considered. Furthermore, selecting the appropriate element type and mesh size that offer precise results within a reasonable computational time frame is crucial in the simulation procedure. This model takes into consideration both geometrical and material nonlinearities. Various element types were tested to determine the most suitable ones for simulating the behavior of CFT rectangular and circular columns. It was found that shell elements are ideal for modeling the steel tube, while solid elements are best suited for representing the concrete core [
15,
16]. These element types were subsequently utilized in the finite element analyses. The steel tube is represented using 8-node shell elements, while the concrete core is modeled using 20-node solid elements. To account for symmetry, only half of the column is modeled, as depicted in
Figure 2. In all cases, this Figure is showing only half of each section, cut along the symmetry plane mentioned in the text. The meshes demonstrate the “refined discretization” described, especially at the base region. A refined discretization is applied to the base region of the column, with the height equal to the width of the cross-section for rectangular sections or to the diameter for circular sections. The nodes on the symmetry plane are appropriately constrained, while the remaining nodes are free to displace in any direction. To optimize computational efficiency, only one half of the column is modeled, resulting in a vertical plane of symmetry that bisects the column. The nodes situated along this symmetry plane are assigned specific movement limitations to accurately simulate the behavior of the entire column. These constraints generally encompass a) the prohibition of movement perpendicular to the symmetry plane, which ensures that the model functions as if the other half of the column were included; b) allowance for free movement parallel to the symmetry plane, which permits the model to deform naturally along the symmetry plane; and c) restriction on rotation around axes parallel to the symmetry plane, which prevents unrealistic twisting that would not occur in a complete column. To solve the nonlinear equations of motion, the Newmark step-by-step time integration method in conjunction with Newton–Raphson iterations is employed in ATENA [
14]. A suitable constraint model has been used for the interface where the relative movement between the concrete core and the steel jacket is minimal. Similar discretization and modeling of CFTs with rectangular and circular sections have been applied by Skalomenos et al. [
15] and Serras et al. [
16].
In this study, the following parameters have been investigated:
Concrete grades: C20, C30, C40, and C50, with concrete compressive strength, fc, equal to 20 MPa, 30 MPa, 40 MPa, and 50 MPa, respectively;
Steel grades: S200, S300, S400, and S500, with steel yield stress, fy, equal to 200 MPa, 300 MPa, 400 MPa, and 500 MPa, respectively;
The diameter of the circular sections, D: 200 mm, 300 mm, 400 mm, 500 mm and 600 mm;
The width of rectangular sections, B: 200 mm, 330 mm, 470 mm, and 600 mm;
The height-to-width ratio of rectangular sections, H/B: 1.00 and 1.20;
The diameter to thickness ratio of circular sections, D/t: 40, 55, 70, 85, and 100;
The width to thickness ratio of rectangular sections, B/tf: 40, 60, 80, and 100;
The width thickness to height thickness ratio of rectangular sections, tf/tw: 1.00 and 1.20;
The member length to section diameter ratio, L/D: 3.0, 6.0, 9.0 and 12.0.
Taking into account the aforementioned combinations of parameters for composite members with circular sections, suitable ones have (4 concrete grades) × (4 steel grades) × (5 diameters) × (5 diameter-to-thickness ratios) × (4 length-to-diameter ratios) = 1600 different models. Additionally, considering the aforementioned combinations of parameters for composite members with rectangular sections, suitable ones have (4 concrete grades) × (4 steel grades) × (4 widths) × (2 height-to-width ratios) × (4 width-to-thickness ratios) × (2 width thickness to height thickness ratios) × (4 length-to-diameter ratios) = 4096 different models. Therefore, the produced databank has 1600 + 4096 = 5696 finite element models, which have been examined under concentric or eccentric axial loads and pure flexural load.
3. Validation of Finite Element Results
In this Section, the validity of the aforementioned finite element analyses of circular and rectangular CFTs is examined. More specifically, some finite element models have been appropriately modified to reproduce the behavior of specific composite beams, columns, and beam–columns, taken from the pertinent literature.
Figure 3 depicts the response of concrete-filled steel tubes under pure bending employing experimental results from Ref. [
17] and finite element analysis.
Figure 4 shows the response of CFT beam–columns through experimental results from Ref. [
18] and finite element analysis.
4. Axial Load-Bearing Capacity of Concrete-Filled Steel Tubes
4.1. Rectangular Composite Sections
For the rectangular composite section, the area for the concrete core can be computed from
and the area for the steel tube is given by
The strict contribution of the concrete core to the maximum axial load can be found by
while the corresponding contribution of steel tube is given by
To quantify the composite action for the aforementioned contribution of the discrete materials, the following empirical expressions can be expressed for the concrete core
where
ac and
a1–
a3 are parameters that are presented in
Table 1. Similarly, the effective axial force for the steel tube is given by
where b
s and b
1–b
3 are parameters that are presented in
Table 1. It should be mentioned that B and H are the rectangular section’s width and height, respectively, and t
f and t
w are the steel tube section thicknesses for the flange and web, respectively. Thus, according to this study, the predicted strength for axial load results from
where parameter
c takes into account the effect of composite column length and can be
where
c1 and
c2 are the parameters that are presented in
Table 1. Therefore, from Equations (5)–(8), one has
Table 2 presents data from the experimental works of various researchers. Furthermore, the evaluation of maximum axial load using Equation (9) is also examined. It is found that the proposed empirical Equation (9) can predict the axial load capacity in a simple, yet effective and accurate manner. This is also obvious in
Figure 5, where the predicted axial forces are close to the experimental ones.
4.2. Circular Composite Sections
For the circular composite section, the area for the concrete core can be computed from
and the area for the steel tube is given by
The strict contribution of the concrete core to the maximum axial load can be found by
while the corresponding contribution of steel tube is given by
To quantify the composite action for the aforementioned contribution of the discrete materials, the following empirical expressions can be expressed for the concrete core
where
and
–
are parameters that are presented in
Table 3. Similarly, the effective axial force for the steel tube is given by
where b
s and b
1–b
3 are parameters that are presented in
Table 3. It should be mentioned that D is the circular section’s diameter and t is the steel tube section thickness. Thus, according to this study, the predicted strength for axial load results from
where parameter c takes into account the effect of composite column length and can be computed by
where
and
are parameters that are presented in
Table 3. Therefore, from Equations (14)–(17), one has
Table 4 presents data from the experimental works of various researchers. Furthermore, the evaluation of maximum axial load using Equation (18) is also examined. It is found that the proposed empirical Equation (18) can predict the axial load capacity in a simple, yet effective and accurate manner. This is also obvious in
Figure 6, where the predicted axial forces are close to the experimental ones.
5. Flexural Bearing Capacity of Concrete-Filled Steel Tubes
In this Section, the flexural capacity of CFTs is examined. More specifically, the bearing capacity of concrete-filled steel tubes with rectangular and circular sections is investigated under the action of pure bending.
5.1. Rectangular Composite Sections
For the rectangular composite section, the plastic section modulus for the concrete core results from
and the plastic section modulus for the steel tube is given by
The contribution of the concrete core to the maximum pure bending moment can be found by
while the corresponding contribution of steel tube is given by
To quantify the composite action for the aforementioned contribution of the discrete materials, the following empirical expressions can be expressed for the concrete core
where
ac and
a1–
a3 are parameters that are presented in
Table 5. Similarly, the effective axial force for the steel tube is given by
where b
s and b
1–b
3 are parameters that are presented in
Table 5. It should be mentioned that B and H are the rectangular section’s width and height, respectively, and t
f and t
w are the steel tube section thicknesses for the flange and web, respectively, as shown in
Figure 1. Thus, according to this study, the maximum pure bending can be estimated by
where parameter
c takes into account the effect of composite column length and can be computed by
where
c1 and
c2 are parameters that are presented in
Table 5. Therefore, from Equations (23)–(26), one has
Table 6 presents data from the experimental works of various researchers. Furthermore, the evaluation of the maximum bending moment using Equation (27) is also examined. It is found that the proposed empirical Equation (27) can predict the flexural capacity in a simple, yet effective and accurate manner. This is also obvious in
Figure 7, where the predicted values of the maximum bending moment are close to the experimental ones.
5.2. Circular Composite Sections
For the circular composite section, the plastic section modulus for the concrete core results from
and the plastic section modulus for the steel tube is given by
The contribution of the concrete core to the maximum pure bending moment can be found by
while the corresponding contribution of steel tube is given by
To quantify the composite action for the aforementioned contribution of the discrete materials, the following empirical expressions can be expressed for the concrete core
where a
c and a
1–a
3 are parameters that are presented in
Table 7. Similarly, the effective axial force for the steel tube is given by
where b
s and b
1–b
3 are parameters that are presented in
Table 7. It should be mentioned that D and t are the steel tube’s diameter and thickness, respectively, as shown in
Figure 1. Thus, according to this study, the maximum pure bending can be estimated by
where parameter c takes into account the effect of composite column length and can be computed by
where c
1 and c
2 are parameters that are presented in
Table 7. Therefore, from Equations (23)–(26), one has
Table 8 presents data from the experimental works from previous studies where the evaluation of the maximum bending moment using Equation (36) is also shown. It is found that the proposed empirical Equation (36) can effectively predict the flexural capacity. This is also obvious in
Figure 8, comparing Equation (36) to experimental results.
6. Concrete-Filled Steel Tubes under Combined Axial Force–Bending Moment
In everyday engineering practice, structural members are under the action of combined axial force–bending moment loading conditions. Thus, the reliable description of the members appears to be vital. This Section examines the behavior of concrete-filled steel tubes under the action of combined axial force–bending moment, i.e., the behavior of CFT beam–columns.
Lai and Varma [
41] and Lai et al. [
42], have elaborated on the development of slenderness limits that are utilized to classify concrete-filled tubular (CFT) columns into the categories of compact, noncompact, or slender, contingent upon the governing slenderness ratio (width-to-thickness b/t or D/t ratio, λ) of the steel tube in detail in their publications. The purpose of this does not elaborate on the use of these ratios. Their studies have shown that one of the most important parameters that affect the behavior of CFT beam–columns is the relative strength ratio, ξ, which can be defined as
In the following, this parameter is used to form the axial force–bending moment curve of CFT beam–columns, i.e., the P-M curve.
6.1. Rectangular Composite Sections
For the rectangular composite section, the relative strength ratio is given by (see Equations (1), (2) and (37))
For specific values of axial load, P, and relative strength ratio, ξ, the corresponding bending moment of the P-M curve can be defined by the following empirical relation
where m
1–m
7 are parameters that are presented in
Table 9.
The applicability and accuracy of the proposed empirical Equation (39) are shown in
Figure 9, comparing this relation to experimental results that have been taken from previous studies. The corresponding data for the studies examined here appear in
Table 10.
6.2. Circular Composite Sections
For the circular composite section, the relative strength ratio is given by (see Equations (1), (10) and (37))
For specific values of axial load, P, and relative strength ratio, ξ, the corresponding bending moment of the P-M curve can be also defined from empirical Equation (41)
where, in this case, the m
1—m
7 parameters are presented in
Table 11.
The applicability and accuracy of the proposed empirical Equation (39) for composite circular sections are shown in
Figure 10. More specifically, this figure depicts the proposed axial force–bending moment interaction curves for various circular CFTs beam–columns in comparison to the experimental results that have been taken from previous studies. The corresponding data for these studies examined here, appear in
Table 12.
7. Discussion
This study examined the behavior of concrete-filled steel tubes, with rectangular and circular sections, under pure axial/compressive loads, pure bending moments, and combined axial load-bending moment loads. From the statistical analysis of finite element results of
Section 2, the following empirical expressions have been proposed:
Equation (9): the maximum compressive bearing capacity for composite rectangular sections;
Equation (18): the maximum compressive bearing capacity for composite circular sections;
Equation (27): the maximum flexural bearing capacity for composite rectangular sections;
Equation (36): the maximum flexural bearing capacity for composite circular sections;
Equation (39): the normalized axial force–bending moment interaction curve for composite rectangular sections;
Equation (41): the normalized axial force–bending moment interaction curve for composite circular sections.
It is important to note that the previously mentioned equations were among the simplest equations that accurately represented the numerical data exhibiting concave curves in every direction. These equations were obtained through the utilization of Table Curve 3D software (Version 5.01, Systat, Richmond, CA, USA) after extensive testing of approximately 8000 mathematical equations. The selection of these equations was based on the criterion for achieving the minimum absolute residual error using the Pearson VII limit, specifically minimizing the sum of ln[√(1 − residual2)].
Equations (9), (18), (27) and (36) take into account the ratio of external section dimensions to the tube thickness, e.g., B/tf, H/tw (rectangular section) and D/t (circular section). The proposed empirical equations can be used for the following range: 30 ≤ (B/tf or H/tw or D/t) 120. Furthermore, these equations also consider the material strength ratio, fy/fc, which can range: 3.5 ≤ fy/fc ≤ 18.
Moreover, Equations (39) and (41) describe the axial force–bending moment curves for rectangular and circular sections, respectively, and they are expressed in normalized form, i.e., P/Pmax − M/Mmax. They take into account the relative strength ratio (see Equation (37)), which can range: 0.20 ≤ ξ ≤ 2.0.
The suggested equations eradicate the necessity of the existing section classification to estimate the compressive strength of the steel tube section caused by local buckling. It is worth noting that these empirical equations were derived from an advanced finite element model, capable of effectively simulating local buckling phenomena (refer to [
15,
16] for further details).
As detailed in
Section 1, the Finite Element Method (FEM) was utilized to simulate the system being studied. The findings from this analysis are depicted graphically throughout the manuscript, offering a visual summary of the principal results. In particular, the figures illustrate the outcomes of the FEM analysis. These visual representations provide an extensive overview of the system’s behavior under varying conditions, as anticipated by the FEM model. While the section on formula derivation emphasizes the theoretical framework, it is crucial to recognize that the FEM results serve to directly inform and validate these mathematical formulations. The interaction between the theoretical model and the FEM analysis represents a significant strength of this research.
8. Conclusions
A simple and effective model was developed to estimate the bearing capacity of the concrete-filled rectangular and circular tubes under axial load, pure bending, and combination axial load-bending moment load. Using appropriate constitutive models for steel and concrete, the finite element method is used to build an expanded database from which the ultimate behaviors of the CFT beams, columns, and beam–columns were determined by statistical analysis. This database took into account the local buckling of the steel tube, the confinement of the concrete core, and the likelihood of slippage between the steel tube and the concrete core. It is also shown that the suggested method’s ability to predict analytically the capacity of the rectangular and circular CFT beam, columns, and beam–columns is practical and reliable by comparing the proposed empirical expressions for the bearing capacity of composite beams, columns, and beam–columns with numerous experimental test results. Finally, the proposed empirical expressions eliminate the necessity of the section classification according to traditional structural codes by introducing a continuous function to estimate the compressive strength of the steel tube section caused by local buckling, taking into account that these equations were derived from an advanced finite element model, capable of effectively simulating local buckling phenomena.