Virtual Work Method Plastic Analysis: 4-Step Kinematic Mechanism Guide (2026)

1. Introduction to Virtual Work in Plastic Structural Analysis

The kinematic assessment of indeterminate structural systems requires an efficient, mathematically rigorous procedure to evaluate prospective collapse modes. Applying the virtual work method plastic analysis provides civil and structural engineers with a direct, work-energy approach to determine ultimate collapse load factors without solving tedious simultaneous differential equations.

Originating from variational mechanics, the virtual work method plastic analysis evaluates a structure at the threshold of plastic collapse. When sufficient plastic hinges form to transform a statically indeterminate structure into a kinematic mechanism with one or more degrees of freedom, the system can undergo virtual incremental displacements without changes in internal stress.

By setting the incremental work performed by external loads equal to the internal energy dissipated across all active plastic hinges, engineers can immediately calculate an upper bound collapse load factor. When systematically combined with upper bound search techniques, utilizing the virtual work method plastic analysis ensures that complex multi-bay frames, continuous bridge girders, and grillage systems are analyzed with speed, physical clarity, and absolute precision.

2. The Principle of Virtual Work in the Plastic Regime

The principle of virtual work states that for any structural system in static equilibrium, the total virtual work done by external forces over a set of compatible virtual displacements equals the internal virtual work done by internal stresses over the corresponding virtual strains.

                           +--------------------------------------------------+
                           |       PRINCIPLE OF VIRTUAL WORK AT COLLAPSE      |
                           +------------------------+-------------------------+
                                                    |
                                         W_ext = W_int
                                                    |
                      +-----------------------------+-----------------------------+
                      |                                                           |
                      v                                                           v
       +------------------------------+                            +------------------------------+
       |     EXTERNAL WORK (W_ext)    |                            |    INTERNAL WORK (W_int)     |
       +------------------------------+                            +------------------------------+
       | • Point Loads: $\sum P_i \delta_i$|                       | • Dissipated purely at       |
       | • Distributed Loads:         |                            |   plastic hinges:            |
       |   $\int w(x) \delta(x) dx$   |                            |   $\sum M_{p,j} |\theta_j|$  |
       | • Rigid segments: $\delta=0$ |                            | • Rigid bars: work = 0       |
       +------------------------------+                            +------------------------------+

2.1 Classical Virtual Displacements vs Plastic Mechanisms

In elastic analysis, virtual displacements $\delta u$ are purely hypothetical, infinitesimal mathematical perturbations imposed on an elastic body. In contrast, in the virtual work method plastic analysis, the displacement field corresponds to an actual kinematic mechanism of rigid body links rotating about localized plastic hinges.

Because the material between plastic hinges remains rigid (under the rigid-plastic idealization), no strain energy or internal deformation work accumulates along the member spans:

$$\delta W_{int,\text{ members}} = \int_V \sigma_{ij} \, \delta \var\epsilon_{ij} \, dV = 0$$

Internal work occurs entirely within the discrete plastic hinge zones through finite localized rotations.

2.2 Internal Plastic Work Dissipation in Yielded Hinges

At every active plastic hinge $j$, the bending moment equals the plastic moment capacity $M_{p,j}$. As the rigid segments rotate through a virtual angle $\theta_j$, the internal work dissipated is the scalar product of the plastic moment and the absolute hinge rotation angle:

$$W_{int} = \sum_{j=1}^{m} M_{p,j} \cdot |\theta_j|$$

Because plastic work dissipation represents irreversible thermodynamic plastic dissipation, $W_{int}$ is strictly positive ($W_{int} > 0$). The absolute value $|\theta_j|$ must be used regardless of whether the hinge rotation is clockwise or counterclockwise.

2.3 External Virtual Work Computation for Concentrated and Distributed Loads

External work represents the energy introduced by external loads undergoing virtual deflections:

1. For Concentrated Point Loads ($P_i$):
   $$W_{ext} = \sum_{i=1}^{n} P_i \cdot \delta_i$$
   where $\delta_i$ is the virtual displacement collinear with load $P_i$.

2. For Uniformly Distributed Loads ($w$):
   $$W_{ext} = \int_{0}^{L} w(x) \cdot \delta(x) \, dx = w \cdot A_{\delta}$$
   where $A_{\delta}$ is the total area under the virtual displacement diagram along the loaded span. For a triangular displacement profile of base $L$ and peak virtual deflection $\delta_0$:
   $$W_{ext} = w \cdot \left( \frac{1}{2} L \delta_0 \right)$$

3. The Four-Step Kinematic Analysis Algorithm

Executing the virtual work method plastic analysis follows an organized four-step engineering sequence:

4. Mathematical Formulations and Hinge Kinematics

4.1 Determining the Number of Independent Mechanisms ($I = N – R$)

