---
title: Visco-Hyperelastic Kelvin–Maxwell Model
url: https://www.emergentmind.com/topics/visco-hyperelastic-kelvin-maxwell-model
type: topic
---

# Visco-Hyperelastic Kelvin–Maxwell Model

The visco-hyperelastic Kelvin–Maxwell model is a thermodynamically consistent constitutive framework for simulating materials that exhibit both finite-strain hyperelasticity and rate-dependent viscous dissipation, generalizing the classical Kelvin–Maxwell solid to large deformations. The model underpins a broad class of nonlinear viscoelastic solids, fluids, and structured continua, serving as a foundation for advanced computational mechanics, the modeling of polymers, biological tissues, and complex fluids. Multiple kinematic and thermodynamic formalisms exist: multiplicative decompositions, additive splits of energies, structure-tensor dynamics, and explicit operator-splitting approaches, all constrained by physical requirements of objectivity, frame invariance, positive-definiteness, and dissipation [1708.08473, 2101.10439, 1712.04257, 2508.14518, 2508.05043, 2104.12399].

## 1. Constitutive and Thermodynamic Frameworks

The canonical finite-strain Kelvin–Maxwell model is formulated through either a multiplicative decomposition of the deformation gradient (e.g., $F = F_e F_i$ [1708.08473]), evolution of strain-like internal variables [2508.14518, 2101.10439], or structure tensors relaxing towards equilibrium [2104.12399]. The constitutive free energy typically decomposes into isochoric (shear) and volumetric components, with the hyperelastic branch defined by a Helmholtz free-energy density (e.g., Mooney–Rivlin, neo-Hookean, or generalized quadratic forms),
\[
\psi(C_e) = c_1(I_1^e - 3) + c_2(I_2^e - 3)
\]
for the multiplicative model [1708.08473] or
\[
\Psi(C, C^v) = \Psi^\mathrm{eq}(C) + \Psi^\mathrm{neq}(C, C^v)
\]
for the series-spring approach [2101.10439]. Lagrangian, Eulerian, and mixed frameworks are utilized, with balance laws written in symmetrizable-hyperbolic conservative form to guarantee finite-speed signal propagation and short-time existence of smooth solutions [2104.12399].

The evolution of viscoelastic (non-equilibrium) internal variables enforces the second law of thermodynamics, typically via a generalized Clausius–Duhem inequality and dissipation potential:
\[
D = -\frac{\partial \Psi}{\partial C^v} : \frac{1}{2} \dot C^v \ge 0.
\]
The kinetic (flow) law is chosen such that the entropy production remains strictly nonnegative.

## 2. Kinematics and Rheological Interpretation

The Kelvin–Maxwell network, in its finite-strain generalization, can be realized as a (i) spring and dashpot in series with large-deformation kinematics [1708.08473, 2101.10439], (ii) equilibrium spring in series or parallel with a generalized Poynting–Thomson (standard-solid) element [2508.14518, 2508.05043], or (iii) a generalized internal variable formalism for structure tensors [2104.12399, 1712.04257]. The kinematic split may use
- Multiplicative (e.g., $F = F_e F_i$)
- Additive strain-like (e.g., $E = E^e + E^\text{v}$ with generalized Hill-strain $E(C)$ [2508.14518])
- Relaxing conformation tensors (e.g., $C^v$, $A_i$)
  
The Maxwell branch is governed by an evolution law for the inelastic or viscous component ($F_i$, $C^v$, $E^\text{v}$, or $A_i$), typically a first-order flow rule proportional to the driving stress or configurational force:
\[
\dot C^v = \frac{1}{\tau} (C - C^v)
\]
or, in multiplicative representation,
\[
L_p := \dot F_i F_i^{-1} = \frac{1}{\eta} F_e^{-1} S_e F_e^{-T} 
\]
where $\tau$ is a relaxation time and $\eta$ is the effective viscosity [1708.08473, 2101.10439, 2508.14518].

For parallel/serial branches or multiple Voigt elements, internal variable evolution becomes a rank-one perturbed ODE system, efficiently integrable via Sherman–Morrison–Woodbury strategies at $\mathcal{O}(M)$ computational cost, crucial for large numbers of non-equilibrium processes [2508.14518].

## 3. Governing Equations and Constitutive Updates

The state variables include the deformation tensor ($F$ or $C=F^T F$), viscoelastic (internal) variables, density, temperature, and possibly energy flux. The strong form of the model consists of:
- Mass, momentum, and (if non-isothermal) energy and entropy conservation [2104.12399]
- Structure-tensor or conformation-tensor evolution equations with upper-convected derivatives [2104.12399, 1712.04257]
- Constitutive laws for the stresses:
  \[
  S = 2\,\frac{\partial\Psi^\mathrm{eq}}{\partial C} + 2\,\frac{\partial\Psi^\mathrm{neq}}{\partial C}
  \]
  or, in Eulerian terms, Cauchy or Kirchhoff stresses via the Piola transformation.

