The Effects of Bond–Slip Laws on the Debonding Failure and Behavior of Flexural Strengthened RC Slabs in Hybrid FRP Retrofit Systems

The hybrid retrofit system using FRP and concrete overlay applied on the top of slabs has proven effective in strengthening and overcoming logistical constraints, compared with conventional strengthening techniques using externally bonded composite materials to the underside of the slabs. Nevertheless, the performance of retrofitted slabs is governed by debonding failure due to the low bond strength between CFRP and concrete overlay. Thus, this study investigates the behavior of flexural strengthened slabs with FRP retrofit systems and the effect of bond–slip laws on debonding failure. Firstly, two full-scale RC slabs with and without a retrofit system were tested in a four-point bending setup as the control specimens. Then, the same retrofitted slab was simulated by utilizing the commercial program ABAQUS. A sensitivity analysis was conducted to consider the influence of bond–slip laws to predict the failure mechanism of the retrofitted slabs based on load–deflection relationships. The results showed that the strengthened slab enhanced the load-carrying capacity by 59%, stiffness by 111%, and toughness by 29%. The initial stiffness of 0.1K0 and maximum shear stress of 0.13τmax, compared with the corresponding values of Neubauer’s and Rostasy’s bond–slip law, can be used to simulate the global response of the retrofitted slab validated by experiment results.


Introduction
Strengthening structures with external bonded materials has become one of the popular choices for rehabilitation and upgrading existing reinforced concrete (RC) structures [1,2]. Researchers and engineers are particularly interested in restoring and strengthening partially damaged RC structures to minimize the impact on environmental deterioration and cost savings [3][4][5]. Fiber-reinforced polymer (FRP) materials have been widely used in the field of external reinforcement in recent decades because of their lightweight, high tensile strength, and non-corrosion resistance [6][7][8][9][10][11][12]. New strengthened techniques and design guidelines have emerged for RC structures using FRP composite materials [13][14][15][16].
The conventional application in strengthening RC slab techniques focuses on attaching FRP to their tensile zone to maximize the tensile strength of these composite materials, which have been widely gaining acceptance in practice [17][18][19]. However, the premature delamination of FRP remains the main disadvantage of the strengthening techniques, as it prevents strengthened slabs from reaching their ultimate load-carrying capacity. In general, sectional and debonding failures are two widespread possible failure modes of FRP-strengthened slabs [20][21][22]. In sectional failure, either the failure of compressive concrete or the rupture of FRP can be predicted through strain compatibility at its limiting strain [23]. In debonding failure, it is more difficult to determine due to the complicated failure mechanism and the factors contributing to problems such as concrete cracks and stress concentration at the FRP-concrete interface. Due to this, the behavior of externally bonded FRP in contact with a concrete surface has attracted considerable interest from the civil engineering community [24][25][26]. In addition, obtaining a well-prepared concrete surface for strengthening FRP on the underside of the RC slab can be a more difficult challenge, or it may not be possible due to logistical constraints.
Consequently, a hybrid retrofit system consists of carbon fiber reinforced polymer (CFRP) and concrete overlay applied on the top surface of the RC slab to overcome the shortcomings of traditionally strengthened solutions rather than exploit the high tensile strength of CFRP, as shown in Figure 1. Previous studies have demonstrated the effectiveness of the retrofit system by increasing the flexural strength and ductility of existing RC slabs. Nevertheless, the results also revealed that debonding failure between the concrete overlay and CFRP was one of the major causes of slabs losing their ultimate load-carrying capacity. Transferring interfacial stress from the slab to FRP or vice versa via a specialized adhesive layer is more efficient than from FRP to overlay members due to the characteristics of the retrofit system, leading to the premature debonding failure of the FRP/overlay interface [27][28][29]. A sound understanding of the behavior of FRP-concrete interface needs to be developed, with a particular interest in FRP/overlay interface for safe in practical design work. strengthened slabs [20][21][22]. In sectional failure, either the failure of compressive concrete or the rupture of FRP can be predicted through strain compatibility at its limiting strain [23]. In debonding failure, it is more difficult to determine due to the complicated failure mechanism and the factors contributing to problems such as concrete cracks and stress concentration at the FRP-concrete interface. Due to this, the behavior of externally bonded FRP in contact with a concrete surface has attracted considerable interest from the civil engineering community [24][25][26]. In addition, obtaining a well-prepared concrete surface for strengthening FRP on the underside of the RC slab can be a more difficult challenge, or it may not be possible due to logistical constraints. Consequently, a hybrid retrofit system consists of carbon fiber reinforced polymer (CFRP) and concrete overlay applied on the top surface of the RC slab to overcome the shortcomings of traditionally strengthened solutions rather than exploit the high tensile strength of CFRP, as shown in Figure 1. Previous studies have demonstrated the effectiveness of the retrofit system by increasing the flexural strength and ductility of existing RC slabs. Nevertheless, the results also revealed that debonding failure between the concrete overlay and CFRP was one of the major causes of slabs losing their ultimate load-carrying capacity. Transferring interfacial stress from the slab to FRP or vice versa via a specialized adhesive layer is more efficient than from FRP to overlay members due to the characteristics of the retrofit system, leading to the premature debonding failure of the FRP/overlay interface [27][28][29]. A sound understanding of the behavior of FRP-concrete interface needs to be developed, with a particular interest in FRP/overlay interface for safe in practical design work. In response to this need, the effect of bond-slip laws on the debonding failure of flexural strengthened RC slabs with a hybrid FRP retrofit system is investigated. To the best of our knowledge, previous studies were limited in evaluating the behavior of the externally bonded FRP on a concrete surface without considering cases of CFRP located between two concrete layers of the strengthened slab with a retrofit system [30][31][32][33][34][35][36][37]. There was still considerable uncertainty and difficulty in determining the behavior of CFRP-toconcrete joints due to their intrinsic bond. Insight into an understanding of the bond-slip law between concrete and CFRP using a cohesive approach is a suitable solution for predicting strengthened structural responses [36,38]. Shear stress at the FRP-concrete interface of cohesive elements, in terms of magnitude and distribution, could be a believable explanation for debonding failure [39,40]. It can be difficult to precisely predict the variations in stresses in concrete structures, whereas the load-deflection curve is far less sensitive to crack locations and sizes. The finite element method (FEM) can reasonably predict the interfacial stresses and the delamination load corresponding to deflection more economically than laboratory tests [41]. A finite element (FE) model of the FRP flexuralstrengthened slab can provide a comprehensive understanding of the various bond-slip parameters associated with the CFRP/concrete interface and be applicable in developing robust predictive equations for practice designs [42]. By combining FEM through ABAQUS software and the experimental program, the iterative adjustment of the bondslip parameters can be easily carried out to simulate the experimental load-displacement curve and failure modes, resulting in reduced experiment time and cost [43].  In response to this need, the effect of bond-slip laws on the debonding failure of flexural strengthened RC slabs with a hybrid FRP retrofit system is investigated. To the best of our knowledge, previous studies were limited in evaluating the behavior of the externally bonded FRP on a concrete surface without considering cases of CFRP located between two concrete layers of the strengthened slab with a retrofit system [30][31][32][33][34][35][36][37]. There was still considerable uncertainty and difficulty in determining the behavior of CFRP-to-concrete joints due to their intrinsic bond. Insight into an understanding of the bond-slip law between concrete and CFRP using a cohesive approach is a suitable solution for predicting strengthened structural responses [36,38]. Shear stress at the FRP-concrete interface of cohesive elements, in terms of magnitude and distribution, could be a believable explanation for debonding failure [39,40]. It can be difficult to precisely predict the variations in stresses in concrete structures, whereas the load-deflection curve is far less sensitive to crack locations and sizes. The finite element method (FEM) can reasonably predict the interfacial stresses and the delamination load corresponding to deflection more economically than laboratory tests [41]. A finite element (FE) model of the FRP flexuralstrengthened slab can provide a comprehensive understanding of the various bond-slip parameters associated with the CFRP/concrete interface and be applicable in developing robust predictive equations for practice designs [42]. By combining FEM through ABAQUS software and the experimental program, the iterative adjustment of the bond-slip parameters can be easily carried out to simulate the experimental load-displacement curve and failure modes, resulting in reduced experiment time and cost [43].
In this study, the behavior of a strengthened RC slab with a hybrid FRP retrofit system is described with particular attention to debonding failure. The effectiveness of the hybrid retrofit system in improving the flexural carrying capacity, stiffness, and toughness is investigated and compared with the full-scale experimental results. A three-dimensional FE model is developed with various bond-slip models for evaluating their accuracy by comparing numerical predictions with experimental measurements. Numerical models are also analyzed for sensitivity to the bond-slip parameters to assess their effects on the debonding failure. Following the validation of the model, the various parameters that strongly influence the behavior of retrofitted slabs are identified and discussed.

