Fact-checked by Grok 2 weeks ago

Yield surface

In continuum mechanics, the yield surface is a mathematical boundary in multi-dimensional stress or strain space that delineates the onset of plastic deformation in a material, separating the region of purely elastic response from that of inelastic behavior under applied loads. It is typically formulated as a function f(\sigma, \kappa) = 0 in stress space, where \sigma represents the stress tensor and \kappa denotes internal state variables accounting for the material's deformation history, ensuring that stresses inside the surface produce only elastic strains while those on or outside trigger yielding. This surface is inherently convex for metals exhibiting linear elasticity unaffected by prior plastic flow, a property rooted in stability principles and Drucker’s postulate, which guarantees unique and physically consistent plastic strain increments. Common yield criteria define the shape of the yield surface to approximate experimental observations under multiaxial loading. The , based on the distortional , posits yielding when the equivalent von Mises \sigma_{VM} = \sqrt{\frac{1}{2}[(\sigma_1 - \sigma_2)^2 + (\sigma_1 - \sigma_3)^2 + (\sigma_2 - \sigma_3)^2]} equals the uniaxial yield \sigma_Y, resulting in an elliptical surface in principal space that effectively captures ductile metal behavior. In contrast, the Tresca criterion relies on the maximum , with yielding occurring when \tau_{\max} = \sigma_Y / 2, yielding a in space suitable for conservative design in brittle or conservative failure predictions. These criteria enable prediction of plastic flow directions via the normality rule, where the plastic strain increment is proportional to the outward normal of the surface at the current point, as per the associated flow rule in classical plasticity theory. The yield surface evolves during plastic deformation through hardening mechanisms, reflecting the material's increasing resistance to further yielding. Isotropic hardening expands the surface uniformly in all directions, preserving its shape while raising the yield stress, as observed in many polycrystalline metals under monotonic loading. Kinematic hardening, conversely, translates the surface without changing its size, modeling the where reverse yielding occurs at lower stresses, common in cyclic loading scenarios. Combined isotropic-kinematic models account for both, providing accurate simulations for complex applications like metal forming and , where the surface's shape—often smooth but potentially featuring vertices—directly influences the material's anisotropic response and failure modes.

Introduction

Definition and Concept

In the theory of elastoplasticity, the yield surface delineates the transition from reversible elastic deformation to irreversible plastic deformation in materials subjected to multiaxial stress states. It is represented as a five-dimensional hypersurface embedded in the six-dimensional space of the stress tensor components, defining the locus of stress states at which yielding commences. This hypersurface encapsulates the material's capacity to withstand loads without permanent change, with its geometry reflecting the onset of plasticity under combined normal and shear stresses. Stress points interior to the yield surface produce exclusively responses, allowing full recovery of the original shape upon unloading. Points on the surface mark the initiation of yielding, where permanent deformation begins alongside . Exterior points are inadmissible in ideal models, as they would necessitate flow rates. The yield surface is characteristically , ensuring a unique direction for increment via the normality rule and maintaining thermodynamic consistency in the material's response. Common visualizations facilitate comprehension of this abstract geometry. Under plane stress assumptions, the yield surface reduces to a closed curve, or yield locus, in the two-dimensional plane of principal stresses σ₁ and σ₂. In full three-dimensional analyses, its projection onto the deviatoric plane—often the octahedral plane perpendicular to the hydrostatic axis—reveals the influence of shear components on yielding. The Haigh-Westergaard stress space offers a for the full five-dimensional surface, separating hydrostatic effects from deviatoric distortions. Yielding itself constitutes a form of instability, wherein a critical combination triggers the onset of permanent deformation, marking the limit of stability. This underpins elastoplasticity, the constitutive framework for materials displaying both recoverable behavior and non-recoverable plastic flow under sufficient loading. The yield surface's form is frequently expressed through invariants, which succinctly capture its rotational invariance and dependence on volumetric and distortional measures.

Historical Development

