Next Article in Journal
Damage Identification of Turnout Rail through a Covariance-Based Condition Index and Quantitative Pattern Analysis
Previous Article in Journal
The Application of Bamboo in the Railway Industry: A Sustainable Solution for Track Construction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application and Experimental Validation of Seven-Degree-of-Freedom Beam Element for Girder Bridges during Deck Construction

1
Civil Engineering Department, University of Louisiana at Lafayette, Lafayette, LA 70503, USA
2
Ameren Corp., St. Louis, MO 63103, USA
3
Department of Civil, Computer and Electrical Engineering, St. Louis University, St. Louis, MO 63103, USA
*
Author to whom correspondence should be addressed.
Infrastructures 2023, 8(12), 175; https://doi.org/10.3390/infrastructures8120175
Submission received: 3 November 2023 / Revised: 4 December 2023 / Accepted: 6 December 2023 / Published: 8 December 2023

Abstract

:
During bridge deck construction, the deck finishing machine and the fresh concrete often produce large vertical loads and torsional moments acting on the bridge girder system. In some cases, these loads can cause excessive vertical deflection and transverse rotation in the bridge girders, leading to many maintenance and safety problems, such as changes in deck thickness and local and global instabilities during construction. To minimize the potential problems caused by deck construction, the AASHOTO LRFD Bridge Design Specification requires consideration of these torsional moments during the design procedure, and a detailed three-dimensional finite element analysis may be conducted. However, for bridge girders with open-section thin-walled sections, only the solid or shell element can be used to recognize the warping of the girder since the torsional warping effect is not included in the classical beam element. In this research, a warping degree of freedom was added to a beam element, and a three-dimensional beam element with seven degrees of freedom (7-DOF) at each node was derived as an alternative method for analyzing girder bridges during deck construction. A computer program based on the 7-DOF beam element was also developed in MATLAB. To assess the 7-DOF beam element, one bridge was selected to measure the transverse rotation, vertical deflection, and stress of the exterior girder and the first interior girder during deck construction. Also, three full-scale numerical models using solid elements, classical three-dimensional beam elements, and 7-DOF beam elements were created based on the geometries and loads of the experimental bridge. A comparative study was conducted by comparing the results from the numerical models and experimental monitoring data to evaluate the 7-DOF beam element. The results showed that the 7-DOF beam element had excellent behavior in analyzing the girder bridges under construction load, especially in the torsional analysis of bridge girders. Also, unlike the solid element model, which also provided reasonable results, the 7-DOF beam element model can compute the internal forces of the cross-sections along the bridge, which allows the 7-DOF beam element to be an alternative approach for design and research requiring less modeling effort and computational complexity.

1. Introduction

Girder bridges are the most common type of small- and medium-span bridges, consisting of several similarly sized longitudinal girders spaced uniformly across the bridge’s width. The deck slab is often designed to span transversely between the girders and extend past the exterior girder. This design increases the deck’s width without the need for additional girders. Deck slab construction typically employs an overhang formwork system and a deck finishing machine. The overhang formwork is usually supported by steel brackets placed at intervals of 90 cm (3 ft) to 120 cm (4 ft) along the exterior girder. These brackets support the weight of fresh concrete, the deck finishing machine, and other construction loads [1,2,3]. Figure 1 shows a typical formwork system and construction loads acting on bridge girders during deck construction. These overhang brackets are often attached to the top flange of the exterior girder and provide support against the bottom of the girder’s web. The deck finishing machine, which creates significant loads during construction, is typically located on the edge of the overhang.
During bridge deck construction, a primary concern is the torsional moment exerted on the exterior girder [4,5,6,7]. The bridge shows less stiffness as the upper concrete deck slab does not contribute as effectively to the structure’s rigidity. As the construction loads are transferred through the overhang bracket formwork system to the exterior girder, significant torsional moments are usually produced and sometimes can result in excessive transverse rotation of the girder [8]. Such excessive rotation can lead to an uneven deck thickness, causing maintenance challenges during the bridge’s service life [8]. Additionally, local and global instabilities are potential issues, particularly for steel plate girders with slender webs when subjected to eccentric loads during deck construction [9,10,11,12]. These challenges can sometimes lead to bridge failures in steel and concrete bridge systems [8].
To mitigate potential issues arising from deck construction, the AASHTO LRFD Bridge Design Specification mandates the consideration of both lateral bending stress and the effects of bracket reactions during the design phase [13]. However, analyzing the lateral–torsional behaviors of bridge girders, especially those with open-section thin-walled beams like I-beams, remains challenging due to warping effects. The Guidelines for Steel Girder Bridge Analysis (2019) present a method that employs an equivalent torsion constant [14]. This method combines the St. Venant torsional stiffness for the open-section thin-walled beam, considering warping fixity at both ends of a specified unbraced length. A limitation of this approach is its assumption of the girder being simply supported, necessitating the calculation of the equivalent torsional constant for every variation in unbraced length or changes in the girder’s cross-sectional properties.
Another common approach to account for warping effects utilizes the finite element method. The AASHTO LRFD Bridge Design Specification and The Guidelines for Steel Girder Bridge Analysis advocate for the finite element method. Specifically, a three-dimensional model with shell or solid elements is suggested to address warping torsion. This is because beam elements cannot recognize the torsional warping degree of freedom [14]. However, implementing finite element modeling with shell or solid elements requires extensive expertise in structural analysis and significant modeling efforts.
To incorporate the warping effect in the analysis of structures with thin-walled open sections, V.Z. Vlasov formulated a torsional theory that includes the impact of restraint warping [15]. In this torsional theory, the warping constant and bi-moment were introduced to account for torsional restraint warping, sectional deformation, and the influence of shear strain on normal stress distribution. Moreover, several studies have been conducted to provide a deeper understanding of the mechanical behavior associated with torsional warping [16,17,18,19,20,21,22,23]. Specialized finite beam elements have been developed for specific applications based on the torsional warping theory. S.H. Zhang et al. introduced a beam element that accounted for the impacts of warping, distortion, and shear lag in box beams [24]. This element features two nodes, each with nine degrees of freedom. Later, in 2017, F. Cambronero-Barrientos proposed a three-dimensional, three-node beam element [25]. It has five degrees of freedom per node and is designed to incorporate effects like shear lag, torsion, distortion, and both homogeneous and non-homogeneous distribution of normal stress.
In addition, various general-purpose beam elements accounting for warping effects have been developed and incorporated into several finite element software packages (e.g., ABAQUS/CAE [26], ANSYS [27], and ADINA [28]). Typically, these software solutions introduce a seventh degree of freedom (7-DOF) at nodal points to encompass the warping effect. However, in certain scenarios, the analysis executed by these software tools might provide inconsistent results, particularly for beams undergoing non-uniform torsion [29]. Expanding on this, Sabat and Kundu explored advanced modeling techniques in ANSYS for complex torsional analysis, enhancing the accuracy of predictions [30]. Similarly, Li et al. proposed improvements in warping analysis in ABAQUS, focusing on the shear lag effect for a thin-walled box-section beam [31].
Some other researchers also developed beam elements for general purposes. For instance, E.J.Sapountzakis formulated a beam element with a 14 × 14 stiffness matrix based on the boundary element method (BEM) [32]. This beam element accounted for warping and shear deformation effects, catering to nodal load vectors of any given homogeneous or composite cross-section. Moreover, Yau and Kuo introduced a new 14 × 14 stiffness matrix for a beam element to analyze the warping effect on the buckling of the I-beams [33]. Their approach leaned on the Vlasov torsional theory and operated under the assumption that the I-beam was perfectly straight, constructed from three slender flat plates, and that the section remained bi-symmetric. Bernuzzi and Simoncelli investigated the interaction between axial forces, bending moments, and bi-moments in the linear elastic range, especially focusing on members with mono-symmetric cross-sections. Their research evaluated the impact of various transformation matrices on analysis outcomes, highlighting how the selection of these matrices significantly affects the results [34]. Furthermore, extensive research has been undertaken to innovate numerical analysis methods for beam warping analysis, as evidenced by several studies in this field [35,36,37,38,39].
In this study, a three-dimensional beam element with seven degrees of freedom (7-DOF) was derived, grounded in the kinematics of beams under non-uniform torsional loads and the theory of potential energy. Compared with the classical three-dimensional beam element with six degrees of freedom in each node, the twist angle per unit length was taken as an additional degree of freedom to represent the torsional warping effect. Based on the proposed 7-DOF beam element, a finite element computer program was developed in this study. In order to evaluate the 7-DOF beam element in analyzing the bridge girders during deck construction, one bridge located in the state of Illinois (USA) was monitored and finite element analyses were conducted using different elements. Four finite element models were established using a solid element, a classical beam element, a B32OS element from ABAQUS/CAE, and the 7-DOF element. The finite element results were compared with the experimental data, and they showed that the 7-DOF frame element developed in this study yielded accurate and reliable results. Additionally, the 7-DOF beam element achieved this with significantly reduced modeling effort and computational time compared to other methods.

