Abstract

This study aims to investigate the hysteretic behavior of H-shaped honeycombed steel web composite columns with rectangular concrete-filled steel tube flanges (STHCCs). Taking the shear span ratio (λs), axial compression ratio (n), steel ratio of section (α), aspect ratio of section (D/B), yield strength of steel tube (fyfk), and compressive strength of concrete (fck) as the main parameters, we designed 22 full-scale STHCCs. By comparing the load-displacement curves between test and simulation, the rationality of finite element modeling method was verified. The quasi-static analysis of 22 specimens was carried out, and the influence regularity of different variables on the hysteretic behavior, skeleton curves, ductility, energy dissipation, resistance degradation, and stiffness degradation of STHCCs was obtained. The results show that the hysteresis curves of all the specimens show full shuttle shape and strong energy dissipation capacity. λs, α, and fyfk have great influence on the bearing capacity of skeleton curves. With the increase of α and fyfk, the initial stiffness of the specimens gradually increases. The stiffness degradation rate of the specimens gradually slows down, and the energy dissipation coefficient gradually decreases by increasing λs, α, and fyfk, but energy dissipation capacity is still at a high level. The resistance degradation of specimens increases gradually by increasing λs, α, fyfk, and D/B. The ductility of specimens gradually increases by increasing n, α, and fck. The maximum bulging deformation and maximum stress of specimens appear at the column foot. The trilinear skeleton curve model and restoring force model of STHCCs are established by statistical regression.

1. Introduction

STHCCs, short for H-shaped honeycombed steel web composite columns with rectangular concrete-filled steel tube flanges, are a kind of novel composite members formed by connecting two concrete-filled rectangular steel tube flanges with a honeycombed steel web. The core concrete can effectively prevent and delay the local buckling of the steel tube and the instability failure of the specimens; the steel tube also has certain restraint action on the concrete, which results in the core concrete under a state of triaxial compression; and the compressive strength and the ability to resist deformation of the core concrete can be significantly improved.

Since the twentieth century, a lot of research has been carried out regarding concrete-filled steel tube (CFST) columns and honeycombed columns at home and abroad. Sakino et al. [1] designed and fabricated 114 short CFST columns with different parameters of the shape of the steel tube, the tensile strength of the steel tube, the diameter-thickness ratio of the steel tube, and the strength of the concrete, and a formula for the ultimate axial bearing capacity was proposed. The mechanical behavior of concrete-filled double-layer circular steel tube columns was studied by Essopjee and Dundu [2]. The tests of 32 specimens subjected to axial compression were conducted, and the bearing capacity formula of concrete-filled double-layer circular steel tube columns was proposed. Anupriya et al. [3] carried out experimental research on honeycomb beams with stiffeners and without stiffeners, analyzed the failure mode of honeycomb beams under concentrated load, and also proposed a method for calculating the compressive strength of honeycomb beams webs after buckling. The nonlinear buckling analysis of 30 STHCCs with different parameters were carried out by Ji et al. [4]. According to the finite element (FE) simulation results, the formula for calculating the stability bearing capacity of STHCCs was derived. In order to investigate the axial compression behavior of short STHCCs with different parameters, an experimental study of 16 short STHCCs was designed and carried out by Ji et al. [5], and the calculation formula of the axial bearing capacity for this type of short column was established. In 2020, Ji et al. [6] investigated the eccentric compression performance of 17 STHCCs by ABAQUS software and derived the calculation formula of eccentric compression bearing capacity for STHCCs. Low-cycle repeated loading tests on CFST beam-column structures were conducted by Varma et al. [7], the FE model of the structure was established, and the reliability of the model was verified, which provided a basis for the calculation of the bearing capacity of CFST structures. Gajalakshmi and Helena [8] took the diameter-thickness ratio of steel tube as a variable to conduct low-cycle repeated loading tests on four circular CFST columns, analyzed the variations of hysteretic curves for CFST columns under different loading levels, and proposed a simplified damage accumulation equation of CFST columns. Low cyclic loading tests on recycled CFST columns and ordinary CFST columns were carried out by Chen et al. [9] and Tang et al. [10]. The influence of axial compression ratio, steel tube strength, steel tube thickness, and other parameters on the seismic behavior of CFST columns was discussed, and a fitting formula for skeleton curve based on Boltzmann mathematical model was proposed. The seismic capacity of 16 ultralarge high strength concrete-filled circular steel tube (HCFTST) columns under cyclic loading was analyzed by Wang et al. [11], and a strength model to predict the bending moment bearing capacity was proposed. Five specimens of multicavity CFST specially shaped columns were designed and fabricated by Yin et al. [12]. Through the combination of quasi-static and numerical simulation, the failure mode, hysteresis characteristic, stiffness degradation, ductility, and energy dissipation capacity of each specimen were obtained, and the effect of various parameters on the mechanism of force was explained. The calculation method of characteristic points’ values for skeleton curve for solid-web steel reinforced concrete T-shaped column was proposed by Liu et al. [13], and the restoring force model was established, which provided the basis for elastic-plastic seismic response analysis of this type of structure.

