Numerical Method for Creep Analysis of Strengthened Fatigue-Damaged Concrete Beams

: Fatigue-damaged concrete improves the load-bearing capacity of components by increasing the cross section. However, the creep performance of damaged components after the repair has received less attention. Thus, this study establishes a constitutive creep model of strengthened fatigue-damaged concrete on the basis of damage mechanics and numerically simulates the strengthened component. The accuracy of the proposed model is veriﬁed by conducting creep tests on fatigue-damaged concrete beams. According to the numerical simulation results, increasing the section height profoundly affects the ability to control their creep deﬂection. The incremental creep deﬂection of the beams with a strengthened section height of 50, 100, and 150 mm loaded for 365 days decreased by 0.107, 0.228, and 0.326 mm, respectively, compared with the unstrengthened damaged beam. Moreover, this reinforcement method excellently controls the deﬂection of the damaged components under a negative bending moment. The model can forecast the creep deformation of undamaged components or damaged components after being strengthened, which facilitates structural maintenance and decision-making about reinforcement.


Introduction
Fatigue and creep of concrete structures are specific time-varying properties during the entire life cycle of engineering structures, primarily studied independently [1], with less coupled consideration. The fatigue damage phenomenon is evident in pier and bridge projects, which are often strengthened by increasing the cross section. Because of the difference between the creep coefficients of new and old concrete after reinforcement, the deflection of the components is continuously increasing. This phenomenon has been confirmed in the related literature [2][3][4]. Due to the effect of fatigue damage, the deflection and stress redistribution of strengthened components may be severe. Therefore, studying the long-term deformation of strengthened fatigue-damaged components is crucial and valuable. The mechanism for the long-term deformation can be analyzed in terms of the micro-difference in the creep of new and old concrete and the degree of fatigue damage.
The long-term performance of concrete components strengthened by the increased cross section method was first noticed in the 1980s [5,6]. Wang et al. investigated the effect of the shrinkage and creep on prestressed new-old concrete composite beams experimentally and revealed that the concrete strain, deflection, and section curvature along the beam height increased gradually with time [7]. In engineering, increasing the cross section is often used to improve the structural bearing capacity and stiffness, so it is essential to identify the effect of increasing the cross section on the creep of components. Yang et al. considered deformation coordination as a calculation condition and presented the influence of shrinkage and creep by introducing the cross section cooperative work coefficient into the classical calculation formula for concrete structure design [8]. Other works have reported where E can be considered as the elasticity modulus of the damaged material, and E is the Young's modulus.
• The effective stress theory The Cauchy stress on the initial cross section of the micro-element of the undamaged material is σ. The stresses on the initial cross section of the micro-element differ between undamaged and damaged materials. The effective stress σ on this cross section can be expressed by the Equation (2).
• The equivalent strain theory.
It is assumed that the solid material is an isotropic homogeneous elastomer and that the object is originally in a stress-free natural state. The equivalent strain theory holds that the deformation state of damaged materials under uniaxial or multiaxial stress state can be expressed by the original undamaged material intrinsic law. As shown in Figure 1, the Cauchy stress is replaced by the effective stress in the intrinsic relationship equation. merical simulation analysis of damaged components can be realized by continuously calling the USERMAT subroutine.

Constitutive Fatigue Model
The damage variable ( D ) is used to define the degree of damage to materials in continuum damage mechanics, and it is shown in Equation (1) [25]. The constitutive relationship of materials under coupled damage can be achieved through effective stress theory and equivalent strain theory.
where E can be considered as the elasticity modulus of the damaged material, and E is the Young's modulus. • The effective stress theory The Cauchy stress on the initial cross section of the micro-element of the undamaged material is σ. The stresses on the initial cross section of the micro-element differ between undamaged and damaged materials. The effective stress σ on this cross section can be expressed by the Equation (2).
The equivalent strain theory.
It is assumed that the solid material is an isotropic homogeneous elastomer and that the object is originally in a stress-free natural state. The equivalent strain theory holds that the deformation state of damaged materials under uniaxial or multiaxial stress state can be expressed by the original undamaged material intrinsic law. As shown in Figure 1, the Cauchy stress is replaced by the effective stress in the intrinsic relationship equation.  When the stress is replaced by effective stress, the constitutive relation of undamaged materials can express the strain of damaged materials [25]. Equation (3) defines the constitutive relation of coupled damage as: When the stress is replaced by effective stress, the constitutive relation of undamaged materials can express the strain of damaged materials [25]. Equation (3) defines the constitutive relation of coupled damage as: where σ ij and ε ij are the stress and strain tensors, respectively, δ ij indicates Kronecker delta, σ kk is volume stress, ν represents Poisson's ratio, and E is the elastic modulus. The concrete part of the component deteriorates under fatigue loading. This study employs the model for the fatigue damage to concrete based on damage mechanics proposed by Mai et al. [26] to define the damage state and damage distribution after fatigue loading: where D is the damage variable, .
D denotes the rate of change in the damage variable, ε represents the strain tensor, m(D) indicates the material parameters, S(ε) is the damage driving source determined by the strain tensor, h(f ) is a continuous and increasing function such that h(0) = 0 and h(1) = 1, and · + stands for the positive scalar function, namely x + = 1 2 {x + |x|}. Equation (5) defines m(D) as: where m 1 , m 2 , and m 3 are the material characteristic parameters determined through experimentation. Section 5.1 presents their specific values. Equation (6) defines h( f ) as: where α, n 1 , and n 2 indicate the material characteristic parameters measured by experimentation. Their specific values are listed in Section 5.1. Equations (7) and (8) can express f and S(ε), respectively: where W is the initial damage threshold.

Constitutive Creep Model with Fatigue Damage
The stresses on the initial cross section of the micro-element differ between undamaged and damaged materials. Considering the influence of damage, effective stress σ is used to replace the stress on the damaged material on the basis of the assumption of strain equivalence. The effective stress on damaged materials in the solution process is completed in Section 2.1. The principal covariance structure considering the initial damage is applied to the solution process of the numerical simulation.
This study adopts the creep prediction model proposed by Zhu for the numerical simulation of the creep of the concrete beam [27]. Assuming that the stress in each time period is varying linearly: where C(t, τ) is in the loading period of τ and the moment of the creep degree of t, f 1 , g 1 , p 1 , γ 1 , f 2 , g 2 , p 2 , and γ 2 represent the parameters determined by the standard creep test of concrete. Section 5.1 displays the specific values of these tabulated parameters.

