3.1 Axial Loading, Normal Stress & Statically Indeterminate Members

3.1 Axial Loading, Normal Stress & Statically Indeterminate Members

Axial loading represents one of the most fundamental stress states encountered in machine design and structural mechanics. Whether sizing tie rods, hydraulic cylinder piston shafts, preloaded fasteners, or structural struts, mechanical engineers must evaluate both stress states and axial deflections under static, thermal, and geometric constraints.

---

1. Uniform Axial Normal Stress & Saint-Venant's Principle

When a prismatic bar with a constant cross-sectional area $A$ is subjected to collinear, concentric axial forces $P$ acting through the centroid of its cross-section, the internal resultant normal force produces normal stress $\sigma$ across the transverse plane:

$$\sigma = \frac{P}{A}$$

By sign convention:

• Tensile forces produce positive normal stress ($\sigma > 0$), resulting in longitudinal elongation and lateral Poisson contraction.
• Compressive forces produce negative normal stress ($\sigma < 0$), causing longitudinal shortening and lateral expansion.

          P <--- [============ Prismatic Bar ============] ---> P
                      |<- Cross-Section A ->|
                      |   sigma = P / A     |

Saint-Venant's Principle

Equation $\sigma = P/A$ strictly describes an idealized, uniform stress distribution. Near points of concentrated load application, geometric transitions, or boundary fixtures, localized stress fields exhibit pronounced gradients. According to Saint-Venant's Principle, the localized stress disturbances caused by statically equivalent load configurations attenuate rapidly with distance from the boundary.

At a longitudinal distance approximately equal to the largest transverse dimension (such as diameter $d$ or width $w$) from the point of load introduction, the stress field becomes essentially uniform and equals $P/A$. For machine components whose length $L \gg d$, the uniform stress approximation holds accurately over the vast majority of the member's volume.

---

2. Hooke's Law, Axial Deformation & Stiffness

Within the linear elastic regime of isotropic materials, uniaxial stress and engineering strain $\epsilon$ relate directly through Hooke's Law:

$$\sigma = E \epsilon$$

where $E$ is Young's Modulus (modulus of elasticity). Expressing engineering strain as the ratio of total elongation $\delta$ to original gauge length $L$ ($\epsilon = \delta / L$) and substituting into the normal stress relationship yields the foundational axial deformation formula:

$$\delta = \frac{P L}{A E}$$

From a structural and mechanical dynamics perspective, axial members function as linear elastic springs. We define axial stiffness $k$ and axial flexibility (compliance) $f$ as:

$$k = \frac{P}{\delta} = \frac{A E}{L} \quad [\text{N/m} \text{ or } \text{lb/in}]$$ $$f = \frac{1}{k} = \frac{\delta}{P} = \frac{L}{A E} \quad [\text{m/N} \text{ or } \text{in/lb}]$$

Segmented Bars with Intermediate Loadings

For non-prismatic stepped members, or bars subjected to multiple external axial loads along their span, total axial deflection is evaluated by piecewise superposition of segments with constant internal force $P_i$, area $A_i$, and elastic modulus $E_i$:

$$\delta = \sum_{i=1}^{n} \frac{P_i L_i}{A_i E_i}$$

To apply this equation correctly, engineers must apply the method of sections, cutting between each load point and support to determine the true internal axial normal force $P_i$ from local static equilibrium.

---

3. Variable Cross-Sections & Distributed Axial Loads

When cross-sectional area $A(x)$, internal axial force $P(x)$, or modulus of elasticity $E(x)$ vary continuously along the longitudinal coordinate $x$, the elongation of an infinitesimal slice $dx$ is $d\delta = \frac{P(x) dx}{A(x) E(x)}$. Integrating over the total length $L$ yields:

$$\delta = \int_{0}^{L} \frac{P(x)}{A(x) E(x)} \, dx$$

Self-Weight Elongation of a Vertical Bar