2. The 7-DOF Beam Element

The major difference between the 7-DOF beam theory and classical beam theory is the assumption of the warping effect for the beam section. The classical beam theory basically assumes that a cross-section orthogonal to the x-axis always remains plane and keeps its shape during deformation. It adopts the St-Venant torsional theory, which assumes that the cross-section is free of warping or the warping along the beam is constant [40,41]. However, for an open thin-walled cross-section such as an I-beam that is widely used in structural design and construction, the warping effect is often dominant during the torsional analysis, especially when the beam is under significant torsional load [25].
The 7-DOF beam element considers warping deformation as an additional degree of freedom during the analysis. Generally, the rotation caused by the torsional moment reduces because the torsional moment is resisted by both St-Venant torsional stress and warping stress. A typical 7-DOF beam element is shown in Figure 2.
The force applied at each element and the corresponding deformation are written as Equations (1) and (2):
F = [ F x 1 F y 1 F z 1 M x 1 M y 1 M z 1 M w 1 F x 2 F y 2 F z 2 M x 2 M y 2 M z 2 M w 2 ] T
U = [ u x 1 u y 1 u z 1 θ x 1 θ y 1 θ z 1 θ w 1 F x 2 F y 2 F z 2 θ x 2 θ y 2 θ z 2 θ w 2 ] T
where M w is the warping moment and θ w is the angle of twist per unit length.

2.1. Deformation and Kinematics of 7-DOF Beam Element

When the beam is restrained at one or more locations along the span, the section no longer remains flat, and additional deformation in the axial direction is induced. Figure 3 shows the general beam deformation in each coordinate plate.
By adding warping-induced axial deformation, the deformation on each coordinate axis can be written as Equations (3)–(5):
u x x , y , z = w x x + z θ x x y θ z x + ω ( y , z ) θ x , x
u y x , y , z = w y x z θ x x
u z x , y , z = w z x + z θ x x
where u x , u y , and u z are the accumulated deformation in the x , y , and z directions, respectively, and ω ( y , z ) is the warping function, which basically is a shape function defining the axial deformation of the cross-section from the rotation component.
As a result of including warping, the strain at a beam cross-section can be described with the following equations (Equations (6)–(8)):
ε x x = u x , x = w x , x + z θ y , x y θ z , x + ω ( x , y ) θ x , x x
γ x y = u x , y + u y , x = w y , x θ z + ( ω , y z ) θ x , x
γ x z = u x , z + u z , x = w z , x + θ y + ( ω , z + y ) θ x , x
Equations (6)–(8) indicate that the warping of the beam section not only produces additional axial strain but also induces extra shear strains on the cross-section. The stress at a certain point of the cross-section can be simply computed based on Hook’s Law.

2.2. Potential Energy and Stiffness Matrix

