Strain Transformation 2D

Strain Transformation 2D

Strain Transformation Theory

Introduction

Strain transformation is the process of converting strain components from one coordinate system to another coordinate system rotated by an angle. This is essential in engineering analysis when strain measurements or calculations are performed in one orientation, but results are needed in a different orientation, or when determining principal strains and their directions.

This tool transforms a known plane strain state (ε_x, ε_y, γ_xy) in the (x,y) coordinate system to the strain state (ε_x', ε_y', γ_x'y') in a rotated (x',y') coordinate system, and calculates principal strains and stresses.

Figure 1: Strain transformation from (x,y) coordinates to (x',y') coordinates rotated by angle θ

1. Plane Strain State

For a point in a loaded structure under plane stress conditions (free surface where σ_zz = τ_xz = τ_yz = 0), the complete state of strain is defined by three independent components:

• ε_x = normal strain in the x-direction (dimensionless, typically expressed in microstrain με)
• ε_y = normal strain in the y-direction (dimensionless, typically expressed in microstrain με)
• γ_xy = engineering shear strain in the xy-plane (dimensionless, typically expressed in microstrain με)

Note: Engineering shear strain γ_xy is used throughout this tool, which is related to tensor shear strain by γ_xy = 2ε_xy. Engineering shear strain represents the total change in angle between initially perpendicular line elements, consistent with standard mechanics of materials practice.

2. Strain Transformation Equations

When the coordinate system is rotated counterclockwise by an angle θ from the (x,y) axes to the (x',y') axes, the strain components transform according to the following equations:

ε_x' = ε_x cos²θ + ε_y sin²θ + γ_xy sinθ cosθ

ε_y' = ε_x sin²θ + ε_y cos²θ - γ_xy sinθ cosθ

γ_x'y' = -2(ε_x - ε_y) sinθ cosθ + γ_xy (cos²θ - sin²θ)

These equations can be derived using double-angle trigonometric identities. Alternatively, they can be expressed using the identities cos²θ = (1 + cos2θ)/2, sin²θ = (1 - cos2θ)/2, and sinθcosθ = sin2θ/2:

ε_x' = (ε_x + ε_y)/2 + (ε_x - ε_y)/2 cos2θ + (γ_xy/2) sin2θ
ε_y' = (ε_x + ε_y)/2 - (ε_x - ε_y)/2 cos2θ - (γ_xy/2) sin2θ
γ_x'y'/2 = -(ε_x - ε_y)/2 sin2θ + (γ_xy/2) cos2θ

Physical Interpretation:

• Normal strains (ε_x', ε_y') in the rotated axes depend on both normal and shear strains in the original axes
• Shear strain (γ_x'y') in the rotated axes similarly depends on the complete original strain state
• The transformation preserves strain invariants (quantities that remain constant under rotation)
• At certain angles θ, the shear strain γ_x'y' becomes zero - these are the principal strain directions

3. Principal Strains

The principal strains (ε₁, ε₂) are the maximum and minimum normal strains that occur at a point. They occur on perpendicular planes where shear strain is zero. Principal strains can be calculated directly from the original strain state without transformation using Mohr's Circle for strain.

Mohr's Circle Method:

Mohr's circle is a graphical and analytical tool for finding principal strains. The circle is constructed in (ε, γ/2) space with:

• Center: (ε_avg, 0) where ε_avg = (ε_x + ε_y) / 2
• Radius: R = √[(ε_x - ε_y)² / 4 + (γ_xy / 2)²]

The principal strains are located at the rightmost and leftmost points on the circle:

ε₁ = ε_avg + R (maximum principal strain)

ε₂ = ε_avg - R (minimum principal strain)

Expanded form:

ε_1,2 = (ε_x + ε_y) / 2 ± √[(ε_x - ε_y)² / 4 + (γ_xy / 2)²]

4. Principal Strain Angle

The angle φ from the x-axis to the ε₁ (maximum principal strain) direction is found from the geometry of Mohr's circle:

tan(2φ) = γ_xy / (ε_x - ε_y)

Solving for φ:

φ = (1/2) arctan[γ_xy / (ε_x - ε_y)]

Important Notes:

• The factor of 2 in tan(2φ) arises because angles on Mohr's circle are doubled relative to physical space
• The arctangent function returns values in the range (-90°, +90°), giving the principal angle range (-45°, +45°)
• At the principal strain directions, shear strain is zero (pure normal strains)
• The maximum shear strain is γ_max = 2R = ε₁ - ε₂, occurring at 45° from the principal directions

Special Cases:

• Isotropic strain state: If ε_x = ε_y and γ_xy = 0, the principal angle is undefined (all directions are principal). The tool sets φ = 0° by convention.
• Pure shear: If ε_x = ε_y but γ_xy ≠ 0, then φ = ±45°
• Zero shear: If γ_xy = 0, then φ = 0° (the x-y axes are already principal axes)

5. Stress Calculations (Optional)

If material elastic properties are known, stresses can be calculated from the measured or calculated strains. This tool supports both isotropic and orthotropic materials under plane stress conditions (free surface where σ_zz = 0).

Orthotropic Stress-Strain Relations:

For an orthotropic material (such as fiber-reinforced composites or rolled metals with directional properties) with principal material axes aligned with the x-y coordinate system:

σ_x = [E_x / (1 - ν_xyν_yx)] (ε_x + ν_yxε_y)

σ_y = [E_y / (1 - ν_xyν_yx)] (ε_y + ν_xyε_x)

τ_xy = G_xy γ_xy

Where:

• E_x = elastic modulus in the x-direction
• E_y = elastic modulus in the y-direction
• G_xy = shear modulus in the xy-plane
• ν_xy = Poisson's ratio (lateral strain in y-direction due to stress in x-direction)
• ν_yx = Poisson's ratio (lateral strain in x-direction due to stress in y-direction)

Reciprocal Relationship: The Poisson's ratios must satisfy ν_xy/E_x = ν_yx/E_y, so ν_yx = ν_xy(E_y/E_x).

Isotropic Materials: If E_x = E_y = E and ν_xy = ν_yx = ν, the material is isotropic. The shear modulus is G_xy = E / [2(1 + ν)]. This tool allows E_y to be left blank for isotropic materials (assumes E_y = E_x).

Microstrain Convention: Strain gages output strain in microstrain (με) where 1 με = 10^-6. The stress equations include a factor of 10^-6 to convert from microstrain to true dimensionless strain:

σ_x = [E_x / (1 - ν_xyν_yx)] (ε_x + ν_yxε_y) × 10^-6
σ_y = [E_y / (1 - ν_xyν_yx)] (ε_y + ν_xyε_x) × 10^-6
τ_xy = G_xy γ_xy × 10^-6

6. Principal Stresses

Once stresses are calculated, the principal stresses (σ₁, σ₂) can be determined using Mohr's circle for stress:

σ_1,2 = (σ_x + σ_y) / 2 ± √[(σ_x - σ_y)² / 4 + τ_xy²]

The principal stress angle:

φ_stress = (1/2) arctan[τ_xy / (σ_x - σ_y)]

Important: For isotropic materials, principal stress directions coincide with principal strain directions (φ_stress = φ_strain). For orthotropic materials, principal stress and strain directions are generally different.

7. Applications

Strain transformation is used in many engineering applications:

• Experimental stress analysis: Converting strain gage measurements to principal strains and stresses
• Structural analysis: Evaluating strain and stress at different orientations to check failure criteria
• Composite materials: Transforming strains between structural coordinates and material principal axes
• Fatigue analysis: Finding maximum principal strains/stresses for fatigue life predictions
• Design verification: Ensuring strains in critical directions remain within allowable limits

8. Summary

The complete strain transformation procedure implemented by this tool:

1. Input: Known strain state (ε_x, ε_y, γ_xy) and rotation angle θ
2. Transform strains: Calculate (ε_x', ε_y', γ_x'y') using transformation equations
3. Calculate principal strains: ε_1,2 = ε_avg ± R using Mohr's circle
4. Determine principal angle: φ = (1/2) arctan[γ_xy / (ε_x - ε_y)]
5. Calculate stresses (optional): If elastic properties provided, compute σ_x, σ_y, τ_xy
6. Calculate principal stresses (optional): σ_1,2 and φ_stress using Mohr's circle for stress

This rigorous framework, based on continuum mechanics and elasticity theory, provides complete characterization of strain and stress states under coordinate transformation.