The development of yield surface theory began in the with empirical observations on metal deformation. In 1864, Henri Tresca proposed a yield criterion based on the maximum , derived from experiments on and processes, marking the first formal description of yielding as a limit state in solids. This idea was further explored in the 1870s by Adhémar Jean Claude Barré de Saint-Venant, who applied Tresca's maximum concept to analyze flow in metals under torsion and extension, emphasizing its role in distinguishing from regimes. Advancements in the early shifted toward energy-based interpretations for isotropic metals. In 1913, introduced the distortion energy hypothesis, positing that yielding occurs when the deviatoric reaches a equivalent to uniaxial tension, providing a more accurate prediction for ductile materials under multiaxial loading. Heinrich Hencky clarified this in 1924 by formulating the equivalent in terms of the second of the deviatoric tensor, offering a physical link to octahedral and solidifying the criterion's theoretical foundation. These works facilitated a transition from purely empirical maximum rules to invariant-based descriptions invariant under coordinate rotations. For geomaterials like soils and rocks, pressure-dependent criteria emerged earlier. laid the groundwork in 1773 with a friction-based failure model for granular media, expressing shear resistance as a of plus , initially applied to retaining walls and landslides. Otto Mohr extended this in 1900 by introducing a graphical envelope of Mohr's circles to represent the nonlinear failure boundary in shear- space, enabling better visualization of triaxial test data. In the 1950s, William Prager, collaborating with Daniel Drucker, generalized the Mohr-Coulomb model into a smooth conical yield surface in principal space, incorporating pressure sensitivity for soils while approximating von Mises behavior under low confinement. Concrete-specific yield surfaces addressed triaxial compression and tension differences in the mid-20th century. In 1958, Borislav Bresler and Karl S. Pister developed a using five independent strength parameters from uniaxial, biaxial, and triaxial tests, capturing the material's asymmetric response under multiaxial loading. This was refined in 1974 by Kurt J. Willam and Edwin P. Warnke, who proposed a three-parameter smooth surface in deviatoric plane, calibrated to match experimental meridians in triaxial compression and tension for normal-weight . Anisotropic extensions recognized directional strength variations in processed s. In the , Ludwik Burzyński introduced a of effort, decomposing into longitudinal, transverse, and components weighted by direction-dependent limits, applicable to rolled metals with differing tensile and compressive strengths. More recently, in the , Frédéric Barlat and colleagues advanced this with the Yld2004-18p criterion, a plane-stress orthotropic model using linear stress transformations and 18 parameters to accurately predict loci in sheet metals from simulations and r-value data. A pivotal event in the 1970s was the integration of yield surfaces into finite element methods for computational plasticity, enabling numerical simulation of large-scale elasto-plastic boundary value problems in engineering structures.

Theoretical Foundations

Representation in Stress Space

The stress state in a continuum is represented by the symmetric \boldsymbol{\sigma}, which possesses six independent components in : the normal stresses \sigma_{11}, \sigma_{22}, \sigma_{33} and the shear stresses \sigma_{12}, \sigma_{23}, \sigma_{31}. This defines a six-dimensional stress space, but for yield surface analysis, the representation is typically reduced to the three principal stresses \sigma_1 \geq \sigma_2 \geq \sigma_3 by diagonalizing the tensor, yielding a three-dimensional principal stress space where the axes align with the principal directions. This simplification eliminates shear components and focuses on the eigenvalues of \boldsymbol{\sigma}, enabling geometric interpretation of yielding without loss of generality for isotropic materials. To facilitate visualization and analysis, the principal stress space is often mapped into Haigh-Westergaard coordinates, a cylindrical system that decouples hydrostatic and deviatoric effects. The coordinate along the hydrostatic axis is \xi = I_1 / \sqrt{3}, where I_1 = \sigma_1 + \sigma_2 + \sigma_3 is the first stress invariant representing volumetric stress. The radial coordinate in the deviatoric plane is \rho = \sqrt{2 J_2}, with J_2 as the second invariant of the deviatoric stress tensor capturing shear distortion. The angular position in this plane is given by the Lode angle \theta, which ranges from $0to\pi/3 and distinguishes triaxial tension (\theta=0), compression (\theta=\pi/3$), and shear states in the octahedral plane. This parameterization highlights the cylindrical symmetry of many yield surfaces around the hydrostatic axis. Yield surfaces in this space are visualized through projections that reveal their shape and orientation. The \pi-plane, perpendicular to the hydrostatic axis at constant \xi, projects the deviatoric section as a closed , illustrating the dependence of strength on the angle for fixed hydrostatic stress. Meridional planes, which include the hydrostatic axis and a specific \theta (e.g., \theta = [0](/page/0) for axisymmetric ), display the surface's cross-section to show pressure-shear interactions, such as conical or parabolic forms in pressure-sensitive materials. These techniques aid in comparing experimental data with theoretical predictions across loading paths. The general yield condition is formulated as f(\boldsymbol{\sigma}) = 0, where f defines the boundary between and regimes in space. To exploit the separation of volumetric and distortional behaviors, \boldsymbol{\sigma} is decomposed into hydrostatic pressure p = -I_1 / 3 (positive in ) and deviatoric tensor \mathbf{s} = \boldsymbol{\sigma} + p \mathbf{I}, yielding f(p, \mathbf{s}) = 0. This additive decomposition, rooted in the properties of \boldsymbol{\sigma}, allows functions to incorporate mean effects independently of . Convexity of the yield surface f(\boldsymbol{\sigma}) \leq 0 is a fundamental requirement in rate-independent , ensuring a unique closest point projection during plastic loading and compliance with the maximum , which maximizes plastic work for irreversible deformation. This , derived from Drucker's stability postulate, implies that the normal to the surface defines the direction of plastic strain increments via the flow rule. Non-convex surfaces can permit ambiguous tangent hyperplanes, leading to unstable flow, multiple solutions in return mapping algorithms, and shear band localization in numerical simulations.

