Abstract
The interface slip between shape steel and concrete is one of the biggest problems in numerical simulation of composite structures, this paper aims to investigate interfacial bond-slip behavior between steel and steel fiber reinforced concrete. For this purpose, a series of push-out tests on 15 circular section specimens which were designed and fabricated with different steel fiber ratio, interface embedded length and concrete cover thickness have been done. Experimental results show that average bond strength versus free end slip curves of each specimen is similar, and then a basic average bond strength-free end slip constitutive model was proposed in this paper. Meanwhile, considering the distribution regularities of steel strains, local bond strength and relative slip along the interface embedded length at different load levels were obtained through the mechanical derivation, the results show that the maximum local bond strength was located near the free end, and the local bond strength away from free end tended to be steady. Based on the basic average bond strength-free end slip constitutive model, two position functions were proposed to describe the local bond-slip constitutive model considering the variation of positions. Two position functions are of great significance for further research on SSFRC non-linear properties by FEM.
Keywords
Introduction
Steel reinforced concrete (SRC), which provides the advantages of high strength, high stiffness and pronounced seismic resistance (Tong et al., 2016, 2017, 2018), is in the family of composite structures. It has been widely applied in modern civil engineering (Elmy and Nakamura, 2017; Tsai et al., 2018), especially in high-rise buildings and areas of high earthquake requirement (Boyd et al., 1995; Hsu et al., 2004; Paultre et al., 2010; Shi et al., 2012; Wu et al., 2018). Choi et al. (2016) pointed out that more than half of high-rise buildings over 200 m were completed in 2014, which employed SRC composite structures. The paper of Xue et al. (2018) showed that SRC frame with irregular section columns is suitable for the construction in the high seismic intensity zone.
Choi et al. (2018) pointed out that the lower angles of inclination for the diagonal reinforcement, when accorded to the detailing requirements specified in ACI 318 (2014), result in reinforcement congestion and construction difficulties. The existence of both rebar cages and reinforcing steel causes the conflict between rebar cages and reinforcing steel in space (Ye et al., 2013). The space conflict leads to two construction problems: (1) inconvenience of forming rebar cages, which results in the discontinuity of rebars and stirrups and (2) difficulty of casting concrete, as shown in Figure 1(a), which leads to poor quality of concrete casting. Presently, the common methods used to solve these problems include opening at the flange or web for stirrups to close up (Figure 1(b)) (Yang et al., 2016) or using special connectors for fixing rebars to the shape steel (Figure 1(c)) (Nzabonimpa and Hong, 2018; Seo and Varma, 2017). However, such solutions not only increase construction difficulty and project investment, but also lower the integral bearing capacity of SRC composite structures. Construction difficulties of steel reinforced concrete. (a) Poor quality of concrete, (b) Openings at steel, and (c) Couplers.
Steel fiber reinforced concrete (SFRC) has been widely applied in buildings and bridges, which is made by mixing a certain ratio of steel fiber into plain concrete (Lim and Nawy, 2005; Yoo and Moon, 2018). Harajli (2006) studied the effect of confinement using steel fibers and indicated that confining the concrete with steel fiber reinforcement increased the bond strength and reduced the bond degradation. The experimental results of Chao et al. (2009) indicated that the bridging effect can effectively enhance post-cracking tensile capacity and limit crack width, which was provided by these irregular steel fibers in SFRC after cracking. Li et al. (2018) carried out an experimental investigation on the stress-strain behavior of SFRC subjected to uniaxial monotonic and cyclic tension, which showed that the determinant of concrete failure mode changes from tensile cracks to shear cracks with increasing fiber volume fraction. In order to solve the construction problems of SRC composite structures, achieving the goal of avoiding the inconvenient procedure of forming rebar cages and improving the pouring quality of concrete, the steel and steel fiber reinforced concrete (SSFRC) structures are adapted by using steel fibers instead of rebar cages. In the existing literature, many scholars have confirmed the feasibility of using steel fibers to replace shear reinforcement bars or even rebar cages in reinforced concrete structures (Boita et al., 2017; Dinh et al., 2010; Ding et al., 2011; Ghazizadeh and Cruz-Noguez, 2018; Hung and Hu, 2018; Valipour et al., 2016; Yoo et al., 2017)
The composite action of SRC is provided by the force transfer between steel and the concrete, which is obtained by relying on bonding properties (Tao and Yu, 2012). Roeder et al. (1999) evaluated bond stress in SRC, and came to the conclusion that bond stress is vital to the response of SRC structures, and shear connectors are required if the bond stress demand exceeds that capacity. Zheng et al. (2016) obtained local bond strength through the mechanical derivation. The paper of Wang et al. (2017) presented the results of bond strength in SRC columns after fire exposure, and the outcome showed that the reduction in bond strength was as high as 54.2% due to fire exposure. The experimental test carried out by Liu et al. (2018) established calculation formula of bond strength for steel reinforced recycled concrete. Further, the bridging effect provided by steel fibers enhances the post-crack tensile properties of concrete, and improves the toughness of concrete and the bonding properties (Tóth et al., 2019).
Similar to steel reinforced concrete, the basis to ensure the joint work of the two materials is the bonding properties between steel and SFRC (Hamad and Abou Haidar, 2016; Morishita and Tomii, 1982). The transfer of force across the interface between steel and SFRC plays a major role in controlling the response of SSFRC structures (Chen and Teng, 2001). The constitutive relation analysis is of great significance to numerical simulation of mechanical properties. The paper of Meng et al. (2016) pointed out that the constitutive model can reflect the bonding mechanism and provides references for theoretical analysis and engineering applications. Yuan et al. (2018) investigated the interfacial bond behavior between basalt fiber reinforced polymer sheet (BFRP) and SFRC. In addition, the analytical bond strength models and interfacial bond-slip models were proposed. Song et al. (2021) studied the interfacial bond-slip behavior between the corrugated steel plate and concrete, and proposed a bond-slip constitutive model of the corrugated steel and reinforced concrete. Constitutive model could provide a certain basis for establishing the relevant theoretical calculation methods. Yu et al. (2021) analyzed bond-slip behaviors self-stressing steel slag concrete filled steel tubular (SSSCFST) columns, a bond-slip constitutive model of SSSCFST columns was established, and the distribution law of the longitudinal compressive stress, inter-facial bond stress, and relative slip were obtained. The understanding of the bond properties between steel and SFRC is the basis for the capacity analysis of SSFRC structure.
Currently, design methods for reinforced concrete members are mainly based on the rectangular section. Mohammed et al. (2019) put forward a simplified new approach for designing RC specimens with circular section which converts circular section specimens to square section specimens with equivalent area, the test results showed that the RC specimens with circular section and its equivalent square specimens have similar nonlinear performance. The paper of Mousa et al. (2019) evaluated the applicability of crack-control models, which were based on rectangular section, to RC bridge members with a circular section, and it was considered that the formula can be applied to circular section after redefining parameters. Same as RC structure, there are few design methods of SRC based on circular section. Furthermore, buildings with traditional style become more and more popular, which mostly use circular section specimens. Therefore, 15 SSFRC specimens were fabricated in this paper to propose formulas based on circular section instead of rectangular section.
The experiment research conducted herein focused on evaluating the average bond strength and local bond strength along the interface at different positions of steel embedded in SFRC through the standard push-out tests. Based on the basic constitutive model of average bond strength versus free end slip, two position functions were proposed to establish the constitutive model of the local bond strength versus relative slip along the interface. Thus, it gives a complete description of the local bond strength at any given free end slip about different positions along the interface, which is of great significance for further research on SSFRC non-linear properties by FEM.
Experimental program
Materials
Steel and steel fibers
Information of steel fibers.
Concrete
Mix proportion of steel fiber reinforced concrete and compressive strength of concrete blocks.
*Note: FA-Fine aggregate; CA-Coarse aggregate.
Fabrication of specimens
Design parameters of specimens.
*Note: L is the length of shape steel.

