 Open access peer-reviewed chapter - ONLINE FIRST

# Perspective Chapter: Simplified Matrix Calculation for Analysis of Girder-Deck Bridge Systems

Written By

Alvaro Gaute-Alonso and David Garcia-Sanchez

Submitted: December 17th, 2021 Reviewed: December 22nd, 2021 Published: February 10th, 2022

DOI: 10.5772/intechopen.102362

From the Edited Volume

## Applied Methods in Design and Construction of Bridges, Highways and Roads - Theory and Practice [Working Title]

Dr. Khaled Ghaedi

Chapter metrics overview

View Full Metrics

## Abstract

In the design of girder-deck bridge systems, it is necessary to determine the cross-sectional distribution of live loads between the different beams that make up the cross section of the deck. This article introduces a novel method that allows calculating the cross-sectional distribution of live loads on beam decks by applying a matrix formulation that reduces the structural problem to 2 degrees of freedom for each beam: the deflection and the rotation of the deck slab at the center of the beam’s span. To demonstrate the proposed method, the procedures are given through three different examples by applying loads to a bridge model. Deflection, bending moment, and shear force of the bridge girders are calculated and discussed through the given examples. The use of the proposed novel method of analysis will result in significant savings in material resources and computing time and contributes in the minimization of total costs, and it contributes in the smart modeling process for girder bridge behavior analysis allowing to feed a bridge digital twin (DT) model based on Inverse Modeling holding the latest updated information provided by distributed sensors. The presented methodology contributes also to speed up real-time decision support system (DSS) demands.

### Keywords

• girder bridge decks
• optimized matrix method
• structural grillage models

## 1. Introduction

Girder bridge decks are a structural typology commonly used in the design of road and railway bridges, and therefore, any optimization in their calculation has a depth impact on the project phase. As explained in , the design process relies solely on the designers’ experience, intuition, and ingenuity resulting in a depth cost material, time, and human effort. It is very common to use structural grillage models [2, 3, 4, 5, 6, 7, 8, 9, 10] to calculate the cross-sectional distribution of live loads between the different beams that make up the cross section of the deck. Another way to deal with the design and calculation of such decks would be to apply different formulations contained in the bridge design standards that allow approximating the cross-sectional distribution of the bending moment and shear stress caused by live loads through what is known as the load distribution factor “LDF” [11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23]. The “LDF” associated with each case study is conditioned by the type and number of beams, their spacing and length, as well as the existence or not of transverse diaphragms that bring transverse rigidity to the deck.

However, the LDF does not allow determining the distribution of bending stresses on all beams; therefore, the design is oversized. The need for a method to calculate the cross-sectional distribution across all types of beams, without resorting to complex structural grillage models or finite element models in specific structure calculation programs, or to approximate parametric methods based on the “LDF,” is one of the authors’ motivations for the development of the research work that has given rise to this article.

## 2. Traditional methods for girder bridge deck analysis

Structural grillage models began to be used for the analysis of cross-sectional distribution on beams in the 1960s. These models divide the beam deck into longitudinal and transverse beams (Figure 1). Longitudinal beams are responsible for providing the longitudinal bending stiffness of the deck, considering as many longitudinal beams as beams conform to the analyzed beam deck. The structural section of each of the longitudinal beams shall be the result of the section composed of the beam analyzed and the effective depth of the contributing deck with that beam [24, 25]. The cross-sectional distribution of the structural model is provided by the cross-beams and the torsional stiffness of the longitudinal beams. The structural section of the cross-beams corresponds to a rectangular section with an equivalent depth to the slab thickness and a depth according to the discretization used in the grillage models. The analysis of such structural models involves the use of specific structure calculation programs that provide the computing power necessary for the resolution of the proposed matrix problem. Figure 1.Grillage discretization: (a) girder bridge deck cross section; (b) structural grillage model.

The use of structural grillage models allows obtaining the structural response of girder bridge decks to live loads, adequately manifesting the cross-sectional distortion of the girder bridge deck and the distribution of stresses on each of the longitudinal beams that make up the deck. However, such models involve complex, time-consuming analysis that involves the need to use specific structure calculation programs, making it necessary to use simplified methods for the start of the design process.

The concept of LDF was first introduced using empirical formulas at the American Association of State Depthway Officials (AASHO) in 1931 ; these methods propose calculating the cross-sectional distribution on beam bridge decks roughly. The LDF is calculated from a series of formulations that parametrically treat the calculation of the percentage of bending moment and shear stress supported by the most requested longitudinal beam. The parameters that condition the calculation of the LDF are the depth of the beam, the span length, the spacing, the number of beams, the position of the load, and the beam position. This method provides an approximate value of the maximum bending stresses on the beams but is not able to reproduce the cross-sectional distribution of longitudinal bending between all the beams that make up the deck.