Invariants and Yield Functions

In , the yield surface is mathematically described using scalar invariants of the tensor, which provide a coordinate-independent representation of the state. These invariants allow the formulation of yield functions that define the boundary between and deformation in a manner consistent with physical principles. The primary invariants used are the first invariant I_1 and the second and third deviatoric invariants J_2 and J_3, derived from the \sigma. For isotropic materials, the function depends solely on these invariants, ensuring the criterion is independent of the . The first invariant I_1 captures the hydrostatic component of the state and is defined as the of the tensor: I_1 = \operatorname{tr}(\sigma) = \sigma_{kk} = \sigma_1 + \sigma_2 + \sigma_3, where \sigma_1, \sigma_2, \sigma_3 are the principal es. This represents the volumetric or hydrostatic pressure p = -I_1 / 3 (positive in compression), which influences yielding in pressure-sensitive materials like soils and rocks but is often negligible for metals. The deviatoric tensor s is then obtained by subtracting the hydrostatic part: s = \sigma - (-p) \mathbf{I} = \boldsymbol{\sigma} + p \mathbf{I}, where \mathbf{I} is the tensor. The second deviatoric J_2 measures the intensity of distortion and is given by J_2 = \frac{1}{2} s_{ij} s_{ij} = \frac{1}{6} \left[ (\sigma_1 - \sigma_2)^2 + (\sigma_2 - \sigma_3)^2 + (\sigma_3 - \sigma_1)^2 \right], with J_2 \geq 0. A common is the octahedral \tau_{\mathrm{oct}} = \sqrt{(2/3) J_2}, which relates to the distortion energy in yielding. The third deviatoric invariant J_3 accounts for the directional in the deviatoric plane and is defined as the of the deviatoric tensor: J_3 = \det(s). It introduces dependence on the Lode angle \theta, related to the Lode parameter \mu = (2 \sigma_2 - \sigma_1 - \sigma_3)/(\sigma_1 - \sigma_3), where -1 \leq \mu \leq 1. The Lode angle influences tension-compression , with \theta = 0 corresponding to uniaxial tension (\mu = -1) and \theta = \pi/3 to uniaxial compression (\mu = 1); this effect is pronounced in materials exhibiting different strengths under triaxial loading. For isotropic hardening , the function takes the general form f(\sigma, k) = f(I_1, J_2, J_3, k) = 0, where k is a strength parameter that evolves with equivalent \varepsilon^p, such as k = k(\varepsilon^p). Many criteria, like the von Mises function, depend only on J_2, but inclusion of J_3 allows for more accurate modeling of shear-induced . These invariants ensure frame-indifference, as they are objective scalars unchanged under rotations, and respect material symmetry groups, such as , by transforming appropriately under symmetry operations. This objectivity is fundamental to constitutive modeling in , guaranteeing that the yield surface description is independent of the observer's frame.

Isotropic Yield Criteria for Metals

Tresca Yield Criterion

The Tresca yield criterion, proposed by French engineer Henri Tresca in to analyze the flow of solids in machine components under small forces, posits that plastic yielding in ductile materials initiates when the maximum reaches a critical value equal to half the uniaxial yield stress. This maximum shear stress theory assumes that failure is governed by shear rather than normal stresses, reflecting the underlying mechanism of dislocation motion in metals. The formulation of the criterion in terms of principal stresses \sigma_1 \geq \sigma_2 \geq \sigma_3 is given by yielding occurring when the maximum principal equals k = \sigma_y / 2, where \sigma_y is the yield strength in uniaxial : \frac{1}{2} \max(|\sigma_1 - \sigma_2|, |\sigma_2 - \sigma_3|, |\sigma_3 - \sigma_1|) = k Equivalently, this simplifies to \sigma_{\max} - \sigma_{\min} = \sigma_y for the onset of yield under . In the space of principal stresses, the yield surface forms an infinite coaxial with the hydrostatic pressure axis (where \sigma_1 = \sigma_2 = \sigma_3), exhibiting a regular hexagonal cross-section in the deviatoric plane that remains invariant under hydrostatic loading. This geometry arises directly from the piecewise linear nature of the maximum difference between principal stresses, resulting in six flat faces corresponding to the pairs of principal stress differences. The simplicity of the Tresca criterion facilitates analytical solutions in design and closely matches experimental observations for and uniaxial in ductile metals, where dominates deformation. It provides a conservative estimate in many practical loading scenarios, aligning well with the insensitivity of metal yielding to hydrostatic . Despite these strengths, the criterion's sharp corners in the deviatoric plane lead to discontinuities, causing non-unique normal vectors and thus ambiguous plastic flow directions under associated flow rules at those points. Furthermore, it overpredicts the strength in certain biaxial states relative to experimental for ductile metals, as the faceted extends beyond observed yield loci in those orientations. Compared to the von Mises criterion, which yields a smoother elliptical cross-section, the Tresca surface's angular facets can complicate numerical implementations.

