Nonlinear Buckling Analysis of Thin-Walled Box Beams considering Shear Lag

In this work, a general geometric nonlinear model of straight thin-walled box beams (STBBs) under combined eccentric and axial loads is established. In order to accurately reflect the behavior of STBB, the additional shear lag warping is added to enrich the displacement field. It is necessary to define the section shape function to describe the local section deformation. (erefore, extension, bending, torsion, distortion, and shear lag effects are expressed by the generalized coordinate method. Based on the stability of transverse unconstrained box beam theory, meaningful higher-order solutions can be obtained by defining a set of coupled deformation modes. (e equilibrium equation is discretized by the Galerkin method, and the Newton–Raphson incremental method is used to derive and solve the nonlinear governing equations. On this basis, the analytical expression of stiffness matrix is established. For solving the stability problem, the effectiveness of the proposed method is verified by comparing the calculation results of shell element (Ansys) with other theories. Numerical examples even show that the proposed method can not only get the influence of shear lag but also obtain the variation of lateral buckling of the beam model.


Introduction
As a kind of important load-bearing structural member, performance and design are a very significant aspect of thinwalled elements because they minimize the weight and cost of the structure. e classical Vlasov beam theory, St. Venant torsion theory, and the thin-walled beam theory put forward several calculation formulas including transverse deformation and displacement, according to the research of different open and closed cross sections [1]. On the account of the kinematics hypothesis that the shear strain was ignored in the central line of the cross section, classical formulas about distortion theory were proposed on the basis of Vlasov theory. Correia et al. discussed the buckling behavior and strength of glass fiber reinforced polymer (GFRP) under uniform compression and conducted experimental, numerical, and analytical studies, focusing on local buckling [2]. Based on the theory of small strain and medium rotations, Ascione et al. established the kinematics model of composite thin-wall beam. e effects of geometrical imperfections and shear deformation on the members were fully considered, and the prebuckling behavior, the global buckling, and local buckling analyses were carried out [3][4][5][6].
In general, the deflection of box beam under self-weight and load may be quite different when the shear lag effect is ignored in the traditional beam theory. Reissner was the first one to use the variational energy method to analyze the biaxial symmetric rectangular box girder and the quadratic parabola as the displacement function of shear lag buckling [7]. Later, many scholars adopted quadratic, cubic, or even quadrilateral parabolas as the displacement function of shear lag on the basis of Reissner's theory [8][9][10][11][12]. Based on these methods, Luo and Wu selected three independent generalized displacements of top, bottom, and side cantilever plates, respectively, to represent different shear lag warping displacements between box girder plates with trapezoidal sections [8]. According to the energy method, Lin and Zhao described the nonuniform longitudinal displacement of flanges by using the higher-order polynomial of the infinite term and further quantified the shear lag effect of thin-walled bending members such as T-beams and box beams [13]. In order to solve the problem of plastic deformation of steel box girder in the limit condition, Lin and Zhao proposed a method for modeling the behavior of inelastic shear lag [14]. Lertsima et al. used shell element modeling for the entire simply supported box girder and studied the stress concentration problem at the flange under the influence of shear lag by using the three-dimensional finite element method [15]. In the same way, Sa-nguanmanasak treated the influence of finite element mesh on shear lag in detail by using the method of multi-degree-offreedom extrapolation and derived the empirical formula to calculate the stress concentration coefficient including the shear lag effect [16]. en, Zhang put forward generalized formulas, which included the geometric properties of sections and generalized internal forces concerning shear lag moments and bending moments [17]. Chen et al. obtained a closed solution of shear lag based on the variational method of energy and carried out harmonic analysis on the closed beam [18]. In addition, Li improved the PM-DP analysis method, which can be used to correct the defects of shear lag warping deformation of the top cantilever flange [19]. e good structural performance of thin-walled box girder has always been the focus of research by many scholars. Considering the coupling of bending, torsion, shear deformation, and warping shear stress, Wang proposed a new type of nonshearing locking thin-wall section space beam element [20]. In view of distortion, a new one-dimensional theory was proposed by Law and Kim, in which the static and dynamic studies on thin-walled closed beams with general cross sections were carried out [21]. e research of Xu and Ye showed that the single box beam had no obvious pure distortion, and the influence of distortion shear stress on shear deformation was analyzed [22]. Razaqpur and Li developed a new finite element model of curved thin-walled multicompartment box girder by combining the generalized Vlasov theory with the FEA technique, which could simulate the extension, bending, torsion, torsional warping, distortion, distortional warping, and shear lag effects of thin-walled box beam [23]. In order to calculate the section deformation, Choi and Jang adopted the beam frame method of two-dimensional Euler beam for simulation and developed the FEA model of quadrilateral thin-walled beam with variable section using the high-order beam theory [24]. It was approximated using the Ritz-Galerkin method's independent descriptions of stress and displacement fields. According to the previous work, Tan further analyzed nonlinear lateral buckling of thin-walled box beam's webs with the unequal thickness by considering the coupling of the bending, torsion, and distortion [25]. e deformation of box girder section under eccentric load and axial force is also studied. Genoese and Pietraszkiewiczestablished the linear kinematic equations for a thin-walled beam of heterogeneous isotropic material [26]. Opoka and Pietraszkiewicz showed the complete boundary value problem (BVP) for geometrical nonlinearity of thin shells, which was the stress action and bending of the reference surface of shells [27]. Coupling both material anisotropy and geometric nonlinearity, Vo presented a general geometric nonlinear model of arbitrary spatial composite thin-walled beams under different loading [28]. Miranda and Gutiérrez developed a new generalized beam theory (GBT), which modified the kinematics formulations of shearing deformation GBT put forward by Silvestre and Camotim in the early period [29]. Considering the effects of axial warping effect and material anisotropy on the surface shear strain in the shell, Zhang proposed an improved structural mechanics model of composite beams with single closed thin-walled section [30]. Meanwhile, Sahraei extended the shear deformation theory and FEA method with computational efficiency to the calculation of the transverse torsional buckling capacity of single symmetric I-beams under general load [31]. Moreover, Cambronero-Barrientos identified a beam element with three nodes, each of which has five degrees of freedom, which exceeded six degrees of freedom of traditional three-dimensional beam [32]. And the influence of homogeneous and nonhomogeneous torsion, distortion, and shear lag had been considered in the normal stress distribution at the same time. en, Meftah investigated the elastic lateral buckling distortional deformation of thinwalled box girder under the combined action of bending and axial loads [33].
In this paper, shear lag effect is considered as a separate mechanical property in the displacement deformation. So, the displacements of straight thin-walled box beams (STBBs) are divided into the sum of the partial displacements of six deformation fields. e main purpose of this work is to consider the influence of shear lag on nonlinear buckling performance of box girder under the combined action of an eccentric load and axial force.

Theoretical Model for Buckling Analysis
In this section, the main symbols and conventions of symbols are given, then the displacement field are described, and the kinematic equations are given on the basis of references [12,21]. erefore, the stiffness matrix is finally obtained.
It is assumed that the displacement field is the sum of the modes of each displacement deformation. Here, since the shear forces Q z and Q y are unknown, it is essential to modify it through the orthogonal warping shear function ω ζz and ω ζy when the shear lag effect is analyzed.

Basic Assumptions and Considerations.
e present assumptions about the structural model are as follows: (1) e material of the model is linear elastic and homogeneous. e influence of the thickness of each wall on warping is ignored. (2) e strain on the contour of each plane is zero, that is, the plane of STBB cannot expand in the tangential direction. (3) e member buckles under the nonextensional mode, so the curvature and central longitudinal strain of the xz plane remain zero during the instability process. is means that members are supposed to buckle under combined axial and bending loads.  Figure 1 shows the global Cartesian coordinate system (x, y, z), where the x-axis is the same as the longitudinal direction of STBB. (s, n, x) is the local coordinate system. For each wall, s passes through the whole cross section along the center line of each plate, and n is the coordinate of the local normal line. e generalized displacement theory of STBB is based on the assumption of deformation modes (in Figures 2 and 3) corresponding to the extension, bending, torsion, shear lag, and distortion. ese modes can be separated into a group of coupled axial and transverse displacement fields, each of which can be decomposed into the functions of the plane and axial coordinates. In the following, it makes sense to note that the displacement modes including global and local modes are described in detail.
Extensional mode, which only produces constant displacement u in the direction of length, is the simplest kind of displacement.
Two bending modes generate a displacement and a rotation, respectively. e rigid solid displacements (v and w) are in the plane of cross section and the rotations (θ y and θ z ) are about the y-and z-axes.
Shear lag deflections considering shear forces, the generalized displacement of shear lag deformations (ζ z and ζ y ), and shear lag warping displacement functions (ω ζz and ω ζy ) are introduced, which are proportional to shear forces Q z and Q y , respectively. erefore, the additional displacements caused by the shear lag are u ζz � − ω ζz ·ζ z and u ζy � − ω ζy ·ζ z . In the analysis of shear lag effect of thin-walled box beam, reference [12] interprets the quadratic parabola theory as a reasonable model of warping displacement function mode, the specific form of which will be detailed in the following part.
Torsion mode (no transverse deformation occurs in the cross section) generates homogeneous and nonhomogeneous torsional warping ω a and ω at , which are related to the normal stress σ x � Eω a ·d 2 θ x /dx 2 and the shear stress τ sx � G·zω at /zs·d 3 θ x /dx 3 , respectively. erefore, the displacement field contains a rigid solid rotation θ x in regard to the shear center (x 0 , y 0 ) of the cross section and two longitudinal displacements-the first-and third-order derivatives of torsion rotation are multiplied by the two warping functions, respectively. Nonhomogeneous torsional warping ω at is caused by the shear stresses. For the calculation of longitudinal displacement, the influence of nonhomogeneous torsion on normal stresses can be ignored.
Distortion mode is the transverse deformation of the cross section. e difference between homogeneous and nonhomogeneous distortion is not discussed in traditional theories. For distortion mode, transverse displacements are (η·χ, ψ s ·χ) in the plane of the cross section. And in terms of free warping or nonfree warping, homogeneous distortion and nonhomogeneous distortion can be obtained. Considering that the thickness of the beam girder compared with other dimensions is very small, the compression of the box wall along the direction of the wall thickness was neglected currently. en, the longitudinal displacement at any point of the beam can be given by the homogeneous and nonhomogeneous distortion warping ω d and ω dt multiplied, respectively, by the first-and third-order derivatives of distortion rotation χ.
In the analysis of deformation modes, the transverse and longitudinal displacements are introduced as (1c)

Structural Mechanisms.
In order to illustrate the general method of obtaining deformation modes of STBB under combined loads, the geometric dimension of STBB is defined as follows: b is the width; h is the height; L denotes the entire beam length; and t i denotes the thickness of each wall which is much smaller than other dimensions. It is well known that the shear lag effect inevitably increases the deflection of the box girder. e shear lag deformation state is described by generalized displacement, that is, additional deflection caused by shear lag effect. e shear lag warping displacement function is not a continuous function throughout the beam that directly affects the size and distribution of cross section stress. erefore, many authors studied and proved these interpolation of functions [8][9][10][11][12]. e shear lag warping stress σ ωζz � − Eω ζz ω ζz ′ (x) (without the influence of axial force) of the cross section of the thin-walled structure should meet the self-balancing condition, that is, When the thin-walled box beams are bent, the normal stress of the cross section must satisfy the axial equilibrium condition, which means A ω ζz dA � 0. After many comparisons, it is reasonable to select the quadratic parabola as the shear lag warping displacement function. In this paper, the quadratic parabola expression of the selected shear lag warping displacement function is

Mathematical Problems in Engineering
where d is the additional uniform displacement function modified by the shear lag warping function and A is the cross section area; this method is valid and correct for this paper. e microelement of generalized distortion displacement in STBB wall is shown in Figure 4. In accordance with any point N 0 in the middle of plane, the longitudinal, normal, and circumferential displacements of the plane at any point N 1 can be obtained by ignoring the compression of the wall in the direction of the plane thickness. In the light of middle shear strain c xs � (zψ s z χ /zx) + (− zωdχ ′ /zs) � 0, the   distortion warping can be expressed as ω d (s) � 0 s ψ s (s)ds. From the self-equilibrium condition of the warping stress, In the case that the box girder has no torsion but only distortion deformation, the interposition pure torsion force system 1/Ω is added at both ends of the microsection dx of STBB and the circumferential displacement ψ s δχ of the distorted microsection is taken as the virtual displacement. en, the virtual work done by the unit pure torque force on the virtual displacement is zero: [(1/Ω)ψ s (s)δ χ ]ds � 0, that is, ψ s (s)ds � 0. rough the frame angular point displacement and the continuity of bending moment, the coefficients can be obtained in reference [21]: where L side is the length of each wall of the cross section and α is a coefficient to satisfy the static force equilibrium. e displacement of STBB is given by

Mathematical Problems in Engineering
in which i � 1∼4 are the wall numbers. e functions Z i (s), Y i (s), h i (s), ω i α (s), and η(s, n) in above equations are expressed as follows:

Expressions of Strains and Stresses.
In order to explain the strain tensor of STBB, the normal and shear strains containing the linear and nonlinear strains are expressed as follows: In the case of normal strain, the expression according to the linear and nonlinear conditions can be denoted by εi · l and εi · nl: e shear strain c i xs can be represented as   Mathematical Problems in Engineering e stresses are carried out by multiplying the normal and shear strain by E (Young's modulus constant) and G (shear modulus constant), respectively.
2.5. Variational Formulation. According to the energy variational principle, the total system potential energy is composed of external potential and internal strain energy (W and U). erefore, the first-order expression of the total potential energy of the system is Substitute longitudinal and shear strain equations (10b)-(10d) into the elastic energy equation (13) to get the internal strain energy expression: where N represents axial force, B ω represents the bimoment, M ωaωd is the high-order coupling bimoment; Q ω ζ z ′ is the generalized shear force; M 1 to M 8 are high-order torques caused by torsion and distortion; M ω a ω ζ z ′ is the shear lag coupling torque; and M ω d ω ζ z ′ is the shear lag distortion coupling torque.
and detailed expressions of these forces and torques can be found in Appendix A. e general formulas about of STBB nonlinear and the buckling analysis are obtained when it is subjected to external forces of combined axial load F x and eccentric load F z . In the initial transverse buckling, the external energy is given: in which W r (s, r) stands for the vertical displacement. In consideration of the classical quadratic kinematics hypothesis, the change of external energy can be expressed as where e y and e z are eccentricities of the concentrated load F z relative to the shear center. According to equations (14) and (17), the virtual displacements and their derivatives can be used to express the change of the total potential energy. Due to the generalized variational principle, balance equations with reference to a first-order variation zu(x), zv(x), zw(x), zθ(x), zχ(x), and zζ(x) can be written as 8 Mathematical Problems in Engineering

e Discrete Equilibrium Problem.
For nonlinear analysis, the governing equations are simplified to generalized coordinates by the Ritz variational method. Based on the Newton-Raphson method, an incremental-iterative method is proposed to solve the nonlinear equilibrium equations. And the buckling load is first determined according to the singularity condition of the tangential stiffness matrix determinant. e discrete kinematic equation can be used to analyze the performance of simply supported structures. e displacement modes adopt the approximate functions that are consistent with the boundary conditions of STBB: where u 0 , v 0 , w 0 , θ 0 , χ 0 , and ζ z0 are related to displacement amplitudes. In addition, these functions can also represent exact solutions of nonlinear equilibrium equations. rough Galerkin's approximation method and integration, the algebraic coupling equilibrium equations of displacement modulus can be obtained: in which M 0 � F z L/4 represents the maximal bending moment when the eccentric load is applied at the middle length of the beam. e detailed expressions of P i algebraic system is shown in Appendix B.

Tangential Stiffness Matrix for Lateral
Buckling. e beam deflects under the action of axial and bending loads. If the beam is loaded to a certain extent, its cross section rotates around the axis of the beam, at which time lateral buckling may occur. is is transverse buckling, and the buckling load is the critical load. Corresponding to the prebuckling state, the deflection is called the "fundamental state." On account of the Jacobian Matlab function, the tangent stiffness matrix [K t ] is automatically defined by the equilibrium equations (20a)-(20f ). In the case of nonlinear highcoupling high-order model and considering the prebuckling state, it can be reasonably assumed that each displacement (u 0 , v 0 , w 0 , θ 0 , χ 0 , ζ z0 ) is approximately equal to (0, 0, w 0 , 0, 0, 0) by analyzing the geometric matrix. Here, the tangent stiffness matrix of the orthogonal displacement modes is Finally, after eliminating the determinant of the tangent stiffness matrix, we get polynomial equations containing F x and F z . In this paper, the roots of the equation yield to transverse buckling loads in different buckling modes (extension, bending, torsion, distortion, and shear lag modes).

Application and Numerical Examples
e goal of the section is to verify the current theoretical model of STBB buckling behavior so as to demonstrate the accuracy and practical application of the formulas. At the same time, the results are compared with those calculated by FEA (Ansys software), Meftah, and Tan. e analysis procedure for the nonlinear buckling analysis of STBB is the same as the general nonlinear analysis process. In the Ansys software, the STBB is idealized by the 4-node 181 isoparametric quadrangle shell elements ( Figure 5). First, a "perfect" model is established, and the sparse matrix solver is used to conduct eigenvalue buckling analysis, which can predict the critical buckling force and provide a basis for the nonlinear buckling analysis. e grid can generate 3888 nodes and 3840 elements. Meanwhile, the grid is verified to discretize these patterns accurately. In the following analysis, perturbation is introduced in the geometry by adding geometric imperfections using the UPGEOM command. e geometric nonlinear analysis of structures with geometric perturbation is conducted using the arclength method, where the postbuckling behavior can be obtained. e size of the perturbations is usually a small percentage of a relative structure, and the relevant scaling factor used in this case is 0.001. Finally, the bifurcation point of the buckling load displacement curve is found out, which is the critical buckling load.
In Ref. [33], Meftah studied the lateral buckling of thinwalled box girder under the combined action of bending and axial loads. And in the stress calculation, the nonlinear higher-order terms were omitted. Similarly, in Ref. [25], Tan further analyzed nonlinear lateral buckling of STBB's webs with the unequal thickness, but the main difference was that the geometrically nonlinear high-order coupled terms (wθ 2 , wχ 2 , w 2 θ. . .) had been taken into account. Both of them did not consider the shear deformation and the shear lag effect. Figure 6 shows a simply supported beam under combined eccentric and axial load. e critical load F z imposed in the middle of beam length takes action on the intersection of the upper flange and the center line of the main web, and the axial load F x is applied on one end of the upper flange. e geometrical and material properties are F x � 2.5 kN, E � 210 GPa, t 1 � t 3 � 4 mm, t 2 � 8 mm, and t 4 � 5 mm. Poisson's ratio is 0.3. e length ranges from 1.0 to 8.0 m. In addition, in order to illustrate the influence of the lengthwidth ratio of cross section on nonlinear lateral buckling of beams, five different groups of cross sections are selected in this paper. Table 1 shows the detailed cross section dimension of the models.  illustrate the comparison between the present and other three theories under different lengths and cross section ratios. e parameters of the beam can indicate the sufficient effectiveness and accuracy of the method in this paper. Tables 2-4 report percentage differences of the displacements (Δu i , Δw i , Δv i ) obtained by theoretical methods and those obtained by FEA. Where, i � 1∼3 represents the present, Meftah, and Tan's methods, respectively.

Example 1: Variation of Length and Cross Section.
As can be seen from Figure 7, the numerical values of the longitudinal displacement are almost between the values of Meftah and FEA. It is easy to see from the picture that Meftah's axial displacements are independent of the change of cross section. At the same length, the longitudinal displacements of the other three methods increase gradually with the decrease of cross section ratio (b/h). It should be noted that the relative errors of present method are the smallest while those of Meftah's method are the largest in Tables 2-4. It is shown that ignoring the shear lag effect will lead to an overestimation of percentage of the longitudinal displacements at any length. With reference to Tan's beams, the values Δu 1 and Δu 3 are the closest, when the aspect ratio of cross section is 1. is is due to the geometrical nonlinear effect caused by bending, torsion, ,distortion and shear lag coupling, which leads to a decrease in the axial stiffness of STBB. It is confirmed that the axial displacements of the present theory are sufficiently close to those predicted by FEA. Especially in IV and V, the values of present method are closest to those of FEA when the length is under 7 and 8 m (the smallest error Δu 1 � 6.6%). is shows that adding additional shear lag displacement is effective for the prediction of longitudinal displacement.
It can be seen from Figure 8 that the vertical displacement deformation trend is completely corresponding to the rigid body performance, which is true for all models. e results of the present theory are confirmed to be closer to those of FEA, which supports the validity of this theory. Comparing the present theory with the other two, the results are satisfactory and can be considered as a supplement to the shear lag. In the same case, the values of FEA are always greater than those of other methods, with the lowest values of Meftah. Because the three theoretical approaches do not take the vertical compressibility of webs into account under vertical load and more nonlinear factors are considered, the vertical displacement decreases. According to the relative errors of vertical displacements from Table 3, when the beam length is less than 2 m, the shear lag effect has little influence on the vertical displacement of STBB. When the beam tends to be slender, the percentage of vertical displacements is smaller (Δw 1 are almost all less than 10%). e present method can well capture the effect of shear deformation on the beam, so it has high prediction accuracy.
In Figure 9, we note that the transversal displacement distribution curve obtained by the present theory is in good agreement with the results obtained by the FEA, which indicates that the proposed beam theory can accurately   predict the structural behavior of STBB under different cross sections. It shows that the postbuckling equilibrium paths are stable. In I (b/h � 2/1), the change of transversal displacement is tiny and the curve is relatively stable when the length is more than 3 m. Again, the same thing happens in III when the length is greater than 5 m. In the other three cases, the curve can be considered as increasing almost linearly. It can be seen that the five displacement modes without considering the shear lag effect are consistent with the prediction of the transverse displacement by present theory. In Table 4, the percentages of transverse displacement of present method are closer than those of the other methods, indicating that the shear lag has a greater impact on the large ratio of the cross section. e results show that although six modes need to deal with more complex deformation relations to obtain the characteristics of transverse, the beam under different lengths and cross sections can be calculated arbitrarily and accurately. Furthermore, as expected, Tables 2-4 point out that the relative displacement errors of Meftah are the largest. e shear lag displacement has a significant influence on the displacements of the long beam (L > 3 m), but a slight  influence on the short beam (L ≤ 3 m). e results highlight that the influence of shearing deformation on displacement is great.
In Figure 10, the present results confirm that the proposed methodology agrees well with that of the shell analysis. When the length-width ratio is greater than 1, Meftah's torsional rotations are always bigger than other values, but when the length-width ratio is less than 1, those of FEA become the maximum. Because of ignoring the effect of shear deformation according to Meftah's theory, torsional stiffness is underestimated. When the length is more than 3 m, the torsional rotations start to become more sensitive. Interestingly, when the length between 4 m and 6 m, the torsional rotations grow slowly and the curve flattens out.
e calculated values of Tan are always slightly larger than those of Meftah. And the variation of the two is more moderate with the length. It is easy to see from Figure 10 that the present torsional rotations are smaller than those of other methods, especially when the length-width ratio is less than 1. e torsion rotations of FEA are slightly different from those of the present model, even in the distribution pattern. e reason is that the shear deformation effect is considered in present model (6 modes), and the shear stress terms not only include the generalized shear lag deformation, torsion, and distortion shear stress terms but also include the coupling torque of torsion and distortion. e variation trends of distortion rotation of the three methods are in good agreement with each other in Figure 11. e present distortion curve is similar to that of Tan. e results of this theory (6 models) are satisfactory, especially when b/h ≤ 1, and the curve value is always slightly smaller than that of Tan. In the case of the same length, the distortional rotations of IV case (b/h � 3/4) are smaller than those of other cross section states. In fact, shear effect reduces the distortion value of the short beam relative to the nonshearable theory. e distortion analysis with shear lag effect shows a conservative behavior compared with a more comprehensive distortion analysis of nonlinear terms. Figure 12 shows the shear lag warping displacement. It can be seen that the curve trend under the different cross sections is in good agreement within this length range. In general, the shear lag warping displacements gradually increase with the increase of length. For 2 ≤ L ≤ 6, the relative rate of change is very small (ζ z min � 3.9287E − 4m, ζ z max � 7.8926E − 4m). In the same case, the larger the length-width ratio is, the larger the shear lag warping displacement is. We can observe that the curve almost tends to a flat line in I (b/h � 1/2). In the case of long beams, the shear lag effect has a significant effect on the buckling behavior. From Figures 7-12, it can be observed that the influence of shear deformation on buckling behavior is more significant when the length increases from 4 m. Due to the influence of defects, the mechanical properties of STBB by the present method are obviously lower than those of FEA. When the cross section length-width ratio is 3/4, the geometric nonlinearity of STBB is optimal. erefore, the theory presented in this paper is valid for the thin-wall aspect ratio assumed in the next example.

Example 2: Variation of Height and ickness.
In this numerical example, the effects of beam length, web thickness, and load height on lateral resistance of STBB are studied. Figure 13 shows three different load positions.
Here, the material and geometrical properties are the same as the previous example. In example 1, it can be seen that when the aspect ratio of cross section is b/h � 3/4, the  Tables 5-7. It can be seen from the three tables that all of the models have good consistency. e critical moments of the FEA model are almost bigger than those obtained by linear and nonlinear analysis. is is because the FEA model (generally speaking, a membrane shell with restricted movement) is stiffer than theoretical model. In order to verify the correlation of the buckling phenomenon, four rigid transverse webs with different thickness are carried out by the FEA method. e results of FEA analysis confirm the correlation of grid deformation.
From the three tables, it is obvious that the critical moments decrease with the increase of beam length. And at the same length, the critical moments increase with the increase of web thickness. e results show that the influence of length on the critical moments is always greater than that of web thickness. is fact is clear and useful to all theories. Furthermore, it can be obtained that the Meftah, Tan, and present nonlinear numerical calculations of the critical buckling moment are in good agreement with the FEA results. Tan's results and present linear results are more conservative than the nonlinear buckling moment calculated by FEA. It should be noted that the nonlinear solution of Meftah is the smallest, followed by the nonlinear solution of Tan's method. at is to say, the more the nonlinear coupling terms are considered, the higher the calculation accuracy is. It should be noted that Meftah and Tan's solutions Table 3: e percentage differences of the vertical displacement (Δw i ).     under the same conditions are lower than those of the present approach; this is due to the generalized eigenvectors of the classical models and the initial assumptions-essentially without taking the shear deformation into account, which is very important for the different thicknesses of the web model. e nonlinear critical moments obtained by present method is closest to those obtained by FEA. And the linear solution of present method is larger than the linear solution of Tan's method; it can be seen that considering the effect of shear lag effect can more comprehensively supplement the mechanical properties of STBB. e difference between linear and nonlinear relative error is very impressive. e linear and nonlinear relative errors of Tan and the present method can reach about 8% and 8.63% at the maximum, respectively. From the relative errors Δ 1 , Δ 2 , and Δ 4 , it can be seen that the numerical and analytical solutions are consistent and highly nonlinear. e relative error Δ 2 is almost the smallest value, while Δ 5 is the biggest one. In this case, by comparing the moments of STBB before and after buckling, it can be seen that only the nonlinear solution is more accurate.
By comparing the three tables, it can be seen that when the load is applied to the top plane of STBB, the bending moment is the smallest. And when the load is applied to the bottom plane, the beam resistance to lateral buckling is the greatest. erefore, we can think that the lateral buckling strength depends on the load height parameter. In the postbuckling state, the shear lag effect term is added to maximize the resistance of STBB to lateral buckling. In general, as the web becomes thicker, the relative error becomes larger overall. Regardless of beam length and the location where the concentrated load was applied, M cr values for nonuniform web thickness (t 1 � 5 mm and t 2 � 8 mm) are higher than those for the same web thickness (t 1 � t 2 � 6.5 mm). With the increase of the length, the difference of the critical moments (for both uniform and unequal webs) becomes more obvious, and the unequal web of beams show better mechanical properties. Buckling moment is mainly affected by the web resistance. It is clear that the phenomenon is due solely to coupling-induced geometrical nonlinear effect.

Conclusion
A new kinematics model is used to analyze the stability of straight thin-walled box beam (STBB). Based on the conventional motion variables, the shear lag modes are added to obtain that the beam section contains six modes of extension, bending, torsion, distortion, and shear lag under the action of shear deformation. e analytical formulas take the deformation, load height, and unequal web thickness into account under combined axial and vertical loads. In this theory, six sets of coupled equilibrium differential equations are derived to describe the longitudinal and transversal deformation through six displacement fields. In particular, when the straight thin-walled box beam is subjected to combined axial and vertical loads, these formulas are naturally applicable to capture bending, torsional, distortion, and shear lag responses. Present theory puts forward a theoretical basis for studying the mechanical properties of STBB and the ways to be affected by adding shear lag effect. Unlike traditional Meftah and Tan's theories, present theory shows that the beam response is largely dependent on the shear lag. Various numerical results are compared with those of Meftah, Tan, and shell finite element (Ansys).
With regard to the different cross sections and lengths, the following conclusions are given: (i) When the cross section ratio is greater than 1, the geometric variation of cross section has little influence on the axial length. (ii) e shear lag displacement is conducive to the prediction of longitudinal displacement. (iii) For the long beam, the shear lag warping displacement can reflect the vertical displacement more accurately. (iv) Generalized shear lag deformation increases the predicted torsional stiffness of the beam.
(v) e distortion under the coupling effect of the six modes (adding shear lag effect) is more conservative than that of the five modes. (vi) e shear lag warping displacement increases with the increase of cross section ratio (b/h).
For the different load positions and the thicknesses of webs, the following results are reached: (i) e nonlinear critical moments of present method are closest to those of FEA. (ii) e relative error Δ 2 is almost the smallest value, while Δ 5 is the biggest one. (iii) e lateral buckling strength depends on the load height parameter. (iv) M cr values for nonuniform web thickness (t 1 � 5 mm and t 2 � 8 mm) are higher than those for the same web thickness (t 1 � t 2 � 6.5 mm).
e effect of shear lag can more comprehensively supplement the mechanical properties of STBB, and the nonlinear solution is more accurate. In the postbuckling state, the shear lag effect term is added to maximize the resistance of STBB to lateral buckling. e calculation results are in good agreement with the results of finite element and other theoretical analysis. e above observation results show that the geometry has a great influence on the nonlinear buckling of STBB. For these reasons, further development of this research should focus on the relationship between the material properties of beams and the modes. e developed analytical method is sufficiently efficient and accurate.

Data Availability
All data and models generated or used during the study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.