Effect of Low-Level Cyclic Loading on Bond Behavior of a Steel Bar in Concrete with Pre-Existing Damage

Understanding the bond behavior of steel rebar in concrete is important in order to determine the performance of a reinforced concrete structure. Although numerous studies have been carried out by many researchers to develop a robust model for numerical analysis, no consensus has been reached as the bond behavior depends on hysteresis. In this study, the bond behavior of a steel bar in concrete with pre-existing damage is investigated under low-level cyclic loading. Based on the experimental bond stress and slip curve, a numerical model for finite element analysis to simulate the effect of low-level cyclic loading is proposed. The results from the numerical analysis show good agreement with the experimental data, including accumulated damage on stiffness and strength throughout entire load cycles.


Introduction
Bonding behavior of steel bars in concrete has been studied by several researchers in recent decades. In 1960s and 1970s, most of experimental studies focused on the bond behavior under monotonic loading with different factors as variables: rebar deformation [1], boundary condition [2], deformation dimensions [3], and cover thickness and concrete strength [4]. Analytical models needed to use refined numerical method to determine reinforced concrete behavior were also reported [4,5]. The effect of hysteresis on the bond behavior of a steel bar in concrete was first reported in the early 1970s [6], and a number of studies reported in the 1980s [7][8][9] laid groundwork for a numerical model. However, there has been no single model that reflects all factors that affect the bond stress-slip behavior of steel rebar in concrete due to difficulties in measuring bond stress and bond slip, and unpredictable hysteresis and boundary conditions. Under cyclic loading, the degradation of bonding behavior is more severe than that under monotonic cyclic loading. The degradation becomes more significant with larger slip value at which reversal of loading occurs. Ismail and Jirsa [7] found that application of load cycles with constant peak stress can yield gradual deterioration in bonding even with a small number of repetitions. Eligehausen et al. [8] carried out an experimental study to investigate the bond stress-slip relationship under cyclic loading with an extensive list of variables. They showed that the envelope of bond stress-slip hysteresis loops for repeated loading lies below the curve from monotonic loading, as shown in Figure 1. Gan listed the most important parameters under cyclic loading as follows: strain (or stress) range, type and rate of loading (strain rate), and the value of maximum imposed bond stress [10]. Based on the test results, the bond stress-slip models were proposed according to two basic approaches: (1) microscopic approach reflecting the local stress state between concrete Due to low-level cyclic loading, such as the moving service loads or small earthquakes, a cumulative damage to the bond between steel rebar and concrete can occur. Since this accumulated damage can compromise the seismic performance of reinforced concrete structures, its effect on bond behavior of steel bar in concrete needs to be understood. In the present study, this effect is investigated by constant amplitude cyclic bond tests under displacement control with a fixed number of loading cycles. In addition, a numerical model for finite element analysis to simulate the effect of low-level cyclic loading is derived based on the experimental result.

Materials
The steel rebar used in this study was manufactured by a local company (Dongkuk Steel Mill Co., Ltd., Seoul, Korea) according to specifications provided in Korean Industrial Standard KS D 3504 [16]. Figure 2 shows a sample of the rebar. The nominal diameter of the rebar was 19 mm and graded as SD400 (minimum yield strength of 400 MPa, minimum tensile strength of 560 MPa). The chemical composition of steel bar was tested according to KS D1801 [17]. Table 1 shows the analysis results, which confirms that the chemical composition of the steel bar satisfies the specifications per KS D 3504. Concrete was prepared by a local ready mixed concrete manufacturer. The strength of concrete was designed as 27 MPa, the slump was 120 mm, the water-cement ratio was 48.2%, the maximum size of aggregate was 25 mm, and the air volume was 4.5%. The concrete composition is shown as Table 2. Due to low-level cyclic loading, such as the moving service loads or small earthquakes, a cumulative damage to the bond between steel rebar and concrete can occur. Since this accumulated damage can compromise the seismic performance of reinforced concrete structures, its effect on bond behavior of steel bar in concrete needs to be understood. In the present study, this effect is investigated by constant amplitude cyclic bond tests under displacement control with a fixed number of loading cycles. In addition, a numerical model for finite element analysis to simulate the effect of low-level cyclic loading is derived based on the experimental result.

