Effect of Concrete Age and Creep on the Behavior of Concrete-Filled Steel Tube Columns

The influence of concrete age and creep on the ultimate axial loading capacity of concrete-filled steel tube (CFST) columns is experimentally and numerically investigated. After validation of numerical models, a parametric study is conducted and the results are used to formulate empirical formulas for predicting the ultimate axial load-bearing capacity of the columns. Formulas are also proposed for predicting both the composite creep and aging coefficients of the CFST columns, which consider the confinement action of steel tubes on concrete. Then, the proposed formulas are validated independently by comparing their predictions with existing test results performed by other researchers. The comparisons show that the empirical formulas have the potential to be used in the practical design of CFST columns.


Introduction
Concrete-filled steel tube (CFST) columns have been extensively used in civil engineering structures over the past few decades.Extensive experimental and theoretical investigations on CFST components have been conducted and are well documented in the literature.The experimental work includes tests reported by O'Shea and Bridge [1] and Giakoumelis and Lam [2], who conducted research on the bearing capacity of short CFST columns.Zhong [3] and Zha [4] also reported their works on short and long CFST columns.Han et al. [5] and Wang [6] also studied property of the CFST under long-time load, while these results cannot be implemented conveniently considering the new code (GB50936-2014).Although significant work has been conducted on the properties of matured CFST columns, less emphasis has been placed on the influence of concrete age.Tan and Qi [7] have worked systematically by testing the effect of creep on columns under axially applied eccentric compressive load.Terry et al. [8] and Luo et al. [9] tested the creep of CFST columns and plain concrete in order to investigate the influence of steel tube reinforcement.Ichinose et al. [10] also reported a series of similar tests to obtain creep coefficients.Naguib and Mirmiran [11] and Liu et al. [12] developed an algorithm for the time-dependent behavior of CFST using the rate of flow method and the double power law functions for the creep of concrete.Shrestha and Chen [13] worked on the aging and creep coefficients of confined concrete, but the effect of steel tube was excluded in their formulas.
This study investigates the effect of concrete age and creep by experimental test and FEM modeling.Then, empirical formulas are proposed by treating a CFST column as a structural component made of a single material through homogenization.

Simulation and Methods
After pouring concrete into the steel tube, the elastic modulus of concrete increases with age.At the early stage, the steel tube is subjected to higher values of load than that specified by the design code.Thus, the ultimate bearing capacity of the CFST column is lower than the designed value.In this section, laboratory tests and FEM model are used to obtain the time-dependent load-bearing capacity.An empirical formula is established to provide a simple form of calculation for estimating the age-related ultimate load-bearing capacity.

FEM Modeling of CFST Columns.
The numerical model is implemented in ABAQUS eight-node quadrilateral inplane continuum shell elements which are used for the steel tube whereas the concrete is modeled employing eight-node linear brick element.The bilinear material model is adopted for the steel, with tangent modulus being 5% of the elastic modulus of 2.065 MPa.The material model used for the concrete is the damaged plasticity model, and hardening of concrete at different ages follows the curves proposed by Yi et al. [14], as shown in Figure 1.Other mechanical and geometrical parameters of the column are given in Table 1.FE meshes of both concrete and steel are shown in Figure 2. The contact between the steel tube and concrete is assumed to be perfect; here the nodes on the interface between solid and shell elements are tied.
Different meshing schemes are tested in order to assess the convergence of the FEM model, and the details are listed in Table 2, where  Ring,  Radius, and  Axial are, respectively, the numbers of elements in the circumferential, radial, and axial directions. Section and  Total are the total number of elements on a cross section and of the whole column, respectively. 3 d and  28 d are the ultimate axial loads of the column at ages of 3 days and 28 days.
Figure 3 shows terms of load against relative displacement of the different meshing FE model for concrete at the age 3 of and 28 days.With the increase in the number of elements, from M1 to M6, Table 2 shows the convergence of the ultimate load.When the elements number is more than that of M4, the difference between the different meshing schemes is nonsignificant.Therefore, meshing scheme M4 will be used in the following numerical simulations.
Laboratory tests have been conducted (see Figure 4) on the column specified in Table 1 in order to verify the FEM model.Load-bearing capacity tests are conducted at ages equal to 3, 7, 14, 21, 28, and 49 days.
The experimental test results of the time-dependent elastic modulus at different ages are shown in Table 3, and the load-bearing capacity of the CFST column is shown in Figure 5 as a function of the considered ages.Comparison between experimental and numerical results shows that the FE model can clearly predict the age-dependent load-bearing capacity.

