Research Article  Open Access
Shiyong Jiang, Weilai Yao, Jin Chen, Tao Cai, "Finite Element Modeling of FRPStrengthened RC Beam under Sustained Load", Advances in Materials Science and Engineering, vol. 2018, Article ID 7259424, 16 pages, 2018. https://doi.org/10.1155/2018/7259424
Finite Element Modeling of FRPStrengthened RC Beam under Sustained Load
Abstract
External bonding of FRP laminates to the tension soffit of concrete members has become a popular method for flexural strengthening. However, the longterm field performance of FRPstrengthened RC members under service conditions is still a concern, and more work needs to be done. Based on concrete smearedcrack approach, this paper presents a finiteelement (FE) model for predicting longterm behavior of FRPstrengthened RC beam, which considers the timedependent 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. The interfacial edge stresses in adhesive layer were determined through appropriate mesh refinement near the plate end, and their timedependent characteristics were investigated. The 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 longterm behavior of the FRPstrengthened beam and help us better understand the complex changes in the stress state occurring over time.
1. Introduction
Using fiberreinforced 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–3]. Relatively few studies have been performed on the behavior of FRPstrengthened RC structures under longterm services conditions. In fact, under sustained load, concrete exhibits clear creep response, and the viscous property of adhesive is remarkable [4–7]. Related to the kind of fiber used and fiber content, some types of FRP products may also present noticeable viscoelasticity [8–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]. Thus, 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 longterm response of strengthened members is not fully explored, and further research is essential in order to enhance the safety and design of FRPstrengthened members [14, 15].
So far, previous studies have proposed some analytical methods to investigate timedependent characteristic of RC structure strengthened by external reinforcement [11, 16–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 longterm 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 tensionstiffening effect appropriately and accurately. Thus, 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 timedependent 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 timedependent 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]. The 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. Those specific simplifications mentioned above may also be avoided, which makes the FE method more generally applicable. However, to date, studying the timedependent performance by FE approach has not been extensively conducted. One of the likely reasons should be that reasonably developing timedependent constitutive models of materials can be a challenging task.
The major thrust of this paper is to propose a finite element model to investigate the longterm response of FRPstrengthened RC beam under sustained load. The timedependent characteristics of all components including the evolution of concrete property along with time were modeled as userdefined 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.
2. LongTerm Deformation Monitoring
Timedependent deformation of a FRPstrengthened RC beam under sustained load was monitored for 200 d. The sectional dimension was 250 × 400 mm with a 30 mm thickness. The 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.
The test beam was cast by commercial concrete. The 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). The reported material strengths above are measured mean values.