Test setup and sketch of specimens. (a) Test setup and (b) Failure mode of specimens.
All the specimens were cast vertically. During the process of casting, the upper part and lower part of the shape steel extended the SFRC 50 and 20 mm, respectively. All specimens were demolded after 3 days, and then were moved to the curing room under natural conditions for 28 days.
Test setup and measurement scheme
A series of push-out tests of SSFRC were carried out by using the test setup as Figure 2 shown. The lower part of concrete supported on special loading plate, which bore the load applied by the TYA-2000 type electric hydraulic test machine, is called concrete loading end. During the test, the supporting load on the upper part was applied to the shape steel. The upper part of the concrete is guaranteed to stay in an unconstrained condition, which is called concrete free end. An H-typed hole, with slightly larger size compared with the inner shape steel, was dug in the middle of the special loading plate. Two displacement transducers were applied to measure the free end slip (S). The free end slip is defined as the mean value of the two displacement transducers. Due to the uneven distribution of steel strain along the interface embedded length, strain gauges were attached to the steel flange along the interface embedded length before casting. The layout of strain gauges is shown in Figure 3. Lay out of strain gauges.
Prior to the loading, the upper and lower surfaces of all specimens were pure smooth to provide a level surface and to ensure even distribution of the applied loads. In the process of loading, the load control was adopted until the ultimate load with the rate of 1 kN/s firstly, and then displacement control method was adopted with the rate of 0.01 mm/s (Wu et al., 2019). At the final stage of loading, end loading when load stayed constant or slip value reached 20 mm.
Experimental results and discussion
Average bond strength-free end slip curve
Figure 4 demonstrates the P-S curves for 15 specimens, from which we can see that the P-S curves are very similar. The entire P-S curves can be summarized into four stages, which are made of the elastic stage, yield stage, steep descending stage and gentle descending stage. During the test, the whole process of interface failure is very complex. Defining the yield point and yield load is helpful to judge the initial damage and initial relative slip of the interface. Geometric drawing method, equal energy method, and R. Park method were used to calculate the yield point on the load-slip curves at the free end. The yield load (Py) at the free end is defined as the average value of two or three methods which have similar results. Load-free end slip curves. (a) No. 1–8 specimens and (b) No. 9–15 specimens.
The average bond strength at the interface between steel and SFRC from test results is defined by equation below:
According to equation (1), the average bond strength Average bond strength-free end slip curves. (a) No. 1–8 specimens and (b) No. 9–15 specimens.
Average bond strength-free end slip
constitutive relationship
Obviously, the A typical 
In the model, the elastic stage (OA section), yield stage (AB section), steep descending stage (BC section), and gentle descending stage (CD section) were defined. Point A, B, C and D are the characteristic points of the
In the OA section, average bond strength increased elastically to
In the AB section, a quadratic relationship between average bond strength and free end slip was defined. In this stage, the hardening of average bond strength lasts until ultimate average bond strength
In the BC section, a hyperbolic relationship between average bond strength and free end slip was defined. In this stage, average bond strength descends steeply until
Values of characteristic points.
The average bond strength versus free end slip curve for SSFRC can be expressed as the equations below:
Statistical regression of characteristic points
Considering the steel fiber ratio, concrete cover thickness, interface embedded length, the concrete strength, and the statistical regression equations for characteristic average bond strength
The comparison of the test and fitting values is shown in Figure 7. It can be seen in Figure 7 that the fitting values are in good agreement with the test values. It needs to be mentioned that the specimen C2-4-20 was subjected to accidental eccentricity during the loading, which led to abnormal performance. Comparison of test values and fitting values of characteristic average bond strength. (a) Comparison of 
The characteristic free end slip is mainly affected by interface embedded length. As shown in Figure 8, through fitting the relationship between characteristic free end slip and interface embedded length, equations for characteristic free end slip Sy, Su, Sr, and Sf are proposed below: Test values and fitting curves of characteristic free end slip. (a) Fitting curve of Sy, (b) Fitting curve of Su, and (c) Fitting curve of Sr.

