Theory of Evolving Natural Configurations

• Theory of Evolving Natural Configurations is a continuum mechanics framework that treats natural configurations as dynamic, evolving reference states influenced by dissipative processes.
• It employs a multiplicative decomposition of the deformation gradient and integrates Helmholtz free energy with dissipation functions to satisfy the second law of thermodynamics.
• The approach enables unified modeling of viscoelastic, plastic, and rod behaviors using both strain-space and stress-space formulations, validated against experimental data.

The theory of evolving natural configurations (ENC), developed by Eckart and expanded by Rajagopal, is a thermodynamically consistent framework for modeling dissipative processes in solids. Its central tenet is that a material point can possess multiple evolving local natural configurations—intermediate, not generally globally compatible or stress-free—and that the sequence of configurations is described via kinematic decompositions which separate the elastic from the dissipative response. In this setting, the natural configuration is not fixed: in viscoelastic and inelastic bodies, the configuration a body would assume after removing all external forces can itself change over time, and the evolution of the natural configuration is governed by the second law of thermodynamics (Singh et al., 7 Aug 2025).

1. Conceptual basis and terminology

Multiple natural configurations arise due to microstructural changes during inelastic processes such as plasticity, viscoplasticity, and viscoelasticity. In one standard formulation, material states are represented by a sequence of configurations: initial or undeformed, intermediate natural, and current. The intermediate configuration may be natural and stress-free, but it may also be incompatible; in related viscoelastic formulations, an evolving natural configuration is locally achieved by instantaneously elastically unloading the current configuration and does not, in general, yield a stress-free state for viscoelastic materials (Singh et al., 6 Feb 2025).

A central shift introduced by ENC is the replacement of a globally fixed natural configuration by an evolving one. In the rod theory of Rajagopal and Rodriguez, this is described as a generalization of the concept of natural configuration from elasticity, where it is globally fixed, to viscoelastic bodies, where it evolves according to dissipative mechanisms and entropy production (Rajagopal et al., 2022). This makes the theory suitable for constitutive settings in which elasticity, viscosity, yielding, and evolving microstructure must be treated within a single framework.

The cited formulations therefore treat the natural configuration as a dynamical object rather than a passive reference. This suggests that ENC is best understood as a constitutive theory of evolving internal geometry, rather than as a mere reparameterization of classical elasticity.

2. Kinematics, energetics, and thermodynamic structure

The core kinematic device is a multiplicative decomposition of the deformation gradient, written in representative forms as

$$\mathbf{F}=\mathbf{F}^e\mathbf{F}^i \qquad\text{or}\qquad \mathbf{F}=\mathbf{F}^e\mathbf{F}^v.$$

Here, the elastic part maps from an intermediate configuration to the current configuration, while the inelastic or viscous part maps from the reference configuration to the intermediate one. Associated measures include $\mathbf{C}^e$, $\mathbf{C}^i$, $\mathbf{D}^e$, and $\mathbf{D}^i$ (Singh et al., 6 Feb 2025).

The thermodynamic structure is built from a Helmholtz free energy and a dissipation function. In the finite-deformation multiple-natural-configurations framework, constitutive modeling employs a Helmholtz free energy $\psi(\mathbf{E},\mathbf{E}^i)$ for elastic response and a dissipation rate $\xi(\mathbf{E}^i,\mathbf{D}^i)$ for inelastic response, with total stress power

$\dot{W}=J\,\mathbf{T}:\mathbf{D}=\dot{\psi}+\xi.$

In the large-strain viscoelastic formulation, the basic thermodynamic constraint is written as

$$\xi=\mathbf{S}:\dot{\mathbf{E}}-\rho_0\dot{\psi}\ge 0.$$

A constrained maximization principle—described as maximum dissipation or maximum entropy production—then yields the evolution equations for the internal variables (Singh et al., 6 Feb 2025, Singh et al., 7 Aug 2025).