Retrofitting Mechanism
In some cases, strengthening the underside of the reinforced concrete slabs may not be possible due to logistical constraints and hindrances by other utilities. A hybrid retrofit system combining the tensile strength of FRP and compressive strength of the concrete overlay was applied on top of the existing RC slabs to improve their strength and ductility, as shown in Figure 1. The overlay thickness of the retrofit system has a critical role in pulling the neutral axis toward the overlay zone, and FRP holds tension at failure. For this case, the retrofitting mechanism for the hybrid FRP system was estimated as recommended by ACI 440.2R [44], as shown in Figure 2. Assume that the ultimate concrete strain (ε cu ) is 0.003, and steel yields at yield stress of f y . Without considering the tension force of concrete, the force equilibrium is calculated as follows: In this study, the behavior of a strengthened RC slab with a hybrid FRP retrofit system is described with particular attention to debonding failure. The effectiveness of the hybrid retrofit system in improving the flexural carrying capacity, stiffness, and toughness is investigated and compared with the full-scale experimental results. A three-dimensional FE model is developed with various bond-slip models for evaluating their accuracy by comparing numerical predictions with experimental measurements. Numerical models are also analyzed for sensitivity to the bond-slip parameters to assess their effects on the debonding failure. Following the validation of the model, the various parameters that strongly influence the behavior of retrofitted slabs are identified and discussed.

Retrofitting Mechanism
In some cases, strengthening the underside of the reinforced concrete slabs may not be possible due to logistical constraints and hindrances by other utilities. A hybrid retrofit system combining the tensile strength of FRP and compressive strength of the concrete overlay was applied on top of the existing RC slabs to improve their strength and ductility, as shown in Figure 1. The overlay thickness of the retrofit system has a critical role in pulling the neutral axis toward the overlay zone, and FRP holds tension at failure. For this case, the retrofitting mechanism for the hybrid FRP system was estimated as recommended by ACI 440.2R [44], as shown in Figure 2. Assume that the ultimate concrete strain (εcu) is 0.003, and steel yields at yield stress of fy. Without considering the tension force of concrete, the force equilibrium is calculated as follows: These internal forces are computed as follows: An evaluation of strain compatibility conditions can determine FRP strain at strength limit (εF) as follows: The moment capacity equilibrium condition can be defined as follows:

