Skip to main content

Beam Grid Example

Beam Grid Example

1. Beam grid

The beam grid is a generally horizontal structure consisting of longitudinal beams, called "main beams", associated with transverse beams. Beam grid modeling allows the calculation of a large number of structures in civil engineering, it is mainly used for the following works:

  • the bridge deck with prefabricated beams made of prestressed concrete with bonded wires, or reinforced concrete girders;

  • the mixed-frame deck when it has more than two beams;

  • the ribbed slab deck;

  • by extrapolation, the slab type deck.

This modeling has the advantage of imposing the main directions of the reduction elements. Indeed, the plate calculation produces principal directions of moments that often deviate from the directions of the reinforcement, and the complex system of deviatoric moments or mxy torsional moments must be taken into account in the strength calculation.

Whereas, for a wireframe model, when the directions of the reinforcement are fixed, according to the geometry of the structure and the recommendations, we will be able to deduce the directions of the beams system in our model. The beam grid remains a simple model to model and above all simple to use. Consequently, we are often led to model a concrete slab (bridge one-way joist slab) by cross girders (transverse beams).

1.1 The GUYON-MASSONNET method

In the 1950s, the construction of multi-beam bridges required calculations tools. The method developed by Yves Guyon [Calculation of Wide Bridges with Multiple Girders Joints with Cross Girders], then improved by Charles Massonnet [Contribution and Complements to the Calculation of Multi Girder Bridges] is the most widely used, it has been compared to model tests or measurements on real bridges and has always been found to be in good or very good agreement with experimental measurements. The method uses two fundamental parameters whose data can be used in a FEM.

The data:

1.- the flexural stiffnesses of the beams and cross girders:

BP=E.IP

BE=E.IE

2.- the torsional stiffnesses of the beams and cross girders:

CP=E.KP

CE=E.KE

The method assumes that all the beams are identical. However, in order to make building easier, the edge beams can be different!

The deck jacking operations require the construction of an end cross girder:

Moreover, this is an isostatic bridge, and even if the method was extended to hyperstatic continuous bridges by the rule of thumb which consists in calculating a fictitious span with the same elastic deflection under a concentrated load applied in the middle of the span, what about the stresses on supports? Can the moments on the piles be calculated using a method based on the deformation calculation?

Modern calculation devices have favored geometrical and architectural criteria in construction and have taken many works out of the method's application field. The FEMs seem to abolish the perimeter of calculation of the structures, and the difficulty often lies in the choice of the data and the model and the exploitation of the results.

Note. - In reinforced concrete it is preferable to use reducing elements rather than stress tensors.

1.2 Purpose of this chapter

The simplicity of exploiting the results makes the beam grid a widely used calculation model. In this chapter, we will discuss data elaboration and the transformation of the results for the effects of the skew and the implementation of hyperstatic spans. These computational difficulties are discussed through the following examples:

  1. Straight isostatic bridge deck with a 100 degree skew.

  2. Bridge deck with 2 hyperstatic straight spans with a 100 degree skew.

  3. Straight isostatic bridge deck with a 70 degree skew.

  4. Straight isostatic bridge deck with a 50 degree skew.

  5. Straight 2 span hyperstatic bridge deck with a 50 degree skew.

1.3 Beam characteristics

As for many statically indeterminate structures, the results are subject to the relative stiffness of each beam grid. The load distribution on the main beams is favored when the cross girders are stiff in relation to the main beams.

1.3.1 Flexural beam stiffness 

The main beams are usually surmounted by a continuous one-way joist slab and the beam grid is made up of a single plan, therefore the sections of each structural element should be well defined.

The geometry of the compression flanges is defined in EN 1992-1-1 art. 5.3.2.1.

with:

1.3.2 Torsional beam stiffness 

Due to the lack of discontinuity between the beams, we cannot strictly speak of torsion, it is rather a plate bending. However, in order to remain consistent with the beam logic, a torsional inertia is determined for each element. By considering the beam free to rotate around its center of torsion, the stiffness is underestimated. Therefore, we will use for rectangular elements a shape coefficient K = 1/3, whatever the slenderness of the part is: J=ab33(*).