von Mises Yield Criterion

The von Mises yield criterion, also known as the maximum distortion energy criterion, posits that yielding in isotropic ductile materials occurs when the elastic distortion energy reaches the same value as that at yield in uniaxial tension. This hypothesis was originally proposed by in 1913, deriving the criterion from the decomposition of total strain energy into dilatational (volumetric) and distortional (shear) components, with yielding governed solely by the latter. The formulation of the is expressed in terms of the second deviatoric invariant J_2, where yielding initiates when \sqrt{3 J_2} = \sigma_y, with \sigma_y denoting the uniaxial . Equivalently, in principal coordinates \sigma_1, \sigma_2, \sigma_3, the takes the form \frac{1}{\sqrt{2}} \sqrt{ (\sigma_1 - \sigma_2)^2 + (\sigma_2 - \sigma_3)^2 + (\sigma_3 - \sigma_1)^2 } = \sigma_y. This defines the effective von Mises \sigma_{VM} = \sqrt{3 J_2}, which must not exceed \sigma_y for elastic behavior. Geometrically, the von Mises yield surface in principal stress space forms a right circular cylinder with its axis aligned along the hydrostatic stress direction (the line \sigma_1 = \sigma_2 = \sigma_3), reflecting independence from hydrostatic pressure. In the deviatoric plane (π-plane), perpendicular to the hydrostatic axis, the cross-section is a circle of radius \sigma_y / \sqrt{3}, making the basic form insensitive to the third deviatoric invariant J_3. This cylindrical shape ensures a smooth, convex surface that bounds allowable stress states. The criterion offers advantages in its differentiability, which facilitates unique associated flow rules in models, and its close alignment with experimental data for ductile metals under various multiaxial states, outperforming prismatic alternatives in most cases. However, it neglects hydrostatic pressure influences on yielding, assuming equal strengths in and , which limits applicability to materials sensitive to mean , such as certain polymers or rocks. Extensions to the basic isotropic form include kinematic hardening, where the cylindrical yield surface maintains a constant radius but translates along the space in the direction of plastic straining, accounting for the in cyclic loading.

Huber Criterion

The criterion, proposed by Polish engineer Maksymilian Tytus in 1904, establishes yielding in ductile materials when the octahedral \tau_\mathrm{oct} reaches a critical value of \frac{\sqrt{2}}{3} \sigma_y, where \sigma_y is the uniaxial yield . This formulation arises from considering the distortional as the measure of material effort, specifically tying to shear on octahedral planes. Mathematically, it corresponds to the second deviatoric invariant satisfying J_2 = \frac{1}{3} \sigma_y^2, marking the onset of plastic deformation under multiaxial loading. Although later recognized as equivalent to the von Mises-Hencky criterion due to identical predictive outcomes, the Huber criterion retains distinct recognition in for its early emphasis on shear-based mechanisms. Huber's work predated von Mises' publication and Hencky's 1924 contributions, yet it was independently developed and applied in Eastern European contexts, often cited separately to highlight its foundational role in energy-based yield hypotheses. In principal stress space, the yield surface defined by the Huber forms a right circular aligned with the hydrostatic , identical in geometry to the von Mises surface, with a of \sigma_y / \sqrt{3}. Its projection onto the octahedral —a plane normal to the direction—appears as a , reflecting the isotropic of the and simplifying of distortion under complex stresses. The octahedral planes correspond to the {111} planes in cubic crystals, which serve as primary slip planes in face-centered cubic metals like aluminum and , providing a physical rationale for the macroscopic distortion energy limit. The criterion's strength is its ability to link microscopic slip processes in metallic crystals to observable bulk yielding, offering predictive accuracy for ductile materials under proportional loading without needing detailed microstructural data. However, it shares the von Mises limitations, assuming no dependence on hydrostatic pressure, which restricts its applicability to pressure-sensitive materials like soils or rocks.

Pressure-Dependent Yield Criteria

Mohr–Coulomb Yield Surface

