Plastic Design Timber Masonry: Non-Linear Limit Analysis Guide (2026)

1. Introduction to Inelastic Behavior in Masonry and Timber Structures

While plastic limit theorems were originally developed for ductile structural steel, modern engineering has successfully adapted these methods for traditional materials like unreinforced masonry and structural timber. Applying inelastic timber-masonry mechanics principles allows structural engineers to evaluate ultimate load capacity, seismic survivability, and historical stability without relying on unrealistic elastic assumptions.

                      MATERIAL CONSTITUTIVE COMPARISON
     Masonry (No Tension)                  Timber (Ductile in Compression)
          Stress (σ)                                Stress (σ)
              ^                                         ^
              |                                  +ft ---+        /-- Brittle Tension
              |                                         |       /
      fc -----+---. (Crushing)                          |      /
              |   |                                     |     /
              |   |                                     +----+----------------> Strain (ε)
              |   |                                     |   /
              +---+-------------> Strain (ε)            |  /
              0                                  -fc ---+-/--------- Plastic Crushing
              (Zero Tension Capacity)                   |/

In unreinforced masonry, non-linearity arises from the material’s near-zero tensile strength, making stability governed by geometry rather than stress limits. Conversely, in structural timber, non-linearity occurs via ductile compressive yielding parallel to the grain and plastic deformation in mechanical fasteners. Utilizing timber nonlinear design alongside masonry limit analysis provides a unified engineering methodology for sustainable, resilient, and historic structures.

By mastering inelastic timber-masonry mechanics techniques, engineers can assess heritage arch bridges, vaulted churches, and multi-storey mass timber frames with high precision. Employing inelastic timber-masonry mechanics ensures accurate structural assessments across historical and modern projects.

2. Masonry Limit Analysis and Heyman’s Classical Formulations

2.1 The Three Fundamental Postulates of Rigid-No-Tension Masonry

In 1966, Professor Jacques Heyman established the modern foundation of masonry limit analysis by formulating three core assumptions:

1. Zero Tensile Strength: Masonry units (stone or brick) and lime mortar have negligible tensile strength ($f_t = 0$). Tensile stresses cause immediate joint opening without resistance.
2. Infinite Compressive Strength: Because real masonry stresses under service loads are typically only 5% to 10% of stone crushing strength $f_c$, the material is idealized as having infinite compressive capacity ($f_c o \infty$).
3. No Sliding Failure: Interlocking friction between masonry blocks is sufficiently high to prevent shear sliding along mortar beds ($\tan \phi \ge \mu_{min}$).

Under these three postulates, masonry behavior becomes purely geometric. Failure occurs when cracks develop through the full thickness of the masonry, forming physical hinge mechanisms.

                    MASONRY HINGE KINEMATICS
        Extrados Hinge (Tension Inside, Compression at Top Edge)
                   | |
                 /     \   <--- Intrados Opening / Crack
                +-------+
                |   *   |  <--- Hinge Point on Extrados
                +-------+

2.2 Thrust Line Analysis and Geometric Safety Factors

The internal line of resistance (thrust line) represents the locus of the resultant compressive force vector across every joint. For stability, the thrust line analysis must satisfy two fundamental conditions:

• Static Admissibility: The thrust line must be in static equilibrium with all external gravity and applied loads.

• Geometric Admissibility: The thrust line must remain entirely contained within the masonry geometry: $y_{intrados}(x) \le y_{thrust}(x) \le y_{extrados}(x)$.

The Geometric Factor of Safety ($FoS_{geom}$) is defined as the ratio by which arch ring thickness $t$ can be scaled down while still accommodating a valid thrust line:

$$FoS_{geom} = \frac{t_{actual}}{t_{minimum\_containment}}$$

2.3 Kinematic Hinging and Mechanism Formation in Arches

A typical single-span masonry arch is statically indeterminate to the third degree. To transform into an unstable collapse mechanism, four plastic hinges must develop in alternating order between the intrados and extrados: $N_{hinges} = D_s + 1 = 3 + 1 = 4\text{ hinges}$.

3. Non-Linear Mechanics and Inelasticity in Timber Structures

3.1 Orthotropic Yield Criteria: Compressive Crushing vs Tensile Splitting

