Strain Gage to Strain State

×

Gage-to-State

Strain Gage Rosette to Strain State Theory

Introduction

Strain gages are electrical resistance devices that measure normal strain in a single direction. When the directions of principal strains are unknown (as is common in complex structures), a single gage is insufficient to characterize the complete state of strain at a point. A strain gage rosette consists of three gages arranged in a tight pattern at known angles, allowing determination of the complete plane strain state and subsequent calculation of principal strains and stresses.

This derivation shows the mathematical process used by this tool to convert three measured strains (ε_A, ε_B, ε_C) at known angles (θ_A, θ_B, θ_C) into:

1. Engineering strain components (ε_xx, ε_yy, γ_xy)
2. Principal strains (ε₁, ε₂) and their orientation (φ)
3. Stresses (σ_xx, σ_yy, τ_xy) if material properties are known

Figure 1: Three-element strain gage rosette showing gage orientation angles θ_A, θ_B, θ_C measured counterclockwise from the x-axis

1. Plane Strain State

For a point on the surface of a loaded structure (plane stress condition, with σ_zz = τ_xz = τ_yz = 0), the state of strain is completely defined by three independent quantities:

• ε_xx = normal strain in the x-direction (dimensionless)
• ε_yy = normal strain in the y-direction (dimensionless)
• γ_xy = engineering shear strain in the xy-plane (dimensionless)

Note: The engineering shear strain γ_xy is related to the tensor shear strain ε_xy by γ_xy = 2ε_xy. Engineering shear strain represents the total angular distortion, while tensor shear strain is half this value. This tool uses engineering shear strain throughout, consistent with standard mechanics of materials practice.

2. Strain Transformation Equation

Consider a strain gage oriented at an angle θ (measured counterclockwise from the positive x-axis). The normal strain ε_θ measured by this gage is related to the strain state components by the strain transformation equation:

ε_θ = ε_xx cos²θ + ε_yy sin²θ + γ_xy sinθ cosθ

Derivation: This equation arises from coordinate transformation of the strain tensor. Starting with the strain tensor in the xy-coordinate system:

[ε] =

εxx	εxy
εxy	εyy

  where ε_xy = γ_xy / 2

Transforming to a rotated coordinate system (x', y') at angle θ using the rotation matrix R:

R =

cosθ	sinθ
−sinθ	cosθ