The research on STHCCs is mainly aimed at axial compression and stability, and few studies have been done on the hysteretic behavior of this kind of composite columns. Based on the research results of axial compression test of STHCCs, the quasi-static analysis of 22 STHCCs under low cyclic loading is further carried out by ABAQUS software. The influence of different parameters on the hysteretic behavior, skeleton curve, ductility, energy dissipation capacity, resistance attenuation, and stiffness degradation of STHCCs is investigated. The hysteretic behavior and failure form of STHCCs are obtained, and the restoring force model of such columns is established.

2. Specimen Overview

Taking the shear span ratio (λs), axial compression ratio (n), steel ratio of section (α), aspect ratio of section (D/B), yield strength of steel tube (fyfk), and compressive strength of concrete (fck) as the main parameters, we design 22 full-scale STHCCs. The schematic diagram and specific parameters of the specimens are shown in Figure 1 and Table 1. Considering the confinement effect of steel tube on concrete, pitch-height ratio, hole-height ratio, and slenderness ratio, the nominal axial compression ratio (n) [14] of the specimens is calculated as follows:where refers to the conversion slenderness ratio of the columns, fsc refers to the equivalent strength of CFST section, Asci denotes the section area of single-limb CFST, Nu denotes the design value of axial compression bearing capacity of the columns, φ represents the stability coefficient of CFST axial compression specimens, and s/ and d/ represent the pitch-height ratio and the hole-height ratio, respectively.

3. Finite Element Model

3.1. Constitutive Model for Materials
3.1.1. Constitutive Model for Steel

According to the metal plasticity theory, the stress-strain curves of steel are generally determined by the tensile test of standard specimens. Here, the ideal bilinear elastoplastic constitutive model (CM) is adopted as CM of steel.

3.1.2. Constitutive Model for Concrete

Mander et al. [15], Teng et al. [16], Han [17], and Pagoulatou et al. [18] have successively given the CM of confined concrete; the CM of unconfined concrete has been given in Code for Design of Concrete Structures (GB50010-2010) by China [19]; and the comparison between them is shown in Figure 2. Through comparative analysis, the CM of confined concrete proposed by Han et al. [17] is adopted as CM of concrete. It is known that concrete is anisotropic material with uneven internal structure and complex mechanical behavior, so the influence of plastic damage on concrete is considered during the FE modeling process.

3.2. Modeling Process and Boundary Conditions

The three-dimensional solid element with eight nodes (C3D8R) was adopted to simulate steel tubes, concrete, and steel webs, and the FE models of STHCCs were established by ABAQUS software [20], as shown in Figure 3. The contact interaction between steel tube and concrete was set as normal contact and bond slip contact in tangential direction, thus shear stress can be transmitted between steel tube and concrete [21]. The friction coefficient (μ) was set as 0.3 [22].

A reference point was set at the center of the upper surface of STHCCs as shown in Figure 3, and the reference point was coupled with the corresponding specimen surfaces, which could ensure that the stress on the specimen surface was uniform in the process of axial compression, and could avoid the bias of the members as well. The displacement in X, Y directions and the rotation in three directions of the reference point were restricted (Ux = Uy = URx = URy = URz = 0), while the bottom of STHCCs was completely fixed. In order to ensure the quality of meshing and the accuracy of calculation, 80 mm was taken as the mesh size of FEM, and the mesh shape was mainly hexahedron.