To ensure all possible failure modes are evaluated during a virtual work method plastic analysis, engineers calculate the number of independent collapse mechanisms using Neal and Symonds’ formula:

$$I = N – R$$

where:

• $N$ = Number of potential plastic hinge critical locations (fixed supports, beam-column knees, under concentrated loads, and span locations of maximum moment).

• $R$ = Degree of static indeterminacy (redundancy) of the structure.

• $I$ = Number of fundamental, independent kinematic mechanisms.

Any complex collapse mechanism can be synthesized by linear combinations of these $I$ independent mechanisms.

       TOTAL MECHANISMS (T) = Independent Mechanisms (I) + Combined Mechanisms (C)

4.2 Kinematic Compatibility and Instant Centers of Rotation

In framed structures with multiple members, the kinematic displacement field must satisfy rigid-body kinematics. For any member $AB$ rotating about a fixed or instantaneous center with virtual rotation $\theta$, the virtual transverse displacement at distance $x$ from the pivot is:

$$\delta(x) = x \cdot \theta$$

At an internal continuous joint connecting two segments with rotations $\theta_1$ and $\theta_2$, the relative hinge rotation $\theta_{rel}$ is:

$$\theta_{rel} = |\theta_1 \pm \theta_2|$$

depending on the relative sense of rotation.

Read the AISC Manual for Structural Steel Calculations and Plastic Analysis

5. Comprehensive Step-by-Step Worked Engineering Examples

We now apply the virtual work method plastic analysis to three foundational structural engineering problems.

5.1 Example 1: Propped Cantilever Under Mid-Span Point Load

Consider a prismatic propped cantilever beam $AB$ of span $L$, fixed at support $A$ and roller-supported at $B$, carrying a vertical point load $P$ at mid-span $C$ ($x = L/2$).

    A                               C                               B
    |===============================*===============================O
    ///                           P \/                            /\
    <------------- L/2 -------------><------------- L/2 ------------>

Step 1: Structural Indeterminacy and Hinge Identification

• Degree of static indeterminacy: $R = 1$ (Fixed base has 2 reactions, roller has 1; $3 – 2 = 1$).

• Potential plastic hinge locations ($N = 2$): Fixed end $A$ and point load location $C$. (Support $B$ is a simple roller and cannot sustain a bending moment, so no hinge forms there).

• Number of independent mechanisms: $I = N – R = 2 – 1 = 1$.

Step 2: Postulate Kinematic Mechanism

Let the beam rotate about fixed end $A$ by virtual angle $\theta$.

• Virtual downward displacement at point load $C$:
  $$\delta_C = \left(\frac{L}{2}\right) \theta$$

• Right segment $CB$ rotates about roller $B$ by angle $\theta_B$:
  $$\delta_C = \left(\frac{L}{2}\right) \theta_B \implies \theta_B = \theta$$

• Hinge rotation at fixed support $A$:
  $$\theta_A = \theta$$

• Relative hinge rotation at load point $C$:
  $$\theta_C = \theta + \theta_B = \theta + \theta = 2\theta$$

Step 3: Formulate Virtual Work Equations

1. External Virtual Work ($W_{ext}$):
   $$W_{ext} = P \cdot \delta_C = P \left( \frac{L}{2} \theta \right) = \frac{P L \theta}{2}$$

2. Internal Plastic Work Dissipation ($W_{int}$):
   $$W_{int} = M_p |\theta_A| + M_p |\theta_C| = M_p (\theta) + M_p (2\theta) = 3 M_p \theta$$

Step 4: Solve for Collapse Load ($P_c$)

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

$$\frac{P_c L \theta}{2} = 3 M_p \theta \implies P_c = \frac{6 M_p}{L}$$

In contrast, the elastic first-yield point load is $P_y = \frac{32 M_y}{11 L} \approx 2.91 \frac{M_y}{L}$. The plastic collapse capacity ($P_c = 6.00 \frac{M_p}{L}$) demonstrates an over $100\%$ capacity increase due to combined shape factor and moment redistribution.

5.2 Example 2: Fixed-Ended Beam Under Uniformly Distributed Load (UDL)

Consider a beam of span $L$ and uniform plastic capacity $M_p$, fixed against rotation at both ends $A$ and $B$, subjected to a uniform load $w$ per unit length.

            w (kN/m)
    ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
    A                                                          B
    |==========================================================|
    ///                                                      ///
    <--------------------------- L ---------------------------->

Step 1: Potential Hinges and Mechanism Geometry

• Degree of static indeterminacy: $R = 2$ (Symmetrical vertical loading).

• Potential hinges form at fixed ends $A$ and $B$, and by symmetry at span centerline $C$ ($N = 3$).

• Independent mechanisms: $I = 3 – 2 = 1$.

