Hubbry Logo
Alternative stress measuresAlternative stress measuresMain
Open search
Alternative stress measures
Community hub
Alternative stress measures
logo
7 pages, 0 posts
0 subscribers
Be the first to start a discussion here.
Be the first to start a discussion here.
Alternative stress measures
Alternative stress measures
from Wikipedia

In continuum mechanics, the most commonly used measure of stress is the Cauchy stress tensor, often called simply the stress tensor or "true stress". However, several alternative measures of stress can be defined:[1][2][3]

  1. The Kirchhoff stress ().
  2. The nominal stress ().
  3. The Piola–Kirchhoff stress tensors
    1. The first Piola–Kirchhoff stress (). This stress tensor is the transpose of the nominal stress ().
    2. The second Piola–Kirchhoff stress or PK2 stress ().
  4. The Biot stress ()

Definitions

[edit]

Consider the situation shown in the following figure. The following definitions use the notations shown in the figure.

Quantities used in the definition of stress measures

In the reference configuration , the outward normal to a surface element is and the traction acting on that surface (assuming it deforms like a generic vector belonging to the deformation) is leading to a force vector . In the deformed configuration , the surface element changes to with outward normal and traction vector leading to a force . Note that this surface can either be a hypothetical cut inside the body or an actual surface. The quantity is the deformation gradient tensor, is its determinant.

Cauchy stress

[edit]

The Cauchy stress (or true stress) is a measure of the force acting on an element of area in the deformed configuration. This tensor is symmetric and is defined via

or

where is the traction and is the normal to the surface on which the traction acts.

Kirchhoff stress

[edit]

The quantity,

is called the Kirchhoff stress tensor, with the determinant of . It is used widely in numerical algorithms in metal plasticity (where there is no change in volume during plastic deformation). It can be called weighted Cauchy stress tensor as well.

Piola–Kirchhoff stress

[edit]

Nominal stress/First Piola–Kirchhoff stress

[edit]

The nominal stress is the transpose of the first Piola–Kirchhoff stress (PK1 stress, also called engineering stress) and is defined via

or

This stress is unsymmetric and is a two-point tensor like the deformation gradient.
The asymmetry derives from the fact that, as a tensor, it has one index attached to the reference configuration and one to the deformed configuration.[4]

Second Piola–Kirchhoff stress

[edit]

If we pull back to the reference configuration we obtain the traction acting on that surface before the deformation assuming it behaves like a generic vector belonging to the deformation. In particular we have

or,

The PK2 stress () is symmetric and is defined via the relation

Therefore,

Biot stress

[edit]

The Biot stress is useful because it is energy conjugate to the right stretch tensor . The Biot stress is defined as the symmetric part of the tensor where is the rotation tensor obtained from a polar decomposition of the deformation gradient. Therefore, the Biot stress tensor is defined as

The Biot stress is also called the Jaumann stress.

The quantity does not have any physical interpretation. However, the unsymmetrized Biot stress has the interpretation

Relations

[edit]

Relations between Cauchy stress and nominal stress

[edit]

From Nanson's formula relating areas in the reference and deformed configurations:

Now,

Hence,

or,

or,

In index notation,

Therefore,

Note that and are (generally) not symmetric because is (generally) not symmetric.

Relations between nominal stress and second P–K stress

[edit]

Recall that

and

Therefore,

or (using the symmetry of ),

In index notation,

Alternatively, we can write

Relations between Cauchy stress and second P–K stress

[edit]

Recall that

In terms of the 2nd PK stress, we have

Therefore,

In index notation,

Since the Cauchy stress (and hence the Kirchhoff stress) is symmetric, the 2nd PK stress is also symmetric.

Alternatively, we can write

or,

Clearly, from definition of the push-forward and pull-back operations, we have

and

Therefore, is the pull back of by and is the push forward of .

Summary of conversion formula

[edit]

Key:

Conversion formulae
Equation for
(non isotropy)
(non isotropy)
(non isotropy) (non isotropy)

See also

[edit]

References

