1. Introduction
Considering the development of bridge structures using light thin-wall, long-span, wide-web spacings and wide cantilever plates, problems related to the spatial mechanics caused by thin-walled structures have become increasingly prominent [
1,
2,
3]. When the thin-walled box girder is in a flexural state, the flexural function of the wing plate is no longer applicable to the assumption of the plane section owing to the influence of shear deformation; the resulting warping phenomenon is called the “shear lag effect” [
4,
5,
6].
The analytical method based on the principle of energy variation is a research approach rooted in the theory of elasticity. Among various methods employed for analyzing the shear lag effect in box girders, the energy variation method stands out as the most commonly utilized. This method involves treating the difference between the longitudinal displacements of the girder body and the longitudinal displacements of the wing plates as unknowns. By applying the minimum potential energy theory to resolve the governing differential equations, the deflection of the girder body and the transverse distribution of normal stresses can be obtained. Wei et al. [
7] expressed the shear lag deflection function using the maximum shear deformation difference function, and combined it with the energy variation method to investigate the shear lag effect in curved box girders. Wang et al. [
8] assumed the shear lag deflection functions of box beams on the cantilever, top, and bottom plates. Based on spline finite points and the principle of energy variation, they discretized the structure using B-spline curve equations, thus revealing the variation law of shear lag effect in box girders. Chen et al. [
9], on the foundation of energy variation theory, introduced spatial girder element functions and derived the Euler girder matrix equation. By considering the influence of shear deformation, they derived the stiffness matrix of Timoshenko beams. Through extensive exploration by numerous scholars, the domain of energy variation methods has been further expanded. From rectangular box sections to curved box girder sections, and then to box girders utilizing composite materials, it has been discovered that employing the energy variation method to solve for the shear lag effect is indeed feasible. This paper also aims to investigate the shear lag effect in thin-walled box girders using the energy variation method. The key issue lies in accurately determining the form of the deflection function, as it directly affects the precision of solving for the shear lag effect.
Currently, the flexural displacement function is primarily based on artificial assumptions and has been considered quadratic parabola, cubic parabola, quadric parabola, cosine function, and catenary-type [
10,
11,
12,
13,
14] functions. Lin et al. [
15] reported that the shear lag effect of thin-walled box girders was essentially caused by shear deformation. By studying the shear flow of a box girder under a vertical load, they proposed a flexural displacement function with shear deformation. Penalba et al. [
16]. represented the shear strain of a steel–concrete composite girder with a cosine function and established a mechanical analysis model of a steel–concrete composite girder considering the influence of shear deformation and interface slip between the concrete and steel plate [
17]. The proposed model was verified by comparing the obtained results with solutions of elastic analysis and finite-element results. Zhang et al. [
18] introduced boundary correction factors to the cantilever plate to reflect the influence of different boundary constraints on the shear lag deflection. By utilizing spatial finite-element software to compute the results, the values of the boundary correction factors for the cantilever plate were determined. Many domestic and international scholars have also established sophisticated bridge finite-element models, verifying their simulation results through dual validation against experiments [
19,
20,
21,
22].
In the literature, the theoretical methods for flexural analysis of thin-walled box girders mainly focus on the shear lag effect of the box girders, whereas research on flexural and shear deformation is lacking [
23]. In particular, the theoretical correlation between the elementary beam theory and spatial mechanics analysis theory of box girders under the influence of shear deformation of each plate is limited. In this study, the flexural displacement function of a box girder under the influence of shear deformation of each plate is proposed, and the spatial mechanics analysis theory of box girders is improved. This paper will also conduct a comparative analysis of the deflection stress calculation results and propose a cantilever displacement correction for the bending deflection function of thin-walled box girders considering the effects of shear deformation in each plate. A solid model of a single-box three-chamber box girder will be established using MIDAS/FEA NX 2021 finite-element software to analyze the factors influencing the shear lag effect under different loading conditions. The theoretical solution for the shear lag effect in the bent box girder will be provided based on the energy variation method to analytically derive the deflection function. Factors affecting the shear lag effect will be introduced to modify the formula for the shear lag effect in a single-box three-chamber box girder, addressing the workload issue in internal force calculations for box girders.
The innovative aspects of this article include theoretically deriving the deflection function of thin-walled box girders under bending, influenced by the shear deformation of each plate, due to uneven shear flow transfer. This also leads to a theoretical analysis of wide-span box girders under the effect of uneven shear flow transfer. The theoretical solutions and influencing factors of the shear lag effect in bent box girders are solved using the energy variation method. Using Origin, an empirical formula for the shear lag effect in a single-box three-chamber girder is fitted, addressing the problem of calculating internal forces in box girders.
2. Analysis of Flexural Displacement Function of Flexural Box Girder under Shear Deformation of Plates
Considering the main girder with a single box and three chambers as the object of analysis, the force diagram of a thin-walled box girder under an arbitrary vertical load distribution is shown in
Figure 1. The section coordinates and dimensions of the box girder are shown in
Figure 2, where the center of the cross-section is the origin of the coordinate system, and the main inertial axes of the section are the coordinate axes.
In thin-walled box girders, stress is assumed to be uniformly distributed along the wall thickness direction. For the study of thin-walled box girders, the position of any point on the section can be represented by the coordinates of the flow curve along the centerline of the thickness of the box wall, and the section of the thin-walled box girder can be simplified. A simplified section of the thin-walled box girder is shown in
Figure 3, wherein the right-hand flow coordinate system is formed using
,
, and
;
represents the displacement of the longitudinal
-axis of the girder body, and
represents the displacement of the circumferential
-axis of the girder wall.
The following basic assumptions were made according to the flexural theory of a box girder [
24,
25]:
(1) Under small-deformation conditions, the contour projection of the box girder section does not deform.
(2) The strain at the midpoint of each plate is zero; that is, for each plate of the box girder, .
(3) The shear deformation outside the plane of each plate of the box girder is not considered, that is, .
In
Figure 3, 1 is the center point of the top plate of the inner chamber, 2 is the apex of the inner web, 3 is the apex of the outer web, 4 is the outer edge point of the cantilever plate, 5 is the center point of the bottom plate of the inner chamber, 6 is the bottom point of the inner web, 7 is the bottom point of the outer web, 8 is the intersection point of the centroidal
-axis and inner web, and 9 is the intersection of the centroidal
-axis and outer web.
To study the equilibrium state of the thin-walled member, microplates with a certain thickness were obtained from the cross-section of the thin-walled box girder, and the thin-walled microplates are shown in
Figure 4. As the structure is in a flexural state, normal and shear stresses exist in the microplates; according to the principle of shear stress reciprocity,
. In addition, because the thickness of the microplate in section
is constant, the shear force
of the unit length along the cross-section of the wall plate can be determined. This is called the shear flow; its flow direction is the same as that of
.
According to the mechanical equilibrium conditions, the resultant force of the thin-walled element shown in
Figure 4 is zero along the
-axis:
By simplifying Equation (1), we obtain
The basic assumption does not consider the out-of-plane shear deformation of each plate of the box girder, and
indicates that each plate is in a state of plane stress. According to the theory of elasticity, the following can be obtained:
The stress distribution in the thin-walled box girder under shear deformation is not clear; therefore, it cannot be calculated using the microelement equilibrium Equations (1)–(5). To solve this problem, we assume that the displacement function is the same as that of the elementary girder to obtain the flexural shear stress of the box girder using the microelement equilibrium equation; the assumed flexural displacement is then modified.
Under the assumption of a plane section, the equation for the flexural normal stress
of the box girder is as follows [
26]:
where
represents the bending moment of the box girder (N·m),
represents the distance from the transverse spindle of the box girder (m), and
represents the inertia distance of the box girder section with respect to the
-axis (m
4).
According to the flexural theory of thin-walled structures with openings [
27], Equation (6) is substituted into Equation (2) to simplify the expression of the shear flow [
28]:
where
represents the flexural shear flow and
represents the distance of the cross-section of the net area with respect to the
-axis.
Because the section of the box girder is a statically indeterminate structure
, to determine the initial shear flow at any starting point, a cut can be set at any section of the box girder to convert it into a statically indeterminate open-section problem for obtaining the initial value of the shear flow
. By superimposing the open-section shear flow
and additional shear flow
, the closed-section shear flow
can be obtained as follows:
The cross-sectional size of the single-box three-chamber box girder is shown in
Figure 5a; the shear center of the section is influenced by the bending moment
. Virtual incisions were made on the sections of each closed-box chamber. The locations of the virtual incisions are shown in
Figure 5b. The calculation equation for the shear flow in the open section is as follows:
Figure 5c presents the distribution of the shear flow
in the open sections of the three chambers in a single box girder. The arrow indicates the direction of the actual shear flow, which is obtained by combining the shear flow direction and the integral direction of the static surface moment.
According to a previous study [
29], if the incision selected is at the top and bottom plates on the sectional symmetry axis, the calculation of the shear flow in statically indeterminate structures can be obtained from the deformation compatibility conditions.
The deformation compatibility condition of chamber
is expressed as
where
represents the basic shear flow in the open section,
represents the additional shear flow in chamber
, and
represents the additional shear flow in the adjacent chamber.
The additional shear flow
in each chamber is solved according to the deformation compatibility condition, as shown in
Figure 5d. By superimposing
Figure 5c,d, the shear flow in the section of a single box and three chambers is obtained, and the corresponding shear flow distribution is shown in
Figure 5e.
Considering the top plate #1 of the box girder as an example, its internal shear deformation was analyzed, and its shear flow equation was obtained as follows:
where
represents the distance between the centerline of the top plate and the
-axis of the center of the box girder section, and
represents the additional shear flow in chamber 2.
When the box girder is bent, the displacement of the top plate along the perimeter
is zero, and
according to the basic assumption; thus,
. According to Equations (5) and (11),
where
represents the shear modulus of elasticity.
The equation for the longitudinal displacement of the top plate of the box girder #1 can be obtained by introducing
into Equation (12) as follows:
where
represents the longitudinal displacement of the center of the top plate of the box girder (starting point of coordinate
).
The expression of the longitudinal displacement at the inner web apex 2 (see
Figure 3) can be obtained via the simplified Equation (13) as follows:
if
where
represents the width of plate #1.
By substituting Equations (14) and (15) into Equation (13) and replacing
representing the horizontal position in the equation with the
coordinate, we obtain
Considering center point 5 of the bottom plate of the box girder (see
Figure 3) as the integral starting point of plate #4 (bottom plate of inner chamber), the equation for longitudinal displacement of plate #4 can be obtained as follows:
where
.
Similarly, the longitudinal displacement equations of the other wall plates can be obtained. The longitudinal displacement equations for the top and bottom plates are as follows:
According to Equation (18), the flexural displacements of the top plate of the inner chamber (plate#1), top plate of the outer chamber (plate#2) and cantilever plate (plate#3), bottom plate of the inner chamber (plate#4), bottom plate of the outer chamber (plate#5) can be represented by , , , and , respectively.
As shown in
Figure 5, the distribution of the shear flow on the box girder web is a quadratic function. By analyzing the shear flow on the box girder web, the shear flows
and
at points 8 and 9 (see
Figure 3) of the box girder web can be obtained as follows:
where
represents the web thickness and
represents the area of the wing plate corresponding to the open-section web.
According to Equations (19) and (20), the flexural shear flow at any point of the box girder web can be expressed as follows:
The web shear deformation equation obtained from Equations (5) and (21) is as follows:
According to the rigid peripheral assumption, if
is the vertical flexural displacement of the girder body, the relationship between the flexural angle of the box girder and the circumferential displacement of the web plate is as follows:
Equation (22) is substituted into Equation (23) and
is replaced with
to represent the position of any vertical point on the web. After simplification of the integral, the equation for the flexural longitudinal displacement at any point on the box girder web is obtained as follows:
By simplifying Equations (19), (20) and (24), , , , and are obtained. By combining Equations (18)–(20) and (24), the flexural longitudinal displacements of the box girder represented by and are obtained, the equations for which are presented below.
Inner chamber top and bottom plates and inner web:
Outer chamber top and bottom plates, cantilever plate, and outer web plate:
where
;
.
and
; the equation is as follows:
Equations (27) and (28) indicate that the flexural longitudinal displacement function of the box girder includes the longitudinal displacements and of the inner and outer webs at the -axis of the centroid and the intersection point. To accurately analyze the flexural longitudinal displacement function of the box girder under the influence of shear deformation of each plate, the longitudinal displacements and should be analyzed correctly.
Because the sum of the stresses on the box girder cross-section in the flexural state is zero, the following can be obtained:
According to Hooke’s law and Equation (3), the normal stress
at any point on the section of the flexural box girder can be expressed as follows:
By substituting Equations (25) and (26) into Equation (30), Equation (31) can be obtained as follows:
Equation (31) is simplified, and
decomposes into the following:
The static moment of the
-axis at the center of the box girder section as well as the integral results of Equations (32) and (33) are zero. If the sum of the stresses is zero, that is, if Equation (31) holds, the sum of Equations (34)–(38) is zero. According to this relationship, the following can be obtained:
By integrating both sides with respect to
z, the following is obtained:
where
is the integration constant. Because the girder body does not expand in the axial direction under a bending load,
;
is the shear lag warpage function that satisfies the axial balance parameter of Equation (31).
Here, represents the cross-sectional area of the top plate, represents the cross-sectional area of the bottom plate, and represents the cross-sectional area of the cantilever.
According to Equations (25), (26) and (40), under the influence of the shear deformation of the plate, the neutral axis of the box girder section no longer overlaps with the centroidal axis, and the deviation distance from the centroidal axis is .
According to the aforementioned description, after
and
are substituted into Equations (25) and (26) and the axial balance parameter
is substituted into
and
, the flexural displacement function of the flexural box girder under the influence of the shear deformation of each plate is expressed as follows:
In Equation (44),
is given as follows:
where
and
are the axial balance coefficients of the web, which are determined by the shear distribution law of the web and the ratio of
to
.
3. Analysis of Flexural Displacement Function of Box Girder via Energy Variational Method
According to the flexural displacement function of the flexural thin-walled box girder under the shear deformation of each plate, the energy variational principle [
30,
31] was applied to analyze the flexural displacement function of the thin-walled box girder. Using Equation (44), the longitudinal strain equation for each point of the thin-walled box girder section can be obtained as follows:
The shear strain of each wing plate is expressed as follows:
The shear strain of each web is expressed as follows:
According to the energy variational principle, the positive strain energy of each plate of the box girder under flexural deformation can be expressed as follows:
The flexural shear strain energy of each plate of the box girder is given as
By substituting Equations (46)–(48) into Equations (49) and (50), the total strain energy of the box girder under flexural deformation can be obtained as follows:
where
represents the cross-sectional warping moment of inertia,
represents the shear warping moment of the wing plate, and
represents the shear buckling moment of the web. The parameter equation of the shear warping sections defined in this study is as follows:
where
represents the area of the top plate,
represents the area of the bottom plate,
represents the area of the cantilever plate,
represents the area of the unilateral web, and
represents the full-section area of the box girder.
The potential energy
owing to the external load of the thin-walled box girder is expressed as follows:
From Equations (51) and (60), the total bending potential energy of the box girder can be obtained as follows:
According to the principle of the minimum potential energy, the stationary value
is substituted into Equation (61), yielding
Integration by parts, i.e.,
, is used to solve (62), and the following is obtained:
According to the principle of the variational method and Equation (63), the control differential equation for the flexural box girder can be obtained as follows:
The resulting boundary conditions are as follows:
By simplifying Equations (64) and (65), we obtain the following:
Here,
, and
and
are shear parameters.
By solving Equation (63), the general solution for the shear displacement is obtained as follows:
where
is the particular solution related to the shear force
, and
and
are the integration constants related to the boundary conditions.
According to Hooke’s law and Equation (61), the stress at any point on the box beam section can be expressed as follows when the influence of the shear deformation of the wing plates and webs is considered.
The equation for the shear lag coefficient
is as follows:
where
represents the actual stress on the section, and
represents the stress calculated using the elementary beam theory.
By substituting Equations (6) and (70) into Equation (71), the shear lag coefficients can be expressed as follows:
The bending curvature of the girder body is expressed as follows:
According to the boundary conditions, the integral constants and were solved to determine the expression of the shear displacement, and the theoretical solution of the flexural displacement function of the box girder was obtained by substituting the expression of the shear displacement into Equations (71)–(73).
5. Correction of Cantilever Plate Displacement for Flexural Displacement Function
Analyzing the finite-element numerical results of the box girder flexural space revealed that when the cantilever of the box girder cantilever plate was longer, the shear lag coefficients of the cantilever plate were significantly smaller than those of the top plate; the difference between the shear lag coefficients at the edge of the cantilever plate and those at the top plate was even larger. To reduce the effect of the stress concentration, the theoretically calculated and finite-element values of the shear lag coefficients in the section at position
of the wing plate were compared; the results are shown in
Figure 10.
As shown in
Figure 10, there was a large difference between the finite-element numerical simulation values at the center of the cantilever plate and the top plate of the box girder and the theoretically calculated values in this study. The shear lag coefficients at the edge of the cantilever plate were smaller than the shear lag coefficients at the center of the top plate, mainly because the edge of the cantilever plate and the center of the top plate had different constraints. Because the cantilever plate at the cantilever end of the constraints did not exist, the top plate of the box chamber on both sides of the web connection was constrained by the different boundary conditions manifested in the difference of the force characteristics. Therefore, it is necessary to optimize the flexural displacement function derived in this study by considering the effect of the shear deformation of each plate and correcting the differences caused by the different boundary constraints. Accordingly, the degree of buckling of the cantilever plate that is corrected using the correction parameter is denoted as
; Equation (45) is corrected as follows:
The final form of the flexural displacement function of the bent box girder under the influence of the shear deformation of each plate can be obtained by substituting Equation (74) into Equation (39) and applying the axial equilibrium condition. The detailed process is not described herein.
Owing to the influence of the aforementioned parameter
, the axial equilibrium parameter in Equation (42) can be corrected as follows:
In this study, the energy variational method was utilized to analyze the flexural displacement function of the box girder. The expressions of the cross-sectional parameters Equations (49)–(59) in the established control differential equations for the corrected cross-sectional parameters, considering the effect of the correction parameter
, are as follows:
Here, represents the area of the top plate, represents the area of the bottom plate, represents the area of the cantilever plate, represents the area of one side of the web, and represents the full cross-sectional area of the box girder.
The equations for the cross-sectional parameters not listed are similar to Equations (53)–(57).
To obtain more accurate values of the shear lag coefficients, the value of the correction parameter
of the cantilever plate must be determined by comparing the results of the numerical analysis calculations obtained via MIDAS/FEA NX and those obtained via the theoretical method for different values of
, as well as the calculation of the upper-left wing plate of the finite-element model of a single-box, three-chamber box girder (
). As shown in
Figure 10, the parameter
must be >1; thus, it was set as 1.0, and increased to 1.7 in steps of 0.1. The shear lag coefficients of the flexural box girder were corrected and calculated, and the results are presented in
Table 3.
As shown in
Table 3, when the value of the cantilever modification parameter
is 1.6, according to the flexural displacement function derived in this study, the calculated shear lag coefficient values of the wing plate on the box girder at each point optimally coincide with the numerical simulation results of MIDAS/FEA NX and fully reflect the distribution law, indicating that the shear lag coefficient of the cantilever plate is smaller than that of the top plate. To more intuitively examine the influence of the correction parameter
, the shear lag coefficient values of the wing plate on the box girder at the
section were compared, as shown in
Figure 11.
Comparing the stress distributions shown in
Figure 10 and
Figure 11 reveals that when the cantilever plate correction parameter
had a value of 1.6, the stresses on the upper-wing plate of the box girder calculated via the proposed theoretical method and the results obtained via the MIDAS/FEA NX numerical simulation agreed well, and the relative difference in the data was <2.00%. The proposed theoretical method can fully reflect the different constraint conditions of the cantilever plate and the top plate, as well as the transverse distribution of the stresses along the wing plate.
7. Conclusions
For the thin-walled box girder under the action of the shear deformation of each plate, the flexural displacement function of the wing and web plates of the box girder cannot satisfy the assumption of the plane section; the longitudinal flexural displacement function in the transverse direction of the wing plate is a quadratic function, and the height of the web plate is the sum of the primary and cubic functions.
According to the theoretical analysis of the axial equilibrium conditions of the thin-walled box girder subjected to bending under the influence of the shear deformation of each plate, the actual neutral axis of the box girder under the influence of shear deformation no longer coincides with the horizontal axis of the box girder cross-section, having an offset distance of .
The measured values of the model test, the simulated values of the finite-element model, and the calculated values of the theoretical method were compared, revealing sufficient agreement among the results. The majority of the shear lag coefficients had a relative error of <10.00%, confirming that the flexural displacement function of the flexural box girder under the influence of the shear deformation of the various plates derived in this study is reliable.
Owing to the different constraint conditions at the edge of the cantilever plate and at the center of the top plate, the cantilever plate and top plate exhibit differences in their two-dimensional force characteristics and boundary constraints. The length of the cantilever plate should be appropriately reduced in engineering design to weaken its influence on the shear lag effect of the box girder.
The article investigates the theoretical solution of shear hysteresis effects in bent box girders based on the energy method, along with the influencing factors of shear hysteresis effects. An empirical formula for shear hysteresis effects in single-box three-chamber box girders is fitted using Origin, addressing the issue of workload in calculating internal forces of box girders. In future developments, the following issues will continue to be addressed:
During the derivation of the deflection function for thin-walled box girders, local constraints such as bulkheads and fillet welds were not considered. In future work, it is necessary to conduct a thorough analysis of the effects of these local constraints.
As box girder structures continue to evolve towards wider cantilevered plates, larger flange-rib spacings, and lighter thin-walled materials, the nonlinear effects in geometric structure and material properties will become increasingly significant. This will become a major research area for thin-walled box girders in the future, requiring further study to refine the theoretical calculations for thin-walled box girders.