Let span centerline $C$ deflect vertically by $\delta_0 = (L/2)\theta$.

• Hinge rotation at support $A$: $\theta_A = \theta$

• Hinge rotation at support $B$: $\theta_B = \theta$

• Central hinge rotation at $C$: $\theta_C = \theta + \theta = 2\theta$

Step 2: Virtual Work Computation

1. External Virtual Work ($W_{ext}$):
   The virtual displacement field is symmetric and triangular. The work done by UDL $w$ is the load multiplied by the area of the displacement triangle:
   $$W_{ext} = w \cdot \left( \frac{1}{2} \cdot L \cdot \delta_0 \right) = w \cdot \left( \frac{1}{2} L \left(\frac{L}{2}\theta\right) \right) = \frac{w L^2 \theta}{4}$$

2. Internal Plastic Work ($W_{int}$):
   $$W_{int} = M_p |\theta_A| + M_p |\theta_C| + M_p |\theta_B| = M_p(\theta + 2\theta + \theta) = 4 M_p \theta$$

Step 3: Solve for Collapse UDL ($w_c$)

$$W_{ext} = W_{int} \implies \frac{w_c L^2 \theta}{4} = 4 M_p \theta$$

$$w_c = \frac{16 M_p}{L^2}$$

Total collapse load on the span is $W_c = w_c L = \frac{16 M_p}{L}$.

5.3 Example 3: Two-Span Continuous Beam with Variable Section Capacity

Consider a continuous beam over three supports $A, B,$ and $C$. Span $AB = L_1 = 6\text{ m}$ with $M_{p1} = 120\text{ kN}\cdot\text{ m}$, and Span $BC = L_2 = 4\text{ m}$ with $M_{p2} = 80\text{ kN}\cdot\text{ m}$. A point load $P_1 = 2P$ acts at mid-span of $AB$, and $P_2 = P$ acts at mid-span of $BC$.

           2P                              P
           |                               |
    A      v       B                       v       C
    O--------------O-----------------------*-------O
    /\             /\                              /\
    <---- 6m -----><-------------- 4m ------------->
       Mp = 120                      Mp = 80

Mechanism Candidate 1: Collapse of Span AB

• Hinges form at mid-span of $AB$ and over central support $B$. (Support $A$ is simple, so $\theta_A$ produces zero work).

• At support $B$, the plastic hinge forms in the weaker adjacent member: $M_{p,B} = \min(120, 80) = 80\text{ kN}\cdot\text{ m}$.

• Hinge at mid-span $AB$: capacity $M_{p,mid} = 120\text{ kN}\cdot\text{ m}$.

• Virtual displacement at $2P$: $\delta_1 = 3\theta$.

• Rotations: $\theta_{mid} = 2\theta$, $\theta_B = \theta$.

$$W_{ext} = (2P) \cdot (3\theta) = 6 P \theta$$

$$W_{int} = M_{p1} (2\theta) + M_{p,B} (\theta) = 120(2\theta) + 80(\theta) = 320 \theta$$

$$6 P_1 \theta = 320 \theta \implies P = \frac{320}{6} = 53.33\text{ kN}$$

Mechanism Candidate 2: Collapse of Span BC

• Hinges form at mid-span of $BC$ and over support $B$.

• $M_{p,mid} = 80\text{ kN}\cdot\text{ m}$, $M_{p,B} = 80\text{ kN}\cdot\text{ m}$.

• Virtual displacement at $P$: $\delta_2 = 2\phi$.

• Rotations: $\theta_{mid} = 2\phi$, $\theta_B = \phi$.

$$W_{ext} = P \cdot (2\phi) = 2 P \phi$$

$$W_{int} = 80(2\phi) + 80(\phi) = 240 \phi$$

$$2 P_2 \phi = 240 \phi \implies P = \frac{240}{2} = 120.0\text{ kN}$$

Governing Collapse Load:

By the upper bound theorem, the governing collapse load is the minimum value:

$$P_c = \min(53.33\text{ kN}, 120.0\text{ kN}) = 53.33\text{ kN}$$

Span $AB$ collapses first at $P = 53.33\text{ kN}$ while span $BC$ remains intact.

6. Advanced Mechanism Combination Techniques

In complex multi-story and multi-bay structural frames, evaluating every individual permutation of plastic hinges manually becomes intractable. The virtual work method plastic analysis overcomes this limitation through the method of combining mechanisms.

1. Identify all independent beam mechanisms, sway mechanisms, and joint mechanisms.
2. Select mechanism combinations that cancel common plastic hinge rotations at heavily loaded joints (since canceling a hinge reduces $W_{int}$ while maintaining or increasing $W_{ext}$).
3. Formulate the combined virtual work equation directly:
   $$W_{ext,\text{ combined}} = \sum W_{ext,k}$$
   $$W_{int,\text{ combined}} = \sum W_{int,k} – 2 \sum M_{p,j} |\theta_{\text{ canceled}}|$$
