Fatigue safety level provided by Brazilian design standards for a prestressed girder highway bridge

São Carlos, SP, Brasil Abstract: It is well known that highway bridges are subjected to fatigue as they work under live loads with different frequencies and amplitude. The safety level for fatigue required by Brazilian codes is still unknown, especially for prestressed concrete girders. Also, current studies on fatigue reliability of bridges only evaluate bending. This work assesses the fatigue safety level provided by Brazilian design standards for a concrete highway bridge, using weigh-in-motion (WIM) data of an important federal Brazilian highway, BR-381 (Fernão Dias Highway). The Palmgren-Miner rule is considered to evaluate the service life and reliability indexes, from the fatigue point of view, of prestressed girders designed according to Brazilian codes. Using limited and complete prestressing levels, different traffic volumes are considered. It is found that the fatigue safety levels of longitudinal and transverse reinforcements are larger than the ones recommended


INTRODUCTION
Fatigue occurs due to stress fluctuations from in-service variable loads that, after a certain number of cycles, may lead to fracture of the material. From the point of view of the loads, fatigue is a Service Limit State -SLS and, from the resistance point of view, it has a character of Ultimate Limit State -ULS. In highway bridges, the number of stress cycles from the traffic determines the fatigue life of the structure. From the design point of view, fatigue is verified considering live load models.
The current Brazilian live load model as in NBR 7188 [1] is shown in Figure 1. This model was originally adopted from the Germany code, as stated by Pfeil [2]. By that time, the Germany live load model was an analogy of a war tank and a platoon crossing the bridge. As the current Brazilian code still uses this model, the procedure developed by Rüsch to estimate the load effect on bridge slabs are still valid and can be applied for the Brazilian live load model. It can be seen from Figure 1 that the contact area of the wheel with the pavement is 0.1 m 2 . According to Brazilian code, the live load model must be weighted by a coefficient related to impact, number of lanes and additional safety coefficients as in NBR 8681 [3]. For fatigue verification in girders, it requires use of a coefficient ψ 1,fad = 0.5 for spans lengths up to 100 meters. Carneiro et al. [4] evaluated the live load model as in NBR 7188 [1], from the unlimited fatigue life point of view, using traffic data from two HS-WIM stations in Brazil (BR-381 and BR-290) and proposed a new fatigue live load model for unlimited fatigue life. The authors indicated that the bias factor of the Brazilian model can vary a lot and may not correspond to the infinite fatigue life approach. Other studies on the assessment of current Brazilian live load model for bridges using real traffic data, are related only to ULS, as seen in Portela [5]. Thus, the reliability indexes of concrete bridges designed with the current Brazilian live load model, especially for prestressed girders, are still unknown in terms of fatigue.
With that in mind, this work assesses the fatigue safety level provided by Brazilian design standards in relation to a concrete highway bridge, using weigh-in-motion (WIM) data of an important federal Brazilian highway, BR-381 (the same data used in Carneiro et al. [4]). The fatigue service life and fatigue reliability indexes for a design life of 50 years are evaluated for prestressed concrete girders designed according to Brazilian codes, in terms of bending moment, shear force and torsion. The cumulative linear damage method, also known as the Palmgren-Miner rule, is considered. For the probabilistic analysis, the Latin Hypercube Sampling -LHS simulation method and the First-Order Reliability Method -FORM, are used. It is important to mention that this paper does not assess the safety of an existing bridge, which real monitoring is more appropriate, but it presents a methodology to evaluate the fatigue safety level provided by design standards.

Literature review of fatigue reliability in concrete bridges
Crespo-Minguillón and Casas [6] present a probabilistic model for fatigue analysis in prestressed bridges. The model uses the S-N curves and the Palmgren-Miner rule to define fatigue strength of reinforcing and prestressing steel. Using traffic data, they assessed the reliability of a slab bridge in terms of bending. They concluded that from the fatigue life point of view the bridge failure was unlikely to occur.
Rodrigues et al. [7] evaluated the fatigue reliability indexes of two-girder concrete bridges in Brazil. Only shortspan (7 m, 10 m and 13 m) bridges were considered in terms of ULS and fatigue. It is important to note, however, that this structural solution, widely used in Brazil over the past years, is no longer the common application in current bridge design. They used traffic data from a static weighing station in São Paulo, from 2005, measured for 204 days. The work concluded that the safety levels, especially in relation to fatigue, are lower than desired. The study only performed bending analysis and the methodology was based on the work of Crespo-Minguillón and Casas [6].
Wassef et al. [8] calibrated the partial safety factors for SLS and fatigue in relation to concrete bridges in the United States. They used data from fifteen different WIM stations around the country, measured for one year. The authors used the infinite fatigue life approach and considered only girder bending.
Yan et al. [9] presented a methodology for assessing the fatigue reliability of short span bridges subjected to real traffic. The authors applied the methodology on a girder and determined the fatigue reliability indexes of longitudinal reinforcements.
Other works like Junges [10], Wang et al. [11], Mankar et al. [12] and [13] presented methodologies and casestudies to assess the fatigue reliability of monitored bridges. Junges [10] measured stress in two Brazilian reinforced concrete bridges and concluded that flexural reinforcement failure is unlikely to occur. Wang et al. [11] also evaluated flexural reinforcement of concrete girders. Mankar et al. [12] and [13] evaluated concrete bridge slabs in terms of fatigue in concrete and reinforcement, respectively. In Mankar et al. [12] the authors used test data to present stochastic models of S-N curves for concrete in compression. In an extensive literature review, no work was found on fatigue reliability assessment of stirrups.