Note. - K is a shape coefficient studied by Mr. Caquot. He indicates that K can be evaluated with a very good approximation by the formula 1k=1+1m23.560-0.56m-1m+12 where m is the ratio ab with a≥b… according to “Reinforced Concrete Form” by  R. Chambaud et P. Lebelle.

(*) Important: The SETRA (operating society for transport and automobile repairs)/CEREMA (center for studies and expertise on risks, environment, mobility and development) guides, the PRP 75, the guide relating to independent-span and post-tensioned prestressed beam viaducts and the guide relating to viaducts with prefabricated beams made of prestressed concrete with bonded wires all offer a more sophisticated approach - it will be necessary to refer to them in the context of a real project.

1.4 Model geometry 

1.4.1 Model characteristics

The model characteristics are:

  • Span: 15.0m;

  • Width: 9.4m;

  • 9 main T-beams with a rib of 0.80x0.40m;

  • A 25cm thick one-way joist slab.

Cross section of the deck

Longitudinal section of the isostatic deck

1.4.2 Geometrical characteristics of the sections

The characteristics of the sections are:

  • Common main beams:

  • Section:

  • Flexural inertia:

  • Torsional inertia:

  • Common cross girders:

The model is made with a cross girder every 150cm:

  • Section:

  • Flexural inertia:

  • Torsional inertia:

The FE model presented below is made with ROBOT v2019 -Autodesk:

1.4.3 Note on torsional inertia

By default, the software calculates the torsional inertia. The table below compares the different calculations:

J[m4]

Software

J=∑ab3/3

J=∑ J_Caquot

Beam

0.0191936

0.0261760

0.0199764

Cross girder

0.0069919

0.0078125

0.0069646


In order to be able to compare the finite element method to the Guyon-Massonnet method we will use the torsional inertias from the software.

1.5 Loading of the model

Each beam of the different geometries is loaded with:

  • A mid-span point load of 100kN:

  • A 10kN/m distributed load along each beam:

1.6 Model results

1.6.1 Comparison of simple geometry with the Guyon-Massonnet method

The mid-span moments per beam for a 100-degree isostatic deck are (for e=y):

Beams y (m)

Guyon-Massonnet method


FEM model


K

Point load

Miso=P∙l4

Distributed load

Miso=q∙l28

Point load


Distributed load


Edge

0.405

150 kN.m

112 kN.m

169 kN.m

99 kN.m

-3.375

0.244

91 kN.m

68 kN.m

121 kN.m

65 kN.m

-2.250

0.188

70 kN.m

52 kN.m

107 kN.m

52 kN.m

-1.125

0.170

63 kN.m

47 kN.m

102 kN.m

48 kN.m

Center

0.166

62 kN.m

46 kN.m

101 kN.m

47 kN.m


It should be noted that the application of the Guyon-Massonnet method results in greater stresses on the edges than the FE model.

This can be explained by the fact that the torsional inertia of the beams joined by the one-way joist slab is underestimated. Indeed, if the torsional inertia of the main beams is higher, the strain of the edge beam will be taken off, the load being better distributed over the entire deck.

Below is a study in which the torsional inertia of the main beams Kp is affected by a coefficient (noted f_Kp) varying from 0.10 to 3.00. The graph presents the comparison of the moments obtained by applying the Guyon-Massonnet method to the two end beams with e=y:

The different torsional inertias used are listed below:

[m4]


J = Σab3/3


J=Σ
J_Caquot


Software

FEM = G-M

[y=b; e=b]


FEM = G-M

[y=3b/4; e=3b/4]


Kp

0.0261760

0.0199764

0.0191936

0.0282102

0.0258026


These results indicate the importance of the evaluation of torsional stiffness in the calculation of moments by the Guyon-Massonnet method. The presence of the one-way joist slab tends to increase the main beams’ torsional inertia.

1.6.2 Consideration of the geometry of the edge beams and end cross girders

The mid-span moments per beam for a 100-degree isostatic deck are (for e=y):