[edit]
Revisions and contributorsEdit on WikipediaRead on Wikipedia
from Grokipedia
Alternative stress measures in refer to a family of tensorial quantities that quantify the internal forces and moments within a deforming body, particularly under finite strains where the standard may not suffice for mathematical convenience or physical interpretation. These measures differ primarily in the reference configuration they employ—either the current (deformed) state or the initial (undeformed) state—and are formulated to ensure objectivity and compatibility with associated strain measures in the development of constitutive relations. The (often denoted as σ\sigma) serves as the foundational "true stress" measure, defined as per unit area in the current deformed configuration, making it directly and physically intuitive for traction calculations on deformed surfaces. It is symmetric and work-conjugate to the rate of deformation tensor, but its use in large-deformation problems can complicate formulations to the evolving geometry. To address these challenges, alternative measures like the first Piola-Kirchhoff stress (also called nominal or engineering stress, denoted PP) relate forces in the deformed state to areas in the reference configuration, providing a non-symmetric two-point tensor useful for equilibrium equations in the material description. The second Piola-Kirchhoff stress (SS), a symmetric material tensor, pulls back the Cauchy stress to the reference frame via S=JF1σFTS = J F^{-1} \sigma F^{-T} (where JJ is the of the deformation FF), rendering it work-conjugate to the Green-Lagrange strain and ideal for models. Another prominent alternative is the Kirchhoff stress (τ=Jσ\tau = J \sigma), a weighted version of the Cauchy stress that maintains conjugacy with the rate of deformation while simplifying objective rate formulations in spatial descriptions. These measures collectively enable robust numerical implementations in finite element analysis and theoretical advancements in nonlinear elasticity, ensuring conservation laws and stability in simulations of complex deformations.

Background

Role of Stress Measures

In , stress is defined as a measure of the internal forces acting within a deformable body, quantified as the average force per unit area across a surface within the material. This concept captures the distribution of forces that arise due to external loads or internal interactions during deformation, enabling the analysis of how materials respond to mechanical influences such as tension, compression, or shear. The origins of stress measures trace back to the 18th and 19th centuries, with foundational contributions from Leonhard Euler, who in 1776 introduced integral forms of principles and concepts of in fluids, laying early groundwork for understanding force distributions in continua. advanced this significantly between 1823 and 1827 by formulating the general mathematical theory of the stress tensor, initially for small deformations but encompassing finite-strain frameworks through his development of deformation tensors using material coordinates. In the , Clifford Truesdell and others unified and extended these ideas into a comprehensive rational , refining stress measures for finite deformations and broadening their application in and physics. A key distinction exists between true stress and engineering stress. True stress accounts for the instantaneous cross-sectional area in the deformed configuration, providing a physically realistic representation of force transmission in the current state of the material. Engineering stress, conversely, uses the original undeformed cross-sectional area as the reference, simplifying calculations for small strains but diverging from actual conditions in large deformations. The represents the primary true stress measure in this context. To comprehend alternative stress measures, it is essential to appreciate the role of the deformation gradient, a tensor that maps infinitesimal line elements from the reference to the current configuration. Stress measures must transform appropriately under this gradient to maintain consistency and objectivity, ensuring that descriptions of internal forces remain invariant under rigid body motions and accurately reflect material behavior in finite strain scenarios.

Cauchy Stress Tensor