The transformed strain tensor is [ε'] = R [ε] R^T, and the normal strain in the x'-direction (the gage direction) is the (1,1) component:

ε_θ = ε'₁₁ = ε_xx cos²θ + ε_yy sin²θ + 2ε_xy sinθ cosθ

Substituting ε_xy = γ_xy / 2:

ε_θ = ε_xx cos²θ + ε_yy sin²θ + γ_xy sinθ cosθ

This is the fundamental equation that relates the measured strain at any angle θ to the unknown strain state (ε_xx, ε_yy, γ_xy).

3. Rosette Transformation Matrix Formulation

A three-element rosette provides three measured strains (ε_A, ε_B, ε_C) at three known angles (θ_A, θ_B, θ_C). Applying the transformation equation to each gage:

ε_A = ε_xx cos²θ_A + ε_yy sin²θ_A + γ_xy sinθ_A cosθ_A

ε_B = ε_xx cos²θ_B + ε_yy sin²θ_B + γ_xy sinθ_B cosθ_B

ε_C = ε_xx cos²θ_C + ε_yy sin²θ_C + γ_xy sinθ_C cosθ_C

This is a system of three linear equations with three unknowns, which can be written in matrix form:

εA

εB

εC

 = 

cos2θA	sin2θA	sinθA cosθA
cos2θB	sin2θB	sinθB cosθB
cos2θC	sin2θC	sinθC cosθC

 

εxx

εyy

γxy

Or more compactly: {ε_measured} = [T] {ε_state}, where [T] is the transformation matrix.

Important: The third column uses sinθcosθ (NOT 2sinθcosθ). This directly corresponds to the engineering shear strain γ_xy in the strain state vector. Some references show 2sinθcosθ, but that formulation is for the tensor shear strain ε_xy = γ_xy/2. Since this tool uses engineering shear strain throughout, the factor of 2 is omitted from the transformation matrix.

This is implemented in output.php (lines 100, 104, 108), where T_i,3 = sinθ_i cosθ_i.

4. Matrix Inversion Method

To solve for the unknown strain state, we invert the transformation matrix:

{ε_state} = [T]⁻¹ {ε_measured}

For a general 3×3 matrix, the inverse is calculated using the cofactor method:

[T]⁻¹ = (1 / det[T]) [C]^T

where det[T] is the determinant and [C]^T is the transpose of the cofactor matrix.

Step 1: Calculate Determinant

For the 3×3 transformation matrix [T] with elements T_ij:

det[T] = T₁₁(T₂₂T₃₃ − T₂₃T₃₂) − T₁₂(T₂₁T₃₃ − T₂₃T₃₁) + T₁₃(T₂₁T₃₂ − T₂₂T₃₁)

This is implemented in output.php line 116. A non-zero determinant indicates the three gage angles are linearly independent (not all colinear or degenerate), which is required for a valid rosette configuration.

Step 2: Calculate Cofactor Matrix

The cofactor C_ij is (−1)^i+j times the determinant of the 2×2 minor matrix obtained by deleting row i and column j:

C₁₁ = T₂₂T₃₃ − T₂₃T₃₂
C₁₂ = −(T₂₁T₃₃ − T₂₃T₃₁)
C₁₃ = T₂₁T₃₂ − T₂₂T₃₁

C₂₁ = −(T₁₂T₃₃ − T₁₃T₃₂)
C₂₂ = T₁₁T₃₃ − T₁₃T₃₁
C₂₃ = −(T₁₁T₃₂ − T₁₂T₃₁)

C₃₁ = T₁₂T₂₃ − T₁₃T₂₂
C₃₂ = −(T₁₁T₂₃ − T₁₃T₂₁)
C₃₃ = T₁₁T₂₂ − T₁₂T₂₁

These are calculated in output.php lines 119-129.

Step 3: Compute Inverse Matrix Elements

The inverse matrix elements are the transposed cofactors divided by the determinant:

T⁻¹_ij = C_ji / det[T]

(Note the index swap in C_ji due to the transpose operation.) This is implemented in output.php lines 132-142.

Step 4: Solve for Strain State

Matrix multiplication yields the strain components:

ε_xx = T⁻¹₁₁ε_A + T⁻¹₁₂ε_B + T⁻¹₁₃ε_C
ε_yy = T⁻¹₂₁ε_A + T⁻¹₂₂ε_B + T⁻¹₂₃ε_C
γ_xy = T⁻¹₃₁ε_A + T⁻¹₃₂ε_B + T⁻¹₃₃ε_C

This is the complete solution for the engineering strain state from rosette measurements (output.php lines 151-153).

5. Principal Strains via Mohr's Circle

Once the strain state (ε_xx, ε_yy, γ_xy) is known, the principal strains (maximum and minimum normal strains) and their orientation can be determined using Mohr's circle for strain.

Mohr's Circle Construction:

Mohr's circle is a graphical and analytical method for finding principal values and orientations. For strain, the circle is plotted in (ε, γ/2) space, with:

• Center: (ε_avg, 0) where ε_avg = (ε_xx + ε_yy) / 2
• Radius: R = √[(ε_xx − ε_yy)² / 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)

Substituting the expressions for ε_avg and R:

ε_1,2 = (ε_xx + ε_yy) / 2 ± √[(ε_xx − ε_yy)² / 4 + (γ_xy / 2)²]

This is implemented in output.php lines 161-166.

Principal Angle:

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

tan(2φ) = γ_xy / (ε_xx − ε_yy)

Solving for φ:

φ = (1/2) arctan[γ_xy / (ε_xx − ε_yy)]

This gives the counterclockwise angle from the x-axis to the maximum principal strain direction (output.php line 169).

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°), which corresponds to the principal angle range (−45°, +45°)
• At ε₁ and ε₂ directions, shear strain is zero (pure normal strains)
• The maximum shear strain is γ_max = 2R, occurring at 45° from the principal directions