To derive the stiffness matrix for the 7-DOF beam element, the potential energy principle was applied. When external forces are applied on the beam element and perform work, the energy is stored in the form of stress and elastic deformation [42,43]. This means the work done by the external forces must be equal to the potential energy stored in the beam. The total potential energy of a linear elastic body can be defined using Equation (9):
Π = U W
where U is the strain energy and W is the total work done by the applied forces, which can be written as Equations (10) and (11):
U = 1 2 E V ε x x 2 d V + 1 2 G V γ x y 2 d V + 1 2 G V γ x z 2 d V
W = F T × U
For the conservative structural system, the kinematically admissible deformations that correspond to the equilibrium state extremize the total potential energy. If the extremum is a minimum, the equilibrium state is stable [44]. Therefore, the minimum potential energy (MPE) can be used to find relations between the applied forces and deformations, which generates a stiffness matrix. The MPE for the 7-DOF beam system can be described using Equation (12).
Π U = 0
By solving the MPE, the relation between the applied force and displacement can be written as Equation (13).
F = K e u
The element stiffness matrix K e is a 14 by 14 symmetric square matrix.

3. Field Bridge Experiment

One bridge located in Lawrence County (Lawrence County, IL, USA) was selected for this research to validate the 7-DOF frame element model. The instrumented bridge is a three-span continuous bridge with six plate girders. The overall span length of the bridge is 157 m (516 ft). The steel plate girders are 167.64 cm (66 in) in depth with 50.8 cm (20 in) wide flanges and are spaced at 2.18 m (86 in). The overhang width is 1.12 m (44 in). A framing plan of the bridge, girder elevations, and details of the cross-frames is shown in Figure 4, Figure 5 and Figure 6.
The deck pouring of the experimental bridge was divided into two sequences. As shown in Figure 7, the first sequence covered from both abutments to 34.4 m (112.75 ft) along the side spans and the middle portion of the second span. In this stage, the screed machine moved from the south abutment to the north abutment along the bridge span. The second deck pour placed the fresh concrete on the rest of the bridge three days after the first pour, with the screed machine moving in the same direction. The fresh concrete was poured at night for both sequences to control early cracking due to temperature effects.

3.1. Bridge Instrumentation

The exterior girder (G1) and first interior girder (G2) of the first span of this bridge were instrumented and monitored for deformation and strain during the first concrete deck pour in the southern side span (shown in Figure 8). Three sections (S1, S2, and S3) were identified for the installation of the sensors. At each section, tilt sensors were installed on the middle and bottom of girder web on both exterior and interior girder. Prisms were installed at the bottom flange of each section in order to measure the deflection of the girders using total stations. Two stain gauges were also placed on the bottom flanges in S2, as shown in Figure 9.
The transverse rotation of the bridge girders was measured using dual-axial (CXTLA02) tilt sensors. The tilt sensors have a twenty-degree maximum range in both directions. The sensitivity is 1% and the resolution is 0.03 deg rms. Also, an open aluminum box was used to hold the tilt sensor for protection, as shown in Figure 10a. Foil strain gauges (CEA-06-125UN-350) were installed on the bottom of the girder flange to measure the strain changes induced by deck construction loads. The strain gauges have a sensitivity of 350 ± 0.3% ohms, a resolution of 1 μm/m, and a range greater than 5000 μm/m. A protective coating was also used to avoid the effect of moisture and protect the strain gauges from damage (Figure 10b). The deflection of the bridge girders during deck construction was measured using the total station, as shown in Figure 11.

3.2. Experimental Results

The readings of sensors were collected at thirty-second intervals during the first pouring sequence of the experimental bridge. The tilt sensors measured both the inward (positive) and outward (negative) transverse rotation of the beam. However, due to the small measuring interval and the vibration of the screed machine on the top of the bridge, the experimental data from tilt sensors showed fluctuations during the whole monitoring period, making it difficult to recognize the changes in transverse rotation that occurred during deck construction. Therefore, the exponential moving average (EMA) method, a weighted moving average method used for smoothing time-series data, was used to reduce these fluctuations in the experimental data.
The EMA method was also used to calculate the value of rotation. Generally, the maximum rotation for a specific section occurred when the screed machine and fresh concrete moved right at the section. Figure 12 shows the exterior girder rotation at the bottom of the web in section S1. The results indicate that a maximum rotation value of 0.08 degrees occurred when the screed machine moved above section S1.
The collected strain data from the bottom of the G1 and G2 girders was analyzed. Figure 13 illustrates the stress time-series data, indicating the behavior of the bridge girders, specifically the bottom of both the G1 and G2 girders, throughout the deck construction phase. It is observed that G1-Bottom experienced a series of fluctuations before establishing a steady increase, closely mirroring the trend exhibited by G2-Bottom. This behavior can be attributed to the fact that G1, being an exterior girder, endures more torsional moments during the deck construction. As the screed machine operated over the S1 section, the stress values peaked at 61.7 MPa (8.95 ksi) and 47.2 MPa (6.85 ksi) for the G2 girders.
Figure 14 illuminates the time-dependent vertical deflection profiles for three distinct sections: S1, S2, and S3. Similar to the stress and rotation, peak deflection is noted when the screed machine is directly positioned over the S1 section. The recorded deflections for the sections are 5.93 cm (2.33 inches), 7.02 cm (2.76 inches), and 7.68 cm (3.02 inches), respectively.

4. Numerical Modeling

