Finite Element Model-Based Weight-Over Process Philosophy for Bridge Loading Capacity Evaluation and Rating Factor Estimation

Existing rating methods estimate bridge loading capacity and demand from secondary actions due to live loads in the primary structural components. In these methods, uniaxial yielding stress is traditionally used to detect component capacity using either stress quantities or shear-moment actions to compute the capacity demand of the bridge. *ese approximations can lead to uncertainties in load capacity estimation. *is article presents the weight-over process (WOP), a novel computer-aided approach to bridge loading capacity evaluation based on tonnage and rating factor estimation. WOP is expected to capture different forms of failure in a more general manner than existing methods. In WOP, a bridge finite element model (FEM) is discretized into many sections and element sets, each containing a single material type, and each assigned a suitable 3D failure criterion. *en, factored gross vehicle weights (GVWs) are incrementally imposed on the bridge FEM with those predefined ultimate unfavored loading scenarios in a manner similar to proof load testing.WOP code runs nonlinear analysis at each increment until a stopping criterion is met. Two representative bridges were selected to confirmWOP’s feasibility and efficacy. *e results showed that WOP-predicted values at the interior girders were between those of the conventional AASHTO and the nondestructive testing (NDT) strain measurement methods. *at may put WOP in a favorable zone as a new method that is less conservative than AASHTO but more conservative than real NDT testing.


Introduction
Highway bridges are vital infrastructure that must be evaluated periodically to ensure their safe operation under traffic loading. Over their lifetime, bridges age, deteriorate, and sustain damage, ultimately becoming structurally deficient. e load-carrying capacity index is a comprehensive integrity assessment tool that engineers can use for bridge rating. e service live load that can be safely carried over a bridge is called the load rating (LR) and is expressed as a rating factor (RF). For a consistent index of load-carrying capacity, all bridges are rated by a standard set of heavy trucks or vehicular live loads (LLs). In the USA, the LR must be biennially executed based on American Association of State Highway Transportation Officials guidelines [1] and Federal Highway Administration rules [2]. e load and resistance factor rating (LRFR) method provides the RF formula for new and existing bridges as the ratio of specific LL resistance to design LL, including its dynamic effects and associated uncertainties. e RF quantity can be interpreted as the number of a particular LL that can safely cross the bridge.
ere are three conventional approaches to practical LR: AASHTO analytical LR, experimental LR through either diagnosis load testing (LT) or proof LT, and finite element model-(FEM-) based LR. AASHTO LR is based on a simplified one-dimensional (1D) beam-line analysis and its prescribed provisions to rate for shear and moment separately as indicated in AASHTO's Manual for Bridge Evaluation [3,4]. e experimental approach to LR uses either diagnosis LT to confirm the AASHTO analytical assumptions or proof LT to prove a bridge's capability to carry its full dead load plus some magnified LL based on in situ LT [5][6][7][8].
e third conventional approach to LR is an advanced, reliable, and nondestructive testing (NDT) method that uses a validated bridge FEM constructed through structural model updating techniques. e common FEM-based LR can computationally facilitate the rating process in which it transforms the bridge superstructure to several equivalent 1D composite beams in the background, similar to the conventional AASHTO LRFR method [9][10][11][12]. FEM-based LR is much more demanding than the common AASHTO method since it properly simulates 3D loading distribution inside composite members of reinforced concrete slabs as well as steel girders, consequently impacting the bridge loading capacity evaluation (BLCE) process [13]. FEM is also capable of considering the reserve capacity of the structure model up to 120 and 185% for simply supported and continuous bridges, respectively [14]. Conventionally, FEM-based LR through the allowable stress design estimates RF by comparing the allowable stress levels at different composite critical sections to those spots' values when a unique moving LL passes through a traffic lane. en, the factored allowable stress of material and stress due to the moving LL are compared to calculate RF. FEM-based LR through LRFR calculates the ultimate capacity when the steel girder and the reinforced concrete section above it reaches a complete plastic section (PS) in which all 1D stresses are equal to the yielding points. e FEM stress levels at critical locations are automatically integrated over the composite girder cross-section as well as the effective tributary deck width to estimate maximum demands during LL traveling on a bridge model for LR [15,16]. Although there are other ways to impose LL for BLCE and RF estimation, mostly based on analyst preference [17], both approaches deploy the uniaxial yielding criteria of each steel girder when it reaches the yielding points [11].
is article proposes a novel algorithm called the weightover process (WOP) to rate bridge structures based on a calibrated FEM. Instead of evaluating the capacity and demand of the bridge composite girder in terms of either internal stress levels or internal actions (moments and shear forces) due to LL, WOP uses direct gross vehicle weight (GVW) to estimate capacity and demand. Since WOP considers tonnage instead of a secondary mechanism, it is expected to capture various failure mechanisms more realistically than transferring LL to associated secondary actions. First, the bridge 3D FEM is simulated and segmented into many slices. Each FEM slice resides in a spatial station along the bridge span axis, and then each slice is discretized with many element sets. WOP is programmed so that all elements that attain a single material property reside in a single element set. Second, suitable 3D failure criteria based on material type are assigned to those sets to detect the bridge's ultimate capacity. ird, WOP directly imposes LL on a bridge FEM with several prescribed loading patterns called the ultimate unfavored loading scenarios (UULSs). Finally, nonlinear analysis is performed while each UULS is incrementally imposed on each traffic lane until a PS is initiated. is process goes over all UULSs for all traffic lanes step by step, and the minimum RF for each girder is declared as its actual RF.