The , denoted as σ\boldsymbol{\sigma}, is a symmetric second-order tensor that represents the state of stress in a continuous body, defined as the force per unit area acting on a surface in the current (deformed) configuration. It arises from Cauchy's fundamental theorem, which posits that the traction vector t\mathbf{t} on any surface with unit normal n\mathbf{n} is given by t=σn\mathbf{t} = \boldsymbol{\sigma} \mathbf{n}, allowing the tensor to fully characterize internal forces across arbitrary orientations. Physically, the components σij\sigma_{ij} of the describe the stress on a surface element in the deformed body: the diagonal elements σii\sigma_{ii} (no sum) represent normal stresses, while the off-diagonal elements σij\sigma_{ij} (i ≠ j) denote shear stresses, quantifying the force in the ii-direction per unit area perpendicular to the jj-direction. This interpretation aligns with the tensor's role in capturing the directional nature of forces, such as and , in the instantaneous state. Key properties of the include its symmetry, σ=σT\boldsymbol{\sigma} = \boldsymbol{\sigma}^T or σij=σji\sigma_{ij} = \sigma_{ji}, which follows from the balance of and ensures no net on material elements in the absence of body couples. It is also objective, transforming under rotations R\mathbf{R} as σ=RσRT\boldsymbol{\sigma}^* = \mathbf{R} \boldsymbol{\sigma} \mathbf{R}^T, preserving its frame-indifference in describing physical stress states. Additionally, the trace of the tensor, tr(σ)=σkk\operatorname{tr}(\boldsymbol{\sigma}) = \sigma_{kk}, relates to the hydrostatic p=13tr(σ)p = -\frac{1}{3} \operatorname{tr}(\boldsymbol{\sigma}), representing the isotropic part of the stress that governs volumetric changes. The Cauchy stress tensor satisfies the local balance of linear momentum, expressed as the Cauchy equation of motion: σ+ρb=ρa,\nabla \cdot \boldsymbol{\sigma} + \rho \mathbf{b} = \rho \mathbf{a}, where σ\nabla \cdot \boldsymbol{\sigma} is the divergence of the tensor, ρ\rho is the current mass density, b\mathbf{b} is the body force per unit mass, and a\mathbf{a} is the material acceleration. This equation derives from applying Newton's second law to an infinitesimal volume element in the Eulerian frame. The advantages of the lie in its intuitive representation of forces relative to the deformed geometry, making it particularly suitable for analyses in , where it directly incorporates pressure and viscous effects, and in small-strain , where deformations do not significantly alter areas or volumes. However, in large-deformation scenarios, its dependence on the current configuration can complicate constitutive relations with reference-based strain measures.

Motivations for Alternatives

Objectivity Requirements

In , the principle of objectivity, also known as material frame indifference, stipulates that the mechanical behavior of a must be independent of the observer's choice of reference frame, particularly under superposed motions such as translations and rotations. This ensures that stress measures reflect intrinsic properties rather than artifacts of the . Translations do not affect tensor components, but rotations require specific transformation rules to maintain physical consistency across observers. Mathematically, a second-order tensor T\mathbf{T} defined in the spatial (current) configuration is deemed objective if, upon a superposed rigid rotation represented by an orthogonal tensor Q\mathbf{Q} (with detQ=1\det \mathbf{Q} = 1), it transforms as T=QTQT.\mathbf{T}^* = \mathbf{Q} \mathbf{T} \mathbf{Q}^T. This rule preserves the tensor's representation of physical quantities like stress under changes in observer orientation. The principle extends to constitutive relations, requiring that response functions yield the same physical outcomes regardless of rigid superpositions. Formulated by Walter Noll in , this axiom forms a of modern , ensuring the invariance of balance laws and constitutive equations under Euclidean transformations. Noll's framework emphasized that objectivity arises from the independence of mechanical descriptions from arbitrary rigid motions, enabling a unified for deformable media. The σ\boldsymbol{\sigma} inherently satisfies objectivity, transforming as σ=QσQT\boldsymbol{\sigma}^* = \mathbf{Q} \boldsymbol{\sigma} \mathbf{Q}^T, since it is defined with respect to force per unit area in the current (deformed) configuration. However, in finite deformation analyses, this spatial reference poses significant challenges: the current configuration evolves with loading, making it difficult to directly link stress to fixed coordinates or initial where constitutive are typically defined and experimentally measured. This mismatch complicates the of models, integration in numerical schemes like finite elements, and interpretation of results relative to undeformed states. Alternative stress measures address these issues by referring forces to the reference (undeformed) configuration while adhering to objectivity through appropriate pull-back operations, allowing constitutive modeling in a fixed frame without loss of frame indifference. Without such measures, analyses in finite strains risk inconsistencies when deformations include large rotations. Non-objective stress measures, such as the (nominal) stress—which approximates force per undeformed area without accounting for deformation-induced changes—fail the transformation rule, typically transforming as T=QT\mathbf{T}^* = \mathbf{Q} \mathbf{T} rather than QTQT\mathbf{Q} \mathbf{T} \mathbf{Q}^T. This leads to observer-dependent predictions, where the computed stress varies with the orientation of the simulation frame even for identical physical loading. In finite element simulations of structures undergoing rotations, such measures introduce spurious coupling between rotation and stress, resulting in artificial dissipation, inaccurate equilibrium paths, or divergent solutions under incremental loading.