Consider a vertical bar of length $L$, uniform area $A$, mass density $\rho$, and specific weight $\gamma = \rho g$ suspended from its top end and hanging freely under gravity. Total member weight is $W = \gamma A L$.

           ///////// Fixed Support
           |       |
           |       |  ^ x (measured from free bottom tip)
           |       |  |
           |   dx  |  |
           |       |  |
           |_______|
             Free Tip (x = 0)

At a distance $x$ measured upward from the free bottom tip, the internal normal force supports only the weight of the material below that section:

$$P(x) = \gamma A x$$

The total elongation under self-weight is:

$$\delta = \int_{0}^{L} \frac{\gamma A x}{A E} \, dx = \frac{\gamma}{E} \left[ \frac{x^2}{2} \right]_{0}^{L} = \frac{\gamma L^2}{2 E} = \frac{W L}{2 A E}$$

Key Principle: The axial elongation of a hanging prismatic bar due entirely to its own weight is precisely one-half (50%) of the elongation that would occur if its entire weight $W$ were concentrated as a point load at the bottom free end.

Linearly Tapered Circular Shaft

For a circular shaft subjected to axial tension $P$, tapering linearly from diameter $d_1$ at $x = 0$ to $d_2$ at $x = L$, the local diameter is $d(x) = d_1 + \left(\frac{d_2 - d_1}{L}\right) x$. The cross-sectional area is $A(x) = \frac{\pi [d(x)]^2}{4}$. Substituting into the integral formulation and integrating results in:

$$\delta = \frac{4 P L}{\pi E d_1 d_2}$$

---

4. Statically Indeterminate Axial Systems

A structural system is statically indeterminate when the number of unknown support reactions and internal forces exceeds the number of independent equations of static equilibrium available from rigid-body mechanics. For 1D axial systems, statics provides only one non-trivial global equilibrium equation: $\sum F_x = 0$.

To resolve indeterminate reactions, engineers must enforce geometric compatibility of deformations alongside material constitutive relations.

The Standard 4-Step Solution Architecture

1. Static Equilibrium: Draw free-body diagrams (FBD) of the entire structure and key subsections. Formulate the global axial equilibrium equation: $\sum F_x = 0$
2. Displacement Compatibility: Formulate geometric relationships dictating the physical displacement constraints of the structure (e.g., total displacement between rigid end walls must equal zero, or relative displacements of parallel bars connected by a rigid pin must remain planar).
3. Force-Displacement (Constitutive) Relations: Express the physical displacement $\delta_i$ of each individual member or segment in terms of its internal normal force $P_i$, material properties, and thermal changes: $\delta_i = \frac{P_i L_i}{A_i E_i} + \alpha_i L_i \Delta T_i$
4. Synthesis & Simultaneous Solution: Substitute the force-displacement expressions into the compatibility equation to eliminate displacements, creating a closed system of linear algebraic equations to solve for unknown support reactions and member stresses.

Compound (Composite) Columns

In parallel structural configurations, such as a concrete core inside a structural steel sleeve or parallel tie rods sharing an axial load $P_{\text{total}}$, both materials deform by the identical longitudinal amount:

$\delta_1 = \delta_2 = \delta$

$\frac{P_1 L}{A_1 E_1} = \frac{P_2 L}{A_2 E_2} \implies P_2 = P_1 \left( \frac{A_2 E_2}{A_1 E_1} \right)$

Substituting into global equilibrium $P_{\text{total}} = P_1 + P_2$ demonstrates that the applied load distributes in direct proportion to the relative axial stiffness ($AE$) of each member:

$$P_1 = P_{\text{total}} \left( \frac{A_1 E_1}{A_1 E_1 + A_2 E_2} \right), \quad P_2 = P_{\text{total}} \left( \frac{A_2 E_2}{A_1 E_1 + A_2 E_2} \right)$$

---

5. Thermal Effects & Restrained Thermal Stress

