Next Article in Journal
A Two-Stage Automatic Container Code Recognition Method Considering Environmental Interference
Next Article in Special Issue
Effect of Subsequent Subgrade on Seismic Response of the High-Speed Railway Track–Bridge System
Previous Article in Journal
The Use of Virtual Reflectance Transformation Imaging (V-RTI) in the Field of Cultural Heritage: Approaching the Materiality of an Ancient Egyptian Rock-Cut Chapel
Previous Article in Special Issue
Investigation of the Mechanical Features of Steel–Concrete Composite Girder Rigid Frame Bridges with V-Shaped Piers during Construction Stages
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Shear Lag Effect of Ultra-Wide Box Girder under Influence of Shear Deformation

1
School of Transportation and Geomatics Engineering, Shenyang Jianzhu University, Shenyang 110168, China
2
School of Resources and Civil Engineering, Northeastern University, Shenyang 110819, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(11), 4778; https://doi.org/10.3390/app14114778
Submission received: 23 March 2024 / Revised: 13 May 2024 / Accepted: 28 May 2024 / Published: 31 May 2024
(This article belongs to the Special Issue Bridge Structural Analysis)

Abstract

:
The objective of this study was to determine the reasonable flexural functions of ultra-wide box girders, reveal the mechanism of the shear lag effect, and improve the analysis theory of ultra-wide box girders. Considering a single-box three-chamber thin-walled box girder as an example, starting from an uneven transfer of shear flow, the flexural displacement function of the curved box girder under the influence of shear deformation of each plate was derived according to the flexural theory of a thin-walled box girder, balance equation of a thin-walled microelement plate, and theory of plane stress. The energy variational method was used to analyze the flexural displacement function, providing a theoretical solution for the shear lag effect of the curved box girder. A displacement correction of the cantilever plate displacement function was performed by comparing the calculation results for the shear lag coefficients. The results indicated that under the shear deformation of each plate, the flexural displacement functions of the wing and web plates of the box girder no longer satisfy the assumption of plane section. The flexural displacement function is a quadratic function of the transverse wing plate, and the web height is the sum of the first- and third-order functions. The theoretical calculation results agree sufficiently well with the experimental results, proving that the flexural displacement function of the box girder under the influence of the shear deformation of each plate is reliable.

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 n , s , and z ; u represents the displacement of the longitudinal z -axis of the girder body, and v represents the displacement of the circumferential s -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 ε s at the midpoint of each plate is zero; that is, for each plate of the box girder, v / s = 0 .
(3) The shear deformation outside the plane of each plate of the box girder is not considered, that is, τ n z = 0 , σ n = 0 .
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 x -axis and inner web, and 9 is the intersection of the centroidal x -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, τ s z = τ z s = τ . In addition, because the thickness of the microplate in section d z is constant, the shear force q = τ t 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 s .
According to the mechanical equilibrium conditions, the resultant force of the thin-walled element shown in Figure 4 is zero along the z -axis:
σ z + σ z z d z t d s σ z t d s + q + q s d s d z q d z = 0 .
By simplifying Equation (1), we obtain
σ z z t + q s = 0 .
The basic assumption does not consider the out-of-plane shear deformation of each plate of the box girder, and τ n z = 0 , σ n = 0 indicates that each plate is in a state of plane stress. According to the theory of elasticity, the following can be obtained:
u z = ε z
v s = ε s
γ s z = u z + v s .
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 σ z of the box girder is as follows [26]:
σ z = M z I x y ,
where M z represents the bending moment of the box girder (N·m), y represents the distance from the transverse spindle of the box girder (m), and I x represents the inertia distance of the box girder section with respect to the x -axis (m4).
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]:
q s = M z I x S y ,
where q s represents the flexural shear flow and S y represents the distance of the cross-section of the net area with respect to the y -axis.
Because the section of the box girder is a statically indeterminate structure q ¯ , 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 q ¯ . By superimposing the open-section shear flow q 0 and additional shear flow q i , the closed-section shear flow q can be obtained as follows:
q s = q 0 + q i .
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 M z . 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:
q 0 = M z I x S y = M z I x y t d s .
Figure 5c presents the distribution of the shear flow q 0 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 i is expressed as
i q 0 t d s + q i i d s t q k i . k d s t = 0 ,
where q 0 represents the basic shear flow in the open section, q i represents the additional shear flow in chamber i , and q k represents the additional shear flow in the adjacent chamber.
The additional shear flow q i 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:
q s = M z I x y s t s + q 2 ,
where y s represents the distance between the centerline of the top plate and the x -axis of the center of the box girder section, and q 2 represents the additional shear flow in chamber 2.
When the box girder is bent, the displacement of the top plate along the perimeter s is zero, and v z = 0 according to the basic assumption; thus, v / z = 0 . According to Equations (5) and (11),
γ s z = u s = τ s z G = q s G t = M z y s G I x s + q 2 G t ,
where G 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 s into Equation (12) as follows:
u s , z = u 0 s , z + u i s , z = M z y s 2 G I x s 2 + q 2 G t s + u 01 z ,
where u 01 ( z ) represents the longitudinal displacement of the center of the top plate of the box girder (starting point of coordinate s ).
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:
u 02 z = M z y s 2 G I x b 1 2 + u 01 z
if
g 01 z = u 02 z u 01 z = M z y s 2 G I x b 1 2 ,
where b 1 represents the width of plate #1.
By substituting Equations (14) and (15) into Equation (13) and replacing s representing the horizontal position in the equation with the x coordinate, we obtain
u x , z = u 02 z 1 x 2 b 1 2 g 1 z β 1 z 1 x b 1
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:
u x , z = u 06 z 1 x 2 b 1 2 g 1 z β 1 z 1 x b 1 ,
where β 1 z = q 2 G t b 1 .
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:
u x , y , z = u 2 z y s b 1 2 x 2 g z 2 q 2 t b 1 x g z Plate # 1 u 6 z + y x b 1 2 x 2 g z 2 q 2 t b 1 x g z Plate # 4 u 3 z y s 4 b 2 2 b 1 + 2 b 2 x 2 g z 2 q 1 t 2 b 2 b 1 + 2 b 2 x g z Plate # 2 u 7 z + y x 4 b 2 2 b 1 + 2 b 2 x 2 g z 2 q 1 t 2 b 2 b 1 + 2 b 2 x g z Plate # 5 u 3 z y s b 3 2 b 1 + 2 b 2 + b 3 x 2 g z Plate # 3
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 u 2 z , u 3 z , u 6 z , and u 7 z , 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 q 8 z and q 9 z at points 8 and 9 (see Figure 3) of the box girder web can be obtained as follows:
q 8 z = q 08 z + q 1 q 2 = M z I x y s A s 1 + M z 2 I x y s 2 t w + q 1 q 2
q 9 z = q 09 z + q 1 = M z I x y s A s 2 + M z 2 I x y s 2 t w + q 1 ,
where t w represents the web thickness and A s i 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:
q y , z = q 8 , 9 z M z 2 I x y 2 t w .
The web shear deformation equation obtained from Equations (5) and (21) is as follows:
γ s z = u s + v z = q z G t w = q 8 , 9 z G t w M z 2 G I x y 2 .
According to the rigid peripheral assumption, if w z 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:
w z = v z .
Equation (22) is substituted into Equation (23) and s is replaced with y 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:
u ( y , z ) = y w ( z ) + ( q 08 , 09 ( z ) G t w + q i G t w ) y M ( z ) 6 G I x y 3 + u 8 , 9 ( z ) .
By simplifying Equations (19), (20) and (24), u 2 z , u 3 z , u 6 z , and u 7 z are obtained. By combining Equations (18)–(20) and (24), the flexural longitudinal displacements of the box girder represented by u 8 z and u 9 z are obtained, the equations for which are presented below.
Inner chamber top and bottom plates and inner web:
u x , y , z = y w z + α 1 y 1 3 y 3 + η 1 x , y g z + u 8 z
Outer chamber top and bottom plates, cantilever plate, and outer web plate:
u x , y , z = y w z + α 2 y 1 3 y 3 + η 2 x , y g z + u 9 z
where α 1 = 2 y s A s 1 t w + y s 2 + 2 q 1 q 2 t w ; α 2 = 2 y s A s 2 t w + y s 2 + 2 q 1 t w .
η 1 x , y and η 2 x , y ; the equation is as follows:
η 1 x , y = y s b 1 2 x 2 2 q 2 t b 1 x Plate # 1 y x b 1 2 x 2 2 q 2 t b 1 x Plate # 4
η 2 x , y = y s 4 b 2 2 b 1 + 2 b 2 x 2 2 q 1 t 2 b 2 b 1 + 2 b 2 x Plate # 2 y s b 3 2 b 1 + 2 b 2 + b 3 x 2 Plate # 3 y x 4 b 2 2 b 1 + 2 b 2 x 2 2 q 1 t 2 b 2 b 1 + 2 b 2 x Plate # 5 .
Equations (27) and (28) indicate that the flexural longitudinal displacement function of the box girder includes the longitudinal displacements u 8 z and u 9 z of the inner and outer webs at the x -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 u 8 z and u 9 z 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:
F N = A σ z d A = 0 .
According to Hooke’s law and Equation (3), the normal stress σ ( x , y , z ) at any point on the section of the flexural box girder can be expressed as follows:
σ x , y , z = E ε z = E u z .
By substituting Equations (25) and (26) into Equation (30), Equation (31) can be obtained as follows:
F N = E A 1 y w z + α 1 y 1 3 y 3 + η 1 ( x , y ) g z + u 8 z d A + E A 1 y w z + α 2 y 1 3 y 3 + η 2 x , y g z + u 9 z d A = 0 .
Equation (31) is simplified, and F N decomposes into the following:
F N 1 = w z E A y d A
F N 2 = α i g z E A y d A
F N 3 = 1 3 g z E A y 3 d A
F N 4 = E g z A 1 η 1 x , y d A
F N 5 = E g z A 2 η 2 x , y d A
F N 6 = u 8 z E A 1 d A
F N 7 = u 9 z E A 2 d A .
The static moment of the x -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:
A 1 u 8 z + A 2 u 9 z = A 1 η 1 x , y d A + A 2 η 2 x , y d A + 1 3 A y 3 d A g z .
By integrating both sides with respect to z, the following is obtained:
A 1 u 8 z + A 2 u 9 z = d g z + C ,
where C is the integration constant. Because the girder body does not expand in the axial direction under a bending load, C = 0 ; d is the shear lag warpage function that satisfies the axial balance parameter of Equation (31).
d = d e + d w .
d e = 2 y s A t y x A b b 1 2 + 4 b 2 2 + A t A b 2 q 2 t b 1 + 2 q 1 t 2 b 2 + y s A c b 3 2 3 A
d w = y x 3 A b y s 3 A t + A c + t w 2 y x 4 y s 4 3 A
Here, A t represents the cross-sectional area of the top plate, A b represents the cross-sectional area of the bottom plate, and A c 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 d × g z .
According to the aforementioned description, after u 8 z and u 9 z are substituted into Equations (25) and (26) and the axial balance parameter d is substituted into η 1 x , y and η 2 x , y , the flexural displacement function of the flexural box girder under the influence of the shear deformation of each plate is expressed as follows:
u x , y , z = y w z + α i y 1 3 y 3 + η ζ x , y g z .
In Equation (44), η ζ x , y is given as follows:
η ζ x , y = y s b 1 2 x 2 2 q 2 t b 1 x + k 1 d Plate # 1 y x b 1 2 x 2 2 q 2 t b 1 x + k 1 d Plate # 4 y s 4 b 2 2 b 1 + 2 b 2 x 2 2 q 1 t 2 b 2 x b 1 + k 2 d Plate # 2 y x 4 b 2 2 b 1 + 2 b 2 x 2 2 q 1 t 2 b 2 x b 1 + k 2 d Plate # 5 y s b 3 2 b 1 + 2 b 2 + b 3 x 2 + k 2 d Plate # 3 k 1 d Inner   web k 2 d Outer   web .
where k 1 and k 2 are the axial balance coefficients of the web, which are determined by the shear distribution law of the web and the ratio of A 1 to A 2 .

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:
ε z = u z = y ω z + α y 1 3 y 3 + η ζ x , y g z = y ω z + η x , y g z .
The shear strain of each wing plate is expressed as follows:
γ x z = u x = η ζ x , y x g z .
The shear strain of each web is expressed as follows:
γ y z = u s + v z = u y + ω z = η x , y y g z .
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:
U 1 = 1 2 σ z ε z d A d z = E 2 ε z 2 d A d z .
The flexural shear strain energy of each plate of the box girder is given as
U 2 = 1 2 τ x z γ x z + τ y z γ y z d A d z = G 2 γ x z 2 + γ y z 2 d A d z
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:
U = 1 2 l A E ε 2 + G γ 2 d A d z = 1 2 l E I x ω 2 + E I η g 2 E I y η ω g + G A x + A y g 2 d z
where I η represents the cross-sectional warping moment of inertia, A x represents the shear warping moment of the wing plate, and A y represents the shear buckling moment of the web. The parameter equation of the shear warping sections defined in this study is as follows:
A x = A η ζ x 2 d A = 1 3 { A t y s 2 b 1 + 2 q 2 / t 2 + 4 b 2 y s + 2 q 1 / t 2 + A b y x 2 b 1 + 2 q 2 / t 2 + 4 b 2 y x + 2 q 1 / t 2 + A c y s 2 b 3 2 } .
A y = A η y 2 d A = α 2 A w + I 4 w α I 2 w
I η = A η x , y 2 d A = α 2 I x + 1 9 I 6 + I ζ 2 3 α I 4 3 α I 1 ζ + I 3 ζ
I y η = A y η x , y d A = α I x 1 3 I 4 + I 1 ζ
I n = A y n d A = y s n A t + A c + y x n A b + 4 t w y s n + 1 + y x n + 1 n + 1 , n = 4 , 6
I n w = A y n d A w = 4 t w y s n + 1 + y x n + 1 n + 1 , n = 2 , 4
I ζ = A η ζ x , y 2 d A = A 1 k 1 2 + A 2 k 2 2 d 2 + A t 8 15 y s 2 b 1 4 + q 2 b 1 t 2 + y s 2 4 b 2 2 2 + 2 q 1 b 2 t 2 + 2 q 3 b 2 t 2 4 3 y s k 1 d b 1 2 + y s k 2 d 4 b 2 2 4 y s q 2 b 1 3 t 8 y s q 1 b 2 3 t 8 y s q 3 b 2 3 t + 4 k 1 d q 2 b 1 t + 4 k 2 d q 1 b 2 t + 4 k 2 d q 3 b 2 t + A b 8 15 y x 2 b 1 4 + q 2 b 1 t 2 + y x 2 4 b 2 2 2 + 2 q 1 b 2 t 2 + 2 q 3 b 2 t 2 + 4 3 y x k 1 d b 1 2 + y x k 2 d 4 b 2 2 4 y x q 2 b 1 3 t 8 y x q 1 b 2 3 t 8 y x q 3 b 2 3 t + 4 k 1 d q 2 b 1 t + 4 k 2 d q 1 b 2 t + 4 k 2 d q 3 b 2 t + A c y s 8 15 y s b 3 2 4 3 k 2 d b 3 2
I n ζ = A y n η ζ x , y d A = 2 3 y s n + 1 A t b 1 2 + 4 b 2 2 + A c b 3 2 + A b b 1 2 + 4 b 2 2 y x y s n + 1 , n = 1 , 3 .
where A t represents the area of the top plate, A b represents the area of the bottom plate, A c represents the area of the cantilever plate, A w represents the area of the unilateral web, and A represents the full-section area of the box girder.
The potential energy V owing to the external load of the thin-walled box girder is expressed as follows:
V = l M z ω z d z .
From Equations (51) and (60), the total bending potential energy of the box girder can be obtained as follows:
= U + V = 1 2 l E I x ω 2 + E I η g 2 2 E I y η ω g + G A x + A y g 2 d z l M z ω z d z .
According to the principle of the minimum potential energy, the stationary value δ = 0 is substituted into Equation (61), yielding
δ = l E I x ω δ ω + E I η g δ g E I y η ω δ g + g δ ω + G A x + A y g δ g + M δ ω d z = l E I x ω E I y η g + M δ ω + E I η g E I y η ω δ g + G A x g + A y g δ g d z
Integration by parts, i.e., a b u d v = u v a b a b v d u , is used to solve (62), and the following is obtained:
δ = l E I x ω E I y η g + M δ ω + E I η g E I y η ω δ g + G A x g + A y g δ g d z = l E I x ω E I y η g + M δ ω + G A x g + A y g E I η g + E I y η ω δ g d z + E I η g E I y η ω δ g 0 l .
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:
E I x ω E I y η g + M = 0
G A x g + A y g E I η g + E I y η ω = 0 .
The resulting boundary conditions are as follows:
E I η g E I y η ω δ g 0 l = 0 .
By simplifying Equations (64) and (65), we obtain the following:
g k 2 g = k 0 Q z
ω = M E I x + B g .
Here, B = I y η / I x , and k = G A x + A y I x E I x I η I y η 2 and k 0 = I y η E I x I η I y η 2 are shear parameters.
By solving Equation (63), the general solution for the shear displacement is obtained as follows:
g z = k 0 C 1 sinh k z + C 2 cosh k z + g * ,
where g * is the particular solution related to the shear force Q z , and C 1 and C 2 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.
σ x , y , z = E y ω z + η x , y g z = y M z I x + η x , y y I y η I x E g z
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:
λ x , y , z = 1 + η x , y y I y η I x E I x g z M z .
The bending curvature of the girder body is expressed as follows:
ω = M z E I x + I y η I x g .
According to the boundary conditions, the integral constants C 1 and C 2 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).

