Next Article in Journal
GPS-Induced Disparity Correction for Accurate Object Placement in Augmented Reality
Previous Article in Journal
Correction: Goo, S.; Yang, H.-D. Restoring Raindrops Using Attentive Generative Adversarial Networks. Appl. Sci. 2021, 11, 7034
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Use of the Effective Width for Simply Supported Generally Loaded Box-Girders

Department of Civil, Construction-Architectural and Environmental Engineering, University of L’Aquila, Piazzale Pontieri, Loc. Monteluco, 67100 L’Aquila, Italy
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(7), 2848; https://doi.org/10.3390/app14072848
Submission received: 5 March 2024 / Revised: 23 March 2024 / Accepted: 25 March 2024 / Published: 28 March 2024
(This article belongs to the Section Civil Engineering)

Abstract

:
A new method for evaluating the effective flange width of simply supported rectangular box-girders in bending is illustrated. It aims to overcome the lack in the literature of a general method able to supply the effective flange width for a general distribution of loads along the beam axis, not restricted, i.e., to the few cases reported in standard codes. The method consists of finding an analytical expression for the effective width under sinusoidal loads of arbitrary wavelength and combining the results in the framework of a Fourier analysis. The model allows investigating the shear-lag effects on displacements and stresses. The analytical formulation proposed is validated against numerical results supplied by refined Finite Element analyses.

1. Introduction

It is well-known that shear effects play an important role on the mechanical response of frames behaving as beam-like structures [1,2,3]. However, they also affect single thin-walled members, especially when the flanges are wide compared to the length of the member. Among these latter structures, the box-girders are particularly prone to the shear-lag phenomenon, which is caused by shear strains in the flanges. This topic is of some technical importance, since box-girders are widespread in civil engineering applications, such as highway bridges, railway bridges and footbridges. Indeed, in box-girders in bending, it happens that longitudinal membrane stresses, developed at the intersection between the flange and web, cause warping of the cross-section and alteration of the normal stresses distribution, with respect to those predicted by the classic Euler–Bernoulli and Timoshenko theories. In particular, the latter assumes that the cross-sections remain planar and predict normal stress constants on the flanges and linear on the webs, while in practice, they are non-uniformly distributed along the flange, with a maximum value at the web intersection, and this value is typically larger than those predicted by the classic Euler–Bernoulli and Timoshenko theories.
The shear-lag effect in box-girders has been largely studied both in the past [4,5] and in recent years [6,7], with the aim of evaluating the maximum value of the normal stresses in the cross-sections. The first person to study the shear-lag effect of rectangular box-girders was Reissner [8] in 1946, by using a variational method based on the principle of minimum potential energy. There is a vast amount of literature on this subject. The shear-lag effect has been studied by analytical methods [9,10,11], numerical methods [12,13,14] and experimental tests [11,15,16]. In particular, Evans and Taherian [17] opted for empirical equations to obtain normal stresses with the finite-stringer method, Wang and Rammerstorfer [18] adopted the finite strip method, Kristek et al. resumed the biharmonic analysis of isolated flanged plates previously applied by [19], and Qin and Liu [20] studied box-girder bridges under uniform load by the Symplectic Elasticity method. Moreover, several researchers adopted the finite-segment method [21,22], the finite element method [23,24,25] and the energy-based variational method [4,26,27].
A simple but efficient method to consider the shear-lag effect of box-girders in bending is the effective flange width method. An effective flange width is defined as the length of the flange in which the normal stresses are assumed constant and equal to the maximum stress arising at the intersection with the web, under the assumption that the cross-section remains planar. This is a simplified approach, but it usually involves reasonable errors in design. The use of the effective width method was originally proposed in 1924 by Von Karman [28] for aeronautical applications. Later on, it was adopted by Sedlacek and Bild [29], who gave a separate determination of the effective widths for bending and torsion, Wang and Rammerstorfer [18], who also investigated the effective width of postcritically loaded stiffened plates, and Wu et al. [30], who proposed an analysis of the ultimate load obtaining formulae for flange effective width ratio. Moreover, the idea of effective flange width has been featured in many design codes and standards, such as Eurocode 4 [31], the AASHTO LRFD Bridge Design Specifications [32], European Committee for Standardization (CEN) [33] and British Standards Institution (BSI) BS 5400 [34], where specific formulas for load type are provided.
A drawback in using effective width lies in the fact that it depends on the load distribution. Therefore, the effective width of a beam loaded by forces uniformly distributed on its whole length is different from that relevant to a segment of distributed loads, or even from that related to a concentrated load, applied, e.g., at different abscissas. Indeed, the standards recalled above give formulas valid for a limited number of different load conditions, since a general formulation is missing in the scientific literature. On the other side, the authors of this paper developed, in [35], a simple analytical model supplying a closed-form expression for a simply supported beam experiencing a sinusoidal deflection of half-wavelength equal to an integer submultiple n of the span. The tool was used to study free vibrations of a box-girder via the adapted classic model (Euler–Bernuolli, Timoshenko), accounting for the shear lag. The idea of this paper is to use such an n-depending effective to analyze static problems. Indeed, since any load can be expanded in a sinus Fourier series, it is possible to associate a different effective width to each of the harmonic components. By superimposing the effects, displacements and stresses, triggered by the original load, can be computed. The primary objective of this work is to evaluate, in each cross-section of the box-girder, the maximum value of the normal and tangential stresses. In particular, in each cross-section, the maximum normal stress occurs at the nodes, while maximum tangential stresses occur at the nodes of the flange and in the middle of the web. The proposed procedure, thanks to its user-friendliness, is intended as a valid alternative to consider the shear-lag effect in bridge designs.
Always in the framework of the Fourier analysis, the evaluation of displacements by using the Euler–Bernoulli beam model with the n-dependent inertia moment is proposed. Displacements by Euler–Bernoulli and Timoshenko beams with integer geometric characteristics are also illustrated for comparison.
Finally, the idea of an average effective width, which by still keeping the necessary dependence on the load shape avoids the need to employ a series expansion, is suggested here. This average effective width is intended as a medium value of the wavelength-dependent effective widths, and it can be obtained by imposing the equivalence of the displacements, or of the normal stresses, in a significant point of the beam.
The paper is organized as follows. In Section 2, the problem is defined, and the analytical expression derived in [35] for the effective flange width is recalled. In Section 3, making use of the reduced geometric characteristics, analytical expressions for both normal and tangential stresses are derived and discussed. In Section 4, always in the framework of the Fourier analysis, the deflections of the beam are evaluated by three different models of increasing refinement. In Section 5, the analytical formula for the average effective width for different load conditions is obtained and compared to the literature. In Section 6, displacements and stresses of a box-girder taken from the literature under different load conditions are evaluated by the present method and numerical results are discussed vs. Finite Element analyses. In Section 7, some conclusions are drawn.

2. Effective Flange Width

