Prediction of Bond-Slip Behavior of Circular/Squared Concrete-Filled Steel Tubes

: Numerous existing formulas predicted the ultimate interfacial bond strength in concrete-ﬁlled steel tubes (CFST) between steel tubes and concrete core without investigating the whole response under push-out load. In this research, four models are proposed to predict the interfacial behavior in CFST including the post-peak branch under the push-out loading test based on 157 circular specimens and 105 squared specimens from the literature. Two models (one for circular and one for squared CFST) are developed and calibrated using artiﬁcial neural network (ANN) and two models (one for circular and one for squared CFST) are developed based on multivariable regression analysis, analysis of variance (ANOVA). The shape of the specimen (circular or squared), diameter of the tube, thickness of the tube, concrete compressive strength, age at the time of testing, and length of the specimen are the main factors considered. These models are then compared to other existing formulas to verify their capability to better predict the ultimate interfacial bond strength. It is found that the ANN model gives better results for most of the considered data. It is also found that ANN models can predict the overall bond-slip response for the considered dataset. In order to simulate the response of any CFST column using ﬁnite element (FE) method, it is vital to have sufﬁcient input data on the overall bond-slip behavior between the interior face of the steel tube and the exterior surface of the concrete core including the post-peak branch. Accordingly, the suggested ANN model is used to generate the required input data related to the cohesive behavior and damage along the interface in ABAQUS model to simulate the response of two circular and two squared CFST columns under concentric compressive load. The results are in good agreement with experimental outcomes. The cohesive criterion and damage interface that are used based on ANN models in FE are found to be sufﬁcient and can be adopted to model CFST columns.


Introduction
Concrete-filled steel tube (CFST) columns are widely used in high-rise buildings, underground infrastructures, and bridges due to their ductility and lower weight compared to reinforced concrete. Their implementation in engineering structures dates back to the 1950s [1]. The use of this composite action of steel and concrete in CFSTs has several advantages over using pure steel or reinforced concrete members. The steel tube acts as formwork for the concrete and provides confinement to the core concrete which enhances the strength and ductility of concrete. On the other hand, the concrete core provides stiffness and compressive strength to the steel tube which reduces the chances of local buckling taking place in steel tubes.
The behavior of the concrete-steel bond in CFSTs has been extensively investigated and tested by many researchers. It has been found that the bond strength is affected by the

