Introduction
To meet future traffic demands, there is a constant need to make transport infrastructure more effective. For railway bridges, this can be achieved by increasing their load-carrying capacity to allow heavier and/or faster trains (where higher dynamic factors increase the static load effect on a bridge) and/or to increase their service life.
This paper describes the assessment of the load capacity (in terms of trains) using a calculation method that has been calibrated using a full-scale test to failure of a redundant concrete railway bridge in Örnsköldsvik (Övik) in northern Sweden. The method uses a nonlinear finite element analysis (NLFEA) with a detailed three-dimensional (3D) solid model using discrete reinforcement elements with both geometrical and material nonlinearities included. The model developed gave a detailed picture of deflections, stresses, and strains in the bridge with increasing train loads, which gave a better correlation between calculated results and the actual behavior of the bridge than the code methods currently in use. Whilst this is thought to be the first use of this method of assessment of the train load capacity of a bridge, the basics of the method are given in FIB (
2008), and examples of advanced finite element modeling of large concrete structures are given in, e.g., Malm (
2009), Schlune et al. (
2009), Richard et al. (
2010), and Bicanic et al. (
2014).
Reports on full-scale tests to failure studying the behavior of bridges are rare. In this paper, some of them are discussed with failures of relevance to this study. A large-scale test (
) investigating the combined shear—torsion—bending of a curved concrete box-girder bridge in California was described by Scordelis et al. (
1977,
1979). Failure was caused by concrete spall from a corner because of high shear and torsion stresses. Two full-scale tests of concrete bridges in Sweden were reported by Plos (
1990,
1995) and Täljsten (
1994). A slab frame bridge with a span of 21 m and a prestressed frame beam bridge with a span of 31 m were both tested to failure, which was because of shear and bending. For the slab beam, a brittle failure occurred when a new shear crack with a low slope emerged at a point load of 4.5 MN. For the prestressed bridge, the failure occurred when one of the beams punched through the end support wall at a point load of 8.45 MN. Neither of the failures were predicted by the available codes. In Switzerland, 89 bridges were studied during their demolition by Zwicky and Vogel (
2000), and by Vogel and Bargähr (
2006). Both sets of authors recommended improved methods for the assessment of real loads and for corrosion damage in prestressed concrete bridges. The tests were also simulated using solid elements by Pimentel et al. (
2007) and using fiber elements capable of taking into account shear effects by Ferreira et al. (
2012). In Switzerland, Fernandez Ruiz et al. (
2007) carried out large-scale tests of the load-carrying capacity of box girder beams with thin webs including posttensioning tendons. The tendons decreased the load-carrying capacity of the compression struts necessary to transform the shear forces.
Details of the Bridge
The bridge studied was a continuous curved reinforced concrete (RC) trough bridge with two spans of 12 m each, designed to carry a single railway line, see Fig.
1. The bridge was designed and built in 1955 and was taken out of service in 2005 because of the building of a new high-speed railway, the Bothnia line. Before demolition, the bridge was loaded to failure to test its ultimate load-carrying capacity as part of the European research project sustainable bridges [
Sustainable Bridges (SB) 2008]. The testing was carried out by applying loads on a beam perpendicular to the bridge (Fig.
2). The embankment south of the bridge was removed before the test. The loads were applied with jacks anchored in the ground beneath the bridge (
Sustainable Bridges (SB)-7.3 2008). The bridge is skewed at an angle of 75° between its longitudinal and transverse directions. It is also curved in the horizontal plane with a radius of 300 m. Fig.
3 is of a longitudinal section showing the steel reinforcement.
To avoid a pure bending failure, for which good, calibrated models already exist, the bridge was strengthened before testing 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 rectangular cross-section of
. They were installed on 100-mm centers using the near-surface mounted reinforcement technique (NSMR) in presawn grooves,
, in the soffit of the bridge (Figs.
4 and
5) (
Sustainable Bridges (SB)-6.3 2007).
Material Properties
The original concrete quality used in 1955 was Swedish K400, with a nominal compressive strength of 40 MPa () measured on 200-mm cubes. This corresponds to a characteristic strength of 31 MPa (lower 5% percentile) and approximately to EC class C28/35. The steel reinforcement is mostly made from 16 and 25-mm diameter bars of quality Ks40 with a nominal yield strength of 400 MPa (58 ksi). Carbon fiber reinforcement Sto FRP Bar M10C has an and a mean tensile strength (363 ksi).
Separate mean concrete properties were determined by testing drilled core samples from each edge beam and the slab which showed that, by the date of the test in 2006, the concrete strength had increased to 68.5 MPa (9,935 psi) corresponding to EC class C55/67. This substantial increase is due the fact that the original cement was coarsely grinded and kept on hydrating and growing in strength after the 28 days when the initial strength was tested.
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 results with it (
FIB 2008). The following properties were assumed, based on the tested material properties: Young’s modulus of elasticity for the concrete
; Poisson’s ratio
; the dilatation angle
; the flow potential eccentricity
; and the biaxial/uniaxial compression plastic strain ratio
and the invariant stress ratio
(
Puurula 2012).
The steel reinforcement bars and the surrounding concrete were modeled together, increasing the nominal virtual stiffness of the steel up to the stress when the concrete cracks, see Fig.
6. This was done in accordance with the results of RILEM Committee 147-FMB “Fracture mechanics to Anchorage and Bond” (
Elfgren and Noghabai 2001,
2002). This procedure increased the stiffness of the calculated load-deflection diagram so that it better followed the curve from the test.
A summary of the material properties is presented in Table
1. Initial characteristic properties are given first based on the original drawings, followed by updated properties based on mean values of the tested samples taken from the bridge following the load test to failure.
Finite Element Model and Calibration with a Full-Scale Field Test
The Örnsköldsvik Bridge was modeled with successively improved models, starting with linear two-dimensional frame models and ending with a nonlinear three-dimensional finite element model (Fig.
7) using Brigade (
2011), which is based on Abaqus software. The calculation models were calibrated with results from the full-scale field test of the bridge (
Puurula 2012). The boundary conditions and the nonlinear material properties during yielding of the steel reinforcement close to failure were deemed important parameters to calibrate.
The model had 1,650 separate structural parts (most of them discrete reinforcement bars), 152,460 elements, 164,003 nodes, and 511,317 variables. Solid elements in the concrete bridge were of type continuum, 3-dimensional, 8-node, reduced integration (C3D8R): 8-node linear brick, reduced integration, and hourglass control. Parameter studies were carried out with different element sizes; elements smaller than 150 mm did not improve the results. Discrete reinforcement bars were modeled as wires, type two-node linear 3D truss elements embedded in the concrete. The CFRP reinforcement bars in their grooves were modeled as perfectly bonded to the surrounding concrete. The steel beam used to introduce the load on the bridge was modeled using shell elements of type S4R: linear quadrilateral, four-node doubly curved shell, reduced integration, and hourglass control. The piles were modeled as springs, each inclined pile as a separate spring with a stiffness in the vertical direction of
and in the horizontal direction of
. The earth pressure
on the East abutment was modeled as
, where
is a coefficient for the soil pressure,
is the weight of the earth, and
is the height from the Earth’s surface (
) (
Puurula 2012). As discussed, the material modeling and the boundary conditions were important parameters in the calibration of the model. The contact between the steel beam and the concrete was modeled with a tie constraint, which does not allow the contact surfaces to move in relation to one another. Another issue was to calculate the overall load-deflection curve of the bridge. Here, it was essential to consider the effect of the surrounding concrete in Fig.
6 and the choice of damage parameters. Dilatation angles between
and 50° gave similar results. Convergence problems were addressed using Riks method (
FIB 2008).
The bridge was tested in July 2006. The load was applied using two jacks on top of a steel beam, which was pulled downwards (Figs.
1 and
2). The monitoring system consisted primarily of strain gauges that were spot-welded to the reinforcement and glued to the CFRP bars and the concrete, an optical laser displacement sensor and linear varying differential transducers (LVDTs)
Sustainable Bridges (SB)-7.3 2008. The load-deflection curve from the final test is given in Fig.
8, which also shows the calculated load-deflection curve and the effect of strengthening. It is shown that the two curves follow one another closely and that the bridge exhibits ductile behavior with a large deflection of the order of 0.1 m before failure. The calculated strains in the steel and CFRP reinforcement and in the concrete also correspond well to the measured values.
At the time of failure, high bond stresses between the concrete and the resin in the outermost groove initiated a bond failure after yielding of the bottom longitudinal steel reinforcement. The bond stresses were calculated to be 11.3 MPa with an alternative, refined model in which the CFRP reinforcement was embedded in epoxy (
Puurula 2012). This is higher than the bond strength of 9.0 MPa for this type of bar (
Nordin and Täljsten 2003). The bond failure lowered the available tensile force at the bottom and increased the inclination of the concrete compression struts which produced higher stresses in the stirrups, as fewer stirrups had to carry the load. These stresses, mostly caused by the vertical shear forces but also, to some extent, by the torsion moment from the loads transferred from the steel beam to the slab, ruptured the stirrups. The torsion moment originates primarily from the load of the steel beam. The edge beams twist outwards because they are nonsymmetrically supported on the bridge slab. The outsides of the edge beams deflect the most because of torsion, which explains why the final failure started in the outermost groove of the CFRP bars (
Puurula et al. 2013). The ultimate load capacity was reached at an applied midspan load of 11.7 MN, see Fig.
9. The inclined failure crack had an inclination of approximately 35° with respect to the 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 the Ph.D. theses of Sas (
2011) and Puurula (
2012) and in Sas (
2012) and Puurula et al. (
2013).
Comparison with Codes
The interaction between the shear force and the bending moment, as a function of a unit load
, is presented in Fig.
10. Because the dead load was already acting on the bridge, its effect [the hatched parts of the diagrams in Figs.
10(b and c)] has not been considered in the analysis. Calculations based on the initial characteristic concrete strength,
, and the tested mean value,
, are given in Sas (
2011), Sas et al. (
2011) and some of the results are summarized in Table
2.
The three codes predict the shear force capacity and the ultimate load capacity of the bridge in a conservative manner. The ratio between the predicted value
and the test result
varies between 0.31 (EC2, with a concrete strut inclination of
), 0.65 (CSA,
), 0.66 (ACI,
), and 0.78 (EC2 with minimum value
). Because EC2 makes use of the variable angle truss model, both minimum and maximum capacities were estimated. They are shown in Table
2. The reason for the differences is the way the shear truss mechanism is applied. Codes ACI and CSA permit the use of the concrete’s contribution to the shear capacity, whereas in EC2, this is not permitted. Another factor that is responsible for some of the differences is that the codes differ in their treatment of the concrete compression strut inclination. In the ACI code, it is assumed to be fixed at 45°, which is quite a conservative assumption. The Canadian code is a simplified version of the modified compression field theory (MFCT), e.g., Collins and Mitchell (
1991) and Bentz et al. (
2006). Here, the inclination is determined iteratively from the cross-sectional equilibrium and is dependent on factors such as the crack spacing, concrete material properties and the average tensile and compressive strains over the cracked sections. In the simplified code definitions, the inclination depends primarily on the longitudinal strain
calculated at the middle of the cross-section. However, this definition has been calibrated using data obtained from beams with only steel reinforcement, whereas the bridge was also strengthened with CFRP reinforcement, a 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 European code compensates for the omission of the concrete contribution to the shear capacity by adopting a crack angle that most closely matches the angle observed in the test. In addition, the truss model used in EC2 is a transparent geometrical method and the change in the tensile longitudinal strains because of the addition of the strengthening can be easily incorporated into the analysis to obtain cross-sectional equilibrium. As the longitudinal force in the tensile chord of the truss increases, the crack angle is reduced; therefore, the assumed crack has to bridge more stirrups to obtain equilibrium. In this way, EC2 predicts an increase of the shear force capacity after strengthening. However, disregarding the concrete contribution leads to conservative estimates of the shear capacity.
When the initial characteristic compressive strength of the concrete is used, the discrepancies between code predictions and test results are even larger (Table
2). This underlines the importance of using actual tested values in assessment of structures.
Assessment of Load-Carrying Capacity for a Train Load
A very simple first estimate of the load capacity of the bridge can be obtained by dividing the actual failure load, , by the originally designed-for axle load of 250 kN. This gives the load capacity of the bridge as . On a 12-m span, there is room for one locomotive with six axles, spaced at 1.6 m centers. This results in the approximation that the strengthened bridge may carry the equivalent of locomotives.
However, by using the finite element model presented previously, now including the soil pressure at both ends of the bridge, it is possible to assess the capacity of the bridge more accurately. In the calculations using Brigade, the train load is increased successively and, at the same time, the uniformly distributed soil pressure caused by the increasing train load is also increased.
The Swedish railways bridge design code, BV Bridge (
2003), contains the BV 2000 loading model, which represents trains with axle loads of 330 kN and is more representative of the type of traffic actually using the bridge. For design and assessment purposes, the actual train loading can be represented by a uniformly distributed characteristic line load of
where 6.4 m is the distance over which the four-axle loads act in two bogies connected to one another. This gives a mean value of the train load of
with a standard deviation of
if the characteristic value has a probability of being exceeded by 3.0% of the trains. As the width of the bridge is 2.9 m, the design loading can be expressed as
. The uniformly distributed soil pressure is calculated as
, where
= train load applied on the bridge divided by the deck area and 0.34 is the coefficient for soil pressure. Additionally, permanent soil pressure, self-weight of the bridge ballast, sleepers, and rail are included in the calculation (Puurula
2012).
Based on results from the model, the following comparisons can be made:
•
bridge without strengthening,
•
two edge beams strengthened with 9 FRP bars as tested during the load test to failure, and
•
bridge slab strengthened with transverse NSM Sto FRP Bar M10C on 150 mm centers. The bars are parallel to the supports, contrary to the original steel reinforcement bars which are at an angle to the supports and perpendicular to the bridge center line (Fig.
11).
In Fig.
12, vertical displacements are presented for the unstrengthened bridge at the train load level
. This load level corresponds to
the modern design loading, representing trains with an axle load of 330 kN or approximately
the original design requirement to cater for an axle load of 250 kN. The deflection in the middle of the bridge slab reached a value of 0.1028 m, when the maximum stresses in the reinforcement in the slab were far beyond their yield limit (Fig.
13). Hence, according to the model, the failure of the original, unstrengthened bridge would start in the deck slab.
The influence of the two potential strengthening alternatives is illustrated in Figs.
14–18. Comparable load-deflection curves for the midpoint of the slab are given in Fig.
14. It is shown that strengthening the edge beams only contributes to a small stiffening of the slab. Strengthening of the slab, however, has a substantial effect on the slab stiffness. The strengthening of the slab also considerably increases the overall load-carrying capacity of the bridge. Failure deflections in the order of 0.1 m are large but are indicative of ductile behavior despite the introduction of nonductile CFRP reinforcement.
The effect of the strengthening of the slab is further demonstrated in Fig.
15, in which load-stress curves are given for the transverse bottom slab steel reinforcement, showing that the CFRP strengthening of the slab decreases the steel tensile stress. In the case in which the slab is strengthened with transverse CFRP bars, much of the tensile force is taken by the carbon reinforcement especially after the steel reinforcement started to yield at a load level of
. The CFRP reinforcement reaches its ultimate failure stress at a load level of
(or
). After initial failure, the stress in the carbon bar decreases and the stresses in the steel reinforcement start to increase, now in the hardening part of the stress-strain curve (Fig.
16). Just strengthening the edge beams does not have any effect on the stresses in the tensile steel reinforcement in the slab.
Load-strain curves for the transverse bottom reinforcement in the middle of the slab for increasing train load are shown in Fig.
16. The steel reinforcement was located squarely in the slab to minimize the distance to the edge beams, whilst the CFRP reinforcement was placed slightly skewed, to be parallel with the supports, see Fig.
11. This is the reason why the steel reinforcement carries larger strains than the carbon reinforcement. After the failure of the carbon bars, the steel reinforcement takes over the load from the carbon bars and the strains in steel reinforcement increase rapidly for the same load level.
Plastic strains in the concrete at the failure load as seen from underneath are shown in Fig.
17 for the strengthened slab. It gives an indication of the crack pattern. In Fig.
18, load-strain curves are given for the longitudinal reinforcement in the centre of the East edge beam when the slab is strengthened transversely with CFRP bars. The steel tensile reinforcement in the edge beams begins to yield at the load level
.
In summary, the transverse carbon strengthening chosen was quite appropriate because the strengthened slab can sustain , and the steel reinforcement in the edge beams can sustain before yielding starts. The rupture of the carbon at is very brittle and is to be avoided. In this case, the yielding of the steel reinforcement in the edge beams leads to a ductile failure. This can, therefore, be predicted to be the ultimate limit of the bridge’s load-carrying capacity.
Summary and Conclusions
This paper has described how the load-carrying capacity of a reinforced concrete railway bridge can be calculated for different strengthening alternatives using a nonlinear three-dimensional finite element model with discrete reinforcement bars which was calibrated using a full-scale test to failure of a redundant 50 year old bridge. The strengthening of the bridge was successful as the carbon fiber bars increased both the stiffness of the bridge and its bending moment capacity.
The three-dimensional nonlinear calculation method with discrete reinforcement was used to model the bridge in an integrated way, including twist and deflections in all directions. The calculation model was first calibrated with test results from a full-scale loading to failure of a railway bridge in Örnsköldsvik in northern Sweden. The bridge behavior with increasing load and its ultimate load-carrying capacity were closely predicted using the calculation method described. The focus was on the assessment of the load-carrying capacity during nonlinear behavior of the materials rather than investigating other failure modes such as intermediate crack (IC) induced debonding, anchorage failure, peeling, or delamination of the carbon fiber reinforcement.
After adjustment and verification of the model, train loading was applied to the bridge using statistical mean values of actual material properties obtained from testing samples taken from the bridge after the conclusion of the load test. The load-carrying capacity of the bridge was then calculated for the bridge before and after strengthening with near-surface mounted reinforcement (NSMR) consisting of carbon fiber reinforced polymers (CFRP). The method can be used to ensure a ductile behavior after strengthening and to avoid brittle failure modes which may be caused by rupturing carbon fibers. Bridge sections for strengthening should be designed so that ductility is maintained by ensuring that the FRP bars do not rupture before extended yielding of the steel reinforcement.
The assessment of the different strengthening alternatives shows that the bridge, after strengthening the deck slab with CFRP bars, would have failed under a load 6.5 times the current maximum axle load of 330 kN, whilst the unstrengthened bridge would have failed under a load 4.7 times the axle load of 330 kN and approximately 6.2 times the originally designed-for axle load of 250 kN. The calculations show that ordinary reinforced concrete bridges of the type studied are ductile and often have a reserve capacity, which may be utilized after a proper assessment procedure. This extra capacity is seldom obvious when a standard evaluation is run using nominal material properties and code models.
Therefore, to increase our knowledge and to use our existing bridges to their best potential, more bridges of different types should be tested to failure in a planned way, when they are going to be demolished. In this way we can (1) analyze the real behavior of typical bridges, (2) develop appropriate strengthening strategies, and (3) determine the real safety of our bridges—with and without strengthening.