Understanding of the stresses induced in beams by bending loads took many years to develop. Galileo worked on this problem, but the theory as we use it today is usually credited principally to the great mathematician Leonard Euler (1707–1783). As will be developed below, beams develop normal stresses in the lengthwise direction that vary from a maximum in tension at one surface, to zero at the beam’s midplane, to a maximum in compression at the opposite surface. Shear stresses are also induced, although these are often negligible in comparision with the normal stresses when the length-to-height ratio of the beam is large. The procedures for calculating these stresses for various loading conditions and beam cross-section shapes are perhaps the most important methods to be found in introductory Mechanics of Materials, and will be developed in the sections to follow. This theory requires that the user be able to construct shear and bending moment diagrams for the beam, as developed for instance in Module 12.
A beam subjected to a positive bending moment will tend to develop a concave-upward curvature. Intuitively, this means the material near the top of the beam is placed in compression along the \(x\) direction, with the lower region in tension. At the transition between the compressive and tensile regions, the stress becomes zero; this is the neutral axis of the beam. If the material tends to fail in tension, like chalk or glass, it will do so by crack initiation and growth from the lower tensile surface. If the material is strong in tension but weak in compression, it will fail at the top compressive surface; this might be observed in a piece of wood by a compressive buckling of the outer fibers. We seek an expression relating the magnitudes of these axial normal stresses to the shear and bending moment within the beam, analogously to the shear stresses induced in a circular shaft by torsion. In fact, the development of the needed relations follows exactly the same direct approach as that used for torsion: 1. Geometrical statement: We begin by stating that originally transverse planes within the beam remain planar under bending, but rotate through an angle \(\theta\) about points on the neutral axis as shown in Figure 1. For small rotations, this angle is given approximately by the \(x\)-derivative of the beam's vertical deflection function \(v(x)\) (The exact expression for curvature is \[\dfrac = \dfrac<[1 + (dv/dx)^2]^<3/2>>.\] This gives \(\theta \approx dv/dx\) when the squared derivative in the denominator is small compared to 1.): \[u = -y v_\] where the comma indicates differentiation with respect to the indicated variable (\(v_ \equiv dv/dx\)). Here \(y\) is measured positive upward from the neutral axis, whose location within the beam has not yet been determined. 2. Kinematic equation: The \(x\)-direction normal strain \(\epsilon_x\) is then the gradient of the displacement: \[\epsilon_x = \dfrac = -yv_\] Note that the strains are zero at the neutral axis where \(y = 0\), negative (compressive) above the axis, and positive (tensile) below. They increase in magnitude linearly with \(y\), much as the shear strains increased linearly with \(r\) in a torsionally loaded circular shaft. The quantity \(v_ \equiv d^2v/dx^2\) is the spatial rate of change of the slope of the beam deflection curve, the "slope of the slope." This is called the curvature of the beam. 3. Constitutive equation: The stresses are obtained directly from Hooke’s law as \[\sigma_x = E\epsilon_x = -y Ev_\] This restricts the applicability of this derivation to linear elastic materials. Hence the axial normal stress, like the strain, increases linearly from zero at the neutral axis to a maximum at the outer surfaces of the beam. 4. Equilibrium relations: Since there are no axial (\(x\)-direction) loads applied externally to the beam, the total axial force generated by the normal \(\sigma_x\) stresses (shown in Figure 2) must be zero. This can be expressed as \(\sum F_x = 0 = \int_A \sigma_x dA = \int_A -y Ev_ dA\) which requires that \(\int_A y dA = 0\) The distance \(\bar\) from the neutral axis to the centroid of the cross-sectional area is \(\bar = \dfrac<\int_A y dA><\int_A dA>\) Hence \(\bar = 0\), i.e. the neutral axis is coincident with the centroid of the beam cross-sectional area. This result is obvious on reflection, since the stresses increase at the same linear rate, above the axis in compression and below the axis in tension. Only if the axis is exactly at the centroidal position will these stresses balance to give zero net horizontal force and keep the beam in horizontal equilibrium. The normal stresses in compression and tension are balanced to give a zero net horizontal force, but they also produce a net clockwise moment. This moment must equal the value of \(M(x)\) at that value of \(x\), as seen by taking a moment balance around point \(O\): \(\sum M_O = 0 = M + \int_A \sigma_x \cdot y dA\) \[M = \int_A (y Ev_) \cdot y dA = Ev_ \int_A y^2 dA\] The quantity \(\int y^2 dA\) is the rectangular moment of inertia with respect to the centroidal axis, denoted \(I\). For a rectangular cross section of height \(h\) and width \(b\) as shown in Figure 3 this is: \[I = \int_^ y^2 b dy = \dfrac\] Solving Equation 4.2.4 for \(v_\), the beam curvature is \[v_ = \dfrac\] 5. An explicit formula for the stress can be obtained by using this in Equation 4.2.3: \[\sigma_x = -y E \dfrac = \dfrac\] The final expression for stress, Equation 4.2.7, is similar to \(\tau_ = Tr/J\) for twisted circular shafts: the stress varies linearly from zero at the neutral axis to a maximum at the outer surface, it varies inversely with the moment of inertia of the cross section, and it is independent of the material’s properties. Just as a designer will favor annular drive shafts to maximize the polar moment of inertia \(J\), beams are often made with wide flanges at the upper and lower surfaces to increase \(I\).
Example \(\PageIndex<1>\): Cantilevered T-beam Consider a cantilevered T-beam with dimensions as shown in Figure 4, carrying a uniform loading of \(w N/m\). The maximum bending moment occurs at the wall, and is easily found to be \(M_ <\max>= (wL)(L/2)\). The stress is then given by Equation 4.2.7, which requires that we know the location of the neutral axis (since \(y\) and \(I\) are measured from there). The distance \(y\) from the bottom of the beam to the centroidal neutral axis can be found using the "composite area theorem" (see Exercise \(\PageIndex<1>\)). This theorem states that the distance from an arbitrary axis to the centroid of an area made up of several subareas is the sum of the subareas times the distance to their individual centroids, divided by the sum of the subareas( i.e. the total area): \(\bar = \dfrac<\sum_i A_i \bar_i>\) For our example, this is \(\bar = \dfrac<(d/2)(cd) + (d + b/2)(ab)>\) The moments of inertia of the individual parts of the compound area with respect to their own centroids are just \(ab^3/12\) and \(cd^3/12\). These moments can be referenced to the horizontal axis through the centroid of the compound area using the "parallel axis theorem" (see Exercise \(\PageIndex\)). This theorem states that the moment of inertia \(I_\) of an area \(A\), relative to any arbitrary axis \(z'\) parallel to an axis through the centroid but a distance \(d\) from it, is the moment of inertia relative to the centroidal axis \(I_z\) plus the product of the area \(A\) and the square of the distance \(d\): \(I_ = I_z + A d^2\) For our example, this is The moment of inertia of the entire compound area, relative to its centroid, is then the sum of these two contributions: \(I = I^ <(1)>+ I^\) The maximum stress is then given by Equation 4.2.7 using this value of \(I\) and \(y = \bar/2\) (the distance from the neutral axis to the outer fibers), along with the maximum bending moment \(M_<\max>\). The result of these substitutions is \(\sigma_x = \dfrac<(3d^2c + 6abd + 3ab^2)wL^2>\) In practice, each step would likely be reduced to a numerical value rather than working toward an algebraic solution.1>
In pure bending (only bending moments applied, no transverse or longitudinal forces), the only stress is \(\sigma_x\) as given by Equation 4.2.7. All other stresses are zero (\(\sigma_y = \sigma_z = \tau_
Transverse loads bend beams by inducing axial tensile and compressive normal strains in the beam's \(x\)-direction, as discussed above. In addition, they cause shear effects that tend to slide vertical planes tangentially to one another as depicted in Figure 8, much like sliding playing cards past one another. The stresses \(\tau_
Figure 8: Shearing displacements in beam bending.
Figure 9: Shear and bending moment in a differential length of beam. where \(b\) is the width of the beam at \(y, \xi\) is a dummy height variable ranging from \(y\) to the outer surface of the beam, and \(A'\) is the cross-sectional area between the plane at \(y\) and the outer surface. Using \(dM = V\ dx\) from Equation 4.2.8 of Module 12, this becomes \[\tau_
Example \(\PageIndex<2>\) Since the normal stress is maximum where the horizontal shear stress is zero (at the outer fibers), and the shear stress is maximum where the normal stress is zero (at the neutral axis), it is often possible to consider them one at a time. However, the juncture of the web and the flange in I and T beams is often a location of special interest, since here both stresses can take on substantial values. Consider the T beam seen previously in Example \(\PageIndex\), and examine the location at point \(A\) shown in Figure 11, in the web immediately below the flange. Here the width \(b\) in Equation 4.2.12 is the dimension labeled \(c\); since the beam is thin here the shear stress \(\tau_\) will tend to be large, but it will drop dramatically in the flange as the width jumps to the larger value a. The normal stress at point \(A\) is computed from \(\sigma_x = My/I\), using \(y = d − y\). This value will be almost as large as the outer-fiber stress if the flange thickness b is small compared with the web height \(d\). The Mohr’s circle for the stress state at point \(A\) would then have appreciable contributions from both \(\sigma_x\) and \(\tau_\), and can result in a principal stress larger than at either the outer fibers or the neutral axis. This problem provides a good review of the governing relations for normal and shear stresses in beams, and is also a natural application for symbolic-manipulation computer methods. Using Maple software, we might begin by computing the location of the centroidal axis: Here the ">" symbol is the Maple prompt, and the ";" is needed by Maple to end the command. The maximum shear force and bending moment (present at the wall) are defined in terms of the distributed load and the beam length as2>
> ybar := ((d/2)*c*d) + ( (d+(b/2) )*a*b )/( c*d + a*b );
For plotting purposes, it will be convenient to have a height variable Y measured from the bottom of the section. The relations for normal stress, shear stress, and the first principal stress are functions of Y; these are defined using the Maple “procedure” command:
> V := w*L; > M := -(w*L)*(L/2);
For plotting purposes, it will be convenient to have a height variable Y measured from the bottom of the section. The relations for normal stress, shear stress, and the first principal stress are functions of Y; these are defined using the Maple “procedure” command:
> sigx := proc (Y) -M*(Y-ybar)/Iz end; > tauxy := proc (Y) V*Q(Y)/(Iz*B(Y) ) end; > sigp1 := proc (Y) (sigx(Y)/2) + sqrt( (sigx(Y)/2)^2 + (tauxy(Y))^2 ) end;The moment of inertia Iz is computed as
> I1 := (a*b^3)/12 + a*b* (d+(b/2)-ybar)^2; > I2 := (c*d^3)/12 + c*d* ((d/2)-ybar)^2; > Iz := I1+I2;
The beam width B is defined to take the appropriate value depending on whether the variable Y is in the web or the flange:
> B:= proc (Y) if YThe command "fi" ("if" spelled backwards) is used to end an if-then loop. The function Q(Y) is defined for the web and the flange separately:
> Q:= proc (Y) if Y int( (yy-ybar)*c,yy=Y..d) + int( (yy-ybar)*a,yy=d..(d+b) ) > else > int( (yy-ybar)*a,yy=Y.. (d+b) ) > fi end;Here "int" is the Maple command for integration, and yy is used as the dummy height variable. The numerical values of the various parameters are defined as
> a:=3: b:=1/4: c:=1/4: d:=3-b: L:=8: w:=100:Finally, the stresses can be graphed using the Maple plot command
> plot(,Y=0..3,sigx=-500..2500);The resulting plot is shown in Figure 12.Example \(\PageIndex<3>\) In the previous example, we were interested in the variation of stress as a function of height in a beam of irregular cross section. Another common design or analysis problem is that of the variation of stress not only as a function of height but also of distance along the span dimension of the beam. The shear and bending moments \(V(x)\) and \(M(x)\) vary along this dimension, and so naturally do the stresses \(\sigma_x (x,y)\) and \(\tau_ (x,y)\) that depend on them according to Equation 4.2.7 and 4.2.12. Consider a short beam of rectangular cross section subjected to four-point loading as seen in Figure 13. The loading, shear, and bending moment functions are: The shear and normal stresses can be determined as functions of \(x\) and \(y\) directly from these functions, as well as such parameters as the principal stress. Since \(\sigma_y\) is zero everywhere, the principal stress is \(\sigma_ = \dfrac + \sqrt<(\dfrac)^2 + \tau_^2>\) One way to visualize the x-y variation of \(\sigma_\) is by means of a 3D surface plot, which can be prepared easily by Maple. For the numerical values \(P = 100, a = h = 10, b = 3\), we could use the expressions (Maple responses removed for brevity):3>
> # use Heaviside for singularity functions > readlib(Heaviside); > sfn := proc(x,a,n) (x-a)^n * Heaviside(x-a) end; > # define shear and bending moment functions > V:=(x)-> -P*sfn(x,0,0)+P*sfn(x,a,0)+P*sfn(x,2*a,0)-P*sfn(x,3*a,0); > M:=(x)-> P*sfn(x,0,1)-P*sfn(x,a,1)-P*sfn(x,2*a,1)+P*sfn(x,3*a,1); > # define shear stress function > tau:=V(x)*Q/(Iz*b); > Q:=(b/2)*( (h^2/4) -y^2); > Iz:=b*h^3/12; > # define normal stress function > sig:=M(x)*y/Iz; > # define principal stress > sigp:= (sig/2) + sqrt( (sig/2)^2 + tau^2 ); > # define numerical parameters > P:=100;a:=10;h:=10;b:=3; > # make plot > plot3d(sigp,x=0..3*a,y=-h/2 .. h/2);The resulting plot is shown in Figure 14. The dominance of the parabolic shear stress is evident near the beam ends, since here the shear force is at its maximum value but the bending moment is small (plot the shear and bending moment diagrams to confirm this). In the central part of the beam, where \(a < x < 2a\), the shear force vanishes and the principal stress is governed only by the normal stress \(\sigma_x\), which varies linearly from the beam’s neutral axis. The first principal stress is zero in the compressive lower part of this section, since here the normal stress \(\sigma_x\) is negative and the right edge of the Mohr’s circle must pass through the zero value of the other normal stress \(\sigma_y\). Working through the plot of Figure 14 is a good review of the beam stress formulas.
Figure 14: Variation of principal stress \(\sigma_\) in four-point bending.Exercise \(\PageIndex<1>\) Derive the composite area theorem for determining the centroid of a compound area. \(\bar = \dfrac<\sum_i A_i \bar_i>\)1>
Exercise \(\PageIndex<2>\) (a)-(d) Locate the centroids of the areas shown. 2>
Exercise \(\PageIndex<3>\) Derive the "parallel-axis theorem" for moments of inertia of a plane area: \(I_x = I_ + A \bar^2\) \(I_y = I_ + A \bar^2\) 3>
Exercise \(\PageIndex<4>\) (a)-(d) Determine the moment of inertia relative to the horizontal centroidal axis of the areas shown. 4>
Exercise \(\PageIndex<5>\) Show that the moment of inertia transforms with respect to axis rotations exactly as does the stress: where \(I_x\) and \(I_y\) are the moments of inertia relative to the \(x\) and \(y\) axes respectively and \(I_\) is the product of inertia defined as \(I_ = \int_A xy dA\)5>
Exercise \(\PageIndex<6>\) (a)–(h) Determine the maxiumum normal stress σx in the beams shown here, using the values (as needed) \(L = 25\ in\), \(a = 5 \ in\), \(w = 10\ lb/in\), \(P = 150\ lb\). Assume a rectangular cross-section of width \(b = 1\) in and height \(h = 2\ in\). 6>
Exercise \(\PageIndex<7>\) Justify the statement in ASTM test D790, "Standard Test Methods for Flexural Properties of Unreinforced and Reinforced Plastics and Electrical Insulating Materials," which reads: When a beam of homogeneous, elastic material is tested in flexure as a simple beam supported at two points and loaded at the midpoint, the maximum stress in the outer fibers occurs at midspan. This stress may be calculated for any point on the load-deflection curve by the following equation: \(S = 3PL/2bd^2\) where \(S\) = stress in the outer fibers at midspan, MPa; \(P\) = load at a given point on the load-deflection curve; \(L\) = support span, mm; \(b\) = width of beam tested, mm; and d = depth of beam tested, mm.7>
Exercise \(\PageIndex<8>\) Justify the statement in ASTM test D790, "Standard Test Methods for Flexural Properties of Unreinforced and Reinforced Plastics and Electrical Insulating Materials," which reads: The tangent modulus of elasticity, often called the "modulus of elasticity," is the ratio, within the elastic limit of stress to corresponding strain and shall be expressed in megapascals. It is calculated by drawing a tangent to the steepest initial straight-line portion of the load-deflection curve and using [the expression:] \(E_b = L^3 m/4bd^3\) where \(E_b\) = modulus of elasticity in bending, MPa; \(L\) = support span, mm; \(d\) = depth of beam tested, mm; and \(m\) = slope of the tangent to the initial straight-line portion of the load-deflection curve, \(N/mm\) of deflection.8>
Exercise \(\PageIndex<9>\) A rectangular beam is to be milled from circular stock as shown. What should be the ratio of height to width \((b/h)\) to as to minimize the stresses when the beam is put into bending? 9>
Exercise \(\PageIndex<10>\) (a)-(h) Determine the maxiumum shear τxy in the beams of Exercise \(\PageIndex\), , using the values (as needed) \(L = 25\ in, a = 5\ in, w = 10\ lb/in, P = 150\ lb\). Assume a rectangular cross-section of width \(b = 1\) in and height \(h = 2\ in\).10>
Exercise \(\PageIndex<11>\) Show that the ratio of maximum shearing stress to maximum normal stress in a beam subjected to 3-point bending is \(\dfrac = \dfrac\) Hence the importance of shear stress increases as the beam becomes shorter in comparison with its height. 11>
Exercise \(\PageIndex<12>\) Read the ASTM test D4475, "Standard Test Method for Apparent Horizontal Shear Strength of Pultruded Reinforced Plastic Rods By The Short-Beam Method," and justify the expression given there for the apparent shear strength: \(S = 0.849 P/d^2\) where \(S\) = apparent shear strength, \(N/m^2\), (or psi); \(P\) = breaking load, \(N\), (or lbf); and \(d\) = diameter of specimen, m (or in.).12>
Exercise \(\PageIndex<13>\) For the T beam shown here, with dimensions \(L = 3, a = 0.05, b = 0.005, c = 0.005, d = 0.7\) (all in \(m\)) and a loading distribution of \(w = 5000 N/m\), determine the principal and maximum shearing stress at point \(A\). 13>
Exercise \(\PageIndex<14>\) Determine the maximum normal stress in a cantilevered beam of circular cross section whose radius varies linearly from \(4r_0\) to \(r_0\) in a distance \(L\), loaded with a force \(P\) at the free end. 14>
Exercise \(\PageIndex<15>\) A carbon steel column has a length \(L = 1\ m\) and a circular cross section of diameter \(d = 20\ mm\). Determine the critical buckling load \(P_c\) for the case of (a) both ends pinned, (b) one end cantilevered, (c) both ends pinned but supported laterally at the midpoint. 15>
Exercise \(\PageIndex<16>\) A carbon steel column has a length \(L = 1\ m\) and a circular cross section. Determine the diameter \(d\) at which the column has an equal probablity of buckling or yielding in compression.16>
This page titled 4.2: Stresses in Beams is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by David Roylance (MIT OpenCourseWare) via source content that was edited to the style and standards of the LibreTexts platform.
- Back to top
- 4.1: Shear and Bending Moment Diagrams
- 4.3: Beam Displacements