Fatigue resistance
In order to determine the fatigue strength, the S-N curve is used. This curve is a plot of the magnitude of an alternating stress (S) and the number of cycles to failure (N) for a given material, as shown in Equation 1. Figure 2 depicts the curve of Equation 1 in logarithmic scale. The parameters of reinforcement S-N curves from fib [14] model code are shown in Table 1.
where m and K : material-related constants; N : number of cycles to failure; σ ∆ : stress range.  Note: m 1 , m 2 and N* are indicated in Figure 2 and Rsk σ ∆ is the characteristic stress range.
To determine the fatigue strength of bent bars with mandrel diameter "D" lower than 25ϕ s , where ϕ s is the bar diameter, the model code indicates multiplication of straight bar values by a reduction factor ξ = 0.35 + 0.026D/ϕ s , also proposed by EN 1992-1-1 [16]. The standards, however, do not make clear if the values obtained with the reduction factor should be used for stirrups.
According to Plos et al. [17], the fatigue strength of bent bars is reduced because the steel undergoes plastic deformation through bending, leading to microcracks. The smaller the bending radius of curvature, the smaller is the fatigue strength. The report adds that bending of stirrups is an exception since tests revealed that, in most cases, stirrups fail in the straight part, due to shear; thus, outside the bent zone. This can also be verified in Higgins et al. [18], [19] and Bachman et al. [20]. Thus, Souza et al. [21] indicate that straight bar S-N curves can be used for stirrups. Hillebrand and Hegger [22] and Hillebrand et al. [23] conducted experimental tests on prestressed concrete beams and also verified that stirrups fail in the straight part.
NBR 6118 [15] considers a reduction factor to determine the fatigue strength of stirrups. Because of this, the design of stirrups for bridge girders, following the Brazilian standard, is usually governed by fatigue. This is not the case for the design of longitudinal reinforcements. The standard considers the stress related to 2x10 6 cycles for fatigue design verification. For straight bars (reinforcing steel) and stirrups, with diameters up to 16 mm, these values are Δ fsd,fad = 190 MPa and Δ fsd,fad = 85 MPa, respectively. The value for straight bars is in accordance with the S-N curve from fib [14] (194 MPa for 2x10 6 cycles). The value of 85 MPa corresponds to the application of the reduction factor ξ with D = 3ϕ s to the straight bar value (194 x 0.428 ≈ 85). For the fatigue design verification, the partial safety factor for steel is 1.0 (characteristic stress).

Fatigue Damage
For variable stress amplitudes, the cumulative fatigue damage calculation may be performed according to the Palmgren-Miner rule, as in Equation 2. In this linear damage rule the failure occurs when the damage reaches a value known as Miner damage at failure (DM).
where i q : number of cycles obtained from the load spectra for each stress range amplitude; i N : number of cycles relative to the failure for each stress range amplitude (obtained from the S-N curve).
Theoretically, Miner damage at failure (DM) should be equal to one. However, as it is an empirical rule, critical damage observed in practice is a random variable. For deterministic analysis, DM = 1 is generally considered.

Fatigue service-life estimation
Considering the fatigue damage related to a one-year load spectra and DM = 1 in Equation 2 (deterministic analysis), the fatigue service life, in years, may be estimated based on Equation 3. for one year, where i q and i N are presented in Section 2.3.