The test beam was simply supported and loaded in four point bending. The 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. The 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.71 kN), 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 longterm 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 longterm working performance is unreliable. Vibrating wire gauges were especially used to obtain both instantaneous and longterm strains in some places because of their reliable performance in long period of monitoring. The 1# wire strain gauge (WG1) was bonded on the top surface of the beam. The 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.
The 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 affected by the change of environmental conditions in the laboratory.
(a)
(b)
3. FE Modeling
Based on ABAQUS 6.12 [25], nonlinear numerical analysis was performed to study the timedependent characteristics of concrete beam strengthened with external FRP reinforcement. The 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 effect of concrete was taken into account and the corresponding calculating method was presented, and also, the method for analyzing the timedependent interfacial stress was developed.
3.1. Modeling of Instantaneous Material Property
In this study, concrete cracking was simulated by using the smearedcracked approach. Unlike the discretecrack approach, there is no need to predefine the crack paths in the smearedcracked approach. In addition, aiming to overcome the meshnonobjectivity problem in the conventional smearedcrack approach, the crack band model was used [27].
The concrete was modeled within the framework of the concretedamaged plasticity model provided by ABAQUS. The stressstrain relationship under uniaxial compression recommended by Saenz [28] (Equation (1)) was applied. Under uniaxial tension, the problem involves tensile cracking. The stresscrack opening displacement curve recommended by Hordijk [29] (Equations (2) and (3)) was used. The parameters in Equations (2) and (3) can be further determined through Equations (4) and (5) based on CEBFIP standard [30]. Figure 4 presents the uniaxial mechanical behavior of concrete.where and (in MPa) are the compressive and tensile stress, respectively. The and are the uniaxial compressive and tensile strength, respectively. The is the compressive strain; is the peak compressive strain (0.002) corresponding to the strength ; is the elastic modulus, according to ACI 318 [31], (in MPa). The (in mm) is the crack opening displacement; (in mm) is the crack opening displacement at the complete release of stress and fracture energy; = 3.0 and = 6.93 are constants determined by test; is the energy required to create a unit area of stressfree crack; (in mm) is size of the maximum aggregate and is assumed as 20 mm if no experimental data are available.
Based on Bazant and Planas [32], mesh sensitivity problem can be effectively overcome through applying the crack band model [33], in which the displacement of crack opening equals to the cracking strain accumulating over the width of the crack band .
According to Rots [34], when meshing concrete with fournode plane stress element with four integration points, the crack band width is taken to be , where is the side length of the element [27]. Thus, the tensile stresscrack opening displacement relationship given by Equation (2) can be transformed into a stressstrain curve needed by FE program through Equation (6).
The damage evolution was also considered to better describe the nonlinearity of concrete. The relationship between the damage factor and plastic strain , which does not take stiffness degradation into account was defined. Damage factor can be determined as Equation (7) according to Tao and Chen [35] for both uniaxial tension and compression. The patterns of damage are illustrated in Figure 4.where is the uniaxial tensile or compressive stress and is uniaxial tensile or compressive strength.
In this study, the response of FRPstrengthened 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, adhesiveCFRP plate, and concretesteel. The wrapped CFRP sheets for avoiding premature debonding were reasonably ignored in the FE model based on the relatively low stress level [20].
3.2. Modeling of TimeDependent Material Property
Under sustained load, concrete presents significant creep behavior, and additional stressindependent shrinkage deformation is also considerable. The different timedependent characteristics of materials induce complicated stress redistribution in structure along with time.
The loading age is written as , and the load sustaining period is from to ( > ). The total time period to is divided by a number of discrete periods: to , to ,…, to , to . This paper presents a quasielastic constitutive model of concrete considering creep and shrinkage based on ageadjusted effective modulus method (Equations (8)–(10)). The model is in recursion form, which comforts codes programming and improves calculating efficiency. The detailed derivation process of Equations (8)–(10) can be found in Jiang et al. [26], which is not presented herein for brevity.
When :where is the stress increment of concrete due to the effect of creep and shrinkage from to , reflecting the redistribution of stress within structure; is the elastic modulus at time ; refers to the shrinkage strain at time ; and . The , , and are material parameters for expanding the creep coefficient as Dirichlet series [26] (Equation (12)). The is the aging coefficient [36] (Equation (13)). The refers to the instantaneous stress state under loading. The is the strain increment of concrete during to , which reasonably considers the effects of stress increment induced by creep and shrinkage.
Due to the progress of hydration, concrete properties including strength and elastic modulus develop with time. The above timedependent elastic modulus of concrete was determined by the timedependent concrete strength (concrete strength at time ) according to ACI 318 [31] (). The time interval ( to , to , …, to , …, to ) was set small (less than 5 d) to obtain accurate enough result. The 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 finished. In this way, the aging effect of concrete was considered. This study applied the evolution law of concrete strength proposed by ACI 209 [37].where is the concrete mean compressive strength at 28 d in MPa and is the age of concrete in days. and 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. Thus, 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. The evolution of concrete property along with time is shown in Figure 5.
The above timedependent property of concrete was programmed as FORTRAN code in userdefined subroutine UEXPAN and USFLD in ABAQUS. Concrete creep coefficient 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 FRPstrengthened system, the performance of adhesive is of crucial importance in providing effective 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 affects the strengthening mechanism, which therefore influences the longterm response of strengthened beam. In this research, to describe the decreased shear relaxation modulus of epoxy with time, Maxwell chain was used [20, 21].where is the shear modulus at infinite time, is number of Maxwell elements, refers to the relaxed shear modulus of the Maxwell element, and is material constant. Following Choi [20], , = 2 d, and . 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 remains constant with time [39, 40], and so, the relaxation modulus in axial direction is derived from that in the shear direction or vice versa: [14, 15].
The viscous property of CFRP was determined as follows [41]:where = timedependent strain; = time after loading, in hour; and = applied stress. The ; = 68950 MPa; = 0.0011; = 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.
3.3. Solution Strategy
In the FE model, only a half of the beam was included by taking advantage of symmetry about midspan plane, which effectively reduces calculating time. The concrete, epoxy adhesive, and CFRP plate were modeled by plane stress element CPS4. The steel reinforcement was modeled using truss elements.
During the FE analyses, stepbystep procedures were applied for calculations. Instantaneous behavior of the beam was calculated first, which was set as the initial condition of the calculation of longterm response afterwards. Then, based on the timedependent properties of materials, numerical simulation of longterm 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. The epoxy layer and the FRP plate were each modeled with two layers of elements in their thicknesses. Nevertheless, when discussing the timedependent interfacial stresses at the bonding edge, mesh refinement is needed, which is particularly mentioned in the later section.
4. Results and Discussion
4.1. Verification of the Recommended FE Model
Under instantaneous loading, comparison of FEpredicted 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. The 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.
 