The formulation developed in [35] to obtain the effective flange width of rectangular box-girders undergoing sinusoidal deflections is summarized here.
A double-symmetric box-girder of dimensions b 1 × b 2 and span l , illustrated in Figure 1, is considered. The box-girder comprises four thin-walled elements: the bottom and upper flanges, represented by i = 1 , 3 and having a thickness h 1 = h 3 , and the right and left webs, represented by i = 2 , 4 , and having a thickness h 2 = h 4 . The box-girder is referenced to the cross-section inertia principal axes x , y at the centroidal axis z 0 , l . The curvilinear abscissa are represented by s i b i h 2 2 , b i h 2 2 along the flanges and by s i b i h 1 2 , b i h 1 2 along the webs. Furthermore, the box-girder is simply supported and torsionally clamped at the ends. Additionally, the cross-section of the box-girder is assumed to remain undeformable within its own plane. In-plane displacements of the cross-section, along the x , y axes, are denoted by u , v respectively and out-of-plane displacements by w.
A simplified analytical model based on the following hypotheses is used: (a) the vertical deflection of the box-girder is sinusoidal; (b) the flange and the web are inextensible in the cross-section plane; and (c) the web is shear undeformable. It is assumed that the box-girder deflects in the y z plane with a sinusoidal law
v z = v ^ sin n π l z
of amplitude v ^ . By virtue of the hypotesis (b), it is u z , s i = 0 . The out-of-plane displacement of the bottom flange is taken as w z , s 1 = w ^ s 1 cos n π l z , with aplitude w ^ s 1 to be determined by equilibrium conditions together with symmetry considerations (the amplitude should be symmetric with respect to the y axis).
The membrane strains on the flange are written, and then, the active stresses in the flange are derived by the elastic law (taking zero, for the moment, the Poisson ratio). By enforcing the indefinite equilibrium in the z direction and integrating ordinary differential equations, the normal stresses are evaluated. Finally, the result of these abscissa-dependent stresses acting on the flange is equated to the product of the maximum stress occurring at the flange–web junction by the unknown effective width b e n . The following result is drawn:
b e n = 2 β tanh b 1 2 β
where
β : = n π l E G
is the wave number. The shear-lag coefficient α n = b e n b 1 is defined as follows:
α n = 2 n π b 1 l E G tanh n π b 1 2 l E G
in which E and G are the Young modulus and the tangential elastic modulus, respectively. This is a ready-for-use formula, giving the effective width as a function of (i) the width-to-length ratio of the box-girder, (ii) the material property ratio E / G = 2 1 + ν , with ν the Poisson coefficient, and (ii) the number n of halfwaves which the beam deflects. For a beam of given geometric and mechanical characteristics, this latter is the only variable. The shear-lag coefficient in Equation (4) is plotted in Figure 2a vs. n, for a fixed b 1 l and different values of ν . It is noticed that the material parameter weakly affects the results. Figure 2b shows the shear-lag coefficient vs. n, for different values of b 1 l . For fixed b 1 l , the graph makes it possible to associate with each n the corresponding value of b e n b 1 and therefore the known b 1 , the relevant b e n .

3. State of Stress

The statics of a simply supported box-girder subjected to vertical loads is studied here, with the purpose of evaluating the maximum normal and tangential stresses.

3.1. Load and Solicitation

The load p z (positive upwards) is expanded in Fourier sine series as
p z = n p ^ n sin n π l z
where the Fourier coefficient p ^ n , representing the amplitude of the n-th harmonic, is given by
p ^ n = 2 l 0 l p z sin n π l z d z
To evaluate the shear force T z and the bending moment M z , the equilibrium equations T z + p z = 0 , M z + T z = 0 (shear force positive counterclockwise, bending moment positive when inducing upward concavity) must be integrated, with the boundary conditions M 0 = M l = 0 . Here, a dash denotes differentiation with respect to z. Since the beam is isostatically constrained, this problem admits a unique solution, namely
T z = n T ^ n cos n π l z M z = n M ^ n sin n π l z
where the solicitation-amplitudes T ^ n , M ^ n are related to the load-amplitudes by
T ^ n = l n π p ^ n
M ^ n = l n π 2 p ^ n
As an example, three cases are dealt with ahead (loads negative, since downward), namely, in Figure 3, (i) the uniformly distributed load p z = p 0 d = const ; (ii) the segment of uniform load p 0 s distributed on an interval ξ Δ 2 , ξ + Δ 2 , i.e., p z = p 0 s W ξ Δ 2 , ξ + Δ 2 , where W · is the window function; and (iii) the concentrated load P applied at the abscissa ξ , i.e., p z = P δ z ξ , where δ · is the Dirac function. In the three cases, the Fourier coefficients read, respectively,
p ^ n = 2 p 0 d 1 1 n n π
p ^ n = 4 p 0 s n π sin Δ n π 2 l sin ξ n π l
p ^ n = 2 l P sin ξ n π l
By using Equation (8), the relevant solicitation amplitudes are evaluated.
Both normal and tangential stress fields are obtained. These stresses are obtained here by using the De Saint Venant theory, in particular, the Navier theory is adopted for the normal stresses and the Jourawsky theory for the tangential ones.

3.2. Stresses

The normal stresses σ z z , s and the tangential stresses τ z s z , s , respectively, induced by the bending moment M z and the shear force T z , are evaluated according to the technical theory of the beam, which extends the De Saint Venant results to beams loaded along their axis. Consequently, the Navier and the Jourawsky formulas are assumed to hold. However, to account for the shear-lag effects, reference must be made to the effective cross-section, in which the true width b 1 is replaced by the effective width b e n = b e n , as illustrated in Figure 4. Due to the dependence on n, it is no possible to define an effective width of the flange valid for any loads, but reference has to be made to the single harmonic component of solicitation, i.e., M ^ n sin n π l z and T ^ n cos n π l z (which, indeed, trigger a sinusoidal deflection of the beam in n halfwaves). For each of these harmonics, the Navier and Jourawsky formulas are applied to a different geometry of the cross-section, and the stresses are added together to obtain
σ z z , s = n σ ^ n s sin n π l z
τ z s z , s = n τ ^ n s cos n π l z
where
σ ^ n s = M ^ n I n y , τ z s s = T ^ n I n h s S n * s
Here,
I n = 2 b 2 h 1 2 2 b e h 1 + 1 12 h 2 b 2 2 h 1 3 + 1 12 b e h 1 3
is the second-order inertia moment (with respect to the x axis) of the effective rectangular box-girder associated with the n-th harmonic of the solicitation, i.e., to flanges of width b e n = b e n . Similarly, S n * s is the static moment of the portion of effective cross-section cuts by a chord at the abscissa s, reading
S n * s 1 , 2 = s ˜ 1 h 1 b 2 h 1 2 , s ˜ 1 0 , b e n h 2 2 h 1 b 2 h 1 b e 2 h 2 4 h 2 8 b 2 2 4 b 2 h 1 + 4 h 1 2 + s 2 2 2 , s 2 b 2 h 1 2 , b 2 h 1 2 ,
respectively, on the flange and on the web. Here, s ˜ 1 = s 1 b 1 b e n 2 is a new translated abscissa.
The maximum (minimum) normal stress occurs at the top (bottom) of the cross-section, where, by using Equation (8b) in the Navier formula, it assumes the value
σ ^ n , m a x = l n π 2 b 2 2 I n p ^ n
The relative maximum (in modulus) tangential stresses occur (i) on the flange, at the junction with the web ( τ ^ n , m a x f l ), i.e., at s ˜ 1 = b e n 2 h 2 , and (ii) on the web at the midpoint ( τ ^ n , m a x w ), i.e., at s 2 = 0 . By using Equation (8a) in the Jourawsky formula, they assume the values
τ ^ n , m a x f l = l n π b 2 h 1 b e n 2 h 2 4 I n p ^ n τ ^ n , m a x w = l n π 1 4 I n h 1 b 2 h 1 b e 2 h 2 h 2 + b 2 2 h 1 2 2 p ^ n