Probabilistic analysis for fatigue
Crespo-Minguillón and Casas [6] consider the Weibull distribution to represent fatigue strength (number of cycles to failure) randomness. The expressions for the probability density function, . , n n n n f n exp n n u n u n u n where α : shape parameter; u : scale parameter; 0 n : lower limit (this value is always 0.0 for strictly positive variables, such as the fatigue strength).
Considering 0 n = 0, the relationship between moments (mean and variance) and parameters is given by Equations 6 and 7, where Γ() represents the Gamma function.
For N µ e 2 N σ determined experimentally, the shape parameter α can be determined iteratively by Equation 8, as reported by Melchers and Beck [24]. With the value of α , the scale parameter u may be determined from Equation 6.
For the fatigue safety assessment, based on Palmgren-Miner rule, the limit state function, as mentioned by Crespo-Minguillón and Casas [6], is given by: (9) where X : vector of random variables involved (the random variables are indicated in Table 6, as described in Sections 4 and 5); DM : Miner damage at failure, which is the resistance variable; 1/N i : elementary damage due to each stress range cycle from the load spectra; ∑(1/N i ): damage due to the load spectra for the reference period considered (50 years or 100 years, normally, for the design service life of bridges), which is the load effect variable.
According to Weibull distribution, Crespo-Minguillón and Casas [6] present S-N curves for reinforcing steel (straight bars) and prestressing steel (posttensioning curved tendons) for a 50% confidence level, as indicated in Table  2. In this paper, several simulations of the random variables are performed and different values of ∑(1/N i ) are calculated using the S-N curves from Table 2. Then, the probability distribution, moments and parameters of the dependent variable ∑(1/N i ) are obtained.
As reported by Crespo-Minguillón and Casas [6], the randomness of DM comes from the randomness of S-N curves. The authors present parameters for Miner damage variable in tests characterized by constant-amplitude stress cycles for reinforcing steel (straight bars) and prestressing steel (posttensioning curved tendons) for different S-N intervals (DM i in Table 3).  For a spectrum of variable amplitude stress cycles, the parameters of DM for the limit state function (9) can be obtained with the weighting of DM i related to the damage for each interval of stress range, as reported by Crespo-Minguillón and Casas [6] and Rodrigues et al. [7]. Therefore, in order to acquire the parameters of DM for (9), related to all stress ranges, the following procedure can be applied: • For each simulation of the random variables (Table 6), the percentages of ∑(1/N i ) corresponding to each interval of stress range (limits in Tables 2 and 3) are obtained; • At the end of all simulations, the average percentages of ∑(1/N i ) corresponding to each interval for all simulations are obtained; • By simulation, for each interval, several DM i values are generated, compatible with the respective weighting (average percentage); • With the DM i values for all the intervals, the probabilistic distribution and the parameters of the DM variable are evaluated.
In this paper, the simulations are performed using the Latin Hypercube Sampling -LHS simulation technique, which permit to obtain reasonable results with a reduced number of simulations. With the parameters of the variables of the limit state function, the FORM method is used, within StRAnD (Beck [31]), to calculate the reliability indexes. The LHS and FORM can be verified in Melchers and Beck [24] and Nowak and Collins [32].

Description of the system and WIM station
Regarding the relevance of actual traffic, this work uses vehicle records obtained from a high-speed weigh-inmotion (HS-WIM) station, as shown in Figure 3. The system is installed on road lanes, and vehicles are registered with no need to stop or to lower their speeds. In general, the system consists of lines of piezoelectric sensors, inductive loops, temperature sensors and a device to collect and analyze the records. Inductive loops detect vehicles, measure the distance between axles and speed, while piezoelectric sensors are responsible for weighting.
In this work, the data from the BR-381 (Fernão Dias Highway, which connects the states of São Paulo and Minas Gerais) station is used (same data used in Carneiro et al. [4]). The system is installed in two road lanes in the same direction (Minas Gerais to São Paulo), as indicated in Figure 3c. For each vehicle that passes over the sensors, the system provides information such as day, hour, lane, speed, total length, total weight, axle spacing and axle weight. The system was installed in July 2015 and the station remained in operation until 2019. The HS-WIM system avoids evasion, as drivers do not notice it. In addition, due to the continuous process of obtaining data (24 hours a day, on consecutive days), the system also yields the real frequency of vehicle occurrences, which is essential for fatigue analysis. These particularities give the technology a great advantage over static weighing stations. In the static weighing process, overloaded vehicles can trace alternative routes, the multiple presence occurrences are unknown, and the data is not collected 24 hours a day, in days and months in a row.  [4], [5]; b) weighting sensor before grout in BR-381 [4], [5], [33]; c) final appearance of the pavement in BR-381 [5], [33] (position sensors, which detect when vehicles change lanes, was not installed).
As WIM sensors are strongly influenced by temperature variation, the system must undergo periodic calibrations. In these measurements, a truck of known weight passes through the sensors at varying speeds, at different times of the day. From July 2015 to August 2017, five calibrations were conducted on the WIM system in BR-381 (July 2015, October 2015, February 2016, February 2017 and May 2017). The error, in relation to the total weight of the trucks, is around 10%. In Brazil, considering the reduced number of WIM stations, the use of this system is not common.