The Mohr–Coulomb yield criterion describes in frictional materials such as soils and rocks when the on a potential reaches a limiting value that increases linearly with the normal on that . The basic formulation is given by |\tau| = c + \sigma \tan \phi, where \tau is the , \sigma is the normal (positive in ), c is the , and \phi is the internal friction angle. Yielding occurs when the Mohr circle representing the stress state becomes tangent to this linear envelope, identifying the most critical . In terms of principal stresses (with \sigma_I \geq \sigma_{II} \geq \sigma_{III}, positive), the criterion for the dominant pair is \sigma_I - \sigma_{III} = (\sigma_I + \sigma_{III}) \sin \phi + 2c \cos \phi. In three dimensions, the full yield surface is defined by the maximum over all three principal stress pairs, ensuring the most severe condition governs . This criterion derives from Mohr's graphical representation of stress states using circles in the \sigma-\tau plane, building on Coulomb's 1776 hypothesis of frictional resistance along discrete failure planes. Mohr's 1900 analysis showed that shear failure initiates on planes oriented at an angle of $45^\circ + \phi/2 to the major principal direction, where the Mohr circle is tangent to the \tau = c + \sigma \tan \phi. This geometric construction captures the influence of confining pressure on , distinguishing it from pressure-insensitive criteria for metals. The associated flow rule implies at angle \psi = \phi, reflecting volume increase during in dense frictional materials. In principal stress space, the Mohr–Coulomb yield surface forms an open-ended irregular , with its axis aligned along the hydrostatic compression direction (isotropic stress line). The cross-section in the deviatoric (π) plane is an irregular , where the distance from the center varies with the angle, being smaller in triaxial extension than in triaxial compression due to the pressure sensitivity. This asymmetry highlights the criterion's ability to model differing tensile and compressive strengths, with uniaxial \sigma_c = 2c \cos \phi / (1 - \sin \phi) typically much larger than tensile strength. The criterion excels in geomechanics for capturing pressure-dependent yielding, tension-compression asymmetry, and associated in granular media like soils and rocks, making it a for stability analyses in and . However, its non-smooth corners at the hexagonal vertices pose numerical challenges in finite element simulations, often requiring regularization. Additionally, the lack of dependence on the third J_3 (intermediate principal effect) leads to overestimation of strength under certain triaxial states, as it assumes equal influence from all deviatoric paths. Variants with a tension cutoff \sigma_{III} = -T (where T = 2c \tan(45^\circ - \phi/2) is the positive tensile strength magnitude) cap the in , rendering the tensile limit independent of mean while preserving the open compressive apex.

Drucker–Prager Yield Surface

The Drucker–Prager yield surface serves as a , conical approximation to the Mohr–Coulomb yield criterion, enabling efficient numerical implementations in simulations of pressure-dependent materials such as soils and rocks. Developed by William Prager in 1952 as an extension of Daniel C. Drucker's stability postulate, the criterion leverages stress invariants to describe yielding in a form amenable to finite element methods (FEM). The yield function is expressed as f = \sqrt{J_2} + \alpha I_1 - k = 0, where I_1 denotes the first stress invariant, J_2 the second deviatoric stress invariant, and the material parameters \alpha and k are calibrated to Mohr–Coulomb cohesion c and friction angle \phi: \alpha = \frac{2 \sin \phi}{\sqrt{3} (3 - \sin \phi)}, \quad k = \frac{6 c \cos \phi}{\sqrt{3} (3 - \sin \phi)}. This formulation aligns the conical surface with the Mohr–Coulomb hexagon in triaxial compression states. Geometrically, the surface appears as a right circular in principal stress space, with its apex at I_1 = k / \alpha along the hydrostatic stress line \sigma_1 = \sigma_2 = \sigma_3. Depending on the calibration, the may be circumscribed (outer tangent) or inscribed (inner tangent) to the Mohr–Coulomb hexagonal pyramid, bounding the yield region conservatively. Advantages of this criterion include its full differentiability, supporting associated flow rules for plastic strain computation without issues at non-smooth vertices, and its incorporation of hydrostatic pressure effects in a simple invariant-based framework ideal for FEM. Limitations stem from the assumed circular deviatoric cross-section, which implies uniform and can lead to overprediction (circumscribed case) or underprediction (inscribed case) of shear capacity relative to the faceted Mohr–Coulomb surface in certain stress paths. For two-dimensional simulations under plane strain conditions, variants adjust the parameters to mimic three-dimensional behavior; for instance, the effective is modified as \tan \beta = \frac{2 \sin \phi}{1 - \sin \phi} in non-dilatant flow to match plane strain yielding.

Yield Criteria for Concrete and Geomaterials

Bresler–Pister Yield Surface