Experimental Program
The reference slab had a 2440 mm total length, 2290 mm clear span, and 130 × 900 mm 2 cross-section. The tensile reinforcement was spaced at 185 mm spacing with five These internal forces are computed as follows: An evaluation of strain compatibility conditions can determine FRP strain at strength limit (ε F ) as follows: The moment capacity equilibrium condition can be defined as follows:

Experimental Program
The reference slab had a 2440 mm total length, 2290 mm clear span, and 130 × 900 mm 2 cross-section. The tensile reinforcement was spaced at 185 mm spacing with five No.13 bars (φ12.7 mm) at the bottom of the RC slab. Transverse steel consisted of six No.10 bars (φ9.5 mm) spaced at 305 mm center on the center, as shown in Figure 3. The yield stress of the reinforcement used for the RC slab was 400 MPa. The concrete employed for the RC slab had a compressive strength of 27 MPa after 28 days. Table 1 provides the mechanical properties and dimensions of the reference slab.
No.13 bars (ϕ12.7 mm) at the bottom of the RC slab. Transverse steel consisted of six No.10 bars (ϕ9.5 mm) spaced at 305 mm center on the center, as shown in Figure 3. The yield stress of the reinforcement used for the RC slab was 400 MPa. The concrete employed for the RC slab had a compressive strength of 27 MPa after 28 days. Table 1 provides the mechanical properties and dimensions of the reference slab.

h (mm) b (mm) As (mm 2 ) d (mm) f′c (MPa) γc (kg/m 3 ) fy (MPa) Es (GPa)
130 900 632.5 80 27 2400 400 200 A hybrid retrofit system was constructed utilizing a combination of CFRP wet layup laminate with a concrete overlay, as shown in Figure 1. The compressive strength of the concrete overlay was 50 MPa at 28-day age. CFRP had a strength of 600 MPa and an elastic modulus of 40 GPa. The mechanical properties and dimensions used in the FRP retrofit system are given in Table 2. Two full-scale RC slabs with and without a hybrid retrofit system were tested in a four-point bending setup with two concentrated loads, as shown in Figure 4. Before measuring the deflection data, the LVDTs were reset to remove self-weight deflection. The load cell capacity applied in the tests was 5000 kN. The applied load and deflection at the midspan section data were recorded during the experimental procedure.  A hybrid retrofit system was constructed utilizing a combination of CFRP wet layup laminate with a concrete overlay, as shown in Figure 1. The compressive strength of the concrete overlay was 50 MPa at 28-day age. CFRP had a strength of 600 MPa and an elastic modulus of 40 GPa. The mechanical properties and dimensions used in the FRP retrofit system are given in Table 2. Two full-scale RC slabs with and without a hybrid retrofit system were tested in a four-point bending setup with two concentrated loads, as shown in Figure 4. Before measuring the deflection data, the LVDTs were reset to remove self-weight deflection. The load cell capacity applied in the tests was 5000 kN. The applied load and deflection at the mid-span section data were recorded during the experimental procedure.

Theoretical Analysis
The experimental data were checked with the predicted analytical moment capacity and deflection at the mid-span due to the applied load of the reference and strengthened RC slabs, as recommended by ACI 440.2R [44] and ACI 318M [45]. Based on the structure,

Theoretical Analysis
The experimental data were checked with the predicted analytical moment capacity and deflection at the mid-span due to the applied load of the reference and strengthened RC slabs, as recommended by ACI 440.2R [44] and ACI 318M [45]. Based on the structure, a self-weighted equivalent distributed load is computed as follows: The distributed load due to the self-weight of the reference slab and strengthened slab is determined as follows: W c = 2.81 N/mm and W s = 3.46 N/mm. The prediction of the load-carrying capacity corresponding to deflection at mid-span due to the applied load for the reference RC slab is given in Table 3. Table 3. Analysis of the reference RC slab.

Calculation Reference Slab
Compressive concrete block depth a = The RC slab and CFRP laminate were assumed to be in perfect bond until the load-carrying capacity was reached. The prediction of the load-carrying capacity corresponding to deflection at mid-span for the CFRP-retrofitted slab due to the applied load is shown in Table 4. Table 4. Analysis of the retrofitted RC slab.

Calculation Retrofitted Slab
Compressive concrete block depth

Finite Element Mesh
The simulations were conducted using one-quarter of the control and strengthened slab specimens, based on specimen symmetry. The symmetric pane was simulated for the xand z-axis by restricting translation in directions 1 and 3. The concrete slab and concrete overlay with reduced integration were constructed with eight-node solid elements (C3D8R). Two-node linear truss elements (T3D2) and four-node shell elements (S4R) were used to model the reinforcements and CFRP laminate, respectively. A mesh size of 20 mm was suggested based on preliminary mesh refinement studies, ranging from 10 mm to 25 mm, for balancing accuracy and computational cost. The embedded function in ABAQUS/CAE 2022 was implemented to simulate concrete-steel bonding. The models were studied using static analysis in ABAQUS/standard. The models under monotonic displacement loads surveyed the behavior of a slab subjected to a four-point bending test. The steel bar was discretized as a rigid part to prevent manufactured stress concentrations under the loading points, as shown in Figure 5.