4. Experimental Verification of FE Models

In order to verify the rationality of the above-mentioned modeling method, numerical simulation analysis on 14 specimens selected from [5] was carried out, and the specific parameters of the 14 specimens are listed in Table 2. The load-displacement curves of specimens obtained by simulation analysis are shown in Figure 4. Compared with the existing test curves, the load-displacement curves obtained by numerical simulation were in good agreement with the test results. The simulation results of the axial compression bearing capacity () of the specimens were extracted correspondingly and compared with the existing test results () subsequently, and the maximum error was 5.58%, as shown in Figure 5. The error was within a reasonable range, which indicated that the modeling method was applicable to such composite columns.

5. Loading Procedure

In the process of quasi-static analysis, the axial load (N) is applied firstly to the top of the columns, then a force-displacement hybrid control program is used for lateral loading, and the loading procedure involved a force control step and a displacement control step as shown in Figure 6. In the force control stage, the increment of lateral load (P) was ±1/3 times the yielding load (Py). The loading procedure was changed into the displacement control stage, and the increment of displacement was 0.4 times the peak displacement (Δm). Each displacement load step was repeated twice. When the lateral load of the specimens dropped below 85% of the peak load, the specimens were regarded as failure [23].

6. Parameter Analysis

6.1. Hysteresis Curves

The hysteresis curves of 22 STHCCs under low cyclic loading can be obtained by FE analysis as shown in Figure 7. The elasticity, elastic-plastic, plastic, ultimate bearing capacity, energy dissipation capacity, and ductility of the specimens can be reflected from the hysteresis curves [24]. It is indicated in Figure 7 that the hysteretic curves of all the specimens are relatively full, and there is no obvious pinching phenomenon from the early stage to the late stage. α has a great influence on the hysteresis curves of STHCC specimens; with the increase of α, the hysteresis curve gradually becomes fuller, and the ultimate displacement of the specimens gradually increases. After reaching the peak load, the decrease of the bearing capacity of the specimens gradually slows down, and the ductility becomes better.

6.2. Skeleton Curves

The skeleton curves can intuitively reflect the force condition of STHCC specimens during the whole process of quasi-static simulation analysis, and the skeleton curves of 22 STHCCs are shown in Figure 8. It is noteworthy that the skeleton curves can be roughly divide into three stages: elastic deformation stage, elastic-plastic stage, and load decline stage, which are similar to the load-displacement curves of the specimens obtained by axial compression tests. The skeleton curves of the specimens show a linear trend at the initial loading stage; at the same time, the specimens show elastic deformation. With the increase of horizontal load, the skeleton curves change into nonlinear gradually, and elastic-plastic deformation appears for the steel tubes. The plastic damage of concrete continues to accumulate until the load reaches the peak load. After the peak load, the skeleton curves begin to decline. It is demonstrated that the specimens show excellent bearing capacity.

With the increase of λs, the initial stiffness of the specimens and the bearing capacity gradually decrease; it can be found that the magnitude of the decline of skeleton curves slows down gradually, as shown in Figure 8(a). Figures 8(b) and 8(c) show that the influence of n (n < 0.4) and fck (fck < 60 MPa) on the bearing capacity of the specimens is not obvious and that the slope of the curves is similar. Figures 8(d) and 8(e) show that with the increase of α and fyfk, the bearing capacity and deformation capacity increase obviously. Figures 8(f) and 8(g) show that with the increase of D/B, the slope of the descending stage of the skeleton curves decreases gradually.

6.3. Establishment of Skeleton Curve Model and Restoring Force Model
6.3.1. Simplified Skeleton Curve Model

In order to obtain the characteristics of STHCC skeleton curves, the skeleton curves should be normalized firstly [25], as shown in Figure 9.

It can be found from Figure 9 that the skeleton curves conform to the trilinear law [26, 27], which can be simplified as a calculation model, as shown in Figure 10, and the specific simplified rules are as follows:(1)A trilinear model consists of three stages: elastic stage, elastic-plastic stage, and decline stage, where the yield point, peak point, and limit point are considered.(2)The values of characteristic points for skeleton curve model are the average positive direction and negative direction values of characteristic points for each specimen.(3)It is considered that the specimens are at the elastic stage before yield point; hence, the elastic stiffness is selected as the stiffness. When the load value raises between the yield point and the peak point, the specimen is considered to be at the elastic-plastic stage. After the peak point, there comes a declining stage, in which the stiffness of the specimens remains at the same value, and the equation of the skeleton curve is shown in Table 3.