4. Displacements

To evaluate the deflections of the beam, three different models are used, of increasing refinement, and always in the framework of the Fourier analysis. The models are
  • The total Euler–Bernoulli model ( E B t ), in which the shear lag is ignored, so that the cross-section is the integer; this model is used only for comparison purposes, since it is expected it overestimates the stiffness of the beam;
  • The total Timoshenko model ( T M t ), in which the shear lag is still ignored, but the shear deformability of the web is accounted for by taking a shear area approximately equal to the web area, namely A * = A w e b ; also this model is expected to give better results of the former, but they are still unsatisfactory;
  • The effective Euler–Bernoulli model ( E B e ), in which the shear-lag is included by exploiting the concept of wavelength-dependent effective width; this model is believed to be sufficiently accurate.
The shear lag could also be introduced in the Timoshenko model as, e.g., done in [35] in a dynamic context. There, it was proven that the Timoshenko model is more accurate than the Euler–Bernoulli model, when considering modes of short wavelength. In statics, however, since regular distribution of equiverse forces, which are of interest in bridge analysis, admit Fourier series in which higher harmonics bring a small contribution, the difference between the two models can be checked to be negligible, as proven ahead in the numerical analysis.
The E B t model is considered first. The equation of the elastic curve and boundary conditions of the simply supported beam read
E I v z = n p ^ n sin n π l z v 0 = v l = 0 v 0 = v l = 0
where v z is the transverse displacement at the abscissa z, I is the second-order inertia moment of the integer cross-section and use has been made of Equation (5) to express the load. After integration, the following displacement field is found:
v z = l 4 π 4 E I n 1 n 4 p ^ n sin n π l z
Concerning, the E B e model, a different inertia moment must be taken for any n, so that the nth problem reads
E I n v n z = p ^ n sin n π l z v n 0 = v n l = 0 v n 0 = v n l = 0
with v n z the response to the nth harmonic load component. From these equations, v n ( z ) = v ^ n sin n π l z is found, with
v ^ n = p ^ n E I n l n π 4
so that the response to p z , i.e., n v ^ n sin n π l z , is
v z = l 4 π 4 E n 1 n 4 I n p ^ n sin n π l z
The solutions in Equations (17) and (20) are similar, but while the inertia I is the same in each term of the former series, the inertia I n decreases with the number n in each term of the latter series. Therefore, the E B e is softer than the E B t model.
Finally, the T M t model is addressed. The equilibrium equations in terms of the (independent) transverse displacement v z and rotation φ z are (see, e.g., [35])
G A * v φ = n p ^ n sin n π l z E I φ G A * v φ = 0
with the boundary conditions
v 0 = v l = 0 φ 0 = φ l = 0
and symbols already introduced. These equations are exactly solved by the series as follows:
v z = n v ^ n sin n π l z φ z = n φ ^ n cos n π l z
whose coefficients are found to be
v ^ n = p ^ n E I l n π 4 1 + E I G A * n π l 2 φ ^ n = p ^ n E I l n π 3
Equation (24) shows that G A * < amplifies the deflection, when compared to the E B e model (which corresponds to G A * = ).

5. Average Effective Width

The theory developed so far explained that the effective width b e n is strictly related to the load shape as well as the wavelengths of its harmonic components. However, one might want to obtain an approximate formula which, still keeping the necessary dependence on the load shape, avoids the need to employ a series expansion, which entails the simple but tedious addition of a larger number of contributions. In this perspective, the idea of an average effective width, obtained as a medium value of the wavelength-dependent effective widths, would be desirable.

5.1. Displacement Equivalence

To obtain such an average value, the following strategy is followed: (a) for a given load p z , the deflection at the midspan (or at another significant point) is computed via the proposed algorithm, through the series expansion Equation (20), i.e.,
v l 2 = l 4 π 4 E n 1 n 4 I n p ^ n sin n π 2 ;
(b) such a the displacement is equated to v ¯ l 2 ,experienced by a beam whose cross-section has flanges of average effective width b ¯ e , i.e., average inertia moment I ¯ ; (c) from the equality v ¯ l 2 = v l 2 , the unknown b ¯ e is derived.
The procedure is explained in detail for the case of uniformly distributed load p z = p 0 d = const , for which p ^ n is given by Equation (9a). Since the midspan displacement of a beam of inertia moment I ¯ , under the same load, is v ¯ l 2 = 5 384 p 0 d l 4 E I ¯ , by enforcing the equality, the unknown I ¯ is derived as
I ¯ = 5 768 π 5 n 1 I n 1 1 n n 5 sin n π 2
Such an expression is still too complicated to draw b ¯ e , so that another strong simplification is introduced. It is assumed (only to the scope at hand) that the inertia moment of the reduced cross sections is almost equal to that of the two flanges (i.e., by neglecting the web contributions), whose area is concentrated in their centroids. Accordingly, I ¯ 1 2 b ¯ e h 1 b 2 h 1 2 and I n = 1 2 b e n h 1 b 2 h 1 2 so that the previous equation provides
b ¯ e b 1 = 5 768 π 5 n 1 1 n n 5 α n sin n π 2
in which b e n = α n b 1 has been introduced. By remembering Equation (4) for α n , it is concluded that b ¯ e b 1 depends on the aspect ratio b 1 l and (weakly) on the Poisson ratio. By fixing the latter, a plot of α ¯ = b ¯ e b 1 vs. b 1 l can be constructed (see Figure 5a), which is ready for use. Indeed, from known b ¯ e , the (now not approximated) inertia moment I ¯ can be computed, and stresses and displacements evaluated, without resorting to series expansions. The accuracy of this simplified method is tested in the next section.
By repeating the procedure for the concentrated load applied at the midspan (Figure 5a), it is found that
b ¯ e b 1 = 1 96 π 4 n 1 n 4 α n sin 2 n π 2
and for a segment of load of extension Δ centered at the midspan (Figure 5b),
b ¯ e b 1 = 1 1536 π 5 8 l 3 Δ 4 l Δ 3 + Δ 4 l 4 n 1 n 5 α n sin Δ n π 2 l sin 2 n π 2