Note. Theor./exp. in the above table refers to the ratio of result predicted by FE model and test result, which reflects the predicting error. 
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. The timedependent midspan deflection is shown in Figure 6(a). The timedependent 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 longterm 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. The 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 longterm 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. The results (FE1) predicted by the FE model agree well with the experimental results. The predicted timedependent 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. The timedependent strains of top concrete, concrete near compression bars, are also predicted with good accuracy as shown in Figure 6(b). Besides, The 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. The reason is that the maturity of the used commercial concrete has been good enough after 76 d of aging (the time that the longterm test commenced), as shown in Figure 5.
(a)
(b)
The result of FE1 and FE2, respectively, refers to the case with and without considering the evolution of concrete property. The additional deflection () is defined as the delayed deflection occurring in load sustaining period. The total additional deflection () is the delayed deflection after 200 d of sustaining load. The ration between and 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 longterm deformation increases continuously and smoothly. Whereas, from the comparison above, the accuracies of the FE results are favorable enough.
The 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 bendingshear section, propagating vertically first but then obliquely afterwards. As shown in Figure 8, the FE results appropriately simulate this phenomenon. The 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)). The 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.
(a)
(b)
(a)
(b)
(c)
The instantaneous cracks are marked by color grey; the further development of cracks in longterm are marked by color red.
The 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 CFRPstrengthened beams and the other three are from Tan and Saha [43]. Sobuz et al. [42] studied the effects of strengthening scheme and level of sustained load on longterm 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), FBC1L, FBC2L, and FBC3L. Tan and Saha [43] tested the timedependent responses of nine GFRPstrengthened specimens. Three typical specimens are chosen herein, namely, GB340, GB349, and GB359. 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. The aging effect of concrete was considered based on the aforementioned method in Section 3.2. For the specimens from Sobuz’s test [42], the ultimate creep coefficient and ultimate shrinkage strain of concrete were 1.87 and 418 με, which were determined through 3.2.1.1 Compressive creep test and 3.2.1.2 Drying shrinkage test in Sobuz’s paper [42]. For the specimens from Tan’s study [43], the specific 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.
(a)
(b)
4.2. TimeDependent Interfacial Edge Stresses
The bonding interfaces are critical regions that transfer stresses between concrete and FRP. One of the most concerned issues of RC beams reinforced by FRP is the premature debonding at the plate edges or near cracks, where the stresses (both shear and vertical normal) in the adhesive layer are relatively high. In this paper, numerical investigations were performed to study the timedependent edge stresses of adhesive interface through the proposed FE approach. The FE model was established on the basis of the test beam in the experimental study of this paper. The responses of specimen under different levels of applied load F (23.54 kN, 47.08 kN, and 70.62 kN) were simulated. Since the wrapped CFRP sheets near the end of beam were not modeled, their influences 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 fine mesh. Following the recommendation of Teng et al. [24], in the vertical direction, for the adhesive layer, the top part (near concreteadhesive interface) and the bottom part (near adhesiveFRP plate interface) were, respectively, meshed by two layers of elements (the height of each element is 0.1 mm). The 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. The CFRP plate has the same mesh pattern as the adhesive layer. A graded and matching mesh was adopted for the concrete. The 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 findings 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. This aspect cannot be captured by existing approximate analytical solutions [22, 23] developed with the assumption that stresses are invariant across the adhesive thickness.
(a)
(b)
(a)
(b)
AC interface refers to adhesivetoconcrete interface, PA interface refers to platetoadhesive interface, and MA refers to the horizontal section of the adhesive layer obtained by a horizontal cut at midthickness.
In this research, the timedependent characteristic of interfacial edge stress at point A (in AC interface) is especially discussed, which is one of the most important factors causing debonding failures. The changes of interfacial shear and vertical normal stresses with time are shown in Figure 13. In this case, only the creep of adhesive in shear direction is considered (both FRP and concrete are modeled as timeinvariant materials). Figures 13(a) and 13(b) show the viscous flow of epoxy in shear direction leading to significant 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 wellknown long creep period of concrete [21]. For the results shown in Figure 14, the axial 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. The 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 influence the edge stresses.
(a)
(b)
(a)
(b)
The combined effect 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. This corresponds to the findings 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 effect by adhesive creep counteracts the increasing effect by concrete creep on the edge stresses to some extent. Therefore, the edge stresses present lower increases with time, comparing to the case that concrete creep works separately. As the relaxation effect is evident enough, which is larger than the increasing effect 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 significant; thus, there is an ongoing increase for the edge stresses afterwards. Besides, since the relaxation effect of adhesive creep in axial direction is smaller, its counteraction to the increasing effect of concrete creep on edge stress is less evident comparing to that of adhesive creep in shear direction.
(a)
(b)
5. Conclusions
In this study, the challenges associated with the timedependent analysis of RC members strengthened by bonding FRP externally have been discussed, and a finite element model has been developed. The recommended FE model simultaneously considers several aspects of the strengthened structures that the existing analytical methods do not take into account: the timedependent characteristics of all components including the aging effect of concrete, the cracking behavior, and tension stiffening of concrete. These aspects have been considered through the constitutive laws of materials. Userdefined subroutines embedded in the main FE program were used to model the timedependent material properties. Unlike some existing analytical methods, specific simplifications were not included, which may further improve the general applicability of this FE approach.
The 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 reflect the environmental effect, the FEpredicted 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.
The analyses of interfacial stresses can be carried out through the same FE model, which only requires mesh refinement near the plate end, rather than to build a quite different analytical model from the beginning. The analyses show that concrete creep causes a significant increase in the shear and vertical normal stresses in the adhesive layer at the edges. This phenomenon may initiate debonding failure, which is of critical importance for the design of FRPstrengthened beams. Whereas, the viscoelastic adhesive may present a favorable effect on alleviating such problem, since it has been presented that the viscous flow of adhesive can relieve the edge stresses.
The FE model developed in this study may set a theoretical platform to better analyze the timedependent characteristics of FRPstrengthened concrete structures. Further aspects of the structural behavior including effects 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 findings and concepts of this study.
Data Availability
The data in Figures 3 and 6 used to support the findings of this study may be released upon application to the corresponding author (email: 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/(ASCE10900268(2006)10:6(474)). These prior studies (and datasets) are cited at relevant places within the text as references [42, 43].
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors are grateful for the financial support from the Chongqing Science and Technology Committee (Project nos: cstc2015shmszx30006 and cstc2011ab0043).
References
 S. Hong, “Effects of the amount and shape of carbon fiberreinforced polymer strengthening elements on the ductile behavior of reinforced concrete beams,” Mechanics of Composite Materials, vol. 50, no. 4, pp. 427–436, 2014. View at: Publisher Site  Google Scholar
 S. Hong, “Effect of intermediate crack debonding on the flexural strength of CFRPstrengthened RC beams,” Mechanics of Composite Materials, vol. 50, no. 4, pp. 523–536, 2014. View at: Publisher Site  Google Scholar
 M. Rezazadeh, H. Ramezansefat, and J. Barros, “NSM CFRP prestressing techniques with strengthening potential for simultaneously enhancing load capacity and ductility performance,” Journal of Composites for Construction, vol. 20, no. 5, Article ID 04016029, 2016. View at: Publisher Site  Google Scholar
 S. H. Aboubakr, U. F. Kandil, and M. R. Taha, “Creep of epoxyclay nanocomposite adhesive at the FRP interface: a multiscale investigation,” International Journal of Adhesion and Adhesives, vol. 54, pp. 1–12, 2014. View at: Publisher Site  Google Scholar
 Y. Jeong, J. Lee, and W. S. Kim, “Modeling and measurement of sustained loading and temperaturedependent deformation of carbon fiberreinforced polymer bonded to concrete,” Materials, vol. 8, no. 2, pp. 435–450, 2015. View at: Publisher Site  Google Scholar
 Y. Jeong, M. M. Lopez, and C. E. Bakis, “Effects of temperature and sustained loading on the mechanical response of CFRP bonded to concrete,” Construction and Building Materials, vol. 124, pp. 442–452, 2016. View at: Publisher Site  Google Scholar
 I. Costa and J. Barros, “Tensile creep of a structural epoxy adhesive: experimental and analytical characterization,” International Journal of Adhesion and Adhesives, vol. 59, pp. 115–124, 2015. View at: Publisher Site  Google Scholar
 J. Jeon, J. Kim, and A. Muliana, “Modeling timedependent and inelastic response of fiber reinforced polymer composites,” Computational Materials Science, vol. 70, pp. 37–50, 2013. View at: Publisher Site  Google Scholar
 A. Sorzia, “Modelling of creep and stress relaxation test of a polypropylene microfibre by using fractionexponential kernel,” Modelling and Simulation in Engineering, vol. 2016, Article ID 3823047, 7 pages, 2016. View at: Publisher Site  Google Scholar
 J. Tanks, K. Rader, S. Sharp, and T. Sakai, “Accelerated creep and creeprupture testing of transverse unidirectional carbon/epoxy lamina based on the stepped isostress method,” Composite Structures, vol. 159, pp. 455–462, 2017. View at: Publisher Site  Google Scholar
 S. Hong and S. K. Park, “Longterm behavior of fiberreinforcedpolymerplated concrete beams under sustained loading: analytical and experimental study,” Composite Structures, vol. 152, pp. 140–157, 2016. View at: Publisher Site  Google Scholar
 D. Trost, Stress Relaxation of Concrete Structures and Relaxation Tests of Very Old Concrete, Springer, Berlin, Germany, 1982.
 K. W. Nasser and A. M. Neville, “Creep of old concrete at normal and elevated temperature,” ACI Journal Proceedings, vol. 64, no. 2, 1967. View at: Publisher Site  Google Scholar
 E. Hamed and M. A. Bradford, “Creep in concrete beams strengthened with composite materials,” European Journal of MechanicsA/Solids, vol. 29, no. 6, pp. 951–965, 2010. View at: Publisher Site  Google Scholar
 E. Hamed and M. A. Bradford, “Flexural timedependent cracking and postcracking behaviour of FRP strengthened concrete beams,” International Journal of Solids and Structures, vol. 49, no. 13, pp. 1595–1607, 2012. View at: Publisher Site  Google Scholar
 A. R. Marí, E. Oller, J. M. Bairán, and N. Duarte, “Simplified method for the calculation of longterm deflections in FRPstrengthened reinforced concrete beams,” Composites Part B: Engineering, vol. 45, no. 1, pp. 1368–1376, 2013. View at: Publisher Site  Google Scholar
 E. Oller and A. R. Marí, “Longterm bond stresses and debonding failure of FRPstrengthened RC cracked members,” Composites Part B: Engineering, vol. 52, pp. 30–39, 2013. View at: Publisher Site  Google Scholar
 K. Hadjazi, Z. Sereir, and S. Amziane, “Creep response of intermediate flexural cracking behavior of reinforced concrete beam strengthened with an externally bonded FRP plate,” International Journal of Solids and Structures, vol. 9495, pp. 196–205, 2016. View at: Publisher Site  Google Scholar
 A. K. ElSayed, R. A. AlZaid, A. I. AlNegheimish, A. B. Shuraim, and A. M. Alhozaimy, “Longterm behavior of wide shallow RC beams strengthened with externally bonded CFRP plates,” Construction and Building Materials, vol. 51, pp. 473–483, 2014. View at: Publisher Site  Google Scholar
 K. K. Choi, M. M. Reda Taha, M. J. Masia, P. L. Shrive, and N. G. Shrive, “Numerical investigation of creep effects on FRPstrengthened RC beams,” Journal of Composites for Construction, vol. 14, no. 6, pp. 812–822, 2010. View at: Publisher Site  Google Scholar
 P. Meshgin, K. K. Choi, and M. M. R. Taha, “Experimental and analytical investigations of creep of epoxy adhesive at the concrete–FRP interfaces,” International Journal of Adhesion and Adhesives, vol. 29, no. 1, pp. 56–66, 2009. View at: Publisher Site  Google Scholar
 S. Benyoucef, A. Tounsi, K. H. Benrahou, and E. A. Adda Bedia, “Timedependent behavior of RC beams strengthened with externally bonded FRP plates: interfacial stresses analysis,” Mechanics of TimeDependent Materials, vol. 11, no. 34, pp. 231–248, 2007. View at: Publisher Site  Google Scholar
 C. Zhang and J. Wang, “Viscoelastic analysis of FRP strengthened reinforced concrete beams,” Composite Structures, vol. 93, no. 12, pp. 3200–3208, 2011. View at: Publisher Site  Google Scholar
 J. G. Teng, J. W. Zhang, and S. T. Smith, “Interfacial stresses in reinforced concrete beams bonded with a soffit plate: a finite element study,” Construction and Building Materials, vol. 16, no. 1, pp. 1–14, 2002. View at: Publisher Site  Google Scholar
 ABAQUS, ABAQUS 6.12 User’s Manual, 2012.
 S. Jiang, W. Yao, J. Chen, and S. Tao, “Time dependent behavior of FRPstrengthened RC beams subjected to preload: experimental study and finite element modeling,” Composite Structures, vol. 200, pp. 599–613, 2018. View at: Publisher Site  Google Scholar
 G. M. Chen, J. G. Teng, and J. F. Chen, “Finite element modeling of intermediate crack debonding in FRPplated RC beams,” Journal of Composites for Construction, vol. 15, no. 3, pp. 339–353, 2010. View at: Publisher Site  Google Scholar
 L. P. Saenz, “Discussion of equation for the stressstrain curve of concreteby Desayi, P. and Krishan, S,” ACI Journal, vol. 61, no. 9, pp. 1229–1235, 1964. View at: Google Scholar
 D. A. Hordijk, “Local approach to fatigue of concrete,” Civil Engineering and Geosciences, vol. 1, no. 4, pp. 442–448, 1991. View at: Google Scholar
 CEBFIP, CEBFIP Model Code 2010, Ernst & Sohn Publishing House, Berlin, Germany, 2013.
 ACI Committee 318, Building Code Requirements for Structural Concrete and Commentary, American Concrete Institute, Farmington Hills, MI, USA, 2014.
 Z. P. Bazant and J. Planas, Fracture and Size Effect in Concrete and Other Quasibrittle Materials, CRC Press, Boca Raton, FL, USA, 1998.
 Z. P. Bažant and B. H. Oh, “Crack band theory for fracture of concrete,” Matériaux Et Construction, vol. 16, no. 3, pp. 155–177, 1983. View at: Publisher Site  Google Scholar
 J. G. Rots, “Computational modeling of concrete fracture,” Ph.D. thesis, Delft University of Technology, Delft, The Netherlands, 1988. View at: Google Scholar
 Y. Tao and J. F. Chen, “Concrete damage plasticity model for modeling FRPtoconcrete bond behavior,” Journal of Composites for Construction, vol. 19, no. 1, Article ID 04014026, 2015. View at: Publisher Site  Google Scholar
 Z. P. Bazant, “Prediction of creep effects using ageadjusted effective modulus method,” ACI Journal Proceedings, vol. 69, no. 4, 1972. View at: Publisher Site  Google Scholar
 ACI Committee 209, Prediction of Creep, Shrinkage, and Temperature Effects in Concrete Structures, American Concrete Institute, Farmington Hills, MI, USA, 2002.
 C. W. Feng, C.W. Keong, Y.P. Hsueh, Y.Y. Wang, and H.J. Sue, “Modeling of longterm creep behavior of structural epoxy adhesives,” International Journal of Adhesion and Adhesives, vol. 25, no. 5, pp. 427–436, 2005. View at: Publisher Site  Google Scholar
 W. N. Findley, J. S. Lai, K. Onaran, and R. M. Christensen, “Creep and relaxation of nonlinear viscoelastic materials,” Journal of Applied Mechanics, vol. 44, no. 2, p. 364, 1977. View at: Publisher Site  Google Scholar
 Y. F. Wen, R. F. Gibson, and J. L. Sullivan, “Prediction of momentary transverse creep behavior of thermoplastic polymer matrix composites using micromechanical models,” Journal of Composite Materials, vol. 31, no. 21, pp. 2124–2145, 1997. View at: Publisher Site  Google Scholar
 W. N. Findley, “Mechanism and mechanics of creep of plastics,” Society of Plastic Engineering Journal, vol. 16, no. 1, pp. 57–65, 1960. View at: Google Scholar
 H. R. Sobuz, E. Ahmed, N. Mohamed Sutan, N. M. Sadiqul Hasan, M. Alhaz Uddin, and M. Jahir Uddin, “Bending and timedependent responses of RC beams strengthened with bonded carbon fiber composite laminates,” Construction and Building Materials, vol. 29, pp. 597–611, 2012. View at: Publisher Site  Google Scholar
 K. H. Tan and M. K. Saha, “Longterm deflections of reinforced concrete beams externally bonded with FRP system,” Journal of Composites for Construction, vol. 10, no. 6, pp. 474–482, 2006. View at: Publisher Site  Google Scholar
 O. Rabinovich and Y. Frostig, “Closedform highorder analysis of RC beams strengthened with FRP strips,” Journal of Composites for Construction, vol. 4, no. 2, pp. 65–74, 2000. View at: Publisher Site  Google Scholar
 E. Hamed and Z. T. Chang, “Effect of creep on the edge debonding failure of FRP strengthened RC beams–a theoretical and experimental study,” Composites Science and Technology, vol. 74, pp. 186–193, 2013. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Shiyong Jiang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.