4. Experimental and Case Studies

To validate the displacement function for the flexural strength of the box girders under the influence of the shear deformation in each plate, the results of the theoretical method used in this study were compared with finite-element calculations and experimental measurements for the physical model of an organic-glass, single-box, three-chamber, simply supported box girder [32].
Figure 6 presents the structural dimensions, measurement-point arrangement, and loading effects of the aforementioned Plexiglas simply supported box girder test model [32]. The box girder test model had a span of L = 1000   mm . Its elastic modulus and Poisson’s ratio, which were measured from a test stretch, were E = 2600   MPa and ν = 0.385 , respectively. The boundary condition was simulated by supporting steel columns of the same gauge underneath the bottom plates corresponding to all the web ends of the girder. Loads were applied on top of all the webs of the L / 2 span section of the single-box, three-compartment, simply supported box girder using a loading system. A computer-controlled voltage jack with a load value of P / 4 = 0.4   kN was used for the centralized force. Photographs of the loading system are shown in Figure 7.
The MIDAS/FEA NX 2021 spatial finite-element software was used to build a solid simulation model of a single-box, three-chamber, simply supported box girder. The cross-section of the finite-element model was a single-box, three-chamber box section. The dimensions and arrangement of the measurement points are shown in Figure 8 and the space finite- element model are shown in Figure 9. Model tests used 120 resistance strain gauges to measure the longitudinal strain of the box girder under load. A total of 24, 15, and 9 strain gauges were, respectively, arranged on the cross-section, top plate, and bottom plate, and they were affixed to the surface of the box girder using structural adhesive. The measurement points of the finite-element model and locations of the measurement points of the experimental model had a correspondence; the span of the model was L = 30000   mm , and the material was C30 concrete with a Young’s modulus of E = 3.0 × 10 4   MPa and Poisson’s ratio of ν = 0.2 . In the finite-element model of the three-chamber box girder, three-dimensional hexahedral solid elements were used to simulate the main girder. The bridge had a total of 234,967 nodes and 157,395 solid units; the centralized load was applied at the top of the 4 webs of the cross-section of the finite-element model of the L / 2 span, and the value of the load was P / 4 = 360   kN .
According to the similarity principle [34], the cross-sectional stiffnesses corresponding to the finite-element numerical model and the experimental model were similar, as well as their cross-sectional geometries; the load of the finite-element model was proportional to that of the model test, and the model boundary conditions were identical. The results of the theoretical method used in this study are presented in Table 1 and Table 2 in comparison with the finite-element calculation results and the measured values of the model test reported in [32].
As shown in Table 1 and Table 2, the theoretically derived flexural displacement function of the thin-walled box girder subjected to bending under the influence of the shear deformation of each plate in this study agreed sufficiently with the measured values of the test and the finite-element numerical simulation values. Most key points of the shear lag coefficient had a relative error of <10.00%. There are several reasons for the occurrence of result discrepancies: firstly, differences in the material properties between the experimental model and the material used in the finite-element numerical simulation; secondly, inherent errors in the experimental setup; thirdly, differences in boundary constraints between the cantilevered plate and the top plate in the model.

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 y = 1450   cm 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:
η ζ x , y = y s b 1 2 x 2 2 q 2 t b 1 x + k 1 d   Inner   chamber   top   plate y s 4 b 2 2 b 1 + 2 b 2 x 2 2 q 1 t 2 b 2 b 1 + 2 b 2 x + k 1 d Outer   chamber   top   plate y s b 3 2 b 1 + 2 b 2 + b 3 x 2 φ + k 1 d   Cantilever   plate y x b 1 2 x 2 2 q 2 t b 1 x + k 1 d Inner   chamber   bottom   plate y x 4 b 2 2 b 1 + 2 b 2 x 2 2 q 1 t 2 b 2 b 1 + 2 b 2 x + k 1 d Outer   chamber   bottom   plate k 1 d Inner web   k 2 d Outer   web
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:
d e = 2 y s A t y x A b b 1 2 + 4 b 2 2 + A t A b 2 q 2 t b 1 + 2 q 1 t 2 b 2 + y s φ A c b 3 2 3 A .
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:
A x = A η ζ x 2 d A = 1 3 { A t y s 2 b 1 + 2 q 2 / t 2 + 4 b 2 y s + 2 q 1 / t 2 +                     A b y x 2 b 1 + 2 q 2 / t 2 + 4 b 2 y x + 2 q 1 / t 2 + A c φ y s 2 b 3 2 }
I ζ = A η ζ x , y 2 d A = A 1 k 1 2 + A 2 k 2 2 d 2 + A t 8 15 y s 2 b 1 4 + q 2 b 1 t 2 + y s 2 4 b 2 2 2 + 2 q 1 b 2 t 2 + 2 q 3 b 2 t 2 4 3 y s k 1 d b 1 2 + y s k 2 d 4 b 2 2 4 y s q 2 b 1 3 t 8 y s q 1 b 2 3 t 8 y s q 3 b 2 3 t + 4 k 1 d q 2 b 1 t + 2 k 2 d q 1 b 2 t + 2 k 2 d q 3 b 2 t + A b 8 15 y x 2 b 1 4 + q 2 b 1 t 2 + y x 2 4 b 2 2 2 + 2 q 1 b 2 t 2 + 2 q 3 b 2 t 2 + 4 3 y x k 1 d b 1 2 + y x k 2 d 4 b 2 2 4 y x q 2 b 1 3 t 8 y x q 1 b 2 3 t 8 y x q 3 b 2 3 t + 4 k 1 d q 2 b 1 t + 2 k 2 d q 1 b 2 t + 2 k 2 d q 3 b 2 t + A c φ y s 8 15 y s b 3 2 4 3 k 2 d b 3 2
I n ζ = A y n η ζ x , y d A = 2 3 y s n + 1 y s n + 1 A t b 1 2 + 4 b 2 2 + A c φ b 3 2 + A b b 1 2 + 4 b 1 2 y x y s n + 1 , n = 1 , 3
Here, A t represents the area of the top plate, A b represents the area of the bottom plate, A c represents the area of the cantilever plate, A w represents the area of one side of the web, and A 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 ( x 0 ). 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 y = 1450   cm 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.