Materials and Methods
is section presents the conventional AASHTO RF philosophy and essential materials for calculating RF under different loading conditions. It then introduces the proposed WOP to estimate RF based on GVW.

AASHTO Rating Factor Method.
ree AASHTO design philosophies-allowable stress rating (ASR), load factor rating (LFR), and LRFR-have been developed over recent decades. Regardless of the type of design philosophy, the RF of a bridge component can be evaluated according to AASHTO's Manual for Bridge Evaluation [4] as follows: where C is the capacity, DC is the dead load effect, DW is the wearing surface effect, P is the permanent load, LL is the live load, and IM is the dynamic load allowance. e load factors c DC , c DW , c P , and c LL represent dead load factor, wearing surface load factor, permanent load factor, and LL factor, respectively. All prescribed load and resistance factors can be found in the AASHTO code. For an inventory level with a high reliability index (β � 3.5), the LL coefficient is c LL � 1.75, and for an operating level with a lower reliability index (β � 2.5), c LL � 1.35. e dead load factor is c DC � 1.25, and the wearing surface load factor is c DW � 1.50 for the strength limit state criteria. e capacity C for strength limit state can be estimated as where ϕ c , ϕ s , ϕ, and R n are the condition factor, system factor, resistance factor, and nominal member resistance specified by the LRFD code, respectively [1]. e condition factor ϕ c provides a reduction in strength due to the uncertainty in member resistance ranging from 0.85 to 1.00. e system factor ϕ s and LRFD resistance factor ϕ are specified by code based on their limit states. Last but not least, the critical parameter in the RF formula is the dynamic allowance factor (IM). e IM indicates how much inertia could escalate the responses of the bridge model to truckinduced vibrations. e dynamic amplification factor (DAF) due to the LL speed on the bridge model can be expressed in terms of IM as follows: where S dyn is the peak value of the dynamic structural response and S sta is the static structural response caused by the same action(s) or under the same load combination situation, estimated based on the most crucial of the understudied bridge points. AASHTO prescribes using HL-93 design truck loading for IM estimation only. e AASHTO formula for IM is [1] 2 Advances in Civil Engineering where D E is the minimum depth of Earth cover above the structure (ft). e IM value is estimated based on bridge surface conditions as follows: 11% for fine, smooth asphalt; 22% for intermediate clearness; and 33% at most for a very bad surface. AASHTO suggests using IM � 0.33 to attain high reliability when no in situ experimental testing or numerical modeling is performed. While some have reported that IM values might be much higher than the AASHTO-prescribed ones, up to IM � 2.80, it depends on bridge type, structural geometry, road surface condition, and truck speed [18]. Engineering judgment and in situ testing using suitable design codes can help analysts pick a fair IM value, either from guidelines or from experimental setup during RF estimation [19]. e RF obtained by (1) can be utilized to estimate the total ultimate safe weight capacity (RT) of the bridge in tonnage based on the AASHTO [4] formula (6A.4.4.4-1) as follows: where GVW is the gross vehicle weight of the rating truck, or LL, weight. e essential parameters to estimate RF through (1) are the capacity and demand of individual bridge composite girders under LL that must be estimated with acceptable accuracy. AASHTO code implicitly rates bridge superstructures on a component basis, meaning that a superstructure is broken down into several composite girders and each is rated independently based on its associated tributary slab width [14]. ose capacity and demand terms are expressed as internal stress values in the ASR code or internal actions as shear and moment forces in the LFR and LRFR codes; their implementations are thoroughly described in the AASHTO code.

