Stress to Strain 3D

Stress to Strain 3D

3D Stress to Strain Conversion Theory

Introduction

This tool converts three-dimensional stress states to corresponding strain states for isotropic materials. The tool performs both material-independent stress analysis (stress invariants, von Mises stress, principal stresses) and material-dependent strain calculations when elastic properties are provided. This is essential for structural analysis, experimental validation, and failure prediction in engineering applications.

Figure 1: Three-dimensional stress tensor showing six independent stress components

1. Three-Dimensional Stress State

The complete stress state at a point in a three-dimensional body is described by the stress tensor, a symmetric 3×3 matrix with six independent components:

• σ_xx = normal stress in the x-direction
• σ_yy = normal stress in the y-direction
• σ_zz = normal stress in the z-direction
• τ_xy = shear stress in the xy-plane (τ_xy = τ_yx by symmetry)
• τ_yz = shear stress in the yz-plane (τ_yz = τ_zy by symmetry)
• τ_zx = shear stress in the zx-plane (τ_zx = τ_xz by symmetry)

Tensor Notation: This tool uses double-subscript tensor notation throughout for clarity and consistency with continuum mechanics literature. Positive normal stresses are tensile; positive shear stresses follow the sign convention where they act in the positive coordinate direction on the positive face of an element.

2. Material-Independent Stress Analysis

Several important stress measures can be calculated directly from the stress tensor without requiring material properties. These are useful for evaluating stress states, comparing results, and applying failure theories.

2.1 Mean (Hydrostatic) Stress

The mean stress (also called hydrostatic stress) represents the average of the three normal stress components:

σ_m = (σ_xx + σ_yy + σ_zz) / 3

The mean stress is associated with volumetric deformation (expansion or contraction) and does not contribute to material yielding in most ductile metals (von Mises criterion).

2.2 Stress Invariants

Stress invariants are scalar quantities that remain constant under coordinate system rotation. They are fundamental properties of the stress tensor used to calculate principal stresses and evaluate failure criteria.

First stress invariant (I₁):

I₁ = σ_xx + σ_yy + σ_zz

I₁ is the trace of the stress tensor and equals three times the mean stress (I₁ = 3σ_m). Units: stress (Pa, PSI, etc.)

Second stress invariant (I₂):

I₂ = σ_xxσ_yy + σ_yyσ_zz + σ_zzσ_xx - τ_xy² - τ_yz² - τ_zx²

Units: stress² (Pa², PSI², etc.)

Third stress invariant (I₃):

I₃ = σ_xxσ_yyσ_zz + 2τ_xyτ_yzτ_zx - σ_xxτ_yz² - σ_yyτ_zx² - σ_zzτ_xy²

I₃ is the determinant of the stress tensor. Units: stress³ (Pa³, PSI³, etc.)

2.3 Von Mises Stress

The von Mises stress (also called equivalent stress) is a scalar measure of the stress state used in the von Mises yield criterion for ductile materials. It represents the stress intensity for distortion (shape change) independent of hydrostatic stress.

σ_vm = √[(σ_xx - σ_yy)² + (σ_yy - σ_zz)² + (σ_zz - σ_xx)² + 6(τ_xy² + τ_yz² + τ_zx²)] / √2

Yield Criterion: According to the von Mises criterion, yielding occurs when σ_vm ≥ σ_yield. This is one of the most widely used failure criteria for ductile metals.

2.4 Octahedral Stresses

Octahedral stresses act on the octahedral plane, which is equally inclined to all three principal axes. They are used in some failure theories and provide insight into the stress state.

Octahedral normal stress:

σ_oct = σ_m = (σ_xx + σ_yy + σ_zz) / 3

Octahedral shear stress:

τ_oct = (1/3)√(σ_xx - σ_yy)² + (σ_yy - σ_zz)² + (σ_zz - σ_xx)² + 6(τ_xy² + τ_yz² + τ_zx²)

The octahedral shear stress is related to von Mises stress by: τ_oct = (√2/3)σ_vm

2.5 Principal Stresses

Principal stresses (σ₁, σ₂, σ₃) are the eigenvalues of the stress tensor. They represent the normal stresses on planes where shear stress is zero, and are ordered such that σ₁ ≥ σ₂ ≥ σ₃.

Principal stresses are found by solving the characteristic equation:

λ³ - I₁λ² + I₂λ - I₃ = 0

This cubic equation is solved using Cardano's method with the trigonometric solution for real roots (the typical case for physical stress states):

Step 1: Calculate reduced coefficients:

p = I₂ - I₁²/3
q = -2I₁³/27 + I₁I₂/3 - I₃

Step 2: Calculate discriminant to verify real roots:

Δ = q²/4 + p³/27

For Δ ≤ 0 (real roots, typical case):

Step 3: Trigonometric solution:

m = 2√(-p/3)
θ = (1/3)arccos[3q/(pm)]

σ₁ = m cos(θ) + I₁/3 (maximum principal stress)

σ₂ = m cos(θ - 2π/3) + I₁/3 (intermediate principal stress)

σ₃ = m cos(θ - 4π/3) + I₁/3 (minimum principal stress)

2.6 Maximum Shear Stress

The maximum shear stress occurs on planes oriented 45° to the principal stress directions:

τ_max = (σ₁ - σ₃) / 2

This is used in the Tresca (maximum shear stress) yield criterion: yielding occurs when τ_max ≥ σ_yield/2.

3. Three-Dimensional Isotropic Constitutive Relations

When elastic material properties are provided, strains can be calculated from stresses using the generalized Hooke's law for isotropic materials. An isotropic material has the same elastic properties in all directions.