## 3. Proposed method for the study of the cross-sectional load distribution on a girder bridge deck

The authors propose the use of a method that allows obtaining the cross-sectional distribution of live loads on girder bridge decks without using empirical formulas of LDFs or complex structural models involving the use of specific structure calculation programs. The proposed method is based on using a virtual model that reflects the transverse stiffness of the slab deck, supported on a series of springs that provide the flexural stiffness Eq. (1) and torsional stiffness Eq. (2) of the longitudinal beams that make up the girder bridge deck (Figure 2). Figure 2.Proposed method model for cross-sectional distribution on a girder bridge deck.
Kv,i=48·EIiL3E1
Kt,i=2·GJiLE2

where EIi = longitudinal bending stiffness of beam “i”; GJi = longitudinal torsional stiffness of beam “i”; L = distance between bridge supports.

The proposed method considers 2 degrees of freedom for each longitudinal beam on the bridge deck: (1) the deflection and (2) the rotation of the deck slab at the center of the beam’s span. Figure 2 represents the structural model scheme for a girder bridge deck composed of three longitudinal beams. The matrix approach that solves the structural problem of the cross-sectional distribution of live loads between the different beams that make up the deck is raised in matrix Eqs. (3), (4), (5). The stiffness of the springs that represents the longitudinal bending and the torsional stiffness of the longitudinal beams that make up the model corresponds to the stiffness equivalent to the center of the span.

Ke=EIeLe3126Le6L4Le2126Le6Le2Le2126Le6Le2Le2126Le6Le4Le2=K11,eK12,eK21,eK22,eE3
KBi=kv,i00kt,iE4
P1Mt,1P2Mt,2P3Mt,3=K11,1+KB1K12,10000K21,1K22,1+K11,2+KB2K12,20000K21,2K22,2+KB3×fv,1θ1fv,2θ2fv,3θ3E5

Where EIe = transverse bending stiffness of element “e” of the upper deck slab; Le = length of element “e” of the upper deck slab; fv,i = vertical displacement experienced by longitudinal beam “i”; Ɵi = transverse rotation of the bridge deck over the longitudinal beam “i.”

### 3.1 Loads applied in the center of the span of the longitudinal beams

The maximum longitudinal bending stress in each of the longitudinal beams that make up the bridge deck corresponds to the application of a point load in the center of the beam span. The distribution of the maximum bending moment and maximum shear stress in the different longitudinal beams is obtained by Eqs. (6) and (7), respectively.

Mfmáx,n=Q·L4·fv,ni=1Nfv,iE6
Qmáx,n=Q·L2·fv,ni=1Nfv,iE7

where Q = point load value; L = girder bridge span length.

Example 3.1.1.

In the study of the structural behavior of the bridge deck represented in Figure 3, it is desired to know the distribution of the bending moment and the shear force in each of the longitudinal beams generated by the application of the following load states: (a) a point load of 300 KN in the center of the span of beam 1, and (b) a point load of 300 KN in the center of the span of beam 2. The girder bridge is made up of four longitudinal beams 1.05 meters deep (two end beams and two central beams) and an upper slab 0.25 meters thick and 16.83 meters wide. Considering the effective width of the upper slab in each of the longitudinal beams, the inertia to the longitudinal bending of the end and central beams takes the value of 0.1445 m4 and 0.1904 m4, respectively. The longitudinal torsional inertia of the longitudinal beams takes the value of 4.3·10−3 m4. The spacing between longitudinal beams is 5.13 meters, while the distance between the support devices of each longitudinal beam is 25 meters. Both the upper slab and the longitudinal beams are made of structural concrete whose modulus of elasticity reaches 35,000 MPa. Figure 3.Girder bridge deck: (a) side view of load states 1 and 2; (b) front view of load state 1; (c) front view of load state 2.

The stiffness of the springs on which the upper slab of the bridge deck rests, and which simulates the bending and torsional stiffness of the longitudinal beams (Figure 4), is calculated as follows:

End beams

Kv,end=48·EIendL3=48·3.5·107KNm2·0.1445m425m3=15,537KNm
Kt,end=2·GJendL=2·3.5·107KNm22·1+0.2·4.3·103m425m=5,017KN·m
KB,end=kv,end00kt,end=15,537KNm005,017KN·m

Central beams

