#### Abstract

Fiber-reinforced polymers (FRP) are being increasingly used for the repair and strengthening of deteriorated or unsafe concrete structures, including structurally deficient concrete highway bridges. The behavior of FRP strengthened concrete bridge girders, including failure modes, failure loads, and deflections, can be determined using an analytical finite element modeling approach, as outlined in this paper. The differences in flexural versus shear FRP strengthening and comparison with available design guidelines are also beneficial to design professionals. In this paper, a common AASHTO type prestressed concrete bridge girder with FRP wrapping was analyzed using the ANSYS FEM software and the ACI analytical approach. Both flexural and shear FRP applications, including vertical and inclined shear strengthening, were examined. Results showed that FRP wrapping can significantly benefit concrete bridge girders in terms of flexure/shear capacity increase, deflection reduction, and crack control. The FRP strength was underutilized in the section selected herein, which could be addressed through decrease of the amount of FRP and prestressing steel used, thereby increasing the section ductility. The ACI approach produced comparable results to the FEM and can be effectively and conveniently used in design.

#### 1. Introduction

Fiber-reinforced polymers (FRP) are recent innovations in structural engineering, as compared to concrete, steel, and wood. These materials have certain advantages over the traditional materials, such as high stiffness-to-weight and strength-to-weight ratios, corrosion resistance, and constructability. Civil engineering applications of FRP sheets include rehabilitation or restoration of the strength of a deteriorated structural member, retrofitting or strengthening a sound structural member to resist increased loads, and correction of design or construction errors. FRP wrapping is increasingly being used with structural members made of reinforced concrete, prestressed concrete, and masonry.

One area in which FRP sheets are being used quite frequently in recent years is the strengthening of structurally deficient or damaged concrete bridges. As is widely known, a significant percentage of the bridges in the US is structurally deficient. Deficiency in bridges can be caused by design flaws, deterioration due to environmental impact, increase in service loads, and accidental vehicular impacts [1]. Traditional techniques to strengthen structural members include externally bonded steel plates, steel or concrete jackets, and external posttensioning (ACI) [2]. Labor and equipment costs to install FRP systems are often lower than traditional techniques, and installation is easier in areas with limited access. FRP systems also provide better aesthetics in many cases.

FRP wrapping of concrete bridge girders and columns can improve flexural, shear, corrosion, seismic, and impact resistance (Figure 1). A state-of-the-art report by ACI Committee 440 [2] provides guidance for the selection, design, and installation of FRP wrapping systems for externally strengthening concrete structures, based on experimental research, analytical work, and field applications. The ACI report outlines design procedure for flexure, shear, and axial force and combined axial and bending forces, which are considered to be conservative. ACI 440 also mentions areas that require additional research, including theoretical modeling of FRP strengthened structures. Besides experimental and field tests, finite element models (FEM) of FRP strengthened structural members are an alternate and economic avenue to determine their behavior.

The objective of this study [3] was to model a prestressed concrete AASHTO type IV bridge girder using the finite element software ANSYS [4], to analyze flexural and shear strengthening with FRP, and to compare the results obtained with ACI 440 provisions. AASHTO type IV bridge girders are popular for moderate spans (24 to 36 m) and are widely used in the USA (Figure 2). This study provides important information on the FEM procedure for FRP strengthened prestressed concrete girders. The results obtained are very helpful to designers and researchers in understanding practical and cost-effective design procedure for flexure and shear strengthening of prestressed concrete bridge girders with FRP wrapping.

