Technical Resource

Isotropic, Kinematic or Combined Hardening?

A practical guide to the Bauschinger effect, yield-surface evolution, cyclic plasticity and the choice of hardening model in nonlinear FEA.

Metal Plasticity Nonlinear FEA Cyclic Loading Material Modeling

Plasticity & Hardening

What happens after a metal first yields?

Have you ever wondered what the Bauschinger effect really means, how it is represented in finite element analysis, and why it can matter to an engineering prediction?

Once plastic deformation occurs, a metal does not necessarily respond identically when the loading direction is reversed. The evolution of the yield condition is commonly represented through isotropic hardening, kinematic hardening, or a combination of both.

Before going deeper into the equations, use the interactive demonstration below to see how the stress state, backstress, yield locus and hysteresis response evolve during cyclic loading.

Interactive Demonstration

Explore the hardening response

Apply cyclic strain and switch between capped isotropic, kinematic and combined hardening. During active plastic yielding, the current stress state lies on the active J₂ yield ellipse; during elastic loading and unloading it remains inside the yield locus.

Stress State

0.0 MPa

Backstress (α)

0.0 MPa

Material Status

Neutral / Idle

📋 Uniaxial Specimen Setup

Deviatoric Stress Space Plane (s₁-s₂)

Cyclic Stress-Strain Curve

🏆 Material Configuration

Young's Modulus E (MPa)
Initial Yield strength σy₀ (MPa)
Saturation Limit σsat (MPa)
Kinematic Modulus Hk (MPa)
Isotropic Modulus Hi (MPa)

⚙️ Hardening Formulation Mode

1️⃣ Pure Kinematic Capped Hardening 2️⃣ Pure Isotropic Capped Hardening 3️⃣ Combined Capped Hardening (Hk + Hi)

🛑 Drag Handle Left/Right To Apply Strain

Min ε (%)

Max ε (%)

Active Displacement (ε): 0.0000 %

📖 Principal Constitutive Formulations