Kv,c=48·EIendL3=48·3.5·107KNm2·0.1904m425m3=20,472KNm
Kt,c=2·GJendL=2·3.5·107KNm22·1+0.2·4.3·103m425m=5,017KN·m
KB,c=kv,end00kt,end=20,472KNm005,017KN·m

The mechanical characteristics of the beam elements that simulate the transverse distribution provided by the upper slab (Figure 4) are calculated as follows:

I=112·20m·0.25m3=2.6·102m4
L=5,13m
K=E·IL3126L6L4L2126L6L2L2126L6L2L2126L6L4L2
K=81,020KNm207,800KN207,800KN710,690KN·m81,020KNm207,800KN207,800KN355,340KN·m81,020KNm207,800KN207,800KN355,340KN·m81,020KNm207,800KN207,800KN710,690KN·m
K=K11K12K21K22

The global stiffness matrix used in the proposed method for the analysis of the transverse distribution of live loads in girder bridge decks is obtained as follows:

KG=K11+KB,endK12K21K22+K11+KB,c00000000K1200000000K2100000000K22+K11+KB,cK12K21K22+KB,end

The transversal distribution of each load state is obtained by planting the compatibility between loads and displacements:

300KN0000000=K11+KB,endK12K21K22+K11+KB,c00000000K1200000000K2100000000K22+K11+KB,cK12K21K22+KB,end×xfv,1θ1fv,2θ2fv,3θ3fv,4θ4

Deflections (Figure 5) and rotations

Maximum bending moments

Mfmáx,n=Q·L4·fv,ni=1Nfv,i=1,588KN·m499KN·m25KN·m187KN·m

Maximum shear forces

Qmáx,n=Q·L2·fv,ni=1Nfv,i=127KN40KN2KN15KN

00300KN00000=K11+KB,endK12K21K22+K11+KB,c00000000K1200000000K2100000000K22+K11+KB,cK12K21K22+KB,end×xfv,1θ1fv,2θ2fv,3θ3fv,4θ4

Deflections (Figure 6) and rotations

Maximum bending moments

Mfmáx,n=Q·L4·fv,ni=1Nfv,i=566KN·m859KN·m478KN·m28KN·m

Maximum shear forces

Qmáx,n=Q·L2·fv,ni=1Nfv,i=45KN69KN38KN2KN

### 3.2 Loads applied at a distance “x” from one of the supports of the longitudinal beams

The proposed method is applicable to calculate the structural response of a girder bridge deck to the application of a vertical load in any cross section. Using the Maxwell-Betti reciprocity theorem , for the application of a point load of value “Q” at a distance “x” from one of the two support points of a longitudinal beam, the deflection at the center of the span is obtained by applying the formulation reflected in Eq. (8). Likewise, the distribution of the maximum bending moment and maximum shear stress in the different longitudinal beams is obtained by Eqs. (9) and (10), respectively.

fcl,Qx,i=fv,i·sinπ·xLE8
Mfmáx,n=Q·Lx·xL·fv,ni=1Nfv,i/xL/2E9
Qmáx,n=Q·LxL·fv,ni=1Nfv,i/xL/2E10

Example 3.2.1.

In the study of the structural behavior of the bridge deck analyzed in Example 3.1.1., it is intended to know the deflection in the center of the span of the longitudinal beams, as well as the distribution of the maximum bending moment and the maximum shear force in each of the longitudinal beams generated by the application of the following load states: (a) a point load of 300 KN at a distance equivalent to L/3 from one of the supports of beam 1, and (b) a point load of 300 KN at a distance equivalent to L/3 from one of the supports of beam 2 (Figure 7). Figure 7.Point load of 300 KN at a distance equivalent to L/3 from one of the supports of a beam.

Deflection in the center of the spam of the longitudinal beams

fcl,Qx,i=fv,i·sinπ·xL=13.12mm4.12mm0.20mm1,55mm

Maximum bending moments

Mfmáx,n=Q·Lx·xL·fv,ni=1Nfv,i=1,412KN·m443KN·m22KN·m167KN·m

Maximum shear forces

Qmáx,n=Q·LxL·fv,ni=1Nfv,i=169KN53KN3KN20KN

Deflection in the center of the spam of the longitudinal beams

fcl,Qx,i=fv,i·sinπ·xL=4.12mm6.25mm3.48mm0.20mm

Maximum bending moments

Mfmáx,n=Q·Lx·xL·fv,ni=1Nfv,i=503KN·m764KN·m425KN·m25KN·m

Maximum shear forces