Experimental and analytical studies of Carbon FRP (CFRP) retrofitted prestressed concrete girders showed that the use of CFRP can result in increased moment and shear resisting capacities. Di Ludovico et al. [5] investigated five prestressed concrete I-shaped girders, which were designed according to ANAS (Italian Transportation Institute). One of them was undamaged, two of them were predamaged, and two of them were CFRP retrofitted. Four-point loading using two hydraulic jacks was applied up to theoretical yielding. It was suggested that fiber debonding between the girder and CFRP may take place in the undamaged girder, which has to be prevented to restore full flexural capacity. It was found that the U-shaped CFRP laminates experienced 1% strain before debonding. Cementitious mortar had to be placed between concrete and the repair zone to ensure perfect bonding. Both stiffness and flexural moment capacity of girders increased due to CFRP wrapping. Cerullo et al. [6] conducted an investigation on a real decommissioned bridge girder repaired with externally bonded CFRP. Repair of one AASHTO type III concrete girder was carried out from a damaged bridge. Before CFRP application, cracks on the damaged girder were mapped and flexural behavior was examined by elastic load test. The shear strength deficiency of the girder was found by investigating the crack pattern and spalling of concrete. Horizontal cracks, caused by flexural loading, could successfully be repaired using CFRP. The Canadian Standards Association (CSA) code was found to be conservative for design. The CFRP retrofitting took only three days, which was found to be very convenient. Petty et al. [7] tested eight simply supported I-shaped 42-year-old prestressed concrete bridge girders retrofitted with CFRP for ultimate shear capacity. A single concentrated load was applied near the critical shear load location. Only 54% of the tensile capacity of CFRP could be utilized due to the variation in cross section of the girder. The strut-and-tie model was found to be effective in predicting the shear capacity of the girders. The AASHTO LRFD method was found to be more accurate than the ACI approach for shear capacity prediction of bridge girders. The girder shear capacity increased by 36% due to the vertical U-shaped strips.

Several design guides, standards, and manufacture’s guidelines are available worldwide for the design and analysis of FRP strengthening systems for concrete structures. Some of these provisions are based on theoretical models, while others are based on experimental work [8]. In the USA, the primary design/analysis source is the guidelines published by Committee 440 of the American Concrete Institute. In ACI 440, the design recommendations are based on limit state method and strength/serviceability requirements. Additional load factors are applied to the contribution of the FRP reinforcement. These factors were determined based on statistical evaluation of variability in mechanical properties, predicted versus full-scale test results, and field evaluations. The following failure modes are considered: () crushing of concrete before yielding of the reinforcing steel; () yielding of tension steel followed by rupture of the FRP laminate; () yielding of the tension steel followed by concrete crushing; () shear/tension delamination of the concrete cover; and () debonding of the FRP from the concrete substrate [8]. Although Mode 2 is preferred in order to fully utilize the strengths of both the prestressing steel and the FRP laminate, laboratory testing to date has only shown this mode as a possibility for beams with proper detailing. In most cases, Mode 4 has controlled the failures. Additional information on this failure mode is available in the literature [9–12].

#### 2. Materials and Methods

ANSYS Parametric Design Language (APDL) 14.5 was used herein to model a typical AASHTO type IV I-girder. ANSYS is capable of predicting the nonlinear behavior of FRP strengthened prestressed girders. A simply supported typical interior girder with a span of 24 m was considered.

Three types of ANSYS elements were used in the model. The Solid 65 element was used to create 3D models of concrete. This element is capable of simulating concrete cracking in tension and crushing in compression. It has eight nodes and three degrees of freedom at each node: translations in the nodal , , and directions. This element is also capable of simulating plastic deformation and creep. Link 180 element (uniaxial tension-compression) was used to model the prestressing steel and shear reinforcements. It has two nodes with three degrees of freedom at each node: translations in the nodal , , and directions. This element is capable of rotation, large deflection, and large strain. FRP wrapping was modeled using the Shell 41 element. This element has four nodes and three degrees of freedom at each node: translations in the nodal , , and directions. Shell 41 is a 3D element having membrane stiffness but no bending stiffness and has variable thickness, stress stiffening, and large deflection options.

Four sets of FEM real constants (constitutive properties) were used. Set 1 for the Solid 65 elements had the smeared reinforcement material numbers, volume ratio, and orientation angle entered as zero, as reinforcement was modeled with Link 180 elements. Set 2 was used for the Link 180 element to represent longitudinal prestressing steel. Set 3 was used for Link 180 element to represent the stirrups, and Set 4 was used for Shell 41 elements to model the epoxy adhesive. Set 5 was used for the FRP modeling with Shell 41 elements. Other parameters, such as element -axis rotation, elastic foundation stiffness, and added mass, were entered as zero, as they were not applicable for this model.

The material properties were defined by five models for concrete, mild steel rebar, prestressed steel, epoxy, and FRP. Parameters needed to define the material models are shown in Table 1.

