Published: 30 June 2017

Free vibration analysis of thin-walled curved box girder considering shear lag deformation

Heng Cai1
Hailin Lu2
Hongyin Yang3
1, 2, 3School of Resource and Civil Engineering, Wuhan Institute of Technology, Wuhan, China
Corresponding Author:
Hailin Lu
Views 58
Reads 19
Downloads 1467


In order to investigate dynamic shear lag effect, based on energy variation method and Hamilton principle, vibration governed differential equations of curved box girder are deduced by assuming different shear lag warping displacement modes, the explicit solution to bending frequency of curved box girder is worked out with Galerkin method. The results of a numerical example show that shear lag warping displacement functions have a limit influence on vibration frequency; second degree parabola or catenary are the appropriate shear lag warping displacement function of box girder; Compared the theoretical calculating results with numerical results of ANSYS, the error between them is only 0.57 %, they agree very well with each other, which demonstrates the correctness and reliability of theoretical deduction.

1. Introduction

Thin-walled box girder with wide flange is widely applied in the engineering construction for its merits such as great stiffness, light weight and so on. While shear lag effect [1] is a prominent problem in box girder, it means that normal stress of flange is non-uniform distributed along width direction because of flange shear deformation. This phenomenon will affect the safety and durability of box girder bridges [2]. Energy variation method, based on principle of minimum potential energy, can not only acquire closed-form solution [3, 4] but also get the real peak stress at the ribs, therefore it is widely applied in calculating shear lag.

Bridges will bear vehicle load, wind load, and sometimes seismic action when they are in service, under the impact of vehicle load and vertical seismic excitation, box girder considered curvature will create strong coupling effect of bending, torsion and shear lag effect, which will make it difficult to calculate its internal force and deformation under dynamic loads, while vibration frequency is an important parameter to reflect the dynamic characteristics of bridges especially when General code for design of highway Bridges and culverts [5] is used to calculate impact factor. Shear lag warping displacement mode is the premise and basis of energy variation method, because it can influence solution accuracy. In this paper, shear lag warping displacement modes are respectively assumed as parabola series [6-8], cosine function, exponential function and catenary, governed differential equations of curved box girder are deduced by energy variation method and Hamilton principle and the explicit solution is also acquired, meanwhile, shear lag warping displacement modes with high precision are studied.

2. System potential energy

In order to simplify mechanic model, assumptions are as following:

(1) The damp is not considered and the material is linear elastic.

(2) The flanges don’t conform to plane cross-section assumption while the webs do.

(3) Shear lag effect is only considered in vertical bending deformation.

(4) Transverse deformation and shear deformation in out-plane are neglected, they mean εy=εz=0, γyz=γxz=0.

The cross section of curved box girder is shown as Fig. 1.

Fig. 1Cross section type of curved box girder

Cross section type of curved box girder

According to the assumptions above, the longitudinal shear lag warping displacement is shown as Eq. (1):


where u(x,t) is maximum difference shear angle of flange and ζ(y) is warping function, they are assumed to be as following:

(1) Parabola series:

ζ(y)=1-ybi n, (n=2, 3, 4, 5).

(2) Cosine function:


(3) Exponential function:


(4) Catenary:

ζy=2-coshηy, η=ln2+3bi,

where bi (i=1, 2) are the half flange width and cantilever slab width respectively.

The total potential energy of curved box girder is shown as Eq. (6):


where Uw, Uf are web and flange energy, and Ut, Ur are free torsion and restricted torsion potential energy respectively, w, ϕ are vertical deflection and torsion angle, r is curvature radius, hi (i=1, 2) are the distances from cross section center to the upper and bottom flange. Besides, E, G are Young’s modulus and shear modulus, It, Iw are torsional moment inertial and warping moment inertial.

The vibration energy of curved box girder can be given as Eq. (7):


where l, m¯ are length and linear mass density.

3. Vibration differential equations

The generalized potential energy function is shown as Eq. (8):


Under the condition of small deformation, there is no external force, based on Hamilton principle t1t2δ(U-T)dt=0, the first variation of Eq. (8) can be got, it is shown as Eq. (9):

Fu-ddxFu'=0, d2dx2Fw''-ddxFw'-ddtFw˙=0, d2dx2Fϕ''-ddxFϕ'+Fϕ=0.

Based on Eq. (9), the vibration differential equations can be shown as Eq. (10), and the corresponding boundary conditions are shown as Eq. (11):

Fu'δu 0l=0,
Fϕ'-ddxFϕ''δϕ 0l=0, Fϕ''δϕ' 0l=0,
Fw'-ddxFw''δw 0l=0, Fw''δw' 0l=0, Fw.δw 0l=0,

where Is, Is¯ are inertia moment and generalized inertia moment of flange, Is¯=(I1/b12+I2/b22)b2, the number in bracket and the “pie” mean the derivative of x, dot is the derivative of time t, α, β, η are shear lag warping constants.