Temperature variations cause volumetric changes in solid materials. In an unconstrained, homogeneous bar subjected to a uniform temperature change $\Delta T = T_{\text{final}} - T_{\text{initial}}$, the unrestricted thermal deformation $\delta_T$ is:

$$\delta_T = \alpha L \Delta T$$

where $\alpha$ is the linear coefficient of thermal expansion ($1/^\circ\text{C}$ or $1/^\circ\text{F}$). Because no external forces resist this free expansion, the internal stress remains identically zero ($\sigma = 0$).

Fully Constrained Member

If an axial member is rigidly held between two unyielding, immovable walls, thermal expansion is completely restrained. The total deformation must equal zero:

$$\delta_{\text{total}} = \delta_T + \delta_P = \alpha L \Delta T + \frac{P L}{A E} = 0$$

Solving for the internal reaction force $P$:

$$P = -A E \alpha \Delta T$$

The resulting thermal normal stress $\sigma_T$ is:

$$\sigma_T = \frac{P}{A} = -E \alpha \Delta T$$

Critical Takeaway: Thermal normal stress in a fully constrained member depends exclusively on the elastic modulus $E$, thermal expansion coefficient $\alpha$, and temperature change $\Delta T$. It is completely independent of cross-sectional area $A$ and length $L$.

Initial Clearance (Gap) Problems

When a bar is separated from a rigid boundary by an initial clearance gap $\delta_{\text{gap}}$, free expansion occurs until contact is established. Normal stress develops only when $\delta_T > \delta_{\text{gap}}$:

$$\delta_T + \delta_P = \delta_{\text{gap}} \implies \alpha L \Delta T + \frac{P L}{A E} = \delta_{\text{gap}}$$ $$\sigma = -E \left( \alpha \Delta T - \frac{\delta_{\text{gap}}}{L} \right)$$

---

6. Stress Concentrations at Geometric Discontinuities

Sudden alterations in geometry—such as holes, keyways, shoulders, and fillet radii—disrupt the uniform trajectory of internal stress trajectories, causing localized stress amplification. The severity of this peak is defined by the theoretical (geometric) stress concentration factor $K_t$:

$$K_t = \frac{\sigma_{\max}}{\sigma_{\text{nom}}}$$

          P <--- [===== ( O ) =====] ---> P
                      w   d   w
                 <- Net Area = (w - d) * t ->

Net Area vs. Gross Area Definition