Timber is a natural cellular orthotropic composite exhibiting distinct failure modes depending on grain orientation:

• Parallel to Grain Compression: Under longitudinal compression, wood fibers undergo micro-buckling and cellular crushing, exhibiting ductile plastic yield plateaus ($\epsilon_p \ge 1.5\%$).

• Parallel to Grain Tension: Under longitudinal tension, timber fractures in a sudden, brittle manner governed by linear elastic fracture mechanics.

• Perpendicular to Grain Compression: Yielding occurs by gradual compaction of the wood lumen cells, providing substantial ductility.

Applying plastic design timber masonry formulations captures these distinct directional response regimes.

3.2 Johansen’s Yield Model for Timber Dowel Connections

In structural timber framing, system ductility is achieved primarily through mechanical connections (dowels, bolts, and nails). Johansen’s European Yield Model (EYM) defines three primary plastic failure modes for dowels in single shear:

• Mode I: Pure crushing of the wood matrix along the fastener length.

• Mode II: Wood crushing accompanied by the formation of one plastic bending hinge in the dowel.

• Mode III: Formation of two plastic bending hinges in the dowel, maximizing connection ductility.

The design shear capacity $R_v$ for Mode III under Eurocode 5 is:

$$R_v = 1.15 \sqrt{2 M_{y,d} f_{h,d} d} + \frac{F_{ax,Rk}}{4}$$

3.3 Moment Redistribution in Continuous Glulam and CLT Beams

In continuous multi-span glulam and CLT beams, compressive fiber yielding over internal supports allows plastic moment redistribution of up to 10%–15% under Eurocode 5 provisions. Integrating plastic design timber masonry allows optimized member sizing in mass timber construction.

4. Modeling Methodologies: Macro-Elements vs Equivalent Frame Models

Engineers evaluate complex masonry and timber systems using specialized non-linear strategies:

4.3 Distinct Element Method (DEM) and Block Discretization

For historic dry-joint or lime-mortared stone masonry vaults, continuum models often fail to predict block detachment, sliding shear keys, and stone rotation. The Distinct Element Method (DEM) models the masonry assembly as a discontinuous collection of rigid or deformable polyhedral blocks interacting through frictional, non-linear contact interfaces. Contact normal forces are governed by zero-tension penalty springs with compressive crushing cutoffs, while shear interactions obey the Coulomb friction criterion with dilation angle $\psi$. DEM simulations provide high-fidelity tracking of multi-ring arch separation, spandrel wall delamination, and barrel vault spreading under dynamic seismic shaking or foundation differential settlements.

In heavy mass timber construction, hybrid analytical models couple non-linear ductile fastener springs with anisotropic continuum timber members. When glued laminated timber (glulam) or cross-laminated timber (CLT) shear walls undergo lateral racking, non-linear kinematic models capture anchor bracket hold-down yielding, shear angle bracket deformation, and panel-to-panel spline friction, giving engineers comprehensive insight into total hysteretic energy dissipation.

5. Comprehensive Step-by-Step Worked Numerical Example

Let us execute a complete plastic design timber masonry calculation to find the collapse load of a semicircular masonry arch subjected to an asymmetric point load.

5.1 Geometry and Loading of a Semicircular Masonry Arch

Consider a symmetrical semicircular stone masonry arch ring:

• Mean centerline radius: $R = 4.0\text{ m}$.

• Ring thickness: $t = 0.40\text{ m}$, Arch width: $b = 1.0\text{ m}$.

• Masonry unit self-weight: $\gamma_m = 22\text{ kN/m}^3$.

• Arch dead weight per unit arc length: $w_d = 22 \times 0.40 \times 1.0 = 8.8\text{ kN/m}$.

• Applied live point load: Vertical load $P$ acting at quarter-point ($\theta = 135^\circ$, $x = 1.172\text{ m}$ from left support $A$).

5.2 Identification of the Four-Hinge Collapse Mechanism

Under asymmetric point loading $P$ at $\theta = 135^\circ$, the arch forms four alternating hinges:

• $H_1$: $\theta_1 = 180^\circ$ (Left Intrados)

• $H_2$: $\theta_2 = 135^\circ$ (Load Point Extrados)