Qmáx,n=Q·LxL·fv,ni=1Nfv,i=60KN92KN51KN3KN

### 3.3 Loads applied at any point on the bridge deck

The proposed method is applicable to calculate the transverse response of girder bridge decks to the application of a vertical load at any point of the cross section. As in any matrix calculation, if the vertical load acts on a section of slab between longitudinal beams, the degrees of freedom of the structural model are locked and the reactions in the locked degrees of freedom are calculated (rigid step). Subsequently, the degrees of freedom are released and loaded with the reactions obtained in the previous step to obtain the vertical displacement of each of the longitudinal beams that make up the bridge deck (flexible step).

Example 3.3.1.

In the study of the structural behavior of the bridge deck analyzed in Example 3.1.1., it is intended to know the deflection in the center of the span of the longitudinal beams, as well as the distribution of the maximum bending moment and the maximum shear force in each of the longitudinal beams generated by the application of the load states described in Figure 8. Figure 8.Load state generated by the actuation of a four-axle truck 70 + 110 + 100 + 100 KN: (a) plan view; (b) front view.

Rigid step (Figure 9).

The load generated by the action of the truck axles on longitudinal beams 1 and 2 is obtained considering the compatibility between loads and movements in the upper slab between longitudinal beams 1 and 2.

VLB1MLB1Qi0Qi0VLB2MLB2=K11,1K12,1K21,1K22,1+K11,200000000K12,200000000K21,200000000K22,2+K11,3K12,3K21,3K22,3·00v2Ɵ2v3Ɵ300

Load on longitudinal beams 1 and 2

Q=1.46·QKN1.09·QKN·m0.54·QKN0.70·QKN·m

Flexible step (Figure 10).

1.46Q1.09Q0.54Q0.70Q0000=K11+KB,endK12K21K22+K11+KB,c00000000K1200000000K2100000000K22+K11+KB,cK12K21K22+KB,endxfv,1θ1fv,2θ2fv,3θ3fv,4θ4

Truck axle 1

Deflection in the center of the spam of the longitudinal beams

fcl,Qx,i=fv,i·sinπ·xL=2.39mm1.20mm0.21mm0.30mm

Maximum bending moments

Mfmáx,n=Q·Lx·xL·fv,ni=1Nfv,i=268KN·m135KN·m24KN·m34KN·m

Maximum shear forces

Qmáx,n=Q·LxL·fv,ni=1Nfv,i=31KN16KN3KN4KN

Truck axle 2

Deflection in the center of the spam of the longitudinal beams

fcl,Qx,i=fv,i·sinπ·xL=4.27mm2.15mm0.38mm0.53mm

Maximum bending moments

Mfmáx,n=Q·Lx·xL·fv,ni=1Nfv,i=468KN·m236KN·m41KN·m59KN·m

Maximum shear forces

Qmáx,n=Q·LxL·fv,ni=1Nfv,i=39KN20KN3KN5KN

Truck axle 3

Deflection in the center of the spam of the longitudinal beams

fcl,Qx,i=fv,i·sinπ·xL=3.69mm1.86mm0.33mm0.46mm

Maximum bending moments

Mfmáx,n=Q·Lx·xL·fv,ni=1Nfv,i=409KN·m206KN·m36KN·m51KN·m

Maximum shear forces

Qmáx,n=Q·LxL·fv,ni=1Nfv,i=41KN21KN4KN5KN

Truck axle 4

Deflection in the center of the spam of the longitudinal beams

fcl,Qx,i=fv,i·sinπ·xL=3.43mm1.73mm0.30mm0.43mm

Maximum bending moments

Mfmáx,n=Q·Lx·xL·fv,ni=1Nfv,i=385KN·m194KN·m34KN·m48KN·m

Maximum shear forces

Qmáx,n=Q·LxL·fv,ni=1Nfv,i=45KN22KN4KN6KN

Deflection and efforts generated by the action of the four-axle truck

fv,i=13.78mm6.94mm1.22mm1.72mm
Mfmáx,i=1,323KN·m667KN·m117KN·m165KN·m
Qmáx,i=138KN70KN12KN17KN

## 4. Conclusions

The proposed method for the analysis of the cross-sectional distribution of live loads on girder bridge decks allows determining the cross-sectional distribution in different configurations of girder bridge decks without the need to resort to complex calculation models that involve depth computing power and excessive analysis time. The method is also applicable for modern synthetic materials, such plastic composites, self-repair. The simplicity of the method allows an easy integration into optimal bridge design strategies  or more heuristic approaches [29, 30, 31, 32, 33] to challenge today’s competitive world intelligently.