Filtering process and traffic characteristics
Even with calibrations, the system might present incorrect measurements: weights greater than the maximum traction capacity, or smaller than the self-weight, are some examples. Furthermore, several data are unnecessary for the study of live load for bridges. Passenger vehicles and light trucks, for instance, can be neglected, in most analyses, due to irrelevant load effect on the bridge. Thus, WIM data needs to be filtered before it can be used.
The filtering criteria need to be defined considering the characteristics of the Brazilian fleet and may vary according to the analysis need, in other words ULS, SLS or fatigue. Among the 14 criteria used in Wassef et al. [8], one filter eliminates trucks with total weight less than 90 kN. As reported by Laranjeiras [34], trucks weighing less than 70 kN accumulate fatigue damage that is irrelevant to the safety of the structure, even if they occur many times.
The following list presents, in order, the filters that were applied to the station considered. The filters were set after a research of National Department of Transportation resolutions (DNIT [35]- [37]) and truck manufacturers catalogs. Vehicles that fit any of the filters have been excluded. 1-GVW ≤ 62 kN (GVW is the total gross vehicle weight). This limit is set based on the maximum physical limit of GVW of lightweight two axle truck (veículo semi-leve in Portuguese). This criterion has the objective of removing irrelevant light trucks. 2-Pi ≤ 22 kN, where "Pi" is the axle weight. This value represents the steering axle load of maximum GVW of lightweight two axle trucks. 3-Pd > 320 kN, where "Pd" is the weight referring to the double tandem. This limit is set based on the heaviest axle weight of a double tandem found in the catalogs of truck manufactures in Brazil. This filter also contributes to eliminate records with GVW greater than the maximum traction capacity of the vehicle. For the sake of comparison, the legal limit for double tandem according to Brazilian regulations is 170 kN. 4-di ≤ 0.92 m, where "di" is the distance between axles. In Brazil, the minimum wheel diameter for a tandem was found to be 0.87 meters. Adding the wheel spacing of 0.05 meters, the final minimum distance between axles is 0.92 meters. 5-C > 36 m, where "C" is the total length of the vehicle. In Brazil, the maximum length allowed is 30 meters. On top of this value, an error margin was considered. 6-C > 15.4 m and GVW ≤ 104.3 kN. This filter screens out the unrealistic long and light trucks, and possible considerations of two trucks in following position as a single unit. The length is defined based on the 14-meter trucks (common in Brazil), with a margin of error of 10%. The GVW is defined based on the sum between 42.3 kN (axle with lower carrying capacity of a 14-meter truck) and 62 kN (minimum GVW for the other truck according to Filter 1). 7-Pi > 180 kN. Among all catalogs that were investigated, the maximum carrying capacity (physical limit) of a single axle was found to be 180 kN. This filter also contributes to eliminate records with GVW greater than the maximum traction capacity of the vehicle. For the sake of comparison, Brazilian legislation established the weight limit of single axles in 100 kN. 8-GVW ≥ 1.1∑Pi or GVW ≤ 0.9∑Pi, where "∑Pi" represents the sum of the weights of the axles. This criterion considers a margin of error of 10% on total weights. 9-∑di > C, where "∑di" represents the sum of axle spacings. This filter eliminates incorrect data of the wheelbase or the truck length. 10-C < 5 m. This filter eliminates incorrect truck length data. 11-V > 170 Km/h, where "V" is the vehicle speed. As in Brazil there are records of trucks traveling at 160 Km/h, the value of 170 Km/h was set as a limit. It should be mentioned that the occurrence of such events is rare. 12-P1 > 100 kN, where "P1" is the weight of the first axle (front). According to truck manufacturers catalogs this value is the maximum carrying capacity of the first axle. 13-GVW > 1500 kN. This limit is set based on the maximum carrying capacity found in catalogs, which corresponds to a nine-axle truck.
To determine the traffic load effect, WIM station data from September 2016 to May 2017 (273 days), is taken into account. It was found that use of a longer period of data does not significantly change the results.
Before filtering, the average daily traffic is 13292. After filtering, only 27% of the entire sample is considered as valid trucks. The first and second filters are responsible for filtering out 68% and 3.5% of the records, respectively. The other filters together eliminate only 1.5% of the records. It is worth mentioning that Filter 8 did not exclude any records. This indicates that the measurements of weights (total and individual per axle) are consistent. Filter eleven removes only one truck per month on average. Filters 3 and 7 eliminate less than one truck per day on average (0,3 for Filter 3 and 0,8 for Filter 7).
In this paper, the records are analyzed using Microsoft Excel spreadsheets. After filtering, the Average Daily Truck Traffic (ADTT) is 3655, where 83% of trucks are on the right lane (slow lane). The ADTT varies from 1079 to 5103. Table 4 shows some statistics for twelve most frequent classes, which represents 95% of the data after filtering, and the respective legal weights (5% of tolerance). It is worth pointing out, nevertheless, that a possible reduction in the values of Filters 1 and 2 would change the statistics of two-axle truck 2C (DNIT class). The maximum GVW obtained are higher than the limits imposed by Brazilian law.
The maximum GVW obtained was 1193 kN, corresponding to nine-axle truck 3M6 (DNIT class). Figure 4 shows the histogram for GVW for class 3M6 and Figures 5 to 7 show the histogram for weights for first axle, double tandem and triple tandem, respectively, for class 3M6. Figure 8 shows the histogram for weights for a single axle for class 3I3 (DNIT class), which is the heaviest class with single axles in Table 4. Based on the histograms presented, it is noted that the maximum weights obtained do not reach the values of filters 3, 7, 12 or 13.
To assess the sensitivity of Filter 1 in relation to the safety analysis, the value of 62 kN was changed to 90 kN, as used in Wassef et al. [8]. It was found that the results presented in Section 5 did not change, that is, the value of 62 kN is shown to be sufficient and may even be increased to eliminate a larger number of trucks, unnecessary for the analysis.