Three finite element models using solid elements, classic beam elements, and 7-DOF beam elements were created and compared with the in-field measurement to evaluate the accuracy of these models. To ensure that all three models were comparable, the loads applied in the models were based on where the maximum transverse rotation occurred in section S1. In this case, the fresh concrete covered from the south abutment to the location of section S1, and the screed machine was right at section S1. Figure 15 shows the location of the fresh concrete and screed machine used in the finite element analysis. The mesh size of the girders was determined to be six inches for all three models.
The solid element model for the experimental bridge was developed in ABAQUS/CAE. The C3D8R element, which is widely used in stress and displacement analysis, was used in the modeling of bridge girders. The cross-frames were modeled using the C3D10 element due to the irregular shape and connection of the components. The C3D10 element is a second-order tetrahedral three-dimensional stress element with ten nodes in each element, which can generate meshes for complex geometry. Figure 16 shows the full-scale finite element model for the experimental bridge in ABAQUS/CAE. The bridge girders and cross-frame were modeled based on the design dimensions. A single-angle K-frame without a top truss was used in this bridge at all non-support locations (shown in Figure 17a). For the support location, an X-frame was used with both top and bottom trusses (shown in Figure 17b). Welded connections between the girder webs and connection plates of the cross-frames were simulated using the tie connection.
The classical beam element model was created using the SAP2000 v14 software. All the structural components including girders and cross-frames were modeled using classical three-dimensional beam elements. In order to obtain a better comparison, the cross-frames were modeled based on the same geometrical dimensions as used in the solid element model. Rigid beam elements were also used to connect beams to cross-frames to maintain the same cross-frame depth as specified in the design. Figure 18 shows the classical beam element model that was built in SAP2000.
A finite element computer program was developed in MATLAB in order to evaluate the 7-DOF beam element developed in this study. The general procedures for finite element analysis using 7-DOF beam elements are shown in Figure 19. The program required the input of the basic bridge information such as girder dimensions, span lengths, overhang width, girder spacings, details of lateral bracing and mesh size, etc. Then, all the information for generating local and global stiffness matrices was created automatically. The local stiffness matrix and coordinate transformation matrix were 14 × 14 matrices. The size of the global stiffness matrix was larger than that in the classical beam element model but much smaller than the solid element model, depending on the total number of nodes. The program outputs included deformations, internal forces, reaction forces, and stresses.
The 7-DOF beam element model was similar to the classical beam element model. These two models were able to share the same node coordinates and element information. However, much attention must be paid when generating the force vector and boundary conditions since the warping degree of freedom was added to the stiffness matrix. The 7-DOF beam element created by MATLAB is shown in Figure 20. Both the girders and cross-frames were modeled using the 7-DOF beam element developed in this study with a mesh size of six inches.
The loads applied on the 7-DOF beam element model were the same as those on the classical beam element model since they were beam elements in both cases. In terms of programming, the analysis using the 7-DOF beam element did not require a total rewrite of the finite element program since only minor modifications were required for the stiffness matrix, boundary conditions, and force vectors.

5. Results

5.1. Rotation and Deflection

The torsional behaviors of the bridge girder are considered the major difference between the classical and 7-DOF beam element models. Figure 21 shows a comparison of transverse rotation on the exterior and first interior girder among three numerical models and experimental data in section S1. The experimental data and solid element model indicated that the middle and bottom of the girder web had very similar transverse rotations for both the exterior girder (G1) and the first interior girder (G2). The 7-DOF element model and solid element model showed very close transverse rotation for the exterior girder—0.087 and 0.073 degrees, respectively, which agreed with the experimental data of 0.087 degrees for the middle of the web and 0.080 degrees for the bottom of web. However, the classical beam element model provided a rotation value of 0.46 degrees for the exterior girder, which was around six times larger than the experimental data and 7-DOF beam element model. Only minor differences were observed among the finite element models and experimental data for the interior girder. This is because no torsional moment was applied to the interior girder, and the rotation resulted from relative deflections among girders as well as the bending of the bracing systems.
A comparison between rotations from the numerical models and experimental data for section S2 is shown in Figure 22. The classical beam model shows a significantly large rotation value of 3.16 degrees for the exterior girder, which is twenty-six times larger than the experimental data (0.12 degrees). This is due to the absence of bracing in section S2, and the classical beam element cannot account for the warping constraint at the adjacent bracings. Furthermore, the 7-DOF beam element and solid element model observed rotations of 0.10 and 0.11 degrees, respectively, which are reasonably close to the experimental data, with a difference of less than 15%. For the interior girder, both the numerical analysis results and experimental data showed rotations of less than 0.10 degrees, as no torsional moment was directly applied to the interior girder.
Figure 23 shows the rotations of both the exterior and the first interior girder in section S3. The 7-DOF beam model yielded a rotation value of 0.09 degrees, which was comparable to the results from the solid element model and experimental data, at 0.11 degrees and 0.12 degrees, respectively. However, the rotation value from the classical beam element model was notably off, showing 0.50 degrees, which is 500% greater than both the 7-DOF beam element result and the experimental data. For the interior girder in section S3, the 7-DOF beam element model indicated a rotation of 0.06 degrees, showing less than a 15% difference when compared to the solid element model and experimental data, both at 0.07 degrees. However, a 40% difference was observed when comparing the results from the classical beam element model to the others.
The comparison among numerical models, which use the 7-DOF beam and other elements, as well as experimental data, demonstrates that the 7-DOF beam element is capable of analyzing the torsional behaviors of girder bridges with thin-walled open sections during deck construction. The transverse rotations observed for both the exterior and the first interior girders across three different sections were reasonable when assessed using the 7-DOF beam element. Moreover, using the 7-DOF beam element requires less computational time and effort in modeling.
Figure 24 presents a comparison of vertical deflection among the numerical models and experimental data. Since warping stiffness does not significantly influence the bending behaviors of the bridge girder, the 7-DOF beam element model exhibits nearly identical deflections to the classical beam element model in the vertical direction across all three sections. Additionally, only slight differences are observed between the 7-DOF beam element and the solid element model, with the maximum variation being 2.3%. The experimental data are approximately 30% lower than the results derived from the 7-DOF beam element. This discrepancy can be attributed to the numerical model’s oversight of the added stiffness provided by the formwork system. In general, the 7-DOF beam element aligns with the other numerical methods when analyzing vertical deflection.

5.2. Stress