Concrete
A model of concrete plastic damage was used with a failure mechanism of compressive crushing and tensile cracking. The stress-strain curve for concrete under uniaxial compression and tension was calculated by using the CEB-FIP model [46]. Compressive strength was used to estimate the tensile strength (fct) and elastic modulus (Ec) of concrete, as recommended by ACI 318M [45].
Fracture energy (Gcr) is an inelastic parameter associated with the softening part of the curve. In mode I, fracture energy is defined as the area under the softening curve.

Concrete
A model of concrete plastic damage was used with a failure mechanism of compressive crushing and tensile cracking. The stress-strain curve for concrete under uniaxial compression and tension was calculated by using the CEB-FIP model [46]. Compressive strength was used to estimate the tensile strength (f ct ) and elastic modulus (E c ) of concrete, as recommended by ACI 318M [45].
Fracture energy (G cr ) is an inelastic parameter associated with the softening part of the curve. In mode I, fracture energy is defined as the area under the softening curve. Figure 6 illustrates the concrete softening curve under uniaxial tension. A linear relationship between the tension damage and crack opening (δ) is assumed to identify tensile damage. The tension damage variable ranges from zero (undamaged material) to one (total loss of strength). Under uniaxial compression, a linear response is observed until reaching the initial value (σc0). In the plastic regime, the response is generally characterized by stress hardening, followed by strain softening beyond the ultimate compressive stress [47].

Reinforced Steel
It is assumed that the reinforcement transmits force axially, and the most common perfectly linear-elastic model for reinforcement was used. Steel's Poisson's ratio is 0.3, and the other mechanical properties are listed in Table 2.

CFRP
For strengthening, CFRP laminate is supposed to be linear and exert isotropic behavior. CFRP failure criteria were not considered for this study because tensile stress in CFRP is too low, compared with its ultimate strength, and therefore CFRP laminate hardly ruptures for these types of structures. The material properties of CFRP, including its thickness, elastic modulus, and ultimate strength, as shown in Table 2, are provided by manufacturers.

FRP-to-Concrete Interface Model
The FRP-to-concrete interface was modeled using cohesive zones. In this study, a bilinear traction-separation model was employed to depict the bonding characteristics of linear adhesive at the interface, as shown in Figure 7 [49]. Under uniaxial compression, a linear response is observed until reaching the initial value (σ c0 ). In the plastic regime, the response is generally characterized by stress hardening, followed by strain softening beyond the ultimate compressive stress [47].

Reinforced Steel
It is assumed that the reinforcement transmits force axially, and the most common perfectly linear-elastic model for reinforcement was used. Steel's Poisson's ratio is 0.3, and the other mechanical properties are listed in Table 2.

CFRP
For strengthening, CFRP laminate is supposed to be linear and exert isotropic behavior. CFRP failure criteria were not considered for this study because tensile stress in CFRP is too low, compared with its ultimate strength, and therefore CFRP laminate hardly ruptures for these types of structures. The material properties of CFRP, including its thickness, elastic modulus, and ultimate strength, as shown in Table 2, are provided by manufacturers.

FRP-to-Concrete Interface Model
The FRP-to-concrete interface was modeled using cohesive zones. In this study, a bilinear traction-separation model was employed to depict the bonding characteristics of linear adhesive at the interface, as shown in Figure 7 [49]. facturers.

FRP-to-Concrete Interface Model
The FRP-to-concrete interface was modeled using cohesive zones. In this study, a bilinear traction-separation model was employed to depict the bonding characteristics of linear adhesive at the interface, as shown in Figure 7   The initiation and evolution of damage in the available traction-separation model were initially assumed to be a linear-elastic behavior [50,51], which can be written in a matrix form as follows: Figure 7. Bilinear traction-separation response.
The initiation and evolution of damage in the available traction-separation model were initially assumed to be a linear-elastic behavior [50,51], which can be written in a matrix form as follows: For specifying the damage initiation, quadratic nominal stress and maximum nominal stress between the CFRP and concrete interface were employed. The damage criterion for quadratic stress may be depicted as follows: The damage criterion for maximum stress can be described as follows: For specifying the damage evolution, the influence of power law and the Benzeggagh-Kenane (BK) fracture criteria on the structure behavior was further evaluated [52]. The power law criterion may be given by: The BK fracture criterion may be defined as follows: In this study, the quantities G f n , G f s , and G f t were assumed to be equal. For FE modeling, first, a perfect FRP-to-concrete bond was assumed to reflect the response of the flexural strengthened RC slab using the ABAQUS program. Then, the interaction behavior of the FRP-to-concrete bond in the strengthened slab applied to the retrofit system was evaluated using the existing constitutive models, as shown in Table 5 [31][32][33][34][35][36]. Table 5. Existing bond-slip models.

Analysis Procedure and Flowchart
A three-dimensional FE model using ABAQUS/CAE 2022 software was performed to simulate the behavior of retrofitted slabs. In the initial phase, an accurate FE model of the reference slab was assessed by comparison with the experimental results. In the second phase, a retrofitted slab model was developed based on the reference slab model, with the addition of FRP laminate and concrete overlay. The well-known bond-slip models were used to simulate the FRP-to-concrete interface to accurately simulate the behavior of the retrofitted slab validated by experiment results. If the available bondslip models were not appropriate, the bond-slip parameters were revised to show good agreement with the experimental data. A flowchart of the analysis procedure for the FE model is shown in Figure 8.

Experimental Analysis
The experimental data from the reference slab showed a good agreement with those of the theoretical analysis. The RC slab behaved linearly and elastically until the first crack appeared at 17 kN of the applied load. The observed result revealed a higher load-carrying capacity of 43.5 kN than the theoretical prediction of 39.9 kN. The ultimate loading (50.8 kN) was 27.3% greater than the predicted load. The linearly elastic-perfectly plastic behavior of the steel can be the reason for the lower strength of the prediction because strain-hardened steel shows higher tensile strength in practice. The mid-span deflection reached a maximum value of 72 mm, as shown in Figure 9. The slope of the load-deflection curve within the elastic limits depicts the stiffness (K) of the slab. The stiffness of the RC slab was determined as 10.1 kN/mm. Calculating the toughness of a structure can be accomplished by integrating the load-deflection curve. Accordingly, the reference slab had a toughness of 3129.2 kNmm.
For the retrofit slab, the cracks emerged and expanded most strongly at the mid-span zone-between the loading points, as shown in Figure 10. The retrofitted slab showed a load-carrying capacity of 69.2 kN at 14 mm mid-span deflection. This was 9% higher than the predicted load of 63 kN. The structure lost carrying capacity due to the micro-buckling of CFRP laminate at the applied load of 32 kN. Then, the applied load continued to increase until CFRP laminate was delaminated at 76.9 kN, corresponding to 22 mm of the mid-span deflection. After delamination, the global response of the concrete slab was equivalent to that of the reference concrete slab. Its failure load and mid-span deflection were 49.2 kN and 81 mm, respectively, as shown in Figure 9. It is possible to define the stiffness of the retrofitted slab at 21.3 kN/mm, which was 2.11 times that of the reference slab. The retrofitted slab had a toughness of 4034.1 kNmm, 1.29 times greater than the reference slab. Compared with theoretical analysis, Table 7 summarizes the experimental results for load-carrying capacity corresponding to the mid-span deflection of the slabs.

Experimental Analysis
The experimental data from the reference slab showed a good agreement with those of the theoretical analysis. The RC slab behaved linearly and elastically until the first crack appeared at 17 kN of the applied load. The observed result revealed a higher loadcarrying capacity of 43.5 kN than the theoretical prediction of 39.9 kN. The ultimate loading (50.8 kN) was 27.3% greater than the predicted load. The linearly elastic-perfectly plastic behavior of the steel can be the reason for the lower strength of the prediction because strain-hardened steel shows higher tensile strength in practice. The mid-span deflection reached a maximum value of 72 mm, as shown in Figure 9. The slope of the load-deflection curve within the elastic limits depicts the stiffness (K) of the slab. The stiffness of the RC slab was determined as 10.1 kN/mm. Calculating the toughness of a structure can be accomplished by integrating the load-deflection curve. Accordingly, the reference slab had a toughness of 3129.2 kNmm.
For the retrofit slab, the cracks emerged and expanded most strongly at the mid-span zone-between the loading points, as shown in Figure 10. The retrofitted slab showed a load-carrying capacity of 69.2 kN at 14 mm mid-span deflection. This was 9% higher than the predicted load of 63 kN. The structure lost carrying capacity due to the micro-buckling of CFRP laminate at the applied load of 32 kN. Then, the applied load continued to increase until CFRP laminate was delaminated at 76.9 kN, corresponding to 22 mm of the mid-span deflection. After delamination, the global response of the concrete slab was equivalent to that of the reference concrete slab. Its failure load and mid-span deflection were 49.2 kN and 81 mm, respectively, as shown in Figure 9. It is possible to define the stiffness of the retrofitted slab at 21.3 kN/mm, which was 2.11 times that of the reference slab. The retrofitted slab had a toughness of 4034.1 kNmm, 1.29 times greater than the reference slab. Compared with theoretical analysis, Table 7 summarizes the experimental results for load-carrying capacity corresponding to the mid-span deflection of the slabs.

Numerical Analysis
In this part of the study, the FE model of the retrofitted slab was used to examine the failure mode and delamination load after obtaining a good agreement between the FE model of the reference slab and the experimental results using the commercial ABAQUS program. A perfect bond model and well-known bond-slip models were applied to simulate the slab behavior through the load-deflection relationship. The material properties of the experimental structure were used to calculate the input parameters for the bondslip models, as shown in Table 6. The FE models were employed to evaluate the suitability of bond-slip laws based on the observed results. Most existing bond-slip models can predict the load-carrying capacity of the retrofitted slab when errors in the prediction of FE models are less than 7%. Nevertheless, the bond-slip models could not accurately predict the mid-span deflection corresponding to the load-carrying capacity and stiffness, with errors higher than 17% and 32%, respectively. The behavior of the well-known bond-slip models also showed a perfect bond between CFRP and concrete, even though delamination failure was observed in the experiment, as shown in Figure 11. These discrepancies

Numerical Analysis
In this part of the study, the FE model of the retrofitted slab was used to examine the failure mode and delamination load after obtaining a good agreement between the FE model of the reference slab and the experimental results using the commercial ABAQUS program. A perfect bond model and well-known bond-slip models were applied to simulate the slab behavior through the load-deflection relationship. The material properties of the experimental structure were used to calculate the input parameters for the bondslip models, as shown in Table 6. The FE models were employed to evaluate the suitability of bond-slip laws based on the observed results. Most existing bond-slip models can predict the load-carrying capacity of the retrofitted slab when errors in the prediction of FE models are less than 7%. Nevertheless, the bond-slip models could not accurately predict the mid-span deflection corresponding to the load-carrying capacity and stiffness, with errors higher than 17% and 32%, respectively. The behavior of the well-known bond-slip models also showed a perfect bond between CFRP and concrete, even though delamination failure was observed in the experiment, as shown in Figure 11. These discrepancies

Numerical Analysis
In this part of the study, the FE model of the retrofitted slab was used to examine the failure mode and delamination load after obtaining a good agreement between the FE model of the reference slab and the experimental results using the commercial ABAQUS program. A perfect bond model and well-known bond-slip models were applied to simulate the slab behavior through the load-deflection relationship. The material properties of the experimental structure were used to calculate the input parameters for the bond-slip models, as shown in Table 6. The FE models were employed to evaluate the suitability of bond-slip laws based on the observed results. Most existing bond-slip models can predict the loadcarrying capacity of the retrofitted slab when errors in the prediction of FE models are less than 7%. Nevertheless, the bond-slip models could not accurately predict the mid-span deflection corresponding to the load-carrying capacity and stiffness, with errors higher than 17% and 32%, respectively. The behavior of the well-known bond-slip models also showed a perfect bond between CFRP and concrete, even though delamination failure was observed in the experiment, as shown in Figure 11. These discrepancies may have resulted from the improper estimations of the bond-slip parameters involving fracture energy, the initial stiffness of the interface elements, maximum interfacial stress, damage initiation criteria, and mixed-mode failure criteria. The FEM results of the load-carrying capacity corresponding to the mid-span deflection and stiffness within the elastic limits of the strengthened slab, compared with its experimental counterpart, are summarized in Table 8. may have resulted from the improper estimations of the bond-slip parameters involving fracture energy, the initial stiffness of the interface elements, maximum interfacial stress, damage initiation criteria, and mixed-mode failure criteria. The FEM results of the loadcarrying capacity corresponding to the mid-span deflection and stiffness within the elastic limits of the strengthened slab, compared with its experimental counterpart, are summarized in Table 8.

Bond-Slip Analysis
The existing bond-slip models tended to overestimate the bond strength of the CFRPto-concrete interface. In this part of the study, a sensitivity analysis of the bond-slip input parameters was performed to compare its results with the structural behavior obtained from the experimental data. The bond-slip model of Neubauer and Rostasy was used to illustrate the numerical results from the strengthened slabs. Initially, the CFRP and concrete bonding interface was assumed to fail using quadratic traction. A cohesive property coefficient of η = 1 was applied to depict the influence of the opening and sliding failure modes on fracture energy. The authors investigated the sensitivity of the results in the interfacial fracture energy between CFRP and concrete. For this case, a change was made from 0.01Gf to Gf for the interfacial fracture energy. There was no evidence that Gf had a notable effect on the numerical results, as shown in Figure 12.