In mechanical engineering calculations, nominal stress $\sigma_{\text{nom}}$ must be carefully defined according to the reference standard used in stress concentration charts (e.g., Peterson's Stress Concentration Factors):

• Gross Area: $A_{\text{gross}} = w \cdot t$ (unnotched width $\times$ thickness).
• Net Area: $A_{\text{net}} = (w - d) \cdot t$ (reduced cross-section at the hole centerline).

For a flat plate of width $w$ containing a central through-hole of diameter $d$, standard machine design charts define nominal stress based on the net area:

$$\sigma_{\text{nom}} = \frac{P}{A_{\text{net}}} = \frac{P}{(w - d) t}$$

For an infinitely wide plate ($d/w \to 0$), theoretical elasticity yields $K_t = 3.00$ at the hole boundary. As the hole diameter becomes a significant fraction of plate width, $K_t$ based on net area decreases toward $2.00$.

Ductile vs. Brittle Material Response

• Ductile Materials (Static Loads): Under static tensile loading, localized stress concentrations exceeding yield strength undergo localized plastic redistribution without compromising the gross structural capacity. Consequently, $K_t$ is typically omitted in static failure analysis of ductile steels, except where low-temperature embrittlement or triaxial constraint applies.
• Brittle Materials & Cyclic Fatigue: In brittle materials (cast iron, hardened tool steels) and in all components subjected to dynamic or cyclic fatigue loading, stress concentrations cannot be ignored. Local notches act as primary crack initiation sites, requiring notch sensitivity adjustments ($K_f = 1 + q(K_t - 1)$).

---

7. Worked Step-by-Step Indeterminate Example

Problem Statement

A stepped structural shaft is rigidly fixed between two immovable, unyielding walls at $A$ and $C$. Segment $AB$ is manufactured from structural steel, and segment $BC$ is manufactured from an aluminum alloy.

   Wall A                                                Wall C
   |//////|=====================|===============|//////|
   |//////|     Segment AB      |  Segment BC   |//////|
   |//////|       (Steel)       |  (Aluminum)   |//////|
   |//////|=====================|===============|//////|
          x = 0                B                x = 1000 mm
                                |--> P = 120 kN

Given Component Data:

• Segment $AB$ (Steel): Length $L_1 = 400\text{ mm}$, diameter $d_1 = 50\text{ mm}$, $E_1 = 200\text{ GPa}$, $\alpha_1 = 12 \times 10^{-6}/^\circ\text{C}$.
• Segment $BC$ (Aluminum): Length $L_2 = 600\text{ mm}$, diameter $d_2 = 30\text{ mm}$, $E_2 = 70\text{ GPa}$, $\alpha_2 = 23 \times 10^{-6}/^\circ\text{C}$.
• An external axial force $P = 120\text{ kN}$ is applied at junction $B$, directed horizontally to the right (toward wall $C$).
• Simultaneously, the operating temperature of the entire assembly increases by $\Delta T = 30^\circ\text{C}$.

Objective: Determine the reactions at support walls $A$ and $C$, the displacement of junction $B$, and the peak normal stress in each material segment.

---

Step 1: Geometric & Stiffness Parameters

First, compute the cross-sectional areas and axial flexibilities of each segment:

$A_1 = \frac{\pi d_1^2}{4} = \frac{\pi (50\text{ mm})^2}{4} = 1963.50\text{ mm}^2 = 1.96350 \times 10^{-3}\text{ m}^2$

$A_2 = \frac{\pi d_2^2}{4} = \frac{\pi (30\text{ mm})^2}{4} = 706.86\text{ mm}^2 = 7.06858 \times 10^{-4}\text{ m}^2$

Compute the mechanical flexibility $f = \frac{L}{AE}$ for each segment:

$f_1 = \frac{0.400\text{ m}}{(1.96350 \times 10^{-3}\text{ m}^2)(200 \times 10^9\text{ Pa})} = 1.01859 \times 10^{-9}\text{ m/N}$

$f_2 = \frac{0.600\text{ m}}{(7.06858 \times 10^{-4}\text{ m}^2)(70 \times 10^9\text{ Pa})} = 1.21261 \times 10^{-8}\text{ m/N}$

Compute the unrestrained free thermal expansion of each segment:

$\delta_{T1} = \alpha_1 L_1 \Delta T = (12 \times 10^{-6}/^\circ\text{C})(0.400\text{ m})(30^\circ\text{C}) = 1.440 \times 10^{-4}\text{ m} = 0.1440\text{ mm}$

$\delta_{T2} = \alpha_2 L_2 \Delta T = (23 \times 10^{-6}/^\circ\text{C})(0.600\text{ m})(30^\circ\text{C}) = 4.140 \times 10^{-4}\text{ m} = 0.4140\text{ mm}$

Total unrestrained thermal expansion:

$\delta_T = \delta_{T1} + \delta_{T2} = 0.1440\text{ mm} + 0.4140\text{ mm} = 0.5580\text{ mm} = 5.580 \times 10^{-4}\text{ m}$

---

Step 2: Static Equilibrium

Let $R_A$ be the reaction exerted by the wall at $A$ on the bar, defined positive directed to the left (inducing tension in segment $AB$). Let $R_C$ be the reaction exerted by the wall at $C$ on the bar, defined positive directed to the left (inducing compression in segment $BC$).

Summing horizontal forces across the entire assembly:

$\sum F_x = -R_A + P - R_C = 0 \implies R_A + R_C = P = 120,000\text{ N}$

Using the method of sections, evaluate internal normal forces in each segment:

• Segment $AB$: $N_1 = R_A$ (positive = tension).
• Segment $BC$: $N_2 = R_A - P = -R_C$ (negative = compression).

---

Step 3: Displacement Compatibility

Because walls $A$ and $C$ are rigid and unyielding, the net total elongation of the assembly must equal zero:

$\delta_{\text{total}} = \delta_1 + \delta_2 = 0$

Expressing each segment's deformation as the sum of mechanical and thermal components:

$\delta_1 = N_1 f_1 + \delta_{T1} = R_A f_1 + \delta_{T1}$

$\delta_2 = N_2 f_2 + \delta_{T2} = (R_A - P) f_2 + \delta_{T2}$

Substituting into the compatibility equation:

$R_A f_1 + \delta_{T1} + (R_A - P) f_2 + \delta_{T2} = 0$

$R_A (f_1 + f_2) - P f_2 + \delta_T = 0$

Solving explicitly for $R_A$:

$R_A = \frac{P f_2 - \delta_T}{f_1 + f_2}$

---

Step 4: Numerical Evaluation

Calculate the numerator and denominator:

$P f_2 = (120,000\text{ N})(1.21261 \times 10^{-8}\text{ m/N}) = 1.45513 \times 10^{-3}\text{ m} = 1.45513\text{ mm}$

$P f_2 - \delta_T = 1.45513\text{ mm} - 0.55800\text{ mm} = 0.89713\text{ mm} = 8.9713 \times 10^{-4}\text{ m}$

$f_1 + f_2 = 1.01859 \times 10^{-9} + 1.21261 \times 10^{-8} = 1.31447 \times 10^{-8}\text{ m/N}$

$R_A = \frac{8.9713 \times 10^{-4}\text{ m}}{1.31447 \times 10^{-8}\text{ m/N}} = 68,250.6\text{ N} \approx 68.25\text{ kN}$

From static equilibrium, evaluate $R_C$:

$R_C = P - R_A = 120.00\text{ kN} - 68.25\text{ kN} = 51.75\text{ kN}$

---

Step 5: Internal Forces, Stresses & Junction Displacement

1. Internal Forces:
   • Steel Segment: $N_1 = +68.25\text{ kN}$ (Tension).
   • Aluminum Segment: $N_2 = -R_C = -51.75\text{ kN}$ (Compression).
2. Normal Stresses: $\sigma_{\text{steel}} = \frac{N_1}{A_1} = \frac{+68,250.6\text{ N}}{1963.50\text{ mm}^2} = +34.76\text{ MPa} \quad (\text{Tension})$ $\sigma_{\text{alum}} = \frac{N_2}{A_2} = \frac{-51,749.4\text{ N}}{706.86\text{ mm}^2} = -73.21\text{ MPa} \quad (\text{Compression})$
3. Displacement of Junction $B$ ($u_B$): Since end $A$ is fixed, the displacement of node $B$ equals the elongation of segment 1: $u_B = \delta_1 = N_1 f_1 + \delta_{T1} = (68,250.6)(1.01859 \times 10^{-9}) + 1.440 \times 10^{-4} = 6.952 \times 10^{-5}\text{ m} + 1.440 \times 10^{-4}\text{ m}$ $u_B = +2.135 \times 10^{-4}\text{ m} = +0.2135\text{ mm} \quad (\text{Displaced to the right toward } C)$

Verification: $\delta_2 = N_2 f_2 + \delta_{T2} = (-51,749.4)(1.21261 \times 10^{-8}) + 4.140 \times 10^{-4} = -6.275 \times 10^{-4} + 4.140 \times 10^{-4} = -2.135 \times 10^{-4}\text{ m} = -0.2135\text{ mm}$. Total assembly deformation: $\delta_1 + \delta_2 = +0.2135\text{ mm} - 0.2135\text{ mm} = 0$, confirming exact equilibrium and compatibility satisfaction.