Numerical Modeling of Creep Deformation with Fatigue Damage
We develop the constitutive creep model for damaged concrete using the USERMAT subroutine. The consistent tangent operator matrix is calculated for each step of the iterative process on the basis of this constitutive model. The stress increment can be deduced according to the matrix and the strain increment in the last step. Figure 2 illustrates the flow chart for the USERMAT subroutine.
where ( ) , C t τ is in the loading period of τ and the moment of the creep degree of t, 1 f , 1 g , 1 p , 1 γ , 2 f , 2 g , 2 p , and 2 γ represent the parameters determined by the standard creep test of concrete. Section 5.1 displays the specific values of these tabulated parameters.

Numerical Modeling of Creep Deformation with Fatigue Damage
We develop the constitutive creep model for damaged concrete using the USERMAT subroutine. The consistent tangent operator matrix is calculated for each step of the iterative process on the basis of this constitutive model. The stress increment can be deduced according to the matrix and the strain increment in the last step. Figure 2 illustrates the flow chart for the USERMAT subroutine. Numerical simulation is divided into two main processes: introducing the initial damage and analyzing the creep effect.
(a) Introducing initial damage Introducing fatigue damage is simulated by coding and calling the USERMAT subroutine. The coding consists of two parts:

•
Calculating the stress increment using the strain increment ( ε Δ ) and updating the stress in ANSYS;  Numerical simulation is divided into two main processes: introducing the initial damage and analyzing the creep effect.
(a) Introducing initial damage Introducing fatigue damage is simulated by coding and calling the USERMAT subroutine. The coding consists of two parts:

•
Calculating the stress increment using the strain increment (∆ε) and updating the stress in ANSYS; • Calculating the Jacobi matrix [P] = ∂∆σ ∂∆ε = [Q] −1 E according to Equations (17) and (20). The state variable USTATEV in the user material subroutine of ANSYS records and saves the damage variable (D) during the operation, which is called and updated in the subsequent operation steps. After coding, the numerical simulation analysis of the damaged materials can be realized by continuously calling USERMAT during the ANSYS runtime.
Essentially, the stress and the damage variable are updated under coupling the stressstrain and damage fields. The specific process is as follows: In the numerical simulation, the time (t) is divided into smaller periods, and for each ∆t, the corresponding equilibrium equation is established to find the increment of stress in this period. Figure 3 shows this process. Stress history plays a crucial role in creep deformation [29]. Generally, the cumulative creep deformation at a specific time can be realized by calculating and storing the stress values of each element in the finite element model at all times. Zhu [30] proposed an improved explicit algorithm with variable step size and an implicit algorithm.

•
As a precondition for ensuring the calculation accuracy and saving computing time, this study adopts a jumping algorithm [28]. (b) Creep simulation The concrete creep model describes the strain growth under constant stress. In the numerical simulation, the time (t) is divided into smaller periods, and for each Δt, the corresponding equilibrium equation is established to find the increment of stress in this period. Figure 3 shows this process. Stress history plays a crucial role in creep deformation [29]. Generally, the cumulative creep deformation at a specific time can be realized by calculating and storing the stress values of each element in the finite element model at all times. Zhu [30] proposed an improved explicit algorithm with variable step size and an implicit algorithm. The process is as follows [30]: The creeping time is divided into several periods,  The process is as follows [30]: The creeping time is divided into several periods, namely ∆t 1 , ∆t 2 , . . . , ∆t n , where ∆t n = t n − t n−1 ; ∆ σ i represents the effective stress increment in step i, and Equation (10) calculates the creep strain at any time.
If ∆t n is small enough, C(t, τ) can be replaced by C(t, t n−0.5 ) at the midpoint. Thus, Equation (10) can be rewritten as: The creep strain increment is also given by: Equation (9) can be rewritten as: where ϕ i (τ) = ( f i + g i τ −p i ). Equation (12) can be written as: where , and the recursive formula for the iterative operator can be defined as: where ω i,1 = ∆ σ 0 ϕ i (t 0 ), and ω i,n is the iterative operator. The creep of the concrete structure can be analyzed by storing ω i,n without recording the stress history, which remarkably saves the storage. For the material constitutive model defined by the user subroutine, it is necessary to calculate the stress increment (∆ σ n ) through the strain increment (∆ε n ) imported by the main program. Without considering the temperature strain, Equation (16) expresses the relationship between the stress increment and the strain increment: where E is the elastic modulus at the midpoint age, and ∆ε e n indicates the increment of the elastic strain.
[Q] and ∆ε e n are expressed by Equations (17) and (18), respectively: where ∆ε e n indicates the increment of the elastic strain: where Combining Equation (16) with Equation (18) yields: where [P] = [Q] −1 E, and [P] is the consistent tangent operator matrix. E is expressed by Equation (20).

Simulation Program
In order to study the creep effect of fatigue-damaged concrete structures after being strengthened, we performed numerical simulations under five different working conditions classified in Table 1.  Figure 4 illustrates the cross-sectional dimensions and reinforcement of the components.

Simulation Program
In order to study the creep effect of fatigue-damaged concrete structures after being strengthened, we performed numerical simulations under five different working conditions classified in Table 1.  Figure 4 illustrates the cross-sectional dimensions and reinforcement of the components. The beam has a rectangular section with a length of 1700 mm and a width of 150 mm; its height depends on the height of the post-cast section, and the concrete cover is 25 mm. The precast section is the existing concrete with dimensions 150 mm × 200 mm. The width of the post-cast part is 150 mm, and its height is set according to the purpose of the study: 0, 50, 100, or 150 mm. The stirrup size of the post-cast part is consistent with that of the precast. Two HRB400 steel bars with a diameter of 8 mm are arranged at the upper part of the section, and two HRB400 steel bars with a diameter of 16 mm are arranged at the lower part of the section under a positive bending moment. The beam has a rectangular section with a length of 1700 mm and a width of 150 mm; its height depends on the height of the post-cast section, and the concrete cover is 25 mm. The precast section is the existing concrete with dimensions 150 mm × 200 mm. The width of the post-cast part is 150 mm, and its height is set according to the purpose of the study: 0, 50, 100, or 150 mm. The stirrup size of the post-cast part is consistent with that of the precast. Two HRB400 steel bars with a diameter of 8 mm are arranged at the upper part of the section, and two HRB400 steel bars with a diameter of 16 mm are arranged at the lower part of the section under a positive bending moment.