Weight-Over Process Method.
is section introduces WOP as a novel FEM-based LR method. e concept of WOP is to substitute the capacity and demand of a bridge (1) in terms of vehicular tonnage rather than secondary actions (stress, shear, moment).

Stage (1): Construct Calibrated FEM.
e first step is to construct a detailed FEM of the bridge that comprises all primary structural components and then calibrate it through structural model updating techniques. Abaqus ® was used to simulate the FEM of both bridge samples presented in this work and perform nonlinear static WOP [20]. e concrete damage plasticity model was used to detect both compressive crushing and tensile cracking at all concrete regions [21][22][23].
e steel material stress-strain behavior was modeled by the nonlinear curve based on Rasmussen's proposed relationship, which is superior to the conventional Ramberg-Osgood model [24]. Four steel alloys of 316L, Grade-30, Grade-42, and Grade-50 were used by the FEM, and their characteristics were chosen from the American Society for Testing and Materials (ASTM) A570 standards [25,26]. Full composite action between the concrete slab and girders was imposed by a surface-to-surface contact algorithm, while the multipoint constraint algorithm was used to connect bridge components. In both samples' FEMs, the concrete body was modeled with a 3D 8-node iso-parametric element (C3D8R) capable of modeling plastic deformation. e reinforcement rebars were modeled as a 3D wire truss element member (T3D2) plus embedded-region constraints inside the concrete medium. e embedded-region constraint can satisfy a perfect bond between the concrete medium and reinforcement members to accurately simulate a unified reinforcement concrete object. e steel cross beams and girders were modeled with eight nodal points, six degrees of freedom, and a reduced integration shell element (S8R).
e Surface-Based Contact element type along with transverse nonlinear springs (defined by the Abaqus Connector Sections option) was used to better simulate the boundary conditions (bearing systems) along both the gravitational and transverse directions [27,28]. e surfacebased interactive elements with both normal behavior ("hard contact") and tangential behavior (penalty method) were defined to express the gravitational interaction of the superstructure over the substructure. e mesh sensitivity analysis was carried out to obtain an acceptable mesh size in which all initial 25 natural frequencies varied by less than 1%. It is required to activate geometric, boundary condition, and material nonlinearities during FEM implementations for BLCE and the LR process. e overall 3D FEM of the first bridge was simulated by 82318 nodes, while the second sample attained 375522 nodes.
After simulating the initial bridge FEM based on the asbuilt blueprint, it is required to perform structural model updating to validate it; model updating techniques are used to make a fair representative of the real bridge from its in situ information [29][30][31][32]. For bridge structures, the most important parameters that must be updated are boundary conditions and material properties using many available advanced techniques [33][34][35][36]. Nonlinear material models must be incorporated to develop WOP, so it can capture the yielding of the bridge in the plastic region. Material nonlinearity, geometric nonlinearity, and large displacement effects must be activated to make the FEM a more realistic representative.
e second stage of WOP is divided into two major substages. First, the FEM of the bridge superstructure is segmented and meshed such that all primary structural components attain the same mesh size along the bridge span axis. For instance, in the case of a reinforced concrete slab and steel girders system, the FEM must be meshed with the same elemental length for rebars, concrete slab, and girders.
is step is essential to better capture the yielding occurrence of the composite girder for the LR process. erefore, the FEM of each composite girder is equally segmented along its span length (e.g., X-direction) into several composite sections (N s ), and each section contains reinforced concrete elements and girder elements as shown in Figure 1.