6. Discussions

Based on the experimental model of an acrylic single-box three-chamber box beam and a solid finite-element model established using MIDAS/FEA NX, stress and deflection at key points of sections L / 2 and L / 4 of the simply supported box beam under different loading conditions were analyzed. The results of the normal stress and beam deflection obtained by solving the deflection function using the energy variational method were compared and analyzed against the finite-element numerical simulation results and experimental measurements. By comparing the stress results, significant discrepancies were found at the edge of the cantilevered plate when using the energy variational method to solve the deflection function. After correcting for the displacement of the cantilevered plate in consideration of the shear deformation effects on each plate, the results were compared with finite-element analysis. The conclusion is as follows:
(1) Comparing the experimental measurements, finite-element simulation values, and theoretical method calculations, the data results show good agreement, with the vast majority of data discrepancies within 10.00%. This confirms the reliability of the deflection function derived in this paper for the bending box beam considering the shear deformation effects on each plate.
(2) Due to different boundary constraints at the edge of the cantilevered plate and at the center of the top plate, the cantilevered plate and the top plate exhibit differences in two-dimensional stress characteristics and boundary constraints. By calculating the boundary conditions of the cantilevered plate for the deflection function of the bending box beam considering the shear deformation effects on each plate, when the cantilevered plate correction parameter is φ = 1.6 , the theoretical method aligns most closely with the finite-element simulation results, with a relative discrepancy of less than 2.00%.

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 d × g z .
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.