As shown in Figure 25, the stress analyzed using the 7-DOF beam element was compared with the other numerical models and experimental data at section S1, where the strain gauges were installed on the bottom flange of the exterior girder (G1) and first interior girder (G2). Stresses of 63.9 MPa (9.27 ksi) in the exterior girder and 61.7 MPa (8.95 ksi) in the interior girder were computed using the 7-DOF beam element. The solid element and classical beam element models gave similar stress values with less than a 1% difference compared with the 7-DOF beam element model. The stresses collected from the experimental bridge were slightly smaller than the numerical analysis result. The differences between the 7-DOF beam element model and experimental data were 4% and 20% for the exterior girder and interior girder, respectively. It is also important to note that the warping effect can change the stress distribution on the cross-section according to Equations (6)–(8). For the experimental bridge, since the primary stress on the cross-section was due to the bending of the girder, the warping stress had a very limited effect on the stress calculation. However, if the bridge girder is under a large torsional moment, the stress induced by the warping effect can significantly affect the stress distribution of the cross-section.

6. Discussion

This comparative study confirmed the accuracy of the 7-DOF beam element in analyzing girder bridges during deck construction. This element stands out for its inclusion of warping as an integral degree of freedom, enhancing the realism and accuracy of structural simulations. This approach marks a significant deviation from traditional methodologies employed in other studies, which often focus on warping beam elements.
The purpose of the 7-DOF beam element in this study is to provide a method for the detailed analysis of complex structural behaviors, especially in scenarios involving girder bridges. It includes a comprehensive set of degrees of freedom, encompassing warping to simulate real-world structural responses more accurately. This is in contrast to the warping beam elements used in other studies, which are primarily designed to address warping deformations in structural members and are less sophisticated.
In terms of accuracy, the 7-DOF beam element demonstrates a high degree of precision, particularly in predicting torsional behavior and warping effects. Its advanced formulation allows for a more detailed representation of real-world conditions, which is a notable improvement over the simpler warping beam elements of other studies. These traditional warping elements, while effective for modeling warping deformations, may not capture the full spectrum of structural behaviors, especially under complex loading conditions.
Validation was achieved by comparing the results from the 7-DOF beam element model with those from the solid element models and experimental data. The classical beam element was found to be inadequate in accurately predicting the torsional behavior of the bridge under torsional moments. This inadequacy largely stemmed from the classical beam element’s assumption that the girder experienced no restrained warping. In contrast, the 7-DOF beam element introduces warping as an additional degree of freedom. This inclusion furnishes the girder with added stiffness to counteract the applied torsional moment. The results from the 7-DOF beam element, both in terms of deformation and stress, aligned closely with those from the solid element model and experimental data. Thus, a key advantage of employing the 7-DOF beam element is its capability to precisely analyze torsional behaviors without the need for intricate solid element models. This efficiency translates to reduced time spent on numerical modeling.
Another significant benefit of the 7-DOF beam element is its lower computational complexity. As shown in Table 1, the experimental bridge analysis using both the 7-DOF and classical beam elements required significantly fewer nodes and elements compared to the solid element model. The computational demands for the solid element model were substantially higher, reaffirming the time-saving benefits of the 7-DOF model in analyzing complex bridge structures.
Additionally, the 7-DOF beam element model is compatible with existing finite element programs that use classical beam elements, facilitating its integration. This compatibility means that existing node and element data can be used without the need to generate new data, requiring only modifications to the stiffness matrix and force vector.
For bridge design, particularly in limit state design, understanding internal forces along the girder is essential. Although solid and shell element models provide accurate stress and deflection data, they do not directly yield internal force information for each cross-section, often necessitating additional post-processing. The 7-DOF element model addresses this gap, offering a reliable method for both torsional analysis and the computation of internal forces across every cross-section.

7. Conclusions

In this research, the warping effect was included as an addition to the classical beam theory in the kinematics of the beam, and equations for computing the strain of the beam cross-sections were derived. With the application of total potential energy theory, a 7-DOF beam element, which had seven degrees of freedom at each node, was derived. A corresponding finite element computer program was also developed in MATLAB. The following conclusion can be drawn after comparing the 7-DOF beam element model with the solid element model, the classical beam element model, and the field monitoring results from the experimental bridge.
(a)
The 7-DOF beam element can accurately evaluate the torsional behaviors of bridge girders under construction load. The results were validated using the results from the solid element model and experimental data.
(b)
In terms of stress and vertical deflection, the 7-DOF beam element had the same behavior as the other numerical methods. However, according to the strain equations for the 7-DOF beam, the warping effect can lead to re-distribution of the stress on the cross-section, especially when the torsional moment applied to the structure is significant.
(c)
The classical beam element that is widely used in commercial software packages failed to compute the transverse rotation of bridge girders when the torsional moment was applied. The transverse rotation based on the results from the classical beam element was often larger than experimental data due to the lack consideration of warping stiffness during the torsional analysis.
(d)
With similar accuracy, the 7-DOF element model requires less modeling effort and computing time compared with the solid element model. For the experimental bridge in this study, the size of global stiffness matrix reduced to 1% when switching the solid element model to the 7-DOF beam element model.
(e)
Unlike the solid and shell element models, with which it is difficult to determine the internal forces directly, the 7-DOF beam element can compute the internal force of cross-sections along the bridge girder for design purposes.
(f)
The 7-DOF beam element model does not require changes in the nodal and element information of the classical beam element model. Therefore, it is possible to convert a classical beam element model into the 7-DOF beam element without much effort.
(g)
The numerical analysis using the 7-DOF beam element can be an alternative approach to the solid element and shell element for bridge analysis, especially when detailed information on stress distribution on the cross-section is not required and internal forces need to be generated.

Author Contributions

Experimental design and methodology, L.H., M.A. and R.H.; experimental testing and data analysis, L.H. and M.A.; numerical simulation, L.H. and M.A.; programming, L.H.; writing—original draft preparation, L.H.; writing—review and editing, L.H., M.A. and R.H.; project administration, R.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Illinois Department of Transportation, grant number ICT-R27-179.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

This publication is based on the results of ICT-R27-179, Effectiveness of Exterior Beam Rotation Prevention Systems for Bridge Deck Construction—Phase II. ICT-R27-179 was conducted in cooperation with the Illinois Center for Transportation; the Illinois Department of Transportation, Office of Program Development; and the US Department of Transportation, Federal Highway Administration.

Conflicts of Interest

Author Md Ashiquzzaman was employed by the company Ameren Corp. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Yang, S.; Helwig, T.; Klingner, R.; Engelhardt, M.; Fasl, J. Impact of Overhang Construction on Girder Design; University of Texas at Austin: Austin, TX, USA, 2009. [Google Scholar]
  2. Hui, L.; Hraib, F.; Gillis, B.; Vicente, M.; Hindi, R. A Simplified Method to Minimize Exterior Girder Rotation of Steel Bridges during Deck Construction. Eng. Struct. 2019, 183, 84–93. [Google Scholar] [CrossRef]
  3. Hui, L.; Hraib, F.; Vicente, M.; Hindi, R. Effectiveness of Exterior-Beam Rotation-Prevention Systems for Bridge Deck Construction—Phase II; Illinois Center for Transportation: Springfield, IL, USA, 2019. [Google Scholar]
  4. Ashiquzzaman, M.; Hui, L.; Schmeltz, J.; Merino, C.; Bozkurt, B.; Ibrahim, A.; Lindquist, W.; Hindi, R. Effectiveness of Exterior Beam Rotation Prevention Systems for Bridge Deck Construction; Illinois Center for Transportation: Springfield, IL, USA, 2016. [Google Scholar]
  5. Ashiquzzaman, M.; Calvo, C.M.; Hui, L.; Ibrahim, A.; Lindquist, W.; Hindi, R. Effectiveness of Different Bracing Systems to Prevent Exterior Girder Rotation during Bridge Deck Construction. Eng. Struct. 2017, 142, 272–289. [Google Scholar] [CrossRef]
  6. Ashiquzzaman, M.; Hui, L.; Ibrahim, A.; Lindquist, W.; Thomson, M.; Hindi, R. Effect of Inconsistent Diaphragms on Exterior Girder Rotation During Overhang Deck Construction. Structures 2016, 8, 25–34. [Google Scholar] [CrossRef]
  7. Ashiquzzaman, M.; Hui, L.; Hindi, R.; Ibrahim, A.; Lindquist, W. Validation of Field Exterior Girder Rotation in Non-Skewed Bridge Due to Construction Loads with FE Analysis. In Maintenance, Monitoring, Safety, Risk and Resilience of Bridges and Bridge Networks; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
  8. Fasl, J.D. The Influence of Overhang Construction on Girder Design; University of Texas at Austin: Austin, TX, USA, 2008. [Google Scholar]
  9. Shokouhian, M.; Shi, Y. Flexural Strength of Hybrid Steel I-Beams Based on Slenderness. Eng. Struct. 2015, 93, 114–128. [Google Scholar] [CrossRef]
  10. Winkler, R.; Kindmann, R.; Knobloch, M. Lateral Torsional Buckling Behaviour of Steel Beams—On the Influence of the Structural System. Structures 2017, 11, 178–188. [Google Scholar] [CrossRef]
  11. Gupta, V.K.; Okui, Y.; Nagai, M. Development of Web Slenderness Limits for Composite I-Girders Accounting for Initial Bending Moment. Struct. Eng./Earthq. Eng. 2006, 23, 229s–239s. [Google Scholar] [CrossRef]
  12. Jáger, B.; Dunai, L.; Kövesdi, B. Lateral-Torsional Buckling Strength of Corrugated Web Girders—Experimental Study. Structures 2022, 43, 1275–1290. [Google Scholar] [CrossRef]
  13. American Association of State Highway and Transportation Officials. AASHTO LRFD Bridge Design Specifications; American Association of State Highway and Transportation Officials: Washington, DC, USA, 2017. [Google Scholar]
  14. American Association of State Highway and Transportation Officials. Guidelines for Steel Girder Bridge Analysis, G13.1-2019; American Association of State Highway and Transportation Officials: Washington, DC, USA, 2019. [Google Scholar]
  15. Vlasov, V.Z. Thin-Walled Elastic Beams; National Technical Information Service: Haifa, Israel, 1984.
  16. Jönsson, J. Distortional Theory of Thin-Walled Beams. Thin-Walled Struct. 1999, 33, 269–303. [Google Scholar] [CrossRef]
  17. Kim, N.-I.; Kim, M.-Y. Exact Dynamic/Static Stiffness Matrices of Non-Symmetric Thin-Walled Beams Considering Coupled Shear Deformation Effects. Thin-Walled Struct. 2005, 43, 701–734. [Google Scholar] [CrossRef]
  18. Pavazza, R. Torsion of Thin-Walled Beams of Open Cross-Section with Influence of Shear. Int. J. Mech. Sci. 2005, 47, 1099–1122. [Google Scholar] [CrossRef]
  19. Roberts, T.M.; Al-Ubaidi, H. Influence of Shear Deformation on Restrained Torsional Warping of Pultruded FRP Bars of Open Cross-Section. Thin-Walled Struct. 2001, 39, 395–414. [Google Scholar] [CrossRef]
  20. Sapountzakis, E.J.; Mokos, V.G. Warping Shear Stresses in Nonuniform Torsion of Composite Bars by BEM. Comput. Methods Appl. Mech. Eng. 2003, 192, 4337–4353. [Google Scholar] [CrossRef]
  21. Prokić, A. Thin-Walled Beams with Open and Closed Cross-Sections. Comput. Struct. 1993, 47, 1065–1070. [Google Scholar] [CrossRef]
  22. Gruttmann, F.; Sauer, R.; Wagner, W. Shear Stresses in Prismatic Beams with Arbitrary Cross-Sections. Int. J. Numer. Meth. Engng. 1999, 45, 865–889. [Google Scholar] [CrossRef]
  23. Mentrasti, L. Torsion of Box Girders with Deformable Cross Sections. J. Eng. Mech. 1991, 117, 2179–2200. [Google Scholar] [CrossRef]
  24. Zhang, S.H.; Lyons, L.P.R. A Thin-Walled Box Beam Finite Element for Curved Bridge Analysis. Comput. Struct. 1984, 18, 1035–1046. [Google Scholar] [CrossRef]
  25. Cambronero-Barrientos, F.; Díaz-del-Valle, J.; Martínez-Martínez, J.-A. Beam Element for Thin-Walled Beams with Torsion, Distortion, and Shear Lag. Eng. Struct. 2017, 143, 571–588. [Google Scholar] [CrossRef]
  26. Dassault Systèmes. Abaqus Analysis User’s Manual (6.14); Dassault Systèmes: Vélizy-Villacoublay, France, 2014. [Google Scholar]
  27. Thompson, M.K.; Thompson, J.M. ANSYS Mechanical APDL for Finite Element Analysis; Elsevier Science: San Diego, CA, USA, 2017; ISBN 978-0-12-813110-7. [Google Scholar]
  28. ADINA R&D Inc. ADINA User’s Manual; ADINA R&D Inc.: Watertown, MA, USA, 2012. [Google Scholar]
  29. Murín, J.; Kutiš, V. An Effective Finite Element for Torsion of Constant Cross-Sections Including Warping with Secondary Torsion Moment Deformation Effect. Eng. Struct. 2008, 30, 2716–2723. [Google Scholar] [CrossRef]
  30. Sabat, L.; Kundu, C.K. Flexural-Torsional Analysis of Steel Beam Structures Using ANSYS. Mater. Today Proc. 2023. [Google Scholar] [CrossRef]
  31. Li, X.; Wan, S.; Mo, Y.L.; Shen, K.; Zhou, T.; Nian, Y. An Improved Method for Analyzing Shear Lag in Thin-Walled Box-Section Beam with Arbitrary Width of Cantilever Flange. Thin-Walled Struct. 2019, 140, 222–235. [Google Scholar] [CrossRef]
  32. Sapountzakis, E.J.; Mokos, V.G. 3-D Beam Element of Composite Cross Section Including Warping and Shear Deformation Effects. Comput. Struct. 2007, 85, 102–116. [Google Scholar] [CrossRef]
  33. Yau, J.D.; Kuo, S.-R. Geometrical Stiffness of Thin-Walled I-Beam Element Based on Rigid-Beam Assemblage Concept. J. Mech. 2012, 28, 97–106. [Google Scholar] [CrossRef]
  34. Bernuzzi, C.; Simoncelli, M. The Transformation Matrix in the 7DOFs Beam Formulation. Thin-Walled Struct. 2023, 190, 110951. [Google Scholar] [CrossRef]
  35. Ramaprasad, S.; Dattaguru, B.; Singh, G. Exact Solutions for Thin-Walled Composite Open Section Beams Using a Unified State Space Coupled Field Formulation. Mech. Based Des. Struct. Mach. 2023, 51, 5983–6007. [Google Scholar] [CrossRef]
  36. Lebastard, M.; Couchaux, M.; Bureau, A.; Hjiaj, M. Lateral-Torsional Buckling of Beams with Warping Restraints at Supports. ce/papers 2021, 4, 2262–2270. [Google Scholar] [CrossRef]
  37. Lebastard, M.; Couchaux, M.; Santana, M.V.B.; Bureau, A.; Hjiaj, M. Elastic Lateral-Torsional Buckling of Beams with Warping Restraints at Supports. J. Constr. Steel Res. 2022, 197, 107410. [Google Scholar] [CrossRef]
  38. Justín, M.; Juraj, H.; Mehdi, A.; Vladimír, K.; Juraj, P.; Stephan, K. Extension of the FGM Beam Finite Element by Warping Torsion. Strojnícky Časopis J. Mech. Eng. 2019, 69, 57–76. [Google Scholar] [CrossRef]
  39. Piotrowski, R.; Szychowski, A.; Vičan, J. Elastic Critical Resistance of the Simple Beam Grillage Resulting from the Lateral Torsional Buckling Condition: FEM Modelling and Analytical Considerations. Materials 2023, 16, 1346. [Google Scholar] [CrossRef] [PubMed]
  40. Burlon, A.; Failla, G.; Arena, F. Exact Stochastic Analysis of Coupled Bending–Torsion Beams with in-Span Supports and Masses. Probabilistic Eng. Mech. 2018, 54, 53–64. [Google Scholar] [CrossRef]
  41. Megson, T.H.G. Torsion of Solid Sections. In Aircraft Structures for Engineering Students; Elsevier: Amsterdam, The Netherlands, 2022; pp. 75–95. ISBN 978-0-12-822868-5. [Google Scholar]
  42. Mohammadi, E.; Hosseini, S.S.; Rohanimanesh, M.S. Elastic Lateral-Torsional Buckling Strength and Torsional Bracing Stiffness Requirement for Monosymmetric I-Beams. Thin-Walled Struct. 2016, 104, 116–125. [Google Scholar] [CrossRef]
  43. Roddis, W.M.K.; Kriesten, M.; Liu, Z. Torsional Analysis for Exterior Girders; University of Kansas Center for Research, Inc.: Lawrence, KS, USA, 1999; 128p. [Google Scholar]
  44. Cook, R.D. Concepts and Applications of Finite Element Analysis, 4th ed.; Wiley: New York, NY, USA, 2001; ISBN 978-0-471-35605-9. [Google Scholar]