4. Bending vibration frequency

According to Eq. (11), the displacement field functions can be given as Eq. (12) [1].


Take derivative of Eq. (12) and substitute them into Eq. (10), the results are shows as Eq. (13):

EIwrπ4l4+EI+GItrπ2l2w0+EIwπ4l4+GItπ2l2+EIr2ϕ0+αEI s rπlu0=0,

In the Eq. (13), there exist nonzero solution to a homogeneous linear equation of w0, ϕ0, u0, therefore the coefficient determinant is zero, it is shown as Eq. (14):



K11=EI+EIwr2π4l4+GItr2-m¯ω2, K12=EIwrπ4l4+EI+GItrπ2l2, K13=αEIsπ3l3,
K22=EIwπ4l4+GItπ2l2+EIr2, K23=αEIsrπl, K33=βEIsπ2l2+ηGIs¯b2, K21=K12,
K31=K13, K32=K23.

The bending frequency of curved box girder is acquired based on Eq. (14).

5. Numerical example

The material properties and structural dimensions of curved box girder are shown in the literature [9], in this paper, dimensions are shown as Fig. 1, mass density ρ= 1180 kg/m3, structural parameters: t1=t2=tw=8 mm, b1=b2=h= 100 mm, curvature radius r= 2.5 m, curvature angleθ= 30°. The material parameters are as following: Young’s modulus E= 3000 MPa, Poisson’s ratio μ= 0.385, the results are shown as Table 1.

In order to verify the correctness of the theoretical formula, SHELL 63 element is adopted to found the spatial shell element model of curved box girder in ANSYS. Considering boundary conditions, Ux, Uy, Uz and Uy, Uz are restricted at the junction between web and bottom flange, while Ux, Uz, and Uz are restricted in the other side, vibration model and frequency are shown as Fig. 2.

It can be seen from Table 1 and Fig. 2 that vibration frequency under different shear lag warping displacement modes are almost the same, based on principle of minimum potential energy, the corresponding vibration frequency is inevitably minimum when the energy is also minimum, therefore second degree parabola and catenary are the appropriate function of the shear lag warping displacement mode, which can better reflect the spatial mechanical behavior of curved box girder.

Besides, the calculating results in this paper agree well with the ANSYS numerical analysis, the error between them is only 0.57 %. Compared with the bridge code [5], the results in this paper have higher precision.

Table 1Bending vibration frequency (Hz)

Displacement modes
Parabola series
Cosine function
Exponential function
Chinese bridge code [5]

Fig. 2Analysis results of ANSYS

Analysis results of ANSYS

6. Conclusions

The following conclusions are acquired:

(1) The theoretical results in this paper agree well with finite element analysis, which can verify the correctness and reliability of theoretical formulas.

(2) Shear lag warping displacement modes have a limit influence on vibration frequency of curved box girder.

(3) Second degree parabola and catenary can better reflect the spatial mechanical behavior of curved box girder, it is advised to preferentially adopt second degree parabola or catenary when shear lag effect is analyzed by energy variation method.


  • Xiang Haifan, Fan Lichu Advanced Theory of Bridge Structures. Second Edition, China Communication Press, Beijing, 2008.
  • Cai H., Lu H. L., Tang Z. 3D finite element analysis on shear lag effect of curved box girder under multi-dimensional seismic excitation. Vibroengineering Procedia, Vol. 10, 2016, p. 293-297.
  • Luo Qizhi Theory and Model Test Studies of the Shear Lag in Thin Walled Box Girders Based on Energy Principle. Hunan University, Hunan, 2005.
  • Zhou W. B., Jiang L. Z., Liu Z. J., et al. Closed-form solution to thin-walled box girders considering effects of shear deformation and shear lag. Journal of Central South University (English Edition), Vol. 19, Issue 9, 2012, p. 2650-2655.
  • JTGD60. General Code for Design of Highway Bridges and Culverts. China Communication Press, Beijing, 2015.
  • Zhang Yuanhai, Lin Lixia Initial parameter method for analyzing shear lag effect of thin-walled box girders. Engineering Mechanics, Vol. 30, Issue 8, 2013, p. 205-211.
  • Lin Pengzhen, Liu Fengkui, Ji Wei Analysis on shear lag effect of concrete box beam by variational principle. Journal of the China Railway Society, Vol. 35, Issue 2, 2013, p. 93-98.
  • Zhang Yuhong The Experimental Study and Theoretical Analysis in Shear Lag Effect on Thin-Walled Box-Girder. Lanzhou Jiaotong University, Lanzhou, 2011.
  • Lu Hailin, Cai Heng, Tang Zheng Analysis on shear lag effect of thin-walled curved box girder under moving load. Railway Engineering, Vol. 9, 2016, p. 14-18.

About this article

07 May 2017
08 May 2017
30 June 2017
Mechanical vibrations and applications
thin-walled curved box girder
shear lag warping displacement modes
vibration frequency
energy variation
finite element