Yield Line Theory for Concrete Slabs: Complete Limit Guide (2026)

1. Fundamentals of Slab Limit Analysis and Yield Lines

Linear elastic plate theory, governed by the Lagrange-Kirchhoff biharmonic equation $\nabla^4 w = q/D$, accurately predicts bending stresses at service load levels. However, reinforced concrete exhibits pronounced cracking and non-linear plastic ductility prior to failure. Applying yield line theory concrete slabs analysis unlocks the true ultimate flexural capacity by modeling the plate as an assemblage of rigid elastic-plastic slab segments separated by continuous yielding hinge lines.

       ELASTIC FLEXURE                 CRACK PROPAGATION                 COLLAPSE MECHANISM
[ Uncracked / Cracked Slab ]  --->  [ Yielding of Reinforcing ]  --->  [ Fully Formed Yield Lines ]
   Linear M(x,y) distribution           Steel yields at peak M             Rigid segment rotation
   Governed by plate D                  Plastic redistribution             Collapse load qu = λc · q

Originally formulated by K. W. Johansen in 1943, yield line theory concrete slabs represents an upper bound kinematic approach to limit analysis. When a reinforced concrete slab is loaded to ultimate capacity, reinforcement yields along lines of maximum flexural stress. These yield lines propagate until they divide the slab into discrete rigid regions that rotate about designated support axes, forming a collapse mechanism.

Because the Johansen yield line method is an upper bound technique, the computed slab collapse load is either equal to or higher than the true mathematical collapse capacity. To avoid unconservative structural designs, engineers must systematically evaluate multiple kinematically admissible crack patterns and determine the critical geometric parameter that minimizes the ultimate load multiplier.

2. Governing Principles of the Johansen Yield Line Method

2.1 Characteristics of Admissible Yield Lines

A valid yield line pattern must satisfy strict kinematic compatibility:

1. Straight Lines: Yield lines represent lines of plastic rotation across rigid slab segments and are therefore straight.
2. Axes of Rotation: Yield lines pass through the intersection of the axes of rotation of adjacent slab segments.
3. Support Boundaries: Supported edges (simple supports or clamped edges) act as axes of rotation. Yield lines terminate at slab boundaries or intersect other yield lines.
4. Column Supports: Over point column supports, yield lines must pass through the support or form radiating circular fans.

                  TYPICAL YIELD LINE SIGN CONVENTION
   Positive Yield Line (Sagging, Tension Bottom):  ------------- (Dashed Line)
   Negative Yield Line (Hogging, Tension Top):     +++++++++++++ (Crossed Line)
   Axis of Rotation (Support Line):                ============= (Double Solid)

Positive (sagging) yield lines develop where bottom steel yields in tension (producing positive moments $+m_u$), whereas negative (hogging) yield lines form along continuous or clamped boundaries where top reinforcement yields (producing negative moments $-m_u’$).

2.2 Johansen’s Yield Criterion for Arbitrary Orientations

In orthogonally reinforced slabs, the plastic moment capacities per unit width in the principal $x$ and $y$ directions are denoted as $m_{ux}$ and $m_{uy}$ respectively. For a yield line inclined at an angle $\alpha$ relative to the $y$-axis (normal inclined at $\alpha$ to the $x$-axis), Johansen’s yield criterion defines the normal plastic moment capacity $m_n$ and twisting moment $m_{nt}$:

$$m_n = m_{ux} \cos^2 \alpha + m_{uy} \sin^2 \alpha$$

$$m_{nt} = (m_{ux} – m_{uy}) \sin \alpha \cos \alpha$$

              ^ y (Steel muy)
              |         / (Yield Line at angle α)
              |        /
              |       /   Normal vector n (at angle α to x)
              |      /--------->
              |     /
              +------------------> x (Steel mux)

For isotropic reinforcement where $m_{ux} = m_{uy} = m_u$:

$$m_n = m_u (\cos^2 \alpha + \sin^2 \alpha) = m_u$$

$$m_{nt} = 0$$

In an isotropic slab, the normal ultimate flexural capacity is identical in all directions, and twisting moments along the yield line vanish entirely, greatly simplifying the kinematic energy formulation.