Material Model #1 was defined for concrete. The Solid 65 element requires linear isotropic and multilinear isotropic concrete material properties. The multilinear isotropic material uses the Von Mises failure criterion along with the Willam and Warnke [13] model to define the concrete failure. The concrete compressive uniaxial stress-strain relationship was obtained using the following equations from MacGregor [14]:where is modulus of elasticity; is stress at any strain *ε*; and *ε*° is strain at the ultimate compressive strength, .

The multilinear relationship requires the first point of the curve to be defined by the user and satisfy Hooke’s Law. The curve is used to help with convergence of the nonlinear solution algorithm. Figure 3(a) shows the concrete stress-strain relationship, based on MacGregor [14]. Linearity changes at 0.3 and at a point defined by and *ε*° = 0.003, indicating traditional crushing strain for unconfined concrete. The intermediate points were calculated using (1). The shear transfer coefficients for open and closed cracks were determined using Wolanski [15] as a basis. The uniaxial cracking stress was based on the modulus of rupture, , as shown in the following equation [16]:

**(a) Concrete**

**(b) Prestressing steel**

**(c) FRP**

The other variables in the concrete model were left to default.

Vertical or inclined U-stirrups of 9.5 mm diameter were spaced at 457 mm center throughout the span. Material Model #2 for shear steel was bilinear isotropic, requiring the yield stress and the tangent modulus of the steel to be defined (Table 1). The prestressing steel consisted of 28 tendons that were 12 mm in diameter. The tendons were arranged in two layers in the bottom flange of the girder, as shown in Figure 4(a). They were lumped and modeled as concentrated at the center of the 28 tendons, connecting the nodes at 127 mm from the bottom surface of the girder. Material Model #3 for prestressing steel, as shown in Figure 3(b), was multilinear isotropic following the Von Mises failure criteria, using the following [15]:where is strain in prestressing steel and is stress in prestressing steel.

**(a) Actual**

**(b) Modified**

**(c) Meshed**

The assumed thickness of epoxy adhesive was 0.51 mm. Material Model #4 for the epoxy was linear isotropic (Table 1). The one-layer CFRP sheet thickness was assumed as 1 mm and was represented with Model #5 with linear orthotropic properties, as shown in Figure 3(c). The assumed values for the CFRP and epoxy are common and can be found at manufacturer’s websites.

AASHTO type IV girder section is not quite adaptable for FEM meshing and reinforcement modeling. Therefore, some conversion was made with the converted shape having the same area and height as the original girder and almost the same moment of inertia. Figure 4 shows the actual, converted, and FEM meshed girder sections, and Table 2 shows the comparison of the cross-sectional properties. The area was meshed with 76 mm mesh size, based on the nodes needed to model reinforcement and the aspect ratio for the elements. The meshed area with 140 nodes was extruded for the 24 m span length, as shown in Figure 5(a). The volume was divided into 53 elements along the longitudinal direction to allow for modeling of the stirrups. The prestressing steel model is shown in Figure 5(b). A two-layered section was used for the epoxy and FRP models (Figures 5(c) and 5(d)).

**(a) Concrete**

**(b) Stirrup and prestressing steel**

**(c) Flexural FRP**

**(d) Shear FRP**

Five different girders with FRP and load configurations were investigated herein through separate modeling. Configurations 1 and 2 (controls) were girders without any FRP application and flexural or shear failure mode, respectively. Configuration 3 assumed one layer of FRP for flexural strengthening throughout the bottom of the girder bottom flange. Configuration 4 for shear strengthening involved vertical U-wrap FRP throughout the girder span, while Configuration 5 had U-wrap FRP as 45-degree inclined strips throughout the span. The shear cracks in beams are developed by diagonal tension. Therefore, it is expected that the 45-degree U-wraps will be more efficient in arresting the shear cracks as the orientation of the fiber is perpendicular to the crack in this configuration.

For the FRP strengthened and control girders, the flexural and shear designs were based on the AASHTO LRFD Bridge Specifications [17] and ACI 440 provisions [2].

##### 2.1. Nonlinear Analysis