5.2. Normal Stresses Equivalence

It is also possible to obtain an average value of the effective width by using the equivalence of the normal stresses instead of the displacements (Section 5.1). For this purpose, the following strategy, similar to that pursued in Section 5.1, is followed: (a) for a given load p z , the maximum normal stress at the midspan (or at another significant point) is computed via the proposed algorithm, through the series expansion Equation (10a), i.e.,
σ z l 2 , b 2 2 = l 2 b 2 2 π 2 n 1 n 2 I n p ^ n sin n π 2
(b) such a normal stress is equated to σ ¯ z l 2 , b 2 2 ,experienced by a beam whose cross-section has flanges of average effective width b ¯ e , i.e., average inertia moment I ¯ ; (c) from the equality σ ¯ z l 2 , b 2 2 = σ z l 2 , b 2 2 , the unknown b ¯ e is derived.
In the three cases, uniformly distributed load, segment of load centered at the midspan and concentrated load applied at the midspan, the normal stresses σ ¯ z l 2 , b 2 2 read, respectively,
σ ¯ z l 2 , b 2 2 = p 0 d l 2 b 2 16 I ¯ σ ¯ z l 2 , b 2 2 = p 0 s Δ 2 l Δ b 2 16 I ¯ σ ¯ z l 2 , b 2 2 = P l b 2 8 I ¯
In the three cases, it is found, respectively, that
b ¯ e b 1 = 1 16 π 3 n 1 1 n n 3 α n sin n π 2 b ¯ e b 1 = 1 32 π 3 Δ 2 l Δ l 2 n 1 n 3 α n sin Δ n π 2 l sin 2 n π 2 b ¯ e b 1 = 1 8 π 2 n 1 n 2 α n sin 2 n π 2
Figure 5 shows the average shear-lag coefficient α ¯ = b ¯ e b 1 vs. b 1 l for the uniformly distributed load, concentrated load at the midspan and segment of distributed load, centered at midspan, for different extensions Δ l = 1 10 , 1 4 , 1 2 , by both the displacement equivalence (dashed line) and the stresses equivalence (continuous line). Results obtained by the displacement equivalence does not appear able to take into account the loading conditions, therefore it is considered to be not sufficiently accurate.
The shear-lag coefficient for the two load cases, concentrated force at midspan and uniformly distributed load, have been studied by the BS5400 [34] and by a variational analysis in [27], where the graph of the shear-lag coefficient b e b 1 vs. b 1 2 l is obtained numerically. Figure 6 shows results by BS5400, by [27] and by the present study all together, for the two load cases mentioned above. They are very close to each other; however, it is important to underline that our study, unlike the other studies mentioned, provides an analytical formula for the relationship, and therefore, it is applicable to different load conditions, and this character of generality constitutes an added value.
In support of the above, the calculation of the average shear-lag coefficient using the equivalence of the normal stresses for two load conditions not present in the literature is proposed here: a concentrated force applied at ξ = l 4 and a segment of distributed load, centered at the abscissa ξ = l 4 . By following the procedure already described for the simplest loading conditions, in the new load cases considered, the normal stresses σ ¯ z l 4 , b 2 2 read, respectively,
σ ¯ z l 4 , b 2 2 = 3 32 P l b 2 I ¯ σ ¯ z l 4 , b 2 2 = p 0 s Δ 3 l 2 Δ b 2 32 I ¯
So it is found, respectively, that
b ¯ e b 1 = 3 32 π 2 n 1 n 2 α n sin 2 n π 4 b ¯ e b 1 = 1 64 π 3 Δ 3 l 2 Δ l 2 n 1 n 3 α n sin Δ n π 2 l sin 2 n π 4
Figure 7 shows the average shear-lag coefficient α ¯ vs. b 1 l for the concentrated load applied at ξ = l 4 and for the segment of the distributed load, centered at the abscissa ξ = l 4 , for different extensions Δ l = 1 10 , 1 4 , 1 2 .

6. Numerical Results

The theory so far developed is applied, for validation purposes, to perform the static analysis of sample simple-supported box-girders in bending, for which the analytical results are compared with those supplied by Finite Element analyses. A case study is considered, consisting of a reinforced concrete box-girder bridge taken from the literature [36]. Three different load cases are analyzed: (i) a uniformly distributed load p 0 d (DL); (ii) a segment of uniformly distributed load p 0 s (SL) on a length Δ = l 4 , centered at the abscissa ξ = l 4 and (iii) a concentrated load P (CL) applied at ξ = l 4 .
The FEM analysis utilized a commercial software, namely SAP2000. The structural representation of the box bridge involves the assembly of thin-shell elements. To realize the boundary conditions of the analytical model, the in-plane displacements of all mesh nodes are constrained at the end cross-sections. Additionally, at one end, the longitudinal displacement of two points along the x-axis is restricted and ’diaphragm constraints’ at every cross-section of the mesh are introduced to simulate the non-deformability of all the cross-sections of the girder, allowing for unrestricted flexural warping.
Regarding the application of loads in the FE model, it is important to note that it exclusively employs thin-shell elements. Consequently, loads can only be applied at joints as concentrated loads using the ‘Assign joint loads’ command. The selected joints for applying forces are located at the intersections of webs (elements 2, 4) and the top flange (element 3) as shown in Figure 8. This choice is justified by the usual procedure carried out in the bridge analysis, in which the loads applied to the flange are “brought to the nodes” by performing first a “local effect analysis”, in which the nodes are fictitiously supported, and successively a “global effect analysis”, in which the reactive forces, changed in sign, are applied to the nodes, freed from constraints.
Chosen p 0 d = p 0 s = p 0 = 100 kN / m , P = 100 kN , the loads DL, SL and CL are expressed as sinusoidal Fourier series and depicted in Figure 9 for N = 20 , 50 , 100 . To determine an appropriate value for N, Fourier analyses are conducted for increasing values of this number. The Fourier series with N = 20 inadequately approximates the loads, particularly in the SL case and, even worse, in the CL case (see Figure 9). For the last case, in particular, it was deemed useful to carry out a convergence analysis of P ˜ = 0 l p z d z (Figure 9d). Consequently, N = 50 is ultimately selected for the analysis.

6.1. Case Study