Empirical
Formula of the Age-Related Load-Bearing Capacity.Parametric studies on a variety of CFST columns presented in Table 4 are conducted to derive an empirical formula for the load-bearing capacity of the CFST columns, whose numerical results are presented in Table 5.
Supposing that the age-strength relationship of concrete in the steel tube has the same form as when it is unconfined, according to ACI-209 [15], the following equation can be established: where the two unknown constants in the equation can be estimated by fitting FE results presented in Table 5.Then,  = 1.243 and  = 0.977 are obtained with the maximum error of 3.4%.Therefore, the ultimate load-bearing capacity of the CFST column is obtained as follows: The previously presented equation is validated through its application on the specimen defined in Table 1, and a comparison between experimental and numerical results is shown in Figure 6.(exp) and (FEM) are obtained applying (2) to the experimental and numerical results, respectively, shown in Figure 5.The predicted ultimate loading capacities obtained from the formula are consistent with those obtained from the experimental tests and FE modeling.

Influence of Creep on CFST Columns
In this section, a composite creep coefficient of CFST columns is proposed, in order to consider steel confinement, and the factor of the effective aging coefficient is obtained by regression on the FEM results.

Creep Coefficient of the CFST Column.
Yi et al. [14] proposed the age-adjusted effective modulus method, also known as BT theory, which represents an effective method to analyze the long-term effect of concrete creep.An aging coefficient was introduced in order to consider the effect of  concrete aging on the ultimate creep, and the integral equation for creep was transformed into an algebraic equation by employing the integral mean value theorem.When stress varies with time, the corresponding total strain at time  can be expressed as follows: where    () is the strain of concrete when creep is with no steel confinement at age  days,  0 1 ( 0 ) is the stress of the concrete at  0 days,   ( 0 ) is the elastic modulus of concrete at  0 days, and (,  0 ) is the creep coefficient at  when loaded at  0 days.
Applying the integral mean value theorem on formula (3) leads to Consider where  1 () is the stress of concrete at age of  days.  is the age-adjusted effective modulus. is the aging coefficient of concrete that is calculated as follows: The following calculation is based on the assumption of plane section and superposition principle of creep regardless of the toroidal shrinkage and strain.The axial strains of concrete and steel tube at age of  0 are denoted as  0 1 and  0 1 , respectively, whereas  0 1 and  0 1 represent the related stresses.The creep of concrete with no steel confinement is    =  0  (,  0 ) at the age of .As a result of the effect of the steel tube confinement, the concrete creep is reduced to   1 , as shown in Figure 7.
It is clear that the strain increment of the steel tube is the actual concrete creep.Hence, where   1 is the axial creep strain of the CFST column,   1 is the axial strain of steel tube due to concrete creep,    is the creep strain of the concrete with no steel confinement,   1 is the reduced strain of the concrete considering the steel tube confinement,  0 1 is the strain of the steel tube when loaded at  0 days, and  0 1 is the strain of the concrete when loaded at  0 days.
The total strain of CFST, or concrete in the steel tube, is derived as follows: From ( 8) and ( 9), it is obtained that and by considering ( 4) and ( 5) where   1 is the stress change of the concrete due to creep.In spite of extensive microcracks that will be induced in the core as reported by previous researchers and splitting cracks that can occur at 30% to 40% stress [16][17][18], the cracks could close if they are subject to an outer sustained compression  load at an early age.The confining stress would delay the formation of splitting cracks and restrict the widening of splitting cracks [19,20], and the confinement effect appears based on the splitting cracks formation [21].There will be no or little splitting cracks in the core concrete under the maximum sustained load after seven stages at an early age that is 40% bearing capacity of the specimens in the experiment and this had been contained naturally in the total creep.Therefore, the effect of splitting cracks on the creep was not considered separately here. is the area ratio of the steel with respect to that of concrete ( =   /  ), whereas  is the hoop effect coefficient ( =     /    ).
The equilibrium force condition in the axial direction satisfies the following equations: where   1 is the stress change of steel tube due to creep.Not considering the creep of the steel tube and taking steel as elastic modulus value a constant give and (11) becomes letting  =   /  ( 0 ) in ( 14), the change in stresses on concrete and steel tube because of the concrete creep is expressed as follows: Consider the following equations: Equations ( 16) are introduced into (15a) and (15b) to obtain the total creep strain   1 of the CFST column, as follows: Here the creep prediction function (,  0 ) employed the equation from ACI-209 [15].Then, the composite creep coefficient of the CFST,   (,  0 ) =   1 / 0  , is obtained as follows: (,  0 ) = ( −  0 )  where  cu is the ultimate creep of concrete;  la is the influence coefficient of loading age;  ℎ is influence coefficient of average thickness or the ratio of volume and area of the component;   is influence coefficient of the fine aggregate content;   is influence coefficient of the concrete slump;   is influence coefficient of the ambient humidity;   is influence coefficient of the air content.In this study,   (,  0 ) represents the aging coefficient that introduce the effect of not only concrete age but also the steel tube confinement when the CFST column is considered as a single homogeneous material, which is distinguished by (,  0 ), the aging coefficient of concrete.
Experimental test and FEM creep analysis of the CFST column are conducted in order to determine the aging coefficient,   (,  0 ).