Total Stress Function:
  σ = E · (ε − εᵖ)
   [E: Young's Modulus, εᵖ: Plastic Strain]

The Yield Criterion:
  f(σ, α) = |σ − α| − σ_y ≤ 0
   [The active hardening response is capped by the configured σ_sat limit]

Capped Kinematic Hardening Evolution:
  dα = H_k · dεᵖ
   [Driven by H_k Modulus; translates yield surface center]

Capped Isotropic Hardening Evolution:
  σ_y = σ_y₀ + H_i · ε̄ᵖ
   [Driven by H_i Modulus; expands radius via accumulated ε̄ᵖ]

🧠 Academic Theory: Uniaxial Deviatoric Stress Space Mapping

The deviatoric stress tensor isolates pure shape distortion by removing hydrostatic pressure components, defined as s_ij = σ_ij − ⅓ σ_kk δ_ij . Under a standard pure uniaxial loading state tensor configuration (σ, 0, 0) , the hydrostatic pressure simplifies directly to p = ⅓σ. Subtracting this from the main diagonal yields the principal deviatoric coordinates: s₁ = ⅔σ and s₂ = s₃ = −⅓σ .

Consequently, any standard uniaxial testing path is geometrically constrained to satisfy s₂ = −½ s₁. For J₂ / von Mises plasticity, with s₃ = −(s₁ + s₂), the yield locus in the ordinary s₁–s₂ Cartesian projection is (s₁−a₁)² + (s₁−a₁)(s₂−a₂) + (s₂−a₂)² = σ_y² / 3 , where a₁ = ⅔α and a₂ = −⅓α. Therefore the J₂ yield locus appears as a rotated ellipse in this s₁–s₂ plot. During active plastic yielding the current stress point lies on this ellipse; during elastic loading/unloading it lies inside it. The configured σ_sat is retained as the cap for the hardening response.

Isotropic Hardening

Expansion of the yield surface

In isotropic hardening, the center of the yield surface remains fixed while its size increases with accumulated plastic deformation.

σ_y = σ_y0 + H_i ε̄^p Simplified linear isotropic hardening representation

In a one-dimensional interpretation, further yielding moves farther away from zero stress in both tension and compression. This can be an effective engineering approximation when plastic loading is mainly monotonic and reverse yielding is not important.

Kinematic Hardening

Translation of the yield surface

Kinematic hardening introduces a backstress, α. Instead of simply growing, the yield surface translates in stress space as plastic deformation develops.

dα = H_k dε^p Simplified linear kinematic hardening representation

The corresponding one-dimensional yield condition can be written as:

f(σ,α) = |σ − α| − σ_y ≤ 0

Translation of the yield surface changes the stress required to produce yielding after a reversal in loading direction.

Bauschinger Effect

Why can reverse yielding occur earlier?

Consider a material with an initial yield stress of 200 MPa. After tensile plastic deformation, suppose a positive backstress of 20 MPa has developed.

Forward tensile yield: approximately +220 MPa

Reverse compressive yield: approximately −180 MPa

The material has therefore not simply become stronger in every direction. Its previous plastic deformation has shifted the yield surface. This reduction in reverse yield stress is the essential idea behind the Bauschinger effect.

Engineering Relevance

Why does this matter in FEA?

The hardening formulation becomes especially important when the loading path contains substantial unloading and reverse plasticity.

Application	Hardening consideration
Predominantly monotonic plastic loading	Isotropic hardening may often be sufficient.
Bending and springback	Reverse yielding can influence residual stresses and final geometry.
Bending-unbending over tooling	The material experiences repeated changes in plastic-flow direction.
Low-cycle fatigue	Cyclic plastic strain and hysteresis behavior depend strongly on the constitutive model.
Thermo-mechanical cycling	Repeated thermal expansion and constraint can produce alternating plastic deformation.
Ratcheting	Accumulation of plastic strain over repeated cycles requires an appropriate cyclic model.

Combined Hardening

Real materials can both translate and expand

Real metals do not necessarily behave as purely isotropic or purely kinematic systems. Combined hardening allows the yield surface to both translate and change size.

Isotropic contribution

Controls evolution of the yield-surface size with accumulated plastic deformation.

Kinematic contribution

Controls translation of the yield surface and directional memory of previous plastic loading.

Combined response

Provides greater flexibility for representing cyclic stress-strain behavior.

Advanced Constitutive Modeling

Beyond basic hardening laws

The simplified hardening laws above are useful for understanding the core mechanisms, but demanding engineering applications may require more advanced constitutive formulations.

Depending on the material, temperature, strain rate and loading history, the response may involve:

Nonlinear isotropic and kinematic hardening
Cyclic hardening and cyclic softening
Combined nonlinear plasticity formulations
Viscoplasticity and strain-rate-dependent behavior
Creep and high-temperature time-dependent deformation
Damage initiation and damage evolution
Ductile and brittle failure
Thermo-mechanical constitutive behavior
Custom constitutive models and user-material implementations

MakeCoTech can support advanced nonlinear material modeling from model selection and parameter calibration through implementation, verification and finite element application.

This includes plastic hardening and softening, cyclic plasticity, viscoplasticity, creep, damage and failure modeling for structural and thermo-mechanical simulations.

Model Calibration

The constitutive parameters must represent the loading path

A sophisticated constitutive law is useful only if its parameters are meaningful. A monotonic tensile test may be sufficient for a simple isotropic plasticity model, while cyclic models generally require additional loading-unloading or tension-compression data.

The objective is not to use the most complicated model available, but to select a model that captures the physical mechanisms that materially affect the engineering prediction.

Key Takeaways

The loading path should guide the hardening model

Isotropic hardening expands the yield surface.
Kinematic hardening translates the yield surface.
Kinematic translation provides a mechanism for representing the Bauschinger effect.
Combined hardening can represent both translation and expansion.
Cyclic applications may require nonlinear hardening formulations.
Viscoplasticity, creep, damage or failure may also be important depending on the application.
Material parameters should be calibrated against representative experimental data.

Advanced Engineering Support

Need a material model beyond standard plasticity?

MakeCoTech provides support in nonlinear finite element analysis, constitutive modeling, material calibration, cyclic plasticity, viscoplasticity, creep, damage and failure.

Discuss Your Application