Figure 1. Typical formwork system and construction loads during deck construction. (a) Deck overhang and construction equipment. (b) Overhang bracket. (c) Loads at location of screed machine.
Figure 1. Typical formwork system and construction loads during deck construction. (a) Deck overhang and construction equipment. (b) Overhang bracket. (c) Loads at location of screed machine.
Infrastructures 08 00175 g001
Figure 2. The 7-DOF beam element.
Figure 2. The 7-DOF beam element.
Infrastructures 08 00175 g002
Figure 3. Deformation field in beam element.
Figure 3. Deformation field in beam element.
Infrastructures 08 00175 g003
Figure 4. Framing plan of the bridge.
Figure 4. Framing plan of the bridge.
Infrastructures 08 00175 g004
Figure 5. Girder elevation of the bridge.
Figure 5. Girder elevation of the bridge.
Infrastructures 08 00175 g005
Figure 6. Cross-frame of the bridge.
Figure 6. Cross-frame of the bridge.
Infrastructures 08 00175 g006
Figure 7. Deck pouring sequence of the bridge.
Figure 7. Deck pouring sequence of the bridge.
Infrastructures 08 00175 g007
Figure 8. Predefined sections for instrumentation.
Figure 8. Predefined sections for instrumentation.
Infrastructures 08 00175 g008
Figure 9. Location of tilt sensors, strain gauge, and prism for the experimental bridge: (a) location of sensors at S1 and S3, (b) location of sensors at S2.
Figure 9. Location of tilt sensors, strain gauge, and prism for the experimental bridge: (a) location of sensors at S1 and S3, (b) location of sensors at S2.
Infrastructures 08 00175 g009
Figure 10. Tilt sensor and strain gauge: (a) Tilt sensor with open aluminum box. (b) Strain gauge with protective coating.
Figure 10. Tilt sensor and strain gauge: (a) Tilt sensor with open aluminum box. (b) Strain gauge with protective coating.
Infrastructures 08 00175 g010
Figure 11. Total station and prism: (a) Total station. (b) Prism.
Figure 11. Total station and prism: (a) Total station. (b) Prism.
Infrastructures 08 00175 g011
Figure 12. Exterior girder rotation at the bottom of the web on section S1.
Figure 12. Exterior girder rotation at the bottom of the web on section S1.
Infrastructures 08 00175 g012
Figure 13. Strain gauge data at section S2.
Figure 13. Strain gauge data at section S2.
Infrastructures 08 00175 g013
Figure 14. Vertical deflection at (a) section S1, (b) section S2, and (c) section S3.
Figure 14. Vertical deflection at (a) section S1, (b) section S2, and (c) section S3.
Infrastructures 08 00175 g014
Figure 15. Location of fresh concrete and screed machine used in finite element analysis: (a) fresh concrete location; (b) loads at section S1.
Figure 15. Location of fresh concrete and screed machine used in finite element analysis: (a) fresh concrete location; (b) loads at section S1.
Infrastructures 08 00175 g015
Figure 16. Solid element model created using ABAQUS/CAE.
Figure 16. Solid element model created using ABAQUS/CAE.
Infrastructures 08 00175 g016
Figure 17. Cross-frame modeled in ABAQUS/CAE: (a) cross-frame at non-support location; (b) cross-frame at support locations.
Figure 17. Cross-frame modeled in ABAQUS/CAE: (a) cross-frame at non-support location; (b) cross-frame at support locations.
Infrastructures 08 00175 g017
Figure 18. Classical beam element model built in SAP2000.
Figure 18. Classical beam element model built in SAP2000.
Infrastructures 08 00175 g018
Figure 19. Programming procedure for 7-DOF frame element model.
Figure 19. Programming procedure for 7-DOF frame element model.
Infrastructures 08 00175 g019
Figure 20. 7-DOF beam element model created in MATLAB.
Figure 20. 7-DOF beam element model created in MATLAB.
Infrastructures 08 00175 g020
Figure 21. Rotation comparison among different numerical models and experimental data for S1 section.
Figure 21. Rotation comparison among different numerical models and experimental data for S1 section.
Infrastructures 08 00175 g021
Figure 22. Rotation comparison among different numerical models and experimental data for S2 section.
Figure 22. Rotation comparison among different numerical models and experimental data for S2 section.
Infrastructures 08 00175 g022
Figure 23. Rotation comparison among different numerical models and experimental data for S3 section.
Figure 23. Rotation comparison among different numerical models and experimental data for S3 section.
Infrastructures 08 00175 g023
Figure 24. Vertical deflection comparison among different finite element models for S3 section.
Figure 24. Vertical deflection comparison among different finite element models for S3 section.
Infrastructures 08 00175 g024
Figure 25. Stress comparison among different finite element models.
Figure 25. Stress comparison among different finite element models.
Infrastructures 08 00175 g025
Table 1. Comparison of numerical models.
Table 1. Comparison of numerical models.
Solid Element Model (SEM)Classical Beam Element Model7-DOF Beam Element Model
Number of Nodes1,348,2886640 (0.49% of the nodes of SEM)6640 (0.49% of the nodes of SEM)
Number of Elements690,1527164 (1.04% of the nodes of SEM)7164 (1.04% of the nodes of SEM)
DOF at Each Node367
Global Stiffness Matrix Size4,044,86439,840 (0.98% of the nodes of SEM)46,480 (1.15% of the nodes of SEM)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hui, L.; Ashiquzzaman, M.; Hindi, R. Application and Experimental Validation of Seven-Degree-of-Freedom Beam Element for Girder Bridges during Deck Construction. Infrastructures 2023, 8, 175. https://doi.org/10.3390/infrastructures8120175

AMA Style

Hui L, Ashiquzzaman M, Hindi R. Application and Experimental Validation of Seven-Degree-of-Freedom Beam Element for Girder Bridges during Deck Construction. Infrastructures. 2023; 8(12):175. https://doi.org/10.3390/infrastructures8120175

Chicago/Turabian Style

Hui, Li, Md Ashiquzzaman, and Riyadh Hindi. 2023. "Application and Experimental Validation of Seven-Degree-of-Freedom Beam Element for Girder Bridges during Deck Construction" Infrastructures 8, no. 12: 175. https://doi.org/10.3390/infrastructures8120175

Article Metrics

Back to TopTop