Shake-Down Analysis Structural Guide: Cyclic Plasticity & Melan (2026)

1. Introduction to Structural Shakedown Under Variable Repeated Loads

When civil and industrial structures undergo variable, cyclic, or moving load patterns, classical plastic limit analysis based on monotonic loading overestimates structural safety. A continuous beam, railway bridge, or offshore platform subjected to fluctuating service actions may fail at loads well beneath the static plastic collapse threshold. Performing a rigorous shake-down analysis structural assessment identifies the precise boundary between long-term safe stabilization and progressive structural degradation.

                    STRUCTURAL RESPONSE UNDER CYCLIC LOADS
    Load Level (P)
        ^
        |               /================== Static Plastic Collapse (Pu)
        |              /
        |             +-------------------- Incremental Collapse / Ratcheting
        |            /
        |           +---------------------- Plastic Shakedown / Alternating Plasticity
        |          /
        |         +------------------------ Elastic Shakedown Limit (Ps)
        |        /
        |       +-------------------------- Initial Yield Limit (Py)
        |      /
        +-----+--------------------------------------------------------> Cycles (N)
              0

When loads exceed the elastic limit but remain below the shakedown threshold, the structure undergoes initial localized plastic deformation during early load cycles. This inelastic straining develops a permanent, self-equilibrating residual stress field. If this residual field allows all subsequent load variations to be resisted entirely elastically, the structure has achieved shakedown (or adaptation).

Understanding the governing theorems behind cyclic limit analysis design protects critical infrastructure against both low-cycle fatigue and progressive geometric instability under arbitrary loading paths. Conducting a methodical shake-down analysis structural evaluation ensures long-term cyclic safety.

2. Inelastic Structural Response Regimes Under Cyclic Actions

2.1 Pure Elastic Domain vs Elastic Shakedown

A structure subjected to fluctuating service loads $P(t)$ bounded within a defined load domain $\Omega$ can exhibit distinct behavioral states:

1. Pure Elastic Response: At every point and at every instant, stress states satisfy $f(\boldsymbol{\sigma}) < 0$. Displacements are entirely reversible.
2. Elastic Shakedown: During initial cycles, plastic strains $\boldsymbol{\epsilon}^p$ develop in high-stress zones. Upon unloading, these plastic deformations lock in a protective residual stress field $\boldsymbol{ ho}$. For all subsequent load cycles within the domain, the combined stress $\boldsymbol{\sigma}_{elastic}(t) + \boldsymbol{ ho}$ strictly satisfies the yield criterion $f \le 0$. Plastic dissipation ceases entirely after finite cycles: $\int_0^\infty \dot{\boldsymbol{\epsilon}}^p \, dt < \infty$.

2.2 Alternating Plasticity and Low-Cycle Fatigue

If the load amplitude is excessively large, no residual stress field can prevent yielding under both positive and negative extremes. The structure enters the alternating plasticity regime. In each cycle, the material undergoes closed plastic strain hysteresis loops:

$$\Delta \boldsymbol{\epsilon}^p = \oint \dot{\boldsymbol{\epsilon}}^p \, dt \approx \mathbf{0}$$

Under cyclic plasticity steel conditions, this cyclical reversal produces micro-crack initiation and causes failure via low-cycle fatigue governed by the Coffin-Manson relationship:

$$\frac{\Delta \epsilon_p}{2} = \epsilon_f’ (2N_f)^c$$

where $\epsilon_f’$ is fatigue ductility coefficient, $2N_f$ is number of reversals, and $c$ is the fatigue ductility exponent.

2.3 Incremental Collapse and Ratcheting Mechanics

When asymmetric cyclic loads act in combination with constant sustained gravity loads, the structure may undergo incremental collapse (progressive ratcheting). In this mode, each cyclic excursion contributes a non-zero residual plastic strain increment:

$$\Delta \boldsymbol{\epsilon}^p = \oint \dot{\boldsymbol{\epsilon}}^p \, dt e \mathbf{0}$$

With each successive cycle, total deflections accumulate monotonically until serviceability limits are violated or mechanism collapse occurs.

3. Fundamental Theorems of Shakedown Analysis

3.1 Melan’s Static (Lower Bound) Shakedown Theorem

Formulated by Ernst Melan, the Melan shakedown theorem provides the classic lower-bound static criterion for elastic shakedown:

Theorem Statement: If there exists a time-independent, self-equilibrating residual stress field $\boldsymbol{ ho}(\mathbf{x})$ such that for all possible load combinations within load domain $\Omega$, the sum of the elastic stress $\boldsymbol{\sigma}^e(\mathbf{x}, t)$ and the residual stress $\boldsymbol{ ho}(\mathbf{x})$ nowhere violates the yield condition:
$$f\left(\boldsymbol{\sigma}^e(\mathbf{x}, t) + \boldsymbol{ ho}(\mathbf{x}) ight) \le 0 \quad \forall \mathbf{x} \in V, \, \forall t \ge 0$$
then the structure will eventually shake down elastically under any arbitrary load history within that domain.

Melan’s theorem guarantees that finding any statically admissible residual stress field provides a safe lower-bound estimate of the true shakedown limit factor $\lambda_s$. In practical engineering, executing a cyclic limit analysis assessment provides this safe bound.

3.2 Koiter’s Kinematic (Upper Bound) Shakedown Theorem

Formulated by W. T. Koiter, the kinematic theorem provides an upper-bound approach:

Theorem Statement: A structure will not shake down if there exists any kinematically admissible cycle of plastic strain rates $\dot{\boldsymbol{\epsilon}}^p(\mathbf{x}, t)$ over cycle duration $T$ such that the total external work done by the elastic stress field exceeds the internal plastic dissipation:
$$\int_0^T \left( \int_V \boldsymbol{\sigma}^e(\mathbf{x}, t) : \dot{\boldsymbol{\epsilon}}^p(\mathbf{x}, t) \, dV ight) dt > \int_0^T \left( \int_V D(\dot{\boldsymbol{\epsilon}}^p) \, dV ight) dt$$

3.3 Self-Equilibrating Residual Stress Fields

In frame systems, the residual stress field reduces to residual bending moments $m_r(x)$ satisfying equilibrium with zero external joint loads. The fundamental objective of cyclic limit analysis engineering is determining the optimal residual moment distribution $m_r(x)$ that maximizes the shakedown load multiplier $\lambda_s$.

4. Mathematical Formulation and Optimization Architecture

4.1 Load Domain Polyhedra and Elastic Stress Envelopes

Let a structural frame be subjected to $K$ independent variable load actions $P_k(t) \in [P_k^{min}, P_k^{max}]$. The envelope of extreme elastic bending moments is:

$$M_{j, max}^e = \max_{P \in \Omega} \left[ \sum_{k=1}^K P_k(t) \mathcal{M}_{j,k}^e ight], \quad M_{j, min}^e = \min_{P \in \Omega} \left[ \sum_{k=1}^K P_k(t) \mathcal{M}_{j,k}^e ight]$$

4.2 Linear Programming Formulation for Shakedown Limit Factor

Using Melan’s static theorem, the shakedown load multiplier $\lambda_s$ is determined by solving a Linear Programming (LP) optimization problem:

$$\text{Maximize } \lambda$$

Subject to:
1. $\lambda M_{j, max}^e + m_{r,j} \le M_{p,j}^+$
2. $\lambda M_{j, min}^e + m_{r,j} \ge -M_{p,j}^-$
3. $\mathbf{C}_{eq} \mathbf{m}_r = \mathbf{0}$
4. $\lambda \left( M_{j, max}^e – M_{j, min}^e ight) \le 2 M_{y,j}$

Engineers implement shake-down analysis structural algorithms to solve these inequalities across continuous frameworks.

5. Comprehensive Step-by-Step Worked Numerical Example

Let us evaluate a two-span continuous steel beam subjected to independent moving concentrated loads.

5.1 Two-Span Continuous Steel Beam Configuration

Consider a symmetrical two-span beam over supports $A$, $B$, and $C$, with each span of length $L = 6\text{ m}$.

                 P1 (0 to P)                 P2 (0 to P)
                      |                           |
                      v                           v
          A +---------v---------+ B +-------------v---------+ C
            ^        Span 1     ^   ^          Span 2       ^
            |<------ L=6m ----->|   |<------- L=6m -------->|

• Span lengths: $L = 6.0\text{ m}$.

• Cross-section: Plastic moment capacity $M_p = 180\text{ kNm}$, yield moment $M_y = 180\text{ kNm}$.

• Loads: Two independent vertical concentrated loads $0 \le P_1(t), P_2(t) \le P$ at midspans $D$ and $E$.

5.2 Maximum and Minimum Elastic Bending Moment Envelopes

From elastic analysis of a two-span continuous beam with $L = 6.0\text{ m}$:

• Unit load $P_1 = 1\text{ kN}$: $M_B = -0.5625\text{ kNm}$, $M_D = +1.21875\text{ kNm}$, $M_E = -0.28125\text{ kNm}$

• Unit load $P_2 = 1\text{ kN}$: $M_B = -0.5625\text{ kNm}$, $M_D = -0.28125\text{ kNm}$, $M_E = +1.21875\text{ kNm}$