A rectangular box-girder bridge, made of reinforced concrete, is considered as a case study, as detailed in [36]. The geometric specifications of the box-girder are as follows: length l = 30 m , width b 1 = 6.35 m , height b 2 = 1.75 m , flange thickness h 1 = 0.25 m , and web thickness h 2 = 0.35 m (refer to Figure 1). The elastic modulus is E = 35,654 × 10 3 kN / m 2 , and the Poisson’s coefficient is ν = 0.15 . The mass per unit volume is calculated as ρ = 25,000 9.81 = 2548.42 kg / m 3 .
The Finite Element analysis involved the utilization of a model comprising 2400 shell elements (40 along the directrix and 60 along the axis). This resulted in a total of 2440 degrees of freedom. In this FE model, the side of the mesh is d = 0.5 m , so for the uniformly distributed load (DL), the number of joints, in correspondence with each web, where the force must be applied is n D L = ( l d + 1 ) = 61 , and the modulus of any concentrated force is p 0 × l 2 n D L = 24.59 kN (see Figure 8a). For the segment of uniformly distributed load (SL), a similar but different mode of working, based on the concept of area of influence, is chosen as follows: at n S L = Δ d = 15 joints in correspondence with each web (the first located at z = 4 m ), the concentrated force p 0 × Δ 2 n S L = 25 kN is applied (see Figure 8b). The concentrated load is applied as a concentrated force P 2 = 50 kN at the joint located at z = 7.5 m in correspondence with each web (see Figure 8c).

6.1.1. Stress Analysis

Maximum normal and tangential stresses (Equations (14) and (15)), at different cross-sections, are plotted in Figure 10, Figure 11 and Figure 12 for the uniformly distributed load case (DL). In particular, Figure 10b–Figure 12b show the maximum stresses for different N both by analytical and numerical models (taken as a reference), revealing a fast convergence, also for low n.
In order to set the sufficient number of harmonics to use, an analysis of harmonics for the stresses on z is carried out. Figure 10a–Figure 12a indicate that starting from n = 5 , the influence of the pertinent harmonic becomes negligible. Consequently, it can be inferred that conducting a Fourier analysis with an N = 50 would be futile. As previously elucidated, the decision to employ a Fourier series with N = 50 has been made to effectively approximate the loads, particularly the loaded segment.
Figure 13, Figure 14 and Figure 15 are about the segment of uniformly distributed load (SL). Similar considerations to the previous case can be made.
Figure 16, Figure 17 and Figure 18 are about the concentrated load (CL). In this case, it is particularly important to consider N = 50 instead of smaller N to capture a stress value sufficiently close to the FE model, especially in correspondence with the concentrated load (see Figure 16b) with a maximum error of about 5% for the maximum normal stresses and 12% for the maximum tangential stresses.

6.1.2. Displacement Analysis

The vertical displacement v z by the three analytical models considered (total Euler-Bernoulli, total Timoshenko and effective Euler-Bernoulli) are plotted in Figure 19 with those by the FEM analysis, for the three load cases. As expected, for all three cases analyzed, the E B e model, in which the shear-lag is included by exploiting the concept of wavelength-dependent effective width, is the most consistent with the FEM results with a maximum error of about 0% for the DL, 1.6% for the SL and 2.8% for the CL. Furthermore, E B t and T M t model present respectively a maximum error of about 7% and 3.6% for the DL, 8.7% and 2% for the SL, and 8.9% and 5.3% for the CL.

6.1.3. Use of the Average Effective Width

The goodness of the results obtained by using the effective width b e n has already been tested by comparison with FEM results, so it is now intended to compare the normal stresses obtained by considering the effective width b e n by Equation (2), taken as a reference, with those by considering the average effective width b ¯ e derived both from the displacement equivalence and the normal stress equivalence (see Section 5). For this comparison, the three simple loading conditions already mentioned in Section 5 are considered, in particular, the uniformly distributed load p 0 d = 100 kN / m (DL), the segment of load p 0 s = 100 kN / m of extension Δ = l 10 centered at midspan (SL-centered) and concentrated load P = 100 kN applied at midspan (CL-centered).
Figure 20 shows that for the DL case, the three curves are almost superimposed, while, on the contrary, for the CL-centered case, the two curves related to the average effective widths are notably different; in particular, deriving from the normal stresses equivalence captures the maximum value quite well, while deriving from the displacement equivalence follows the curve better but is not able to approximate the maximum value. The case of SL-centered is an intermediate case. The amount of the difference between the two approximate curves depends on the extension Δ of the segment loaded.

7. Conclusions

A new method has been proposed for the bending static analysis of simply supported, but generally loaded, rectangular box-girders with undeformable cross-section, undergoing shear-lag. The method is based on (i) expanding the general load in the Fourier series; (ii) finding the effective flange width relevant to each harmonic; (iii) computing, for each harmonic, displacements and stresses by using a beam model with effective geometric characteristics; and (iv) combining the results of all the harmonics. As an alternative to the effects combination, approximate formulae giving ready-for-use expressions for an “average effective width” have been derived, based on equivalences at selected points with (i) displacements or (ii) stresses of the more refined model. Results have been compared with much more expensive Finite Element analyses to validate the method presented. The following conclusions, relevant to sample loading conditions, have been drawn.
  • The maximum displacement and stress supplied by the algorithm are in excellent agreement with the FE results, even in the most complex load cases analyzed.
  • The average shear-lag coefficient based on the stress equivalence, gives excellent results in accordance with the literature. Differently from the literature, however, they can be generalized to any load conditions.
  • The displacement-based equivalence generally gives less accurate results. In some circumstances, it describes better the longitudinal trend, while the stress-based equivalence captures the maximum values better.
  • The neglection of the shear deformation of the web, which could be accounted by the Timoshenko beam model, has been proven to not significantly affect the results, as it emerges from the comparison with shear deformable FE.
In conclusion, based on the excellent agreement of the results with FEM analyses, it is believed that the formulae here derived for the effective width flange of simply supported box girders are reliable and usable. Above all, a methodological and innovative algorithm has been developed, easily implementable for any load distributions of interest to the designer.

Author Contributions