Materials
The steel rebar used in this study was manufactured by a local company (Dongkuk Steel Mill Co., Ltd., Seoul, Korea) according to specifications provided in Korean Industrial Standard KS D 3504 [16]. Figure 2 shows a sample of the rebar. The nominal diameter of the rebar was 19 mm and graded as SD400 (minimum yield strength of 400 MPa, minimum tensile strength of 560 MPa). The chemical composition of steel bar was tested according to KS D1801 [17]. Table 1 shows the analysis results, which confirms that the chemical composition of the steel bar satisfies the specifications per KS D 3504. Concrete was prepared by a local ready mixed concrete manufacturer. The strength of concrete was designed as 27 MPa, the slump was 120 mm, the water-cement ratio was 48.2%, the maximum size of aggregate was 25 mm, and the air volume was 4.5%. The concrete composition is shown as Table 2.

Specimen Preparation
As shown in Figure 3, test specimens were designed as a cube having sides dimension equal to 15 times the diameter of the bar (15db, 300 mm) to avoid splitting failure with embedment length of 5db, based on the study by Casanova et al. [18]. Steel bars were cut to equal length of 1800 mm. The length of the steel bar was determined in order to acquire room for LVDT (Linear Variable Differential Transformer) installation and screw thread for mechanical coupler, which was adapted to anchor the rebar during the test. Threads were made at both ends of the steel bar over the length of 60 mm according to the specification set by the coupler manufacturer (Jungwoo B&C Co., Ltd., Paju, Korea).

Specimen Preparation
As shown in Figure 3, test specimens were designed as a cube having sides dimension equal to 15 times the diameter of the bar (15d b , 300 mm) to avoid splitting failure with embedment length of 5d b , based on the study by Casanova et al. [18].

Specimen Preparation
As shown in Figure 3, test specimens were designed as a cube having sides dimension equal to 15 times the diameter of the bar (15db, 300 mm) to avoid splitting failure with embedment length of 5db, based on the study by Casanova et al. [18]. Steel bars were cut to equal length of 1800 mm. The length of the steel bar was determined in order to acquire room for LVDT (Linear Variable Differential Transformer) installation and screw thread for mechanical coupler, which was adapted to anchor the rebar during the test. Threads were made at both ends of the steel bar over the length of 60 mm according to the specification set by the coupler manufacturer (Jungwoo B&C Co., Ltd., Paju, Korea). Steel bars were cut to equal length of 1800 mm. The length of the steel bar was determined in order to acquire room for LVDT (Linear Variable Differential Transformer) installation and screw thread for mechanical coupler, which was adapted to anchor the rebar during the test. Threads were made at both ends of the steel bar over the length of The embedded length (5d b ) was secured in the middle of the reinforcement guided by two PVC pipes of which openings were later sealed with silicone to prevent intrusion of concrete into the pipe during concrete placement, as shown in Figure 3. Other parts of reinforcement were sprayed with WD-40 (Water Displacement 40th formula by WD-40 Company, San Diego, CA, USA) to prevent corrosion. After the concrete placement, all specimens were cured in moist condition for up to 28 days.

Cyclic Load Test Method
The testing machine shown in Figure 4 was used for cyclic load test. The loading device was designed to exert the load on the specimen in two different modes: displacement control or load control. The double acting hydraulic cylinder with 200-tonf capacity at pressure of 700 bar was manufactured by a local company (TDC-20015 by Daejin Hydraulic Machinery Co., Ltd., Busan, Korea). A low-profile load cell with a 20-tonf compression and tension measurement range, manufactured by a local company (Bongshin Loadcell Co., Ltd., Osan, Korea), was used to measure applied load. Two linear displacement transducers (CDP-50M by Tokyo Measuring Instruments Lab Co., Ltd., Tokyo, Japan) were used to measure local slip as shown as LVDT2 and LVDT3 in Figure 4. An additional linear displacement transducer, LVDT1, was installed on the base plate to monitor the global displacement of the specimen holder. All measurements were recorded using a data logger (UCAM-60B by Kyowa Electronic Instruments Co., Ltd., Tokyo, Japan) at one second intervals. The embedded length (5db) was secured in the middle of the reinforcement guided by two PVC pipes of which openings were later sealed with silicone to prevent intrusion of concrete into the pipe during concrete placement, as shown in Figure 3. Other parts of reinforcement were sprayed with WD-40 (Water Displacement 40 th formula by WD-40 Company, San Diego, CA, USA) to prevent corrosion. After the concrete placement, all specimens were cured in moist condition for up to 28 days.

Cyclic Load Test Method
The testing machine shown in Figure 4 was used for cyclic load test. The loading device was designed to exert the load on the specimen in two different modes: displacement control or load control. The double acting hydraulic cylinder with 200-tonf capacity at pressure of 700 bar was manufactured by a local company (TDC-20015 by Daejin Hydraulic Machinery Co., Ltd., Busan, Korea). A low-profile load cell with a 20-tonf compression and tension measurement range, manufactured by a local company (Bongshin Loadcell Co., Ltd., Osan, Korea), was used to measure applied load. Two linear displacement transducers (CDP-50M by Tokyo Measuring Instruments Lab Co., Ltd., Tokyo, Japan) were used to measure local slip as shown as LVDT2 and LVDT3 in Figure 4. An additional linear displacement transducer, LVDT1, was installed on the base plate to monitor the global displacement of the specimen holder. All measurements were recorded using a data logger (UCAM-60B by Kyowa Electronic Instruments Co., Ltd., Tokyo, Japan) at one second intervals. After placing a specimen inside the specimen holder by hand, each specimen was positioned and secured as follows: (a) thin steel sheets (0.5 mm, 1 mm, 2 mm thickness) were used as shim plate to secure the specimen in the middle of the holder as needed; (b) the position of the specimen was then adjusted using the hydraulic cylinder so that the specimen sits at the center between the reaction supports; (c) the initial readings from the LVDT and the load cell were recorded, and the steel bar was fastened with a set of semispherical washer and a rebar coupler at each end. The load was monitored during the fastening so that it remained consistent with the initial load reading.
A monotonic pullout test was performed first by applying displacement load at the rate of 0.03 mm/s until the average local slip value reached 20 mm in order to determine the low-level cyclic loading parameters. From the monotonic pullout test, the average local slip corresponding to the maximum bond stress, Sm, was found to be 1.67 mm. Based on the Sm obtained from the monotonic test, the cyclic load level for this study was determined as (±)0.5 mm, which equals to (±)0.3 × Sm.
Before the commencement of the cyclic load test with the (±)0.5 mm load level, in order to emulate the pre-existing damage in a reinforced concrete structure, the specimen was subjected to three cycles of load in the range of (−)0.5~(+)0.5 tonf, which showed no measurable local displacement, followed by four cycles in the range of (±)1 tonf which resulted in a discernable local displacement. Afterward, the (±)0.5 mm load was applied After placing a specimen inside the specimen holder by hand, each specimen was positioned and secured as follows: (a) thin steel sheets (0.5 mm, 1 mm, 2 mm thickness) were used as shim plate to secure the specimen in the middle of the holder as needed; (b) the position of the specimen was then adjusted using the hydraulic cylinder so that the specimen sits at the center between the reaction supports; (c) the initial readings from the LVDT and the load cell were recorded, and the steel bar was fastened with a set of semi-spherical washer and a rebar coupler at each end. The load was monitored during the fastening so that it remained consistent with the initial load reading.
A monotonic pullout test was performed first by applying displacement load at the rate of 0.03 mm/s until the average local slip value reached 20 mm in order to determine the low-level cyclic loading parameters. From the monotonic pullout test, the average local slip corresponding to the maximum bond stress, S m , was found to be 1.67 mm. Based on the S m obtained from the monotonic test, the cyclic load level for this study was determined as (±)0.5 mm, which equals to (±)0.3 × S m .
Before the commencement of the cyclic load test with the (±)0.5 mm load level, in order to emulate the pre-existing damage in a reinforced concrete structure, the specimen was subjected to three cycles of load in the range of (−)0.5~(+)0.5 tonf, which showed no measurable local displacement, followed by four cycles in the range of (±)1 tonf which resulted in a discernable local displacement. Afterward, the (±)0.5 mm load was applied 10 times under the displacement control, at a rate of 0.03 mm/s, to record bond-slip behavior. Upon completion of 10 load cycles, the pullout behavior of the steel bar was recorded up to the local slip value of 20 mm in the same way as for the monotonic pullout test specimen as mentioned above. The pullout result is shown in Figure 5. 10 times under the displacement control, at a rate of 0.03 mm/s, to record bond-slip behavior. Upon completion of 10 load cycles, the pullout behavior of the steel bar was recorded up to the local slip value of 20 mm in the same way as for the monotonic pullout test specimen as mentioned above. The pullout result is shown in Figure 5.  Figure 6 shows that the ascending branch of the first cycle begins with a shallow slope and the slope sharply increases with the inflection point at the bond stress value of 4.5 MPa. The ascending branch softens gradually and reaches the maximum bond stress of 8.3 MPa at the slip of 0.5 mm. Upon unloading, the load drops rapidly with a little change in slip with the slope of 200 MPa/mm. As the displacement load increases in reverse direction, bond stress increases very slowly up to (−)1.03 MPa, at which point a plateau is observed until the slip value reached zero. It is noted that the plateau is consistent with a frictional resistance; with increase in the displacement, bond stress increases with the inflection point at about (−)3 MPa. At the slip value of (−)0.5 mm, the bond stress reached (−)7.7 MPa, which is 6% less than the maximum bond stress at (+)0.5 mm. As the displacement direction is reversed, the load drops to zero with 0.1 mm change in slip and exhibits a new plateau at bond stress value of (+)0.8 MPa until the slip reaches zero. The reduction in frictional resistance and the slope of the descending branch of the curve in the negative slip range can be attributed to the cumulated damage to the bond between the concrete and steel.  Figure 6 shows that the ascending branch of the first cycle begins with a shallow slope and the slope sharply increases with the inflection point at the bond stress value of 4.5 MPa. The ascending branch softens gradually and reaches the maximum bond stress of 8.3 MPa at the slip of 0.5 mm. Upon unloading, the load drops rapidly with a little change in slip with the slope of 200 MPa/mm. As the displacement load increases in reverse direction, bond stress increases very slowly up to (−)1.03 MPa, at which point a plateau is observed until the slip value reached zero. It is noted that the plateau is consistent with a frictional resistance; with increase in the displacement, bond stress increases with the inflection point at about (−)3 MPa. At the slip value of (−)0.5 mm, the bond stress reached (−)7.7 MPa, which is 6% less than the maximum bond stress at (+)0.5 mm. As the displacement direction is reversed, the load drops to zero with 0.1 mm change in slip and exhibits a new plateau at bond stress value of (+)0.8 MPa until the slip reaches zero. The reduction in frictional resistance and the slope of the descending branch of the curve in the negative slip range can be attributed to the cumulated damage to the bond between the concrete and steel. Materials 2021, 14, x FOR PEER REVIEW 6 of 14 The general shape of the descending branch is retained for the subsequent load cycles (second to tenth), though it is noted that there is no inflection point in the ascending branch. The absence of the inflection point in the ascending branch suggests that the slipresistant mechanism is consolidated to a single mode after the first load cycle. As can be seen in Figure 7a, the maximum bond stress (τm) and the bond stress corresponding to the frictional resistance (τf) show significant degradation (41% and 30%, respectively) in the first three cycles, and approach asymptotic values with final degradation of 57% and 41% at the tenth cycle, with respect to those of the first cycle. Furthermore, it is noted that frictional resistance is always lower in positive direction than in negative direction by 21% on average. Figure 7b also shows asymmetry of the frictional plateau. The directionality of the bond-slip behavior seems to be an inevitable phenomenon as the degradation is not likely to occur uniformly over the lug of the deformed bar.  The general shape of the descending branch is retained for the subsequent load cycles (second to tenth), though it is noted that there is no inflection point in the ascending branch. The absence of the inflection point in the ascending branch suggests that the slip-resistant mechanism is consolidated to a single mode after the first load cycle. As can be seen in Figure 7a, the maximum bond stress (τ m ) and the bond stress corresponding to the frictional resistance (τ f ) show significant degradation (41% and 30%, respectively) in the first three cycles, and approach asymptotic values with final degradation of 57% and 41% at the tenth cycle, with respect to those of the first cycle. Furthermore, it is noted that frictional resistance is always lower in positive direction than in negative direction by 21% on average. Figure 7b also shows asymmetry of the frictional plateau. The directionality of the bond-slip behavior seems to be an inevitable phenomenon as the degradation is not likely to occur uniformly over the lug of the deformed bar. The general shape of the descending branch is retained for the subsequent load cycles (second to tenth), though it is noted that there is no inflection point in the ascending branch. The absence of the inflection point in the ascending branch suggests that the slipresistant mechanism is consolidated to a single mode after the first load cycle. As can be seen in Figure 7a, the maximum bond stress (τm) and the bond stress corresponding to the frictional resistance (τf) show significant degradation (41% and 30%, respectively) in the first three cycles, and approach asymptotic values with final degradation of 57% and 41% at the tenth cycle, with respect to those of the first cycle. Furthermore, it is noted that frictional resistance is always lower in positive direction than in negative direction by 21% on average. Figure 7b also shows asymmetry of the frictional plateau. The directionality of the bond-slip behavior seems to be an inevitable phenomenon as the degradation is not likely to occur uniformly over the lug of the deformed bar.  The effect of cyclic loading on the pullout behavior of the rebar is shown in Figure 5. The cyclic loading reduced the maximum bond stress and corresponding slip by 19% and 8%, respectively, compared to those of the specimen subjected to a simple monotonic loading. It is noted that two curves become essentially identical beyond a slip value of 7 mm, corresponding to the frictional resistance branch. The area under the bond stressslip curve up to the maximum bond stress was 41% smaller in the cyclic loading case than that of the monotonic case. Meanwhile, the area under the entire curve was 14% smaller in the cyclic loading case than that of the monotonic case, indicating that the low-level cyclic loading affects the pre-peak bond-slip behavior more than the post-peak.

