Viscoelasticity Simulator: Creep, Relaxation & Hysteresis in Polymers

simulator intermediate ~10 min
Loading simulation...
τ = 50 s — stress relaxes to 1/e in 50 seconds

With E=100 MPa and η=5000 MPa·s, the Maxwell relaxation time is 50 s. At ε̇=0.1/s (De=5), the response is predominantly elastic with some viscous dissipation.

Formula

σ(t) = σ₀ · exp(-t·E/η)  (Maxwell relaxation)
De = τ · ε̇ = (η/E) · ε̇  (Deborah number)
E*(ω) = E'(ω) + i·E''(ω)  (complex modulus)

Between Solid and Liquid

Polymers occupy a unique middle ground between elastic solids and viscous liquids. Apply a sudden strain to a polymer sample and the stress does not remain constant (as in a perfect solid) nor drops instantly to zero (as in a liquid) — it decays exponentially over a characteristic relaxation time τ. This time-dependent response, called viscoelasticity, governs everything from shoe sole cushioning to asphalt pavement creep.

The Maxwell Model

The simplest viscoelastic model places a spring (elastic element, modulus E) in series with a dashpot (viscous element, viscosity η). Under constant strain, stress relaxes as σ(t) = σ₀·exp(-t/τ) where τ = η/E. The simulation shows this exponential decay in real time and lets you tune both elements to match different polymer behaviors — from fast-relaxing melts (τ ~ ms) to slow-relaxing glasses (τ ~ years).

The Deborah Number

Whether a polymer appears solid or liquid depends on how fast you probe it relative to its relaxation time. The Deborah number De = τ·ε̇ captures this: at high De (fast loading) the material is stiff and elastic; at low De (slow loading) it flows viscously. Silly Putty famously demonstrates this — it bounces like a ball (high De) but flows under its own weight over minutes (low De). The simulation lets you sweep the Deborah number continuously.

Hysteresis and Energy Dissipation

When a viscoelastic material is cyclically loaded, the stress-strain curve forms a hysteresis loop — the area inside represents energy dissipated as heat per cycle. This is quantified by the loss tangent tan(δ) = E''/E', the ratio of loss to storage modulus. High tan(δ) materials make good vibration dampers; low tan(δ) materials make efficient springs. The simulation visualizes both the hysteresis loop and the frequency-dependent modulus spectrum.

FAQ

What is viscoelasticity?

Viscoelasticity describes materials that exhibit both viscous (liquid-like, time-dependent) and elastic (solid-like, reversible) behavior. All polymers are viscoelastic: they store energy like a spring during deformation but also dissipate energy like a dashpot. The balance depends on temperature and loading rate.

What is the Deborah number?

The Deborah number De = τ/t_obs compares the material's relaxation time τ to the observation time. When De >> 1 (fast loading), the material appears elastic; when De << 1 (slow loading), it flows like a liquid. The name references the biblical prophetess Deborah: 'The mountains flowed before the Lord.'

What is the Maxwell model?

The Maxwell model represents a viscoelastic material as a spring (modulus E) in series with a dashpot (viscosity η). Under constant strain, stress relaxes exponentially: σ(t) = σ₀·exp(-t/τ) where τ = η/E. It captures stress relaxation well but predicts instantaneous creep, so real polymers require multi-element models.

How does temperature affect viscoelastic behavior?

Temperature shifts the entire viscoelastic response along the time axis via time-temperature superposition. The WLF equation quantifies this: a 10 K increase near Tg can reduce relaxation time by a factor of 10–100. This principle allows constructing master curves spanning decades of frequency from limited experimental data.

Sources

Embed

<iframe src="https://homo-deus.com/lab/polymer-science/viscoelasticity/embed" width="100%" height="400" frameborder="0"></iframe>
View source on GitHub