Beams 

Moment - Bar model:

beams of identical inertia = basic case

Moment - Bar model: variable inertia beams (edge) and end cross girders


y(m)

Point load

Distributed load

Point load


Impact

Distributed load


Impact

Edge

169 kN.m

99 kN.m

154 kN.m

-9.81%

86 kN.m 

-15.24%

-3.375

121 kN.m

65 kN.m

123 kN.m

1.31%

66 kN.m

1.62%

-2.250

107 kN.m

52 kN.m

107 kN.m

0.60%

53 kN.m

1.15%

-1.125

102 kN.m

48 kN.m

102 kN.m

0.01%

48 kN.m

0.02%

Center

101 kN.m

47 kN.m

100 kN.m

-0.23%

46 kN.m

-0.47%

These results, which are far from calling into question the method’s validity, indicate the importance of the decrease in the stiffness of the edge beams and the influence of the end cross girder.

Note: the torsional inertias retained are given below (software column):

J[m4]


Software

J = Σab3/3


J=ΣJ_Caquot


Beam

0.0191936

0.0261760

0.0199764

Cross girder

0.0069919

0.0078125

0.0069646

End cross girder

0.0232905

0.0318938

0.0232925

Edge beam

J=ΣJ_Caquot

0.0242880

0.0181062

1.6.3 Complex model results

Max. span moment


Beams – y(m)

edge

-3.375

-2.250

-1.125

Center

Isostatic deck with a 70 degree skew

Point load

147 kN.m

120 kN.m

107 kN.m

101 kN.m

99 kN.m

Distributed load

77 kN.m

61 kN.m

51 kN.m


44 kN.m


43 kN.m


Note

The influence of the skew is very small and the common cross girders (hollow core elements) are arranged skewed.

Isostatic deck with a 50 degree skew

Point load

138 kN.m

105 kN.m

98 kN.m

88 kN.m

93 kN.m

Distributed load

72 kN.m


54 kN.m


43 kN.m


38 kN.m


37 kN.m


Note

The common cross girders are laid out straight. The end cross girder, on the obtuse angle side, is in flexion rather than torsion, therefore it reduces the span moments to its detriment.

Statically indeterminate deck (2 spans) with a 100 degree skew

Point load

130 kN.m

105 kN.m

93 kN.m 

90 kN.m

88 kN.m

Distributed load

65 kN.m


51 kN.m


41 kN.m


37 kN.m


38 kN.m


Note

There is no variation on the load distribution between beams caused by the deck being statically indeterminate

Statically indeterminate deck (2 spans) with a 50 degree skew

Point load

132 kN.m

97 kN.m

91 kN.m

82 kN.m

87 kN.m

Distributed load

66 kN.m


48 kN.m


36 kN.m


33 kN.m


31 kN.m


Note

The fact that the deck is statically indeterminate causes little variation in the span moments.

Max. pile moment


Beams – y(m)

edge

-3.375

-2.250

-1.125

Center

Statically indeterminate deck (2 spans) with a 100 degree skew

Point load

-63 kN.m

-40 kN.m

-33 kN.m

-31 kN.m

-30 kN.m

Distributed load

-70 kN.m


-46 kN.m


-41 kN.m


39 kN.m


38 kN.m


Note

No variation on the load distribution between beams caused by the deck being statically indeterminate

Statically indeterminate deck (2 spans) with a 50 degree skew

Point load

-84 kN.m

-53 kN.m

-31 kN.m

-28 kN.m

-27 kN.m

Distributed load

-92 kN.m


-46 kN.m

-37 kN.m

-36 kN.m

-36 kN.m

Note

There is a concentration of moments close to the obtuse angle of the load.

1.7 Conclusion

For a classical bridge deck, i.e. a high stiffness of the beams compared to the cross girders, the difference of the finite element model beam grid compared to other methods is mainly expressed through the consideration of a lower stiffness of the edge beams and the consideration of the end cross girders.

When the skew of the deck is pronounced, only a finite element model can accurately value the reduction elements of the beams.