The Bresler–Pister yield surface provides a failure envelope for under multiaxial states, derived empirically from triaxial and tension-compression tests on plain specimens. Developed by Boris Bresler and Karl S. Pister, the criterion extends the Rankine maximum principal by introducing interactions among the principal stresses to better fit experimental on combined loading. This approach accounts for 's pronounced difference in tensile and compressive capacities, where failure is often governed by tensile cracking in brittle materials like . The original formulation emerged from analyses showing that simple maximum criteria underpredicted strengths in biaxial and triaxial regimes. In principal stress space, with \sigma_1 as the algebraically largest principal stress (typically in tension) and \sigma_2, \sigma_3 the others, the yield function is expressed as: f(\sigma_1, \sigma_2, \sigma_3) = \frac{\sigma_1^2}{f_t^2} + \frac{\sigma_2^2 + \sigma_3^2}{f_c^2} - 2\nu \frac{(\sigma_1\sigma_2 + \sigma_2\sigma_3 + \sigma_3\sigma_1)}{f_t f_c} where yield initiates when f = 1, f_t is the uniaxial tensile strength, f_c is the uniaxial compressive strength (with f_c \gg f_t), and \nu is Poisson's ratio (typically 0.2 for concrete). This equation, adapted from the failure condition for use in associated plasticity models, defines an ellipsoidal surface tilted relative to the hydrostatic axis, opening wider in the compressive octant due to higher f_c. The ellipsoid's asymmetry highlights concrete's tension-weak behavior without relying on explicit hydrostatic pressure terms, distinguishing it from frictional criteria. The criterion's advantages lie in its straightforward quadratic structure, which facilitates analytical solutions and numerical in finite element analyses for design. It effectively models the tension-compression asymmetry using only three material parameters, making it suitable for early computational tools lacking advanced formulations. Notably, it was incorporated into initial design codes and standards for predicting ultimate strengths under combined stresses. However, limitations include its neglect of intermediate principal stress influences, leading to inaccuracies in triaxial extension or shear-dominated states, and its principal stress-based form, which is not fully and requires careful axis alignment. These shortcomings prompted evolution toward more sophisticated models incorporating third stress effects.

Willam–Warnke Yield Surface

The Willam–Warnke yield surface is a three-parameter model designed to predict the onset of yielding in and similar cohesive-frictional materials under triaxial loading, incorporating dependence on the angle to account for the influence of the third deviatoric invariant. This criterion builds briefly on approaches like Bresler–Pister by adding angle effects for a more accurate deviatoric shape. Developed specifically for , it separates the yield into meridional (pressure-dependent) and deviatoric (-dependent) components, enabling better representation of failure modes in , , and . The model was introduced by Karl J. Willam and Ernest P. Warnke in their 1974 work, derived from experimental triaxial test data on specimens to fit both the / meridians and the deviatoric cross-sections. The is expressed in Haigh-Westergaard stress space, where the hydrostatic axis is given by \xi = I_1 / \sqrt{3} (with I_1 the first stress invariant), the deviatoric radius by \rho = \sqrt{2 J_2} (with J_2 the second deviatoric invariant), and the Lode angle \theta ranging from 0° () to 60° (). The yield function takes the form f = \sqrt{a_1 \rho^2 + a_2 \rho + a_3} + m b_1 I_1 - b_2 = 0, where the coefficients a_1, a_2, a_3 define the shape of the deviatoric section as a function of \theta and pressure, while b_1, b_2 and the multiplier m (often related to compressive strength) are calibrated to match uniaxial tensile/compressive strengths and biaxial/triaxial data. This polynomial form for the deviatoric part ensures smoothness and convexity, avoiding sharp corners in the failure envelope. Geometrically, the yield surface manifests as a rounded triangular pyramid in principal stress space, opening irregularly with increasing hydrostatic pressure; the deviatoric cross-section is a smooth, nearly triangular curve that transitions from a more rounded shape in tension to sharper edges in compression, reflecting concrete's differential strength in these regimes. Key advantages of the Willam–Warnke criterion include its ability to capture the nonlinear curvature of the compression and the asymmetry due to the Lode angle, providing improved predictions for high-strength under multiaxial states compared to simpler isotropic models. It is particularly effective for simulating shear failure and hydrostatic sensitivity in structural applications. However, extending to a five-parameter version for greater flexibility in meridian shapes increases complexity, and the model inherently assumes fixed ratios between tensile and compressive meridian curvatures, limiting adaptability to highly variable material data. The criterion has been widely adopted in commercial finite element codes, such as , where it underpins the concrete damaged plasticity model for nonlinear of elements under dynamic and static loads.

Podgórski and Rosendahl Trigonometric Yield Surfaces

