Stress Recovery Explorer

SEAA 4243 Finite Element Method · QR Explorer

Stress Recovery Explorer

Stresses of the text’s cantilever from its nodal displacements, {σ} = [E][B]{û}. Tap the plate to probe a point; switch the display option the way you will in ANSYS.

Your tasks

  1. 4-node 2 × 2 mesh, probe (3.5, 1.3): the lecture point (Plane Stress Part 3, Activity 1). Compare with your hand calculation.
    Check answer

    σxx = 315, σyy = −45, σxy = −311 kN/m²; σ1 = 494, σ2 = −224; von Mises 637. By hand you need the displacements to five decimals: εyy is the difference of v6 = −0.13042 and v9 = −0.13058, and rounding them to three decimals doubles it.

  2. Keep the point and refine: 8-node 16 × 8, then 32 × 16. Does σxx approach beam theory’s 225 kN/m²? Then probe mid-span, (2.0, 1.3).
    Check answer

    No. It settles near 150 kN/m², because 0.5 m from the loaded edge the uniform end traction is not the parabolic shear that beam theory assumes (St-Venant). At mid-span the FE gives 900 kN/m², equal to beam theory. Check against beam theory only away from the ends.

  3. Run the corner study. Does the peak σxx next to the fixed corner converge?
    Check answer

    No: 4 134 → 5 372 → 6 400 → 7 467 kN/m² as the elements shrink, and it passes the beam value. A re-entrant corner of a fully fixed edge is a stress singularity of the elastic model; do not report that peak.

  4. On the 4-node 4 × 2 mesh compare Elemental mean, Unaveraged and Averaged. Where do neighbouring elements disagree most?
    Check answer

    Near the fixed edge, especially at its corners, and wherever the stress changes fastest. The jumps shrink as you refine: their size is a free indicator of the mesh error. ANSYS shows the same three options.

Element

Mesh

Quantity

Display option

Probe point

Corner study

4-node meshes 8 × 4 to 64 × 32; σxx at the centroid of the element in the top-left corner, against beam theory at the same point.