Note that according to the P-S curve and loading system, Sf is defined as 20 mm. The comparison of the test and fitting values is shown in Figure 9. It can be seen in Figure 9 that the fitting values are in good agreement with the test values. Comparison of test values and fitting values of characteristic free end slip. (a) Comparison of Sy, (b) Comparison of Su, and (c) Comparison of Sr.
Comparison of test curves and fitting curves
As the coordinates of characteristic points are determined, then according to equations (2)–(10), the fitting curves of Comparison of test curves and fitting curves of 
Local bond strength and relative slip
Distribution of steel strain along the interface embedded length
Figure 11 presents the steel strain distribution along interface embedded length at different load levels, which were obtained by strain gauges. Distance away from the free end is defined as x. According to Figure 11, it can be evidently seen that the steel strain distribution along the interface embedded length is close to the exponential relationship in the whole process of loading. Therefore, the steel strain distribution can be expressed by the exponential equation: Steel strain distribution along interface embedded length at different load levels. (a) C1-2-20, (b) C1-6-20, (c) C2-2-20, (d) C2-4-20, (e) C2-6-20, (f) C3-2-20, (g) C3-4-20, (h) C3-6-20, (j) C2-2-40, (k) C2-3-40, (l) C2-4-40, (m) C2-1-65, (n) C2-2-65, and (p) C2-3-65.

