Plasticity, hysteresis, and recovery mechanisms in spider silk fibers

Spider silk is a remarkable biomaterial with exceptional stiffness, strength, and toughness stemming from a unique microstructure. While recent studies show that silk fibers exhibit plasticity, hysteresis, and recovery under cyclic loading, the underlying microstructural mechanisms are not yet fully understood. In this work, we propose a mechanism explaining the loading-unloading-relaxation response through microstructural evolution: initial loading distorts intermolecular bonds, resulting in a linear elastic regime. Upon reaching the yield stress, these bonds dissociate and the external load is transferred to the polypeptide chains, which deform entropically to allow large deformations. Unloading is driven by entropic shortening until a traction free state with residual stretch is achieved. Subsequently, the fiber recovers as chains reorganize and bonds reform, locking the microstructure into a new stable equilibrium that increases stiffness in subsequent cycles. Following these mechanisms, we develop a microscopically motivated, energy-based model that captures the macroscopic response of silk fibers under cyclic loading. The response is decoupled into two parallel networks: (1) an elasto-plastic network of inter- and intramolecular bonds governing the initial stiffness and yield stress, and (2) an elastic network of entropic chains that enable large deformations. The model is validated against experimental data from Argiope bruennichi dragline silk. The findings from this work are three-fold: (1) explaining the mechanisms that govern hysteresis and recovery and linking them to microstructural evolution; (2) quantifying the recovery process of the fiber, which restores and enhances mechanical properties; and (3) establishing a predictive foundation for engineering synthetic fibers with customized properties.

💡 Research Summary

This paper investigates the remarkable mechanical behavior of spider dragline silk under cyclic tensile loading, focusing on the observed plasticity, hysteresis, and recovery phenomena. The authors propose that the silk’s response can be understood as the interplay of two parallel microstructural networks: (1) an elasto‑plastic network formed by intermolecular and intramolecular hydrogen bonds, and (2) an elastic network of entropic polypeptide chains.

During the initial loading phase, the hydrogen‑bond network deforms elastically, providing the linear stiffness. When the applied stress reaches a yield value, bonds begin to dissociate, transferring load to the polypeptide chains. These chains then stretch entropically, allowing large deformations akin to rubber elasticity. Upon unloading, the entropic chains contract, but because the bond density has been reduced, the fiber does not return to its original dimensions, leaving a residual stretch (plastic strain).

A subsequent relaxation period enables the broken hydrogen bonds to reform and the chains to reorganize. This “recovery” locks the fiber in a new, slightly elongated configuration, increasing the number of effective cross‑links and the alignment of chains. Consequently, the next loading cycle exhibits higher stiffness and a larger yield stress, a behavior that the authors capture quantitatively.

Mathematically, the total free‑energy density ψ is expressed as the sum of the bond‑network energy ψ_b(λ_b^e) and the chain‑network energy ψ_n(λ_n^e). The bond energy is modeled as ψ_b = (E/2) ln²(λ_b^e), yielding a stress σ_b = E ln(λ_b^e). Plasticity is governed by a yield surface f = σ_b – σ_y(λ_b^p) = 0, where the yield stress σ_y decreases with accumulated plastic stretch λ_b^p, reflecting progressive bond breakage. The chain network follows a classical entropic elasticity model (e.g., neo‑Hookean), providing σ_n. The overall true stress is σ = σ_b + σ_n, and the total stretch λ is decomposed as λ = λ^r λ^l, with λ^r representing the relaxed (post‑recovery) stretch and λ^l the elastic stretch applied during reloading.

The authors calibrate the model against experimental data from Argiope bruennichi dragline silk, including single‑cycle loading to failure and multi‑cycle tests up to ten repetitions. Parameter fitting (elastic modulus E, initial yield stress σ₀, degradation rate k, recovery amplitude η, and relaxation time τ) reproduces: (i) the initial linear modulus and yield point, (ii) the area of hysteresis loops corresponding to dissipated energy, (iii) the progressive increase in stiffness and yield stress across cycles, and (iv) the modest reduction of residual strain during short relaxation intervals.

Key insights from the study are:

1. Plastic deformation originates primarily from the progressive dissociation of hydrogen bonds, not from chain slippage.
2. Hysteresis is a direct manifestation of energy dissipated during bond breakage and reformation.
3. Recovery is driven by bond re‑formation and chain realignment, which effectively “heal” the microstructure and enhance mechanical performance.
4. The parallel‑network framework provides a physically grounded, yet computationally tractable, tool for predicting silk behavior under arbitrary loading histories.

The paper concludes that a microscopically motivated, energy‑based model can faithfully capture the complex cyclic response of spider silk, linking macroscopic observables to underlying molecular events. This understanding opens pathways for designing synthetic fibers that mimic or surpass natural silk by tailoring bond density, chain extensibility, and recovery kinetics to achieve desired stiffness, toughness, and self‑healing capabilities. Future work is suggested to incorporate environmental factors such as humidity and temperature, as well as to extend the model to multiaxial loading and impact scenarios.