Experimental Data-Driven Model
The experimental data obtained in this study suggest no model currently available in the literature would not properly capture the bond-slip behavior of a rebar in concrete subjected to low-level cyclic loading with an existing damage; existing numerical bond-slip models consider strain softening due to large crack formation around rebars. Additionally, initial stiffness degradation due to the pre-existing damage on bond capacity must be addressed in the numerical model. In this section, a new numerical model is suggested for low-level cyclic bond-slip behavior.

Simplification of Low-Level Cyclic Bond-Slip Behavior
The experimental cyclic curves are simplified as in Figure 8 for the suggested model, where the nonlinear loading curve is characterized by two control points (∆ 0 , θ 0 ) and (∆ tr , θ ∆ ). The first control point, ∆ 0 , represents the damaged bond range limit in which a reinforcing steel rebar slides with the minimum bond stress, σ 0 , and is evolved with loading cycles. From the experimental results, it can be reasonably assumed that disturbed and residual regions near the rib of a steel bar exist during each cycle stage as depicted in Figure 9. It is the disturbed region where steel bars slip backward and forward in the damaged region with the minimum bond stress, while the nonlinear loading behavior is due to the onset of the residual region. The peak stress, σ ∆ , at each cycle is also evolved and gradually decreased to represent the degradation on bond strength and stiffness. The bond stress is unloaded as the elastic unloads with the stiffness of E u . It is noted that the experimental results show the unloading stiffnesses for the pulling direction, which is the initial loading direction, and the opposite (pushing) direction are different. The unloaded bond stress is bounded by the minimum bond stress, σ 0 , in the opposite direction and keeps constant until the slip displacement reaches ∆ 0 in the opposite direction to represent the slip with the minimum bond stress in the disturbed region.