Evaluated bridge
This work evaluates the fatigue safety of the bridge girders illustrated in Figure 9. The dimensions were obtained from a real bridge that has a simply supported span of 28 meters and crossbeams only on supports. The girders are designed based on the Brazilian Design Standards, considering the ULS, SLS and Fatigue Limit States. The same reinforcement detailing for all girders is considered, i.e., the greatest rebar areas are considered.
The transverse load distribution is performed according to the Fauchart's method, using the Ftool program (Martha [38]). The method, which can be verified in Stucchi and Skaf [39], is also used for the traffic load effect. For the pavement, a total load of 3.2 kN/m 2 is considered (5 cm of pavement with weight of 24 kN/m 3 and an additional load of 2 kN/m 2 for paving improvements). Two possibilities are considered for the design, one with limited prestressing and another with complete prestressing level. Prestressing losses are considered as 25%. The Additional Impact Coefficient (CIA, in Portuguese) of NBR 7188 [1] is not computed in the design, since the coefficient should be only applied to slabs and crossbeams in the vicinity of the joints. The stirrups design, about the ULS and fatigue, consider the combination of shear force with torsional moment. It is important to observe, however, that the maximum shear and torsion for design are not simultaneous. The section close to the support without web enlargement, indicated in Figure 10, even with shear smaller than the section adjacent to the support, show the greatest stress range for stirrups in fatigue verification and the largest calculated area of transverse reinforcement. The calculated reinforcement area for stirrups for ULS (section indicated in Figure 10) needed to be increased to meet fatigue assessment, that is, the stirrups design was controlled by fatigue. Flexural fatigue analysis did not affect the areas of longitudinal reinforcement.
The NBR 6118 [15] provides two models for transverse reinforcement design based on the truss model. In Model I, the strut angle assumes a fixed value of 45 degrees; whereas in Model II, this angle may vary from 30 to 45 degrees. Furthermore, the portion of the shear force resisted by complementary mechanisms of the truss model differs from one model to another. Based on the standards equations, it is noticeable that Model II is more refined, while Model I is simpler. For the studied bridge, Model II with an angle of 30 degrees led to smaller areas for transverse reinforcement (reduction of around 15%, comparing to Model I). Thus, Model II is adopted in design, to obtain the smallest possible reliability indexes, according to design prescriptions by NBR 6118 [15].
The girder dimensions are indicated in Figure 10, and the calculated reinforcement areas are shown in Table 5. The compressive strength of concrete at 28 days (fck) is 40 MPa for precast beams and 30 MPa for slabs. The yield stress of reinforcing steel (f yk ) is 500 MPa. The strength (f ptk ) and yield stress (f pyk ) of prestressing steel are 1900 and 1710 MPa, respectively. The modules of elasticity for reinforcing and prestressing steel are 210 and 200 GPa, respectively, as indicated by NBR 6118 [15]. The secant modulus of elasticity of concrete is given as a function of f ck , as indicated by NBR 6118 [15], and granite is considered as coarse aggregate. For complete prestressing, it was found that the longitudinal reinforcing would not be necessary. In this case, the ratio of 0.09% is used for reinforcing, which corresponds to 0.5ρ mín , where ρ mín is the minimum ratio of longitudinal reinforcement (reinforcing and prestressing) for beams recommended by NBR 6118 [15]. The value of 0.5ρ mín was obtained based on the minimum ratio for positive reinforcing in one-way prestressed slabs of NBR 6118 [15], as recommended by IBRACON [40].