The 24 m simple span was simulated with a three-point loading setup. For Configurations 1 and 3, the point load was applied at mid-span for maximum bending moment effect, as shown in Figure 6(a). For Configurations 2, 4, and 5, the point load was applied adjacent to the support for maximum shear force effect (Figure 6(b)). A static analysis approach was utilized herein, with simple support displacement constraints provided at the nodes at the two ends. In the first load step, only initial prestrain was applied. Prestressing was defined using ANSYS Command Window [4]. The initial strain was determined from the effective prestress and the modulus of elasticity of prestressing steel. The applied prestress was 690 MPa, selected herein to facilitate convergence. The camber of the model girder due to prestress is shown in Figure 7. In the second load step, the self-weight of the girder was applied. The FEM was set up to examine three different girder behaviors: initial cracking, yielding of the prestressing steel, and flexural/shear failure of the girder. The FRP strengthened beam design was developed in order to avoid the Mode 4 [2] cover delamination failure in order to concentrate on fully achieving the tensile strength of the FRP composite fabric.

**(a) Flexural load**

**(b) Shear load**

##### 2.2. Analytical Calculations

Hand calculation was used herein to determine the expected theoretical capacities of the girders. The analysis is based on the steel and concrete strengths and the cross-sectional properties and steel layout (both prestressed and regular). Provisions from the AASHTO LRFD were utilized for the strength analysis of the control girder [17, Chapter 5]. AASHTO LRFD Specifications currently do not include provisions for the strength determination of FRP strengthened bridge members. Therefore, provisions from the ACI 440 guidelines were used to calculate flexural strengths for the FRP-reinforced girders [2, Chapter 10]. No load factors or resistance factors were applied. Based on the theoretical flexural capacities, equivalent three-point concentrated loads (Figure 6) were calculated from statics. The analytical and FEM load capacities for each girder are compared in Section 3.

#### 3. Results and Discussion

FEM results are presented herein in the form of the values of the equivalent concentrated applied load corresponding to the various loading stages.

##### 3.1. Unstrengthened (Control) Girder for Flexure (Configuration 1)

The flexural load-deflection curve for the control girder showed seven distinct points: effective prestress, addition of self-weight, zero deflection, decompression, initial cracking, steel yielding, and failure (Figure 8). Hand calculations were used herein to validate the FEM results herein, as shown in Table 3. Initial prestress application stage, initial flexural cracks, cracking at yielding load, and cracking at flexural capacity are shown in Figure 9. Nonconvergence of the nonlinear algorithm occurred at a load of 1003 kN, indicating flexural failure. The excessive cracking that occurred at this load stage throughout the entire span is evident in Figure 9(c). The FEM was validated by comparing various parameters with hand calculations and strain compatibility approach. Table 3 shows sample calculations for Configuration 1 girder. The corresponding loads and deflections values were close in general, validating the results from the FEM.

**(a) Prestress and self-weight**

**(b) Initial cracking**

**(c) At prestress yielding**

**(d) At flexural capacity**

##### 3.2. Flexural Strengthening (Configuration 3)

In the flexural FEM, the point load was applied at mid-span, as shown in Figure 6(a). The load-deflection curve for this girder is shown in Figure 10, and various cracking stages are shown in Figure 11. The cracking, yielding, and flexural failure loads increased by 17.65%, 5.14%, and 38%, respectively, due to FRP flexural strengthening, as compared to the control unstrengthened girder (Table 4). The maximum deflection at the failure load decreased by 35% due to FRP wrapping, and localized cracking at the girder ends due to the fact that prestress was reduced. The height of the flexural cracks and the area of concrete subjected to flexural cracks also diminished, as compared to the control girder.

**(a) Prestress and self-weight**

**(b) Initial cracking**

**(c) At prestress yielding**

**(d) At capacity**

The analytical procedure from ACI 440-2R [2] was used herein to predict the flexural capacity of the strengthened girder, corresponding well with the FEM prediction. The detailed equations are quite few and are readily available from the reference. The ACI 440 prediction was close to that from the FEM. The effective failure strain prediction from ACI 440 was 0.008786, again close to the FEM output, as was the maximum FRP stress (Table 4). The flexural failure mode of the girder was governed by yielding of the prestressing steel, followed by concrete crushing and FRP rupture, as per the FEM and ACI 440. This may be inferred from the fact that the tensile strength of the FRP (931 MPa) was greater than the actual failure stress in the FRP (531 MPa). FRP strengthening caused 38% and 24% increases in the flexural load capacity of the girder, according to the FEM and the ACI 440 approach, respectively.