6. Orthotropic Stress-Strain Relations (Plane Stress)

If the material's elastic properties are known, stresses can be calculated from the measured strains. For orthotropic materials (such as fiber-reinforced composites or rolled metals with directional properties), the stress-strain relationships are more complex than for isotropic materials.

Plane Stress Assumption: On the free surface of a structure, σ_zz = τ_xz = τ_yz = 0. The in-plane stresses (σ_xx, σ_yy, τ_xy) are related to in-plane strains by the plane stress constitutive equations.

Orthotropic Constitutive Equations:

For an orthotropic material with principal material axes aligned with the x-y coordinate system, the stress-strain relationships are:

σ_xx = (E_x / (1 − ν_xyν_yx)) (ε_xx + ν_yxε_yy)
σ_yy = (E_y / (1 − ν_xyν_yx)) (ε_yy + ν_xyε_xx)
τ_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 for Poisson's Ratios:

The Poisson's ratios are not independent. From strain energy considerations (symmetry of the compliance matrix), they must satisfy:

ν_xy / E_x = ν_yx / E_y

This means if E_x, E_y, and ν_xy are known, then ν_yx = ν_xy (E_y / E_x).

Isotropic Special Case: If E_x = E_y = E and ν_xy = ν_yx = ν, the material is isotropic, and the shear modulus is related by G_xy = E / [2(1 + ν)]. This tool allows E_y to be left blank, in which case it assumes E_y = E_x (isotropic assumption), as seen in output.php lines 192-198.

Microstrain Conversion:

Strain gages typically output strain in microstrain (με), where 1 με = 10⁻⁶. True (dimensionless) strain is obtained by dividing by 10⁶. The stress equations become:

σ_xx = (E_x / (1 − ν_xyν_yx)) (ε_xx + ν_yxε_yy) × 10⁻⁶
σ_yy = (E_y / (1 − ν_xyν_yx)) (ε_yy + ν_xyε_xx) × 10⁻⁶
τ_xy = G_xy γ_xy × 10⁻⁶

This is exactly as implemented in output.php lines 206-208.

Units Convention: If moduli are input in ksi (1000 psi) or MPa (10⁶ Pa), and strains in με, the resulting stresses will be in ksi or MPa, respectively. This tool automatically handles unit conversions and displays stresses in engineering units appropriate to the selected unit system (KSI for USCS, MPa for SI).

7. Principal Stresses via Mohr's Circle

Once the stress state (σ_xx, σ_yy, τ_xy) is calculated, the principal stresses (maximum and minimum normal stresses) and their orientation can be determined using Mohr's circle for stress.

Important Distinction - Principal Stress vs Principal Strain Directions:

• For isotropic materials: Principal stress directions coincide with principal strain directions (φ_stress = φ_strain)
• For orthotropic materials: Principal stress and strain directions are generally DIFFERENT (φ_stress ≠ φ_strain)

Mohr's Circle for Stress:

Similar to the strain analysis, Mohr's circle for stress is plotted in (σ, τ) space, with:

• Center: (σ_avg, 0) where σ_avg = (σ_xx + σ_yy) / 2
• Radius: R_stress = √[(σ_xx − σ_yy)² / 4 + τ_xy²]

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

σ₁ = σ_avg + R_stress (maximum principal stress)
σ₂ = σ_avg − R_stress (minimum principal stress)

Substituting the expressions:

σ_1,2 = (σ_xx + σ_yy) / 2 ± √[(σ_xx − σ_yy)² / 4 + τ_xy²]

This is implemented in output.php lines 323-328.