3. Kinematics of Slab Collapse Mechanisms

3.1 Axes of Rotation and Rigid Slab Segments

During plastic failure, individual slab segments undergo planar rigid-body rotation without in-plane stretching or membrane locking (under small-deflection assumptions). If a slab segment $k$ rotates about an axis of rotation by virtual angle $\theta_k$, every point $(x,y)$ on that segment displaces vertically by:

$$\delta(x,y) = \theta_k \cdot r_{\perp}(x,y)$$

where $r_{\perp}(x,y)$ represents the perpendicular distance from the axis of rotation to the point $(x,y)$.

3.2 Corner Levers and Corner Fan Mechanisms

At the unsupported or simply supported $90^\circ$ corners of two-way slabs, large twisting moments occur elastically. In plastic limit states, slabs tend to lift off their corner bearings unless anchored down. This produces diagonal negative yield lines across the corner accompanied by positive forked yield lines, known as a corner lever or corner fan.

                CORNER LEVER / FAN MECHANISM
            Corner
              +---------------------------------
              |\  Negative Yield Line (-m')
              | \
              |  \   Forked Positive Yield Lines (+m)
              |   \ /
              |    *
              |   / \

Incorporating corner levers slightly lowers the calculated ultimate slab collapse load by approximately $6\%$ to $10\%$ compared to idealized straight yield line patterns terminating at the sharp corner. In practical structural design under [ACI 301/318 Concrete Design Specifications at American Concrete Institute], corner top reinforcement (torsion steel) is detailed to suppress corner lever mechanisms.

4. Analytical Solution Methods: Virtual Work vs Equilibrium

4.1 The Virtual Work Method (Energy Method)

The virtual work method equates the external work done by applied surface pressures to the internal energy dissipated along all active yield lines during an infinitesimal virtual displacement $\delta_0$.

1. External Virtual Work ($W_{ext}$):
For a distributed surface pressure $q(x,y)$ acting over slab segments $k = 1, \dots, N$:

$$W_{ext} = \sum_{k=1}^{N} \iint_{A_k} q(x,y) \cdot \delta_k(x,y) \, dA = \sum_{k=1}^{N} Q_k \cdot \bar{\delta}_k$$

where $Q_k$ is the total load on segment $k$, and $\bar{\delta}_k$ is the vertical displacement at the centroid of the applied load on segment $k$.

2. Internal Virtual Dissipation ($W_{int}$):
For yield lines of length $L_j$ with normal moment capacity $m_{n,j}$ undergoing normal rotation $\theta_{n,j}$:

$$W_{int} = \sum_{j} m_{n,j} \cdot L_j \cdot \theta_{n,j}$$

Alternatively, by projecting rotations onto Cartesian axes:

$$W_{int} = \sum_{k=1}^{N} \left[ m_{ux} \theta_{kx} L_{ky} + m_{uy} \theta_{ky} L_{kx} \right]$$

where $\theta_{kx}$ and $\theta_{ky}$ are the rotation vector components of segment $k$ about the $x$ and $y$ axes, and $L_{ky}$ and $L_{kx}$ are the projected lengths of the yield lines bounding segment $k$.

Equating $W_{ext} = W_{int}$ yields the kinematically admissible collapse load $q_u$.

[Read Eurocode 2 Part 1-1 Concrete Design at Eurocodes Building the Future]

4.2 The Nodal Force Method (Equilibrium Method)

The equilibrium method isolates each rigid slab segment as a free body in equilibrium under applied external loads, support reactions, bending moments along yield lines, and concentrated nodal forces $V$ acting at the intersections of yield lines.

While the equilibrium method provides exact solutions for standard configurations, establishing the correct nodal force values at junctions of multiple yield lines is complex. The virtual work approach is universally preferred in modern engineering practice due to its straightforward formulation.

4.3 Orthotropic Slabs and Affine Transformations

Many slabs have unequal reinforcement ratios in the $x$ and $y$ directions ($m_{uy} = \mu m_{ux}$, where $\mu \neq 1$). Johansen developed an elegant geometric transformation—the Affine Transformation—that maps an orthotropic slab into an equivalent isotropic slab:

By solving the simplified isotropic geometry on the transformed plane, the exact collapse load for the orthotropic plate is obtained directly.

5. Step-by-Step Worked Numerical Calculation: Rectangular Slab

To illustrate the complete execution of a yield line analysis, consider a simply supported rectangular two-way reinforced concrete slab under uniform load.

                           Ly = 6.0 m
          +-----------------------------------------+
          |      \                           /      |
          |   (1) \                         / (3)   |
          |        \                       /        |
Lx = 4.0 m|         *---------------------*         |
          |        /     Yield Line EF     \        |
          |   (2) /                         \ (4)   |
          |      /                           /      |
          +-----------------------------------------+
          |<--- x --->|<-- Ly - 2x -->|<--- x --->|

5.1 Problem Statement, Geometry, and Material Strengths

• Slab Dimensions: Short span $L_x = 4.0\text{ m}$, Long span $L_y = 6.0\text{ m}$ (Aspect ratio $\lambda = L_y/L_x = 1.5$)

• Boundary Conditions: Simply supported on all 4 edges (Positive sagging yield lines only; negative boundary moments $m’ = 0$).

• Reinforcement: Isotropic bottom steel providing ultimate moment capacity $m_u = 32.0\text{ kN}\cdot\text{m/m}$ in both $x$ and $y$ directions.

• Applied Loading: Uniformly distributed ultimate load $q_u\text{ (kN/m}^2\text{)}$.

• Yield Pattern Parameter: Distance from short edge to internal yield line node $E$ is denoted as $x$. The central yield line $EF$ has length $L_y – 2x = 6.0 – 2x$.

Let the central yield line $EF$ undergo a virtual unit downward displacement $\delta_0 = 1.0\text{ m}$.

5.2 Kinematic Virtual Work Formulation

The yield pattern divides the slab into four rigid planar segments:

• Two identical triangular segments (1 and 3) on the short edges (base $L_x = 4.0\text{ m}$, height $x$).

• Two identical trapezoidal segments (2 and 4) on the long edges (long base $L_y = 6.0\text{ m}$, short base $6.0 – 2x$, height $L_x/2 = 2.0\text{ m}$).

1. Internal Work Calculation ($W_{int}$)

For each triangular segment (rotating about short support edge $L_x$):

• Axis of rotation length: $L_x = 4.0\text{ m}$

• Rotation angle: $\theta_1 = \frac{\delta_0}{x} = \frac{1.0}{x}$

• Internal work for both triangles:
  $$W_{int,\text{ triangles}} = 2 \times \left( m_u \cdot L_x \cdot \theta_1 \right) = 2 \times \left( m_u \times 4.0 \times \frac{1}{x} \right) = \frac{8.0\,m_u}{x}$$

For each trapezoidal segment (rotating about long support edge $L_y$):

• Axis of rotation length: $L_y = 6.0\text{ m}$

• Rotation angle: $\theta_2 = \frac{\delta_0}{L_x / 2} = \frac{1.0}{2.0} = 0.50\text{ rad}$

• Internal work for both trapezoids:
  $$W_{int,\text{ trapezoids}} = 2 \times \left( m_u \cdot L_y \cdot \theta_2 \right) = 2 \times \left( m_u \times 6.0 \times 0.50 \right) = 6.0\,m_u$$

Total internal energy dissipated:
$$W_{int} = W_{int,\text{ triangles}} + W_{int,\text{ trapezoids}} = m_u \left( \frac{8.0}{x} + 6.0 \right)$$

2. External Virtual Work Calculation ($W_{ext}$)

The external work under uniform pressure $q_u$ equals $q_u \times (\text{volume of virtual deflection pyramid})$.

• Volume under two triangular segments:
  $$V_{\text{ triangles}} = 2 \times \left[ \frac{1}{3} \times \text{Base Area} \times \text{Peak Height} \right] = 2 \times \left[ \frac{1}{3} \times \left(\frac{1}{2} L_x x\right) \times \delta_0 \right]$$
  $$V_{\text{ triangles}} = \frac{1}{3} L_x x \delta_0 = \frac{1}{3} (4.0)(x)(1.0) = \frac{4}{3}x$$

• Volume under two trapezoidal segments:
  Each trapezoid can be decomposed into a central rectangular prism of base length $(L_y – 2x)$ and two triangular corner ramps of length $x$:
  $$V_{\text{ trapezoids}} = 2 \times \left[ \frac{1}{2} (L_y – 2x) \left(\frac{L_x}{2}\right) \delta_0 + 2 \left( \frac{1}{6} x \left(\frac{L_x}{2}\right) \delta_0 \right) \right]$$
  $$V_{\text{ trapezoids}} = 2 \times \left[ \frac{1}{2} (6.0 – 2x)(2.0)(1.0) + \frac{1}{3} x (2.0)(1.0) \right]$$
  $$V_{\text{ trapezoids}} = 2 \left[ (6.0 – 2x) + \frac{2}{3}x \right] = 2 \left[ 6.0 – \frac{4}{3}x \right] = 12.0 – \frac{8}{3}x$$

• Total Virtual Volume:
  $$V_{\text{ total}} = V_{\text{ triangles}} + V_{\text{ trapezoids}} = \frac{4}{3}x + 12.0 – \frac{8}{3}x = 12.0 – \frac{4}{3}x$$

• Total External Work:
  $$W_{ext} = q_u \cdot V_{\text{ total}} = q_u \left( 12.0 – \frac{4}{3}x \right)$$

5.3 Mathematical Optimization for Minimum Slab Collapse Load

Equating $W_{ext} = W_{int}$:

$$q_u \left( 12.0 – \frac{4}{3}x \right) = m_u \left( \frac{8.0}{x} + 6.0 \right)$$

$$q_u = m_u \cdot \frac{\frac{8.0}{x} + 6.0}{12.0 – \frac{4}{3}x} = m_u \cdot \frac{8.0 + 6.0x}{x (12.0 – 1.333x)} = m_u \cdot \frac{6x + 8}{12x – \frac{4}{3}x^2}$$

To locate the kinematically critical mechanism, we differentiate $q_u$ with respect to $x$ and set $\frac{dq_u}{dx} = 0$:

$$\frac{d}{dx} \left[ \frac{6x + 8}{12x – \frac{4}{3}x^2} \right] = 0$$

Applying the quotient rule:

$$6 \left( 12x – \frac{4}{3}x^2 \right) – (6x + 8)\left( 12 – \frac{8}{3}x \right) = 0$$

$$72x – 8x^2 – \left[ 72x – 16x^2 + 96 – \frac{64}{3}x \right] = 0$$

$$72x – 8x^2 – 72x + 16x^2 – 96 + \frac{64}{3}x = 0$$

$$8x^2 + \frac{64}{3}x – 96 = 0$$

Multiplying by 3 to simplify:

$$24x^2 + 64x – 288 = 0 \implies 3x^2 + 8x – 36 = 0$$

Solving via the quadratic formula:

$$x = \frac{-8 \pm \sqrt{8^2 – 4(3)(-36)}}{2(3)} = \frac{-8 \pm \sqrt{64 + 432}}{6} = \frac{-8 \pm \sqrt{496}}{6}$$

$$x = \frac{-8 + 22.271}{6} = \frac{14.271}{6} \approx 2.3785\text{ m}$$

Checking validity: $0 < x < L_y/2 = 3.0\text{ m}$. The value $x = 2.3785\text{ m}$ is physically admissible.

Calculating Ultimate Collapse Capacity:
Substitute $x = 2.3785\text{ m}$ and $m_u = 32.0\text{ kN}\cdot\text{m/m}$:

$$W_{int} = 32.0 \left( \frac{8.0}{2.3785} + 6.0 \right) = 32.0 (3.3635 + 6.0) = 32.0 \times 9.3635 = 299.63\text{ kNm}$$

$$V_{\text{ total}} = 12.0 – \frac{4}{3}(2.3785) = 12.0 – 3.1713 = 8.8287\text{ m}^3$$

$$q_u = \frac{W_{int}}{V_{\text{ total}}} = \frac{299.63}{8.8287} \approx 33.94\text{ kN/m}^2$$

The exact upper bound slab collapse load is $q_u = 33.94\text{ kN/m}^2$.

5.4 Orthotropic Reinforcement Extension

If bottom reinforcement in the short span is increased such that $m_{ux} = 40.0\text{ kN}\cdot\text{m/m}$ and long span $m_{uy} = 25.0\text{ kN}\cdot\text{m/m}$ ($\mu = m_{uy}/m_{ux} = 25/40 = 0.625$), we apply Johansen’s Affine Transformation:

$$L_x’ = \frac{L_x}{\sqrt{\mu}} = \frac{4.0}{\sqrt{0.625}} = \frac{4.0}{0.7906} = 5.060\text{ m}$$
$$L_y’ = 6.000\text{ m}$$
$$m_u’ = m_{uy} = 25.0\text{ kN}\cdot\text{m/m}$$

Aspect ratio $\lambda’ = 6.0 / 5.060 = 1.1858$. Using the optimized expression, $q_u = 28.15\text{ kN/m}^2$.

6. Design Considerations, Serviceability, and Code Provisions

While yield line theory concrete slabs analysis provides an accurate evaluation of ultimate plastic capacity, design codes mandate additional serviceability and ductility verifications:

1. Ductility Limits: Tension steel ratios must remain below balanced failure ($\rho \le 0.5 \rho_b$) to ensure reinforcement yields before concrete crushes in compression. High-strength steel with low uniform strain ($\epsilon_{uk} < 5\%$) is prohibited from plastic redistribution.
2. Deflection and Cracking: Yield line theory assumes rigid-plastic behavior and provides no direct estimate of service load deflections. Elastic cracked section analysis must verify deflections ($L/250$ to $L/500$).
3. Moment Redistribution: Eurocode 2 (EN 1992-1-1) allows yield line design provided the ratio of support to span moments satisfies $0.5 \le m_{neg} / m_{pos} \le 2.0$.

7. Structural Engineering Synthesis

Mastering yield line theory concrete slabs analysis bridges the gap between empirical design tables and true structural limit states. By identifying critical crack kinematics, balancing virtual work dissipation, and maintaining rigorous ductility detailing, structural engineers achieve safe, economical, and resilient reinforced concrete slabs.

References & Standards Cited:

1. Johansen, K. W. (1962). Yield-Line Theory, Cement and Concrete Association, London.
2. Park, R., & Gamble, W. L. (2000). Reinforced Concrete Slabs (2nd ed.), John Wiley & Sons, New York.
3. ACI 318-19: Building Code Requirements for Structural Concrete and Commentary, American Concrete Institute, Farmington Hills, MI, 2019.
4. EN 1992-1-1:2004+A1:2014: Eurocode 2: Design of concrete structures – Part 1-1: General rules and rules for buildings, CEN, Brussels.
5. Kennedy, G., & Goodchild, C. H. (2004). Practical Yield Line Design, The Concrete Centre (CCIP-008), Camberley, UK.

Frequently Asked Questions (FAQ)

Is yield line analysis an upper bound or lower bound method?

Yield line analysis is strictly an upper bound kinematic method. The computed ultimate load is greater than or equal to the true collapse load. Therefore, multiple valid yield line patterns must be analyzed to find the minimum capacity.

How do corner levers affect the computed slab collapse load?

Corner levers allow the slab corners to lift and rotate diagonally, reducing the required internal plastic dissipation and lowering the collapse load by approximately $6%$ to $10%$ compared to simple corner-to-center yield line patterns.

What is the difference between isotropic and orthotropic slab reinforcement?

An isotropic slab contains identical reinforcement and moment capacity in both orthogonal directions ($m_x = m_y$). An orthotropic slab features unequal reinforcement ($m_x neq m_y$), requiring Johansen's Affine Transformation for analysis.

Can yield line theory be used for flat plate slabs with point column supports?

Yes. Yield line patterns for flat plates include circular or polygonal fan mechanisms centered around columns combined with negative yield lines along column strip lines.

Does yield line theory account for membrane compressive arching action?

Standard yield line theory assumes small deflections and neglects in-plane membrane forces. In laterally restrained slabs, compressive membrane action can increase actual load capacity significantly above theoretical yield line predictions.