Traffic load effect
For the transverse load distribution on the girders, the trucks are assumed to be in the center of the traffic lanes and the value of 2 meters is adopted for the transverse wheel spacing, as illustrated in Figure 9. For the longitudinal analysis, the bending moment is calculated at midspan, the shear force and the torsional moment are calculated adjacent to the support, and the influence lines were implemented in Microsoft Excel. To determine the shear force and the torsional moment at the critical section (Figure 10), the load effects obtained by the influence line are linearly reduced.
To consider dynamic amplification, the factor presented in Almeida et al. [41] is used. Using trucks with five and six axles, with a total weight of 450 kN, the authors performed dynamic monitoring on Brazilian bridges with spans ranging from 7.5 to 45 meters. Equation 10 presents the obtained dynamic application factor (DAF). For the considered span length, this coefficient varies between 1.11 and 1.37, according to the WIM vehicle speed. . .
The traffic passage along the bridge generates irregular load cycles, with variable frequency and amplitudes, which does not enable straight use of the Palmgren-Miner rule. Then, the Rainflow counting method is used, which makes it possible to obtain individual cycles. Before performing the Rainflow count, however, a computational routine was implemented to obtain the load effects, at the considered sections, due to the traffic passage along the span of the structure. The load effects are obtained both for the single-vehicle passage and multiple presences.
To identify the situation of single and multiple presences, a routine was implemented to the WIM data, in Microsoft Excel spreadsheets, that uses the following information: • truck travel lane; • truck speed; • the time vehicle goes over the sensor (resolution of 0.01 sec); • the overall length of the vehicle's axle group; • bridge span length (simply supported structures).
The single passages of the vehicle on the bridge, i.e., when there are no axles of other trucks at the same time on the bridge, present the occurrence of 85% for the considered span (15% of vehicles are in multiple presence). The "sideby-side" case, i.e., when the front axle spacing between two vehicles, in different lanes, is less than half the length of the axle group of the first vehicle, present the occurrence of 2.35% in relation to the total valid trucks. Hence, for every 100 trucks on average, one is next to the other.
Using the Rainflow counting procedure and calculating the safety levels, it was found that, except for the side-byside case, the other situations of multiple presences can be replaced, without loss to the analysis, by the individual passage of the vehicles. The results presented in Section 5 are the same when considering the different situations of multiple presences (15%) and when considering only the "side-by-side" case (2.35%) as multiple presence. This occurs because, except for the side-by-side case, the position of the second vehicle in the other cases of multiple presences, for the considered span, coincides with the lower ordinates of the influence lines. It also applies to the case of three or more vehicles on the bridge. It is essential to observe, however, that for longer span lengths or continuous bridges, other situations of multiple presences might be relevant.
Applying Rainflow and calculating the safety levels, it was also found that the method can be substituted, without loss to the analysis, by the maximum vehicles load effects at the studied sections. The results, presented in Section 5, are the same when using the Rainflow method and when considering only the maximum vehicle load effects in the considered sections. This occurs because, except for the cycles related to the maximum load effects, the other cycles obtained with Rainflow, accumulate irrelevant damages, since they present, in their majority, small load effects in magnitude. It is important to note, however, that for other bridge scenarios, such as continuous structures, for instance, the cycles obtained by Rainflow can be relevant.
With the traffic load effects, the stress range in the reinforcement and the damage for each girder of the bridge are calculated. Overall, the stress range in the reinforcement is determined according to Equation 11, where max σ and min σ are maximum and minimum stresses, respectively.
The equations for calculating the stress range in the reinforcement are presented in the following sections. Although the present paper considers only prestressed girders, the equations also apply to reinforced elements. In this case, it is enough to eliminate the parcels related to prestressing.

Stress range in prestressing steel
For the stress calculation in longitudinal reinforcement, fib [14] recommends the cracked section to be considered, if for any combination of loads at SLS, cracking occurs. For the girders of the studied bridge, there is cracking due to traffic. Applying Equation 11, the stress range in the prestressed reinforcement,   : maximum bending moment at midspan, given by M g +M q,WIM,max, where M g is the bending moment due to dead loads (M g1 + M g2 + M g3 ) and M qWIM , máx is the maximum bending moment due to traffic (with impact) for the considered period for the WIM data; M g1 : bending moment at midspan due to the self-weight of the precast beam; M g2 : bending moment at midspan due to the self-weight of the slab and barriers; M g3 : bending moment at midspan due to pavement;

