Change specimen geometry, load and fringe constant to visualize isochromatic fringe patterns corresponding to the principal stress difference in real time. Experience the foundation of FEM validation.
Disk (Hertz): $\sigma_y = -\dfrac{2P}{\pi D t},\quad \sigma_x = \dfrac{6P}{\pi D t}$
Denser fringes indicate larger principal stress difference.
What is Photoelasticity Stress Analysis?
🙋
What exactly is photoelasticity? I see colorful fringes in the simulator, but what do they mean?
🎓
Basically, it's a visual stress analysis technique. When you shine polarized light through a transparent, stressed model, the light splits and recombines, creating colorful bands called isochromatic fringes. Each color band corresponds to a specific level of stress difference. In the simulator, try moving the 'Load P' slider—you'll see the fringes get denser, showing where the stress is higher.
🙋
Wait, really? So the colors aren't random? How do we get a number from a color?
🎓
Exactly! Each full cycle from black to black (or a specific color sequence) is one "fringe order," labeled N. The key is that the fringe order is directly proportional to the difference between the two principal stresses ($\sigma_1 - \sigma_2$). For instance, if you change the 'Specimen Shape' from a disk to a beam, the fringe pattern changes completely, revealing the new stress distribution.
🙋
So the thickness and material matter too? What happens if I make the model thicker using the 'Thickness t' control?
🎓
Great question! Yes, thickness is crucial. For the same stress difference, a thicker model creates more "optical retardation," resulting in a higher fringe order N. Try increasing 'Thickness t' in the simulator—you'll see more fringes appear instantly. That's why we need the material's 'Fringe Constant f' to calibrate the system and get actual stress values from the pretty picture.
Physical Model & Key Equations
The fundamental law of photoelasticity relates the observed fringe pattern to the physical stress state in the model. The fringe order N you see is calculated from the principal stress difference, thickness, and a material property.
$$N = \dfrac{(\sigma_1 - \sigma_2)\,t}{f}$$
Where: $N$ = Fringe order (e.g., 1st, 2nd fringe) $\sigma_1, \sigma_2$ = Principal stresses (maximum and minimum normal stress) $t$ = Specimen thickness $f$ = Material fringe constant (a property of the photoelastic material)
For specific geometries, we have analytical stress solutions. For the disk under diametral compression (a classic test), the stresses at the center are given by Hertz's solution.
$$\sigma_y = -\dfrac{2P}{\pi D t},\quad \sigma_x = \dfrac{6P}{\pi D t}$$
Where: $P$ = Applied load $D$ = Disk diameter $t$ = Disk thickness
At the center, $\sigma_1 = \sigma_x$ and $\sigma_2 = \sigma_y$, so the fringe order directly reveals the load. This is what the simulator calculates when you choose the 'Disk' shape.
Frequently Asked Questions
In photoelasticity, the same color repeats each time the order N, which is proportional to the principal stress difference, becomes an integer. Since the order increases as 0 (black) → 1 (red/blue) → 2 (red/blue)..., more fringes appear concentrically in higher stress areas. When counting the order, count the change from black to red as one order.
The material fringe value f is a constant that indicates how sensitively the material exhibits birefringence in response to stress. The smaller f is, the greater the birefringence even for a small stress difference, resulting in a higher order N for the same stress difference. In other words, materials with higher sensitivity produce denser fringes.
Increasing the load P increases the principal stress difference inside the disk, raising the fringe order at the center. Specifically, the number of fringes increases outward from the central black point (order 0), and the concentric fringes become denser. The load and fringe order are proportional.
This simulator is based on theoretical solutions such as Hertzian contact theory, and under ideal conditions, it shows high agreement with actual photoelastic experiments. However, real test specimens involve material inhomogeneity, residual stress, and edge effects, so it is suitable for qualitative trend understanding and basic learning for FEM verification.
Real-World Applications
Structural Component Design: Engineers use photoelasticity to find stress concentrations in complex parts like gear teeth, crane hooks, or load-bearing brackets. By coating the actual part with a photoelastic material or using a scaled model, they visually identify where fringes are densest—the points most likely to fail—and then reinforce those areas.
Biomechanics and Implant Analysis: Researchers analyze stress transfer in bone-implant systems, such as dental implants or hip replacements. A photoelastic model of bone with an embedded implant shows how stress flows around it, helping to design shapes that minimize bone stress-shielding and promote integration.
Validation of Computer Simulations (FEA): Before trusting complex Finite Element Analysis (FEA) software, engineers often validate their digital models against photoelastic experiments. The physical fringe pattern provides a direct, trustworthy benchmark to check if the computer's stress predictions are accurate.
Geotechnical and Rock Mechanics: Scientists model how stress propagates through granular materials or rock layers under load. Photoelastic models using crushed glass or specialized polymers reveal force chains and contact stresses, which is crucial for understanding foundation stability and tunnel safety.
Common Misconceptions and Points to Note
First, you might tend to think "darker fringes = higher stress," but strictly speaking, this is not accurate. Isochromatic fringes represent the principal stress "difference." For example, uniform tension (large σ1, σ2=0) and uniform shear (σ1 and σ2 have the same absolute value but opposite signs) result in completely different stress states even with the same stress difference. Try selecting the "disk" in the simulator and applying a load to observe its center. Fringes appear due to a large stress difference from a combination of tension and compression. On the other hand, in a simple tensile test specimen, applying a load initially shows almost no fringes (because σ2 is near zero, resulting in a small difference). Always remember you are observing the stress "difference," not its "magnitude."
Next, handling the fringe constant "f". This value is material-specific. For instance, it is about 3~5 kN/m·fringe for epoxy resin and about 10~15 kN/m·fringe for acrylic, varying significantly between materials. If you change this value in the simulator, you'll see the number of fringes change drastically even under the same load, right? Knowledge gained from experiments with a known material (e.g., epoxy) in practice cannot be directly applied to designing another material (e.g., aluminum). This is because the difference in the f-value alters the "apparent risk level."
Finally, the fundamental limitation that only "in-plane stresses in 2D" can be evaluated. This method assumes that the stress in the thickness direction, through which light passes, is constant (a state of plane stress). In actual thick components, "3D stress" occurs where the stress state differs between the surface and the interior. Therefore, this method is not suitable for precise evaluation of complex 3D shapes; it's crucial to understand its role as primarily for trend identification or FEM verification.
Select material type (polycarbonate, epoxy, or acrylic) using the material dropdown—each has different stress-optic coefficients (fσ range 0.08–0.12 MPa⁻¹).
Set geometry dimensions (vDNum for diameter/width in mm, vTNum for thickness in mm) and apply load magnitude (sP in kN) via the load input field.
Adjust specimen parameters (vPNum for load position, sPNum for support spacing) and observe real-time fringe order N calculation and principal stress values (σ₁, σ₂, τmax in MPa) updating in the output panel.
Rotate or zoom the 3D fringe pattern visualization to identify isochromatic bands where color transitions indicate equal stress magnitude differences.
Worked Example
Circular polycarbonate disc (diameter 60 mm, thickness 8 mm, fσ = 0.10 MPa⁻¹) under diametral compression load of 5 kN. Using standard disc compression geometry with 50 mm support spacing: calculated principal stress difference (σ₁ − σ₂) = 48 MPa yields fringe order N = 4.8, displayed as alternating red-blue-green isochromatic bands across the disc diameter. Maximum shear stress τmax = 24 MPa peaks at the load contact point.
Practical Notes
High fringe orders (N > 10) indicate stress concentration zones; use finer material thickness or reduce load magnitude for clearer fringe separation in tensile specimen necks or notched geometries.
Isochromatic patterns are independent of material birefringence direction—use isoclinic lines (dark fringe family, separate analysis) to map principal stress orientations simultaneously.
For edge-loaded rectangular beams: maximum fringe order appears near neutral axis centerline; compare experimental wraparound patterns to FEA predictions when validating stress concentration factors.