Implementation relies on efficient time-stepping algorithms:
- Exponential-mapping and iteration-free implicit updates for multiplicative models [1708.08473]
- Explicit, unconditionally stable updates for coupled internal variables via analytical matrix exponentials [2508.14518]
- Heun's predictor–corrector or backward Euler for Lagrangian strain-space models [2508.05043]
- Operator-splitting with finite-volume (Suliciu-type Riemann) solvers for flow problems [1712.04257]
  
These methods preserve frame-invariance, positive-definiteness, symmetry of internal variables, and strict inelastic incompressibility, a critical requirement for large deformations [1708.08473]. Approximation schemes may be unconditionally stable, first- or second-order accurate in time, and iteration-free, with geometric conservation properties exactly satisfied.

## 4. Specialization, Limiting Cases, and Model Variants

The small-strain ($F\approx I+\varepsilon$) limit of the finite-strain Kelvin–Maxwell models reproduces the classical linear viscoelastic solid:
\[
\dot\sigma + \frac{2(c_1+c_2)}{\eta} \sigma = 2(c_1+c_2) \dot\varepsilon
\]
where $\sigma$ is the stress, $c_1$ and $c_2$ are Mooney–Rivlin moduli, and $\eta$ the viscosity [1708.08473, 2101.10439].

Standard-solid (Zener or Poynting–Thomson) generalizations are expressed by additional springs/dashpots in series or parallel, with analytic limits corresponding to Maxwell, Kelvin–Voigt, and pure elastic response [2508.05043]. Model selection is informed by experimental fitting, as in uniaxial extension tests for various polymers.

In non-isothermal extensions, the free energy is augmented by temperature-dependent volumetric terms and the balance laws incorporate a Maxwell–Cattaneo heat flux for finite-speed thermal waves [2104.12399]. Causal properties are preserved, and strictly convex total entropy ensures symmetric–hyperbolic well-posedness.

## 5. Numerical Schemes, Robustness, and Benchmarks

Numerical realization of the visco-hyperelastic Kelvin–Maxwell model employs:
- Operator-splitting, relaxation-based finite-volume schemes for shallow free-surface gravity flows (Saint-Venant–Maxwell), with provable entropy consistency and exact mass/momentum conservation [1712.04257]
- NURBS-based inf-sup stable elements and generalized-alpha schemes for high-fidelity incompressible viscoelasticity with proven energy stability and large-deformation robustness [2101.10439]
- Matrix-exponential integration harnessing the rank-one perturbation structure for multi-process networks, enabling $\mathcal O(M)$ cost per integration step [2508.14518]

Benchmark problems include classical dam break, lid-driven cavity flow, uniaxial extension, and large stretch-rate tests, consistently demonstrating agreement with observed behavior of elastomers and avoidance of the High-Weissenberg Number Problem without breakdown [1712.04257, 2508.05043]. Design parameters are established experimentally (e.g., for VHB-4910 polymers: $\alpha_1=21$, $\alpha_2=20$, $\alpha_3=20$, $\alpha_4=28$, $\alpha_5=4$, $\alpha_6=5$, $\eta=140$ [2508.05043]).

## 6. Hyperbolic Structure, Entropy, and Causality

Properly formulated Kelvin–Maxwell models admit a symmetrizable hyperbolic structure, with strictly convex entropy and all real, distinct characteristic speeds (acoustic, elastic, heat waves), guaranteeing finite-propagation signal speeds and mathematical well-posedness [2104.12399, 1712.04257]. Polyconvex energy functions in the presence of structure tensors ensure stability under large deformations, and dissipation inequalities prescribed at the tensorial evolution level enforce thermodynamic admissibility.

The inclusion of relaxation source terms for stresses and heat-flux, as well as polyconvexity of the free energy in all arguments, ensures short-time existence and uniqueness of solutions even in non-smooth, multi-physics regimes.

## 7. Applications and Extensions

The visco-hyperelastic Kelvin–Maxwell model is central to the modeling of
- Large-deformation polymer mechanics, elastomers, and soft biological tissues across broad strain rates and loading paths
- Complex fluids exhibiting both solid-like and fluid-like responses, including free-surface gravity-driven flows, with Saint-Venant type approximations [1712.04257]
- High-strain-rate and non-isothermal continuum mechanics, underpinning fully causal viscoelastic flow models with energy transport [2104.12399]

Variants extend to multiple relaxation processes (generalized Maxwell), anisotropy via objective structure-tensor evolution, and direct embedding into thermodynamically stable finite element frameworks [2508.14518, 2101.10439].

The robust and efficient numerical algorithms developed—especially non-iterative, explicit time-integration and entropy-consistent finite-volume solvers—yield algorithmic scalability and reliability, facilitating simulations of realistic three-dimensional boundary value problems and supporting quantitative agreement with experimental data [1708.08473, 2508.05043, 2508.14518].

---

**Key references**: [1708.08473], [2101.10439], [1712.04257], [2508.14518], [2508.05043], [2104.12399]

Source: https://www.emergentmind.com/topics/visco-hyperelastic-kelvin-maxwell-model