4. Compute the combined load factor and identify the global minimum.

Consult the Eurocode 3 (EN 1993-1-1) Inelastic Design Provisions

7. Common Forensic Pitfalls and Dynamic Work Corrections

Structural forensic investigations highlight several critical errors when applying the virtual work method plastic analysis:

       M_p
        ^
        |     CORRECT: Sign of Hinge Rotation MUST be positive: W_int = \sum M_p |\theta|
        |     ERROR: Assuming negative work for hogging hinges cancels energy dissipation!
        +----------------------------------------------------------------------------> \theta

1. Incorrect Signs on Internal Work: A frequent error is assigning negative signs to internal hinge work based on hogging/sagging sign conventions. Internal plastic dissipation $M_p |\theta|$ is always positive.
2. Ignoring Movable Plastic Hinges Under UDL: Under distributed loads on propped cantilevers, the maximum sagging moment hinge does not form at mid-span ($0.50L$), but at $x = (\sqrt{2}-1)L \approx 0.414L$. Assuming a mid-span hinge yields an unconservative upper-bound error of approximately $3.0\%$.
3. Overlooking Joint Equilibrium: Failing to verify that non-hinged joints satisfy $|M| \le M_p$ can lead to adopting an invalid collapse mode.

8. Synthesis and Engineering Wrap-Up

The virtual work method plastic analysis stands as one of the most elegant and practical tools in structural engineering mechanics. By transforming complex differential boundary conditions into a simple scalar energy balance between external work and internal plastic dissipation, engineers can rapidly determine the ultimate limit state of ductile structures. When applied with rigorous kinematic consistency, this method guarantees that structural safety assessments rest upon solid physical foundations.

References & Standards Cited

1. AISC (2022). Specification for Structural Steel Buildings (ANSI/AISC 360-22), American Institute of Steel Construction, Chicago, IL.
2. CEN (2005). Eurocode 3: Design of steel structures — Part 1-1: General rules and rules for buildings (EN 1993-1-1), European Committee for Standardization, Brussels.
3. Neal, B. G. (1977). The Plastic Methods of Structural Analysis, 3rd Edition, Chapman and Hall, London.
4. Horne, M. R., & Merchant, W. (1965). The Stability of Frames, Pergamon Press, Oxford.
5. Baker, J. F., Horne, M. R., & Heyman, J. (1956). The Steel Skeleton: Vol 2, Plastic Behaviour and Design, Cambridge University Press.
6. ASCE (2022). Minimum Design Loads and Associated Criteria for Buildings and Other Structures (ASCE/SEI 7-22), American Society of Civil Engineers, Reston, VA.

Frequently Asked Questions (FAQ)

Why is internal plastic work always positive in virtual work equations?

Plastic hinge rotation represents irreversible thermodynamic energy dissipation via plastic yield deformation. Regardless of whether a hinge undergoes clockwise or counter-clockwise rotation, or whether the bending moment is hogging (negative) or sagging (positive), physical work is absorbed by the material. Thus, $W_{int} = sum M_p |theta|$ is always positive.

How is the location of the plastic hinge determined under a distributed load?

Under a distributed load, the exact span hinge forms at the point of zero shear force at collapse, which corresponds to the location of maximum bending moment. In a propped cantilever under UDL, setting $dM/dx = 0$ gives $x = (sqrt{2}-1)L approx 0.414L$ from the propped end, yielding the true minimum collapse load $w_c = (11.656) M_p / L^2$.

What happens if a structural frame possesses more than one degree of freedom at collapse?

If a structure develops more plastic hinges than necessary to form a single-degree-of-freedom mechanism, it forms a **complete collapse mechanism** or an **over-complete mechanism**. The virtual work method remains valid by formulating the virtual displacement field in terms of independent virtual rotation parameters ($theta_1, theta_2$) and minimizing the load factor with respect to each parameter.

Can the virtual work method be applied if columns have smaller $M_p$ than beams?

Yes. At any rigid beam-column knee joint where members of different capacities meet ($M_{p,text{beam}} neq M_{p,text{column}}$), the plastic hinge will form in the weaker member immediately adjacent to the joint. The internal work term must use $M_{p,text{weaker}} |theta|$.

Why are elastic deflections omitted when writing the virtual work equation at collapse?

Under the rigid-plastic idealization of limit analysis, elastic deformations are assumed to be negligible compared to the large plastic mechanism rotations developed at incipient collapse ($delta_{elastic} ll delta_{plastic}$). Neglecting elastic deformations simplifies calculations while maintaining high accuracy for ductile frames.