Stress range in longitudinal reinforcing steel
According to fib [14], the effect of differences in bond behavior of prestressing and reinforcing steel must be considered by multiplying the stress in the reinforcing steel by the factor s η . Thus, the stress range calculation in the longitudinal reinforcing steel, (   In NBR 6118 [15], the same factor s η is recommended both for reinforcing steel and bonded prestressed steel. In this work, however, it is considered only for the reinforcing steel, according to fib [14].

Stress range in stirrups
Applying Equation 11, the stress range in the stirrups, , for the critical section ( Figure 10), is given by Equation 17. As the bridge only have crossbeams on the supports, in addition to the shear force, the torsional moment must also be considered.
where sw A : Cross-sectional area of the stirrups at the section under consideration; s : longitudinal spacing between the stirrups in the evaluated section; θ : inclination angle of the compression diagonals (struts) regarding the longitudinal axis of the girder (without correction, since the parcel tanθ makes de adjustment in the inclination for fatigue). This angle is 45 degrees for Truss Model I and might vary between 30 to 45 degrees for Truss Model II. In this paper, Model II with a 30-degree angle is considered.
V g : absolute value of the shear force at the section under consideration due to dead loads (V g = V g1 + V g2 + V g3 ); V g1 : absolute value of the shear force at the section under consideration due to the self-weight of the precast beam; V g2 : absolute value of the shear force at the section under consideration due to the self-weight of the slab and barriers; V g3 : absolute value of the shear force at the section under consideration due to the pavement; V q : absolute value of the shear force at the section under consideration due to each traffic cycle; V p : absolute value of the shear force at the section under consideration due the prestressing; ( )  Figure 10); V c : parcel of the shear force resisted by the complementary mechanisms of the truss (it is important to note that, as indicated by NBR 6118 [15], a 50% reduction is considered in V c to calculate V 1 and V 2 . In NBR 6118 [42], the reduction was not indicated for Truss Model II). V c = V c0 , for Model I in reinforced elements (not used in this paper); V c = V c1 , for Model II in reinforced elements (not used in this paper); V c = V c0 (1+M 0 /M s,máx ) ≤ 2V c0 , for Model I in prestressed elements (not used in this paper); V c = V c1 (1+M 0 /M s,máx ) ≤ 2V c1 , for Model II in prestressed elements; Fusco [43] explains that the ratio between the decompression moment (M 0 ) and the maximum moment in the region under analysis is a relative measure of the possible cracking degree of the member. Strictly speaking, the author elucidates that the decompression moment could be substituted by the cracking moment r M (Equation 14). Thereby, this work computes r M instead of M 0 and M s,máx = M g + M qWIM , máx , where M g is the nominal moment due to the dead loads (M g1 + M g2 + M g3 ) and M qWIM,máx , is the maximum moment due to traffic (with impact), all referring to midspan.
Regarding torsion, the variation of torsional moment T ∆ assumes an absolute value of T q, as long as there is no inversion in the direction between (T g + T q ) and T g. In the case of an inversion, T ∆ is given by the highest absolute value between (T g + T q ) e T g . T q : torsional moment at the section under consideration from the traffic, associated with V q ; T g : torsional moment at the section under consideration due to the dead loads (T g = T g1 + T g2 + T g3 ); T g1 : torsional moment at the section under consideration due to the self-weight of the precast beams; T g2 : torsional moment at the section under consideration due to the self-weight of the slab and barriers; T g3 : torsional moment at the section under consideration due to the pavement; (two-legged stirrup).

FATIGUE-SAFETY ASSESSMENT
In this paper, fatigue safety is assessed for different traffic volumes. For the studied WIM station, the Average Daily Truck Traffic (ADTT) for two lanes is 3655. It is worth noting that the ADTT computes only heavy vehicles, according to the filter criteria presented. For fatigue-safety assessment, Rodrigues et al. [7] used ADTT = 5000 and Crespo-Minguillón and Casas [6] used ADTT = 6000 (both for two traffic lanes). Wassef et al. [8] used ADTT = 5000. The present paper assesses fatigue safety for three ADTT: 2500, 5000 and 7500. The correction in the traffic volume, therefore, is made by multiplying the damage calculated for the WIM traffic by the ratio ADTT/3655.
All results presented in this Section correspond to bridge girders designed following NBR 6118 [15]. The reliability indexes and fatigue life estimates for stirrups are evaluated using S-N curves for straight bars, as justified by several studies (Section 2.2).

Fatigue service life estimation
For the damage calculation (Equation 2), related to fatigue service life estimation, the S-N curves from fib [14] are utilized, as indicated in Table 1  Based on the Palmgren-Miner model, Jacob and Labry [44] explain that the estimates below 50 years are not adequate for a bridge project, estimates between 50 and a few hundred years can be considered as acceptable or good and estimates greater than 1000 years indicate unlimited fatigue life. These indicators, however, are only preliminary, since the conclusion about fatigue safety should be based, mainly, on the reliability indexes (probabilistic analysis performed in the next Section). This because even if the analysis considers the real spectrum of vehicles traveling over a bridge, the fatigue life calculation is based on deterministic values (characteristic S-N curves; characteristic concrete strength; DM = 1; nominal values for dead loads, reinforcement areas and elements dimensions), which, in reality, present variability.
The fatigue life estimation was performed for all girders. The estimates of longitudinal reinforcements (midspan section) and stirrups (section indicated in Figure 10) of all girders, for all ADTT considered, were found to be higher than 1000 years for both limited and complete prestressing, which indicates infinite fatigue life.