Author Contributions

Conceptualization, methodology, software, and validation, formal analysis, investigation, resources, data curation, Y.L. (Yanfeng Li) and J.X.; writing—original draft preparation, writing—review and editing, visualization, Y.L. (Yanfeng Li) and F.W.; supervision, project administration, funding acquisition, F.W. and Y.L. (Yuanhui Li). All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by National Key R&D Program of China (Nos. 2018YFC0809600 and 2018YFC0809606); Natural Science Foundation of Liaoning Province (No. 2019-MS-265); and Key Laboratory of the Ministry of Education on Safe Mining of Deep Metal Mines (No. 2019SMDM-KF-A01); Shenyang Municipal Science and Technology Planning Project Fund (No. 23-407-3-19).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Zhou, S. Effect of Shear Lag on Structural Behaviours in Indeterminate Box-Girders. J. Civ. Archit. Environ. Eng. 2010, 32, 7–11. [Google Scholar] [CrossRef]
  2. Luo, Z.; Song, Y.; He, S.; Huang, F. Combined Element Method for Shear Lag Effect Analysis of Thin-walled Box Girder. J. Highw. Transp. Res. Dev. 2014, 31, 58–63+75. [Google Scholar] [CrossRef]
  3. Zhang, Y.; Liu, Z.; Wang, C. Rational Method for Computing Polar Inertia Moment of Thin-walled Box Girders with Cantilever Slab. J. China Railw. Soc. 2019, 41, 94–99. [Google Scholar] [CrossRef]
  4. Chen, Y.; Dong, J.; Xu, T.; Xiao, Y.; Jiang, R.; Nie, X. The shear-lag effect of composite box girder bridges with corrugated steel webs and trusses. Eng. Struct. 2019, 181, 617–628. [Google Scholar] [CrossRef]
  5. Zhou, S. Finite beam element considering shear-lag effect in box girder. J. Eng. Mech. 2010, 136, 1115–1122. [Google Scholar] [CrossRef]
  6. Luo, Q.; Wu, Y.; Li, Q.; Tang, J.; Liu, G. A finite segment model for shear lag analysis. Eng. Struct. 2004, 26, 2113–2124. [Google Scholar] [CrossRef]
  7. Wei, C.; Li, B.; Zeng, Q. Transfer Matrix Method Considering Both Shear Lag and Shear Deformation Effects in Non-Uniform Continuous Box Girder. Eng. Mech. 2008, 25, 111–117. [Google Scholar]
  8. Ding, N.; Lin, L.; Qian, Y. An Approximate Method to Analyze the Effect of Shear Lag and Shear Deformation of Box Beam with Varying Depth. J. Railw. Sci. Eng. 2011, 8, 14–18. [Google Scholar]
  9. Chen, D.; Xu, Z.; Wang, Y. Derivation of Stiffness Matrix of Variable Section Beam Element Based on Energy Variation Method. In Proceedings of the 2019 Exchange Meeting on New Materials, Technologies and Their Engineering Applications in Civil Engineering, Beijing, China, 17 May 2019; Zhengzhou University: Zhengzhou, China, 2019; pp. 1639–1645. [Google Scholar]
  10. Zhang, Y.; Zhang, Y.; Zhang, H. Influence of Lateral Displacement of Flanges on Shear Lag Effect of Box Girder. J. Southeast Univ. (Nat. Sci. Ed.) 2022, 52, 247–253. [Google Scholar] [CrossRef]
  11. Lin, P.; Zhou, S.; Liu, F. Additional Axial Force Analysis Caused by Parabolical Warping Displacement about Shear Lag. Eng. Mech. 2010, 27, 90–93+109. Available online: https://pubs.cstam.org.cn/article/id/em_1213 (accessed on 1 August 2010).
  12. Song, A.; Li, Z.; Zhou, P.; Wang, S.; Su, Q. Analysis of Shear Lag Effect in Thin-walled Box Girder Based on Cosine Warping Displacement Function. Highw. Eng. 2021, 46, 12–18. [Google Scholar] [CrossRef]
  13. Guo, W.; Yang, Z.; Feng, Q.; Dai, C.; Yang, J.; Lei, X. A new method for band gap analysis of periodic structures using virtual spring model and energy functional variational principle. Mech. Syst. 2022, 168, 683–692. [Google Scholar] [CrossRef]
  14. An, H.; Wu, X.; Shi, Y.; Yin, Y. Theoretical analysis of shear lag warping displacement function of box girder. Eng. J. Wuhan Univ. 2020, 53, 692–697. [Google Scholar] [CrossRef]
  15. Lin, P.; Liu, F.; Yang, J.; Ji, W.; Liu, Y. Definition and Application of Shear-lag Warping Displacement Functions for Box Girders. Chin. J. Comput. Mech. 2012, 29, 789–794. Available online: https://pubs.cstam.org.cn/article/id/cjcm_jslxxb201205025 (accessed on 1 October 2010).
  16. Penalba, M.; Ulazia, A.; Ibarra-Berastegui, G.; Ringwood, J.; Sáenz, J. Wave energy resource variation off the west coast of Ireland and its impact on realistic wave energy converters’ power absorption. Appl. Energy 2018, 224, 205–219. [Google Scholar] [CrossRef]
  17. Lu, P.; Li, D.; Zhang, Y.; Zhou, Y. Bending characteristics of composite box girder considering all surfaces shear deformation based on Hamiltonian system. Mech. Based Des. Struct. Mach. 2023, 51, 2583–2606. [Google Scholar] [CrossRef]
  18. Zhang, Y.; Hu, Y.; Lin, L. Analysis on Shear Lag Effect of Thin-Walled Box Girders Based on A Modified Warping Displacement Mode. China Civ. Eng. J. 2015, 48, 44–50. [Google Scholar]
  19. Elshazli, M.T.; Abdulazeez, M.M.; ElGawady, M.; Ibrahim, A. Comprehensive Numerical Modeling of Prestressed Girder Bridges under Low-Velocity Impact. Buildings 2024, 14, 640. [Google Scholar] [CrossRef]
  20. Maleska, T.; Beben, D. Behaviour of soil–steel composite bridges under strong seismic excitation with various boundary conditions. Materials 2023, 16, 650. [Google Scholar] [CrossRef]
  21. Jing, Y.; Zhang, X.; Zhou, Y.; Zhao, Y.; Li, W. Dynamic response and impact force calculation of PC box girder bridge subjected to over-height vehicle collision. Buildings 2023, 13, 495. [Google Scholar] [CrossRef]
  22. Wang, Q.; Dai, R.; Wang, S.; Wang, Y. Prediction Algorithm and Interpretation for Autogenous Shrinkage of Concrete Based on Machine Learning. J. Shenyang Jianzhu Univ. 2023, 39, 1050–1057. [Google Scholar]
  23. Zhou, M.; Zhang, Y.; Lin, P.; Zhang, Z. A new practical method for the flexural analysis of thin-walled symmetric cross-section box girders considering shear effect. Thin-Walled Struct. 2022, 171, 108710. [Google Scholar] [CrossRef]
  24. Huang, J. Torsional Analysis of Thin-Walled Structures: Curved Beams and Diagonally Supported Box Girders; China Railway Publishing House: Beijing, China, 1988. [Google Scholar]
  25. Gupta, T.; Kumar, M. Flexural response of skew-curved concrete box-girder bridges. Eng. Struct. 2018, 163, 358–372. [Google Scholar] [CrossRef]
  26. Shan, Z. Mechanics of Materials, 3rd ed.; Higher Education Press: Beijing, China, 2016. [Google Scholar]
  27. Cai, J.; Moen, C. Elastic buckling analysis of thin-walled structural members with rectangular holes using generalized beam theory. Thin-Walled Struct. 2016, 107, 274–286. [Google Scholar] [CrossRef]
  28. Guo, J. Theory of Box Girder Design, 2nd ed.; People’s Transportation Publishing House: Beijing, China, 2008. [Google Scholar]
  29. Bao, S. Structural Mechanics of Thin-Walled Bars; China Architecture Publishing & Media: Beijing, China, 2016. [Google Scholar]
  30. Dym, C.; Shames, I. Utilixue Blanfenfa; China Railway Publishing House: Beijing, China, 1984. [Google Scholar]
  31. Lai, Z.; Jiang, L.; Zhou, W.; Liu, X.; Yu, J.; Zhang, Y. Analytical soluteon to mapping rail deformation under bridge transverse deformation using energy variational principle. J. Cent. South Univ. 2022, 29, 2654–2664. [Google Scholar] [CrossRef]
  32. Luo, M.; Lin, P. Experimental Study of Shearing Force Lag Effect of Three-cell Box Girder. J. Exp. Mech. 2015, 30, 621–628. [Google Scholar] [CrossRef]
  33. Zhang, K. The Theoretical and Experimental Investigations on the Space Mechanical Behavior of Single Box and Triple Cells Girder under Symmetrical Load. Master’s Thesis, Lanzhou Jiaotong University, Lanzhou, China, 2017; 67p. Available online: https://cdmd.cnki.com.cn/Article/CDMD-10732-1018234535.htm (accessed on 3 June 2017).
  34. Chu, F.; Cao, P.; Wan, L. Similarity Theory and Its Application to Simulation Tests. Min. Technol. 2007, 4, 64–65+78. [Google Scholar] [CrossRef]
