The Bond-Slip Relationship at FRP-to-Brick Interfaces under Dynamic Loading

Interface debonding between fiber reinforced polymers (FPR) and substrates is the principal failure mode for FRP-reinforced structure. To understand the bond–slip relationship at FRP-to-brick interfaces under dynamic loading, the influences of the dynamic enhancement of material performance on the bond–slip curve were studied. Single-lap shear tests under two different loading rates were performed, and the slip distribution curves at different loading stages were fitted to derive the bond–slip relationship. Then a numerical model considering the strain rate effects on materials was built and verified with test results. Further, the influences of brick strength, FRP stiffness and slip rate on the bond–slip relationship were investigated numerically. The research results show that FRP stiffness mainly influences the shape of the bond–slip curve, while brick strength mainly influences the amplitude of the bond–slip curve. The variations of the bond–slip relationship under dynamic loading, i.e., under different slip rates, are mainly caused by the dynamic enhancement of brick strength, and also by the dynamic enhancement of FRP stiffness, especially within a specific slip rate range. The proposed empirical formula considering dynamic FRP stiffness and dynamic brick strength can be used to predict the bond–slip relationship at the FRP-to-brick interface under dynamic loading.


Introduction
Externally bonded fiber reinforced polymer (FRP) is widely used as a proven reinforcement technique for masonry structures [1][2][3][4][5]. Practical applications and numerous studies have found that the debonding of FRP from the surface of substrates (e.g., concrete, clay brick and mortar) is the principal failure mode of reinforced structures [1,[6][7][8][9]. The bondslip relationship between FRP and substrate is the basis for understanding the process of interfacial debonding. The masonry structure is likely to be subjected to dynamic loads such as earthquakes, explosions and impacts. Experimental studies of the FRP-to-substrate interface show that the dynamic interfacial strength is significantly higher than that of static [10][11][12][13], and the dynamic bond-slip relationship is naturally different from that of static [14][15][16][17][18]. An accurate bond-slip relationship can be used for strength calculations of FRP-reinforced structures, helping to design safer and more economical FRP-reinforced solutions. Therefore, the study of the bond-slip relationship at the FRP-to-brick interface under dynamic loading is of great significance.
Almost all the methods for deriving the bond-slip relationship at the FRP-to-substrate interface are based on the pull test [19] (including the single-lap shear test, the double-lap shear test and the beam test) or its numerical simulation. According to the differences in the test data used, the methods for calculating the bond-slip relationship are roughly divided into three: one direct method and two indirect methods.
The direct method is most commonly used. The method directly applies the differential and integral of the axial strain of the FRP to calculate the bond stress and slip, respectively, and then derive the bond-slip relationship by synthesizing the data from different loading stages [20]. Based on the direct method of using strain data, different models have been established to fit the bond-slip relationship. These bond-slip models all conform to the tendency of tensile softening, but with different shapes: (1) bilinear [21,22], e.g., Monti et al.'s model [22]; (2) multilinear [23,24], e.g., Ghiassi et al.'s model [24]; (3) single curve type [25][26][27][28][29][30][31], e.g., Popovics et al.'s model [25] and its development of Nakaba et al.'s model [26] and Savioa et al.'s model [27]; (4) double curve type, e.g., Dai and Ueda's model [32] and Lu et al.'s model [33]; (5) mixed type with line and curve, e.g., Pan and Wu's model [34]. All of these five types of models appear in the study of the bond-slip relationship at the FRP-concrete interface.
Some studies have suggested that there are three main reasons why a consistent, reliable and accurate bond-slip relationship cannot be derived with the direct method. First, the debonding process is difficult to capture due to the highly nonlinear and brittle nature of local fractures [35]. Second, the irregularities in the surface treatment and adhesive layer cause an irregular fluctuation in strain measurements [36]. Third, the random distribution of cracks and aggregates in concrete also induces irregular fluctuations in strain measurements, resulting in a drastic variation in bond stress calculated from strain data [33]. To overcome this issue, Ueda and Dai et al. [35,37,38] ignored the dense strain information, fitted the load-slip curve at the loading end and then derived the bond-slip relationship indirectly based on a simple and rigorous analytical method. To further eliminate the chance errors of the experimental measurements and improve the accuracy of the analysis results, Wu et al. [36,39,40] fitted the slip distribution curves at different loading stages simultaneously to calculate the bond-slip relationship indirectly.
Summarizing the concerns of the above references, it can be found that the research on the bond-slip relationship is mainly focused on the FRP-to-concrete interface, while there is less research on the FRP-to-brick interface. The research on the bond-slip relationship at the FRP-to-brick interface under dynamic loading is even more limited. Besides, due to the harsh requirements of the dynamic loading control on the experimental equipment, most of the studies are on the tensile strength of the interface and lack the study of the bond-slip relationship.
Accordingly, the bond-slip relationship at the FRP-to-brick interface under dynamic loading is investigated in this paper. Firstly, a numerical model of the FRP-to-brick interface was constructed based on the single-lap shear tests under two different loading rates. Secondly, the effects of FRP stiffness and brick strength on the bond-slip relationship were studied numerically. Finally, the dynamic enhancement effect on material performance was integrated to study and analyze the dynamic bond-slip relationship under different slip rates.

Single-Lap Shear Tests
The digital image correlation (DIC) method was adopted to implement full-field strain measurements, and a special fixture was designed to ensure stability during the dynamic loading. The CFRP-to-brick (CFRP, carbon fiber reinforced polymer) single-lap shear tests under two different loading rates were performed on a universal testing machine.

Materials and Specimen Preparation
The bricks used in the tests are standard commercial clay bricks with a size of 240 × 115 × 53 mm 3 . The carbon cloth is a unidirectional woven carbon fiber fabric UT30-30G produced by Toray Japan. The epoxy adhesive HM-180C3P is a two-component fiber adhesive produced by Horse China specifically for structural reinforcement. Their specific material parameters are shown in Table 1. The CFRP-to-brick specimens are shown in Figure 1. The carbon cloth with a width of 40 mm was bonded on the centerline of the brick, which contains a 160 mm long bonded region and a 40 mm long unbonded region. The surface of the bricks was sanded smooth and cleaned previously. Then the unbonded region was covered with tape to prevent impregnation by the adhesive. Finally, the carbon cloth was bonded to the brick by the wet lay-up method.  The CFRP-to-brick specimens are shown in Figure 1. The carbon cloth with a width of 40 mm was bonded on the centerline of the brick, which contains a 160 mm long bonded region and a 40 mm long unbonded region. The surface of the bricks was sanded smooth and cleaned previously. Then the unbonded region was covered with tape to prevent impregnation by the adhesive. Finally, the carbon cloth was bonded to the brick by the wet lay-up method.

Instruments and Test Procedures
With the fixture as shown in Figure 2, the single-lap shear tests were performed on an MTS universal testing machine (MTS, Eden Prairie, MN, America). The loading rate of the MTS universal testing machine with an ultimate load of 100 KN ranges from 1 to 1000 mm/min. Two groups of tests were performed, Group A: 10 mm/min and Group B: 1000 mm/min. Each group contains five specimens.

Instruments and Test Procedures
With the fixture as shown in Figure 2, the single-lap shear tests were performed on an MTS universal testing machine (MTS, Eden Prairie, MN, America). The loading rate of the MTS universal testing machine with an ultimate load of 100 KN ranges from 1 to 1000 mm/min. Two groups of tests were performed, Group A: 10 mm/min and Group B: 1000 mm/min. Each group contains five specimens.  The CFRP-to-brick specimens are shown in Figure 1. The carbon cloth with a width of 40 mm was bonded on the centerline of the brick, which contains a 160 mm long bonded region and a 40 mm long unbonded region. The surface of the bricks was sanded smooth and cleaned previously. Then the unbonded region was covered with tape to prevent impregnation by the adhesive. Finally, the carbon cloth was bonded to the brick by the wet lay-up method.

Instruments and Test Procedures
With the fixture as shown in Figure 2, the single-lap shear tests were performed on an MTS universal testing machine (MTS, Eden Prairie, MN, America). The loading rate of the MTS universal testing machine with an ultimate load of 100 KN ranges from 1 to 1000 mm/min. Two groups of tests were performed, Group A: 10 mm/min and Group B: 1000 mm/min. Each group contains five specimens.   The strain was measured by a two-dimensional digital image correlation (DIC) method combined with a high-speed camera (iX Cameras, Rochford Essex, UK). As shown in Figure 3, the debonding process was recorded by a high-speed camera (ix cameras i-SPEED716) with a maximum frame rate of 500,000 fps. The required light was provided by two 200 W LED lights. The high-speed camera was triggered by a TTL (Transistor-Transistor Logic) pulse from the testing machine. The strain was measured by a two-dimensional digital image correlation (DIC) method combined with a high-speed camera (iX Cameras, Rochford Essex, UK). As shown in Figure 3, the debonding process was recorded by a high-speed camera (ix cameras i-SPEED716) with a maximum frame rate of 500,000 fps. The required light was provided by two 200 W LED lights. The high-speed camera was triggered by a TTL (Transistor-Transistor Logic) pulse from the testing machine.

Test Results and Analysis
The test results revealed that all specimens exhibited the peeling of CFRP and a thin layer of brick, as shown in Figure 4, the fracture of the brick dominated the failure of FRPto-brick interface. The test results are in agreement with the related literature [30,41]. The load-slip curves at the loading end are shown in Figure 5. The loading process can be roughly divided into three stages: linear growth stage, softening stage and stable stage. With the increase of slip, the tensile force first goes through the linear growth stage and the softening stage to reach the ultimate load Fu, then into the stable stage. In the stable stage, the slip further increases, the tensile force remains basically unchanged.

Test Results and Analysis
The test results revealed that all specimens exhibited the peeling of CFRP and a thin layer of brick, as shown in Figure 4, the fracture of the brick dominated the failure of FRP-to-brick interface. The test results are in agreement with the related literature [30,41].
The load-slip curves at the loading end are shown in Figure 5. The loading process can be roughly divided into three stages: linear growth stage, softening stage and stable stage. With the increase of slip, the tensile force first goes through the linear growth stage and the softening stage to reach the ultimate load F u , then into the stable stage. In the stable stage, the slip further increases, the tensile force remains basically unchanged. The strain was measured by a two-dimensional digital image correlation (DIC) method combined with a high-speed camera (iX Cameras, Rochford Essex, UK). As shown in Figure 3, the debonding process was recorded by a high-speed camera (ix cameras i-SPEED716) with a maximum frame rate of 500,000 fps. The required light was provided by two 200 W LED lights. The high-speed camera was triggered by a TTL (Transistor-Transistor Logic) pulse from the testing machine.

Test Results and Analysis
The test results revealed that all specimens exhibited the peeling of CFRP and a thin layer of brick, as shown in Figure 4, the fracture of the brick dominated the failure of FRPto-brick interface. The test results are in agreement with the related literature [30,41]. The load-slip curves at the loading end are shown in Figure 5. The loading process can be roughly divided into three stages: linear growth stage, softening stage and stable stage. With the increase of slip, the tensile force first goes through the linear growth stage and the softening stage to reach the ultimate load Fu, then into the stable stage. In the stable stage, the slip further increases, the tensile force remains basically unchanged.  Comparing the average of the two groups of curves, it can be found that the linear growth rate of both groups is the same, but the linear growth stage of group B is longer so that the ultimate load is larger, which is 1.2 times of that of group A. Besides, Group B also has a longer stable stage and ends up with a larger slip.
To deeply understand the process of interfacial failure, the strain at different loading stages marked in Figure 5 was analyzed. The corresponding strain contours obtained by the DIC method are shown in Figure 6, in which the process of strain develops from the loading end to the free end can be clearly seen. Before the ultimate load Fu, the strain increases while it develops from the loading end to the free end. After the ultimate load Fu, the strain continues to develop towards the free end while the strain at the loading end hardly increases, indicating that the interface at the loading end may have been damaged.  Comparing the average of the two groups of curves, it can be found that the linear growth rate of both groups is the same, but the linear growth stage of group B is longer so that the ultimate load is larger, which is 1.2 times of that of group A. Besides, Group B also has a longer stable stage and ends up with a larger slip.
To deeply understand the process of interfacial failure, the strain at different loading stages marked in Figure 5 was analyzed. The corresponding strain contours obtained by the DIC method are shown in Figure 6, in which the process of strain develops from the loading end to the free end can be clearly seen. Before the ultimate load F u , the strain increases while it develops from the loading end to the free end. After the ultimate load F u , the strain continues to develop towards the free end while the strain at the loading end hardly increases, indicating that the interface at the loading end may have been damaged. Comparing the average of the two groups of curves, it can be found that the linear growth rate of both groups is the same, but the linear growth stage of group B is longer so that the ultimate load is larger, which is 1.2 times of that of group A. Besides, Group B also has a longer stable stage and ends up with a larger slip.
To deeply understand the process of interfacial failure, the strain at different loading stages marked in Figure 5 was analyzed. The corresponding strain contours obtained by the DIC method are shown in Figure 6, in which the process of strain develops from the loading end to the free end can be clearly seen. Before the ultimate load Fu, the strain increases while it develops from the loading end to the free end. After the ultimate load Fu, the strain continues to develop towards the free end while the strain at the loading end hardly increases, indicating that the interface at the loading end may have been damaged.  The strain on the median line marked in Figure 6 was extracted. The strain distribution data and the fitted curves of specimen A4 and B2 are shown in Figure 7a,b, respectively. The strain development at the loading end (X = 160 mm) matches their load-slip relationships shown in Figure 5. The development of strain from the loading end to the free end can be clearly seen in Figure 7. By comparing Figure 7a,b, it can be found that the final strain of B2 (1000 mm/min) is larger than that of A4 (10 mm/min). There are irregular fluctuations in strain measurements, as shown in Figure 7. It is because of the irregularities in the surface treatment and the random distribution of cracks and aggregates in brick. The formula for fitting the curves in the figure will be given in the next section along with the derivation of the bond-slip relationship. The strain on the median line marked in Figure 6 was extracted. The strain distribution data and the fitted curves of specimen A4 and B2 are shown in Figure 7a,b, respectively. The strain development at the loading end (X = 160 mm) matches their load-slip relationships shown in Figure 5. The development of strain from the loading end to the free end can be clearly seen in Figure 7. By comparing Figure 7a,b, it can be found that the final strain of B2 (1000 mm/min) is larger than that of A4 (10 mm/min). There are irregular fluctuations in strain measurements, as shown in Figure 7. It is because of the irregularities in the surface treatment and the random distribution of cracks and aggregates in brick. The formula for fitting the curves in the figure will be given in the next section along with the derivation of the bond-slip relationship.

Bond-slip Relationship
The slip distribution data derived by integrating the strain is shown in Figure 8. To eliminate the chance errors of the experimental measurements and improve the accuracy of the analysis results, the method proposed by Wu et al. [36,39,40] was applied to calculate the bond-slip relationship. The following mathematical function was used to fit the slip distribution curves at different loading stages simultaneously: where s is the slip, x is the distance from the free end, α and β are the shared fitting parameters and x0 is the fitting parameter that varies with different loading stages. The fitted results are shown in Figure 8.

Bond-Slip Relationship
The slip distribution data derived by integrating the strain is shown in Figure 8. To eliminate the chance errors of the experimental measurements and improve the accuracy of the analysis results, the method proposed by Wu et al. [36,39,40] was applied to calculate the bond-slip relationship. The following mathematical function was used to fit the slip distribution curves at different loading stages simultaneously: where s is the slip, x is the distance from the free end, α and β are the shared fitting parameters and x 0 is the fitting parameter that varies with different loading stages. The fitted results are shown in Figure 8. The strain distribution can be derived from the slip distribution [36]: Then, the distribution of the bond stress can be obtained as follows [36]: where E f t f is called the FRP stiffness, E f is the FRP modulus of elasticity and t f is the FRP thickness. The strain distribution can be derived from the slip distribution [36]: Then, the distribution of the bond stress can be obtained as follows [36]: where ff Et is called the FRP stiffness, f E is the FRP modulus of elasticity and f t is the FRP thickness.
Equation (1) can be expressed as [36]: Substituting Equation (4) into Equation (3) yields the bond-slip relationship [36]: Further, parameters such as the maximum bond stress max , the local slip 0 s at the maximum bond stress and the interfacial fracture energy f G can be derived [40]: The bond-slip parameters obtained by using the aforementioned method are listed in Table 2, and the bond-slip curves for Specimens A4 and B2 are shown in Figure 9. It can be found that the maximum bond stress and the corresponding slip are greater for Group B (1,000 mm/min) compared with Group A (10 mm/min), and the curve coverage area i.e., interfacial fracture energy is greater for Group A as well. Equation (1) can be expressed as [36]: Substituting Equation (4) into Equation (3) yields the bond-slip relationship [36]: Further, parameters such as the maximum bond stress τ max , the local slip s 0 at the maximum bond stress and the interfacial fracture energy G f can be derived [40]: The bond-slip parameters obtained by using the aforementioned method are listed in Table 2, and the bond-slip curves for Specimens A4 and B2 are shown in Figure 9. It can be found that the maximum bond stress and the corresponding slip are greater for Group B (1000 mm/min) compared with Group A (10 mm/min), and the curve coverage area i.e., interfacial fracture energy is greater for Group A as well.
It should be noted that the interfacial fracture energy is the energy consumed for FRP peeling from the substrate per unit area of FRP, not the fracture energy of cracks in the substrate material. The fracture energy of cracks in the substrate material is strain-rateindependent [42,43]. However, the corresponding fracture surface under the FRP is not flat, as shown in Figure 4, and it may change with the loading rate [18]. The variation of the fracture surface may be responsible for the dynamic enhancement of the interface, but the specific microscopic mechanism requires further investigation based on more specimens.
By comparing the load-slip curves, strain distribution, slip distribution and bond-slip curves of group A and group B, it is indicated that the interface behavior of FRP-to-brick has a significant dynamic enhancement effect. This may be related to the strain rate effect on the mechanical performance of materials such as FRP and brick under dynamic loading.   It should be noted that the interfacial fracture energy is the energy consumed for FRP peeling from the substrate per unit area of FRP, not the fracture energy of cracks in the substrate material. The fracture energy of cracks in the substrate material is strain-rateindependent [42,43]. However, the corresponding fracture surface under the FRP is not flat, as shown in Figure 4, and it may change with the loading rate [18]. The variation of the fracture surface may be responsible for the dynamic enhancement of the interface, but the specific microscopic mechanism requires further investigation based on more specimens.
By comparing the load-slip curves, strain distribution, slip distribution and bondslip curves of group A and group B, it is indicated that the interface behavior of FRP-tobrick has a significant dynamic enhancement effect. This may be related to the strain rate effect on the mechanical performance of materials such as FRP and brick under dynamic loading.

Numerical Analysis and Validation
To further investigate the bond-slip relationship at the FRP-to-brick interface under various loading rates, a plastic-damage model was established to simulate the deformation and damage of brick in the finite element software ABAQUS. The strain rate effect on the materials was considered, and the debonding process of the FRP-to-brick interface was simulated and compared with the test results.

Numerical Modeling
There are two finite element methods to model the bond behavior: the interface modeling method and the direct modeling method. The method using interface elements (such as CZM, cohesive zone method) is commonly adopted for modeling FRP-reinforce structures due to heavy computational demands. A constitutive law for the interface elements must be obtained previously. This is thus not really a predictive approach for the bond behavior but for structural behavior [44]. The other is the direct modeling method, in which debonding is simulated by modeling the failure of the substrate adjacent to FRP.
As described in Section 2.3, the fracture of the brick dominated the failure of the FRP-tobrick interface. Therefore, the direct modeling method, which assumes a perfect adhesion between FRP and substrate, was used to simulate the FRP-to-brick interface [42][43][44]. The plastic-damage model was used to simulate the deformation and damage of the brick [45]. Thus the slippage of the interface and the damage of the brick are simulated.
A three-dimensional numerical model for the single-lap shear test was constructed according to the dimensions of the specimen and the boundary conditions of the test, as shown in Figure 10. Due to the symmetry of the structure and load, only half of the model was built. The brick was modeled with C3D8 eight-node hexahedral elements, and the FRP was modeled with S4R four-node shell element. The FRP and brick near the interface use the same mesh size, and the corresponding nodes were tied separately, so a perfect adhesion between FRP and brick was fulfilled. The bottom and right side faces of the brick were fixed.
There are two finite element methods to model the bond behavior: the interface modeling method and the direct modeling method. The method using interface elements (such as CZM, cohesive zone method) is commonly adopted for modeling FRP-reinforce structures due to heavy computational demands. A constitutive law for the interface elements must be obtained previously. This is thus not really a predictive approach for the bond behavior but for structural behavior [44]. The other is the direct modeling method, in which debonding is simulated by modeling the failure of the substrate adjacent to FRP.
As described in Section 2.3, the fracture of the brick dominated the failure of the FRPto-brick interface. Therefore, the direct modeling method, which assumes a perfect adhesion between FRP and substrate, was used to simulate the FRP-to-brick interface [42][43][44]. The plastic-damage model was used to simulate the deformation and damage of the brick [45]. Thus the slippage of the interface and the damage of the brick are simulated.
A three-dimensional numerical model for the single-lap shear test was constructed according to the dimensions of the specimen and the boundary conditions of the test, as shown in Figure 10. Due to the symmetry of the structure and load, only half of the model was built. The brick was modeled with C3D8 eight-node hexahedral elements, and the FRP was modeled with S4R four-node shell element. The FRP and brick near the interface use the same mesh size, and the corresponding nodes were tied separately, so a perfect adhesion between FRP and brick was fulfilled. The bottom and right side faces of the brick were fixed. The mechanical behavior of FRP was described with a linear elastic model. The mechanical behavior of clay bricks is similar to that of concrete. There is a stiffness degradation caused by the propagation of micro-cracks and micro-defects, and there is a plastic deformation such as slip and flow associated with the mesoscopic mechanism of materials The mechanical behavior of FRP was described with a linear elastic model. The mechanical behavior of clay bricks is similar to that of concrete. There is a stiffness degradation caused by the propagation of micro-cracks and micro-defects, and there is a plastic deformation such as slip and flow associated with the mesoscopic mechanism of materials deformation. The correct constitutive relationship of clay bricks should be a plastic-damage constitutive that reflects both damage and plastic deformation.

Plastic-Damage Model
The plastic-damage model in ABAQUS was improved and developed by Lee and Fenves [45] based on Lubliner [46]. The model uses the damage combined with plastic strain to simulate the inelastic behavior of clay bricks. Since the stiffness degradation caused by tensile and compressive plastic strains are considered, this model is suitable for clay bricks under complex loading conditions.
The material failure is achieved by stiffness degradation in the plastic-damage model. Isotropic stiffness degradation is assumed. The concepts of effective stress and strain decomposition are employed. Then the Cauchy stress tensor can be expressed as [45]: where the damage factor consists of tensile damage d t and compression damage d c , σ is the effective stress tensor, E 0 is the initial elastic modulus, ε is the strain tensor and ε p is the plastic component of the strain tensor.
The nonassociative flow rule is adopted in the model. The plastic potential function associated with effective stress is in the Drucker-Prager form [45]: where I 1 is the first invariant of the stress tensor, J 2 is the second invariant of the stress deviator, f t0 is the initial uniaxial tensile strength, ξ is the smoothing parameter and ψ is the dilatancy angle. The plastic parameters used in this paper are listed in Table 3, where f b0 /f c0 and K are yield parameters. A detailed explanation of these parameters can be found in the relevant literature [45,47]. The axial compressive strength f c of MU15 clay bricks under quasi-static is 11.4 MPa and the elastic modulus is 7630 MPa. The shape of the compression curve refers to GB50010-2010 [48]. Few experimental studies focused on the tensile properties of clay bricks, the tensile strength (0.25 f c ) and curve shape; refer to the study of D'Altri et al. [49].
To overcome the mesh sensitivity problem due to the softening behavior and stiffness degradation of the plastic-damage model [47,49], regularization of this model was achieved by scaling the fracture energy with equivalent length. A detailed description of this can be found in the relevant literature [44]. Scaled fracture energy is noted to affect the crack speed in some studies [50] because no damping and no friction is included in the interactions between damaged material points. However, the damage results are not affected by the mesh size [51]. In this paper, the material strength is set according to the loading rate, and the simulation results for the bonding behavior are not influenced by the crack speed.

Strain Rate Effects on Materials
The FRP did not fail during interfacial debonding, so this study only considers the strain rate effect on its elastic modulus. Zhang et al. [52] found that the strain rate effect on the elastic modulus of unidirectional CFRP sheets has obvious stage differences, and the dynamic enhancement factor (DIF) is: The compressive strength of bricks is a key parameter that influences interface performance. Zhang et al. [53] carried out the compressive strength tests of clay bricks under different strain rates and obtained the dynamic enhancement factor of compressive strength:

Validation with Test Results
By comparing the load-slip curves at the loading end simulated with different element sizes, as shown in Figure 11, the mesh convergence was analyzed, and 2 mm was adopted by considering both computational efficiency and accuracy of results. The aforementioned numerical model was used to simulate the single-lap shear tests under the loading rates of 10 mm/min and 1000 mm/min, respectively. The calculated load-slip curves at the loading end were compared with the test results, as shown in Figure 12, where the test results are the mean and standard deviation over the common slip range. It can be found that the linear growth stage and the stable stage of the numerical results are in good agreement with the test results, which proves the reliability of the numerical model.

Validation with Test Results
By comparing the load-slip curves at the loading end simulated with different element sizes, as shown in Figure 11, the mesh convergence was analyzed, and 2 mm was adopted by considering both computational efficiency and accuracy of results. The aforementioned numerical model was used to simulate the single-lap shear tests under the loading rates of 10 mm/min and 1000 mm/min, respectively. The calculated load-slip curves at the loading end were compared with the test results, as shown in Figure 12, where the test results are the mean and standard deviation over the common slip range. It can be found that the linear growth stage and the stable stage of the numerical results are in good agreement with the test results, which proves the reliability of the numerical model.  In the same way as analyzing the test data, firstly, the strain distribution data at different loading stages marked in Figure 12 were extracted, as shown in Figure 13a,b. Then the slip distribution data were obtained by integrating the strain, as shown in Figure 13c,d. Finally, the seven slip distribution curves were fitted by Equation (1), and the fitted parameters α and β were substituted into Equation (5)  In the same way as analyzing the test data, firstly, the strain distribution data at different loading stages marked in Figure 12 were extracted, as shown in Figure 13a,b. Then the slip distribution data were obtained by integrating the strain, as shown in Figure 13c,d. Finally, the seven slip distribution curves were fitted by Equation (1), and the fitted parameters α and β were substituted into Equation (5) to derive the bond-slip relationship.
The numerical results of the bond-slip curves were compared with the test results, as shown in Figure 14. As the feature parameters of the bond-slip relationship, the errors of the numerical results for τ max , s 0 and G f relative to the test results were calculated. The errors under the loading rate of 10 mm/min were −3.5%, 0.8% and −3.2%, respectively; the errors of 1000 mm/min were −3.7%, −1.6% and −5.5%, respectively.
The strain rate effect on the material strength used in this paper is a macroscopic mechanical model. The generation and development of cracks are not simulated. Therefore, the dynamic enhancement of the interfacial fracture energy in the numerical simulations is mainly due to the higher brick strength being considered at higher loading rates.
In the same way as analyzing the test data, firstly, the strain distribution data at different loading stages marked in Figure 12 were extracted, as shown in Figure 13a,b. Then the slip distribution data were obtained by integrating the strain, as shown in Figure 13c,d. Finally, the seven slip distribution curves were fitted by Equation (1), and the fitted parameters α and β were substituted into Equation (5) to derive the bond-slip relationship.   The strain rate effect on the material strength used in this paper is a macroscopic mechanical model. The generation and development of cracks are not simulated. Therefore, the dynamic enhancement of the interfacial fracture energy in the numerical simulations is mainly due to the higher brick strength being considered at higher loading rates.

Debonding Analysis
The interface debonding process was analyzed by taking a specimen of 10 mm/min as an example. As shown in Figure 15a, the slip distribution, strain distribution and bond stress distribution at different loading stages can be derived by substituting the fitted parameters α, β and the corresponding x 0 into Equations (1)-(3), respectively. Overall, with increasing loading, the interfacial debonds from the loading end to the free end. The strain and slip of FRP develop toward the free end, while the interfacial shear bond stress concentration region moves toward the free end. The slip shows an exponentially monotonically increasing distribution during debonding. The shear bond stress shows a similar distribution as Gaussian, which also can be found from the shear bond stress contours shown in Figure 15b. The maximum bond stress in the contours corresponds one-to-one with the peaks in the bond stress distribution, as indicated by the dashed line in Figure 15. In addition, the shear stress concentration caused by squeezing occurs near the loading end, which is the result of the combined effect of boundary conditions and loading conditions.
The numerical results of the load-slip and bond-slip relationships are in good agreement with the test results, and the numerical results give a better understanding of the development of interface debonding. Therefore, the numerical model based on the plasticdamage constitutive and the strain rate effects on materials can simulate the FRP-to-brick interface behavior under different loading rates. The interface debonding process was analyzed by taking a specimen of 10 mm/min as an example. As shown in Figure 15a, the slip distribution, strain distribution and bond stress distribution at different loading stages can be derived by substituting the fitted parameters α, β and the corresponding 0 into Equations (1)-(3), respectively. Overall, with increasing loading, the interfacial debonds from the loading end to the free end. The strain and slip of FRP develop toward the free end, while the interfacial shear bond stress concentration region moves toward the free end. The slip shows an exponentially monotonically increasing distribution during debonding. The shear bond stress shows a similar distribution as Gaussian, which also can be found from the shear bond stress contours shown in Figure 15b. The maximum bond stress in the contours corresponds one-to-one with the peaks in the bond stress distribution, as indicated by the dashed line in Figure 15. In addition, the shear stress concentration caused by squeezing occurs near the loading end, which is the result of the combined effect of boundary conditions and loading conditions. The numerical results of the load-slip and bond-slip relationships are in good agreement with the test results, and the numerical results give a better understanding of the development of interface debonding. Therefore, the numerical model based on the plasticdamage constitutive and the strain rate effects on materials can simulate the FRP-to-brick interface behavior under different loading rates.

Parameter Study on Dynamic Behavior
Numerous experimental studies of FRP-substrate interfaces have shown that the FRP stiffness and the substrate strength are the main parameters affecting the bond-slip rela-

Parameter Study on Dynamic Behavior
Numerous experimental studies of FRP-substrate interfaces have shown that the FRP stiffness and the substrate strength are the main parameters affecting the bond-slip relationship [20,29,34,38,41,54]. In this section, firstly, the effects of single variables of FRP stiffness and brick strength on the bond-slip relationship at the FRP-to-brick interface were investigated. Then the strain rate effects on FRP stiffness and brick strength were taken into account together to study the bond-slip relationship under different slip rates.

Effect of FRP Stiffness
In the study of the bond-slip relationship, the product of the elastic modulus E f and thickness t f of FRP is generally called the FRP stiffness. The FRP stiffness is mainly determined by FRP type and thickness and may be affected by the loading strain rate. Therefore, different FRPs (CFRP: carbon fiber reinforced polymer, GFRP: glass fiber reinforced polymer and BFRP: basalt fiber reinforced polymer), number of layers (1ply: one layer, 2ply: two layers and 3ply: three layers) and dynamic increase factor (1.2, 1.4, 1.6 and 1.8) were considered in the parameter study. The uniaxial compressive strength of the bricks was set to 11.4 MPa, and the numerical model was used to calculate the bond-slip relationship with different FRP stiffness. The specific calculation settings and results are shown in Table 4. It should be noted that the E f t f , α, β, τ max , s 0 and G f in Table 4 are the standardized parameters of E f t f , α, β, τ max , s 0 and G f , respectively. Specimen CFRP-1ply was used as a reference, so all the standardized parameters of Specimen CFRP-1ply are 1. The purpose of parameter standardization is to facilitate the observation of parameter sensitivity. The load-slip curves with different FRP stiffness are shown in Figure 16. As the FRP stiffness increases, it can be found that the ultimate load increases significantly, while the final slip decreases significantly with no obvious stable stage in CFRP-3ply. This phenomenon can be explained by the two extreme cases of FRP stiffness. By assuming that the FRP stiffness is infinitely large, i.e., the FRP is rigid, the slip at the loading end is identical to the slip at the free end. All the regions in contact with the FRP are simultaneously stressed, so the ultimate load will be large and the final slip will be small. On the contrary, if the FRP stiffness is infinitesimal, the deformation of the FRP at the loading end cannot develop towards the free end. Only the regions contacted by the FRP at the loading end will be stressed, so the ultimate load will be small and the final slip will be large. Moreover, from the point of view of energy conservation, the work done by the tensile force should be equal to the fracture energy of the interface debonding under static loading. If the area of the fracture surface keeps constant, the work done by the tensile force should be fixed, so the ultimate load increases while the final slip decreases. The same phenomenon can be found in the corresponding bond-slip curves, as shown in Figure 17, where the maximum bond stress increases with increasing FRP stiffness while the corresponding slip at the maximum bond stress decreases. the area of the fracture surface keeps constant, the work done by the tensile force should be fixed, so the ultimate load increases while the final slip decreases. The same phenomenon can be found in the corresponding bond-slip curves, as shown in Figure 17, where the maximum bond stress increases with increasing FRP stiffness while the corresponding slip at the maximum bond stress decreases.  A specific understanding of the effects of FRP stiffness on the bond-slip relationship can be found by further quantitative analysis of the bond-slip parameters. The relationships between standardized α, β and standardized FRP stiffness are shown in Figure 18a. As the FRP stiffness increases, α decreases while β increases, and both the standardized parameters of α and β have a power function relationship with the standardized FRP stiffness. The changes in interfacial fracture energy, maximum bond stress and the corresponding slip with FRP stiffness can also be calculated according to Equation (6). As the area of the fracture surface keeps constant, the work done by the tensile force should be fixed, so the ultimate load increases while the final slip decreases. The same phenomenon can be found in the corresponding bond-slip curves, as shown in Figure 17, where the maximum bond stress increases with increasing FRP stiffness while the corresponding slip at the maximum bond stress decreases.  A specific understanding of the effects of FRP stiffness on the bond-slip relationship can be found by further quantitative analysis of the bond-slip parameters. The relationships between standardized α, β and standardized FRP stiffness are shown in Figure 18a. As the FRP stiffness increases, α decreases while β increases, and both the standardized parameters of α and β have a power function relationship with the standardized FRP stiffness. The changes in interfacial fracture energy, maximum bond stress and the corresponding slip with FRP stiffness can also be calculated according to Equation (6). As A specific understanding of the effects of FRP stiffness on the bond-slip relationship can be found by further quantitative analysis of the bond-slip parameters. The relationships between standardized α, β and standardized FRP stiffness are shown in Figure 18a. As the FRP stiffness increases, α decreases while β increases, and both the standardized parameters of α and β have a power function relationship with the standardized FRP stiffness. The changes in interfacial fracture energy, maximum bond stress and the corresponding slip with FRP stiffness can also be calculated according to Equation (6). As shown in Figure 18b, the changes in maximum bond stress and the corresponding slip have the opposite trend and the same magnitude; thus, the interfacial fracture energy remains unchanged. This is consistent with the analysis of the load-slip curves and bond-slip curves.
shown in Figure 18b, the changes in maximum bond stress and the corresponding slip have the opposite trend and the same magnitude; thus, the interfacial fracture energy remains unchanged. This is consistent with the analysis of the load-slip curves and bondslip curves. Therefore, FRP stiffness mainly influences the shape of the bond-slip curve. The maximum bond stress increases with increasing FRP stiffness, and the corresponding slip decreases accordingly; thus the interfacial fracture energy keeps constant.

Effect of Brick Strength
The brick strength mentioned in this paper refers specifically to the axial compressive strength. The quasi-static strength of MU15 clay bricks is 11.4 MPa, and its strength at the strain rate of 200 s −1 is 25.57 MPa calculated by Equation (10). The FRP stiffness is consistent with CFRP-1ply in Table 3 and keeps constant. The bond-slip relationships under nine different brick strength in the range of 11.4 MPa to 25.57 MPa were calculated, the specific calculation settings and results are shown in Table 5. The load-slip curves with different brick strengths are shown in Figure 19. As the brick strength increases, it can be observed that both the ultimate load and the final slip increase significantly, which indicates that the work done by the tensile force increases. Therefore, FRP stiffness mainly influences the shape of the bond-slip curve. The maximum bond stress increases with increasing FRP stiffness, and the corresponding slip decreases accordingly; thus the interfacial fracture energy keeps constant.

Effect of Brick Strength
The brick strength mentioned in this paper refers specifically to the axial compressive strength. The quasi-static strength of MU15 clay bricks is 11.4 MPa, and its strength at the strain rate of 200 s −1 is 25.57 MPa calculated by Equation (10). The FRP stiffness is consistent with CFRP-1ply in Table 3 and keeps constant. The bond-slip relationships under nine different brick strength in the range of 11.4 MPa to 25.57 MPa were calculated, the specific calculation settings and results are shown in Table 5. The load-slip curves with different brick strengths are shown in Figure 19. As the brick strength increases, it can be observed that both the ultimate load and the final slip increase significantly, which indicates that the work done by the tensile force increases. The same phenomenon can be found in the corresponding bond-slip curves, as shown in Figure 20, where both maximum bond stress and the corresponding slip increase with increasing brick strength, making the curve coverage area, i.e., interfacial fracture energy, increase significantly.
The same phenomenon can be found in the corresponding bond-slip curves, as shown in Figure 20, where both maximum bond stress and the corresponding slip increase with increasing brick strength, making the curve coverage area, i.e., interfacial fracture energy, increase significantly.  The relationships between standardized α, β and standardized brick strength are shown in Figure 21a. With the increase in brick strength, α increases significantly, while β decreases slightly, and both the standardized parameters of α and β have a power function relationship with the standardized brick strength. The changes in interfacial fracture energy, maximum bond stress and the corresponding slip with brick strength can also be calculated according to Equation (6). As shown in Figure 21b, variations of the maximum bond stress and the corresponding slip have a similar trend, which indicates that the shape of the bond-slip curve remains basically unchanged. Therefore, the interfacial fracture energy increases significantly with brick strength. The same phenomenon can be found in the corresponding bond-slip curves, as shown in Figure 20, where both maximum bond stress and the corresponding slip increase with increasing brick strength, making the curve coverage area, i.e., interfacial fracture energy, increase significantly.  The relationships between standardized α, β and standardized brick strength are shown in Figure 21a. With the increase in brick strength, α increases significantly, while β decreases slightly, and both the standardized parameters of α and β have a power function relationship with the standardized brick strength. The changes in interfacial fracture energy, maximum bond stress and the corresponding slip with brick strength can also be calculated according to Equation (6). As shown in Figure 21b, variations of the maximum bond stress and the corresponding slip have a similar trend, which indicates that the shape of the bond-slip curve remains basically unchanged. Therefore, the interfacial fracture energy increases significantly with brick strength. The relationships between standardized α, β and standardized brick strength are shown in Figure 21a. With the increase in brick strength, α increases significantly, while β decreases slightly, and both the standardized parameters of α and β have a power function relationship with the standardized brick strength. The changes in interfacial fracture energy, maximum bond stress and the corresponding slip with brick strength can also be calculated according to Equation (6). As shown in Figure 21b, variations of the maximum bond stress and the corresponding slip have a similar trend, which indicates that the shape of the bond-slip curve remains basically unchanged. Therefore, the interfacial fracture energy increases significantly with brick strength.
Being different from FRP stiffness, the brick strength mainly influences the amplitude of the bond-slip curve. Both the maximum bond stress and the corresponding slip increase as the brick strength increases, and the shape of the bond-slip curve remains almost unchanged; thus, the interfacial fracture energy increases significantly. Being different from FRP stiffness, the brick strength mainly influences the amplitude of the bond-slip curve. Both the maximum bond stress and the corresponding slip increase as the brick strength increases, and the shape of the bond-slip curve remains almost unchanged; thus, the interfacial fracture energy increases significantly.

Empirical Formulas and Validation
Based on the effects of FRP stiffness and brick strength on the bond-slip relationship in the above two sections, empirical formulas for α and β were derived by combining the fit results in Figures 18a and 21a (11) Then the empirical formulas for the maximum bond stress, the corresponding slip at the maximum bond stress, and the interfacial fracture energy were obtained with Equation (6): The strain rate effects on both FRP stiffness and brick strength were taken into account. The bond-slip relationships were calculated with different slip rates. The slip rate here is the rate of local slip at the interface, which is equivalent to the loading rate at the loading end. For easy parameter conversion, the slip rate is given in m/s. The settings and calculation results under different slip rates are shown in Table 6, where the slip rate ṡ ranges from 1.67E-6 m/s to 10 m/s. As shown in Figure 22, the strain history of FRP under

Empirical Formulas and Validation
Based on the effects of FRP stiffness and brick strength on the bond-slip relationship in the above two sections, empirical formulas for α and β were derived by combining the fit results in Figures 18a and 21a: Then the empirical formulas for the maximum bond stress, the corresponding slip at the maximum bond stress, and the interfacial fracture energy were obtained with Equation (6): Where the units of α, β, t f and s 0 are mm; the units of, f c , E f and τ max are MPa; and the unit of G f is N/mm. The strain rate effects on both FRP stiffness and brick strength were taken into account. The bond-slip relationships were calculated with different slip rates. The slip rate here is the rate of local slip at the interface, which is equivalent to the loading rate at the loading end. For easy parameter conversion, the slip rate is given in m/s. The settings and calculation results under different slip rates are shown in Table 6, where the slip rate . s ranges from 1.67 × 10 −6 m/s to 10 m/s. As shown in Figure 22, the strain history of FRP under the slip rate of 10 m/s shows that the slope of the linearly rising section of the FRP strain at different locations is almost identical. In this study, the slope of the linearly rising section at X = 160 mm was taken as the strain rate and the corresponding strain rate . ε ranges from 1.96 × 10 −5 s -1 to 118.57 s -1 .
The relationship between the slip rate and the strain rate is initially obtained from the single-lap shear tests. Then, the material parameters for different strain rates and the corresponding slip rates were set in the numerical simulations. The numerical results under different slip rates further validate the relationship between the slip rate and the strain rate, which was expressed by Equation (13)  The dynamic FRP stiffness and dynamic brick strength at a given slip rate can be calculated by combining Equations (9), (10) and (13). Then, the bond-slip parameter at a given slip rate can be calculated by substituting the dynamic FRP stiffness and dynamic brick strength into the empirical formulas (11) and (12). The maximum bond stress and the corresponding slip calculated by the empirical formulas were compared with the test and numerical results. As shown in Figure 23, the prediction results of the empirical formulas are in good agreement with the test and numerical results, indicating the reliability of the empirical formulas for calculating the bond-slip relationship at the FRP-to-brick interface under dynamic loading. the slip rate of 10 m/s shows that the slope of the linearly rising section of the FRP strain at different locations is almost identical. In this study, the slope of the linearly rising section at X = 160 mm was taken as the strain rate and the corresponding strain rate ε̇ ranges from 1.96E-5 s -1 to 118.57 s -1 .  The relationship between the slip rate and the strain rate is initially obtained from the single-lap shear tests. Then, the material parameters for different strain rates and the corresponding slip rates were set in the numerical simulations. The numerical results under different slip rates further validate the relationship between the slip rate and the strain rate, which was expressed by Equation (13).
  10 10 log log 1.0687 s (13) The dynamic FRP stiffness and dynamic brick strength at a given slip rate can be calculated by combining Equations (9), (10) and (13). Then, the bond-slip parameter at a given slip rate can be calculated by substituting the dynamic FRP stiffness and dynamic brick strength into the empirical formulas (11) and (12). The maximum bond stress and the corresponding slip calculated by the empirical formulas were compared with the test and numerical results. As shown in Figure 23, the prediction results of the empirical formulas are in good agreement with the test and numerical results, indicating the reliability of the empirical formulas for calculating the bond-slip relationship at the FRP-to-brick interface under dynamic loading.

Analysis and Discussion
The load-slip curves and bond-slip curves under different slip rates are shown in Figures 24 and 25, respectively. The variations of the curves under different slip rates are basically similar to those with different brick strength; that is, the curve shape remains unchanged and the amplitude increases. This indicates that the influence of the slip rate on the bond-slip relationship is dominated by the brick strength. The specific reasons can be analyzed with the empirical formulas and the dynamic enhancement law of material performance.

Analysis and Discussion
The load-slip curves and bond-slip curves under different slip rates are shown in Figures 24 and 25, respectively. The variations of the curves under different slip rates are basically similar to those with different brick strength; that is, the curve shape remains unchanged and the amplitude increases. This indicates that the influence of the slip rate on the bond-slip relationship is dominated by the brick strength. The specific reasons can be analyzed with the empirical formulas and the dynamic enhancement law of material performance.

Analysis and Discussion
The load-slip curves and bond-slip curves under different slip rates are shown in Figures 24 and 25, respectively. The variations of the curves under different slip rates are basically similar to those with different brick strength; that is, the curve shape remains unchanged and the amplitude increases. This indicates that the influence of the slip rate on the bond-slip relationship is dominated by the brick strength. The specific reasons can be analyzed with the empirical formulas and the dynamic enhancement law of material performance.  On one hand, the comparison between Figures 18b and 21b shows that the maximum bond stress and the corresponding slip are more sensitive to brick strength than FRP stiffness, and the interfacial fracture energy is only affected by brick strength. On the other hand, the dynamic enhancement effect on brick strength is always greater than FRP stiffness, as shown in Table 4 (9) and (10) are shown in Figure 26. In the strain rate range of 20-76 s −1 (the slip rate is 1.71-6.49 m/s), the strain rate sensitivity (slope of the curve) of FRP stiffness is slightly stronger than that of brick strength. In this case, the bond-slip relationship is clearly influenced by FRP stiffness, which is the main reason for the shape variation of the bond-slip curve under the slip rate of 5 m/s. However, as the strain rate increases further On one hand, the comparison between Figures 18b and 21b shows that the maximum bond stress and the corresponding slip are more sensitive to brick strength than FRP stiffness, and the interfacial fracture energy is only affected by brick strength. On the other hand, the dynamic enhancement effect on brick strength is always greater than FRP stiffness, as shown in Table 4, where the DIF f c is 1.735 at the strain rate of 19.6 s −1 , while the DIF E f t f is only 1.028. Noteworthy, the dynamic enhancement laws of FRP stiffness and brick strength in Equations (9) and (10) are shown in Figure 26. In the strain rate range of 20-76 s −1 (the slip rate is 1.71-6.49 m/s), the strain rate sensitivity (slope of the curve) of FRP stiffness is slightly stronger than that of brick strength. In this case, the bond-slip relationship is clearly influenced by FRP stiffness, which is the main reason for the shape variation of the bond-slip curve under the slip rate of 5 m/s. However, as the strain rate increases further and becomes greater than 76 s −1 (the slip rate is greater than 6.49 m/s), the strain rate sensitivity of brick strength increases significantly and becomes much greater than the strain rate sensitivity of FRP stiffness, as shown in Figure 26. The dynamic enhancement effect on brick strength is therefore more pronounced and again dominates the variations in the bond-slip relationship under dynamic loading, as can be seen from the bond-slip curve under the slip rate of 10 m/s. On one hand, the comparison between Figures 18b and 21b shows that the maximum bond stress and the corresponding slip are more sensitive to brick strength than FRP stiffness, and the interfacial fracture energy is only affected by brick strength. On the other hand, the dynamic enhancement effect on brick strength is always greater than FRP stiffness, as shown in Table 4, where the c f DIF is 1.735 at the strain rate of 19.6 s −1 , while the ff Et DIF is only 1.028.
Noteworthy, the dynamic enhancement laws of FRP stiffness and brick strength in Equations (9) and (10) are shown in Figure 26. In the strain rate range of 20-76 s −1 (the slip rate is 1.71-6.49 m/s), the strain rate sensitivity (slope of the curve) of FRP stiffness is slightly stronger than that of brick strength. In this case, the bond-slip relationship is clearly influenced by FRP stiffness, which is the main reason for the shape variation of the bond-slip curve under the slip rate of 5 m/s. However, as the strain rate increases further and becomes greater than 76 s −1 (the slip rate is greater than 6.49 m/s), the strain rate sensitivity of brick strength increases significantly and becomes much greater than the strain rate sensitivity of FRP stiffness, as shown in Figure 26. The dynamic enhancement effect on brick strength is therefore more pronounced and again dominates the variations in the bond-slip relationship under dynamic loading, as can be seen from the bond-slip curve under the slip rate of 10 m/s.  In summary, the variations of the bond-slip relationship under dynamic loading are mainly caused by the dynamic enhancement of brick strength. In other words, the amplitude of the bond-slip curve increases with increasing slip rate, and the shape of the curve almost remains unchanged. However, the bond-slip relationship is also clearly influenced by FRP stiffness especially within a specific slip rate range, where the strain rate sensitivity of FRP stiffness is stronger than that of brick strength.

Conclusions
In this paper, the effects of FRP stiffness, brick strength and slip rate on the bond-slip relationship at FRP-to-brick interfaces are numerically investigated based on the single-lap shear tests with two different loading rates. The conclusions are as follow: 1.
The numerical model based on the plastic-damage constitutive and the strain rate effects on material performance can simulate the FRP-to-brick interface behavior under different loading rates, and the numerical results of the bond-slip relationship are in good agreement with the test results.

2.
FRP stiffness mainly influences the shape of the bond-slip curve. As the maximum bond stress increases with increasing FRP stiffness, the corresponding slip at the maximum bond stress decreases and the interfacial fracture energy remains constant.

3.
Brick strength mainly influences the amplitude of the bond-slip curve. Both the maximum bond stress and the corresponding slip increase as the brick strength increases, and the shape of the bond-slip curve remains almost unchanged; thus, the interfacial fracture energy increases significantly.

4.
The variations of the bond-slip relationship under dynamic loading are mainly a consequence of the dynamic enhancement of brick strength, so the magnitude of the bond-slip curve changes significantly. However, the dynamic bond-slip relationship is also influenced by FRP stiffness especially within a specific slip rate range, where the strain rate sensitivity of FRP stiffness is stronger than that of brick strength.

5.
The empirical formulas considering dynamic FRP stiffness and dynamic brick strength can be used to predict the bond-slip relationship at the FRP-to-brick interface under dynamic loading, and its prediction results are in good agreement with the test and numerical results. Informed Consent Statement: Not applicable.

Data Availability Statement:
The raw/processed data required to reproduce these findings cannot be shared at this time as the data also form part of an ongoing study.