3.1 Material Properties

• E = elastic modulus (Young's modulus) - relates normal stress to normal strain
• ν = Poisson's ratio - relates transverse strain to axial strain (typically 0.25-0.35 for metals)
• G = shear modulus - relates shear stress to shear strain

For isotropic materials, the shear modulus is not independent:

G = E / [2(1 + ν)]

If G is not provided, this tool calculates it automatically from E and ν.

3.2 Normal Strain Components

The three normal strain components are calculated using the 3D isotropic stress-strain relations. Each normal strain depends on all three normal stresses due to the Poisson effect:

ε_xx = (1/E)[σ_xx - ν(σ_yy + σ_zz)]

ε_yy = (1/E)[σ_yy - ν(σ_xx + σ_zz)]

ε_zz = (1/E)[σ_zz - ν(σ_xx + σ_yy)]

Physical Interpretation: Each normal strain has three terms:

• Extension in the direction of the applied stress (σ/E)
• Contraction due to Poisson effect from perpendicular stresses (-νσ/E)

3.3 Shear Strain Components (Tensor Form)

Tensor shear strains (ε_ij) represent half the change in angle between initially perpendicular line elements. They are related to shear stresses by:

ε_xy = τ_xy / (2G)

ε_yz = τ_yz / (2G)

ε_zx = τ_zx / (2G)

Note: For isotropic materials, shear strains are independent of normal stresses (no coupling). This is a key feature of isotropy.

4. Engineering vs. Tensor Shear Strains

Two conventions exist for shear strain, differing by a factor of 2:

• Tensor shear strain (ε_ij): Used in tensor notation and theoretical mechanics. Represents half the change in angle.
• Engineering shear strain (γ_ij): Used in engineering practice and experimental measurements. Represents the total change in angle.

Relationship:

γ_xy = 2ε_xy = τ_xy / G

γ_yz = 2ε_yz = τ_yz / G

γ_zx = 2ε_zx = τ_zx / G

This tool calculates both forms internally but displays engineering shear strains in tables for consistency with experimental practice.

5. Volumetric Strain

The volumetric strain (ε_vol) represents the fractional change in volume of a material element. It is the trace of the strain tensor:

ε_vol = ε_xx + ε_yy + ε_zz

For isotropic materials, volumetric strain can also be expressed in terms of stresses:

ε_vol = (1 - 2ν)(σ_xx + σ_yy + σ_zz) / E = 3(1 - 2ν)σ_m / E

Bulk Modulus: The ratio of mean stress to volumetric strain is the bulk modulus K = E / [3(1 - 2ν)].

Incompressibility: Materials with ν = 0.5 are incompressible (ε_vol = 0 under any stress state). Rubber and many polymers approach this limit.

6. Microstrain Convention

Strains are dimensionless quantities (change in length per unit length), but are typically very small in engineering structures. For convenience, strains are often expressed in microstrain:

1 με = 1 × 10^-6 strain

This convention is standard in experimental stress analysis where strain gages output readings in microstrain. This tool provides strain results in both dimensionless form and microstrain for convenience.

Example: A strain of 0.001234 (dimensionless) equals 1,234 με.

7. Applications

This 3D stress-to-strain conversion tool has numerous applications in engineering analysis:

• Structural analysis: Converting finite element analysis (FEA) stress results to strains for comparison with allowable strain limits
• Experimental validation: Converting measured strains from strain gages to stresses for comparison with analytical predictions or FEA results
• Failure prediction: Calculating von Mises stress, octahedral shear stress, or maximum shear stress for use in yield criteria (von Mises, Tresca, octahedral shear stress theory)
• Material characterization: Verifying constitutive relations and material properties from multi-axial test data
• Design verification: Ensuring critical stress and strain measures remain within safe limits under combined loading
• Fatigue analysis: Determining principal stresses and strains for fatigue life prediction under complex loading

8. Calculation Procedure

This tool implements the following calculation sequence:

1. Input: Six stress components (σ_xx, σ_yy, σ_zz, τ_xy, τ_yz, τ_zx) and optional material properties (E, ν, G)
2. Material-independent stress analysis:
   • Calculate mean stress σ_m
   • Calculate stress invariants I₁, I₂, I₃
   • Calculate von Mises stress σ_vm
   • Calculate octahedral stresses σ_oct, τ_oct
   • Solve characteristic equation for principal stresses σ₁, σ₂, σ₃ using Cardano's method
   • Calculate maximum shear stress τ_max
3. Strain calculations (if E and ν are provided):
   • Calculate or verify shear modulus G = E/[2(1+ν)]
   • Calculate normal strains ε_xx, ε_yy, ε_zz using 3D Hooke's law
   • Calculate tensor shear strains ε_xy, ε_yz, ε_zx
   • Convert to engineering shear strains γ_xy, γ_yz, γ_zx
   • Calculate volumetric strain ε_vol
4. Output: Tables of stress measures, stress invariants, strain components (dimensionless), and strain components in microstrain

9. Summary

This comprehensive 3D stress-to-strain analysis tool provides both material-independent stress characterization and material-dependent strain calculations. The material-independent analysis (stress invariants, von Mises stress, principal stresses) is useful for evaluating stress states and applying failure criteria without requiring material properties. When elastic properties are provided, the tool calculates the complete 3D strain state using isotropic constitutive relations.

The theoretical framework is based on continuum mechanics, elasticity theory, and tensor analysis, following the approaches presented in classical references on experimental stress analysis and advanced mechanics of materials. All calculations use tensor notation for consistency with the engineering mechanics literature.