The trigonometric yield surfaces developed by Podgórski and Rosendahl provide smooth, pressure-sensitive models for the yielding of geomaterials such as and rock, incorporating to capture the influence of the intermediate principal without singularities. These criteria extend traditional linear models like Mohr-Coulomb or Drucker-Prager by blending the effects of principal stresses through sinusoidal terms, resulting in closed, surfaces that better approximate experimental from triaxial tests.111:2(188)) Podgórski's formulation, introduced in the 1980s, derives from failure loci and generalizes to three dimensions for isotropic media. The yield function is expressed in terms of ordered principal stresses σ₁ ≥ σ₂ ≥ σ₃ as: f = \frac{(\sigma_1 - \sigma_3) \left[1 + \sin\left(\frac{\pi}{2} \cdot \frac{(\sigma_2 - \sigma_3)}{(\sigma_1 - \sigma_3)}\right)\right]}{2k} + \alpha (\sigma_1 + \sigma_3) = 0 where k represents the shear yield strength and α is a sensitivity related to the material's friction angle. This sine term smoothly interpolates the contribution of the intermediate stress σ₂, transitioning the deviatoric section from a triangular at low pressures (resembling Tresca) to more rounded forms at higher pressures, avoiding sharp corners that can cause numerical issues in simulations. The derivation starts from biaxial and uniaxial test data to fit the locus, then extends to using stress invariants for convexity assurance.111:2(188)) The Rosendahl variant, developed in the as an extension for 3D yielding in and , modifies the Drucker-Prager with a trigonometric adjustment emphasizing the θ through a cos(3θ) term. This allows the deviatoric plane to vary smoothly from triangular to nearly circular shapes with increasing hydrostatic , capturing axisymmetric in compression-dominated regimes. The features rounded edges and a closed surface in principal stress space, ensuring differentiability for associated rules.111:2(188)) These trigonometric surfaces offer advantages over linear models by eliminating singularities at transitions and providing superior fits to triaxial experimental for geomaterials, where intermediate effects are pronounced. For instance, they predict higher strengths in balanced triaxial compared to Mohr-Coulomb, aligning with observed rock behavior under moderate confinement. However, parameter calibration remains challenging, requiring multiple tests to determine k and α accurately, and these criteria are less commonly implemented in commercial finite element software due to their relative complexity.

Anisotropic and Advanced Yield Surfaces

Burzyński-Yagn Criterion

The Burzyński-Yagn criterion is an early pressure-sensitive yield model that can be extended to anisotropic materials exhibiting direction-dependent strength, often due to processes like rolling. It originates from an energy-based hypothesis considering both distortional and volumetric energies, generalized for by incorporating directional tensile, compressive, and strengths. For the isotropic case, the yield surface is a second-order about the hydrostatic axis, formulated using as $3I_2' = \frac{\sigma_{\text{eq}} - \gamma_1 I_1}{1 - \gamma_1} \frac{\sigma_{\text{eq}} - \gamma_2 I_1}{1 - \gamma_2}, where I_2' is the second deviatoric , \sigma_{\text{eq}} the equivalent , I_1 the first , and \gamma_1, \gamma_2 material parameters controlling tension-compression asymmetry and pressure sensitivity. Anisotropic extensions modify this by direction-dependent coefficients, capturing orthotropy in polycrystalline metals through calibration with uniaxial and tests along principal axes. Proposed by W. Burzyński in , the criterion draws from experimental observations of pressure effects on yielding, deriving limits from uniaxial tension, , and . It has been applied and modified in later works for anisotropic pressure-dependent solids. In principal stress space, the yield surface appears as a distorted surface (e.g., or ), with facets reflecting interactions between normal and yielding in orthotropic directions, highlighting reduced compared to isotropic criteria like von Mises. The criterion's advantages include its ability to capture initial in rolled sheet metals using directional strength measurements and explicit -compression , valuable for textured alloys. It offers a for preliminary where orthotropic effects dominate. However, the quadratic nature can lead to non-smoothness in anisotropic extensions, causing discontinuities in the flow rule at certain points that complicate finite element simulations. It may also inadequately predict yielding under high or equi-biaxial due to its energy-based boundaries. As a foundational model for pressure-sensitive , it influenced subsequent criteria.

Bigoni–Piccolroaz Yield Surface