Principal Stress Angle:

The angle φ_stress from the x-axis to the σ₁ direction is:

φ_stress = (1/2) arctan[τ_xy / (σ_xx − σ_yy)]

This gives the counterclockwise angle from the x-axis to the maximum principal stress direction (output.php lines 339-350, with divide-by-zero handling for special cases).

Special Cases:

• Pure shear stress: If σ_xx = σ_yy but τ_xy ≠ 0, then φ_stress = ±45°
• Isotropic stress: If σ_xx = σ_yy and τ_xy = 0, principal angle is undefined (set to 0° by convention)
• Zero shear: If τ_xy = 0, then φ_stress = 0° (principal axes align with x-y axes)

8. Summary of Complete Analysis Procedure

The complete strain gage rosette analysis procedure implemented by this tool is:

1. Measure: Three strains (ε_A, ε_B, ε_C) at known angles (θ_A, θ_B, θ_C)
2. Construct transformation matrix [T]: Using cos²θ, sin²θ, sinθcosθ for each gage
3. Invert matrix: Calculate det[T], cofactor matrix [C], and [T]⁻¹ = (1/det)[C]^T
4. Solve for strain state: {ε_xx, ε_yy, γ_xy} = [T]⁻¹ {ε_A, ε_B, ε_C}
5. Calculate principal strains: ε_1,2 = ε_avg ± R via Mohr's circle
6. Determine principal strain angle: φ = (1/2) arctan[γ_xy / (ε_xx − ε_yy)]
7. Calculate stresses (optional): If E_x, E_y, G_xy, ν_xy are known, compute σ_xx, σ_yy, τ_xy using orthotropic constitutive equations
8. Calculate principal stresses (optional): σ_1,2 = σ_avg ± R_stress and φ_stress = (1/2) arctan[τ_xy / (σ_xx − σ_yy)]

This rigorous mathematical framework, based on fundamental continuum mechanics and elasticity theory, allows complete characterization of the stress and strain state at a point from three simple strain gage measurements.

9. Common Rosette Configurations

Rectangular (45°) Rosette: θ_A = 0°, θ_B = 45°, θ_C = 90°
This is the most common configuration, providing excellent numerical conditioning and simplicity. The transformation matrix for this configuration has a particularly simple form.

Delta (60°) Rosette: θ_A = 0°, θ_B = 60°, θ_C = 120°
Provides 120° angular separation between gages, often used when space constraints favor a more compact rosette pattern.

Tee (Two-Element) Rosette: θ_A = 0°, θ_B = 90°
Uses only two gages when the principal directions are known to align with the gage axes. Cannot determine shear strain independently; assumes γ_xy = 0 at the known principal orientation.

This tool supports arbitrary rosette angles, allowing analysis of custom gage configurations or cases where the rosette is mounted at an angle to the structural coordinate system.

10. Practical Considerations

Measurement Errors: Strain gage measurements are subject to errors from temperature effects, gage misalignment, adhesive creep, and electrical noise. Proper installation techniques, temperature compensation (using dummy gages or self-temperature-compensated gages), and careful calibration are essential for accurate results.

Singularity Conditions: If the three gage angles are colinear or degenerate (e.g., all at the same angle), the transformation matrix becomes singular (det[T] = 0), and the strain state cannot be determined. This tool will report an error or infinite values in such cases.

Spatial Resolution: The three gages in a rosette occupy a finite area (typically a few millimeters). In regions with steep strain gradients (near crack tips, notch roots, or concentrated loads), the gages may average over a non-uniform strain field. Smaller rosettes or finite element analysis may be required for such cases.

Material Orientation: For orthotropic materials, the stress-strain relationships provided assume the x-y axes align with the material's principal directions. If the rosette is mounted at an arbitrary angle to the material axes, additional coordinate transformation is required. Advanced analysis may involve transforming strains to material coordinates, applying constitutive equations, then transforming stresses back to the measurement coordinate system.

Close