Figure 1. Loading force diagram.
Figure 1. Loading force diagram.
Applsci 14 04778 g001
Figure 2. Coordinates and size diagram of a single-box three-chamber box girder section.
Figure 2. Coordinates and size diagram of a single-box three-chamber box girder section.
Applsci 14 04778 g002
Figure 3. Simplified schematic of the box girder sections.
Figure 3. Simplified schematic of the box girder sections.
Applsci 14 04778 g003
Figure 4. Diagram of the thin-walled element.
Figure 4. Diagram of the thin-walled element.
Applsci 14 04778 g004
Figure 5. Shear flow diagram of a single-box three-chamber girder.
Figure 5. Shear flow diagram of a single-box three-chamber girder.
Applsci 14 04778 g005
Figure 6. Schematic of the structural dimensions, measurement points, and loads of the model (unit: mm).
Figure 6. Schematic of the structural dimensions, measurement points, and loads of the model (unit: mm).
Applsci 14 04778 g006
Figure 7. Loading system used for the model test [32,33].
Figure 7. Loading system used for the model test [32,33].
Applsci 14 04778 g007
Figure 8. Cross-sectional dimensions and measurement-point arrangement of the finite-element model (unit: mm).
Figure 8. Cross-sectional dimensions and measurement-point arrangement of the finite-element model (unit: mm).
Applsci 14 04778 g008
Figure 9. Space finite-element model of plexiglass single-box three-chamber simple support box girder.
Figure 9. Space finite-element model of plexiglass single-box three-chamber simple support box girder.
Applsci 14 04778 g009
Figure 10. Comparison of the shear lag coefficients of the top wing plate of the section at y = 1450   cm for a single-box, three-chamber box girder.
Figure 10. Comparison of the shear lag coefficients of the top wing plate of the section at y = 1450   cm for a single-box, three-chamber box girder.
Applsci 14 04778 g010
Figure 11. Comparison of the flexural stresses in the upper-wing plate of the box girder at y = 1450   cm for a single-box three-chamber box girder at φ = 1.6 .
Figure 11. Comparison of the flexural stresses in the upper-wing plate of the box girder at y = 1450   cm for a single-box three-chamber box girder at φ = 1.6 .
Applsci 14 04778 g011
Table 1. Comparison of the shear lag coefficients at the key points of the top plate of the L / 2 and L / 4 span sections under centralized loading.
Table 1. Comparison of the shear lag coefficients at the key points of the top plate of the L / 2 and L / 4 span sections under centralized loading.
LocationKey PointTheoretical ValueFinite-Element ValueMeasured Value
L / 2 cross-section T 6 0.9720.8460.877
T 5 1.0011.0200.97
T 4 1.0301.6441.098
T 3 0.9731.0661.041
T 2 1.0371.6631.093
T 1 0.9741.0340.842
L / 4 cross-section T 1 0.9951.0251.323
T 2 1.0071.0231.550
T 3 0.9971.0041.005
T 4 1.0040.980.988
T 5 0.9990.980.168
T 6 0.9940.9770.931
Table 2. Comparison of the shear lag coefficients at the key points of the bottom plate of the L / 2 and L / 4 span sections under concentrated loading.
Table 2. Comparison of the shear lag coefficients at the key points of the bottom plate of the L / 2 and L / 4 span sections under concentrated loading.
LocationKey PointTheoretical ValueFinite-Element ValueMeasured Value
L / 2 cross-section B 4 1.1271.1531.093
B 3 0.9050.8920.972
B 2 1.0351.0200.993
B 1 0.8640.8710.884
L / 4 cross-section B 1 0.9690.9651.116
B 2 0.9800.9640.97
B 3 1.0031.0120.996
B 4 1.0331.0410.977
Table 3. Statistics of the shear lag coefficients of the wing plate at the y = 1450   cm section for the cantilever correction.
Table 3. Statistics of the shear lag coefficients of the wing plate at the y = 1450   cm section for the cantilever correction.
x Coordinate (mm)Finite-Element ValueCorrection Parameter φ
1.01.11.21.31.41.51.61.7
−80100.82750.97540.96140.94960.92620.90490.87640.83340.8014
−73400.91070.99290.98550.98010.96650.95510.93980.91650.8994
−66700.98891.01121.00931.00931.00391.00060.99630.98980.985
−60001.04321.03131.03251.0361.03541.03711.03941.04291.0455
−53401.03221.01051.01231.01631.01671.01931.02291.02821.0323
−46701.02080.99470.99681.00121.00231.00551.00991.01651.0215
−40001.01780.98050.98380.98950.99250.99761.00461.01501.0227
−33401.02160.99650.99871.00341.00471.00841.01321.02051.0260
−26701.03581.01391.01591.02031.021231.02441.02871.03521.0400
−20001.04621.02251.02471.02941.03071.03451.03941.04701.0524
−13401.03551.00731.00981.01451.01611.02011.02531.03311.0389
−6701.01380.9880.99070.99560.99751.00161.00711.01541.0216
01.00590.97530.97850.98370.98620.99070.99701.00631.0132
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

Li, Y.; Xie, J.; Wang, F.; Li, Y. Shear Lag Effect of Ultra-Wide Box Girder under Influence of Shear Deformation. Appl. Sci. 2024, 14, 4778. https://doi.org/10.3390/app14114778

AMA Style

Li Y, Xie J, Wang F, Li Y. Shear Lag Effect of Ultra-Wide Box Girder under Influence of Shear Deformation. Applied Sciences. 2024; 14(11):4778. https://doi.org/10.3390/app14114778

Chicago/Turabian Style

Li, Yanfeng, Jiyuan Xie, Fengchi Wang, and Yuanhui Li. 2024. "Shear Lag Effect of Ultra-Wide Box Girder under Influence of Shear Deformation" Applied Sciences 14, no. 11: 4778. https://doi.org/10.3390/app14114778

APA Style

Li, Y., Xie, J., Wang, F., & Li, Y. (2024). Shear Lag Effect of Ultra-Wide Box Girder under Influence of Shear Deformation. Applied Sciences, 14(11), 4778. https://doi.org/10.3390/app14114778

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop