1 Fundamental concepts
Nonlinear elasticity describes the deformation of solid bodies when stress and strain are not linked by a simple proportional law. In this setting, the response may depend on the magnitude of deformation, the path taken to reach it, or the geometry of the body after deformation. The subject extends classical elasticity to situations involving large stretches, strong rotations, and pronounced changes in shape.
A central idea is that a solid can be compared in an undeformed reference state and in a deformed current state. The relationship between these states is described by a deformation map, from which strain, stress, and other quantities are derived. Because the deformation may be large, many formulas used in linear theory must be replaced by finite-deformation measures.
1.1 Deformation and strain
Deformation is the change in position of material points in a body. In nonlinear elasticity, strain is not treated as a simple small quantity; instead, it is measured relative to a reference configuration using tensors that account for finite changes in length and angle. These measures distinguish stretching from rigid motion, allowing rotation to be separated from true deformation.
1.2 Stress measures
Stress in a deformed body can be expressed in more than one form, depending on whether it is measured in the current configuration or referred back to the reference state. Different stress measures are useful for different balances and constitutive laws. In nonlinear elasticity, careful choice of stress measure is essential because finite deformation changes the geometry over which forces act.
1.2.1 Cauchy stress
The Cauchy stress is the true stress acting on surfaces in the current, deformed configuration. It describes force per unit area in the body as it exists at that moment. This measure is often used in physical interpretation because it directly corresponds to the actual internal traction on an oriented surface.
1.2.2 First and second Piola–Kirchhoff stress
The first Piola–Kirchhoff stress relates forces in the current configuration to areas in the reference configuration, while the second Piola–Kirchhoff stress is fully referred to the undeformed state. These measures are convenient in theoretical formulations because they pair naturally with deformation gradients and strain tensors defined in the reference frame.
1.3 Reference and current configurations
The reference configuration is the original state of the body, usually assumed stress-free or nearly so, while the current configuration is the deformed state after loading. Comparing these two configurations allows one to define the deformation gradient and derive kinematic quantities. This distinction is fundamental in finite elasticity because geometry itself changes under load.
1.4 Finite versus infinitesimal deformation
Infinitesimal deformation theory assumes that displacements and rotations are small enough to linearize the governing equations. Nonlinear elasticity removes that restriction, permitting large strains and rotations. The finite-deformation setting is necessary when linear approximations no longer provide accurate predictions, such as in rubberlike materials or heavily loaded structural components.
2 Constitutive theory
Constitutive theory describes how a material responds to deformation through a stress-strain relation. In nonlinear elasticity, these relations are often nonlinear functions of strain measures and may depend on material structure, compressibility, and symmetry. The theory aims to capture both mechanical response and material-specific features in a mathematically consistent way.
2.1 Hyperelastic materials
Hyperelastic materials are elastic solids whose stress response can be derived from a stored energy function. Their behavior is path-independent in the sense that the stress at a given deformation depends only on the current strain state, not on how that state was reached. This framework is especially important for rubbers, gels, and soft tissues.
2.1.1 Strain energy density
The strain energy density is the elastic energy stored per unit reference volume as the body deforms. Once this energy function is specified, stresses follow by differentiation with respect to appropriate strain measures. The form of the energy function determines how the material stiffens, softens, or resists volumetric change.
2.1.2 Material isotropy and anisotropy
An isotropic material has the same response in all directions, whereas an anisotropic material responds differently depending on orientation. Many synthetic elastomers are approximately isotropic, while fiber-based or biological materials often show directional dependence. Anisotropy introduces additional structure into the constitutive law through preferred directions or embedded families of fibers.
2.2 Incompressibility and compressibility
Some elastic materials, especially rubberlike solids and soft tissues, are modeled as nearly incompressible, meaning their volume changes very little under load. Others are compressible and can undergo significant volumetric expansion or contraction. The choice affects the pressure-like terms in the governing equations and has a major influence on the predicted deformation.
2.3 Objective stress-strain relations
An objective constitutive relation is one that does not change under superposed rigid-body motions. This requirement ensures that the material law depends only on physically meaningful deformation and not on the observer’s frame of reference. Objectivity is a key consistency condition in continuum mechanics.
2.4 Material symmetry
Material symmetry describes invariance of the constitutive response under specific transformations, such as rotations that preserve the internal structure of the solid. Symmetry reduces the number of independent material parameters and helps classify materials into isotropic, transversely isotropic, orthotropic, and other categories. It is particularly useful in formulating models for composites and biological tissues.
3 Governing equations
The governing equations of nonlinear elasticity express conservation laws together with constitutive behavior. They determine how a body moves, deforms, and maintains equilibrium under applied forces and constraints. Because deformation can be large, these equations are usually written in tensor form and evaluated in either the reference or current configuration.
3.1 Balance of mass
The balance of mass states that material does not appear or disappear during deformation. In continuum form, it links the density in the reference configuration to the density in the current configuration through the deformation. For incompressible materials, this balance imposes an additional constraint on the allowable deformation.
3.2 Balance of momentum
The balance of momentum expresses Newton’s second law for a deformable body. It relates internal stresses to body forces and inertia. In nonlinear elasticity, this balance can be written for static problems where acceleration vanishes or for dynamic problems where motion evolves over time.
3.2.1 Static equilibrium
In static equilibrium, the body is at rest and the sum of internal and external forces is zero. The resulting equations describe the deformation under sustained loading, such as stretching, compression, or indentation. Many practical elasticity problems are formulated in this regime.
3.2.2 Dynamic motion
Dynamic motion includes time-dependent deformation and inertia. The equations then describe how waves, vibrations, and transient responses propagate through the solid. Large-deformation dynamics are relevant in impact, oscillation, and fast mechanical loading.
3.3 Balance of angular momentum
The balance of angular momentum requires that internal moments be consistent with the absence of spontaneous rotational forces. For ordinary elastic materials without couple stresses, this leads to symmetry conditions on the stress tensor in the appropriate configuration. This balance is essential for the internal consistency of the theory.
3.4 Boundary and initial conditions
Boundary conditions prescribe displacements, tractions, or mixed constraints on the body’s surface. Initial conditions specify the starting deformation and velocity in dynamic problems. Together, they complete the mathematical problem and determine whether a unique solution exists.
4 Strain measures and invariants
Different strain measures are used in nonlinear elasticity to describe finite deformation in mathematically convenient ways. Many constitutive laws are written in terms of tensor invariants, which remain unchanged under rigid rotations. These quantities provide a compact and physically meaningful description of deformation.
4.1 Green–Lagrange strain
The Green–Lagrange strain is a finite strain measure defined relative to the reference configuration. It captures both stretching and quadratic contributions from deformation, making it suitable for large-strain analysis. In the small-deformation limit, it reduces to the familiar linear strain tensor.
4.2 Right and left Cauchy–Green tensors
The right Cauchy–Green tensor measures deformation with respect to the reference configuration, while the left Cauchy–Green tensor describes the deformation in the current configuration. Both are central to finite elasticity because they separate pure stretch from rotation. They are often used to construct strain energy functions.
4.3 Principal stretches
Principal stretches are the maximum and minimum stretch ratios experienced by material line elements along orthogonal directions. They provide an intuitive description of how a body elongates or contracts. Many material models can be expressed directly in terms of these stretch values.
4.4 Invariant-based formulations
Invariant-based formulations express constitutive laws using quantities unchanged by rotation, such as traces or determinants of strain tensors. This approach simplifies the description of isotropic materials and helps ensure objectivity. It also provides a flexible basis for anisotropic models when additional structural directions are included.
5 Material models
Material models specify the relationship between stress and strain for particular classes of solids. In nonlinear elasticity, these models are often chosen to fit experimental data while remaining mathematically tractable. Some are designed for isotropic rubberlike behavior, while others capture fiber reinforcement or tissue anisotropy.
5.1 Neo-Hookean models
Neo-Hookean models are among the simplest hyperelastic laws. They describe rubberlike materials with a relatively simple energy function and are often used for moderate large deformations. Despite their simplicity, they can capture basic nonlinear stiffening.
5.2 Mooney–Rivlin models
Mooney–Rivlin models extend the neo-Hookean form by including additional invariant terms. This gives greater flexibility in fitting experimental data for elastomers. They are widely used when a simple model is needed but more accuracy is desired than the neo-Hookean law provides.
5.3 Ogden models
Ogden models represent the strain energy in terms of powers of principal stretches. This makes them highly adaptable for fitting a wide range of elastic responses, especially for rubber and soft polymers. Their parameterization can capture strong nonlinearity over large deformations.
5.4 Saint Venant–Kirchhoff model
The Saint Venant–Kirchhoff model is a classical nonlinear extension of linear elasticity. It uses the Green–Lagrange strain in a quadratic energy function. Although mathematically convenient, it is less reliable for very large deformations than more modern hyperelastic laws.
5.5 Anisotropic models
Anisotropic models account for direction-dependent response caused by internal structure. They are important for materials with aligned fibers, layered organization, or preferred orientations. Such models often combine isotropic matrix behavior with directional reinforcement.
5.5.1 Fiber-reinforced solids
Fiber-reinforced solids contain stiff embedded fibers that carry load preferentially along certain directions. The fibers can greatly increase strength and stiffness while producing complex nonlinear behavior. These materials are common in composites and many biological systems.
5.5.2 Biological tissue models
Biological tissue models often include anisotropy, near incompressibility, and nonlinear stiffening at large stretches. They are used to describe organs, membranes, tendons, and connective tissues. Because tissues are structurally heterogeneous, model selection typically balances realism and simplicity.
6 Analytical methods
Analytical methods seek exact or approximate solutions to nonlinear elasticity problems. Since fully nonlinear equations are often difficult to solve, researchers use special techniques that exploit symmetry, small parameters, or variational structure. These methods provide insight into deformation patterns and limiting behavior.
6.1 Exact solutions
Exact solutions are closed-form solutions that satisfy the full nonlinear equations without approximation. They are usually available only for special geometries, loading conditions, or material laws. Even when limited in scope, exact solutions are valuable for benchmarking theory and numerical methods.
6.2 Perturbation methods
Perturbation methods approximate solutions by expanding them around a known state using a small parameter. They are useful for weak nonlinearity, small departures from equilibrium, or nearly symmetric problems. Such methods can reveal how nonlinear effects alter stresses, stability, and wave speeds.
6.3 Variational principles
Variational principles formulate elasticity as an extremum problem for an energy functional. This approach is especially natural in hyperelasticity, where equilibrium corresponds to stationary potential energy. Variational methods also provide a foundation for numerical discretization and theoretical existence results.
6.3.1 Energy minimization
Energy minimization states that stable equilibrium configurations tend to minimize total potential energy subject to constraints. This principle helps identify admissible deformations and compare competing states. It is widely used in deriving equilibrium conditions and proving stability.
6.3.2 Principle of virtual work
The principle of virtual work requires that the work of external forces equals the work of internal stresses for all admissible virtual displacements. It provides an equivalent weak form of the equilibrium equations. This formulation is especially useful in finite element analysis.
6.4 Linearization about a deformed state
Linearization about a deformed state studies small disturbances around a finite pre-stress or pre-deformation. It is used to analyze incremental behavior, stability, and wave propagation in loaded bodies. This technique bridges fully nonlinear theory and local response near an equilibrium state.
7 Stability and bifurcation
Stability analysis examines whether an equilibrium deformation persists under small disturbances. Bifurcation occurs when multiple deformation paths become possible from a single loading state. In nonlinear elasticity, these phenomena are common because geometry and material response can interact strongly.
7.1 Material and geometric instability
Material instability arises from the constitutive response of the solid, such as softening or nonconvex energy functions. Geometric instability results from the shape of the body and the way loads are applied, even if the material itself is stable. Both types can act together and produce sudden changes in deformation.
7.2 Buckling and wrinkling
Buckling is a sudden change in shape under compression or other destabilizing loads, while wrinkling is a localized pattern of surface undulations. These effects are often seen in slender structures, thin sheets, and soft layers. Nonlinear elasticity provides the framework for predicting their onset and evolution.
7.3 Loss of ellipticity
Loss of ellipticity refers to a mathematical change in the governing equations that signals the possible formation of localized deformation patterns. It is often associated with material softening, shear banding, or the breakdown of smooth solutions. In practice, it can mark the boundary between stable and unstable mechanical behavior.
7.4 Post-buckling behavior
Post-buckling behavior describes the deformation path after instability has occurred. The structure may settle into a new equilibrium shape, continue deforming through secondary bifurcations, or undergo complex localized changes. Understanding this regime is important for safe design and for interpreting experiments.
8 Wave propagation
Wave propagation in nonlinear elasticity studies how disturbances move through deformable solids when amplitudes are not small. Nonlinearity can change wave speed, shape, and interaction behavior. This subject is relevant to impacts, vibrations, acoustics, and the characterization of materials.
8.1 Finite-amplitude waves
Finite-amplitude waves have displacements large enough that linear wave theory is inadequate. Their speed may depend on amplitude, pre-stress, or deformation state. Such waves can steepen, spread, or interact in ways not seen in infinitesimal theory.
8.2 Surface and bulk waves
Surface waves travel near boundaries, while bulk waves propagate through the interior of the solid. In nonlinear media, both types may be influenced by finite strain and preloading. Their analysis is important in nondestructive testing and material characterization.
8.3 Nonlinear acoustics in solids
Nonlinear acoustics in solids deals with sound propagation when the elastic response is not strictly linear. Harmonic generation, waveform distortion, and amplitude-dependent speed are common features. This field connects elasticity with acoustic measurement techniques.
9 Numerical methods
Numerical methods are essential for solving realistic nonlinear elasticity problems, which often lack closed-form solutions. They approximate the governing equations on a discretized domain and update the solution iteratively. These tools are widely used in engineering simulation and scientific computing.
9.1 Finite element formulation
The finite element method divides the body into small elements and approximates displacement fields with basis functions. It is well suited to complex geometry, heterogeneous materials, and nonlinear constitutive laws. In nonlinear elasticity, the method is usually expressed in a weak form and solved incrementally.
9.2 Incremental-iterative schemes
Incremental-iterative schemes solve nonlinear problems by loading the system in small steps and correcting the solution at each step. Newton-type methods are common because they converge efficiently near a solution. Such schemes are needed when stiffness, geometry, or constraints evolve during deformation.
9.3 Contact and large-deformation analysis
Contact problems involve interaction between deforming bodies or self-contact within a single body. Large-deformation analysis must also track changing boundaries and possible surface sliding or separation. These problems are computationally demanding because the contact conditions and geometry both vary with time.
9.4 Computational challenges
Computational challenges include convergence difficulties, mesh distortion, material incompressibility, and sensitivity to initial guesses. Accurate simulation may require adaptive meshing, stabilization methods, or specialized solvers. These issues become more pronounced when deformations are extreme or when material response is highly nonlinear.
10 Applications
Nonlinear elasticity is applied wherever large deformation or strong material nonlinearity affects mechanical response. It is especially important in soft solids, engineered elastomers, and structures that operate beyond the small-strain range. The theory supports both design and interpretation of experiments.
10.1 Rubber and elastomers
Rubber and elastomers are classic examples of materials requiring nonlinear elastic models. They can sustain large strains and exhibit pronounced stiffening at high extension. Their behavior is central to products such as tires, seals, hoses, and flexible mounts.
10.2 Soft biological tissues
Soft biological tissues often show anisotropy, viscoelastic effects, and near incompressibility, making nonlinear elasticity an important approximation. It is used to model deformation in organs, vessels, ligaments, and skin. The framework helps relate internal structure to mechanical function.
10.3 Structural components
Structural components may experience significant deformation under service loads, during forming, or in extreme events. Nonlinear elasticity helps predict response in arches, shells, beams, and other load-bearing elements when linear assumptions fail. It is also useful in assessing stability and residual deformation.
10.4 Sealing, damping, and vibration isolation
Seals, dampers, and isolators often rely on rubberlike materials whose performance depends on nonlinear deformation. Accurate modeling improves predictions of compression, contact pressure, energy absorption, and vibration transmission. These applications depend on both material law and component geometry.
11 History and development
The development of nonlinear elasticity reflects the broader evolution of continuum mechanics from classical small-strain theory to modern finite-deformation modeling. Early work focused on idealized solids and linear approximations, while later research addressed large deformation, stability, and numerical simulation. The field now integrates analysis, computation, and material modeling.
11.1 Early classical elasticity
Early classical elasticity was dominated by linear theories that described small displacements and small strains. These ideas were successful for many engineering problems and established the basic concepts of stress, strain, and equilibrium. However, they could not fully explain large deformations or strongly nonlinear material responses.
11.2 Emergence of finite elasticity
Finite elasticity emerged to address problems where geometry changes significantly and the linear approximation breaks down. Researchers developed kinematic tools, objective stress measures, and energy-based constitutive laws to handle finite strain. This work laid the foundation for modern nonlinear elasticity.
11.3 Modern computational nonlinear mechanics
Modern computational nonlinear mechanics combines theory with numerical algorithms capable of handling realistic geometries, complex materials, and contact interactions. Advances in computing made it possible to solve previously intractable problems and to compare models with experiments in detail. The field continues to expand through applications in biomechanics, soft matter, and advanced structural analysis.