Distribution of local bond strength along the interface embedded length
In this study, the distribution of local bond strength was obtained through mechanical derivation. The applied load P satisfies the equilibrium equation below:
The distribution of local bond strength was calculated and shown in Figure 12. It is obvious that the local bond strength increases with the growing of applied load, and the maximum local bond strength is located near the free end. Local bond strength distribution along interface embedded length at different load levels. (a) C1-2-20, (b) C1-6-20, (c) C2-2-20, (d) C2-4-20, (e) C2-6-20, (f) C3-2-20, (g) C3-4-20, (h) C3-6-20, (j) C2-2-40, (k) C2-3-40, (l) C2-4-40, (m) C2-1-65, (n) C2-2-65, and (p) C2-3-65.
Distribution of relative slip along the interface embedded length
The relative slip ΔS between steel and SFRC can be expressed as follows:
Substituting equation (11) into equation (14) gives: Distribution of relative slip along interface embedded length at different load levels. (a) C1-2-20, (b) C1-6-20, (c) C2-2-20, (d) C2-4-20, (e) C2-6-20, (f) C3-2-20, (g) C3-4-20, (h) C3-6-20, (j) C2-2-40, (k) C2-3-40, (l) C2-4-40, (m) C2-1-65, (n) C2-2-65, and (p) C2-3-65. Distribution of relative slip along interface embedded length at ultimate load (a) C2-2-20, (b) C2-4-20, and (c) C2-6-20.


Local bond-slip relationship
Local bond-slip relationship of different positions τ= τ(ΔS, x)
The typical Relationship of local bond strength versus relative slip at different positions. (a) C1-2-20, (b) C1-6-20, (c) C2-2-20, (d) C2-4-20, (e) C2-6-20, (f) C3-2-20, (g) C3-4-20, (h) C3-6-20, (j) C2-2-40, (k) C2-3-40, (l) C2-4-40, (m) C2-1-65, (n) C2-2-65, and (p) C2-3-65.
Comparing the features of the set of curves τ= τ(ΔS, x) and
The determination of G(x) and F(x)
The determination of G(x)
Through the statistical regression analysis about the average value of τu(x)/ Determination of G(x). (a) Statistical regression of G(x) and (b) Fitting result of G(x).
The determination of F(x)
By the same procedure as above, position function F(x) can be acquired through the statistical regression analysis about the average value of ΔSu(x)/Su in the position x for 14 specimens. Figure 17(a) shows the procedure of statistical regression analysis, which indicated that a cubic polynomial has the highest degree of fitting. The relevant fitting results are shown in Figure 17(b). Determination of F(x). (a) Statistical regression of G(x) and (b) Fitting result of G(x).
The local bond-slip constitutive model considering the variation of positions
On the basis of the typical The local bond-slip constitutive model.

The characteristic average bond strength
Conclusions
The interfacial local bond strength between steel and SFRC has been investigated in this paper. A total of 15 circular section SSFRC were tested under monotonic load. The following conclusions can be drawn from this paper: (i) There are four characteristic points of the (ii) The fitting results of (iii) It can be seen that the distribution regularities of steel strain along the interface embedded length is close to the exponential relationship in the whole process of loading. According to this and the mechanical derivation, the distributions of steel strain, local bond strength and relative slip along the interface embedded length at different positions were acquired. (iv) The distribution of local bond strength is uneven along the interface embedded length, the maximum local bond strength is located near the free end, and the local bond strength away from free end tends to be steady. (v) Comparing the features of the set of curves τ= τ(ΔS, x) and
Footnotes
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work is supported by the National Natural Science Foundation of China (Grant No.51208175) and Fundamental Research Funds for the Central Universities (Grant No. B200202067 and Grant No. B210203021).