The global form of the second law is also significant. For special Cosserat rods with rate-dependent evolving natural configurations, the constitutive relations are derived by prescribing frame indifferent forms of the Helmholtz energy and the total dissipation rate and requiring that the state variables evolve in a way that maximizes the rate of total entropy production. The framework explicitly distinguishes between the point-wise Clausius-Duhem inequality and the strictly weaker requirement that the rate of total entropy production is non-decreasing (Rajagopal et al., 2023).

3. Constitutive formulations for viscoelastic solids

ENC has been used to revisit nonlinear constitutive models of viscoelastic solids at large strains. Within a Lagrangian framework, a Maxwell model and a Kelvin-Voigt model, together with their associated standard solids—a Zener solid and a Poynting-Thompson solid—have been modeled using the theory. A central conclusion is that a strain-space formulation of the evolving natural configurations is useful in modeling Maxwell-type materials, whereas a stress-space formulation that incorporates a rate of dissipation function in terms of the relevant configurational forces is required for modeling the Kelvin-Voigt type materials (Singh et al., 7 Aug 2025).

In the strain-space formulation, the dissipation is chosen as a quadratic function in the viscous strain rate. In the stress-space formulation, the configurational force is

$$\mathbf{A}^v=-\rho_0\frac{\partial \psi}{\partial \mathbf{E}^v},$$

the dissipation is

$\mathbf{C}^e$0

and a typical evolution equation is

$\mathbf{C}^e$1

The basic Maxwell and Kelvin-Voigt models can be obtained as limiting cases from the derived standard solid models, integration algorithms have been developed, and numerical solutions for a relevant boundary value problem are obtained. The response of the developed models has been compared and benchmarked with experimental data; specifically, the response of the novel Poynting-Thompson model is studied in details and shows a very good match with the existing experimental data obtained from a uniaxial stretching of polymers over a large extent of strain (Singh et al., 7 Aug 2025).

Formulation	Dissipation potential	Primary models
Strain-space	$\mathbf{C}^e$2	Maxwell, Zener
Stress-space	$\mathbf{C}^e$3	Kelvin-Voigt, Poynting-Thompson

This division between strain-space and stress-space is not merely formal. It suggests that the appropriate dissipative variable depends on the rheological architecture being represented.

4. Multiple natural configurations and rheological connections

A major extension of ENC concerns the incorporation of rheological connections at finite deformation. Classical rheological models are mostly developed within a one-dimensional small-strain framework, and one of the key impediments in extending them to three-dimensional finite deformation is determining how series and parallel connections are incorporated into the material models. Within a multiple natural configurations framework, both the series and the parallel connection between two rheological elements can be modeled without changing or introducing new configurations (Singh et al., 6 Feb 2025).

The key distinction is expressed through stress power. For two elements with stress powers $\mathbf{C}^e$4 and $\mathbf{C}^e$5, the finite-deformation postulates are

$\mathbf{C}^e$6

and

$\mathbf{C}^e$7

Accordingly, the difference in a series and a parallel connection is manifested in the ratio of the stress powers expended during the deformations of the associated rheological elements (Singh et al., 6 Feb 2025).

This formulation reorganizes the usual small-strain intuition. In finite deformation, the fundamental distinction is not direct stress or strain superposition, but how the input power is split between branches. Examples developed in this framework include the Standard Linear Solid, an elastic-perfectly plastic material, and elastoplasticity with strain hardening. The account in (Singh et al., 6 Feb 2025) states that the multiple natural configurations framework removes the need for artificial configuration assignments or arbitrary product decompositions outside their natural context and resolves previous mathematical and physical inconsistencies, especially for parallel connections under finite deformation.

5. Rod theories, exact formulations, and long-time behavior

ENC has also been formulated for rod-like bodies. For a planar, inextensible, unshearable, viscoelastic rod, the only nontrivial strain is the curvature

$\mathbf{C}^e$8