##### 3.3. Shear Strengthening (Configurations 2, 4, and 5)

For the shear strengthening configurations, the point load in the FEM was applied at 1.07 m (critical shear location according to AASHTO) from the girder end at one side (Figure 6(b)) to provide critical shear force effect on the girder. The load-displacement curve is presented in Figure 12, and pertinent results are presented in Table 5. Because of the low maximum bending moment near mid-span due to mostly self-weight, flexural failure did not control over shear failure in this FEM. The crack progression in the shear strengthened girders showed reduction in the height of the shear cracks, as compared to the control girder. The ultimate shear capacity of the vertical U-wrapped FRP strengthened girder was about 15% greater from the ACI 440 provisions, as compared to the FEM results. The FEM stress in the prestressing steel at failure was 435 MPa, much below the yielding capacity, as expected in a shear failure situation. The maximum vertical U-wrap FRP strain at failure was 0.0018 from the FEM and 0.002 from ACI 440 calculations, again showing that the shear failure was caused by stirrup yielding and not FRP rupture. The shear load capacity increased by about 13% due to the vertical wrapping. The ACI 440 shear capacity prediction was about 14% greater than the FEM output.

The influence of a 45-degree shear U-wrap FRP on the girder capacity is also presented in Table 5. The shear failure load capacity of the girder with the inclined FRP orientation increased by about 25% over that from the vertical orientation, and about 38% over the control girder shear capacity.

It is noted that 57% of the FRP tensile strength was used in flexural strengthening, while 74% was used in case of shear strengthening. This is because the FRP rupture mode did not control the girder failures herein, and the FRP strength was not fully utilized. FRP application reduced concrete cracking and girder deflection.

FRP strengthening effectively increased the girder load carrying capacity and reduced cracking and deflections. The prestressing steel strain at the nominal strength should be checked to maintain sufficient degree of ductility, according to ACI 440 [2]. Adequate ductility is achieved if the strain in the prestressing steel at the nominal strength is at least 0.013, a bit higher than 0.011 from the FEM for flexural strengthening. To address this, the flexural strength reduction factor has to be decreased to 0.77 in design, according to ACI 440. To increase the girder ductility, the prestressing steel and the FRP areas could be reduced compared to the levels used in this study.

#### 4. Conclusions

The following conclusions may be made based on the results from this study:(1)An AASHTO type IV prestressed concrete bridge girder was successfully modeled using finite element software ANSYS for both flexural and shear FRP strengthening. The approach is simpler and more economic than full scale experimental testing.(2)The FEM results are accurate, based on the validation of the results with hand calculations. Validation parameters included camber due to the initial prestress and after application of the self-weight of the girder, zero deflection, decompression, initial cracking, and failure loads.(3)Excellent contributions from the FRP flexural strengthening to the bridge girder were evident for cracking, yielding, and flexural capacity loads. The flexural failure load increased by almost 38% due to FRP wrapping, and the maximum deflection decreased by 35%. Tension cracks at the girder ends reduced due to FRP, and the height/extent of flexural cracks diminished.(4)The analytical procedure from ACI 440 [2] yielded comparable results to the FEM output for flexural FRP strengthening in terms of the FRP strain, FRP stress, and failure load. The ACI 440 predictions were slightly conservative, as compared to the FEM results.(5)Prestress steel yielding was the dominant failure mode for both the ACI 440 and FEM simulations, followed by concrete crushing. About 57% of the FRP tensile strength was utilized at flexural failure. The amounts of FRP and prestressing steel could be decreased within the design constraints to increase the prestressing steel strain and girder ductility.(6)The FRP shear wrapping also contributed well to the girder shear capacity. The FEM results showed that the shear capacity increased by about 13% due to vertical shear wrapping and 38% due to inclined shear wrapping, respectively. The prestressing steel stress was well below yielding at failure, as expected.(7)The yielding of stirrups was the dominant failure mode in shear, as demonstrated by both the FEM and the ACI 440 outputs. The ACI 440 provisions, however, were found to be a bit on the nonconservative side for shear strengthening, as compared to the FEM outputs. The FRP strains at failure from the two approaches were comparable.(8)Both FEM and ACI 440 analytical techniques may be effectively used to predict the flexural and shear behavior of CFRP strengthened AASHTO type bridge girders.

#### Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.