Bond-Slip Analysis
The existing bond-slip models tended to overestimate the bond strength of the CFRPto-concrete interface. In this part of the study, a sensitivity analysis of the bond-slip input parameters was performed to compare its results with the structural behavior obtained from the experimental data. The bond-slip model of Neubauer and Rostasy was used to illustrate the numerical results from the strengthened slabs. Initially, the CFRP and concrete bonding interface was assumed to fail using quadratic traction. A cohesive property coefficient of η = 1 was applied to depict the influence of the opening and sliding failure modes on fracture energy. The authors investigated the sensitivity of the results in the interfacial fracture energy between CFRP and concrete. For this case, a change was made from 0.01G f to G f for the interfacial fracture energy. There was no evidence that G f had a notable effect on the numerical results, as shown in Figure 12. The initial stiffness of the interface elements was changed from 0.01K0 to K0 to evaluate its effect. It had a substantial influence on the load-deflection curve's slope of the structure. Based on the simulated results, the stiffness of the retrofitted slab decreased with decreasing the initial stiffness but was not proportional. With the initial stiffness of 0.1K0, the global response of the slab within the elastic limits could match the experimental results, as shown in Figure 13. However, the constitutive model exhibited a perfect bond of CFRP-to-concrete interface. As mentioned above, the low bond strength between the concrete overlay and CFRP was the main reason for the delamination failure of the retrofit slab. It was possible to obtain better explanations for the differences between the simulated and experimental results by modifying the maximum shear stress. Therefore, the maximum interfacial stress of the concrete overlay/CFRP was further investigated using sensitivity analysis. An appropriate value for this parameter was crucial for the accurate prediction of CFRP debonding failures. As shown in Figure 14a, the maximum shear stress was increased from 0.1τmax to τmax. The shear stress values had a considerable impact on the debonding failure. The The initial stiffness of the interface elements was changed from 0.01K 0 to K 0 to evaluate its effect. It had a substantial influence on the load-deflection curve's slope of the structure. Based on the simulated results, the stiffness of the retrofitted slab decreased with decreasing the initial stiffness but was not proportional. With the initial stiffness of 0.1K 0 , the global response of the slab within the elastic limits could match the experimental results, as shown in Figure 13. However, the constitutive model exhibited a perfect bond of CFRP-to-concrete interface. The initial stiffness of the interface elements was changed from 0.01K0 to K0 to evaluate its effect. It had a substantial influence on the load-deflection curve's slope of the structure. Based on the simulated results, the stiffness of the retrofitted slab decreased with decreasing the initial stiffness but was not proportional. With the initial stiffness of 0.1K0, the global response of the slab within the elastic limits could match the experimental results, as shown in Figure 13. However, the constitutive model exhibited a perfect bond of CFRP-to-concrete interface. As mentioned above, the low bond strength between the concrete overlay and CFRP was the main reason for the delamination failure of the retrofit slab. It was possible to obtain better explanations for the differences between the simulated and experimental results by modifying the maximum shear stress. Therefore, the maximum interfacial stress of the concrete overlay/CFRP was further investigated using sensitivity analysis. An appropriate value for this parameter was crucial for the accurate prediction of CFRP debonding failures. As shown in Figure 14a, the maximum shear stress was increased from 0.1τmax to τmax. The shear stress values had a considerable impact on the debonding failure. The Figure 13. Sensitivity to initial stiffness of interface elements.
As mentioned above, the low bond strength between the concrete overlay and CFRP was the main reason for the delamination failure of the retrofit slab. It was possible to obtain better explanations for the differences between the simulated and experimental results by modifying the maximum shear stress. Therefore, the maximum interfacial stress of the concrete overlay/CFRP was further investigated using sensitivity analysis. An appropriate value for this parameter was crucial for the accurate prediction of CFRP debonding failures. As shown in Figure 14a, the maximum shear stress was increased from 0.1τ max to τ max . The shear stress values had a considerable impact on the debonding failure. The delamination load of the retrofitted slab was sensitive to changes in the ranges of maximum shear stress from 0.1τ max to 0.15τ max , whereas a perfect bond was shown for values greater than or equal to 0.2τ max . It could also be easily seen that Neubauer's and Rostasy's bond-slip law, along with other equivalent models, overestimated the ultimate shear strength of the concrete overlay-CFRP interface. delamination load of the retrofitted slab was sensitive to changes in the ranges of maximum shear stress from 0.1τmax to 0.15τmax, whereas a perfect bond was shown for values greater than or equal to 0.2τmax. It could also be easily seen that Neubauer's and Rostasy's bond-slip law, along with other equivalent models, overestimated the ultimate shear strength of the concrete overlay-CFRP interface.
According to the simulations, the delamination loads were 66.2 kN and 74 kN, corresponding to 10.1 mm and 48 mm of the mid-span deflection, respectively, occurring with the maximum shear stresses of 0.1τmax and 0.15τmax, respectively. The numerical results also indicated that the debonding failure emerged earlier, at 0.1τmax and later at 0.15τmax, compared with the experimental results. Moreover, the overestimation of the initial stiffness should also be highlighted in this case. The maximum shear stress ranged from 0.1τmax to 0.15τmax, combined with the initial stiffness of 0.1K0; thus, the case was further evaluated. The delamination load and mid-span deflection of the retrofitted slab increased with the increase in the maximum shear stress but was not proportional; however, there were no notable differences in the delamination load after the yield in steel. The global response of the strengthened slab was displayed as a function of the shear stress values. As a result, it was possible to reasonably simulate the behavior of the retrofitted slab and compare it with its experimental counterpart with the initial stiffness of 0.1K0 and maximum shear stress of 0.13τmax, as shown in Figure 14b. In this part of the study, we analyzed the damage initiation criteria as the quadratic and maximum nominal stresses based on the proposed results mentioned above. The quadratic stress criterion is defined by accounting for the effect of traction components, while the maximum stress criterion is established by comparing the traction components to their respective allowable values [54]. Despite having similar responses until the yield in steel, the quadratic and maximum nominal stress criteria exhibited different mid-span deflections at debonding failures of 22 mm and 28 mm, respectively, as shown in Figure 15. Poor results were obtained using the maximum stress criterion for predicting delamination in the retrofitted slabs. There was a significant interaction between stresses for decohesion elements dealing with mixed-mode delamination onset and propagation, as mentioned by Cui et al. [55]. Debonding behavior can emerge before any of the involved traction components reach their respective allowable limits. An accurate prediction of delamination failure requires taking into account the interaction. Thus, the damage initiation criterion applying the quadratic nominal stress was suggested in this case. According to the simulations, the delamination loads were 66.2 kN and 74 kN, corresponding to 10.1 mm and 48 mm of the mid-span deflection, respectively, occurring with the maximum shear stresses of 0.1τ max and 0.15τ max , respectively. The numerical results also indicated that the debonding failure emerged earlier, at 0.1τ max and later at 0.15τ max , compared with the experimental results. Moreover, the overestimation of the initial stiffness should also be highlighted in this case. The maximum shear stress ranged from 0.1τ max to 0.15τ max , combined with the initial stiffness of 0.1K 0 ; thus, the case was further evaluated. The delamination load and mid-span deflection of the retrofitted slab increased with the increase in the maximum shear stress but was not proportional; however, there were no notable differences in the delamination load after the yield in steel. The global response of the strengthened slab was displayed as a function of the shear stress values. As a result, it was possible to reasonably simulate the behavior of the retrofitted slab and compare it with its experimental counterpart with the initial stiffness of 0.1K 0 and maximum shear stress of 0.13τ max , as shown in Figure 14b. In this part of the study, we analyzed the damage initiation criteria as the quadratic and maximum nominal stresses based on the proposed results mentioned above. The quadratic stress criterion is defined by accounting for the effect of traction components, while the maximum stress criterion is established by comparing the traction components to their respective allowable values [54]. Despite having similar responses until the yield in steel, the quadratic and maximum nominal stress criteria exhibited different mid-span deflections at debonding failures of 22 mm and 28 mm, respectively, as shown in Figure 15. Poor results were obtained using the maximum stress criterion for predicting delamination in the retrofitted slabs. There was a significant interaction between stresses for decohesion elements dealing with mixed-mode delamination onset and propagation, as mentioned by Cui et al. [55]. Debonding behavior can emerge before any of the involved traction components reach their respective allowable limits. An accurate prediction of delamination failure requires taking into account the interaction. Thus, the damage initiation criterion applying the quadratic nominal stress was suggested in this case. Here, we describe the effects of damage evolution based on energy. The load-deflection relationship was derived from the power law or BK law, as shown in Figure 16a. Likewise, a study of the cohesive coefficient sensitivity of FE results was also conducted, the result of which is shown in Figure 16b. Nevertheless, these criteria had no effect on the delamination load or global response of the strengthened RC slabs. In brief, these results were consistent with the previous literature regarding the influential factors on debonding failure, which is significantly affected by the maximum interfacial stress [40,56,57] but insensitive to changes in fracture energy and mixed-mode failure criteria [53]. However, disagreement can also be observed once the global response of the strengthened slab with a hybrid retrofit system was found to be impacted by changing the criteria of the stiffness of interface elements and damage initiation. Developing an appropriate bond-slip model will require further investigation to determine the precise values for notable influential factors on debonding failure.