Work Conjugacy with Strain Measures

In , particularly for finite deformations, a stress measure S\mathbf{S} and a strain measure E\mathbf{E} are defined as work-conjugate if the internal virtual work per unit reference volume is expressed as the double contraction S:δE\mathbf{S} : \delta \mathbf{E}, ensuring that the integral over the domain accurately represents the total without additional transformation factors. This conjugacy condition is fundamental for deriving consistent constitutive relations, as it guarantees that the stress power equals the rate of change of the density, W˙=S:E˙\dot{W} = \mathbf{S} : \dot{\mathbf{E}}. The importance of work conjugacy is pronounced in hyperelasticity, where materials store energy path-independently through a potential W(E)W(\mathbf{E}), with the stress obtained as S=WE\mathbf{S} = \frac{\partial W}{\partial \mathbf{E}} only for conjugate pairs; non-conjugate choices lead to inconsistencies in energy conservation. In finite element methods for large-deformation problems, conjugate pairs preserve variational consistency in incremental formulations, promoting by avoiding spurious energy modes and ensuring reliable convergence in nonlinear solvers. For instance, violations of conjugacy in analyses can introduce large errors in critical load predictions, as observed in simulations of compressed elastic solids. Common work-conjugate pairs include the σ\boldsymbol{\sigma} with the spatial rate of deformation d\mathbf{d} (symmetric part of the velocity gradient), suitable for hypoelastic models; the second Piola-Kirchhoff stress S\mathbf{S} with the Green-Lagrange strain E\mathbf{E}, ideal for hyperelastic formulations in the reference configuration; and the first Piola-Kirchhoff stress P\mathbf{P} with variations of the deformation gradient δF\delta \mathbf{F}, often used in total Lagrangian schemes. These pairs facilitate objective and efficient implementation in . An illustrative example arises in modeling hyperelastic materials like rubber under cyclic loading, where non-conjugate pairs—such as those in certain hypoelastic approximations—introduce artificial and , resulting in unphysical residual stresses and energy loss upon unloading, contrary to the reversible nature of hyperelasticity. This underscores the need for conjugacy to maintain physical fidelity in large-strain simulations.

Key Alternative Measures

Kirchhoff Stress

The Kirchhoff stress tensor, denoted as τ\boldsymbol{\tau}, is defined as τ=Jσ\boldsymbol{\tau} = J \boldsymbol{\sigma}, where σ\boldsymbol{\sigma} is the Cauchy stress tensor and J=detFJ = \det \mathbf{F} represents the Jacobian determinant of the deformation gradient F\mathbf{F}, quantifying the local volume ratio between the current and reference configurations. This measure is a spatial tensor, symmetric like the Cauchy stress, and operates in the current (deformed) configuration. As a scalar multiple of the , the Kirchhoff stress shares the same principal directions but is scaled by the factor JJ, which accounts for volumetric changes during deformation. It is an objective stress measure, invariant under superposed rotations, preserving the transformation properties of the Cauchy stress while incorporating the volume scaling. Physically, the Kirchhoff stress adjusts the true stress for both area and volume changes in finite strain analyses, providing a weighted representation of force per unit deformed area that bridges spatial and material descriptions without direct physical interpretation beyond this adjustment. One key advantage of the Kirchhoff stress lies in updated Lagrangian formulations, where it simplifies the derivation of equilibrium equations by absorbing the Jacobian factor, facilitating clearer expressions for stress power as τ:d\boldsymbol{\tau} : \mathbf{d} (with d\mathbf{d} the rate of deformation tensor) and easing the handling of boundary conditions. It is particularly prevalent in metal plasticity models, where plastic incompressibility implies J1J \approx 1, making τσ\boldsymbol{\tau} \approx \boldsymbol{\sigma} and streamlining numerical implementations for isochoric flow in crystal plasticity simulations.