and the natural configuration is determined by a natural curvature $\mathbf{C}^e$9, assumed spatially uniform. Rajagopal and Rodriguez formulate the quasistatic equations of motion as an infinite-dimensional dynamical system, prove global well-posedness, show that for every value of the terminal thrust the equations contain a smooth embedded curve of static solutions, characterize the spectrum of the linearized equations about an arbitrary equilibrium point, and—using a convergence result due to Brunovský and Poláčik—prove that every solution converges to an equilibrium point as time goes to infinity (Rajagopal et al., 2022).

The asymptotic analysis is based on a Lyapunov function. In that model,

$\mathbf{C}^i$0

which ensures that all trajectories approach the set of equilibria. The theory therefore supports not only constitutive specification but also rigorous dynamical-systems analysis of existence, boundedness, spectrum, $\mathbf{C}^i$1-limit sets, and convergence (Rajagopal et al., 2022).

A geometrically exact ENC framework has also been developed for special Cosserat rods with rate-dependent evolving natural configurations. The rod configuration consists of a centerline $\mathbf{C}^i$2 and directors $\mathbf{C}^i$3, with natural configuration fields $\mathbf{C}^i$4 and $\mathbf{C}^i$5, and decompositions

$\mathbf{C}^i$6

Two models arise depending on whether the second law is enforced pointwise or only globally. In the strong form,

$\mathbf{C}^i$7

whereas in the weaker global form,

$\mathbf{C}^i$8

The paper shows that, in contrast to other viscoelastic Cosserat rod models introduced in the past, certain quadratic strain energies in this model yield both solid-like stress relaxation and creep (Rajagopal et al., 2023).

6. Scope, adjacent formulations, and terminological breadth

The literature suggests that the phrase “evolving natural configurations” does not denote a single standardized doctrine across all domains. In continuum mechanics, it names a constitutive framework centered on multiple natural configurations, dissipation, and thermodynamic constraints. In other domains, related language denotes different mathematical programs.

One line of work defines a natural process as “an act, by which a system organizes itself with time,” relates organization and action through

$\mathbf{C}^i$9

and proposes that the development of a system towards states of greater organization is cyclic in nature and thus evolution is a cyclic process (Chatterjee, 2011). Reginatto and Hall, by contrast, consider probabilities $\mathbf{D}^e$0 as evolving configurations, endow them with the Fisher-Rao metric, introduce a canonically conjugate field $\mathbf{D}^e$1, extend the geometry to a Kähler structure, and recover wave functions, the quantum free-particle Hamiltonian, and Hilbert space from the geometry of probabilities in motion (Reginatto et al., 2011).

In cultural and cognitive theory, Context-driven Actualization of Potential and the state-context-property formalism describe change of state as context-driven and potentially nondeterministic, with Hilbert-space representations, superposition, and entanglement used to model concept combination and creativity (Gabora, 2013). In biology, natural autoencoding proposes that evolution proceeds by retaining repeating biological interactions while non-repeatable interactions disappear, with species interaction codes defined as compact descriptions of the repeating interactions of species organisms with their external and internal environments (Cohen et al., 2022). Stability-Driven Assembly, in turn, models stochastic assembly with differential persistence and argues that populations are biased toward longer-lived motifs, realizing evolution as a natural, emergent genetic algorithm driven solely by stability (Adler, 22 Jan 2026). A review of living systems further identifies interconnectivity, plasticity, and interdependency as three main features and states that slow evolutionary dynamics is an emergent phenomenon governed by the replicator-mutator equation as the direct consequence of a constrained variational nonequilibrium process (Domenico, 27 Dec 2025).

This suggests a family of related research programs concerned with evolving structure, dissipation, stability, and information, rather than a single universally shared formalism. Within that broader landscape, the ENC of continuum mechanics remains distinctive for its emphasis on multiple natural configurations, constitutive closure by dissipation or entropy-production principles, and explicit finite-deformation modeling of viscoelasticity, plasticity, and rod mechanics.