Extreme elastic moment envelopes:

• Midspan $D$: $M_{max}^e = +1.21875 P$, $M_{min}^e = -0.28125 P$

• Support $B$: $M_{max}^e = 0$, $M_{min}^e = -1.12500 P$

5.3 Formulation of Residual Moment Distribution

For a single degree of indeterminacy, the residual moment is parameterized by support moment $m_{r,B} = R$. Midspan residual moments are $m_{r,D} = m_{r,E} = 0.5 R$.

5.4 Calculation of Shakedown Factor vs Static Collapse Multiplier

Applying Melan’s static inequalities within this shake-down analysis structural procedure:
1. Midspan $D$: $1.21875 P + 0.5 R \le 180$
2. Support $B$: $-1.12500 P + R \ge -180 \implies R \ge 1.12500 P – 180$

Equating the active bounds:
$$1.21875 P + 0.5(1.12500 P – 180) = 180 \implies 1.78125 P = 270 \implies P_s = 151.58\text{ kN}$$

Comparing with limit states:

• Initial Yield Load: $P_y = \frac{180}{1.21875} = 147.69\text{ kN}$

• Elastic Shakedown Limit: $P_s = 151.58\text{ kN}$

• Static Plastic Collapse Load: $P_u = \frac{6 M_p}{L} = \frac{6 \times 180}{6.0} = 180.00\text{ kN}$

Under variable repeated loading exceeding $P_s = 151.58\text{ kN}$, the beam fails by incremental collapse or alternating plasticity. A thorough shake-down analysis structural review ensures that cyclic loads remain strictly below $151.58\text{ kN}$.

6. Engineering Applications: High-Speed Rail, Bridges, and Pressure Vessels

Shakedown principles are incorporated across engineering sectors:

• High-Speed Railway Bridges (EN 1993-2): Millions of axle load passages require shakedown checks to prevent cumulative track deflections.

• Pavement Geotechnics: Heavy vehicle wheel repetitions cause subgrade rutting, controlled by shakedown limits.

• Pressure Vessels & Piping (ASME Section III): Cyclic thermal gradients combined with internal pressure induce thermal ratcheting (Bree problem).

A comprehensive shake-down analysis structural approach guarantees structural longevity in all these demanding applications.

7. Cyclic Structural Limit Synthesis

Performing a thorough shake-down analysis structural evaluation reveals that safety under monotonic loads does not guarantee survivability under repeated cyclic actions. By calculating residual stress equilibria and bounding the energy dissipated across variable load envelopes, engineers protect vital infrastructure from low-cycle fatigue and incremental collapse.

References & Standards Cited

1. Melan, E. (1936). “Theorie statisch unbestimmter Systeme aus ideal-plastischem Baustoff.” Sitzungsberichte der Akademie der Wissenschaften in Wien, Abt. IIa, 145, 195-218.
2. Koiter, W. T. (1960). “General Theorems for Elastic-Plastic Solids.” Progress in Solid Mechanics, North-Holland, Amsterdam, 1, 165-221.
3. European Committee for Standardization (CEN). (2006). Eurocode 3: Design of Steel Structures – Part 2: Steel Bridges (EN 1993-2). Brussels, Belgium.
4. American Society of Mechanical Engineers (ASME). (2021). ASME Boiler and Pressure Vessel Code, Section III. New York, NY.
5. König, J. A. (1987). Shakedown of Elastic-Plastic Structures. Elsevier, Amsterdam.

Frequently Asked Questions (FAQ)

What is the fundamental physical difference between shakedown and static plastic collapse?

Static collapse occurs when sufficient plastic hinges form simultaneously under a single monotonic load. Shakedown describes long-term stabilization under variable loads, where initial plastic deformations lock in residual stresses keeping all subsequent cycles elastic.

Can a structure fail under cyclic loads even if every load stays below the static collapse capacity?

Yes. If cyclic loads exceed the elastic shakedown limit ($P_s$), the structure can fail through either alternating plasticity (low-cycle fatigue) or incremental collapse (ratcheting).

How does strain hardening affect the shakedown limit?

Strain hardening expands the yield surface, which increases the shakedown limit. Design codes typically neglect hardening to maintain a safe lower-bound estimate.

What is the Bree diagram in cyclic thermo-plasticity?

The Bree diagram maps structural behavior under combined constant mechanical tension and cyclic thermal gradients into Elastic, Shakedown, Alternating Plasticity, and Ratcheting regimes.

Why are linear programming algorithms widely used in shakedown analysis?

Because Melan's static shakedown theorem formulates the search for a self-equilibrating residual stress field as a constrained linear optimization problem.