First Piola-Kirchhoff Stress

The first Piola-Kirchhoff stress tensor, also known as the nominal stress tensor or stress tensor, is a measure that relates forces in the current (deformed) configuration to areas in the reference (undeformed) configuration. It is defined mathematically as P=JσFT\mathbf{P} = J \boldsymbol{\sigma} \mathbf{F}^{-T}, where J=detFJ = \det \mathbf{F} is the of the deformation F\mathbf{F}, and σ\boldsymbol{\sigma} is the . This tensor serves as a pull-back of the Cauchy stress to the reference configuration, facilitating analysis in Lagrangian descriptions of finite deformations. Unlike the symmetric , the first Piola-Kirchhoff stress tensor P\mathbf{P} is generally not symmetric, even when σ\boldsymbol{\sigma} is symmetric, due to the involvement of the inverse of F\mathbf{F}. However, it satisfies the relation PFT=FPT\mathbf{P} \mathbf{F}^T = \mathbf{F} \mathbf{P}^T, reflecting the of σ\boldsymbol{\sigma}. Regarding frame indifference, P\mathbf{P} transforms under a superposed rigid Q\mathbf{Q} as P=QP\mathbf{P}^* = \mathbf{Q} \mathbf{P}, embodying objectivity in a mixed sense as a two-point tensor connecting reference and current frames. Physically, P\mathbf{P} represents the force per unit area in the undeformed configuration, with the force vector expressed in the deformed configuration's basis. In component form, using Cartesian bases ei\mathbf{e}_i for the current configuration and EJ\mathbf{E}_J for the reference, the components are PiJP_{iJ}, where PiJP_{iJ} denotes the i-component of force per unit reference area normal to the J-direction. This interpretation arises from Nanson's formula relating area vectors across configurations, ensuring P\mathbf{P} maps reference normals to current traction vectors appropriately. In applications, particularly finite element analysis of large deformations, P\mathbf{P} is valuable for imposing boundary conditions, as prescribed loads naturally align with the initial in the reference configuration. It enables straightforward computation of equilibrium equations in the Lagrangian framework, such as 0P+b0=ρ0x¨\nabla_0 \cdot \mathbf{P} + \mathbf{b}_0 = \rho_0 \ddot{\mathbf{x}}, where 0\nabla_0 is the reference , b0\mathbf{b}_0 body forces, and ρ0\rho_0 reference density. This utility extends to hybrid stress formulations and variational principles in nonlinear simulations.

Second Piola-Kirchhoff Stress

The second Piola-Kirchhoff stress tensor, denoted as S\mathbf{S}, is defined as the pull-back of the first Piola-Kirchhoff stress tensor P\mathbf{P} by the inverse deformation gradient F1\mathbf{F}^{-1}, yielding S=F1P\mathbf{S} = \mathbf{F}^{-1} \mathbf{P}. Equivalently, it relates to the σ\boldsymbol{\sigma} through S=JF1σFT\mathbf{S} = J \mathbf{F}^{-1} \boldsymbol{\sigma} \mathbf{F}^{-T}, where J=detFJ = \det \mathbf{F} is the . This tensor was originally introduced by Gabrio Piola in 1833 and later formalized by in 1852, providing a foundational measure for stress in finite deformation theory. A key property of S\mathbf{S} is its , S=ST\mathbf{S} = \mathbf{S}^T, which follows directly from the symmetry of the under typical assumptions in . Additionally, S\mathbf{S} is objective in the Lagrangian sense, remaining unchanged under superposed rigid-body rotations on the current configuration, as its components are defined with respect to the fixed reference configuration. Physically, S\mathbf{S} represents a "corrected" engineering stress, interpreting the force per unit undeformed area but resolved in the undeformed material directions, thus accounting for the geometric changes in orientation during deformation. In hyperelastic constitutive models, S\mathbf{S} plays a central role, with the Cauchy stress obtained via the push-forward relation σ=1JFSFT\boldsymbol{\sigma} = \frac{1}{J} \mathbf{F} \mathbf{S} \mathbf{F}^T. This simplifies the expression of laws, as S\mathbf{S} is often derived directly from the strain energy density function as a function of the Green-Lagrange strain tensor. Its advantages include ensuring in total Lagrangian finite element , where equilibrium is enforced in the reference configuration, thereby enhancing computational stability for large-deformation problems without frequent updates.