• $H_3$: $\theta_3 = 60^\circ$ (Right Intrados)

• $H_4$: $\theta_4 = 0^\circ$ (Right Extrados)

5.3 Virtual Work Formulation for Collapse Load Factor

Let virtual rotation $\delta \theta_1$ occur at hinge $H_1$:
1. Downward displacement under load $P$: $\delta v_P = R \sin(45^\circ) \delta \theta_1 = 2.8284 \delta \theta_1$.
2. Left quadrant self-weight work: $W_1 = 27.65\text{ kN}$, $\delta v_{G1} = 1.250 \delta \theta_1 \implies W_{ext, 1} = +34.56 \delta \theta_1\text{ kNm}$.
3. Right segment self-weight work: $W_2 = 82.94\text{ kN}$, $\delta v_{G2} = -0.880 \delta \theta_1 \implies W_{ext, 2} = -72.99 \delta \theta_1\text{ kNm}$.

Equating total virtual work to zero ($\delta W_{int} = 0$):
$$P(2.8284 \delta \theta_1) + 34.56 \delta \theta_1 – 72.99 \delta \theta_1 = 0 \implies 2.8284 P = 38.43 \implies P_{collapse} = 13.59\text{ kN}$$

5.4 Geometric Thrust Line Verification

At $P = 13.59\text{ kN}$, the thrust line touches the arch boundary at four points without crossing it, proving exact collapse equilibrium. Practicing plastic design timber masonry confirms that geometric hinging governs failure.

6. Design Standards, Eurocode 5 / 6 Provisions, and Preservation Engineering

Modern international codes govern non-linear design:
1. Eurocode 6 (EN 1996-1-1 / EN 1998-3): Allows pushover analysis of masonry walls; limits rocking drift to $0.008 h_{\pier}$.
2. Eurocode 5 (EN 1995-1-1): Requires minimum connector ductility $\mu \ge 4.0$ for plastic fastener redistribution.
3. ISCARSAH Guidelines: Prioritizes geometric limit analysis for historic preservation over destructive modifications.

7. Structural Engineering Synthesis

Applying plastic design timber masonry methodologies bridges the gap between historical craft and rigorous modern mechanics. By recognizing that masonry collapses by geometry while timber yields through localized fastener and fiber ductility, engineers can accurately assess, preserve, and build resilient structures that honor material physics.

References & Standards Cited

1. Heyman, J. (1966). “The Stone Skeleton.” International Journal of Solids and Structures, 2(2), 249-279.
2. Johansen, K. W. (1949). “Theory of Timber Connections.” IABSE, 9, 249-262.
3. European Committee for Standardization (CEN). (2004). Eurocode 5: Design of Timber Structures (EN 1995-1-1). Brussels, Belgium.
4. European Committee for Standardization (CEN). (2005). Eurocode 6: Design of Masonry Structures (EN 1996-1-1). Brussels, Belgium.
5. ISCARSAH. (2003). Recommendations for the Analysis, Conservation and Structural Restoration of Architectural Heritage. ICOMOS.

Frequently Asked Questions (FAQ)

Why does classical elasticity often fail when evaluating historic masonry arches?

Elastic analysis assumes equal tensile and compressive capacities. Because masonry has zero tensile capacity, elastic models predict premature tensile failures that do not reflect true geometric collapse capacity.

What is Heyman's Safe Theorem in masonry mechanics?

Heyman's Safe Theorem states that if a thrust line in static equilibrium with loads can be drawn completely within the masonry geometry, the structure will not collapse.

How do timber fasteners provide system ductility in earthquake-resistant frames?

While wood is brittle in tension, steel dowels yield in bending while crushing wood fibers in compression, providing reliable hysteretic energy dissipation.

What is the difference between intrados and extrados hinges in masonry vaults?

An intrados hinge forms when cracks open on the outside (extrados), forcing thrust against the inner face (intrados). An extrados hinge forms when cracks open inside, pushing thrust to the outer boundary.

Can cross-laminated timber (CLT) panels undergo plastic redistribution?

CLT wood panels behave elastically, but metallic bracket connections (hold-downs and angle brackets) yield plastically to provide ductile rocking mechanisms.