3.4 Beam Bending Stresses & Deflection Analysis

3.4 Beam Bending Stresses & Deflection Analysis

Beam flexure is among the most pervasive design considerations in mechanical engineering. Transmission shafts, crane girders, vehicle chassis rails, and machine tool beds must maintain structural strength against bending stresses while simultaneously satisfying tight deflection limits to preserve gear mesh alignment and bearing life.

---

1. Pure Bending & The Flexure Formula

Euler-Bernoulli beam theory (classical beam theory) describes the behavior of slender beams subjected to transverse flexural loads under the following core assumptions:

1. Plane Sections Remain Plane: Cross-sections perpendicular to the longitudinal axis remain plane and perpendicular to the deformed neutral axis after bending (no interlaminar shear deformation, valid for span-to-depth ratios $L/h > 10$).
2. Linear Elastic Homogeneous Material: The beam material is isotropic, homogeneous, and obeys Hooke's Law ($\sigma = E \epsilon$).
3. Prismatic Geometry with Symmetry: The beam possesses a vertical longitudinal plane of symmetry, and loads act strictly in this plane.

                    Compression Zone (sigma < 0)
               ----------------------------------
    - - - - - - - - - - - Neutral Axis - - - - - - - - - - (sigma = 0)
               ----------------------------------
                     Tension Zone (sigma > 0)

Kinematic Derivation & Strain Distribution

When an internal positive bending moment $M$ bends a beam segment into a circular arc of curvature $\kappa = 1/\rho$ (where $\rho$ is the radius of curvature of the neutral surface), longitudinal strain $\epsilon_x$ varies linearly with vertical distance $y$ measured from the neutral axis:

$$\epsilon_x(y) = -\frac{y}{\rho} = -y \kappa$$

Applying 1D Hooke's Law yields the normal stress distribution:

$$\sigma_x(y) = E \epsilon_x(y) = -\frac{E y}{\rho}$$

To satisfy axial equilibrium under pure bending, the net internal normal force must vanish:

$$\int_{A} \sigma_x \, dA = -\frac{E}{\rho} \int_{A} y \, dA = 0 \implies \int_{A} y \, dA = 0$$

Fundamental Theorem: The first moment of area $\int y dA$ vanishes identically if and only if the reference axis passes through the centroid of the cross-section. Thus, the neutral axis of an elastic beam coincides with the cross-sectional centroid.

Enforcing moment equilibrium between the internal stress distribution and applied bending moment $M$:

$$M = -\int_{A} y \sigma_x \, dA = \frac{E}{\rho} \int_{A} y^2 \, dA = \frac{E I}{\rho}$$

where $I = \int_A y^2 dA$ is the second moment of area (area moment of inertia) about the centroidal neutral axis. Substituting $\frac{E}{\rho} = \frac{M}{I}$ back into the stress equation yields the Flexure Formula:

$$\sigma(y) = -\frac{M y}{I}$$

Elastic Section Modulus ($S$)

Maximum tensile and compressive stresses occur at the extreme fibers furthest from the neutral axis ($y = c_{\text{top}}$ and $y = c_{\text{bot}}$). Defining $c = \max(c_{\text{top}}, c_{\text{bot}})$:

$$\sigma_{\max} = \frac{M c}{I} = \frac{M}{S}$$

where $S = I / c$ is the elastic section modulus.

---

2. Parallel Axis Theorem & Built-Up Cross-Sections

For complex or asymmetric cross-sections composed of multiple standard geometric sub-shapes (such as T-beams, channels, or angle irons), the area moment of inertia about the composite neutral axis is evaluated using the Parallel Axis Theorem:

$$I_x = \bar{I}_{x'} + A d^2$$

where $\bar{I}_{x'}$ is the centroidal moment of inertia of the sub-part, $A$ is its cross-sectional area, and $d = |\bar{y}_{\text{composite}} - y_i|$ is the perpendicular distance between the composite neutral axis and the sub-part's centroidal axis.

Systematic Procedure for Built-Up Sections

1. Establish a Common Datum: Typically chosen at the bottom-most or top-most surface.
2. Locate the Composite Neutral Axis (Centroid): $\bar{y} = \frac{\sum A_i y_i}{\sum A_i}$
3. Compute Centroidal Moment of Inertia for Each Component: $\bar{I}_i = \frac{b_i h_i^3}{12}$.
4. Apply Parallel Axis Theorem: $I_{\text{NA}} = \sum_{i=1}^{n} \left( \bar{I}_i + A_i d_i^2 \right)$
5. Evaluate Asymmetric Section Moduli: For sections where the neutral axis is not centered, calculate distinct top and bottom moduli: $S_{\text{top}} = I / c_{\text{top}}$ and $S_{\text{bot}} = I / c_{\text{bot}}$.

---

3. Shear & Bending Moment Diagrams: Equilibrium Relations

Constructing shear ($V$) and bending moment ($M$) diagrams relies on the fundamental differential equilibrium relations established for a beam element subjected to a distributed transverse downward load $w(x)$:

                  w(x) dx [Distributed Load]
                     vvvvv
       V      +-----------------+      V + dV
      ^       |                 |        |
      |-----> |   dx Element    | ------>V
      ( M )   |                 |      ( M + dM )
              +-----------------+

$$\frac{dV}{dx} = -w(x) \iff \Delta V = -\int w(x) \, dx$$ $$\frac{dM}{dx} = V(x) \iff \Delta M = \int V(x) \, dx$$

Practical Diagram Rules

• Slope of $V(x)$: Equal to the negative of load intensity $-w(x)$.
• Jump in $V(x)$: Concentrated downward point load $P$ causes an instantaneous downward vertical step in $V(x)$ of magnitude $P$.
• Slope of $M(x)$: Equal to the shear force $V(x)$.
• Change in Moment $\Delta M$: The area under the shear force diagram between two points equals the net change in bending moment between those points.
• Maximum/Minimum Moment Location: Local maxima and minima of $M(x)$ occur where $V(x) = 0$ or where $V(x)$ changes sign across a concentrated load. Always check supports and free ends too: the absolute maximum can sit at a fixed support (as in the propped-cantilever example below) or over an overhang support.
• Jump in $M(x)$: Concentrated clockwise external moment $M_0$ causes an instantaneous upward step in the bending moment diagram.

---

4. Differential Equation of Beam Deflection & Boundary Conditions

In classical Euler-Bernoulli theory, curvature $\kappa$ approximates the second spatial derivative of transverse deflection $v(x)$ for small slopes ($[v'(x)]^2 \ll 1$):

$$\frac{d^2 v}{dx^2} = \frac{M(x)}{E I}$$

Differentiating with respect to $x$ generates the complete hierarchical family of beam differential equations:

$$\text{Deflection: } v(x)$$ $$\text{Slope: } \theta(x) = \frac{dv}{dx}$$ $$\text{Bending Moment: } M(x) = E I \frac{d^2 v}{dx^2}$$ $$\text{Shear Force: } V(x) = \frac{dM}{dx} = E I \frac{d^3 v}{dx^3}$$ $$\text{Distributed Load: } -w(x) = \frac{dV}{dx} = E I \frac{d^4 v}{dx^4}$$

Standard Boundary Conditions

---

5. Deflection Tables & Superposition for Statically Indeterminate Beams

Because the governing differential equation is linear, the Principle of Superposition states that the net deflection and slope caused by combined loadings equal the algebraic sum of deflections and slopes caused by each individual load acting separately.

Classic Beam Deflection Formulas

A. Cantilever Beams (Length $L$, Fixed at $x = 0$)

• Tip Point Load $P$ (Downward at $x = L$): $v_{\max} = \frac{P L^3}{3 E I}, \quad \theta_{\max} = \frac{P L^2}{2 E I}$
• Uniform Distributed Load $w$ (Over entire span $L$): $v_{\max} = \frac{w L^4}{8 E I}, \quad \theta_{\max} = \frac{w L^3}{6 E I}$
• End Couple $M_0$ (At $x = L$): $v_{\max} = \frac{M_0 L^2}{2 E I}, \quad \theta_{\max} = \frac{M_0 L}{E I}$

B. Simply Supported Beams (Span $L$, Pinned at Ends)

• Center Point Load $P$ (At $x = L/2$): $v_{\max} = \frac{P L^3}{48 E I} \quad (\text{at midspan}), \quad \theta_{\text{ends}} = \frac{P L^2}{16 E I}$
• Uniform Distributed Load $w$ (Over entire span $L$): $v_{\max} = \frac{5 w L^4}{384 E I} \quad (\text{at midspan}), \quad \theta_{\text{ends}} = \frac{w L^3}{24 E I}$
• Intermediate Point Load $P$ (At distance $a$ from left, $b = L - a$ from right): $v(x) = \frac{P b x}{6 L E I} (L^2 - b^2 - x^2) \quad (\text{for } 0 \le x \le a)$

Solving Statically Indeterminate Beams by Superposition

For a propped cantilever beam (fixed at $A$, supported by a roller prop at $B$):

1. Identify the redundant reaction (e.g., roller prop reaction $R_B$).
2. Release the prop to create a primary cantilever beam.
3. Write the deflection compatibility condition at support $B$: $v_B = v_{B,\text{loads}} - v_{B,R_B} = 0$
4. Substitute standard cantilever formulas and solve for $R_B$.

---

6. Moment-Area Theorems

The Moment-Area method provides a powerful geometric-analytical approach for computing slopes and deflections directly from the area and first moment of the $M / (EI)$ diagram.

First Moment-Area Theorem (Change in Slope)

The change in slope between any two points $A$ and $B$ on the elastic curve equals the total area under the $M / (EI)$ diagram between those points:

$$\theta_{B/A} = \theta_B - \theta_A = \int_{A}^{B} \frac{M(x)}{E I} \, dx$$

Second Moment-Area Theorem (Tangential Deviation)

The vertical tangential deviation $t_{B/A}$ of point $B$ on the elastic curve relative to the tangent line drawn from point $A$ equals the first moment of the area under the $M / (EI)$ diagram between $A$ and $B$, evaluated about an axis passing through $B$:

$$t_{B/A} = \int_{A}^{B} \frac{M(x)}{E I} x \, dx = \bar{x}_B \left( \text{Area}_{A \to B} \right)$$

where $\bar{x}_B$ is the horizontal distance from point $B$ to the centroid of the $M/(EI)$ area.

---

7. Composite (Bimetallic) Beams & The Transformed Section Method

When a flexural member consists of two or more distinct materials bonded together (e.g., steel reinforcement plates on wood beams, or bimetallic strips in thermal actuators), the cross-section cannot be treated with the standard flexure formula directly because $E$ varies across depth.

   Original Composite Section:            Transformed Equivalent Section:
        +-----------------+                    +---------------------------------+
        | Material 2 (E2) |                    | Material 1 Equivalent (n * b2)  |
        +-----------------+   ======>          +---------------------------------+
        | Material 1 (E1) |                    |         Material 1 (b1)         |
        +-----------------+                    +---------------------------------+

Modular Ratio ($n$)

Select one material as the baseline reference (Material 1 with $E_1$). The modular ratio $n$ scales the second material's stiffness:

$$n = \frac{E_2}{E_1}$$

Transformation Rule

To model the beam as a homogeneous body of Material 1, replace Material 2 with an equivalent width $b^*$ scaled by $n$ while keeping vertical heights and distances from the neutral axis unchanged:

$$b^* = n \cdot b_2$$

Stress Evaluation in Composite Beams

1. Locate the transformed centroid $\bar{y}_{\text{trans}} = \frac{\sum A_{tr,i} y_i}{\sum A_{tr,i}}$.
2. Compute transformed area moment of inertia $I_{\text{trans}}$ using the parallel axis theorem.
3. Evaluate normal bending stress:
   • In the reference material (Material 1): $\sigma_1(y) = -\frac{M y}{I_{\text{trans}}}$
   • In the transformed material (Material 2), multiply by the modular ratio $n$: $\sigma_2(y) = -n \frac{M y}{I_{\text{trans}}}$

Caution: Stresses in the transformed material exhibit an abrupt jump at the material interface proportional to $n$ because strain remains continuous ($\epsilon = -y/\rho$) across the bonded boundary.

8. Worked Step-by-Step Indeterminate Propped Cantilever Example

Problem Statement

A machine support beam consists of a propped cantilever beam fixed against rotation and translation at support $A$ ($x = 0$) and supported by a frictionless roller prop at support $B$ ($x = L = 6.0\text{ m}$).

   Support A                                               Support B
    |//////|                                                  (O)
    |//////|===================================================|===
    |//////|   w = 12 kN/m (Uniform downward load)             |
    |//////|vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv|
           x = 0                 |                            x = 6 m
                                 v P = 30 kN (at midspan x = 3 m)

Given Loading & Section Parameters:

• Span length: $L = 6.0\text{ m}$.
• Uniformly distributed downward load: $w = 12.0\text{ kN/m}$ across the entire span.
• Concentrated downward point load: $P = 30.0\text{ kN}$ at midspan ($x = 3.0\text{ m}$).
• Beam cross-section: Wide-flange beam with flexural rigidity $E I = 30.0 \times 10^6\text{ N}\cdot\text{m}^2 = 30,000\text{ kN}\cdot\text{m}^2$.
• Section depth: $d = 300\text{ mm}$, yielding $c = d/2 = 150\text{ mm} = 0.150\text{ m}$, with $I = 1.50 \times 10^{-4}\text{ m}^4$.

Objective:

1. Use the method of superposition to calculate the reaction force $R_B$ at the roller prop.
2. Determine the reaction force $R_A$ and clamping moment $M_A$ at fixed support $A$.
3. Construct the internal shear and moment expressions, identifying the location of zero shear and maximum positive span moment.
4. Determine the maximum peak bending stress $\sigma_{\max}$ in the beam.

---

Step 1: Degree of Indeterminacy & Redundant Selection

The beam has four reaction components ($R_A, M_A, H_A$ at fixed wall $A$, and vertical reaction $R_B$ at prop $B$). With negligible axial forces ($H_A = 0$), two vertical equilibrium equations exist ($\sum F_y = 0, \sum M = 0$) for three unknown reactions ($R_A, M_A, R_B$). The structure is statically indeterminate to the 1st degree.

We select the roller prop reaction $R_B$ (directed upward) as the redundant reaction.

---

Step 2: Displacement Compatibility Formulation

Releasing support $B$ yields a primary cantilever beam fixed at $A$. Because support $B$ is rigid, the net vertical deflection at $B$ must equal zero:

$v_B = v_{B,w} + v_{B,P} - v_{B,R_B} = 0$

where:

• $v_{B,w}$ is the downward tip deflection due to uniform load $w$.
• $v_{B,P}$ is the downward tip deflection due to midspan point load $P$.
• $v_{B,R_B}$ is the upward tip deflection due to redundant force $R_B$.

---

Step 3: Deflection Component Computations

Using standard cantilever deflection tables:

1. Deflection from Uniform Load $w$ ($12\text{ kN/m}$): $v_{B,w} = \frac{w L^4}{8 E I} = \frac{(12\text{ kN/m})(6.0\text{ m})^4}{8 \times (30,000\text{ kN}\cdot\text{m}^2)} = \frac{12 \times 1296}{240,000} = \frac{15,552}{240,000} = 0.06480\text{ m} = 64.80\text{ mm}$

2. Deflection from Concentrated Load $P$ at Midspan ($a = 3.0\text{ m}$): For a cantilever loaded at distance $a$, the tip deflection is $v_{\text{tip}} = \frac{P a^2}{6 E I}(3L - a)$: $v_{B,P} = \frac{(30\text{ kN})(3.0\text{ m})^2}{6 \times (30,000\text{ kN}\cdot\text{m}^2)} [3(6.0\text{ m}) - 3.0\text{ m}] = \frac{270}{180,000} [15.0] = \frac{4050}{180,000} = 0.02250\text{ m} = 22.50\text{ mm}$

3. Total Downward Deflection without Prop: $v_{B,\text{down}} = v_{B,w} + v_{B,P} = 64.80\text{ mm} + 22.50\text{ mm} = 87.30\text{ mm} = 0.08730\text{ m}$

4. Upward Deflection from Prop Force $R_B$: $v_{B,R_B} = \frac{R_B L^3}{3 E I} = \frac{R_B (6.0\text{ m})^3}{3 \times (30,000\text{ kN}\cdot\text{m}^2)} = \frac{216 R_B}{90,000} = 0.00240 \cdot R_B \quad [\text{m/kN}]$

---

Step 4: Solve for Prop Reaction $R_B$

Enforcing compatibility:

$0.00240 \cdot R_B = 0.08730\text{ m}$

$R_B = \frac{0.08730}{0.00240} = 36.375\text{ kN}$

---

Step 5: Statics for Support Reactions at Wall $A$

Summing vertical forces:

$\sum F_y = R_A + R_B - w L - P = 0$

$R_A = (12\text{ kN/m})(6.0\text{ m}) + 30\text{ kN} - 36.375\text{ kN} = 72.0 + 30.0 - 36.375 = 65.625\text{ kN}$

Summing moments about point $A$ (defining counterclockwise clamping moment $M_A$ positive):

$\sum M_A = M_A - \frac{w L^2}{2} - P \left(\frac{L}{2}\right) + R_B L = 0$

$M_A = \frac{12(6.0)^2}{2} + 30(3.0) - (36.375)(6.0) = 216.0 + 90.0 - 218.25 = 87.75\text{ kN}\cdot\text{m}$

Thus, the wall exerts a counterclockwise restraining moment of $87.75\text{ kN}\cdot\text{m}$ (producing a negative bending moment $M(0) = -87.75\text{ kN}\cdot\text{m}$).

---

Step 6: Shear & Moment Profile Analysis

• Segment $1$ ($0 < x < 3\text{ m}$): $V(x) = R_A - w x = 65.625 - 12 x$ At $x = 0^+$: $V = +65.625\text{ kN}$. At $x = 3.0^-\text{ m}$: $V = 65.625 - 12(3.0) = +29.625\text{ kN}$.
• Jump at $x = 3.0\text{ m}$: $V(3.0^+) = 29.625 - P = 29.625 - 30.0 = -0.375\text{ kN}$ Observation: Shear force transitions abruptly from $+29.625\text{ kN}$ to $-0.375\text{ kN}$ at midspan. Because $V(x)$ passes through zero at $x = 3.0\text{ m}$, the local maximum positive span bending moment occurs precisely at midspan!
• Segment $2$ ($3\text{ m} < x < 6\text{ m}$): $V(x) = 65.625 - 30.0 - 12 x = 35.625 - 12 x$ At $x = 6.0^-\text{ m}$: $V = 35.625 - 12(6.0) = -36.375\text{ kN} = -R_B$, matching the boundary condition exactly.

Now, evaluate internal bending moments at critical sections:

1. At Fixed Wall $A$ ($x = 0$): $M(0) = -M_A = -87.75\text{ kN}\cdot\text{m}$
2. At Midspan ($x = 3.0\text{ m}$): $M(3.0) = -M_A + R_A (3.0) - \frac{w(3.0)^2}{2} = -87.75 + (65.625)(3.0) - \frac{12(9.0)}{2}$ $M(3.0) = -87.75 + 196.875 - 54.0 = +55.125\text{ kN}\cdot\text{m}$
3. At Prop Support $B$ ($x = 6.0\text{ m}$): $M(6.0) = 0\text{ kN}\cdot\text{m}$

---

Step 7: Maximum Peak Bending Stress

Comparing moment magnitudes:

• Fixed support: $|M_{\text{wall}}| = 87.75\text{ kN}\cdot\text{m}$.
• Midspan positive peak: $|M_{\text{midspan}}| = 55.125\text{ kN}\cdot\text{m}$.

The absolute peak bending moment occurs at the fixed clamping wall $A$ ($M_{\max} = 87.75\text{ kN}\cdot\text{m} = 87,750\text{ N}\cdot\text{m}$).

Applying the flexure formula with extreme fiber distance $c = 0.150\text{ m}$ and $I = 1.50 \times 10^{-4}\text{ m}^4$:

$\sigma_{\max} = \frac{M_{\max} c}{I} = \frac{(87,750\text{ N}\cdot\text{m})(0.150\text{ m})}{1.50 \times 10^{-4}\text{ m}^4} = 87,750,000\text{ Pa} = 87.75\text{ MPa}$

At Support $A$: Tensile normal stress of $+87.75\text{ MPa}$ develops at the top flange (resisting hogging curvature), and compressive normal stress of $-87.75\text{ MPa}$ develops at the bottom flange.