Visco-Hyperelastic Kelvin–Maxwell Model

• Visco-hyperelastic Kelvin–Maxwell model is a constitutive framework that unifies finite-strain hyperelasticity with rate-dependent viscous dissipation under strict thermodynamic principles.
• It employs methodologies such as multiplicative decomposition, additive energy splits, and structure-tensor dynamics to ensure objectivity, positive-definiteness, and dissipation compliance.
• Robust numerical schemes like operator-splitting and matrix-exponential integration enable reliable simulations of polymers, biological tissues, and complex fluids with energy stability and finite propagation speeds.

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 (Shutov, 2017, Liu et al., 2021, Boyaval, 2017, Zhao et al., 20 Aug 2025, Singh et al., 7 Aug 2025, Boyaval et al., 2021).

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$ (Shutov, 2017)), evolution of strain-like internal variables (Zhao et al., 20 Aug 2025, Liu et al., 2021), or structure tensors relaxing towards equilibrium (Boyaval et al., 2021). 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 (Shutov, 2017) or

$$\Psi(C, C^v) = \Psi^\mathrm{eq}(C) + \Psi^\mathrm{neq}(C, C^v)$$

for the series-spring approach (Liu et al., 2021). 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 (Boyaval et al., 2021).

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 (Shutov, 2017, Liu et al., 2021), (ii) equilibrium spring in series or parallel with a generalized Poynting–Thomson (standard-solid) element (Zhao et al., 20 Aug 2025, Singh et al., 7 Aug 2025), or (iii) a generalized internal variable formalism for structure tensors (Boyaval et al., 2021, Boyaval, 2017). 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)$ (Zhao et al., 20 Aug 2025))
• 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$, $\psi(C_e) = c_1(I_1^e - 3) + c_2(I_2^e - 3)$0, $\psi(C_e) = c_1(I_1^e - 3) + c_2(I_2^e - 3)$1, or $\psi(C_e) = c_1(I_1^e - 3) + c_2(I_2^e - 3)$2), typically a first-order flow rule proportional to the driving stress or configurational force: $\psi(C_e) = c_1(I_1^e - 3) + c_2(I_2^e - 3)$3 or, in multiplicative representation,

$\psi(C_e) = c_1(I_1^e - 3) + c_2(I_2^e - 3)$4

where $\psi(C_e) = c_1(I_1^e - 3) + c_2(I_2^e - 3)$5 is a relaxation time and $\psi(C_e) = c_1(I_1^e - 3) + c_2(I_2^e - 3)$6 is the effective viscosity (Shutov, 2017, Liu et al., 2021, Zhao et al., 20 Aug 2025).

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 $\psi(C_e) = c_1(I_1^e - 3) + c_2(I_2^e - 3)$7 computational cost, crucial for large numbers of non-equilibrium processes (Zhao et al., 20 Aug 2025).

3. Governing Equations and Constitutive Updates

The state variables include the deformation tensor ($\psi(C_e) = c_1(I_1^e - 3) + c_2(I_2^e - 3)$8 or $\psi(C_e) = c_1(I_1^e - 3) + c_2(I_2^e - 3)$9), 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 (Boyaval et al., 2021)
• Structure-tensor or conformation-tensor evolution equations with upper-convected derivatives (Boyaval et al., 2021, Boyaval, 2017)
• Constitutive laws for the stresses:

$$\Psi(C, C^v) = \Psi^\mathrm{eq}(C) + \Psi^\mathrm{neq}(C, C^v)$$0

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 (Shutov, 2017)
• Explicit, unconditionally stable updates for coupled internal variables via analytical matrix exponentials (Zhao et al., 20 Aug 2025)
• Heun's predictor–corrector or backward Euler for Lagrangian strain-space models (Singh et al., 7 Aug 2025)
• Operator-splitting with finite-volume (Suliciu-type Riemann) solvers for flow problems (Boyaval, 2017)

These methods preserve frame-invariance, positive-definiteness, symmetry of internal variables, and strict inelastic incompressibility, a critical requirement for large deformations (Shutov, 2017). 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 ($$\Psi(C, C^v) = \Psi^\mathrm{eq}(C) + \Psi^\mathrm{neq}(C, C^v)$$1) limit of the finite-strain Kelvin–Maxwell models reproduces the classical linear viscoelastic solid: $$\Psi(C, C^v) = \Psi^\mathrm{eq}(C) + \Psi^\mathrm{neq}(C, C^v)$$2 where $$\Psi(C, C^v) = \Psi^\mathrm{eq}(C) + \Psi^\mathrm{neq}(C, C^v)$$3 is the stress, $$\Psi(C, C^v) = \Psi^\mathrm{eq}(C) + \Psi^\mathrm{neq}(C, C^v)$$4 and $$\Psi(C, C^v) = \Psi^\mathrm{eq}(C) + \Psi^\mathrm{neq}(C, C^v)$$5 are Mooney–Rivlin moduli, and $$\Psi(C, C^v) = \Psi^\mathrm{eq}(C) + \Psi^\mathrm{neq}(C, C^v)$$6 the viscosity (Shutov, 2017, Liu et al., 2021).

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 (Singh et al., 7 Aug 2025). 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 (Boyaval et al., 2021). 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 (Boyaval, 2017)
• NURBS-based inf-sup stable elements and generalized-alpha schemes for high-fidelity incompressible viscoelasticity with proven energy stability and large-deformation robustness (Liu et al., 2021)
• Matrix-exponential integration harnessing the rank-one perturbation structure for multi-process networks, enabling $$\Psi(C, C^v) = \Psi^\mathrm{eq}(C) + \Psi^\mathrm{neq}(C, C^v)$$7 cost per integration step (Zhao et al., 20 Aug 2025)

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 (Boyaval, 2017, Singh et al., 7 Aug 2025). Design parameters are established experimentally (e.g., for VHB-4910 polymers: $$\Psi(C, C^v) = \Psi^\mathrm{eq}(C) + \Psi^\mathrm{neq}(C, C^v)$$8, $$\Psi(C, C^v) = \Psi^\mathrm{eq}(C) + \Psi^\mathrm{neq}(C, C^v)$$9, $$D = -\frac{\partial \Psi}{\partial C^v} : \frac{1}{2} \dot C^v \ge 0.$$0, $$D = -\frac{\partial \Psi}{\partial C^v} : \frac{1}{2} \dot C^v \ge 0.$$1, $$D = -\frac{\partial \Psi}{\partial C^v} : \frac{1}{2} \dot C^v \ge 0.$$2, $$D = -\frac{\partial \Psi}{\partial C^v} : \frac{1}{2} \dot C^v \ge 0.$$3, $$D = -\frac{\partial \Psi}{\partial C^v} : \frac{1}{2} \dot C^v \ge 0.$$4 (Singh et al., 7 Aug 2025)).

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 (Boyaval et al., 2021, Boyaval, 2017). 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 (Boyaval, 2017)
• High-strain-rate and non-isothermal continuum mechanics, underpinning fully causal viscoelastic flow models with energy transport (Boyaval et al., 2021)

Variants extend to multiple relaxation processes (generalized Maxwell), anisotropy via objective structure-tensor evolution, and direct embedding into thermodynamically stable finite element frameworks (Zhao et al., 20 Aug 2025, Liu et al., 2021).

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 (Shutov, 2017, Singh et al., 7 Aug 2025, Zhao et al., 20 Aug 2025).

---

Key references: (Shutov, 2017, Liu et al., 2021, Boyaval, 2017, Zhao et al., 20 Aug 2025, Singh et al., 7 Aug 2025, Boyaval et al., 2021)