Reliability indexes
In literature, different values are observed for the target reliability index ( t β ) for fatigue. Aspects such as consequences of failures, inspection, repair possibility or reference period may influence the choice of t β . For a reference period of 50 years, fib [45] proposes t β = 3.1 for the calibration of the partial safety factors for fatigue. For the same reference period, ISO 13822 [46] indicates t β = 2.3 in the case of the possibility of inspecting the member subjected to fatigue and t β = 3.1, in the case of no inspection possibility. Thus, the reinforcement in the concrete structures fits in t β = 3.1, as stated by fib [45]. It is worth to mention that fib [45] recommends t β = 3.8 for ULS, corresponding to medium consequences of failure, and t β = 1.5 for SLS, corresponding to irreversible service failure. According to the model code, the value of t β = 3.1 for fatigue corresponds to the ULS assessment in the case of low consequence of failure. For the unlimited fatigue life approach, which is not the case of the present paper, Wassef et al. [8] indicates t β = 1.0.
The reliability indexes of the girders are calculated by evaluation of the limit state function, according to Equation 9, for the longitudinal reinforcement (midspan section) and stirrups (section indicated in Figure 10), for the design service life of 50 years (compatible with the reference period of t β from fib [45]). The reliability indexes are calculated for the girders that presented the highest values of ∑(1/N i ), related to the longitudinal and transverse reinforcement.
The random variables are presented in Table 6, where some of them are Brazilian statistics, as reported by Santiago [25], [26] and Santiago et al. [47], [48]. Other variables consider the statistics utilized in Nowak [27] and Wassef et al. [8]. The statistics for the dead load, obtained from Nowak [27], are computed in this work directly in the loads, i.e., the span length is considered as deterministic. Variables from Section 4 which are not indicated in Table 6 are considered as deterministic. For each random variable, related to longitudinal or transverse reinforcement, 50 values are generated by using the Latin Hypercube Sampling -LHS simulation technique. Crespo-Minguillón and Casas [6] also considers this number of simulations using LHS. Table 6. Random variables (indicated in Figure 10 and Section 4). The damage relative to the WIM data is calculated for each simulation, according to the S-N curves indicated in Table 2, and extrapolated for the intended reference period (50 years) considering the correction of ADTT. With the calculated values, the probability distribution, moments, and parameters of the dependent variable ∑(1/N i ) are obtained, which represents the load effect variable in the limit state function (9). The interior girder showed the highest values of ∑(1/N i ) for the stirrups, and the exterior girder (near to rightmost lane) presented the highest values of ∑(1/N i ) for the longitudinal reinforcement.

Category
The damage calculation is performed using Microsoft Excel, while the simulations and parameters obtaining are performed in Matlab. It is found that the ∑(1/N i ) variable fit the Lognormal distribution. Using ADTT = 5000 and limited prestressing, for example, Figure 11 shows the 50 values of ∑(1/N i ) for stirrups of interior girder in Lognormal probability paper (moments and parameters are indicated in Table 9). After all simulations for limited and complete prestressing, for all ADTT values considered, it is verified that almost all calculated stress ranges are lower than 205 MPa for reinforcing steel (longitudinal and stirrups) and all stress ranges are lower than 165 MPa for prestressing steel. Therefore, the parameters of DM, which is the variable of resistance in the limit state function, can be obtained directly from Table 3 (Weibull distribution) and the procedure presented in Section 2.5, for weighting of DM i , does not need to be applied. Tables 7, 8 and 9 present moments for ∑(1/N i ) and DM according to reinforcement, ADTT and prestress level. Using FORM, and the StRAnD [31] software, the fatigue reliability indexes are obtained, as indicated in Tables 7, 8 and 9. The longitudinal reinforcement (prestressing and reinforcing) and stirrups have shown fatigue reliability indexes higher than the intended value for design ( t β = 3.1). As expected, the reliability indexes decrease with the increase in traffic volume. The increase in ADTT from 5000 to 7500, however, did not significantly change the results. Complete prestressing presents higher safety levels than limited prestressing for both longitudinal reinforcements and stirrups. Although the Brazilian live load model is not in accordance with the unlimited fatigue life approach, as stated by Carneiro et al. [4], the design criteria, which involve partial safety factors and stress limits for fatigue design, guaranteed satisfactory fatigue safety levels for the evaluated bridge. It is important to mention that according to EN 1992-2 [49], a fatigue design verification is generally not necessary for prestressing and longitudinal reinforcing steel, in regions where, under the frequent combination of actions (including prestressing) only compressive stresses occur at the extreme concrete fibres (like complete prestressing of NBR 6118 [15]). The high reliability indexes for complete prestressing from Tables 7 and 8 for longitudinal reinforcements confirm the EN 1992-2 [49] consideration.

CONCLUSIONS
This work assessed the fatigue safety level provided by Brazilian design standards for a concrete highway bridge, using weigh-in-motion (WIM) data of an important federal Brazilian highway. The fatigue service life and fatigue reliability indexes for design life (50 years) were evaluated for prestressed girders. Using limited and complete prestressing levels, different traffic volumes were considered.
It was found that multiple presence trucks can be safely replaced by the analysis of individual trucks, except for side-by-side truck case. In addition, the Rainflow cycle method can be replaced by the maximum vehicles load effects. The design of longitudinal reinforcements (reinforcing and prestressing) and stirrups according to Brazilian codes ensured fatigue reliability indexes higher than values recommended by fib [45], that is, given the design conditions, it is unlikely that these reinforcements will fail due to fatigue. All fatigue life estimates exceeded 1000 years, which indicates unlimited fatigue life. Although the Brazilian live load model does not comply with the unlimited fatigue life approach, as stated by Carneiro et al. [4], the partial safety factors and stress limits for fatigue design guaranteed satisfactory fatigue safety levels for the evaluated bridge. It is important that the methodology presented in this paper be applied to a greater number of concrete bridges, to provide more general assessment of fatigue safety levels provided by Brazilian design standards.