Conclusions
In this paper, a hybrid retrofit system considering the difficulty of accessing and installing CFRP laminates to the underside of RC slabs to enhance the flexural carrying ca- Here, we describe the effects of damage evolution based on energy. The load-deflection relationship was derived from the power law or BK law, as shown in Figure 16a. Likewise, a study of the cohesive coefficient sensitivity of FE results was also conducted, the result of which is shown in Figure 16b. Nevertheless, these criteria had no effect on the delamination load or global response of the strengthened RC slabs. Here, we describe the effects of damage evolution based on energy. The load-deflection relationship was derived from the power law or BK law, as shown in Figure 16a. Likewise, a study of the cohesive coefficient sensitivity of FE results was also conducted, the result of which is shown in Figure 16b. Nevertheless, these criteria had no effect on the delamination load or global response of the strengthened RC slabs. In brief, these results were consistent with the previous literature regarding the influential factors on debonding failure, which is significantly affected by the maximum interfacial stress [40,56,57] but insensitive to changes in fracture energy and mixed-mode failure criteria [53]. However, disagreement can also be observed once the global response of the strengthened slab with a hybrid retrofit system was found to be impacted by changing the criteria of the stiffness of interface elements and damage initiation. Developing an appropriate bond-slip model will require further investigation to determine the precise values for notable influential factors on debonding failure.

Conclusions
In this paper, a hybrid retrofit system considering the difficulty of accessing and installing CFRP laminates to the underside of RC slabs to enhance the flexural carrying ca- In brief, these results were consistent with the previous literature regarding the influential factors on debonding failure, which is significantly affected by the maximum interfacial stress [40,56,57] but insensitive to changes in fracture energy and mixed-mode failure criteria [53]. However, disagreement can also be observed once the global response of the strengthened slab with a hybrid retrofit system was found to be impacted by changing the criteria of the stiffness of interface elements and damage initiation. Developing an appropriate bond-slip model will require further investigation to determine the precise values for notable influential factors on debonding failure.

Conclusions
In this paper, a hybrid retrofit system considering the difficulty of accessing and installing CFRP laminates to the underside of RC slabs to enhance the flexural carrying capacity, stiffness, and toughness of the existing RC slabs was proposed. The proposed retrofitting mechanism proved suitable through a good agreement between the experimental results and theoretical analyses at the mid-span section. The global response and the effect of the bond-slip law in predicting the debonding failure of the retrofitted RC slabs were presented. The FE model using a sensitivity analysis based on the load-deflection relationship was performed to evaluate the parameters influencing the interfacial behavior between CFRP and concrete. Based on the obtained results, the following conclusions could be drawn as follows: The efficiency of the retrofit system was verified through the experimental test of the strengthened slab in the positive moment section. The retrofit system enhanced the load-carrying capacity of the slab by 59%, stiffness by 111%, and toughness by 29%.
The neutral axis of the retrofitted slab was located within the overlay, and CFRP held tension at the mid-span section at the ultimate failure state.
The bond-slip models overestimated the criteria regarding the damage initiation between CFRP and the concrete overlay. The numerical analysis results with the initial stiffness of 0.1K 0 and maximum shear stress of 0.13τ max were compared with the corresponding values of Neubauer's and Rostasy's bond-slip law and showed a good agreement with the experimental data.
The stiffness of structures is notably impacted by the initial stiffness of the interface elements, while the delamination load and failure load can be decided by the maximum interfacial stress. Quadratic nominal stress is recommended as a criterion for damage initiation based on the considerable interaction between stresses.
The global response of the retrofitted slab was not sensitive to changes in the interfacial fracture energy and damage evolution regarding the mix-mode failure criteria and the cohesive coefficient of the interface.
Our study contributes to the evaluation of the impact of bond-slip parameters on the behavior of retrofitted slabs. Nonetheless, further research into an appropriate bond-slip model for retrofitted RC structures should consider the influence of factors regarding materials properties and geometric dimensions.

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