Model Development
Elements solid186 and link180 simulate the concrete and the steel bar, respectively. The reinforcement design in this study ensures sufficient anchorage length, so the reinforcements are firmly anchored in concrete to avoid slippage. Therefore, neither the slippage between new and old concretes nor the slippage between reinforcement and concrete is considered in the modeling process. A fixed hinge support is set at one end of the finite element model, and a rolling hinge support is set at the other. A steel block with high stiffness is set at the support to prevent rapid damage development due to the stress concentration at the support.
Further, there is no slippage between the precast beam and the steel block. The beam adopts the mapping grid, and the grid width of the beam and the steel bar is 25 mm. The precast concrete consists of 3264 elements, and the post-cast concrete comprises 1632 elements; the longitudinal bar consists of 64 elements. Figure 5 displays the finite element model of the composite beam. centration at the support.
Further, there is no slippage between the precast beam and the steel block. The beam adopts the mapping grid, and the grid width of the beam and the steel bar is 25 mm. The precast concrete consists of 3264 elements, and the post-cast concrete comprises 1632 elements; the longitudinal bar consists of 64 elements. Figure 5 displays the finite element model of the composite beam. In this study, both the reinforcement and the concrete are under low stress during the loading process due to the slight creep loading. Under such circumstances, the creep of concrete components is affected more by the concrete and less by the reinforcements. Therefore, a bilinear isotropic strengthening model is adopted for the steel bars. The numerical simulation process is divided into two stages, namely fatigue damage introduction and creep loading, realized by two loading steps. The superposition of the post-cast concrete is realized by the life-and-death element method in ANSYS software. The element corresponding to the superimposed layer is activated after the fatigue loading step but before the creep loading step, and its damage variable is set to zero.

Experimental Procedures
Static loading test, fatigue loading test, and creep experiment were conducted to validate the accuracy and applicability of the numerical model proposed herein. The deflection and strain of the specimens were measured and recorded during testing.

Materials
PO 42.5 Portland cement was used as the binder, and the coarse aggregate was crushed stones with a continuous gradation of 5-16 mm; the fine aggregate was river sand. Further, a similar mix proportion presented in Table 2 was used for the precast beam and the lamination layer of concrete. The designed compressive strength of concrete was 30.0 In this study, both the reinforcement and the concrete are under low stress during the loading process due to the slight creep loading. Under such circumstances, the creep of concrete components is affected more by the concrete and less by the reinforcements. Therefore, a bilinear isotropic strengthening model is adopted for the steel bars. The numerical simulation process is divided into two stages, namely fatigue damage introduction and creep loading, realized by two loading steps. The superposition of the post-cast concrete is realized by the life-and-death element method in ANSYS software. The element corresponding to the superimposed layer is activated after the fatigue loading step but before the creep loading step, and its damage variable is set to zero.

Experimental Procedures
Static loading test, fatigue loading test, and creep experiment were conducted to validate the accuracy and applicability of the numerical model proposed herein. The deflection and strain of the specimens were measured and recorded during testing.

Materials
PO 42.5 Portland cement was used as the binder, and the coarse aggregate was crushed stones with a continuous gradation of 5-16 mm; the fine aggregate was river sand. Further, a similar mix proportion presented in Table 2 was used for the precast beam and the lamination layer of concrete. The designed compressive strength of concrete was 30.0 MPa, and the measured 28-day compressive strength of cubic concrete specimens was 30.9 MPa. The measured elastic modulus of the concrete was also 22.9 GPa. The specifications and types of the reinforcement are as described in the numerical simulation model, and Table 3 lists the main mechanical properties of the steel bars.

Section Design
The cross-sectional design is consistent with the numerical model. According to code GB50010-2010 [31], the stirrups of the precast beam should be extended into the lamination layer to guarantee the overall stability of the composite concrete beam. The internal hoops of the cast-in-place concrete are not closed hoops after this setting. On the other hand, cracks generated by fatigue damage are randomly distributed in the concrete [32]. Although the beam undergoes bending failure under static loading, the oblique section shear failure may occur under fatigue loading [33]. Therefore, the shear capacity of the precast beams is guaranteed by adhering carbon fibers to the inclined section to ensure that they do not undergo shear failure during fatigue loading [34][35][36]. The reinforcement schematic is shown in Figure 6.