Conceptualization, A.L.; methodology, A.L.; software, F.P.; validation, F.P.; formal analysis, A.L. and F.P.; investigation, F.P.; resources, A.L. and F.P.; data curation, F.P.; writing—original draft preparation, F.P.; writing—review and editing, A.L. and F.P.; visualization, F.P.; supervision, A.L.; project administration, A.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Meng, X.; Xiao, F.; Yan, Y.; Chen, G.S.; Ma, Y. Non-destructive damage evaluation based on static response for beam-like structures considering shear deformation. Appl. Sci. 2023, 13, 8219. [Google Scholar] [CrossRef]
  2. Xiao, F.; Zhu, W.; Meng, X.; Chen, G.S. Parameter identification of frame structures by considering shear deformation. Int. J. Distrib. Sens. Netw. 2023, 2023, 6631716. [Google Scholar] [CrossRef]
  3. Luongo, A.; D’Annibale, F.; Ferretti, M. Shear and flexural factors for static analysis of homogenized beam models of planar frames. Eng. Struct. 2021, 228, 111440. [Google Scholar] [CrossRef]
  4. Dezi, L.; Mentrasti, L. Nonuniform bending-stress distribution (shear lag). J. Struct. Eng. 1985, 111, 2675–2690. [Google Scholar] [CrossRef]
  5. Chang, S.T.; Zheng, F.Z. Negative shear lag in cantilever box girder with constant depth. J. Struct. Eng. 1987, 113, 20–35. [Google Scholar] [CrossRef]
  6. Singh, G.J.; Mandal, S.; Kumar, R.; Kumar, V. Simplified analysis of negative shear lag in laminated composite cantilever beam. J. Aerosp. Eng. 2020, 33, 04019103. [Google Scholar] [CrossRef]
  7. Li, X.; Wan, S.; Mo, Y.; Shen, K.; Zhou, T.; Nian, Y. An improved method for analyzing shear lag in thin-walled box-section beam with arbitrary width of cantilever flange. Thin-Walled Struct. 2019, 140, 222–235. [Google Scholar] [CrossRef]
  8. Reissner, E. Analysis of shear lag in box beams by the principle of minimum potential energy. Q. Appl. Math. 1946, 4, 268–278. [Google Scholar] [CrossRef]
  9. Luo, Q.Z.; Tang, Q.; Li, Q.; Liu, G.; Wu, J. Membrane forces acting on thin-walled box girders considering shear lag effect. Thin-Walled Struct. 2004, 42, 741–757. [Google Scholar] [CrossRef]
  10. Lin, Z.; Zhao, J. Least-work solutions of flange normal stresses in thin-walled flexural members with high-order polynomial. Eng. Struct. 2011, 33, 2754–2761. [Google Scholar] [CrossRef]
  11. Lin, Z.; Zhao, J. Modeling inelastic shear lag in steel box beams. Eng. Struct. 2012, 41, 90–97. [Google Scholar] [CrossRef]
  12. Wu, Y.; Liu, S.; Zhu, Y.; Lai, Y. Matrix analysis of shear lag and shear deformation in thin-walled box beams. J. Eng. Mech. 2003, 129, 944–950. [Google Scholar] [CrossRef]
  13. Lertsima, C.; Chaisomphob, T.; Yamaguchi, E. Stress concentration due to shear lag in simply supported box girders. Eng. Struct. 2004, 26, 1093–1101. [Google Scholar] [CrossRef]
  14. Sa-nguanmanasak, J.; Chaisomphob, T.; Yamaguchi, E. Stress concentration due to shear lag in continuous box girders. Eng. Struct. 2007, 29, 1414–1421. [Google Scholar] [CrossRef]
  15. Luo, Q.Z.; Wu, Y.; Tang, J.; Li, Q. Experimental studies on shear lag of box girders. Eng. Struct. 2002, 24, 469–477. [Google Scholar] [CrossRef]
  16. Moffatt, K.R.; Dowling, P.J. British shear lag rules for composite girders. J. Struct. Div. 1978, 104, 1123–1130. [Google Scholar] [CrossRef]
  17. Evans, H.R.; Taherian, A.R. The prediction of the shear lag effect in box girders. Proc. Inst. Civ. Eng. 1977, 63, 69–92. [Google Scholar] [CrossRef]
  18. Wang, X.; Rammerstorfer, F.G. Determination of effective breadth and effective width of stiffened plates by finite strip analyses. Thin-Walled Struct. 1996, 26, 261–286. [Google Scholar] [CrossRef]
  19. Kristek, V.; Evans, H.R. A hand calculation of the shear lag effect in unstiffened flanges and in flanges with closely spaced stiffeners. Civ. Eng. Pract. Des. Eng. 1985, 4, 163–190. [Google Scholar]
  20. Qin, X.; Liu, H. Effective flange width of simply supported box girder under uniform load. Acta Mech. Solida Sin. 2010, 23, 57–65. [Google Scholar] [CrossRef]
  21. Luo, Q.Z.; Li, Q.; Tang, J. Shear lag in box girder bridges. J. Bridge Eng. 2002, 7, 308–313. [Google Scholar] [CrossRef]
  22. Luo, Q.Z.; 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]
  23. Zhou, S.J. Finite beam element considering shear-lag effect in box girder. J. Eng. Mech. 2010, 136, 1115–1122. [Google Scholar] [CrossRef]
  24. Lee, C.K.; Wu, G. Shear lag analysis by the adaptive finite element method: 1. Analysis of simple plated structures. Thin-Walled Struct. 2000, 38, 285–309. [Google Scholar] [CrossRef]
  25. Fafitis, A.; Rong, A.Y. Analysis of thin-walled box girders by parallel processing. Thin-Walled Struct. 1995, 21, 233–240. [Google Scholar] [CrossRef]
  26. Chen, J.; Shen, S.L.; Yin, Z.Y.; Horpibulsuk, S. Closed-form solution for shear lag with derived flange deformation function. J. Constr. Steel Res. 2014, 102, 104–110. [Google Scholar] [CrossRef]
  27. Lin, Z.; Zhao, J. Revisit of AASHTO effective flange-width provisions for box girders. J. Bridge Eng. 2011, 16, 881–889. [Google Scholar] [CrossRef]
  28. Von Karman, T. Die Mittragende Poreite; Springer: Berlin/Heidelberg, Germany, 1924. [Google Scholar]
  29. Sedlacek, G.; Bild, S. A simplified method for the determination of the effective width due to shear lag effects. J. Constr. Steel Res. 1993, 24, 155–182. [Google Scholar] [CrossRef]
  30. Wu, Y.; Yu, S.; Shi, C.; Li, J.; Lai, Y.; Zhu, Y. Ultimate load analysis of thin-walled box beams considering shear lag effect. Thin-Walled Struct. 2004, 42, 1199–1210. [Google Scholar] [CrossRef]
  31. EN 1994-2; EUROCODE 4: Design of Composite Steel and Concrete Structures—Part 2: General Rules and Rules for Bridges. European Committee for Standardization: Brussels, Belgium, 2005.
  32. American Association of State Highway and Transportation Officials. AASHTO LRFD Bridge Design Specifications, 1st ed.; American Association of State Highway and Transportation Officials: Washington, DC, USA, 2008. [Google Scholar]
  33. CEN (European Committee for Standardization). Design of Concrete Structures—Part 1-1: General Rules and Rules for Buildings; Eurocode 2; CEN (European Committee for Standardization): Brussels, Belgium, 2004. [Google Scholar]
  34. British Standards Institution. BS5400, Steel, Concrete and Composite Bridges, Part 3, Code of Practice for Design of Steel Bridges; British Standards Institution: London, UK, 1982. [Google Scholar]
  35. Pancella, F.; Luongo, A. Shear-flexural beam models for free vibration analysis of warpable cross-undeformable box girders. Thin-Walled Struct. 2023, 188, 110837. [Google Scholar] [CrossRef]
  36. Petrangeli, M.; Zechini, A. Il calcolo dei ponti a cassone unicellulare con pareti sottili. G. Genio Civ. 1973, 2, 108–122. [Google Scholar]