Experimental Database
The circular and squared specimens considered in this study were chosen from twelve previous studies [5,6,11,18,21,[24][25][26][27][28][29][30][31]. No treatment along the interfacial length between the concrete core and steel tube was provided such that the interaction between the two surfaces relied on the natural bond. Additionally, no shear connectors or studs were used in the considered specimens. Besides, the additives that replaced aggregates in the concrete mixture for some specimens were considered through the cubic compressive strength parameter, .
The considered variables of specimens were; diameter (D) or width (B) of the circular or squared specimen, respectively, the thickness of the steel tube (t), the interfacial length between the concrete and steel tube (Li), the cubic strength of the concrete mixture ( ), and finally the age of the concrete mixture.
Tables S1 and S2 are included as supplementary documents that summarize the details of circular and squared specimens from literature and used in this study.

Analysis of Variance (ANOVA)
A multivariate linear analysis of variance (ANOVA) was carried out to calculate the parameters that represent the interfacial behavior, namely; , , and α for both circular and square CFSTs based on the aforementioned test results. To achieve the best predictions by applying these equations, it is recommended to use specimens with (Li) less than 1200 mm and the concrete age (T) of not more than one year. The equations proposed for circular specimens are  The average bond stress was commonly used by most of the researchers to indicate the bond strength in the push-out test regardless of the non-uniform interfacial bond stress distribution [4] and was defined as where p f is the maximum load which results in breaking the bond, and A i is the interfacial area between the internal inner surface of the steel tube and the outer surface of the concrete core. However, different definitions of bond strength were used in the literature, such as the peak bond stress or the post-peak residual bond stress. While the slip value (S u ) which corresponds to the ultimate bond strength represents the relative slip between the concrete core and the steel tube. Nevertheless, the slip results include both the relative slip between the infill concrete and the steel tube as well as the compressive deformations of the concrete and steel tube. However, because these compressive deformations have a small effect on the bond mechanism and strength, and because it is difficult to distinguish them from interfacial slip results due to the experimental conditions, only relative slip between infill concrete and steel tube was considered.
After reaching the peak bond strength (τ u ), the push-out curves through the literature were noticed to have three different trends for the post-peak behavior depending on the macro-dimensional deviations [24]. This trend can represent hardening, softening, or constant plasticity as indicated in Figure 1a by trends a, c, and b, respectively. No predictions were found in the literature to capture both the ascending and post-peak branches which form the whole behavior of the push-out test.
In this research, the softening trend is adopted to represent the post-peak behavior of CFSTs, and, as a result, only specimens that followed this trend are chosen for analysis. This post-peak branch is expressed as an exponential function using β and α parameters as presented by Equation (2). This branch is normalized against the ultimate bond strength and shifted to the origin point as shown in Figure 1b, so the equation becomes

Experimental Database
The circular and squared specimens considered in this study were chosen from twelve previous studies [5,6,11,18,21,[24][25][26][27][28][29][30][31]. No treatment along the interfacial length between the concrete core and steel tube was provided such that the interaction between the two surfaces relied on the natural bond. Additionally, no shear connectors or studs were used in the considered specimens. Besides, the additives that replaced aggregates in the concrete mixture for some specimens were considered through the cubic compressive strength parameter, f cu .
The considered variables of specimens were; diameter (D) or width (B) of the circular or squared specimen, respectively, the thickness of the steel tube (t), the interfacial length between the concrete and steel tube (L i ), the cubic strength of the concrete mixture ( f cu ), and finally the age of the concrete mixture.
Tables S1 and S2 are included as supplementary documents that summarize the details of circular and squared specimens from literature and used in this study.

Analysis of Variance (ANOVA)
A multivariate linear analysis of variance (ANOVA) was carried out to calculate the parameters that represent the interfacial behavior, namely; τ u , s u , β and α for both circular and square CFSTs based on the aforementioned test results. To achieve the best predictions by applying these equations, it is recommended to use specimens with (L i ) less than 1200 mm and the concrete age (T) of not more than one year. The equations proposed for circular specimens are and for squared specimens are The τ u , S u , β and α results of ANOVA equations of circular and square CFSTs are listed in Tables S1 and S2 respectively. These results are also presented in Figure 2 for circular specimens and Figure 3 for squared specimens. By applying ANOVA analysis, the mean value (µ), standard deviation (σ), and R 2 of (τ u Exp. /τ u Pred. ) are presented in Table 1. It is found that the ANOVA results compare well with the experimental results and the mechanical-based discussion is presented in Section 4. ANOVA is an efficient parametric method for analyzing experiment data since it is practical and adaptable [32]. However, it is a considerably complex and subtle method to use as there are numerous ANOVA variations, each of which corresponds to a specific experimental situation. As such, it is possible to use the wrong type of ANOVA for a given experimental situation and draw incorrect conclusions from the data [33]. iate regression methods. Therefore, using ANOVA in this present study can be complicated because these parameters are correlated and affected by various factors.Accordingly, it is suggested to present another method of analysis for comparison. Machine Learning (ML) methods study links between inputs and outputs using statistical methods that enable computer systems to learn from a dataset without being explicitly programmed, which can be a good solution to this problem [34]. Several unique ML prediction methods have been used in recent decades. Artificial neural networks (ANNs) are one of them and are discussed further in the next section.     The predictions of this current study were based on specimens of various parameters that had been collected from various literature works then modeled using multivariate regression methods. Therefore, using ANOVA in this present study can be complicated because these parameters are correlated and affected by various factors.Accordingly, it is suggested to present another method of analysis for comparison. Machine Learning (ML) methods study links between inputs and outputs using statistical methods that enable computer systems to learn from a dataset without being explicitly programmed, which can be a good solution to this problem [34]. Several unique ML prediction methods have been used in recent decades. Artificial neural networks (ANNs) are one of them and are discussed further in the next section.

Artificial Neural Network (ANN)
It is vital to mention that several novel machine learning (ML) methods, i.e., artificial neural network (ANN), multi expression programming (MEP), and gene expression programming (GEP), have become very important forecasting and optimizing tools in many civil engineering fields. ML methods were introduced and utilized in predicting and modeling the non-linear behavior of different types of geoengineering materials like soil, concrete, asphalt, and clay which opened the door for revolutionary applications in the field of material science and optimization [35,36]. Additionally, ML methods can describe and recognize unnoticeable non-linear patterns among a high number of variables in extremely complex datasets due to the use of a non-linear activation function [37]. Moreover, ML has the advantage over regression in that the form of the model need not be pre-determined.ML was extensively utilized for predicting the mechanical behavior of concrete and concrete structural elements [38]. ML methods such as GEP were used to predict the ultimate load capacity of CFSTs [39]. On the other hand, ANN proved its ability to predict the mechanical properties of CFSTs such as the ultimate pure bending, the axial compression capacity, and fire performance [40][41][42].
In this study, the Feed-forward error Back-propagation (FB) method is utilized for calibrating the ANN models in this study. To model the behavior of CFST columns, two models are developed: ANN1 and ANN2 for circular and square columns, respectively. To develop the ANN models, 262 testing results datasets are utilized from the literature. These models are trained by utilizing the TRSEQ1 program, which was built on the FB algorithm [43]. The models are structured to have one hidden layer with a sigmoidal activation function.
ANN training can reflect efficient results if both inputs and outputs are normalized and scaled so the whole data vary from 0.1 to 0.9 before modeling. Once ANN models' inputs are scaled, the outputs are also be scaled, so the outputs are denormalized to give the results. Table 2 lists the parameters, maximum and minimum parameters values, and the ANN chosen maximum and minimum values for each parameter.   (12) and (13). However, ANN models' prediction is restricted to be within the actual maximum and minimum values. Equations (12) and (13) are utilized in normalizing and denormalizing the data respectively.
Normalized Parameter n = Parameter n − ANN n min ANN n max − ANN n min (12) Parameter n = (Normalized Value × (ANN n max − ANN n min )) + ANN n min (13) The datasets utilized in ANN modeling are divided into three groups: training datasets, testing datasets, and validation datasets. Training datasets is utilized to optimize the network's connection links weights for different networks with a varying number of hidden nodes from 1 to 10. For each trained network, the statistical accuracy measurements; Average Squared Errors (ASE), Mean Absolute Relative Error (MARE), and coefficient of determination (R 2 ) are calculated and listed. The testing datasets are used to test the accuracy of the trained networks. The trained network with the lowest ASE value of the testing datasets is defined as the network with the optimum number of hidden nodes. After the optimization, both training and testing datasets are combined into one train-all datasets group and utilizedtoretrain the optimum network to get the connection links weights of the train-all model. Combining these two datasetgroups allows the network to use more datasets in model training and capturing the relations between inputs and outputs which helps in increasing the model accuracy and reducing the error. The validation datasets (which were not used in training or testing the model) are utilized to check and validate the model by comparing them with the model's outputs. Figure 4 illustrates a 5-n-4 (5 inputs, n hidden nodes, and 4 outputs) ANN model schematic. However, 5-n-4 ANN model mathematical equations can be written as in Equations (14) and (15).
The developed ANN models are: 1. ANN1 (for circular column): 157 datasets of circular co tasets are divided into 117 training datasets, 30 testing datasets. ANN models' ASEtraining, ASEtesting, ASEtrain-all listed in Figure 5a. Based on the values shown in Figure  hidden nodes for ANN1 is found to be 9; because it ha Hence, ANN1 is denoted by its architecture as 5-9-4. 2. ANN2 (for square column): 105 datasets of circular co tasets are divided into 75 training datasets, 20 testing da tasets. As shown in Figure 5b, the optimum number o found to be 7. So, ANN2 is denoted by its architecture a The developed ANN models are:  Figure 5a. Based on the values shown in Figure 5a, the optimum number of hidden nodes for ANN1 is found to be 9; because it has the lowest ASE testing value. Hence, ANN1 is denoted by its architecture as 5-9-4. 2.
ANN2 (for square column): 105 datasets of circular columns are utilized. The datasets are divided into 75 training datasets, 20 testing datasets, and 10 validation datasets. As shown in Figure 5b, the optimum number of hidden nodes for ANN2 is found to be 7. So, ANN2 is denoted by its architecture as 5-7-4. Figure 5 illustrates how combining testing and training datasets together and including them in the models training helps in reducing the ASE values and increasing the models' accuracy.
The connection links weights for ANN1 and ANN2 to be used in Equations (12) and (13) are shown in Tables 3 and 4, respectively. However, ANN modeling results are listed in Tables S1 and S2. The connection links weights for ANN1 and ANN2 to be used in Equations (12 and (13) are shown in Table 3 and Table 4, respectively. However, ANN modeling re sults are listed in Tables S1 and S2.    For both models, Table 5 lists the training, testing, train-all, and validation statistical accuracy measurements. However, Figure 4; Figure 5 show the experimental datasets versus the ANN prediction results of the validation results for ANN1 and ANN2 models respectively. Both validation statistical accuracy measurements (listed in Table 5) and the validation predicted versus experimental data results (shown in Figures 6 and 7) illustrate that ANN1 and ANN2 predictions are significantly reasonable and acceptable. Excel sheets of the ANN models are provided as supplementary documents.

Comparison of Bond Strength ( ) Predicted by ANOVA and ANN with Existing Models
Results of based on proposed equations in this paper and equations from literature for circular CFSTs of Roeder et al. [7], Xue and Cai [20], Lyu and Han [21], and Chen et al. [19]and equations for squared CFSTs of Parsley et al. [18], Lyu and Han [21], and Xue and Cai [20]are compared. It can be concluded that results by ANOVA and ANN

Comparison of Bond Strength (τ u ) Predicted by ANOVA and ANN with Existing Models
Results of τ u based on proposed equations in this paper and equations from literature for circular CFSTs of Roeder et al. [7], Xue and Cai [20], Lyu and Han [21], and Chen et al. [19] and equations for squared CFSTs of Parsley et al. [18], Lyu and Han [21], and Xue and Cai [20] are compared. It can be concluded that results by ANOVA and ANN proposed models can better represent most of the specimens compared to available equations. Mean values (µ) and standard deviation (σ) results estimated by ANOVA and ANN formulas and by the existing models are presented in Table 6.

Comparisons of τ u , S u , β and α Predicted by ANOVA and ANN for Both Circular and Squared CFSTs
It is found that most predicted values of τ u , S u , β and α by ANN models are closer to the experimental than ANOVA results (Table 1). ANN outperforms the ANOVA and provides more accurate prediction results due to the use of the sigmoidal activation function which allows the ANN to capture the non-linear mechanical behavior of concrete. Figure 8 shows the total interfacial behavior obtained experimentally of one circular and one squared specimen versus the predicted behavior using ANN and ANOVA analysis. These two specimens were previously investigated by Abendeh et al. [26] The results also confirms that the ANN analysis can better represent the bond-slip behavior of CFSTs. Proposed excel sheets for predicting the entire bond-slip behavior of circular and squared CFSTs by using ANN are available as supplementary documents. Figures 9 and 10 present the whole experimental and predicted interfacial behavior (based on ANN models) of 6 circular specimens and 6 squared specimens, respectively. These specimens are selected from the different references adopted in this study to verify the capability of ANN models to represent the whole experimental bond-slip curves. For all predicted curves, the first point, after the ultimate bond strength is reached, on the post-peak curve is selected at S u plus standard deviation. tions, the curve showed a constant plasticity trend.
The difference between predicted and experimental results is considered reasonable since the test results of identical specimens encountered such a difference due to variation of properties of concrete, manufacturing, and previously mentioned irregularities in internal forces. For example, Figures 9 and 10 present the difference in τ u of the three identical t2.5-1.2%-28d specimens, in S of the three identical CTRL-SS1 specimens, in β and α of the two identical specimens of Parsley.    It should be stated that the proposed bond-slip behavior model ( Figure 1) follows what was observed inthe literature; the bond stress increase linearly, till the maximum bond stress value is reached. At this stage, no separation is seen in the interface between the concrete core and steel tubes. Hence, as both materials still work together, chemical adhesion dominates the bond resistance. Furthermore, the micro-interlocking forces increase with an increase in the push-out load. However, the macro-interlocking force does not contribute since no slip occurred between the various materials.
After reaching the peak, three kinds of curves are noted based on the macro-dimensional deviation [21,24]. The macro-interlocking force is produced from the confinement of steel tubes that prevented lateral expansion. It is also directly correlated to the size effect, which is produced by lateral expansion in concrete or the presence of irregularities in the internal surfaces of the steel tubes. Additionally, this force increases with an increase in the separation and relative slip. When the macro-dimensional deviation is large, the increase in friction is higher compared to the reduction in the initial friction. This leads to a hardening trend. On the other hand, if this friction value is lower compared to the reduction, the curve shows a softening trend. In intermediary situations, the curve showed a constant plasticity trend.
The difference between predicted and experimental results is considered reasonable since the test results of identical specimens encountered such a difference due to variation of properties of concrete, manufacturing, and previously mentioned irregularities in internal forces. For example, Figures 9 and 10 present the difference in τ u of the three identical t2.5-1.2%-28d specimens, in S u of the three identical CTRL-SS1 specimens, in β and α of the two identical specimens of Parsley.

Finite Element Modeling of CFST's Columns
In this paper, ABAQUS is used to model previously tested CFST columns under axial loads. The development of such models involves the discretization of components, the selection of the materials' constitutive models, and the approach to defining the interaction along the interfaces.
To develop calibrated finite element models of CFST members, it is vital to quantify and qualify the interaction characteristics along the concrete-steel interface. Abu Shamah and Allouzi [44] defined a surface to surface contact between the core-concrete and the steel tube using a shear stress limit according to Han et al. [45]. After the shear stress limit is reached, the softening behavior is not considered. The importance to identify the concrete-steel interface model was emphasized by Effendi [46]. Interface element and contact analysis were used along the interface and it was found that this method is better than rigid bar elements since rigid bar elements artificially add stiffness to the system. Nguyena et al. [47] investigated the use of the tie and slip contact properties between concrete and steel tubes. It was found that the tie model represented the test outcomes better than the slip model. Neither the tie model nor the slip model represents the actual interaction along the concrete-steel interface. Alatshan and Mashiri [48] defined the friction coefficient along the interface between the concrete and steel tube. The magnitude of this coefficient was reviewed from the literature.
The ANN models developed in this paper are used to define the interaction along the steel-concrete interface of simulated CFST columns.

Experimental Data for Calibration
Four CFST columns are developed using ABAQUS to verify if the proposed ANN models of the concrete-steel interaction can be used to simulate the interface characteristics along the concrete-steel interface of these columns. The details of these columns are presented in Table 7.

Finite Element Discretization
Concrete and steel are modeled as 8-noded 3D elements. The mesh size adopted to model specimens "SA" and "MA" is8 mm. A larger mesh size of 25 mm is used to model specimens "SU-40" and "SU-70" due to the larger sections of these specimens. As supports are modeled as a rigid body, the mesh size has no effect; accordingly, a mesh size of 20 mm is used.

Constitutive Models of Materials
The Concrete Damaged Plasticity model available in ABAQUS [23] is used which is a continuum model to represent the behavior of brittle materials such as concrete. This model involves tension and compression plasticity.
The input data of plasticity of concrete under compression is defined based on Tsai's equation [50] as The strain at the compressive strength f c (in MPa) is taken as [51]. n and r are parameters to control the shape of the stress-strain curve and are taken as E c ε c f c and f c 5.2 − 1.9, respectively [52]. The plasticity of concrete under tension is modeled based on Tsai's equation as . n t and r are parameters to control the shape of the curve. The strain at peak tensile strength ( f t ) is ε t0 .
For the concrete plasticity definition, the dilation angle of 22 • is adopted. The flow potential eccentricity is taken as 0.1 and the ratio of biaxial to uniaxial compressive strength is taken as 1.16. The invariant stress ratio is 0.667 and the viscosity is 0.001.
An elastic-plastic model is used to simulate the steel tube.

Tube-Core Interface
Cohesive behavior is implemented between steel tube and concrete core, and it is represented based on traction versus separation law. Damage is also enforced to simulate the post-peak behavior in cohesion over two steps; damage initiation and damage evolution.
Damage initiation is selected based on the quadratic criterion of contact stress ratios that is represented based on ABAQUS manual [23] as where t n ,t s ,t t are the contact stresses and t o n , t o s ,t o t are the peak values of the contact stresses. When the damage initiation criterion is met, the damage evolution is defined based on tabulated values of damage ratios versus slip displacement after the peak bond strength. The proposed ANN models presented in this paper define the bond-slip behavior of concrete-filled steel tube after bond strength that can be used to define damage evolution. The damage, d, is calculated as where τ is the bond stress at specific slip displacement after the peak bond strength.

Finite Element Results
The axial load versus axial displacement curves of the four axial specimens are estimated using ABAQUS [23]. The measured versus numerically calculated responses of these specimens are shown in Figure 11. As shown in this figure, the simulated models can represent the experimental results very well for columns "SA" and "MA". Note that the length of column "MA" exceeds the maximum length of the data of circular specimens used to develop the ANN models and this can be the reason for the difference between numerical and experimental results. For columns "SU-40" and "SU-70", the results can represent the measured experimental results well and it is expected to get even better results if the steel constitutive model is tri-linear but the details of steel material were not provided in the reference. This verifies that the used cohesive criterion and damage based on the ANN models are sufficient and can be adopted to model CFST's columns under axial load. Buildings 2022, 12, x FOR PEER REVIEW 18 of 22 Figure 11. Load-deflection curves of the simulated models.

Sensitivity Analysis
The sensitivity of the results of finite element modeling to the input parameters is essential to study since it shows the effectiveness of the model and accordingly the outcomes. The main input parameters that are investigated for sensitivity are; the mesh size of steel tube and concrete and the dilation angle of concrete.
The dilation angle of concrete is the angle of internal friction and it controls the amount of plastic volumetric strain. For the definition of the concrete constitutive model, the dilation angle usually is selected based on the value that can represent the experimental response of one specimen and verify it for the other specimens [53]. The loaddeflection results of specimens SA and SU-40 that are simulated with dilation angles of 22° and 30° are presented in Figure 12. It is found the circular column is more sensitive to the dilation angle value compared to the squared section.
A mesh convergence study is performed where each column is simulated with different mesh sizes and the results are compared. Sufficiently small elements are used. Three mesh sizes are adopted to model specimens "SA" and "MA", namely; 8 mm, 15 mm, and 20 mm. Three mesh sizes are used to model specimens "SU-40" and "SU-70", namely; 15 mm, 20 mm, and 25 mm. The results of load-deflection of all considered columns at these different mesh sizes are introduced in Figure 13. It is found that the peak load results are not affected by the different mesh sizes.

Sensitivity Analysis
The sensitivity of the results of finite element modeling to the input parameters is essential to study since it shows the effectiveness of the model and accordingly the outcomes. The main input parameters that are investigated for sensitivity are; the mesh size of steel tube and concrete and the dilation angle of concrete.
The dilation angle of concrete is the angle of internal friction and it controls the amount of plastic volumetric strain. For the definition of the concrete constitutive model, the dilation angle usually is selected based on the value that can represent the experimental response of one specimen and verify it for the other specimens [53]. The load-deflection results of specimens SA and SU-40 that are simulated with dilation angles of 22 • and 30 • are presented in Figure 12. It is found the circular column is more sensitive to the dilation angle value compared to the squared section. Figure 11. Load-deflection curves of the simulated models.

Sensitivity Analysis
The sensitivity of the results of finite element modeling to the input parameters is essential to study since it shows the effectiveness of the model and accordingly the outcomes. The main input parameters that are investigated for sensitivity are; the mesh size of steel tube and concrete and the dilation angle of concrete.
The dilation angle of concrete is the angle of internal friction and it controls the amount of plastic volumetric strain. For the definition of the concrete constitutive model, the dilation angle usually is selected based on the value that can represent the experimental response of one specimen and verify it for the other specimens [53]. The loaddeflection results of specimens SA and SU-40 that are simulated with dilation angles of 22° and 30° are presented in Figure 12. It is found the circular column is more sensitive to the dilation angle value compared to the squared section.
A mesh convergence study is performed where each column is simulated with different mesh sizes and the results are compared. Sufficiently small elements are used. Three mesh sizes are adopted to model specimens "SA" and "MA", namely; 8 mm, 15 mm, and 20 mm. Three mesh sizes are used to model specimens "SU-40" and "SU-70", namely; 15 mm, 20 mm, and 25 mm. The results of load-deflection of all considered columns at these different mesh sizes are introduced in Figure 13. It is found that the peak load results are not affected by the different mesh sizes. A mesh convergence study is performed where each column is simulated with different mesh sizes and the results are compared. Sufficiently small elements are used. Three mesh sizes are adopted to model specimens "SA" and "MA", namely; 8 mm, 15 mm, and 20 mm. Three mesh sizes are used to model specimens "SU-40" and "SU-70", namely; 15 mm, 20 mm, and 25 mm. The results of load-deflection of all considered columns at these different mesh sizes are introduced in Figure 13. It is found that the peak load results are not affected by the different mesh sizes.

Conclusions
Concrete-Filled Steel Tube (CFST) columns are composite columns that are getting significant attention as a replacement for reinforced concrete due to their ductility and lower weight. CFST consists of an outer steel tube that is filled with concrete. Various models available in the literature are limited to predicting the ultimate bond strength along the steel-concrete interface of CFSTs under push-out loading. Nevertheless, there isno prediction of the post-peak behavior in the literature. Accordingly, the behavior along the steel-concrete interface in composite columns is investigated in this paper to propose prediction models of the overall behavior which can be used for engineering practice.
The dataset considered for this investigation includes 157 circular specimens and 105 squared specimens from several references. These specimens had the following key parameters; diameter or width that ranges from 90.85 mm to 254 mm, thickness that ranges from 2.5 mm to 6.6 mm, length that ranges from 190 mm to 1498.6 mm, coreconcrete compressive strength that ranges from 9.11 MPa to 96.43 MPa, and age of concrete when the push-out test was performed that ranges from 28 days to 365 days.
ANN and ANOVA models are proposed in this study to obtain the overall stressslip behavior of the CFSTs. The ANN analysis is shown to be the most appropriate models for anticipating overall interfacial behavior among the major existing equations and the suggested ANOVA formulas herein. The ratios of the experimental to the predicted ultimate bond strength based on ANOVA and ANN formulas have mean values of 1.001

Conclusions
Concrete-Filled Steel Tube (CFST) columns are composite columns that are getting significant attention as a replacement for reinforced concrete due to their ductility and lower weight. CFST consists of an outer steel tube that is filled with concrete. Various models available in the literature are limited to predicting the ultimate bond strength along the steel-concrete interface of CFSTs under push-out loading. Nevertheless, there isno prediction of the post-peak behavior in the literature. Accordingly, the behavior along the steel-concrete interface in composite columns is investigated in this paper to propose prediction models of the overall behavior which can be used for engineering practice.
The dataset considered for this investigation includes 157 circular specimens and 105 squared specimens from several references. These specimens had the following key parameters; diameter or width that ranges from 90.85 mm to 254 mm, thickness that ranges from 2.5 mm to 6.6 mm, length that ranges from 190 mm to 1498.6 mm, core-concrete compressive strength that ranges from 9.11 MPa to 96.43 MPa, and age of concrete when the push-out test was performed that ranges from 28 days to 365 days.
ANN and ANOVA models are proposed in this study to obtain the overall stress-slip behavior of the CFSTs. The ANN analysis is shown to be the most appropriate models for anticipating overall interfacial behavior among the major existing equations and the suggested ANOVA formulas herein. The ratios of the experimental to the predicted ultimate bond strength based on ANOVA and ANN formulas have mean values of 1.001 and 0.999, respectively for circular CFSTs, and 1.067 and 1.004, respectively for squared CFSTs. These mean values represent how ultimate bond strength is well predicted compared to the mean values of the results estimated by the other existing formulas. Moreover, the proposed ANOVA and ANN formulas have a standard deviation of 0.277 and 0.166, respectively, for circular specimens and 0.596 and 0.188, respectively, for squared specimens compared to the standard deviation values of the results estimated by the existing models. Hence, ANOVA and ANN are more appropriate models to predict the ultimate bond strength for circular and squared CFSTs.
The proposed ANN models provide better prediction results than ANOVA due to the use of the sigmoidal activation function which allows the ANN to capture the non-linear mechanical behavior of concrete. ANN models have a lesser standard deviation of the ratios of the experimental to the predicted ultimate bond strength, slip at ultimate bond strength, and the parameters controlling the post-peak curve. R 2 values are found ranging from 0.674 to 0.876 for circular CFSTs and from 0.548 to 0.837 for squared CFSTs, which are high compared to ANOVA results. Consequently, it can be concluded that ANN models are appropriate to anticipate the full behavior of circular and square CFSTs and are better than ANOVA.
The difference between ANN results and experimental outcomes is considered reasonable since the test results of identical specimens had similar differences in their slip-bond response curves. This can be attributed to the variation of properties of concrete, manufacturing, and previously mentioned irregularities in internal forces.
This paper also investigates the ability of proposed ANN models to be used in numerical non-linear finite element models to predict the performance of CFST's columns under axial load. It is concluded that calibrated finite element models can represent the experimental results which may make it beneficial to further study CFST's columns. The Concrete Damage Plasticity model and cohesive-damage interface based on ANN models are sufficient and can be adopted to simulate the behavior of CFST's columns under axial load.
Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/buildings12040456/s1, Table S1. Details of circular CFSTs specimens and results obtained from experiments and the prediction from ANOVA and ANN; Table S2. Details of squared CFSTs specimens and results obtained from experiments and the prediction from ANOVA and ANN. Excel sheets of the ANN models is provided as supplementary documents.

Conflicts of Interest:
The authors declare no conflict of interest.