Advances in Civil Engineering
In each section, several types of bridge materials (N m ) are attributed to the concrete, rebar, and girder elements. Hence, there are (N s * N m ) element sets for each composite girder of the bridge FEM; each set attains spatial characteristics to identify the sets that reside at the same section but contain different material types. Eventually, (N s * N m * N G ) FEM element sets are created by WOP for all bridge composite girders; (N G ) is the number of bridge girders.
In the second substage, appropriate 3D failure criteria are assigned to each element set based on their material types. While the bridge span is being loaded based on the WOP algorithm, some elements at several sets start to pass their yielding points. However, the plastic capacity of a composite girder is reached when all fibers inside a section pass their yielding points and enter the plastic region, and a PS is initiated based on AASHTO's Manual for Bridge Evaluation [4]. Since WOP deals with the 3D FEM of a composite bridge with various materials, appropriate 3D failure criteria must be chosen rather than the conventional uniaxial stress-strain relationship for LR. In general, any suitable failure criteria can be deployed to execute WOP based on material types and engineering judgment [37,38]. Common failure surfaces for ductile and brittle materials are shown in Figure 2, in which 3D FEM stress fields can be compared to their associated uniaxial behaviors to detect PSs, since WOP automatically computes all structural components' interactions (axialshear-flexural) and compares appropriate equivalent stress to material yield stress. In this study, the Von Mises yield criteria are used to capture the 3D yielding of FEM elements of steel components. For resilient and flexible components made with metallic materials, the Von Mises stress, or equivalent stress (σ e ), is commonly used as the failure criteria [39]: where σ 1 , σ 2 , σ 3 , and σ y are the first, second, and third principal stresses and the uniaxial yield stress, respectively, and σ 1 ≥ σ 2 ≥ σ 3 . For structural components that are made with brittle materials such as concrete, the maximum principal stress (σ mps ) is used as the failure criteria hereafter. e maximum principal stress criteria state that yield occurs when the largest principal stress exceeds its correspondent uniaxial yield strength; it is expressed as [40,41] σ mps � where σ yc and σ yt are compressive strength and tensile strength, respectively. is means that the ultimate model capacity is when, at specific FEM section (s), all concrete element sets pass σ mps and all steel element sets pass σ e , as defined in equations (7) and (8), respectively. At this point, the PSs are initiated, and, consequently, the capacity can be calculated.

Stage (3): Ultimate Unfavored Loading Scenarios and
Imposing LLs. In the third stage, several UULSs are prescribed and then imposed sequentially on the bridge FEM to run a nonlinear analysis. e term UULS embodies the failure mechanism that could occur if an ultimate number of LL trucks is placed on the bridge at the most unfavored locations, which produce the highest internal stresses. It is expected that UULS will generate worst-case loading scenarios to quickly initiate yielding on the bridge. For each UULS, all traffic lines are iteratively loaded by LL until a plastic section is initiated as an indication of the ultimate loading capacity of the bridge. Whether the bridge is simply supported or continuous span, several UULSs (N u ) can be imposed on the bridge's FEM to generate the worst loading patterns in the primary structural components. For shear failure, the most unfavorable position is close to the bridge support, and for flexural maximum action, it is the bridge midspan [7,14,15,42]. ere are three UULSs (N u � 3) for a single-span simply supported bridge, as shown in Figure 3. e order in which LLs are imposed on the bridge is indicated by the numbers on the trucks in Figure 3. UULS-1 activates the flexural failure by imposing an LL from the midspan toward the bridge sides, while UULS-2 and UULS-3 activate the shear failure by imposing LLs from the sides toward the midspan. e number of traffic lanes on bridges can be determined by (W/3600) in millimeters, which divides the clear roadway width (W) curbs to a standard lane width [1]. At each traffic lane, an integer number of LL, called the geometric number (N g ), can be placed practically and is defined as follows: where L BL , L LL , and "floor" represent the bridge span length, LL vehicle length, and floor function, which rounds the "N g " value to the closest and smallest integer number. Since in the proof LT an integer number of LL can be imposed on a bridge and WOP replicates the proof LT numerically, "N g " must be an integer and not a decimal in order to accurately represent proof LT scenarios. is geometric number indicates how many LL can be practically imposed on FEM to replicate realistic loading conditions. In order to be conservative when imposing the legal load, it was assumed that the trucks were placed on each lane, bumper to bumper, without any spacing between them, similar to what happens during proof load testing. e way LL is imposed on FEM, illustrated in equation (10), is to artificially replicate proof LT in a numerical setup. For any bridge type, (N u ) numbers of UULSs exist. In each case, the order of imposing LLs is specified by a number associated with UULS-i, as shown in Figures 3-5. At each analysis step, one LL (GVW i ) is magnified by AASHTO factors and then added to the FEM at the ith spot, which is specified in UULS. For the ith UULS, the overall sequential loading pattern is implemented as follows:    (1) requires that bridge capacity be deducted by the associated resistance factors indicated in equation (2). To meet this criterion, FEM materials can be manipulated through WOP. e stressstrain (σ − ε) datasets are adjusted to the reduced stress levels by resistance factors (ϕ c ϕ s ϕ) as follows:

Stage (4): Quantify BLCE. Equation
Also, to estimate the remained capacity of the bridge in the equation (1) numerator, the dead loads, wearing surface loads, and permanent loads (c DC (DC) + c DW (DW) ± c P (P)) are also enlarged, based on prescribed AASHTO load factors, and then imposed on the bridge FEM. e following steps will be taken to implement W(i) on the FEM. As described in stage (3), there are (N u ) numbers of UULSs, depending on whether the bridge is single-span simply supported or continuous span. From each UULS, one critical weight W cr (i) can be estimated by summing all loaded FEM nodal forces of "m" nodes, which are located at "n" pivotal lanes (F n m ) at the moment of PS initiation as follows: where N i t and N i n are the number of loaded traffic lanes and the number of loaded FEM nodes for the ith UULS, respectively. From those (N u ) cases, the primary critical weight (W P cr ) is the minimum weight at which the bridge collapses, as follows: So far, the remained capacity of the bridge in weight/ tonnage quantity has been estimated by (W P cr ), while appropriate AASHTO factors were inserted during WOP implementation to be consistent with the numerator of equation (1). As indicated in equation (1), the vehicular LL capacity and demand are stress quantity in the ASR method or shear-moment actions in the LFR and LRFR methods. In order to attain unit consistency, the demand of vehicular LL in the WOP method is the LL weight; thus, the LL of the denominator of equation (1) must be substituted with a single GVW that has the same unit (tonnage) as the remained capacity W P cr . erefore, the FEM-based RF based on the proposed WOP can be evaluated as follows:

Advances in Civil Engineering
e WOP algorithm is depicted in a flowchart in Figure 7 that shows the stages one by one.