6.3.2. Comparisons of Regression Values and Simulation Values of the Skeleton Curves

The comparisons of calculation results and simulation results of the skeleton curves for 22 specimens are illustrated in Figure 11. It can be found that the calculation results are in good agreement with the simulation results, which indicates that the skeleton curve model obtained by statistical regression can reflect the load-displacement relationships accurately.

6.4. Stiffness Degradation

The stiffness of the specimens is reflected by the secant stiffness (Ki) [10, 28], which could be expressed bywhere Pi denotes peak horizontal load of the first cycle under the i-th control displacement and Δi refers to the corresponding displacement.

The secant stiffness (Ki) versus the horizontal displacement (Δ) relationships of 22 STHCCs are shown in Figure 12. It can be seen from Figure 12 that the stiffness degradation trends of STHCCs are similar. Ki shows a decreasing trend with the increase of Δ, and the curves gradually slow down with the increase of Δ. λs, α, D/B, and fyfk have a significant impact on the secant stiffness, and with the increase of λs, α, and fyfk, the stiffness degradation rate of the specimen gradually slows down. With the increase of D/B, the stiffness degradation rate of the specimen improves, and the stiffness degradation phenomenon becomes more significant.

6.4.1. Degradation Law of Stiffness

It is found that the secant stiffness of the specimens decreases constantly during the whole loading and unloading process. In Figure 13, point 1 represents the i-th horizontal load peak point under forward loading, and point 2 represents the i-th horizontal load peak point under reverse loading. The slope of the line from point 1 to point 2 expresses secant stiffness. The relationships between Ki/K0 and Δi/Δmax can be obtained by statistical regression based on the test results, as shown in Figure 14. The regression formula is expressed as follows:where Ki denotes the secant stiffness of the first stage horizontal load peak under the i-th control displacement and Δi refers to the i-th control displacement.

6.5. Ductility

Ductility is regarded as a main parameter to measure the deformation performance of the specimens [29], and it can be explained by the ductility coefficient (μ), which can be calculated by the following equation:where refer to positive limit displacement and yield displacement, respectively, and denote negative limit displacement and yield displacement, respectively.

6.5.1. Simulation Results

The ductility coefficients of 22 specimens determined by the above-mentioned calculation method are shown in Table 4. The existing research [25, 30, 31] shows that μ of steel reinforced concrete (SRC) columns ranges from 3.50 to 6.00 when λs value is greater than 2. Table 4 shows that the ductility of STHCCs is similar to that of SRC columns, which indicates that the specimens have good deformation behavior.

6.5.2. Shear Span Ratio

The relationships between Pu, μ, and λs of the STHCCs are shown in Figure 15. With the increase of λs, Pu and μ of the STHCCs decrease gradually. When λs changes from 2.08 to 2.45, 2.64, 3.33 in turn, Pu of STHCCs decreases by 17.32%, 5.97%, and 18.92%, respectively, and μ decreases by 9.21%, 1.01%, and 45.83%, respectively.

6.5.3. Axial Compression Ratio

The relationships between Pu, μ, and n of the STHCCs are shown in Figure 16. It can be seen from Figure 16 that when n changes from 0.1 to 0.3, Pu of the STHCCs decreases by 1.92% and 1.81%, respectively, and μ increases by 8.93% and 85.76%, respectively, in which n has a significant effect on the ductility of the STHCCs. With the further increase of n (n > 0.7) [12], μ of the STHCCs decreases obviously, so the limit value of n should be reasonably considered, which ensures that the STHCCs have good ductility [25, 30, 31].

6.5.4. Steel Ratio of Section

The relationships between Pu, μ, and α of the STHCCs are illustrated in Figure 17. It can be seen from Figure 17 that when α increases from 4 to 6, 8, and 10, Pu of the STHCCs increases by 14.86%, 13.32%, and 13.52%, respectively, while μ increases by 2.04%, 37.5%, and 58.73%, respectively, which indicates that the influence of α on the ductility of the STHCCs is significant.