Biot Stress

The Biot stress tensor, often denoted as T\mathbf{T}, is also expressed as the symmetric part of the first Piola-Kirchhoff stress tensor rotated by the deformation, T=12(RP+PR)\mathbf{T} = \frac{1}{2} (\mathbf{R}^\top \mathbf{P} + \mathbf{P}^\top \mathbf{R}), where R\mathbf{R} is the rotation tensor from the polar decomposition of F\mathbf{F}. The Biot stress tensor is generally unsymmetric unless the material is isotropic, in which case it shares symmetry with the underlying Cauchy stress. Physically, it represents the internal force per unit reference area acting in the direction of the current deformed line elements, providing a direct interpretation of stress along extended material fibers without resolving rigid rotations. This fiber-aligned perspective facilitates analysis of deformation work as T:U˙\mathbf{T} : \dot{\mathbf{U}}, where U\mathbf{U} is the right stretch tensor, establishing work conjugacy with the Biot strain measure. Due to its basis in deformed fiber directions, the Biot stress finds niche applications in modeling anisotropic materials where directional reinforcement is prominent. In crystal plasticity, it aids in capturing slip along lattice directions under large strains, as seen in finite-strain formulations for polycrystals. For fiber-reinforced composites, it supports homogenization of effective properties under finite deformations, enabling prediction of transverse loading responses. In , particularly for soft tissues like membranes or cardiovascular structures, it models nonlinear elasticity with fiber orientations, improving fits to experimental data on residual stresses and . These uses stem from seminal work by Biot on incremental deformations, though its complexity limits broader adoption compared to other measures.

Interconversion Relations

Cauchy to First Piola-Kirchhoff

The transformation from the σ\boldsymbol{\sigma} to the first Piola-Kirchhoff stress tensor P\mathbf{P} arises in the context of finite deformation continuum mechanics, where forces and areas are related between the current (deformed) and (undeformed) configurations. This relation ensures that the stress measure accounts for the mapping of surface elements under deformation, preserving the balance of forces in equilibrium equations. The derivation begins with Cauchy's fundamental theorem, which states that the traction t\mathbf{t} on a surface with normal n\mathbf{n} in the current configuration is t=σn\mathbf{t} = \boldsymbol{\sigma} \mathbf{n}, where the magnitude relates to the force per unit deformed area dada. To refer this to the reference configuration, Nanson's formula provides the transformation for area vectors: da=JFTdA\mathbf{da} = J \mathbf{F}^{-T} \mathbf{dA}, where J=detFJ = \det \mathbf{F} is the Jacobian determinant of the deformation gradient F\mathbf{F}, and dA\mathbf{dA} is the reference area vector. Equating the forces, the traction in the reference configuration becomes T=PN\mathbf{T} = \mathbf{P} \mathbf{N}, with P\mathbf{P} satisfying PN=σ(JFTN)\mathbf{P} \mathbf{N} = \boldsymbol{\sigma} (J \mathbf{F}^{-T} \mathbf{N}), leading to the explicit relation P=JσFT\mathbf{P} = J \boldsymbol{\sigma} \mathbf{F}^{-T}. This push-forward/pull-back operation maintains objectivity under rigid body motions. Geometrically, the term JFTJ \mathbf{F}^{-T} represents the cofactor matrix of F\mathbf{F}, which scales areas by the deformation while accounting for orientation changes; it transforms reference areas into deformed ones, ensuring that P\mathbf{P} computes per unit reference area, in contrast to σ\boldsymbol{\sigma}'s use of deformed areas. This is particularly useful for integrating equilibrium over fixed reference domains in numerical simulations. In the example of uniaxial tension along the xx-direction, assume a bar with uniform deformation where F=diag(λ,1/λ,1/λ)\mathbf{F} = \operatorname{diag}(\lambda, 1/\sqrt{\lambda}, 1/\sqrt{\lambda})
Add your contribution
Related Hubs
User Avatar
No comments yet.