Figure 1. Box-girder bridge: cross-section and curvilinear abscissa.
Figure 1. Box-girder bridge: cross-section and curvilinear abscissa.
Applsci 14 02848 g001
Figure 2. Shear lag coefficient α n = b e n b 1 vs. n: (a) for different materials Poisson coefficients ν = 0 , 0.2 , 0.3 and fixed b 1 l = 1 10 ; (b) for different b 1 l = 1 10 , 1 20 , 1 30 and fixed materials Poisson coefficients ν = 0.2 .
Figure 2. Shear lag coefficient α n = b e n b 1 vs. n: (a) for different materials Poisson coefficients ν = 0 , 0.2 , 0.3 and fixed b 1 l = 1 10 ; (b) for different b 1 l = 1 10 , 1 20 , 1 30 and fixed materials Poisson coefficients ν = 0.2 .
Applsci 14 02848 g002
Figure 3. Loading conditions: (a) uniformly distributed load, (b) segment of load, (c) concentrated load; load equally divided between the two webs in the next numerical simulations.
Figure 3. Loading conditions: (a) uniformly distributed load, (b) segment of load, (c) concentrated load; load equally divided between the two webs in the next numerical simulations.
Applsci 14 02848 g003
Figure 4. Normal and tangential stresses induced by the nth harmonic of the load.
Figure 4. Normal and tangential stresses induced by the nth harmonic of the load.
Applsci 14 02848 g004
Figure 5. Average shear-lag coefficient α ¯ = b ¯ e b 1 vs. the b 1 l aspect ratio: (a) uniformly distributed load (DL) and concentrated load at midspan (CLc); (b) segment of distributed load centered at midspan (SLc), for different extensions Δ l = 1 10 , 1 4 , 1 2 ; by the displacement equivalence (d.e.), by the stresses equivalence (s.e.). Poisson ratio ν = 0.15 .
Figure 5. Average shear-lag coefficient α ¯ = b ¯ e b 1 vs. the b 1 l aspect ratio: (a) uniformly distributed load (DL) and concentrated load at midspan (CLc); (b) segment of distributed load centered at midspan (SLc), for different extensions Δ l = 1 10 , 1 4 , 1 2 ; by the displacement equivalence (d.e.), by the stresses equivalence (s.e.). Poisson ratio ν = 0.15 .
Applsci 14 02848 g005
Figure 6. Average shear-lag coefficient α ¯ = b ¯ e b 1 vs. the b 1 2 l aspect ratio by the stresses equivalence: (a) concentrated load applied at the midspan; (b) uniformly distributed load. Poisson ratio ν = 0.15 . VA is for variation analysis by [27].
Figure 6. Average shear-lag coefficient α ¯ = b ¯ e b 1 vs. the b 1 2 l aspect ratio by the stresses equivalence: (a) concentrated load applied at the midspan; (b) uniformly distributed load. Poisson ratio ν = 0.15 . VA is for variation analysis by [27].
Applsci 14 02848 g006
Figure 7. Average shear-lag coefficient α ¯ = b ¯ e b 1 vs. the b 1 l aspect ratio by the stresses equivalence: (a) concentrated load applied at ξ = l 4 ; (b) segment of distributed load, centered at the abscissa ξ = l 4 , for different extensions Δ l = 1 10 , 1 4 , 1 2 . Poisson ratio ν = 0.15 .
Figure 7. Average shear-lag coefficient α ¯ = b ¯ e b 1 vs. the b 1 l aspect ratio by the stresses equivalence: (a) concentrated load applied at ξ = l 4 ; (b) segment of distributed load, centered at the abscissa ξ = l 4 , for different extensions Δ l = 1 10 , 1 4 , 1 2 . Poisson ratio ν = 0.15 .
Applsci 14 02848 g007
Figure 8. FE-Model: (a) DL; (b) SL; (c) CL.
Figure 8. FE-Model: (a) DL; (b) SL; (c) CL.
Applsci 14 02848 g008
Figure 9. Loads (in kN / m and kN ): (a) distributed load (DL) expanded in Fourier sine series; (b) segment of distributed load (SL) expanded in Fourier sine series; (c) concentrated load (CL) expanded in Fourier sine series; (d) convergence analysis of P ˜ vs. N. For (ac), N = 20 (dashed line), N = 50 (thin line), N = 100 (thick line).
Figure 9. Loads (in kN / m and kN ): (a) distributed load (DL) expanded in Fourier sine series; (b) segment of distributed load (SL) expanded in Fourier sine series; (c) concentrated load (CL) expanded in Fourier sine series; (d) convergence analysis of P ˜ vs. N. For (ac), N = 20 (dashed line), N = 50 (thin line), N = 100 (thick line).
Applsci 14 02848 g009
Figure 10. Maximum normal stress (in N / mm 2 ) vs. z: (a) n-th harmonic σ ^ n , m a x sin n π l z , n = 1 , 2 , 3 , 4 , 5 ; (b) by Equation (10a) for N = 1 , 2 , 3 , 4 , 5 , 50 (from the lightest line to the darkest). FEM analysis (bullets). Distributed load (DL).
Figure 10. Maximum normal stress (in N / mm 2 ) vs. z: (a) n-th harmonic σ ^ n , m a x sin n π l z , n = 1 , 2 , 3 , 4 , 5 ; (b) by Equation (10a) for N = 1 , 2 , 3 , 4 , 5 , 50 (from the lightest line to the darkest). FEM analysis (bullets). Distributed load (DL).
Applsci 14 02848 g010
Figure 11. Maximum tangential stress (in N / mm 2 ) vs. z on flange: (a) n-th harmonic τ ^ n , m a x f l cos n π l z , n = 1 , 2 , 3 , 4 , 5 ; (b) by Equation (10) for N = 1 , 2 , 3 , 4 , 5 , 50 (from the lightest line to the darkest). FEM analysis (bullets). Distributed load (DL).
Figure 11. Maximum tangential stress (in N / mm 2 ) vs. z on flange: (a) n-th harmonic τ ^ n , m a x f l cos n π l z , n = 1 , 2 , 3 , 4 , 5 ; (b) by Equation (10) for N = 1 , 2 , 3 , 4 , 5 , 50 (from the lightest line to the darkest). FEM analysis (bullets). Distributed load (DL).
Applsci 14 02848 g011
Figure 12. Maximum tangential stress (in N / mm 2 ) vs. z on web: (a) n-th harmonic τ ^ n , m a x w cos n π l z , n = 1 , 2 , 3 , 4 , 5 ; (b) by Equation (10b) for N = 1 , 2 , 3 , 4 , 5 , 50 (from the lightest line to the darkest). FEM analysis (bullets). Distributed load (DL).
Figure 12. Maximum tangential stress (in N / mm 2 ) vs. z on web: (a) n-th harmonic τ ^ n , m a x w cos n π l z , n = 1 , 2 , 3 , 4 , 5 ; (b) by Equation (10b) for N = 1 , 2 , 3 , 4 , 5 , 50 (from the lightest line to the darkest). FEM analysis (bullets). Distributed load (DL).
Applsci 14 02848 g012
Figure 13. Maximum normal stress (in N / mm 2 ) vs. z: (a) n-th harmonic σ ^ n , m a x sin n π l z , n = 1 , 2 , 3 , 4 , 5 ; (b) by Equation (10a) for N = 1 , 2 , 3 , 4 , 5 , 50 (from the lightest line to the darkest). FEM analysis (bullets). Segment of uniformly distributed load (SL).
Figure 13. Maximum normal stress (in N / mm 2 ) vs. z: (a) n-th harmonic σ ^ n , m a x sin n π l z , n = 1 , 2 , 3 , 4 , 5 ; (b) by Equation (10a) for N = 1 , 2 , 3 , 4 , 5 , 50 (from the lightest line to the darkest). FEM analysis (bullets). Segment of uniformly distributed load (SL).
Applsci 14 02848 g013
Figure 14. Maximum tangential stress (in N / mm 2 ) vs. z on flange: (a) n-th harmonic τ ^ n , m a x f l cos n π l z , n = 1 , 2 , 3 , 4 , 5 ; (b) by Equation (10b) for N = 1 , 2 , 3 , 4 , 5 , 50 (from the lightest line to the darkest). FEM analysis (bullets). Segment of uniformly distributed load (SL).
Figure 14. Maximum tangential stress (in N / mm 2 ) vs. z on flange: (a) n-th harmonic τ ^ n , m a x f l cos n π l z , n = 1 , 2 , 3 , 4 , 5 ; (b) by Equation (10b) for N = 1 , 2 , 3 , 4 , 5 , 50 (from the lightest line to the darkest). FEM analysis (bullets). Segment of uniformly distributed load (SL).
Applsci 14 02848 g014
Figure 15. Maximum tangential stress (in N / mm 2 ) vs. z on web: (a) n-th harmonic τ ^ n , m a x w cos n π l z , n = 1 , 2 , 3 , 4 , 5 ; (b) by Equation (10b) for N = 1 , 2 , 3 , 4 , 5 , 50 (from the lightest line to the darkest). FEM analysis (bullets). Segment of uniformly distributed load (SL).
Figure 15. Maximum tangential stress (in N / mm 2 ) vs. z on web: (a) n-th harmonic τ ^ n , m a x w cos n π l z , n = 1 , 2 , 3 , 4 , 5 ; (b) by Equation (10b) for N = 1 , 2 , 3 , 4 , 5 , 50 (from the lightest line to the darkest). FEM analysis (bullets). Segment of uniformly distributed load (SL).
Applsci 14 02848 g015
Figure 16. Maximum normal stress (in N / mm 2 ) vs. z: (a) n-th harmonic σ ^ n , m a x sin n π l z , n = 1 , 2 , 3 , 4 , 5 ; (b) by Equation (10a) for N = 1 , 2 , 3 , 4 , 5 , 50 (from the lightest line to the darkest). FEM analysis (bullets). Concentrated load (CL).
Figure 16. Maximum normal stress (in N / mm 2 ) vs. z: (a) n-th harmonic σ ^ n , m a x sin n π l z , n = 1 , 2 , 3 , 4 , 5 ; (b) by Equation (10a) for N = 1 , 2 , 3 , 4 , 5 , 50 (from the lightest line to the darkest). FEM analysis (bullets). Concentrated load (CL).
Applsci 14 02848 g016
Figure 17. Maximum tangential stress (in N / mm 2 ) vs. z on flange: (a) n-th harmonic τ ^ n , m a x f l cos n π l z , n = 1 , 2 , 3 , 4 , 5 ; (b) by Equation (10b) for N = 1 , 2 , 3 , 4 , 5 , 50 (from the lightest line to the darkest). FEM analysis (bullets). Concentrated load (CL).
Figure 17. Maximum tangential stress (in N / mm 2 ) vs. z on flange: (a) n-th harmonic τ ^ n , m a x f l cos n π l z , n = 1 , 2 , 3 , 4 , 5 ; (b) by Equation (10b) for N = 1 , 2 , 3 , 4 , 5 , 50 (from the lightest line to the darkest). FEM analysis (bullets). Concentrated load (CL).
Applsci 14 02848 g017
Figure 18. Maximum tangential stress (in N / mm 2 ) vs. z on web: (a) n-th harmonic τ ^ n , m a x w cos n π l z , n = 1 , 2 , 3 , 4 , 5 ; (b) by Equation (10b) for N = 1 , 2 , 3 , 4 , 5 , 50 (from the lightest line to the darkest). FEM analysis (bullets). Concentrated load (CL).
Figure 18. Maximum tangential stress (in N / mm 2 ) vs. z on web: (a) n-th harmonic τ ^ n , m a x w cos n π l z , n = 1 , 2 , 3 , 4 , 5 ; (b) by Equation (10b) for N = 1 , 2 , 3 , 4 , 5 , 50 (from the lightest line to the darkest). FEM analysis (bullets). Concentrated load (CL).
Applsci 14 02848 g018
Figure 19. Vertical displacement (in mm ) evaluated at half side of element 1 (at s 1 = 0 ): (a) uniformly distributed load (DL); (b) segment of uniformly distributed load (SL); (c) concentrated load (CL). Total Euler–Bernoulli (dashed line), total Timoshenko (thin line), effective Euler–Bernoulli (thick line), FEM (bullets).
Figure 19. Vertical displacement (in mm ) evaluated at half side of element 1 (at s 1 = 0 ): (a) uniformly distributed load (DL); (b) segment of uniformly distributed load (SL); (c) concentrated load (CL). Total Euler–Bernoulli (dashed line), total Timoshenko (thin line), effective Euler–Bernoulli (thick line), FEM (bullets).
Applsci 14 02848 g019
Figure 20. Maximum normal stress (in N / mm 2 ) vs. z: (a) DL; (b) SL-centered; (c) CL-centered. Evaluated with b e n (black continuous line), b ¯ e by the normal stresses equivalence (gray line), b ¯ e by the displacements equivalence (black dashed line).
Figure 20. Maximum normal stress (in N / mm 2 ) vs. z: (a) DL; (b) SL-centered; (c) CL-centered. Evaluated with b e n (black continuous line), b ¯ e by the normal stresses equivalence (gray line), b ¯ e by the displacements equivalence (black dashed line).
Applsci 14 02848 g020
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

Pancella, F.; Luongo, A. On the Use of the Effective Width for Simply Supported Generally Loaded Box-Girders. Appl. Sci. 2024, 14, 2848. https://doi.org/10.3390/app14072848

AMA Style

Pancella F, Luongo A. On the Use of the Effective Width for Simply Supported Generally Loaded Box-Girders. Applied Sciences. 2024; 14(7):2848. https://doi.org/10.3390/app14072848

Chicago/Turabian Style

Pancella, Francesca, and Angelo Luongo. 2024. "On the Use of the Effective Width for Simply Supported Generally Loaded Box-Girders" Applied Sciences 14, no. 7: 2848. https://doi.org/10.3390/app14072848

APA Style

Pancella, F., & Luongo, A. (2024). On the Use of the Effective Width for Simply Supported Generally Loaded Box-Girders. Applied Sciences, 14(7), 2848. https://doi.org/10.3390/app14072848

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