Experimental Test and FEM Simulation of the Creep Model during Construction.
A special loading device was designed and manufactured in order to simulate the loading process during construction, as shown in Figure 8. Six columns, divided into two groups, were tested.The three columns in the first group are filled with concrete with recycled aggregates, whereas the other three columns in the second group are filled with concrete with conventional aggregates.The volume percentage of the expansive agent is 0%, 5%, and 10% for the three columns in each group, respectively.The design details of the columns are shown in Table 6.
The ultimate bearing capacity of the columns is designed as 1,100 kN, whose 40% represents the maximum dead load of 350 kN reached at the final stage of loading.The entire loading process is divided into seven stages and the load at each stage is recorded by the force sensor, as shown in Table 7 and Figure 9.An initial load of approximately 50 kN is applied to all columns of the two groups when concrete is three days  old.Dead load increases by approximately 50 kN at each stage with a three-day interval, until it reaches 350 kN.The creep deformation of the columns is measured by excluding the instantaneous displacements from the total displacement, as shown in Figure 10.Instantaneous displacements from the dead load are recorded immediately after such a load is applied, at each stage of the loading process.From Figure 10, the longitudinal deformation of the columns caused by the creep is notable and an increase in the expansion agent results in a larger creep deformation.
Column 4 is modeled by ABAQUS using a user-supplied subroutine of the creep effect.The creep model of ACI-209 [15], which considers the moist-cured condition, is adopted in the subroutine, where the progression method proposed by Gao et al. [22] is used to consider the change in concrete elastic modulus with age.
Creep deformation of column 4 obtained from the FE analysis is compared with the experimental results shown in Figure 11.In general, the FE solutions are consistent with the experimental test results.
After validating the FE analysis, the model is used to evaluate the aging coefficient   (,  0 ) of the columns specified in Table 8; here the composite creep coefficient   (,  0 ) is obtained by ratio between creep and instantaneous displacements from the numerical result, shown in Figure 12.    (,  0 ) can be back calculated employing (18a), (18b), and (18c) with known   (,  0 ) in Figure 12. Figure 13 shows the resulting aging coefficients, which are almost linear for the steel ratio  ranging from 0.05 to 0.20, representing the typical reinforcement ratio of practical designs.Thus, the aging coefficient within this range can be approximately computed as follows: Then, the composite creep coefficient   (,  0 ) of a CFST column can be calculated employing (18a), (18b), (18c), and (19).Predictions by (18a), (18b), and (18c) are compared with the test results from the studies of Terry et al. [8], in order to verify such equations, as shown in Figure 14, where Equ(Terry-1) is the equation result obtained by (18a), (18b), and (18c) related to the experimental test of Terry-1.A quite agreement can be seen, with some differences due to the different value adoptions of the ultimate creep of concrete and the loading duration.