6.5.5. The Compressive Strength of Concrete

Figure 18 shows the relationships between Pu, μ, and fck of the STHCCs. When fck increases from 30 MPa to 40 MPa, 50 MPa, and 60 MPa, Pu of STHCCs increases by 6.31%, 2.34%, and 9.56%, respectively, while μ decreases by 36.54%, 1.99%, and 20.01%, respectively. With the increase of fck, Pu of the specimens increases but the ductility decreases, which indicates that the ductility of the STHCCs is significantly affected by fck.

6.5.6. The Yield Strength of Steel Tube

The relationships between Pu, μ, and fyfk of the STHCCs are shown in Figure 19. When fyfk increases from 235 MPa to 345 MPa, 390 MPa, and 420 MPa, Pu of the STHCCs increases by 27.76%, 22.88%, and 2.86%, respectively, and μ increases by 29.47%, 10.49%, and 6.19%, respectively, which indicates that Pu and deformation capacity of the STHCCs improve with the increase of fyfk.

6.5.7. Section Aspect Ratio

The relationships between Pu, μ, and D/B of the STHCCs are shown in Figure 20. It can be seen from Figure 20 that when D/B increases from 1.6 to 1.73, 2.08, and 2.72, μ of the STHCCs increases by 1.58%, 1.86%, and 34.54%, respectively.

6.6. Energy Dissipation Capacity

The energy dissipation capacity is applied to evaluate the ability to absorb the energy released in the earthquake of the structure [9]. The energy dissipation coefficient (E) is used to evaluate the energy dissipation capacity of STHCC specimens, and the energy dissipation coefficient refers to the ratio of the shadow area to the triangle DFO and the triangle BEO area, as shown in Figure 21. The energy dissipation coefficient can be calculated as follows:

Figure 22 shows the horizontal displacement (Δ) versus the energy dissipation coefficient (E) relationships of the 22 STHCCs under the first reciprocating load of each stage. It can be found from Figure 22 that the energy dissipation capacity values of the 22 STHCCs are substantially similar.

The energy dissipation capacity of STHCCs is shown in Figure 22. With the increase of λs, n, fyfk, and D/B, E and the energy dissipation capacity of the specimens decrease gradually. It can be seen from Figures 22(c) and 22(d) that fck and α have little effect on the energy dissipation capacity of STHCCs. The existing research [32] shows that E of RC columns is about 1.0 and E of SRC columns is about 2.4; E of specimens in this article can reach above 2.0, so it is indicated that the STHCCs have excellent energy dissipation capacity.

6.7. Resistance Degradation

The resistance degradation coefficient (ξ) is used to measure seismic action of STHCC specimens as shown in Figures 23(a), 23(b), 23(d), 23(e), and 23(g). With the increase of λs, n, α, fyfk, and D/B, ξ of the STHCCs increases and the resistance degradation rate gradually slows down. Figure 23(c) shows that with the increase of fck, the resistance degradation of the specimens gradually decreases and the resistance degradation rate accelerates. It can be concluded that within a certain range, the increase of λs, n, α, fyfk, and D/B is conducive to improving the seismic performance of the structure.

6.8. Restoring Force Model

The restoring force model of STHCCs is established by statistical regression, as shown in Figure 24, in which points D and C represent the yield points of positive direction and negative direction. Points E and B represent the peak points. AB and EF represent the descent segment. The horizontal auxiliary line (L1) represents the straight line (P = 0.45PE) [13], and represents the symmetric line of L1 about horizontal coordinate axis.

The specimens are approximately at the elastic stage before yielding, and the paths of loading and unloading process are along the elastic stage of the skeleton curves; more specifically, the OD section is loaded forward, and the OC section is loaded backward after specimens unloading to the O point. It can be found that the loading stiffness and unloading stiffness are both the elastic stiffness of the skeleton curves.

After the specimens yielded, the loading path develops along the skeleton curves. Before reaching the peak load, the unloading and reverse loading are carried out according to the path of point E to the negative yield point C and then develop along the path of point C to point B. The overall loading path is D-E-1-C-B-1′-D.