First Bridge Sample.
To deploy and validate the WOP for a real bridge, a sample bridge in Iowa (FHWA #31690) was selected. As-built blueprints were received from the Iowa Department of Transportation and used to construct the FEM of the bridge. It is a single-span simply supported composite bridge with a reinforced concrete deck, steel plate girders, and cross beams, as shown in Figure 8(a). e bridge has a length L Bl � 18.5(m), a span width L Bw � 9.15(m), two traffic lanes N tl � 2, four girders (N G � 4), and a geometric number N g � 2. e UULS for a single-span bridge is shown in Figure 3 for the most severe load configuration and truck position. e FEM of the first sample and its traffic lane loaded by HL-93 LL are shown in Figures 8(b) and 8(c).   Figure 9. Due to the sample symmetry, two PSs were generated, PS1 and PS2. Although the plastic hinge theoretically happens exactly in the middle of a homogenous girder under flexural failure, the first PS occurs at stations of 0.38 * L Bl and 0.62 * L Bl , spaced 12% from the bridge's middle station due to the composite action of the superstructure, various material yielding behaviors, and 3D modeling.

Sensitivity Analysis.
It is common to deploy FEM in order to perform sensitivity analysis to confirm experimental results or potential outcomes of severe incidents such as damage scenarios, strengthening of components, and corrosion [43] erefore, the four following sensitivity analysis cases are defined and analyzed to compare BLCE and RF. In Figure 10, the force displacement of the WOP nonlinear analysis is shown; the force is the total reaction force of the bridge FEM toward the midspan deflection. Figure 11 shows the estimated RFs for the various cases as well.
(1) Rebar Reduction Strength. e first adjustment in the FEM was to reduce the area of the reinforcement rebars by 50%; results are indicated in the plots by "Rebar Reduction Case." Rebar corrosion happens in reinforced concrete structures due to cover cracking, spalling, and degradation of the outer concrete layers; this causes a reduction in the rebar cross-sectional area at high levels of corrosion [43]. To investigate the rebar strength for the BLCE, only the reinforcement area properties were adjusted; all other parameters were kept constant with baseline FEM.
(2) Increasing Concrete Strength. e second case study increased the concrete strength to investigate its effects on the BLCE and RF of the bridge structures; results are indicated in the plots by "Enhanced Concrete Case." For this sensitivity study, the concrete strength for the bridge deck was enhanced from 20 MPa to 40 MPa in the FEM code to determine its effects on BLCE. It is assumed that all other structural system characteristics remained constant.
(3) Increasing Steel Girder Strength. e third case study increased the steel girder strength to investigate its effects on the BLCE and RF of bridge structures; results are indicated in the plots by "Enhanced Steel Case." For this sensitivity study, the steel strength for the bridge girders was enhanced from steel 316L to steel Grade-42 in the FEM code to determine its effects on BLCE. Again, it is assumed that all other structural system characteristics remained constant.

(4) Settlement Sensitivity Analysis.
e last sensitivity study simulated settlement modeling for one of the interior girders to investigate its effects on BLCE and RF. Results are indicated in the plots by "Single Settlement Case." is is among the most prevalent damage in a bridge structure that happens immediately after construction or during its lifetime. e deflection of bridges is limited to the span length divided by 800 (L/800) for simple and continuous spans and L/300 for cantilever arms [4]. To simulate a severe settlement scenario above the allowable deflection limits, a settlement of δ � L Bl /200 � 9.2(cm) was considered on one side of an interior girder to do a sensitivity study.

Second Bridge Sample.
e second sample is the Powder Mill Bridge (PMB), a three-span continuous composite concrete deck with steel girders located in Massachusetts. e exterior span length is 11.75 m, the interior span is 23.5 m, there are six girders (N G � 6), and the overall bridge length is 47 m, as shown in Figure 12(a).
ere were two main reasons for selecting PMB as the second sample. e first was to test WOP for a continuous bridge model. e second was that its information and experimental RFs are available in the literature [43][44][45][46]. In these references, the three RFs for each bridge girder are the AASHTO conventional LR approach, LR modified by NDT data, and conventional FEM-based LR. e NDT method uses field measurements to adjust and scale AASHTO-proposed distribution factors and, consequently, estimate RFs [43]. is type of NDT data is normally used to calibrate the girder distribution factor and the actual girder loading capabilities based on in situ information from strain measurements. While the NDT result is not perfect, it is considered the most trusted and authentic outcome since it illuminates as-is and real structural behavior in its operational conditions. e FEM of the PMB was remodeled in Abaqus by this article's authors based on detailed geometric information and material strengths presented in the references (see Figure 12(b)). Fourteen UULSs were defined by the authors to impose potential severe LL positioning. e UULSs were defined from three-span continuous influence line concepts; they either generate maximum shear internal actions at girder supports (or bridge pier cap spots) or maximum S, Mises (Avg: 75%) +2.702e+08 +2.200e+08 +2.163e+08 +2.125e+08 +2.088e+08 +2.050e+08 +2.013e+08 +1.975e+08 +1.938e+08 +1.900e+08 +1.863e+08 +1.825e+08 +1.788e+08 +1.750e+08 +1.764e+04 Y X Z S, Min. Principal (Avg: 75%) +4.165e+08 +1.985e+08 +1.948e+08 +1.910e+08 +1.873e+08 +1.836e+08 +1.798e+08 +1.761e+08 +1.723e+08 +1.686e+08 +1.649e+08 +1.611e+08 +1.574e+08 +1.537e+08 +3.110e+06 flexural moment at the midspan of each girder span (s). In general, the overall UULSs must contain all potential severe LL placements, which have the highest internal effects on the bridge's primary structural components. e overall UULSs are separated based on their potential failure modes: shear UULS, shown in Figure 4, and flexural UULS, shown in Figure 5. Rating factors estimated by WOP and those reported in the literature are compared in Figure 13.

Recommendation for Practical Implementation.
To successfully implement the WOP for real-world bridge LR, the first and essential practical step is to perform the structural model updating process to construct a fair representation of the bridge by its FEM as explained in the WOP first stage. Although, in general, any available robust model updating technique can be exploited to perform WOP's first stage, detailed descriptions of the model

12
Advances in Civil Engineering updating techniques used in this article are available in the references [29,30]. After simulating the initial FEM based on the as-built blueprints of a bridge, the analyst must calibrate its FEM either with in situ NDT data or visual inspections based on AASHTO guidelines for field inspection. With the exception of the first WOP step, all other stages can be executed with the developed computer program that imposes the potential UULS and declares RFs with higher accuracy.

Conclusions
is article proposes a methodology called the weight-over process (WOP) to evaluate a bridge's ultimate capacity and its RFs based on actual GVW in excessive traffic scenarios. By considering rightful GVW instead of secondary action estimations, WOP could capture various failure mechanisms with higher accuracy than other conventional approaches.
e WOP process examines all traffic lanes to find the most critical RF for each girder and then declares the minimum RF for each girder to be the dominant RF. With WOP, there is no need to manually find the RF in the negative and positive sections; WOP incrementally places the LLs at the critical spots defined as UULS. Bridge failure mode in the case of excessive traffic loading is also identified through WOP and can be used for the rehabilitation plan, depending on the type of failure. It is well known that all FEM-based LR approaches return higher RFs than the AASHTO method because they consider the reserved material capacity and structural redundancies and also incorporate 3D composite action between different components accurately without any simplification.
For the first bridge sample, several sensitivity analyses were defined by changing the rebar area, concrete strength, steel girder strength, and settlement to find the most critical components in the BLCE and RFs. It was indicated that the settlement phenomenon is the most crucial cause of damage because it decreased RF tremendously relative to the other simulations. For the second bridge sample's exterior girders, WOP returns higher RFs than the experimental NDT method, while for interior girders it returns lower RFs. Except for the first girder, WOP results are much closer to the experimental NDT LR method than those of the conventional FEM-based method, as shown in Figure 13.
Except for small differences at girders 1 and 5, the results in Figure 13 show that WOP demonstrated better performance than the conventional FEM. e NDT data could be considered the best available reference for comparison and validation purposes with WOP; it compares theoretical strains with measured strains and scales the conventional AASHTO rating based on these differences [43,46]. At the interior girders, WOP-predicted values were between those of AASHTO and NDT, which may indicate that it is less conservative than AASHTO but more conservative than the real (NDT) testing.
As shown in Figure 13, WOP returns higher RFs than the AASHTO and NDT LR methods for the exterior girders and lower RFs than the NDT LR method for the interior girders.
is could be related to the nature of the NDT data used in this work. e NDT used measured strain to enhance the conventional AASHTO rating, and it seems that the external girders were not sufficiently stressed during the testing [43]. at could be why WOP presents higher values than NDT. Another issue could be the lack of information about sidewalk barriers and guardrails; their equivalent dead loads were assumed and then approximately imposed by the authors. e lower rating by WOP for the interior girders could be related to how the strain data was scaled for the NDT method or to the possibility that the 3D failure criteria and stress field could capture secondary stress fields during PS capacity estimation more effectively than the conventional 1D approaches. Except for girders 1 and 5, WOP returns were closer to the experimental RFs than those of the conventional FEM-based method. In the case of girders 1 and 5, the deviation from the real experimental RF is very small. Also, the bridge FEM in the second example was resimulated by the authors based on information provided in the references. e main computational cost of WOP is attributed to its implementation for different types of bridges. Difficulties such as the FEM segmentation, UULS definition, and sequential analysis steps could be valid concerns for program developers. Commercial software must be reprogrammed and upgraded to automatically prescribe UULSs from its library and impose sequential LLs on the FE model. Also, an FE model segmentation stage is needed to slice the model along its pivotal axis to myriad spatial stations/elements that require swift processors in practice.

Data Availability
All the datasets and the finite element code used in this article are available via the Data On-Request option. is dataset can be obtained by contacting the first author, Dr. Ali Karimpour (ali-karimpour@uiowa.edu). e experimental dataset for the second example is also available in references [43,46]. Disclosure e contents reflect the views of the authors, who are responsible for the facts and the accuracy of the information presented herein and are not necessarily representative of the sponsoring agencies.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.