The is a for pressure-sensitive materials, offering versatility through angle dependence and extensions to . Proposed in isotropic form by Bigoni and Piccolroaz in 2004, it has been generalized to include anisotropic effects and surfaces with corners, suitable for materials with directional variations. The model separates hydrostatic and deviatoric effects, controlling tension-compression asymmetry and shear response. The yield condition is f(p, q, \theta) = F(p) + \frac{q}{g(\theta)} = 0, where p = -\operatorname{tr}(\sigma)/3 is hydrostatic pressure, q = \sqrt{3 J_2} the von Mises stress, and \theta the Lode angle. The meridian function is F(p) = -M p_c \sqrt{(\phi - \phi_m)[2(1-\alpha)\phi + \alpha]} for \phi = (p + c)/(p_c + c) \in [0,1], and infinite otherwise; the deviatoric function is g(\theta) = 1 / \cos[\beta \pi/6 - (1/3) \cos^{-1}(\gamma \cos 3\theta)]. Parameters include M (friction angle), p_c, c (compression/tension strengths), \alpha, m (meridian shape), \beta, \gamma (deviatoric shape). Anisotropic versions vary parameters across Lode angle sectors or incorporate directional tensors. This ensures a smooth, convex surface for stable simulations. In the deviatoric plane, the section is rounded triangular, with parameters tuning uniaxial / and ratios (e.g., higher \beta toward Tresca). Hydrostatic affects scaling via F(p), capturing strengthening. For meridians, it uses , , and interpolating curves for smooth transitions. The anisotropic extension builds on the trigonometric deviatoric function, adding orthotropy via variable coefficients. Advantages include flexibility to recover criteria like Tresca, Mohr-Coulomb, or von Mises as limits, while ensuring smoothness and convexity via parameter constraints. It supports hardening and applies to foams, soils, and metals with Lode effects. However, up to seven (or more for anisotropic) parameters require extensive tests like triaxial across Lode angles, limiting routine use compared to Drucker-Prager.

Cosine Ansatz (Altenbach-Bolchoun-Kolupaev)

The cosine , developed by Altenbach, Bolchoun, and Kolupaev (2011–2015), is a yield criterion for orthotropic materials like textured composites and polycrystals under complex loading. It uses to model directional yield variations, extending forms for better representation. It builds on early anisotropic ideas, incorporating periodic functions for material symmetries. The function uses invariants: (3I_2')^3 \cdot \frac{1 + c_3 \cos 3\theta + c_6 \cos^2 3\theta}{1 + c_3 + c_6} = \left[ \frac{\sigma_{\rm eq} - \gamma_1 I_1}{1 - \gamma_1} \right]^{6-l-m} \left[ \frac{\sigma_{\rm eq} - \gamma_2 I_1}{1 - \gamma_2} \right]^l \sigma_{\rm eq}^m, where \cos 3\theta = \frac{3\sqrt{3}}{2} \frac{I_3'}{I_2'^{3/2}}, c_3, c_6 shape the deviatoric (with convexity constraints like c_6 \geq (5/12) c_3^2 - 1/3), \gamma_1, \gamma_2 hydrostatic nodes, and l, m powers for meridian curvature. \sigma_{\rm eq} scales to , calibrated from orthotropic tests. This leverages three-fold for orthotropy. Geometrically, it yields a closed, undulating surface in principal stress space, smooth in the deviatoric plane without sharp corners, promoting convex envelopes for flow rules. The cosine terms approximate , ensuring convexity for suitable parameters. Advantages include guaranteed convexity for in finite element analysis, few parameters (c_3, c_6, \gamma_1, \gamma_2, l, m, \sigma_y), efficient orthotropy modeling without heavy , suiting composites. It extends Hill's , improving predictions. Limitations: Primarily for plane stress, less applicable to full 3D with out-of-plane effects; focuses on deviatoric anisotropy over strong pressure dependence, potentially needing combination with Drucker-Prager for geomaterials.

Barlat's Yield Surface

Barlat's yield surfaces are a family of anisotropic criteria for sheet metals, especially in plane stress forming. Evolving from Yld89 (plane stress orthotropy), Yld91 incorporated six parameters from uniaxial tests and r-values for directional yield and strain ratios. Later, Yld2000 enhanced aluminum flexibility, and Yld2004-18p extended to 3D with 18 parameters from sheet tests. The Yld91 function is \phi = |\phi_1' - \phi_2'|^a + |\phi_2' - \phi_3'|^a + |\phi_3' - \phi_1'|^a = 2k, where \phi_i' are principals of one linear of the deviatoric (with eight \alpha_k coefficients for orthotropy), and a = 6 or $8 (crystal-dependent); k relates to shear yield. A second is used in some variants. In plane , it distorts the von Mises circle into an irregular ellipse capturing directional variations. Yld2004-18p uses two transformations for general stresses, yielding convex but irregular surfaces. These excel in simulating earing in for aluminum/steel sheets, reproducing yield loci from r-value data for better formability predictions. Convexity ensures finite element stability. The exponent a needs via experiments or polycrystal models, and variants increase computation. Recent extensions, such as a 2020 model for tension-compression , improve textured modeling. Barlat's models are standard in automotive simulations for defect prediction.