Assessment of The Strengthening of An RC Railway Bridge With CFRP Utilizing A Full-Scale Failure Test and Finite-Element Analysis
Assessment of The Strengthening of An RC Railway Bridge With CFRP Utilizing A Full-Scale Failure Test and Finite-Element Analysis
Abstract: A finite element (FE) model was calibrated using the data obtained from a full-scale test to failure of a 50 year old reinforced
concrete (RC) railway bridge. The model was then used to assess the effectiveness of various strengthening schemes to increase the load-
carrying capacity of the bridge. The bridge was a two-span continuous single-track trough bridge with a total length of 30 m, situated in
Örnsköldsvik in northern Sweden. It was tested in situ as the bridge had been closed following the construction of a new section of the railway
line. The test was planned to evaluate and calibrate models to predict the load-carrying capacity of the bridge and assess the strengthening
schemes originally developed by the European research project called Sustainable bridges. The objective of the test was to investigate shear
failure, rather than bending failure for which good calibrated models are already available. To that end, the bridge was strengthened in flexure
before the test using near-surface mounted square section carbon fiber reinforced polymer (CFRP) bars. The ultimate failure mechanism
turned into an interesting combination of bending, shear, torsion, and bond failures at an applied load of 11.7 MN (2,630 kips). A computer
model was developed using specialized software to represent the response of the bridge during the test. It was calibrated using data from the
test and was then used to calculate the actual capacity of the bridge in terms of train loading using the current Swedish load model which
specifies a 330 kN (74 kips) axle weight. These calculations show that the unstrengthened bridge could sustain a load 4.7 times greater than
the current load requirements (which is over six times the original design loading), whilst the strengthened bridge could sustain a load
6.5 times greater than currently required. Comparisons are also made with calculations using codes from Canada, Europe, and the United
States. DOI: 10.1061/(ASCE)ST.1943-541X.0001116. This work is made available under the terms of the Creative Commons Attribution
4.0 International license, http://creativecommons.org/licenses/by/4.0/.
Author keywords: Bridge; Train load; Failure analysis; Ultimate load-carrying capacity; Shear; Near-surface mounted reinforcement (NSMR);
Carbon fiber reinforced polymers (CFRP); Strengthening; Nonlinear finite element analysis (NLFEA); Structural safety and reliability.
carried out large-scale tests of the load-carrying capacity of box C55/67. This substantial increase is due the fact that the original
girder beams with thin webs including posttensioning tendons. cement was coarsely grinded and kept on hydrating and growing
The tendons decreased the load-carrying capacity of the compres- in strength after the 28 days when the initial strength was tested.
sion struts necessary to transform the shear forces. A concrete damage plasticity model was used in the finite
element calculations for the bridge slab and mid columns. The
model was chosen mostly because other users had obtained good
Details of the Bridge results with it (FIB 2008). The following properties were assumed,
based on the tested material properties: Young’s modulus of elas-
The bridge studied was a continuous curved reinforced concrete ticity for the concrete Ec ¼ 25.4 GPa; Poisson’s ratio ν ¼ 0.167;
(RC) trough bridge with two spans of 12 m each, designed to carry the dilatation angle β ¼ 35°; the flow potential eccentricity
a single railway line, see Fig. 1. The bridge was designed and ε ¼ 0.1; and the biaxial/uniaxial compression plastic strain ratio
built in 1955 and was taken out of service in 2005 because of the fbo =f c ¼ 1.16 and the invariant stress ratio κ ¼ 0.6666 (Puurula
building of a new high-speed railway, the Bothnia line. Before 2012).
demolition, the bridge was loaded to failure to test its ultimate The steel reinforcement bars and the surrounding concrete were
load-carrying capacity as part of the European research project sus- modeled together, increasing the nominal virtual stiffness of the
tainable bridges [Sustainable Bridges (SB) 2008]. The testing was steel up to the stress when the concrete cracks, see Fig. 6. This
carried out by applying loads on a beam perpendicular to the bridge was done in accordance with the results of RILEM Committee
(Fig. 2). The embankment south of the bridge was removed before 147-FMB “Fracture mechanics to Anchorage and Bond” (Elfgren
the test. The loads were applied with jacks anchored in the ground and Noghabai 2001, 2002). This procedure increased the stiffness
beneath the bridge (Sustainable Bridges (SB)-7.3 2008). The bridge of the calculated load-deflection diagram so that it better followed
is skewed at an angle of 75° between its longitudinal and transverse the curve from the test.
directions. It is also curved in the horizontal plane with a radius A summary of the material properties is presented in Table 1.
of 300 m. Fig. 3 is of a longitudinal section showing the steel Initial characteristic properties are given first based on the original
reinforcement. drawings, followed by updated properties based on mean values of
To avoid a pure bending failure, for which good, calibrated the tested samples taken from the bridge following the load test to
models already exist, the bridge was strengthened before testing failure.
with bars of carbon fiber reinforced polymers (CFRP) (Täljsten
et al. 2011). The two edge beams of the bridge were each fitted
with nine Sto FRP Bar M10C with a length of 10 m and a Finite Element Model and Calibration with a
rectangular cross-section of 10 × 10 mm. They were installed on Full-Scale Field Test
100-mm centers using the near-surface mounted reinforcement
technique (NSMR) in presawn grooves, 15 × 15 mm, in the soffit The Örnsköldsvik Bridge was modeled with successively improved
of the bridge (Figs. 4 and 5) (Sustainable Bridges (SB)-6.3 2007). models, starting with linear two-dimensional frame models and
Fig. 1. View of bridge in Örnsköldsvik in northern Sweden prior to testing (adapted from Sustainable Bridges (SB) 7.3 2008; image by
Lennart Elfgren)
Fig. 2. Elevation, plan, and section of the bridge showing the test loading arrangement using a steel beam placed in the middle of one of the two spans
and pulled downwards (adapted from Sustainable Bridges (SB) 7.3 2008)
ending with a nonlinear three-dimensional finite element model bonded to the surrounding concrete. The steel beam used to intro-
(Fig. 7) using Brigade (2011), which is based on Abaqus software. duce the load on the bridge was modeled using shell elements of
The calculation models were calibrated with results from the full- type S4R: linear quadrilateral, four-node doubly curved shell, re-
scale field test of the bridge (Puurula 2012). The boundary condi- duced integration, and hourglass control. The piles were modeled
tions and the nonlinear material properties during yielding of the as springs, each inclined pile as a separate spring with a stiffness in
steel reinforcement close to failure were deemed important param- the vertical direction of kv ¼ 285 MN=m and in the horizontal
eters to calibrate. direction of kh ¼ 71 MN=m. The earth pressure σh on the East
The model had 1,650 separate structural parts (most of them abutment was modeled as σh ¼ k0 · γ · h, where k0 ¼ 0.34 is a co-
discrete reinforcement bars), 152,460 elements, 164,003 nodes, and efficient for the soil pressure, γ ¼ 20 kN=m3 is the weight of the
511,317 variables. Solid elements in the concrete bridge were of earth, and h is the height from the Earth’s surface (m) (Puurula
type continuum, 3-dimensional, 8-node, reduced integration 2012). As discussed, the material modeling and the boundary con-
(C3D8R): 8-node linear brick, reduced integration, and hourglass ditions were important parameters in the calibration of the model.
control. Parameter studies were carried out with different element The contact between the steel beam and the concrete was modeled
sizes; elements smaller than 150 mm did not improve the results. with a tie constraint, which does not allow the contact surfaces to
Discrete reinforcement bars were modeled as wires, type two-node move in relation to one another. Another issue was to calculate
linear 3D truss elements embedded in the concrete. The CFRP the overall load-deflection curve of the bridge. Here, it was essen-
reinforcement bars in their grooves were modeled as perfectly tial to consider the effect of the surrounding concrete in Fig. 6
Fig. 3. Longitudinal section of the bridge showing the layout of the main steel reinforcement and the outline strengthening scheme
Fig. 4. Cross section showing the principal dimensions of the bridge and reinforcement details together with the location of the near-surface mounted
reinforcement (NSMR) with FRP Bar M10C on 100 mm centers
Fig. 5. Installation of the NSMR showing (a) sawing of grooves; (b) filling grooves with epoxy adhesive; (c) grooves following insertion of the CFRP
reinforcement
and the choice of damage parameters. Dilatation angles between (Figs. 1 and 2). The monitoring system consisted primarily of strain
β ¼ 10° and 50° gave similar results. Convergence problems were gauges that were spot-welded to the reinforcement and glued to the
addressed using Riks method (FIB 2008). CFRP bars and the concrete, an optical laser displacement sensor
The bridge was tested in July 2006. The load was applied using and linear varying differential transducers (LVDTs) Sustainable
two jacks on top of a steel beam, which was pulled downwards Bridges (SB)-7.3 2008. The load-deflection curve from the final
out considering the effect of surrounding concrete (data from Elfgren crack had an inclination of approximately 35° with respect to the
and Noghabai 2001, 2002) horizontal axis.
A preliminary description of the load test and the finite element
calculations is given in Puurula et al. (2008) and in more detail in
test is given in Fig. 8, which also shows the calculated load- the Ph.D. theses of Sas (2011) and Puurula (2012) and in Sas
deflection curve and the effect of strengthening. It is shown that (2012) and Puurula et al. (2013).
the two curves follow one another closely and that the bridge ex-
hibits ductile behavior with a large deflection of the order of 0.1 m
before failure. The calculated strains in the steel and CFRP rein- Comparison with Codes
forcement and in the concrete also correspond well to the measured
values. The load capacity of the bridge was calculated using three
At the time of failure, high bond stresses between the concrete major codes from the United States, Canada, and Europe, namely,
and the resin in the outermost groove initiated a bond failure after ACI-318 [American Concrete Institute (ACI) 2011], CSA-A23.3
yielding of the bottom longitudinal steel reinforcement. The bond [Canadian Standards Association (CSA) 2004], and EC2 [European
stresses were calculated to be 11.3 MPa with an alternative, refined Committee for Standardization (CEN) 2004], respectively. The
model in which the CFRP reinforcement was embedded in epoxy background to them is given in ASCE-ACI 445 (1998).
Fig. 7. (a) Model of bridge with nonlinear concrete in the bridge slab and midcolumns and linear concrete in the other parts; failure started in the east
edge beam; (b) nonlinear discrete reinforcement is embedded in the concrete model of the bridge with a perfect bond to the concrete
Fig. 8. Comparison of load-deflection curves for the cases with and data obtained from beams with only steel reinforcement, whereas
without FRP strengthening; 1 = effect on stiffness, 2 = strengthening the bridge was also strengthened with CFRP reinforcement, a
effect material that displays linear elasticity until failure. This might be
one reason why the crack angle estimated by the CSA code does
not correspond to the angle observed in the test.
The interaction between the shear force and the bending The European code compensates for the omission of the concrete
moment, as a function of a unit load P, is presented in Fig. 10. contribution to the shear capacity by adopting a crack angle that
Because the dead load was already acting on the bridge, its effect most closely matches the angle observed in the test. In addition,
[the hatched parts of the diagrams in Figs. 10(b and c)] has not the truss model used in EC2 is a transparent geometrical method
been considered in the analysis. Calculations based on the initial and the change in the tensile longitudinal strains because of the ad-
characteristic concrete strength, fc ¼ 31 MPa, and the tested mean dition of the strengthening can be easily incorporated into the analy-
value, fc ¼ 68 MPa, are given in Sas (2011), Sas et al. (2011) and sis to obtain cross-sectional equilibrium. As the longitudinal force
some of the results are summarized in Table 2. in the tensile chord of the truss increases, the crack angle is reduced;
The three codes predict the shear force capacity and the ultimate therefore, the assumed crack has to bridge more stirrups to obtain
load capacity of the bridge in a conservative manner. The ratio be- equilibrium. In this way, EC2 predicts an increase of the shear force
tween the predicted value PV and the test result PTest varies between capacity after strengthening. However, disregarding the concrete
0.31 (EC2, with a concrete strut inclination of θ ¼ 45°), 0.65 (CSA, contribution leads to conservative estimates of the shear capacity.
θ ¼ 38°), 0.66 (ACI, θ ¼ 45°), and 0.78 (EC2 with minimum When the initial characteristic compressive strength of the con-
value θ ¼ 22°). Because EC2 makes use of the variable angle truss crete is used, the discrepancies between code predictions and test
model, both minimum and maximum capacities were estimated. results are even larger (Table 2). This underlines the importance of
They are shown in Table 2. The reason for the differences is the using actual tested values in assessment of structures.
Fig. 9. Failure cracks with ruptured stirrups in the beams after the maximum load of 11.7 MN in the span close to the South abutment: (a) west beams;
(b) east beams (images by Lennart Elfgren)
Table 2. Capacities Predicted by Codes Using Tested and Nominal Characteristic Material Properties
ACI CSA EC2
Parameter studied fc ¼ 68 MPa f c ¼ 31 MPa fc ¼ 68 MPa f c ¼ 31 MPa f c ¼ 68 MPa fc ¼ 31 MPa
θ 45° 45° 38°a 37°4′a 45° 22° 45° 22°
V s (MN) 1.06 1.06 1.31 1.22 1.02 2.54 0.93 2.32
V c (MN) 1.24 0.83 0.97 0.66 — — — —
V Rd (MN) 2.29 1.89 2.28 1.87 1.02 2.54 0.93 2.32
M Rd (MNm) 6.12 5.62 6.12 5.62 6.12 6.12 5.62 5.62
PV (MN) 7.74 6.25 7.67 6.3 3.64 9.12 3.33 8.33
PV =PTest 0.66 0.53 0.65 0.54 0.31 0.78 0.28 0.71
Note: M Rd = bending moment capacity; PV = vertical load corresponding to the shear capacity; and PV =PTest = ratio of the predicted capacity to the actual load
obtained at the PTest ¼ 11.7 MN; V s and V c = vertical shear forces taken by steel and concrete, respectively; V Rd ¼ V s þ V c is the total shear capacity; and
θ = inclination of the compression struts.
a
This value iteratively calculated.
Fig. 11. Proposed strengthening of the slab with Sto FRP Bar M10C at 150 mm centers, parallel to the supports but not parallel to the steel
reinforcement
Fig. 12. Modeled vertical displacements of the bridge slab for an unstrengthened bridge at the load level q ¼ 304 kN=m2 ¼ 4.71q0 , where
q0 ¼ 64.5 kN=m2 is the mean live load for a train with an axle load of 330 kN. The maximum displacement in the midpoint of the slab is 0,1028 m
Fig. 13. Maximum tensile stress of 629 MPa, well above the yield stress of 440 MPa, is obtained for a load level of q ¼ 304 kN=m2 ¼ 4.71q0 for an
unstrengthened bridge
The effect of the strengthening of the slab is further demon- This paper has described how the load-carrying capacity of a re-
strated in Fig. 15, in which load-stress curves are given for the inforced concrete railway bridge can be calculated for different
transverse bottom slab steel reinforcement, showing that the CFRP strengthening alternatives using a nonlinear three-dimensional fi-
strengthening of the slab decreases the steel tensile stress. In the nite element model with discrete reinforcement bars which was
case in which the slab is strengthened with transverse CFRP bars, calibrated using a full-scale test to failure of a redundant 50 year
much of the tensile force is taken by the carbon reinforcement es- old bridge. The strengthening of the bridge was successful as the
pecially after the steel reinforcement started to yield at a load level carbon fiber bars increased both the stiffness of the bridge and its
of 230 kN=m2 . The CFRP reinforcement reaches its ultimate fail- bending moment capacity.
ure stress at a load level of q ¼ 487 kN=m2 (or 7.54q0 ). After ini- The three-dimensional nonlinear calculation method with dis-
tial failure, the stress in the carbon bar decreases and the stresses in crete reinforcement was used to model the bridge in an integrated
the steel reinforcement start to increase, now in the hardening part way, including twist and deflections in all directions. The calcula-
of the stress-strain curve (Fig. 16). Just strengthening the edge tion model was first calibrated with test results from a full-scale
beams does not have any effect on the stresses in the tensile steel loading to failure of a railway bridge in Örnsköldsvik in northern
reinforcement in the slab. Sweden. The bridge behavior with increasing load and its ultimate
Load-strain curves for the transverse bottom reinforcement load-carrying capacity were closely predicted using the calculation
in the middle of the slab for increasing train load are shown in method described. The focus was on the assessment of the load-
Fig. 16. The steel reinforcement was located squarely in the slab carrying capacity during nonlinear behavior of the materials rather
Fig. 17. Plastic strains in concrete, i.e., the crack pattern of a strengthened slab at the failure load as seen from underneath the slab