Formulation of Bond-Slip Model for Unloading/Reloading
The loading curve is suggested in the form of a second-order polynomial to represent both the pre-damaged initial and subsequential cyclic bond-slip behavior: where three coefficients are determined from the slopes at two control points (∆ 0 , θ 0 ) and (∆ tr , θ ∆ ) in Figure 8. The coefficients are determined for the initial cycle stage as: where ∆ 0 = 0, and two stiffnesses θ 0 and θ ∆ are determined from the experimental data. During the initial cycle, stiffness degradation appears on the loading curve after the stiffness reaches θ ∆ . In the present study, such degradation is modeled as the hypoelastic stress formulation and Equation (1) is replaced if dσ/d∆ > θ ∆ :

Formulation of Bond-Slip Model for Unloading/Reloading
The loading curve is suggested in the form of a second-order polynomial to represent both the pre-damaged initial and subsequential cyclic bond-slip behavior: where three coefficients are determined from the slopes at two control points (∆ , ) and (∆ , ∆ ) in Figure 8. The coefficients are determined for the initial cycle stage as:

Formulation of Bond-Slip Model for Unloading/Reloading
The loading curve is suggested in the form of a second-order polynomial to represent both the pre-damaged initial and subsequential cyclic bond-slip behavior: where three coefficients are determined from the slopes at two control points (∆ 0 , 0 ) and (∆ , ∆ ) in Figure 8. The coefficients are determined for the initial cycle stage as: .
where ∆ lim is the slip displacement limit for the present model, after which the bond strength drops rapidly and the low-level cyclic loading assumption is not valid. After the first cycle the coefficients are differently set to consider the slip with the minimum bond stress in the disturbed region and degradation on strength: where the parameter, ∆ 0 , and the peak bond stress σ ∆ are evolved with loading cycles and extracted from experimental results.
In the present study, the growth-saturation formulation is suggested to be used for the evolution of the damaged bond range limit parameter: where N is the cycle number, which is counted if the slip displacement exceeds ∆ 0 in each direction. The three parameters, γ 0 , α 1 and β 1 in Equation (4) can be calibrated from the experimental data. A similar form of the growth-saturation formulation is also used to represent the monotonically decreasing evolution of the peak stress due to the cyclic degradation of the bond strength: where σ max is the maximum bond stress in the previous loading cycles and N is set to one whenever the maximum is newly updated. The parameters, α 2 and β 2 in Equation (5) are also calibrated from the experimental data. For unloading, an elastoplastic-type incremental bilinear model is used: where . σ is the bond stress increment and .
∆ is the slip displacement increment. The unloading stiffness can be defined asymmetrically for each slip direction from the observation of the experimental results:

Numerical Simulation and Model Validation
The parameters in the suggested numerical model are determined from the experimental results in Section 2. To validate the numerical model, the experiment is numerically simulated. The finite element mesh and boundary conditions are depicted in Figure 10. Plane stress elements are used for ambient concrete and truss elements are used for steel reinforcing bars. For modeling presently developed bond-slip behavior, zero-length link elements are used. Since the concrete material and nonlinear contact behavior in bond-slip regions under a small amount of cyclic slip displacements is readily considered in the numerical model, the concrete is modeled as an elastic one. The slip displacement limit parameter in Equation (2) is ∆ lim = 0.6 mm. From the experimental data analysis, γ 0 = 0.3, α 1 = 0.6, and β 1 = 2.0 in Equation (4), and α 2 = 0.6 and β 2 = 1.85 in Equation (5)  bond-slip regions under a small amount of cyclic slip displacements is readily considered in the numerical model, the concrete is modeled as an elastic one. The slip displacement limit parameter in Equation (2)   Full 10 cycles as in the experiment are simulated and the numerical result is plotted in Figure 11. It can be observed that the result appropriately simulates the experimental bond-slip behavior discussed in Section 2. To evaluate the reproduction capability of the presently suggested numerical model, the first two cycles are closely compared in Figure  12. It is shown that the numerical model realistically represents the pre-damaged bond behavior by providing the initial damage on stiffness. It is also verified that the evolution of the damaged bond-slip boundary, ∆ , and the peak bond stress, ∆ , during the second loading path is appropriately reproduced. Full 10 cycles as in the experiment are simulated and the numerical result is plotted in Figure 11. It can be observed that the result appropriately simulates the experimental bond-slip behavior discussed in Section 2. To evaluate the reproduction capability of the presently suggested numerical model, the first two cycles are closely compared in Figure 12. It is shown that the numerical model realistically represents the pre-damaged bond behavior by providing the initial damage on stiffness. It is also verified that the evolution of the damaged bond-slip boundary, ∆ 0 , and the peak bond stress, σ ∆ , during the second loading path is appropriately reproduced.   To check the validation of the numerical model, relative errors between the experimental, and numerical results on the peak bond stresses in the pulling direction are evaluated and listed in Table 3. The maximum relative error is 7.98% at the sixth cycle and the To check the validation of the numerical model, relative errors between the experimental, and numerical results on the peak bond stresses in the pulling direction are evaluated and listed in Table 3. The maximum relative error is 7.98% at the sixth cycle and the average relative error for all the 10 cycles is 3.21%, showing that the suggested numerical bond-slip model simulates the experiment within the acceptable accuracy. The overall results show that the presently suggested numerical bond-slip model realistically reproduces the experimental results and can be used in simulating low-level cyclic bond-slip behavior with the pre-existing damage on bond capacity due to cyclic service loading.

