Finite Element Modeling of FRP-Strengthened RC Beam under Sustained Load

External bonding of FRP laminates to the tension soffit of concrete members has become a popular method for flexural strengthening. However, the long-term field performance of FRP-strengthened RC members under service conditions is still a concern, and more work needs to be done. Based on concrete smeared-crack approach, this paper presents a finite-element (FE) model for predicting long-term behavior of FRP-strengthened RC beam, which considers the time-dependent properties of all components including the aging effect of concrete. According to the comparison between theoretical predictions and test results, the validity of the FE model is verified. (e interfacial edge stresses in adhesive layer were determined through appropriate mesh refinement near the plate end, and their time-dependent characteristics were investigated. (e results show that creep of concrete and epoxy resin cause significant variations of the edge stresses with time. According to the research in this paper, the FE approach is found to be able to properly simulate the long-term behavior of the FRP-strengthened beam and help us better understand the complex changes in the stress state occurring over time.


Introduction
Using fiber-reinforced polymer (FRP) composites as externally bonded reinforcement is in the forefront of retrofitting technology.Most of previous studies have been rightfully focused on the ultimate limit state behavior [1][2][3].Relatively few studies have been performed on the behavior of FRP-strengthened RC structures under long-term services conditions.In fact, under sustained load, concrete exhibits clear creep response, and the viscous property of adhesive is remarkable [4][5][6][7].Related to the kind of fiber used and fiber content, some types of FRP products may also present noticeable viscoelasticity [8][9][10].As time goes on, structure deformation increases and cracks further develop, which negatively affects the serviceability [11].In addition, it has been shown that the creep deformation of very old concrete can still be substantial [12,13].
us, if a structure is strengthened because the service load will be increased, although the concrete may be aged, creep response caused by the additional load should still receive critical attention.Nevertheless, the long-term response of strengthened members is not fully explored, and further research is essential in order to enhance the safety and design of FRP-strengthened members [14,15].So far, previous studies have proposed some analytical methods to investigate time-dependent characteristic of RC structure strengthened by external reinforcement [11,[16][17][18][19].It can be found that most of them are based on cross section analyses, which satisfy strain compatibility and equilibrium of forces.One of the drawbacks of the analytical method is that it may be difficult to consider some significant characteristics.For instance, the basic mechanism for strengthening is the shear deformation and stress transferring of epoxy adhesive between concrete and external FRP.Normally, epoxy presents evident flow property [20,21], which therefore leads to relaxation of interfacial stress and affects the long-term behavior.Whereas, it is actually not easy to model such shear creep characteristic of adhesive in cross section analyses.Besides, the tension stiffening of concrete needs to be reasonably taken into account when predicting deflections.Nevertheless, for the method of section analysis, it is also a challenging task to model such tension-stiffening effect appropriately and accurately.
us, it can be noted that the above influential factors are usually roughly approximated and even neglected [11,16,19].On the other hand, in order to improve the efficiency of solving process or reduce the difficulty in mathematics, simplifications have been ordinarily introduced in analytical method.Based on axial strain reduction coefficient and curvature reduction coefficient, Hong and Park [11] calculated the time-dependent change of the concrete stress in a simple and empirical way.Marí et al. [16] assumed that the stresses of tension bars and FRP laminate remain constant over time.Although, reasonable theoretical results are still obtained in their studies [11,16], there is no evidence to show those simplifications are suitable for other cases, where the creep responses of materials, the geometry of member, and reinforcement ratio (for steel and FRP) are different.It is therefore implied that the universality of the analytical method applying specific simplification is not favorable enough.Additionally, it should be noted that despite simplifications have been made, large amount of iterative calculations still may be needed to obtain the final results [11].Moreover, the stresses at the adhesive interfaces near the bonding edges are responsible for the debonding failures, and their redistribution with time is an issue worth studying.Nevertheless, when determining and investigating the time-dependent interfacial stresses, cross section analyzing is no longer competent, and a quite different analysis procedure (particularly to interfacial forces) needs to be carried out [22,23].In the related studies, it can be found that although precise and complicated derivations had been conducted, assumption not fully corresponding to the fact was still made in order to make the mathematical process more manageable [22,23].e stresses within the adhesive layer were assumed to be constant through the thickness of the adhesive layer in both studies of Benyoucef et al. [22] and Zhang and Wang [23], which in fact vary strongly across the layer thickness [24].
To eliminate the aforementioned deficiencies, finite element (FE) modeling may be a feasible scheme.Comparing to the analytical method, to study the behavior of FRPstrengthened RC member under sustained load, the FE method may provide a more powerful tool, which calculates numerically and does not operate on the basis of cross section analysis.
ose specific simplifications mentioned above may also be avoided, which makes the FE method more generally applicable.However, to date, studying the time-dependent performance by FE approach has not been extensively conducted.One of the likely reasons should be that reasonably developing time-dependent constitutive models of materials can be a challenging task.
e major thrust of this paper is to propose a finite element model to investigate the long-term response of FRPstrengthened RC beam under sustained load.
e timedependent characteristics of all components including the evolution of concrete property along with time were modeled as user-defined subroutines in ABAQUS 6.12 [25].Besides, load sustaining experiment of a strengthened beam was also conducted.By comparing the experiment and FE results, the validity of the FE model is demonstrated.Also, interfacial stresses were accurately determined by the same FE model, with appropriate mesh refinement near the end of CFRP plate.Additional numerical investigations were carried out to study the changes of edge stresses over time induced by the creep of materials.

Long-Term Deformation Monitoring
Time-dependent deformation of a FRP-strengthened RC beam under sustained load was monitored for 200 d.e sectional dimension was 250 × 400 mm with a 30 mm thickness.e total and effective span was 3300 mm and 2900 mm, respectively.D28 reinforcement (diameter is 28 mm) was used for reinforcement on the compression side, while D14 reinforcement (diameter is 14 mm) was used for reinforcement on the tension side.Under service state, concrete stress is not high, which presents linear creep.In this study, adequate compression bars were used to control instantaneous stress in compressive concrete at lower level simulating the service condition.CFRP plate with length and width of 2600 and 200 mm was bonded on the tension soffit of beam.D8 (diameter is 8 mm) reinforcement at interval of 200 mm was used for the stirrups.To avoid premature debonding near the end of CFRP plate, CFRP sheets were used to wrap around the test beam.Figure 1 provides the detailed information of the test beam.
e test beam was cast by commercial concrete.e 28day cylinder compressive strength of the concrete was 34.0 MPa.Other parameters of concrete are listed in Table 1.Table 2 shows the properties of steel reinforcement.Table 3 gives the properties of epoxy resin (applied as adhesive), CFRP sheet (for avoiding premature debonding), and CFRP plate (for flexural strengthening).
e reported material strengths above are measured mean values.e test beam was simply supported and loaded in four point bending.e sustained load F was 75% of the theoretical ultimate capacity F u (on the basis of cross section analysis) of the same dimensional RC beam without strengthening.e beam was cured for 7 d, then placed indoors for another 69 d before loading.In the test, a dedicated experimental system, made up of vertical reaction frame, screw jack, and steel distribution beam, was designed to apply and sustain the external load, as illustrated in Figure 2. Load was applied to the beam through slowly lifting up the screw jack, meanwhile, its value was monitored by the force sensor at the top of the jack.As the value reached to 0.75 F u (117.71kN), screw jack was locked to prevent lifting or declining, which sustained the load.When the deformation of test beam increased due to creep behavior, the applied load might decrease, then the Jack was unlocked and further lifted up to replenish the force to the designed level (0.75F u ).During the long-term test, the error of applied load was controlled within ±1% of the planned value.
As shown in Figure 1, in the midspan, resistance strain gauges were bonded to the tension, compression steel, and CFRP plate to measure the instantaneous strains only, since their long-term working performance is unreliable.Vibrating wire gauges were especially used to obtain both 2 Advances in Materials Science and Engineering instantaneous and long-term strains in some places because of their reliable performance in long period of monitoring.e 1# wire strain gauge (WG1) was bonded on the top surface of the beam.e 2# wire strain gauge (WG2) tied with the middle compression bar was buried in the specimen when casting concrete.One dial gauge was installed at the bottom midspan point and other two dial gauges installed at pedestals were used mainly to check the symmetrical nature of the loaded beam.e experiment was conducted indoors.Figure 3 shows the recorded temperature and humidity, which presented unavoidable variations with time.It should be noted that the observed creep might be a ected by the change of environmental conditions in the laboratory.

FE Modeling
Based on ABAQUS 6.12 [25], nonlinear numerical analysis was performed to study the time-dependent characteristics of concrete beam strengthened with external FRP reinforcement.e FE model used herein is on the basis of the one developed by Jiang et al. [26]; however, the current FE model in this paper was further updated and is more advanced: the aging e ect of concrete was taken into account and the corresponding calculating method was presented, and also, the method for analyzing the time-dependent interfacial stress was developed.

Modeling of Instantaneous Material
Property.In this study, concrete cracking was simulated by using the smearedcracked approach.Unlike the discrete-crack approach, there is no need to prede ne the crack paths in the smeared-cracked approach.In addition, aiming to overcome the meshnonobjectivity problem in the conventional smeared-crack approach, the crack band model was used [27].
e concrete was modeled within the framework of the concrete-damaged plasticity model provided by  Advances in Materials Science and Engineering ABAQUS.
where σ c and σ t (in MPa) are the compressive and tensile stress, respectively.e f c and f t are the uniaxial compressive and tensile strength, respectively.e ε c is the compressive strain; ε cp is the peak compressive strain (0.002) corresponding to the strength f c ; E 0 is the elastic modulus, according to ACI 318 [31], E 0 4730 f c (in MPa).e w t (in mm) is the crack opening displacement; w cr (in mm) is the crack opening displacement at the complete release of stress and fracture energy; c 1 3.0 and c 2 6.93 are constants determined by test; G F is the energy required to create a unit area of stress-free crack; d a (in mm) is size of the maximum aggregate and is assumed as 20 mm if no experimental data are available.Advances in Materials Science and Engineering Based on Bazant and Planas [32], mesh sensitivity problem can be e ectively overcome through applying the crack band model [33], in which the displacement of crack opening w equals to the cracking strain ε cr accumulating over the width of the crack band h c .
According to Rots [34], when meshing concrete with four-node plane stress element with four integration points, the crack band width is taken to be 2 √ e, where e is the side length of the element [27].
us, the tensile stress-crack opening displacement relationship given by Equation ( 2) can be transformed into a stress-strain curve needed by FE program through Equation ( 6).
e damage evolution was also considered to better describe the nonlinearity of concrete.
e relationship between the damage factor d and plastic strain ε p , which does not take sti ness degradation into account was dened.Damage factor can be determined as Equation ( 7) according to Tao and Chen [35] for both uniaxial tension and compression.e patterns of damage are illustrated in Figure 4.
where σ is the uniaxial tensile or compressive stress and f is uniaxial tensile or compressive strength.In this study, the response of FRP-strengthened RC beam under service state is investigated, therefore, the instantaneous material properties of epoxy adhesive, steel reinforcements, and CFRP plate were modeled as linear elastic.Also, perfect bonds were assumed between adhesiveconcrete, adhesive-CFRP plate, and concrete-steel.e wrapped CFRP sheets for avoiding premature debonding were reasonably ignored in the FE model based on the relatively low stress level [20].

Modeling of Time-Dependent Material Property.
Under sustained load, concrete presents signi cant creep behavior, and additional stress-independent shrinkage deformation is also considerable.e di erent time-dependent characteristics of materials induce complicated stress redistribution in structure along with time.e loading age is written as t 0 , and the load sustaining period is from t 0 to t n (t n > t 0 ).e total time period t 0 to t n is divided by a number of discrete periods: t 0 to t 1 , t 1 to t 2 ,. .., t n−2 to t n−1 , t n−1 to t n . is paper presents a quasielastic constitutive model of concrete considering creep and shrinkage based on age-adjusted e ective modulus method (Equations ( 8)- (10)).e model is in recursion form, which comforts codes programming and improves calculating e ciency.
e detailed derivation process of Equations ( 8)-( 10) can be found in Jiang et al. [26], which is not presented herein for brevity. Δε : equivalent plastic strain (considering damage)

Advances in Materials Science and Engineering
When n 1: where σ * (t i ) is the stress increment of concrete due to the e ect of creep and shrinkage from t i−1 to t i , re ecting the redistribution of stress within structure; E(t i ) is the elastic modulus at time t i ; ε sh (t i ) refers to the shrinkage strain at time t i ; and Δt i t i − t i−1 (i 0, 1, 2, . . ., n − 1, n). e b j , λ j , and m are material parameters for expanding the creep coe cient φ(t, τ) as Dirichlet series [26] (Equation ( 12)).e χ(t n , t n−1 ) is the aging coe cient [36] (Equation ( 13)).e σ 0 refers to the instantaneous stress state under loading.e ε n is the strain increment of concrete during t n−1 to t n , which reasonably considers the e ects of stress increment induced by creep and shrinkage.
Due to the progress of hydration, concrete properties including strength and elastic modulus develop with time.
e above time-dependent elastic modulus of concrete E(t i )(i 0, 1, 2, . . ., n − 2, n − 1, n) was determined by the time-dependent concrete strength f c (t i ) (concrete strength at time t i ) according to ACI 318 [31] (E c 4730 f c ). e time interval (t 0 to t 1 , t 1 to t 2 , . .., t i to t i+1 , . .., t n−2 to t n−1 ) was set small (less than 5 d) to obtain accurate enough result.
e elastic modulus of concrete was assumed to be constant in each time interval, but renewed and updated when calculation of the latest time interval was nished.In this way, the aging e ect of concrete was considered.
is study applied the evolution law of concrete strength proposed by ACI 209 [37].
where f c28 is the concrete mean compressive strength at 28 d in MPa and t i is the age of concrete in days.a 4.0 and b 0.85 are material constant for concrete using ordinary Portland cement.In this study, the test beam was loaded at the age of 76 d after casting concrete.us, the calculated 76 d strength 37.7 MPa and its corresponding elastic modulus 29.0 GPa were used to determine the instantaneous response of the test beam under loading.e evolution of concrete property along with time is shown in Figure 5. e above time-dependent property of concrete was programmed as FORTRAN code in user-de ned subroutine UEXPAN and USFLD in ABAQUS.Concrete creep coe cient and shrinkage strain was calculated using ACI 209 code [37], and during the calculation, the measured average humidity (73%) was used.
In this study, both shear and axial creep of adhesive resin were modeled in the FE approach.For FRP-strengthened system, the performance of adhesive is of crucial importance in providing e ective stress transfer.It has been demonstrated that epoxy adhesive presented evident viscoelastic manner [20,21,38].As a result, the relaxation of interfacial stresses caused by adhesive creep a ects the strengthening mechanism, which therefore in uences the long-term response of strengthened beam.In this research, to describe the decreased shear relaxation modulus of epoxy G a (t) with time, Maxwell chain was used [20,21].
where G u is the shear modulus at in nite time, N G is number of Maxwell elements, G i refers to the relaxed shear modulus of the ith Maxwell element, and τ Gi is material constant.Following Choi [20], N G 1, τ G1 2 d, and G u G 1 /5.In addition, the creep of epoxy in axial direction was also taken into account, since near the midspan, where section moment is large and the instantaneous stress of epoxy resin in axial direction usually reaches a noticeable value, which causes nonignorable creep response.It is assumed that Poisson's ratio of adhesive ] a remains constant with time [39,40], and so, the relaxation modulus in axial direction E a (t) is derived from that in the shear direction or vice versa: [14,15].e viscous property of CFRP was determined as follows [41]: where ε CFRP (t) time-dependent strain; t time after loading, in hour; and σ applied stress.e ε 0 0.379; σ m 68950 MPa; m 0.0011; n 1 0.123 are material constants [41].Similar to concrete, the above mechanical properties of adhesive and CFRP related to time were also programmed in user subroutines UEXPAN and USFLD.6 Advances in Materials Science and Engineering midspan plane, which effectively reduces calculating time.e concrete, epoxy adhesive, and CFRP plate were modeled by plane stress element CPS4.e steel reinforcement was modeled using truss elements.
During the FE analyses, step-by-step procedures were applied for calculations.Instantaneous behavior of the beam was calculated first, which was set as the initial condition of the calculation of long-term response afterwards.en, based on the time-dependent properties of materials, numerical simulation of long-term behavior of the test beam was carried out.In this study, if not otherwise stated, square element with an element size of 10 mm was used to model concrete.Matching element sizes were chosen to represent the steel bars, the epoxy layer, and the FRP plate.e epoxy layer and the FRP plate were each modeled with two layers of elements in their thicknesses.Nevertheless, when discussing the time-dependent interfacial stresses at the bonding edge, mesh refinement is needed, which is particularly mentioned in the later section.

Verification of the Recommended FE Model.
Under instantaneous loading, comparison of FE-predicted results and experimental results of the specimen is shown in Table 4.As shown in Figure 1, the tension steel gauge, compression steel gauge, and FRP plate gauge are labeled as TS, CS, and FP.Several strain gauges were employed at the place of midspan: TS1, TS2, and TS3 were bonded to each of the tension bar; CS1, CS2, and CS3 were bonded to each of the compression bar; and FP1 and FP2 were adjacently bonded to the CFRP plate.In Table 4, the strain of tension bar, compression bar, and CFRP plate is the average value of the measured data by TS1, TS2, TS3; CS1, CS2, CS3; FP1, and FP2, respectively.e strain of the top concrete of test beam was obtained by the 1# wire strain gauge (WG1).Table 4 shows that the test results are reasonably predicted by the FE model.
From Table 4, the compressive strain of top concrete is 277 με under instantaneous loading, indicating the stress level does not exceed 40 percent of strength, thus, applying the linear creep theory for concrete is appropriate in this paper.e time-dependent midspan deflection is shown in Figure 6(a).
e time-dependent strains of top concrete and the concrete near compression bars are presented in Figure 6(b), which were experimentally obtained by the vibrating wire gauges WG1 and WG2, respectively.Distinct additional long-term deformation was observed during load sustaining period.After time passage of 200 d, the deflection, strain of top concrete, and strain of concrete near compression bars were respectively 56.7%, 98.6%, and 59.2% higher than their initial value under instantaneous loading.
e increasing of deformation corresponds to the typical law of RC member under sustained load: quick in early period and slow in later time.From Figure 6(a), it is shown that after sustaining load for 78 d (39% of the total time period), most part of additional long-term deflection (87.2%) is finished.Also, as shown in Figure 6(b), there are no clear increases for the strains of concrete after 78 d of sustaining load.e results (FE1) predicted by the FE model agree well with the experimental results.e predicted time-dependent deflection at the 50, 99, 155, and 200 d after sustaining load is 4.97, 5.24, 5.42, and 5.52 mm, respectively, which is close to the corresponding test values of 5.00, 5.43, 5.52, and 5.59 mm, with differences of only 0.6%, 3.5%, 1.8%, and 1.3%, respectively.e time-dependent strains of top concrete, concrete near compression bars, are also predicted with good accuracy as shown in Figure 6(b).Besides, e predicted results (FE2) without considering the further evolution of concrete property are also displayed in Figure 6.As expected, neglecting the evolution of concrete property leads to overestimated results, however, with slight effect only.e reason is that the maturity of the used commercial concrete has been good enough after 76 d of aging (the time that the long-term test commenced), as shown in Figure 5.
e result of FE1 and FE2, respectively, refers to the case with and without considering the evolution of concrete property.e additional deflection (Δf) is defined as the delayed deflection occurring in load sustaining period.e total additional deflection (Δf total ) is the delayed deflection after 200 d of sustaining load.e ration between Δf and Δf total represents the creep progress of the test beam.
As shown in Figure 6, in this study, the test results of both strain and deflection do not present continuing growth along with time, which is likely to be caused by the fluctuations of environmental humidity and temperature [11].However, effects of climate changing were currently not considered during FE modeling; therefore, the predicted long-term deformation increases continuously and smoothly.Whereas, from the comparison above, the accuracies of the FE results are favorable enough.
e proposed FE model rationally predicts the cracking behavior of concrete.Figure 7(a) shows the experimentally observed cracking pattern of test beam under instantaneous loading.In this study, because adequate stirrups were deployed to prevent the occurrence of shear failure, in the earlier stage of loading, the observed cracks were flexural cracks in vertical direction mainly.With the increasing of load, new cracks gradually emerged in the bending-shear section, propagating vertically first but then obliquely afterwards.As shown in Figure 8, the FE results appropriately simulate this phenomenon.e effect of tension stiffening, which causes the concrete between adjacent cracks to suffer tensile stress, is also simulated by the FE model, as presented in Figure 8(c).During the test, further development of instantaneous cracks during load sustaining period was clearly observed (Figure 7(b)).e combined effect of materials' creep behaviors and decreasing of section stiffness induced by further formation of cracks causes the additional longterm deformation of test beam.
e instantaneous cracks are marked by color grey; the further development of cracks in long-term are marked by color red.
e FE model in this paper has been further validated against the test results in other literatures [42,43].Totally, seven specimens are selected herein: four of them are from Sobuz et al. [42] including one control beam (not strengthened) and three CFRP-strengthened beams and the other three are from Tan and Saha [43].Sobuz et al. [42] Advances in Materials Science and Engineering studied the e ects of strengthening scheme and level of sustained load on long-term response of RC beam strengthened by CFRP.Four representative specimens are selected herein to demonstrate the accuracy of the proposed FE approach, namely, CBC (not strengthened, the control beam), FBC-1L, FBC-2L, and FBC-3L.Tan and Saha [43] tested the time-dependent responses of nine GFRPstrengthened specimens.ree typical specimens are chosen herein, namely, GB3-40, GB3-49, and GB3-59.Table 5 presents the detailed information of these specimens.
According to Sobuz et al. [42] and Tan and Saha [43], the viscoelasticity of FRP plate was ignored.e aging e ect of concrete was considered based on the aforementioned method in Section 3.2.For the specimens from Sobuz's test [42], the ultimate creep coe cient and ultimate shrinkage strain of concrete were 1.87 and 418 με, which were determined through 3.2.1.1Compressive creep test and 3.2.1.2Drying shrinkage test in Sobuz's paper [42].For the specimens from Tan's study [43], the speci c laws of concrete creep and shrinkage have been given in their paper [43], which were, therefore, directly used in the FE simulations.From Figure 9, the proposed FE approach well predicts the test results.

Pure bending section Wrapped CFRP sheets
The top face of test beam   in uences on the localized interfacial stresses were not discussed currently.
As shown in Figure 10, due to the existence of singular points (Point A and B), in order to obtain accurate stresses, the region near the plate end needs a particularly ne mesh.Following the recommendation of Teng et al. [24], in the vertical direction, for the adhesive layer, the top part (near concrete-adhesive interface) and the bottom part (near adhesive-FRP plate interface) were, respectively, meshed by two layers of elements (the height of each element is 0.1 mm).
e rest part of the adhesive layer was evenly meshed by other three layers of elements.In the horizontal direction, a graded mesh was used starting with an aspect ratio of 1 for the minimum height (0.1 mm) elements.
e CFRP plate has the same mesh pattern as the adhesive layer.A graded and matching mesh was adopted for the concrete.e stress distribution near the plate end under instantaneous loading is given in Figures 11 and 12 (applied load F 70.62 kN), whose characteristics correspond to the ndings of previous researches [24,44].It shows that high vertical normal (peeling) and shear stresses exist near the bonding end, which are responsible for the debonding failures widely reported in tests.Also, the normal and shear stresses vary strongly across the thickness of the adhesive layer in the vicinity of the end of the plate.is aspect cannot be captured by existing approximate analytical solutions [22,23] developed with the assumption that stresses are invariant across the adhesive thickness.
AC interface refers to adhesive-to-concrete interface, PA interface refers to plate-to-adhesive interface, and MA refers   13.In this case, only the creep of adhesive in shear direction is considered (both FRP and concrete are modeled as time-invariant materials).Figures 13(a) and 13(b) show the viscous ow of epoxy in shear direction leading to signi cant relaxations for both shear and vertical normal stresses at the plate end.With the increase of sustained load, the level of instantaneous stress in epoxy layer is promoted, which causes more evident creep response and stress relaxation.In Figure 13, the relaxation rate of stress is fast in the earlier stage, but almost slows down to zero afterwards.Such phenomenon is believed to be induced by the shear creep law of adhesive adopted herein, which presents that most of creep deformation develops within relatively short time period compared to the well-known long creep period of concrete [21].For the results shown in Figure 14, the axial Advances in Materials Science and Engineering creep of epoxy is only considered to separately investigate its effect on the edge stress.It is shown that like the shear creep, the axial creep of adhesive also causes relaxations of the edge stresses, however, with a much more inferior effect.e likely reason is that near the plate end, the section moment is small, and there is low axial stress in the adhesive layer;  hence, the axial creep is not obvious, which does not clearly in uence the edge stresses.
e combined e ect of concrete creep and adhesive creep was also investigated, as illustrated in Figure 15.It is found that the creep of concrete individually leads to increases of both interfacial shear and vertical normal stresses along with time, which may eventually lead to premature debonding.is corresponds to the ndings in previous studies [14,15,45], which demonstrates the rationality of the suggested FE approach.However, as the creep of concrete and adhesive occurs at the same time, the relaxation e ect by adhesive creep counteracts the increasing e ect by concrete creep on the edge stresses to some extent.erefore, the edge stresses present lower increases with time, comparing to the case that concrete creep works separately.As the relaxation e ect is evident enough, which is larger than the increasing e ect of concrete creep, the stresses even show decreases in the earlier stage (curve A, B, and C in Figure 15).Following that, the creep of adhesive is almost completed, while the creep of concrete is still signi cant; thus, there is an ongoing increase for the edge stresses afterwards.Besides, since the relaxation e ect of adhesive creep in axial direction is smaller, its counteraction to the increasing e ect of concrete creep on edge stress is less evident comparing to that of adhesive creep in shear direction.Advances in Materials Science and Engineering

Conclusions
In this study, the challenges associated with the timedependent analysis of RC members strengthened by bonding FRP externally have been discussed, and a nite element model has been developed.e recommended FE model simultaneously considers several aspects of the strengthened structures that the existing analytical methods do not take into account: the time-dependent characteristics of all components including the aging e ect of concrete, the cracking behavior, and tension sti ening of concrete.ese aspects have been considered through the constitutive laws of materials.User-de ned subroutines embedded in the main FE program were used to model the time-dependent material properties.Unlike some existing analytical methods, speci c simpli cations were not included, which may further improve the general applicability of this FE approach.e capability and accuracy of the proposed FE model were demonstrated through comparisons of its predictions with the results of test performed by the authors.Although, detailed changes of environmental conditions have not been considered and the average humidity during test was used to generally re ect the environmental e ect, the FE-predicted results are in qualitative agreement with the experimental observations.Besides, test data in other research studies were also used to further verify the FE approach.
e analyses of interfacial stresses can be carried out through the same FE model, which only requires mesh re nement near the plate end, rather than to build a quite di erent analytical model from the beginning.e analyses show that concrete creep causes a signi cant increase in the shear and vertical normal stresses in the adhesive layer at the edges.
is phenomenon may initiate debonding failure, which is of critical importance for the design of FRPstrengthened beams.Whereas, the viscoelastic adhesive may present a favorable e ect on alleviating such problem, since it has been presented that the viscous ow of adhesive can relieve the edge stresses.
e FE model developed in this study may set a theoretical platform to better analyze the time-dependent characteristics of FRP-strengthened concrete structures.Further aspects of the structural behavior including e ects of environment changes (e.g., variations of humidity and temperature), development of debonded regions, and creep rupture need to be further researched, which can be built on the ndings and concepts of this study.
Data Availability e data in Figures 3 and 6 used to support the ndings of this study may be released upon application to the corresponding author (e-mail: jiangshiy@163.com).Previously reported data in Figure 9 were used to support this study and are available at doi: 10.1016/j.conbuildmat.2011.11.006 and doi: 10.1061/(ASCE1090-0268(2006)10:6(474)).ese prior studies (and datasets) are cited at relevant places within the text as references [42,43].

Conflicts of Interest
e authors declare that they have no con icts of interest.

Figure 6 :
Figure 6: Test and FE results of the specimen.(a) De ection at the midspan.(b) Concrete strains.

Figure 10 Figure 11 :
Figure 10: e element re nement near the plate end.

Figure 14 :
Figure 14: e separate e ect of axial creep of adhesive on the edge stresses.(a) Variation of vertical normal stress with time.(b) Variation of shear stress with time.

Figure 15 :
Figure 15: e combined e ect of adhesive creep and concrete creep on the edge stresses.(a) Variation of vertical normal stress with time.(b) Variation of shear stress with time.

Table 2 :
Materials' properties of steel reinforcement.

Table 4 :
Immediate response at midspan of the test beam.Note.eor./exp. in the above table refers to the ratio of result predicted by FE model and test result, which re ects the predicting error.

Table 5 :
Properties of selected specimens for numerical simulation.
Figure 13: e separate e ect of shear creep of adhesive on the edge stresses.(a) Variation of vertical normal stress with time.(b) Variation of shear stress with time.