After the peak load, the loading path develops along the skeleton curve from point E to point 4 (point 5). The unloading path is from point 4 (point 5) to point 2 (point 3), the reverse loading path is from point 2 (point 3) to point 4′ (point 5′), the reverse unloading path is from point 4' (point 5′) to point 2′ (point 3′), and then the forward load is loaded to the path of the skeleton curve. The overall loading path is 4 (5)-2 (3)-6 (7)-4′ (5′)-2′ (3′)-6′ (7′)-4 (5).

6.9. Comparisons between Restoring Force Model and Simulation Results

The hysteretic curves of STHCCs calculated by the restoring force model are shown in Figure 25. By comparing them with the simulation results, we can see that both are in good agreement; therefore, the restoring force model is feasible.

6.10. Destruction Form

The failure modes of STHCCs subjected to low cyclic loading are substantially the same, as shown in Figure 26. The column foot of the double-limb steel tube flange shows obvious outward expansion deformation. Due to the effect of axial force and horizontal force, the honeycomb steel web produces vertical compression deformation and transverse in-plane buckling deformation. The core concrete of the two-limb flange generates vertical compression deformation and bulging deformation of the upper column foot. The maximum stress of the steel tube appears in the inner side of the column foot, and the maximum bulging deformation is concentrated in the upper part of the column foot. The maximum stress of concrete occurs at the column foot and the outer side of the concrete columns. Therefore, it can be concluded that the steel tube exhibits confinement effect on the core concrete, which can effectively improves the strength of concrete and inhibits the damage of concrete. The honeycomb steel web also has a certain support for the CFST inside the flange, which can enhance the constraint effect of the inner steel tube on the concrete, the strength of the concrete, and the bearing capacity of the CFST flange. The core concrete in the steel tube effectively inhibits and delays buckling failure of the steel tube and improves the deformation capacity and hysteretic behavior of the STHCCs.

7. Conclusions

Based on the simplified constitutive model of steel and the nonlinear constitutive model of concrete, 14 STHCCs are numerically simulated by ABAQUS software. By comparison, the maximum error between the simulation and test results is 5.58%, which can verify the rationality of the FE modeling method.

The hysteresis curves of all the specimens show full shuttle shape and strong energy dissipation capacity. α has a great influence on the hysteresis curves of STHCCs; with the increase of α, the ultimate displacement of the specimens gradually increases. After reaching the peak load, the decrease of the bearing capacity of the specimens gradually slows down, and the ductility becomes better. When n < 0.4 and fck < 60 MPa, the effects of n and fck on the bearing capacity of specimens are not remarkable. With the increase of α and fyfk, the bearing capacity and deformation capacity increase obviously. The stiffness degradation rate of specimens gradually slows down, and the E gradually decreases by increasing λs, α, and fyfk, but the energy dissipation capacity is still higher than that of ordinary concrete columns. With the increase of λs, α, fyfk, and D/B, the resistance degradation increases; the degradation rate decreases; and n, α, and fck have great influence on the ductility of specimens.

Based on the hysteresis curves obtained by the quasi-static analysis, the trilinear skeleton curve model and the restoring force model are established, and the hysteresis rules are proposed correspondingly.

By observing the failure mode of STHCCs, the steel tube has significant constraint on the core concrete and improves the compressive strength of concrete effectively. The honeycomb steel web provides effective support for the rectangular CFST flanges, and the core concrete delays the buckling failure of the steel tube validly, which can reflect the good overall hysteretic behavior of the STHCCs.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors are grateful for the financial support received from the Natural Science Foundation of Heilongjiang Province, grant no. LH2020E018; Scientific Research Fund of Institute of Engineering Mechanics, China Earthquake Administration, grant no. 2020D07; Opening Fund for Key Laboratory of the Ministry of Education for Structural Disaster and Control of Harbin Institute of Technology, grant no. HITCE201908; 2021 Social Science Development Research Project of Hebei Province, grant no. 20210301135; Humanities and Social Science Research Project of Higher Education Institutions of Hebei Province, grant no. SQ2021115; General Research Project of Heilongjiang Higher Education Teaching Reform, grant no. SJGY20190096; and Northeast Petroleum University Guided Innovation Fund, grant no. 2020YDL-02.