Section Design
The cross-sectional design is consistent with the numerical model. According to co GB50010-2010 [31], the stirrups of the precast beam should be extended into the lamin tion layer to guarantee the overall stability of the composite concrete beam. The inter hoops of the cast-in-place concrete are not closed hoops after this setting. On the oth hand, cracks generated by fatigue damage are randomly distributed in the concrete [3 Although the beam undergoes bending failure under static loading, the oblique secti shear failure may occur under fatigue loading [33]. Therefore, the shear capacity of precast beams is guaranteed by adhering carbon fibers to the inclined section to ensu that they do not undergo shear failure during fatigue loading [34][35][36]. The reinforcem schematic is shown in Figure 6. The design of carbon fiber reinforcement is according to code GB50367-2013 [37]. C bon-fiber-reinforced plastic (CFRP) is pasted on the shear span of the precast beam, so has a slight influence on the bending bearing capacity of the beam, which can ensure t the damage to concrete primarily concentrates in the pure bending section in the fatig test [38]. The CFRP is U-shaped, is pasted on the whole side and bottom of the beam, a has a thickness of 0.111 mm [39]. The design value of the tensile strength of the CFRP also 1120 N/mm 2 . The CFRP is removed at the end of the fatigue test, and the laminati layer is cast. The purpose of pasting the lamination layer is to improve the load-beari capacity of the section and prevent the shear failure of the beam. The casting size is show in Figure 4, and the materials are consistent with the precast part.

Introducing Fatigue Damage
Firstly, the static loading test was conducted to measure the ultimate bearing capac of the beam so as to determine the fatigue load and the creep load. The static loading t was performed in three phases. In the first phase, each stage was loaded with 3.0 kN befo cracking, and in the second phase, each stage was loaded with 5.0 kN from cracking to The design of carbon fiber reinforcement is according to code GB50367-2013 [37]. Carbon-fiber-reinforced plastic (CFRP) is pasted on the shear span of the precast beam, so it has a slight influence on the bending bearing capacity of the beam, which can ensure that the damage to concrete primarily concentrates in the pure bending section in the fatigue test [38]. The CFRP is U-shaped, is pasted on the whole side and bottom of the beam, and has a thickness of 0.111 mm [39]. The design value of the tensile strength of the CFRP is also 1120 N/mm 2 . The CFRP is removed at the end of the fatigue test, and the lamination layer is cast. The purpose of pasting the lamination layer is to improve the load-bearing capacity of the section and prevent the shear failure of the beam. The casting size is shown in Figure 4, and the materials are consistent with the precast part.

Introducing Fatigue Damage
Firstly, the static loading test was conducted to measure the ultimate bearing capacity of the beam so as to determine the fatigue load and the creep load. The static loading test was performed in three phases. In the first phase, each stage was loaded with 3.0 kN before cracking, and in the second phase, each stage was loaded with 5.0 kN from cracking to the theoretical ultimate load-bearing capacity of the component (72.6 kN). In the third phase, 3.0 kN was loaded until the component failed. The maximum value under the stable state during loading, i.e., 106.0 kN, was taken as the actual ultimate load-bearing capacity of the component.
According to the test results of static loading, the upper load limit for the fatigue loading was considered 65% of the ultimate bearing capacity of the beam, i.e., 68.9 kN, and the stress ratio was taken as 0.1. Fatigue tests were performed on a PMW800 electro-hydraulic pulsation fatigue loading machine with a loading frequency of 4 Hz. During the fatigue loading, linear variable displacement transducers (LVDTs) were employed to measure the midspan deflection of the beams, and strain gauges were used to record the midspan strain of the beams. A DH5922 dynamic data acquisition instrument synchronously collected the deflection and the strain with a sampling frequency of 100 Hz. Figure 7 illustrates the fatigue loading test on the precast beams. measure the midspan deflection of the beams, and strain gauges were used to record th midspan strain of the beams. A DH5922 dynamic data acquisition instrument synchr nously collected the deflection and the strain with a sampling frequency of 100 Hz. Figu 7 illustrates the fatigue loading test on the precast beams.      Figure 8 demonstrates that the midspan deflection of the beam during fatigue loading presents three-phase characteristics. The first phase is the rapid growth of the beam deflection, which accounts for approximately 6.5% of the fatigue life. After that is the stable fatigue evolution phase, in which the beam deflection develops slowly, accounting for 92.5% of the fatigue life. The last phase is the failure stage, where the midspan deflection of the beam develops rapidly into a brittle fracture. Because the fatigue life of the strengthened concrete beams is discrete, and the third phase of fatigue is rapid and difficult to control, three representative working conditions, namely the first phase of the fatigue, the early stage of the second phase, and the middle and late stages of the second phase, are selected as the three initial damage states. Therefore, the loading cycles for fatigue damage degree of 2%, 10%, and 60% are 3000, 15,000, and 90,000, respectively.

Creep Test on Concrete Specimens
The creep coefficient is an essential material parameter for examining the long-term performance of structures. Different stress levels change only the initial strain of concrete, and the creep development has the same pattern, presenting a similar creep coefficient in the linear creep stage of concrete. Standard creep tests were carried out on two specimens, numbered C1 and C2, to obtain the creep coefficient of the material. Moreover, the standard shrinkage test was conducted on specimen SK to assess the shrinkage strain of the concrete. All the above tests were performed at a temperature of 20 ± 2 • C and a relative humidity of 60 ± 5%.
The specimen was loaded by a spring-type creep meter, and the data were recorded by a micrometer and displayed in real time by a pressure sensor. Figure 9 displays the layout of the experiments.
ble fatigue evolution phase, in which the beam deflection develops slowly, accounting for 92.5% of the fatigue life. The last phase is the failure stage, where the midspan deflection of the beam develops rapidly into a brittle fracture. Because the fatigue life of the strengthened concrete beams is discrete, and the third phase of fatigue is rapid and difficult to control, three representative working conditions, namely the first phase of the fatigue, the early stage of the second phase, and the middle and late stages of the second phase, are selected as the three initial damage states. Therefore, the loading cycles for fatigue damage degree of 2%, 10%, and 60% are 3000, 15,000, and 90,000, respectively.

Creep Test on Concrete Specimens
The creep coefficient is an essential material parameter for examining the long-term performance of structures. Different stress levels change only the initial strain of concrete, and the creep development has the same pattern, presenting a similar creep coefficient in the linear creep stage of concrete. Standard creep tests were carried out on two specimens, numbered C1 and C2, to obtain the creep coefficient of the material. Moreover, the standard shrinkage test was conducted on specimen SK to assess the shrinkage strain of the concrete. All the above tests were performed at a temperature of 20 ± 2 °C and a relative humidity of 60 ± 5%.
The specimen was loaded by a spring-type creep meter, and the data were recorded by a micrometer and displayed in real time by a pressure sensor. Figure 9 displays the layout of the experiments. As shown in Figure 9, a jack applied a creep load equal to 30% of the ultimate load, and the pressure was maintained at a changeless level by an adjustable spring. The deformation of the concrete specimen was recorded by dial gauges installed on the opposite surface of the concrete specimens. Therefore, the material parameters t, 1 f , 1 g , 1 p , 1 γ , As shown in Figure 9, a jack applied a creep load equal to 30% of the ultimate load, and the pressure was maintained at a changeless level by an adjustable spring. The deformation of the concrete specimen was recorded by dial gauges installed on the opposite surface of the concrete specimens. Therefore, the material parameters t, f 1 , g 1 , p 1 , γ 1 , f 2 , g 2 , p 2 , and γ 2 in Equation (8) could also be obtained from these tests. Their specific values are listed in Section 5.1.

Creep Test on Strengthened Concrete Beams
After the fatigue damage test was completed, the CFRP of the precast beams was removed, and the post-cast part was poured. After curing, the creep test was carried out employing surcharge creep loading. The composite beams were loaded at three-point with a design value of 5.0 kN. Multiple concrete blocks were used for loading, and the two loading points of the component were uniformly loaded during the loading process. Further, a micrometer was arranged vertically to record the mid-span deflection, and four micrometers were arranged horizontally to record the concrete strain at the corresponding positions. Figure 10 depicts the device for the creeping loading of the composite beams.
ildings 2023, 13, x FOR PEER REVIEW 2 f , 2 g , 2 p , and 2 γ in Equation (8) could also be obtained from these tests values are listed in Section 5.1.

Creep Test on Strengthened Concrete Beams
After the fatigue damage test was completed, the CFRP of the prec removed, and the post-cast part was poured. After curing, the creep test w employing surcharge creep loading. The composite beams were loaded with a design value of 5.0 kN. Multiple concrete blocks were used for loadin loading points of the component were uniformly loaded during the loadin ther, a micrometer was arranged vertically to record the mid-span deflec micrometers were arranged horizontally to record the concrete strain at the positions. Figure 10 depicts the device for the creeping loading of the comp

Material Parameters
The creep deformation of the beams was calculated in the serviceabilit reinforcement and concrete were in a state of low stress. The performance materials markedly affected the creep deformation. The uniaxial compr blocks determined the material parameters in the concrete fatigue damage uniaxial compression creep tests on prisms fitted the creep model.
The stress-strain curve derived from the uniaxial compression test an the cracked concrete obtained from the tensile test were used to determine t age threshold (W), as presented in Table 4. The values of the other param listed in the table following the literature [27,40].

Material Parameters
The creep deformation of the beams was calculated in the serviceability stage, so both reinforcement and concrete were in a state of low stress. The performance of the concrete materials markedly affected the creep deformation. The uniaxial compression tests on blocks determined the material parameters in the concrete fatigue damage model, and the uniaxial compression creep tests on prisms fitted the creep model.
The stress-strain curve derived from the uniaxial compression test and the strain of the cracked concrete obtained from the tensile test were used to determine the initial damage threshold (W), as presented in Table 4. The values of the other parameters are also listed in the table following the literature [27,40]. The creep development curve of the uniaxial compression creep test of the concrete prism is shown in Figure 11. The concrete creep model is fitted according to the standard     Figure 12 delineates the numerical and experimental midspan deflection and midspan deflection increment of the components.

Sustained load (d)
F-0%-100-80-10 in test F-2%-100-80-10 in test F-10%-100-80-10 in test F-60%-100-80-10 in test F-0%-100-80-10 in numerical stimulation F-2%-100-80-10 in numerical stimulation F-10%-100-80-10 in numerical stimulation F-60% -100-80-10in numerical stimulation According to Figure 12, the higher the degree of damage to the precast beam, the more significant the midspan deflection and midspan deflection increment of the composite beam in the initial creep stage. The deflection increment is more suitable for characterizing the creep deflection of the strengthened damaged concrete beams. Comparing the experimental results of each component with the numerical data reveals that both are close in the first 40 days of loading; the experimental values are slightly higher than the numerical ones in the later stages, but the difference in the numerical and experimental deflections is less than 15%. The test results fluctuate to a certain extent due to the influence of envi-       Figure 12 delineates the numerical and experimental midspan deflection and midspan deflection increment of the components.

Sustained load (d)
F-0%-100-80-10 in test F-2%-100-80-10 in test F-10%-100-80-10 in test F-60%-100-80-10 in test F-0%-100-80-10 in numerical stimulation F-2%-100-80-10 in numerical stimulation F-10%-100-80-10 in numerical stimulation F-60% -100-80-10in numerical stimulation According to Figure 12, the higher the degree of damage to the precast beam, the more significant the midspan deflection and midspan deflection increment of the composite beam in the initial creep stage. The deflection increment is more suitable for characterizing the creep deflection of the strengthened damaged concrete beams. Comparing the experimental results of each component with the numerical data reveals that both are close in the first 40 days of loading; the experimental values are slightly higher than the numerical ones in the later stages, but the difference in the numerical and experimental deflections is less than 15%. The test results fluctuate to a certain extent due to the influence of envi- According to Figure 12, the higher the degree of damage to the precast beam, the more significant the midspan deflection and midspan deflection increment of the composite beam in the initial creep stage. The deflection increment is more suitable for characterizing the creep deflection of the strengthened damaged concrete beams. Comparing the experimental results of each component with the numerical data reveals that both are close in the first 40 days of loading; the experimental values are slightly higher than the numerical ones in the later stages, but the difference in the numerical and experimental deflections is less than 15%. The test results fluctuate to a certain extent due to the influence of environmental temperature and humidity on concrete creep. Therefore, the numerical simulation method used to explore the influence of fatigue damage on the creep of the beams, i.e., their midspan deflection, is practical.

Stress on Cross Section
As the concrete was being loaded, its strain was evaluated at a distance of 50, 150, 225, and 275 mm from the bottom of the beam. The concrete strain was measured along the height of the concrete section by pasting dial gauges on its surface. According to the tests by Tang [41] and Fang [13], the strain variation caused by the concrete shrinkage of the old and new concrete of composite components is significant, and it is challenging to completely distinguish the creep strain from the shrinkage strain in tests. Therefore, this study compares only the experimental elastic strain in the midspan section of the composite beams with the simulated data, as shown in Figure 13. ronmental temperature and humidity on concrete creep. Therefore, the numerical simulation method used to explore the influence of fatigue damage on the creep of the beams, i.e., their midspan deflection, is practical.

Stress on Cross Section
As the concrete was being loaded, its strain was evaluated at a distance of 50, 150, 225, and 275 mm from the bottom of the beam. The concrete strain was measured along the height of the concrete section by pasting dial gauges on its surface. According to the tests by Tang [41] and Fang [13], the strain variation caused by the concrete shrinkage of the old and new concrete of composite components is significant, and it is challenging to completely distinguish the creep strain from the shrinkage strain in tests. Therefore, this study compares only the experimental elastic strain in the midspan section of the composite beams with the simulated data, as shown in Figure 13.  Figure 13 shows that the numerical elastic strains calculated without considering the shrinkage effect are closer to the experimental results. The simulation of the creep effect for the fatigue-damaged components is closer to the actual situation, verifying that the developed numerical method is practical.

Fatigue damage
There is a difference in the shrinkage strains between the precast and post-cast parts of the composite beams, leading to additional deflection. This study employs the calculation method of additional deflection proposed by Zhu to correct the deflection [2829]. Figure 14a delineates the variation in the numerical and experimental maximum midspan deflection of the beams with the fatigue cycles. Figure 14b plots the change in the numer-  Figure 13 shows that the numerical elastic strains calculated without considering the shrinkage effect are closer to the experimental results. The simulation of the creep effect for the fatigue-damaged components is closer to the actual situation, verifying that the developed numerical method is practical.

Fatigue damage
There is a difference in the shrinkage strains between the precast and post-cast parts of the composite beams, leading to additional deflection. This study employs the calculation method of additional deflection proposed by Zhu to correct the deflection [28,29]. Figure 14a delineates the variation in the numerical and experimental maximum midspan deflection of the beams with the fatigue cycles. Figure 14b plots the change in the numerical and experimental maximum midspan deflection increment of the beams with the fatigue cycles. Figure 14 demonstrates that the development of the deflection of each precast beam follows similar trends in the fatigue tests and the numerical simulation. The numerical and experimental trends of the maximum midspan deflection increment of the beams are also close.
There is a difference in the shrinkage strains between the precast and post-cast parts of the composite beams, leading to additional deflection. This study employs the calculation method of additional deflection proposed by Zhu to correct the deflection [2829]. Figure 14a delineates the variation in the numerical and experimental maximum midspan deflection of the beams with the fatigue cycles. Figure 14b plots the change in the numerical and experimental maximum midspan deflection increment of the beams with the fatigue cycles. Figure 14 demonstrates that the development of the deflection of each precast beam follows similar trends in the fatigue tests and the numerical simulation. The numerical and experimental trends of the maximum midspan deflection increment of the beams are also close.

Discussion
Many factors affect the shrinkage and creep performance of materials, such as composition [42], temperature, and humidity [43]. The stress on beams and the concrete age at loading are the main factors affecting the creep deformation, but the influence of these factors is not evident. Therefore, it is necessary to discuss the mechanism for the effect of damage on the creep development of beams. Although increasing the section is a standard method for structural reinforcement [44], it is still worth examining the creep of damaged components after being strengthened.

Stress Magnitude
The essential impact of damage on the creep of concrete beams is the degradation of the cross-sectional stiffness, which leads to the difference in the stress distribution of the cross section under the same load. The creep effect depends on the stress state of the beams, so a difference in the cross-sectional stress distribution causes a difference in the creep effect. Hence, this section presents the investigation of the mechanism for the effect of damage on the creep development of the beams through experiments with different sustained loads applied to the undamaged components. The components are described in Section 3.1, and Figure 15 depicts the test results.

Discussion
Many factors affect the shrinkage and creep performance of materials, such as composition [42], temperature, and humidity [43]. The stress on beams and the concrete age at loading are the main factors affecting the creep deformation, but the influence of these factors is not evident. Therefore, it is necessary to discuss the mechanism for the effect of damage on the creep development of beams. Although increasing the section is a standard method for structural reinforcement [44], it is still worth examining the creep of damaged components after being strengthened. The essential impact of damage on the creep of concrete beams is the degradation of the cross-sectional stiffness, which leads to the difference in the stress distribution of the cross section under the same load. The creep effect depends on the stress state of the beams, so a difference in the cross-sectional stress distribution causes a difference in the creep effect. Hence, this section presents the investigation of the mechanism for the effect of damage on the creep development of the beams through experiments with different sustained loads applied to the undamaged components. The components are described in Section 3.1, and Figure 15 Figure 15a shows that the midspan deflection increment of the concrete beams is significantly higher under larger loads because creep correlates closely with the distribution and magnitude of stress on concrete. The undamaged beams are continuously loaded with 5 and 15 kN for 365 days, with midspan deflections of 0.270 and 0.734 mm, respectively. Although the distribution of stress on the beam is not uniform, the stress on the beams under large loads is generally higher than that on the beams under small loads. This is also why the strain of the beams under large loads is generally higher than that of the beams under small loads.
According to a previous study, the creep of concrete components is proportional to the stress when the stress on them is lower than 40% of their strength [45]. The creep coefficient is similar to this case. It can be inferred from Figure 15a that there is a distinct difference in the creep coefficients of the three components: the creep coefficients of the beams under larger loads are relatively smaller. This is because the stresses at various points of the beam constantly change over time and are affected by multiple factors other than the strain effects (see Figure 16). Although the strain is linearly related to the stress under creep, the change in stress does not correlate linearly with the load. As a result, the creep coefficients of the beams under various loads are not identical. The damage causes the redistribution of the stresses on the concrete in the beam cross section. The initial elastic strain and subsequent strains of the fatigue-damaged beams increase to some extent because some part of the concrete withdraws from bearing the load  Although the distribution of stress on the beam is not uniform, the stress on the beams under large loads is generally higher than that on the beams under small loads. This is also why the strain of the beams under large loads is generally higher than that of the beams under small loads.
According to a previous study, the creep of concrete components is proportional to the stress when the stress on them is lower than 40% of their strength [45]. The creep coefficient is similar to this case. It can be inferred from Figure 15a that there is a distinct difference in the creep coefficients of the three components: the creep coefficients of the beams under larger loads are relatively smaller. This is because the stresses at various points of the beam constantly change over time and are affected by multiple factors other than the strain effects (see Figure 16). Although the strain is linearly related to the stress under creep, the change in stress does not correlate linearly with the load. As a result, the creep coefficients of the beams under various loads are not identical.  Figure 15a shows that the midspan deflection increment of the concrete nificantly higher under larger loads because creep correlates closely with th and magnitude of stress on concrete. The undamaged beams are continuously 5 and 15 kN for 365 days, with midspan deflections of 0.270 and 0.734 mm, Although the distribution of stress on the beam is not uniform, the stress under large loads is generally higher than that on the beams under small also why the strain of the beams under large loads is generally higher tha beams under small loads.
According to a previous study, the creep of concrete components is pr the stress when the stress on them is lower than 40% of their strength [45]. efficient is similar to this case. It can be inferred from Figure 15a that ther difference in the creep coefficients of the three components: the creep coeffi beams under larger loads are relatively smaller. This is because the stress points of the beam constantly change over time and are affected by multiple than the strain effects (see Figure 16). Although the strain is linearly related under creep, the change in stress does not correlate linearly with the load. A creep coefficients of the beams under various loads are not identical. The damage causes the redistribution of the stresses on the concrete in th section. The initial elastic strain and subsequent strains of the fatigue-damag crease to some extent because some part of the concrete withdraws from bea The damage causes the redistribution of the stresses on the concrete in the beam cross section. The initial elastic strain and subsequent strains of the fatigue-damaged beams increase to some extent because some part of the concrete withdraws from bearing the load (or the effective area is reduced due to the damage), and the remaining concrete is subjected to increased stress under the same bending moment. This leads to the higher deflection of the damaged beam compared with undamaged beams under creep. It should be emphasized that damage gives rise to different stress states and various rates of stress change. Therefore, it is not reasonable to use undamaged components to anticipate the creep strains of damaged components.

Loading Age
The age of concrete at loading affects its creep properties, and all creep models consider the loading age a crucial parameter. In order to explore the mechanism for the impact of damage on the creep development of concrete beams, here we apply a sustained load to the fatigue-damaged beams of different ages using the component identities presented in Section 3.1. Figure 17 depicts the change in the deflection increment of and the stress on the fatigue-damaged concrete beams of different ages with the loading time. (or the effective area is reduced due to the damage), and the remaining concrete is subjected to increased stress under the same bending moment. This leads to the higher deflection of the damaged beam compared with undamaged beams under creep. It should be emphasized that damage gives rise to different stress states and various rates of stress change. Therefore, it is not reasonable to use undamaged components to anticipate the creep strains of damaged components.

Loading Age
The age of concrete at loading affects its creep properties, and all creep models consider the loading age a crucial parameter. In order to explore the mechanism for the impact of damage on the creep development of concrete beams, here we apply a sustained load to the fatigue-damaged beams of different ages using the component identities presented in Section 3.1. Figure 17 depicts the change in the deflection increment of and the stress on the fatigue-damaged concrete beams of different ages with the loading time. According to Figure 17a, if the additional deflection caused by the difference in the beam shrinkage is not taken into account, the deflection of both damaged and undamaged laminated components under the influence of creep diminishes with age. The deflection of beam A-60%-100-365-10 and beam A-60%-100-3650-10 is 0.252 mm and 0.264 mm at the age of 365 and 3650 days, respectively. It increases by approximately 0.01 mm in 3285 days. For the same material, creep declines with increasing the loading age. This law is also reflected in the creep model used in this paper. A lower creep results in a lower creep strain for the same stress state, leading to a smaller beam deflection. For the same loading age, damage to the beams enlarges their midspan deflection increment, implying that raising stress increases the creep of the beams. Figure 17b demonstrates that the stress on the top of the beam in the post-cast part is lower for the damaged and undamaged beams with longer ages. After 3650 days, the stress on the precast concrete of the undamaged beam is 0.7 MPa, while that on the postcast concrete is -0.7 MPa. The stress on the old concrete is approximately 1.4 MPa greater than that on the new concrete. However, for the precast section, the tensile stress on the beam with a longer age is higher and varies slightly with time. Indeed, the components with a more extended age have lower creep strains, and the smaller strain reduces the change in stress.
Under similar initial stress, the stress on the components with a longer age is always more considerable under the effect of creep. In contrast, at a similar loading age and creep, According to Figure 17a, if the additional deflection caused by the difference in the beam shrinkage is not taken into account, the deflection of both damaged and undamaged laminated components under the influence of creep diminishes with age. The deflection of beam A-60%-100-365-10 and beam A-60%-100-3650-10 is 0.252 mm and 0.264 mm at the age of 365 and 3650 days, respectively. It increases by approximately 0.01 mm in 3285 days. For the same material, creep declines with increasing the loading age. This law is also reflected in the creep model used in this paper. A lower creep results in a lower creep strain for the same stress state, leading to a smaller beam deflection. For the same loading age, damage to the beams enlarges their midspan deflection increment, implying that raising stress increases the creep of the beams. Figure 17b demonstrates that the stress on the top of the beam in the post-cast part is lower for the damaged and undamaged beams with longer ages. After 3650 days, the stress on the precast concrete of the undamaged beam is 0.7 MPa, while that on the post-cast concrete is −0.7 MPa. The stress on the old concrete is approximately 1.4 MPa greater than that on the new concrete. However, for the precast section, the tensile stress on the beam with a longer age is higher and varies slightly with time. Indeed, the components with a more extended age have lower creep strains, and the smaller strain reduces the change in stress.
Under similar initial stress, the stress on the components with a longer age is always more considerable under the effect of creep. In contrast, at a similar loading age and creep, the difference in the strain of the midspan section of the beams directly causes the difference in the stress on them. The strain of the components with a more extended age is lower in the span of the precast part, which is the opposite of the previously mentioned stress variation law. Since the stress on the precast part of damaged beams is coupled with the effect of damage, it is impossible to directly and accurately determine the stress, but the law of undamaged beams can deduce it: the stress increases with the loading age under the creep effect. Nonetheless, the deflection of the beams is relatively small, so the loading age can affect the creep of damaged components and the stress state.
As mentioned above, the loading age is one factor affecting the creep development due to damage to concrete beams. The primary mechanism is that increasing the loading age reduces the degree of creep and the creep coefficient at the same loading time and declines the creep strain in the same stress state, which not only gives rise to the difference in the development of the deflection of the beams under creep but also leads to the difference in the stress state after loading. The effect of creep also depends on the loading age or its creep properties. Although the loading age can be controlled in experimentation or numerical simulation, in practical engineering, the damage is often accompanied by a long loading age, and even the degree of damage is coupled with the loading age. Therefore, the impact of the loading age cannot be neglected in the creep analysis.

Creep Effect of Strengthened Damaged Beams
The discussion in Section 6.1 demonstrates that both the magnitude of the sustained load and the age of the concrete impact the stresses on the strengthened beam section and eventually change the creep effect of the concrete. In fact, when reinforcing the beams by increasing their cross section, the height of the reinforcement significantly impacts the stress state of the component [46] because the height of the reinforcement affects the creep effect of the component. Thus, this section discusses the creep of the damaged beams that were strengthened and subjected to positive and negative bending moments.

Under Positive Bending Moment
This section investigates the effect of reinforcing damaged beams on their creep development from a stress perspective. Thus, different heights of the post-cast layer are set for comparison, as described in Section 3.1. Figure 18 presents the test results. the difference in the strain of the midspan section of the beams directly causes the difference in the stress on them. The strain of the components with a more extended age is lower in the span of the precast part, which is the opposite of the previously mentioned stress variation law. Since the stress on the precast part of damaged beams is coupled with the effect of damage, it is impossible to directly and accurately determine the stress, but the law of undamaged beams can deduce it: the stress increases with the loading age under the creep effect. Nonetheless, the deflection of the beams is relatively small, so the loading age can affect the creep of damaged components and the stress state. As mentioned above, the loading age is one factor affecting the creep development due to damage to concrete beams. The primary mechanism is that increasing the loading age reduces the degree of creep and the creep coefficient at the same loading time and declines the creep strain in the same stress state, which not only gives rise to the difference in the development of the deflection of the beams under creep but also leads to the difference in the stress state after loading. The effect of creep also depends on the loading age or its creep properties. Although the loading age can be controlled in experimentation or numerical simulation, in practical engineering, the damage is often accompanied by a long loading age, and even the degree of damage is coupled with the loading age. Therefore, the impact of the loading age cannot be neglected in the creep analysis.

Creep Effect of Strengthened Damaged Beams
The discussion in Section 6.1 demonstrates that both the magnitude of the sustained load and the age of the concrete impact the stresses on the strengthened beam section and eventually change the creep effect of the concrete. In fact, when reinforcing the beams by increasing their cross section, the height of the reinforcement significantly impacts the stress state of the component [4846] because the height of the reinforcement affects the creep effect of the component. Thus, this section discusses the creep of the damaged beams that were strengthened and subjected to positive and negative bending moments.

Under Positive Bending Moment
This section investigates the effect of reinforcing damaged beams on their creep development from a stress perspective. Thus, different heights of the post-cast layer are set for comparison, as described in Section 3.1. Figure 18 presents the test results.  at a loading period of 365 days. The incremental creep deflection of the beams declines with increasing the height of their reinforcement. The incremental creep deflection of the beams with a strengthened section height of 50, 100, and 150 mm decreases by 0.107, 0.228, and 0.326 mm, respectively, compared with beam H-60%-0-365-10 at a loading period of 365 days. Indeed, increasing the section height makes the damaged components more rigid and weakens the stress state of the section under the same bending moment. Further, the creep strain is slight when the stress is low, reducing the deflection increment of the beams. In addition, the incremental deflection of the strengthened damaged beams is lower than that of the undamaged beams (H-0%-0-365-10). The method of reinforcing damaged beams with laminating post-cast parts offers promising results in terms of controlling their creep deflection; in particular, increasing the section height has a profound effect.
According to Figure 18b, since the stress on the top of the unstrengthened concrete beams coupled with the effect of damage cannot be depicted realistically, the stress on the damaged beams under the effect of creep should be higher than that on the undamaged beam. This figure demonstrates that the stress on the top of the damaged beams declines to some extent during creep after post-cast concrete reinforcement and is lower than that on the undamaged beam (H-0%-0-365-10). Therefore, the reinforcement method of laminating a post-cast concrete part to the damaged beam performs well. The tensile stress on the bottom longitudinal reinforcement of the damaged beams decreases when the height of the reinforcement section increases. However, if the height of the reinforcing post-cast layer reaches 150 mm, the tensile stress on the longitudinal bars of the damaged concrete beams under the effect of creep is similar to that on the undamaged concrete beam (H-0%-0-365-10). Nonetheless, the longitudinal tensile stress on the damaged beam when the height of the reinforcement reaches 150 mm is lower than that on the undamaged beam after 127 days of sustained loading. The method offers an excellent controlling effect for the long-term performance of the strengthened damaged beams.

Under Negative Bending Moment
In engineering practice, components affected by negative bending moments are also typical. Since the damage occurs between the precast and post-cast parts, it is possible to examine whether this damage impacts the creep effect of strengthened components. Section 3.1 defines the component numbers, and Figure 19 depicts the test results.  Figure 19a demonstrates that, under a negative bending moment, fatigue damage increases the midspan deflection of the beams compared with the undamaged beam under the effect of creep. After increasing the cross section, the deflection increment of the damaged components declines to some extent, and the creep deflection increment of damaged component MS-60%-100-365-10 and undamaged component MS-0%-365-10-100 decreases by 21.6% and 32.0%, respectively, after 365 days of sustained loading. Therefore, increasing the cross section better controls the deflection of the damaged components under a negative bending moment.  Figure 19a demonstrates that, under a negative bending moment, fatigue damage increases the midspan deflection of the beams compared with the undamaged beam under the effect of creep. After increasing the cross section, the deflection increment of the damaged components declines to some extent, and the creep deflection increment of damaged component MS-60%-100-365-10 and undamaged component MS-0%-365-10-100 decreases by 21.6% and 32.0%, respectively, after 365 days of sustained loading. Therefore, increasing the cross section better controls the deflection of the damaged components under a negative bending moment.

Conclusions
The study presented in this paper established a simulation for the creep effect of fatigue-damaged concrete components and verified the validity of the model utilizing experimental data so as to analyze the creep of different damaged components under sustained loading. The following are the main conclusions drawn:

•
On the basis of the principle of damage mechanics and strain equivalence, we established a creep structure of damaged concrete and proposed a numerical method for the creep effect of fatigue-damaged concrete components, which is suitable only for fatigue damage components but not for damage components, such as steel corrosion, concrete erosion, etc. The experimental results also validated the proposed model for the creep of damaged concrete components. The model could be used to forecast the creep deformation of undamaged components or damaged components after being strengthened, facilitating structural maintenance and decision-making about reinforcement.

•
The mechanism for the influence of damage on the creep development of concrete beams was analyzed in conjunction with numerical calculations. The damage changed the stiffness of the cross section of the components and varied the stress state of the concrete in the cross section during loading, affecting the creep development of the concrete. Further, the creep of the components developed slowly as the age of the concrete increased.

•
The undamaged beams were continuously loaded with 5 and 15 kN for 365 days, experiencing midspan deflections of 0.270 and 0.734 mm, respectively. When the load was elevated from 5 to 15 kN, the deflection enlarged by roughly 0.46 mm. The deflection of beam A-60%-100-365-10 and A-60%-100-3650-10 was 0.252 mm and 0.264 mm at the age of 365 and 3650 days, respectively. It increased by approximately 0.01 mm in 3285 days. Further, the stress magnitude had a more profound effect on the creep than the loading.

•
The incremental creep deflection of the beams with a strengthened section height of 50, 100, and 150 mm loaded for 365 days decreased by 0.107, 0.228, and 0.326 mm, respectively, compared with the unstrengthened damaged beam. The method of reinforcing damaged beams with laminating post-cast part offered promising results in terms of controlling their creep deflection; in particular, increasing the section height had a profound effect. Moreover, this method excellently controlled the deflection of the damaged components under a negative bending moment.

Prospect
• This work investigated the impact of the creep of cracked concrete under low stress. For the nonlinear creep of damaged components at high stress levels, the initial damage to concrete may differ during the creep process. Thus, the nonlinear creep effect should be further investigated.

•
This study focused on a single sort of damage, i.e., fatigue damage, and its influence on the creep effect of the component. However, in practice, multiple factors affect components' creep. Therefore, further research on creep is required under the combined impact of factors, such as material flaws, load, and the environment. Data Availability Statement: All data, models, and code generated or used during this study are available upon request from the corresponding author.

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