Discussion
Based on the current study, it was found that cyclic loads, which are well below the elastic limit, can cause a significant change in the initial part of the bond stress-slip curve. Such initially degraded bond behavior is especially important in low-to-moderate intensity earthquake analysis and reflected in the proposed numerical model, which is suited for the analysis of old reinforced concrete structures. The following conclusions can be derived: (1) The low-level cyclic load, (±)0.3 × S m can cause significant degradation to steel rebar-concrete bond in the first three cycles, but the rate of degradation diminishes rapidly in subsequent cycles. (2) Damage caused by cyclic loading affects the pullout behavior of the steel rebar.
Measurable difference can be seen in the maximum bond stress, corresponding slip, and slope of descending branch. But the effect of the cyclic loading on the frictional resistance was negligible. (3) A new numerical model has been developed to simulate the effect of low-level cyclic load on the bond-slip behavior. The results obtained from the simulation showed good agreement with that from the experiment.

Implications and Limitations
Bond behavior of steel reinforcing bar in concrete has been studied by several researchers; however, to the best of the authors' knowledge, no attempts have been made to consider the effect of pre-existing damage in the bond. The implications and limitations of this study can be drawn as follows: (1) Since the accumulated damage to the bond between steel rebar and concrete, due to low-level cyclic loading, such as the moving service loads or small earthquakes, can compromise the seismic performance of reinforced concrete structures, its effect on bond behavior of steel bar in concrete need to be understood. In the present experimental study this pre-existing damage effect was verified, and it was discovered that the initial loading behavior is significantly different from the conventional ones. (2) The developed numerical bond-slip model can be embedded into a nonlinear finite element analysis procedure to simulate the realistic effect of low-level cyclic loading with the pre-existing bond stiffness degradation. This can contribute to enhancing the accuracy of performance evaluation of existing deteriorated reinforced concrete structures subjected to low-to-moderate intensity earthquakes. (3) In order to refine the proposed numerical model, further experimental studies are needed to investigate the effects of other parameters, such as rebar size and concrete strength, pertaining to hysteresis effects. Additionally, the proposed numerical model is not applicable to predict structural behavior beyond low-level cyclic loads, as the slip range considered in the model is limited up to the peak bond strength.
Consequently, the softening of bond behavior cannot be represented, and should be addressed in the future development of a high-level cycle model. Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.
Data Availability Statement: Data sharing is not applicable to this article.