Conclusion
In this study, the effect of concrete age and creep on the ultimate loading capacity of CFST columns had been experimentally and numerically investigated.The validated numerical models were employed in order to propose a set of empirical equations that can be used to predict the axial loading capacity of CFST columns with different concrete ages and composite creep coefficients.
According to the analysis of formula derivation and verification by experimental test and FE method, the ultimate bearing capacity considering the age of the CFST column can be calculated using (2).When treating the CFST column as a single material, both the unity creep and the aging coefficient, determined by employing (18a), (18b), (18c), and (19), respectively, can be used for long-term property analysis.
Equations in this study will be simplified to a certain degree on the design and time-dependent analysis of the CFST column.More studies considering the step upload will be conducted in the future.

Symbols
(): The elastic modulus of concrete at time   28  : The elastic modulus of concrete at 28th day   : The elastic modulus of the steel tube   : The elastic modulus of the CFST column   : The area of the cross section of the concrete   : The area of the cross section of the tube : Ratio between cross areas of steel and concrete : Th er a t i ob e t w e e nt h ee l a s t i cm o d u l u so f the steel and the concrete at  time  ut : Th eu l t i m a t eb e a r i n gc a p a c i t yo fC F S T column at  time (days)   : The bearing capacity of hollow steel tube   : The stability of the CFST column  28 : Th eu l t i m a t eb e a r i n gc a p a c i t yo fC F S T column at 28th day   ( 0 ): The elastic modulus of concrete at loading age  0 (,  0 ): The concrete creep coefficient  =   /  : The hoop effect coefficient : The outer diameter of the steel tube : The thickness of the steel tube : The length of the steel tube : The age of the concrete : The eccentricity of the load : The slenderness of the column    (): The total strain of concrete when creep is with no steel confinement at age  days  0 1 ( 0 ): The initial stress of the concrete when loaded at  0 days   ( 0 ): The elastic modulus of concrete at  0 days (,  0 ): The creep coefficient at  when loaded at  0  1 (): The stress of concrete at age  days : The aging coefficient of concrete   : The age-adjusted effective modulus   1 : The axial creep strain of the CFST column   1 : The axial strain of steel tube due to concrete creep    : The creep strain of the concrete with no steel confinement   1 : The reduced strain of the concrete due to the effect of the steel tube confinement  0 1 : The initial strain of CFST when loading  0 1 : The initial strain of the steel tube when loaded  0 1 : The initial strain of the concrete when loaded   1 : The stress change of concrete due to creep   1 : The stress change of steel tube due to creep   (,  0 ): TheunitycreepcoefficientofCFST   (,  0 ): The aging coefficient of CFST.

Figure 1 :
Figure 1: Stress strain relations of concrete at different ages.

Figure 2 :
Figure 2: Finite element mesh of (a) concrete and (b) steel tube.

Figure 3 :
Figure 3: The load curve for different meshes at the (a) age of 3 days and (b) age of 28 days.

Figure 4 :
Figure 4: Load-bearing capacity test of the CFST column.

Figure 5 :Figure 6 :
Figure 5: Load-bearing capacity of the CFST column at different ages.

Figure 10 :
Figure 10: Longitudinal creep deformation curve obtained by FE model.

2
Advances in Materials Science and Engineering

Table 2 :
Meshes parameter and results.

Table 3 :
Elastic modulus of the concrete at different age.

Table 6 :
Geometrical and mechanical properties of tests columns.

Table 7 :
Registered load at each stage (kN).

Table 8 :
[8]ameters of column needed in calculation.isthetime of the load duration;  cu is the ultimate creep coefficient of concrete according to Terry et al.[8]. TL