Mohr's Circle and Principal Stress Calculator

Give this calculator a plane-stress state — the two normal stresses σx and σy and the shear stress τxy acting on a small element — and it returns the principal stresses σ1 and σ2, the maximum in-plane shear stress, the angle you would have to rotate the element to reach the principal planes, and the von Mises equivalent stress used in yield checks. It also transforms the state to any angle you choose and plots the Mohr's circle itself, so you can see why the answers are what they are rather than just reading them off.

Calculator

This calculator runs in your browser. Enable JavaScript for live results — the inputs, formula and worked example below remain fully readable without it.

Inputs this calculator takes, with typical values
InputWhat to enterExample
Normal stress σxPositive in tension, negative in compression, acting on the face whose normal points along x.80 MPa
Normal stress σyPositive in tension, negative in compression, acting on the face whose normal points along y.20 MPa
Shear stress τxyPositive when it acts on the +x face in the +y direction — the standard tensor sign convention.40 MPa
Rotation angle θ for transformed stressesCounter-clockwise rotation of the element from the x axis; set it to 0 to leave the state untransformed.30 °

It returns

  • Maximum principal stress σ₁ — The largest in-plane normal stress. No shear acts on the plane it lies on.
  • Minimum principal stress σ₂
  • Maximum in-plane shear τmax (circle radius R)
  • Average stress (circle centre)
  • Principal plane angle θp — Rotate the element by this angle counter-clockwise and σ₁ acts on the new x face.
  • Von Mises equivalent stress
  • Transformed normal stress σx′ at θ
  • Transformed shear stress τx′y′ at θ

The formula

σ1,2=σx+σy2±(σxσy2)2+τxy2
σx=σx+σy2+σxσy2cos2θ+τxysin2θ
σvM=σ12σ1σ2+σ22

In plain text: σ₁,₂ = (σx + σy)/2 ± √( ((σx − σy)/2)² + τxy² )

  • σ₁, σ₂Maximum and minimum in-plane principal stresses (MPa)
  • σx, σyNormal stresses on the x and y faces of the element (MPa)
  • τxyShear stress on the element faces (MPa)
  • RRadius of Mohr's circle, equal to the maximum in-plane shear stress (MPa)
  • θpCounter-clockwise rotation from x to the plane carrying σ₁ (°)

Plane stress assumes the third principal stress is zero — a free surface, a thin plate loaded in its plane, or the outer fibre of a beam or shaft.

Updated Category Engineering FE & PE Exam Math Verified against published test cases Reading time 12 min

What principal stresses are and why you need them

The stresses you calculate from a load — bending stress from My/I, torsional shear from Tr/J, direct stress from P/A — are tied to the axes you happened to choose. Rotate the element and every one of those numbers changes. The principal stresses are the values you find on the one orientation where the shear stress vanishes entirely, leaving pure tension and pure compression. They are properties of the stress state itself, not of your coordinate choice.

That matters for three practical reasons. Brittle materials crack on the plane of maximum tensile stress, so σ1 tells you where the crack will start and which way it will run. Ductile materials yield on planes of maximum shear, so τmax and the von Mises stress predict yielding. And fatigue life is driven by the principal stress range, not by the component stresses.

Christian Otto Mohr's 1882 construction turns the algebra into a circle. Plot each face of the element as a point (σ, τ) and every possible orientation of that element lies on one circle whose centre is the average normal stress and whose radius is the maximum shear stress. Rotating the physical element by θ moves you 2θ around the circle. Once you see that, most transformation problems become geometry rather than trigonometry.

Reading the formula off the circle

Start with the two quantities that define the circle. The centre sits at (σxy)/2 on the σ axis, always — the sum of the normal stresses is invariant, so no rotation moves it. The radius is R = √( ((σx−σy)/2)² + τxy² ), which is just Pythagoras applied to the point representing the x face: it sits (σx−σy)/2 to the right of the centre and τxy above it.

Everything else follows. The circle crosses the σ axis at centre ± R, and those two crossings are σ1 and σ2. They are the extreme normal stresses, and because they lie on the horizontal axis, the shear there is zero — which is the definition of a principal plane. The top and bottom of the circle are at ±R, so the maximum in-plane shear stress equals the radius. Those planes are a quarter turn around the circle from the principal planes, which is 45° in the real element.

The angle is the part students most often get wrong. tan2θp = 2τxy / (σx − σy), and a plain arctangent loses the quadrant. This calculator uses atan2 with the numerator and denominator kept separate, which resolves it correctly and always returns the angle to the σ1 plane rather than to whichever principal plane happens to be nearer. Halving the result converts circle angle to element angle.

The von Mises stress is not a Mohr's circle quantity but is calculated from the principal stresses: σvM = √(σ1² − σ1σ2 + σ2²) for plane stress. Compare it directly to the yield strength — if it exceeds it, the material yields somewhere in the element.

Worked example: σx = −20, σy = 90, τxy = 60 MPa

This is a classic mechanics-of-materials state: compression in one direction, tension in the other, with substantial shear.

  1. Centre. (−20 + 90) ÷ 2 = 70 ÷ 2 = 35 MPa.
  2. Half-difference. (−20 − 90) ÷ 2 = −110 ÷ 2 = −55 MPa. This is the horizontal offset of the x-face point from the centre; it is negative, so that point sits to the left.
  3. Radius. R = √((−55)² + 60²) = √(3,025 + 3,600) = √6,625 = 81.394 MPa.
  4. Principal stresses. σ1 = 35 + 81.394 = 116.394 MPa. σ2 = 35 − 81.394 = −46.394 MPa. Check the invariant: 116.394 + (−46.394) = 70 = −20 + 90. ✓
  5. Maximum in-plane shear. τmax = R = 81.394 MPa, acting on planes where the normal stress is the average, 35 MPa.
  6. Principal angle.p = atan2(2 × 60, −20 − 90) = atan2(120, −110). The reference angle is arctan(120/110) = 47.490°, and because the point is in the second quadrant, 2θp = 180 − 47.490 = 132.510°. So θp = 66.255° counter-clockwise.
  7. Maximum-shear planes. 66.255 − 45 = 21.255°.
  8. Von Mises. √(116.394² − 116.394 × (−46.394) + (−46.394)²) = √(13,547.6 + 5,399.9 + 2,152.4) = √21,099.9 = 145.26 MPa. Against a 250 MPa mild-steel yield, the element has a factor of 1.72 against yielding.

Notice how much larger σ1 is than either input normal stress. Shear stress always pushes the principal values further apart than the normal stresses you started with, which is why a shaft in combined bending and torsion is more highly stressed than either loading alone suggests.

How to interpret the result

Compare the right quantity to the right property. For a ductile metal, compare the von Mises stress to the yield strength. If you prefer the more conservative Tresca criterion, compare twice the absolute maximum shear stress to the yield strength instead — Tresca and von Mises agree in uniaxial tension and differ by at most 15.5%, with Tresca always the safer of the two. For a brittle material such as grey cast iron or concrete, compare σ1 to the tensile strength, because it fails on the maximum-tension plane.

Do not forget the third principal stress. Plane stress means σ3 = 0, and zero is still a principal stress. When both in-plane principals are tensile, σ3 = 0 is the smallest of the three, so the true maximum shear stress is σ1/2 on an out-of-plane plane, which is larger than the in-plane radius R. The same applies with the signs reversed when both are compressive. This calculator reports both the in-plane value and the absolute value in the steps, and warns you when they differ.

The sign of θp follows the sign convention you used for τxy. A positive τxy acts on the +x face in the +y direction. Many textbooks plot Mohr's circle with the shear axis inverted so that rotations on the circle and in the element go the same way; the numerical answers are identical either way, only the picture flips.

A circle of zero radius is a real answer, not an error. Equal normal stresses with no shear is the hydrostatic state, where every plane is principal and no plane carries shear. That is why a sphere under uniform pressure has no preferred failure plane.

Principal stresses for common loading states

Each row is derived from the same formula with σ representing the applied magnitude. Use it to sanity-check a numerical answer before trusting it.
Stress stateσ₁σ₂τmax (in-plane)θpvon Mises
Uniaxial tension σσ00.5σσ
Uniaxial compression −σ0−σ0.5σσ
Pure shear ττ−ττ45°1.7321τ
Equal biaxial tension σ, σσσ0anyσ
Equal tension and compression σ, −σσ−σσ1.7321σ
Thin cylinder under pressure (2σ, σ)σ0.5σ1.7321σ

The pure-shear row is the origin of the shear yield strength τy = σy/√3 = 0.577σy used in the von Mises criterion; Tresca gives 0.5σy for the same quantity.

Pitfalls that produce a wrong principal stress

  • Using arctan instead of atan2 for the angle. arctan(2τ/(σx−σy)) returns a value between −90° and +90° for 2θp, which silently maps some states onto the σ₂ plane instead of the σ₁ plane. The error is exactly 90° in the element and easy to miss.
  • Halving the wrong thing. The circle angle is 2θ. Compute 2θp first, then divide by two. Halving the input stresses is a different and much more visible mistake; halving at the wrong step is not.
  • Forgetting that σ₃ = 0 is a principal stress. For a biaxial tensile state, the in-plane radius understates the true maximum shear stress, so a Tresca check based on R alone is unconservative.
  • Mixing sign conventions for shear. Some texts define positive shear as acting counter-clockwise on the element. The principal stresses are unaffected because τxy is squared, but the sign of θp flips.
  • Applying plane stress where the body is thick. Deep inside a thick component under multiaxial load, σ₃ is not zero and you need the full three-dimensional principal stress solution, which is a cubic equation rather than a circle.
  • Comparing von Mises stress to ultimate tensile strength. Von Mises predicts the onset of yielding, so it belongs against the yield strength. Fracture in ductile materials involves large plastic strain that the criterion does not model.

Where this fits in a strength check

Mohr's circle sits between the load analysis and the failure check. You compute component stresses first — P/A for axial, Mc/I for bending, Tr/J for torsion, VQ/It for transverse shear — superpose them at the point of interest, then transform. In a shaft under combined bending and torsion, the critical point is the outer surface where bending stress is maximum and transverse shear is zero, giving σx = Mc/I, σy = 0, τxy = Tr/J. Feed those three numbers in and the von Mises output is the number you compare to the allowable.

The same transformation mathematics applies to strain, with strain gauge rosettes as the everyday application: three measured normal strains give you the full plane-strain state, and the identical circle construction returns the principal strains and their direction. It also applies to the second moment of area, which is why a Mohr's circle of inertia gives the principal axes of an angle section — the axes about which the buckling check and the section modulus must be evaluated.

For beams in ordinary bending, the principal stress at the extreme fibre is simply the bending stress and no transformation is needed; the transformation earns its keep at the neutral axis of a thin web, and in any member where torsion is present. If your problem is a deflection limit rather than a strength limit, the deflection check is the governing one instead, and stress transformation does not enter it.

Key terms

Plane stress
A state in which one principal stress is zero. It holds at any free surface, in thin plates loaded in their own plane, and at the outer fibre of beams and shafts — which covers most points where failure actually starts.
Principal plane
An orientation on which the shear stress is zero. There are two in plane stress, always 90° apart in the element and 180° apart on the circle.
Stress invariant
A combination of stress components that does not change with rotation. In plane stress, σx + σy and σxσy − τxy² are both invariant, which gives you two independent checks on any hand transformation.
Von Mises stress
A single equivalent uniaxial stress that produces the same distortion energy as the actual multiaxial state. Comparing it to the yield strength is the standard ductile yield criterion.
Tresca criterion
Yield occurs when the maximum shear stress reaches half the yield strength. Simpler than von Mises and never less conservative, with the largest gap — 15.5% — occurring in pure shear.

Frequently asked questions

Which angle does the calculator return, θp to σ₁ or to σ₂?

It returns the counter-clockwise angle from the x axis to the plane carrying σ₁, the larger principal stress. Because it uses atan2 with the numerator 2τxy and denominator σx − σy kept separate, the quadrant is resolved correctly and you never have to add 90° yourself. The plane carrying σ₂ is at θp + 90°, and the two maximum-shear planes are at θp ± 45°.

Why is my maximum shear stress larger than the circle radius?

Because the third principal stress matters. Plane stress means σ₃ = 0, and when both in-plane principals have the same sign, zero lies outside the pair. The absolute maximum shear stress is then half the spread across all three values — σ₁/2 if both are tensile, |σ₂|/2 if both are compressive — and that shear acts on a plane inclined out of the xy plane. The in-plane radius R is still the maximum shear you would find by rotating within the plane.

How do I get the stress state for a shaft in bending and torsion?

Take the outer surface at the point of maximum bending moment. There σx = Mc/I (or 32M/πd³ for a solid round shaft), σy = 0, and τxy = Tr/J (16T/πd³). Enter those three and read the von Mises output. For a solid circular shaft this reduces to the familiar equivalent-stress expression √(σ² + 3τ²), which you can use to check the calculator by hand.

What sign convention should I use for τxy?

Positive τxy acts on the face whose outward normal points in +x, in the +y direction — the standard tensor convention. The principal stresses and the maximum shear stress do not depend on the choice because τxy is squared, but the sign of θp does: reversing the convention reverses the direction of rotation. If you are following a textbook that plots the shear axis downward, expect its θp to be the negative of the one shown here.

Is von Mises or Tresca the right criterion to use?

Von Mises matches experimental yield data for ductile metals more closely and is what finite-element software reports by default. Tresca is simpler, always equal to or more conservative, and is written into several pressure-vessel codes for that reason. The two agree exactly in uniaxial tension and differ most in pure shear, where Tresca predicts yielding at 0.500σy and von Mises at 0.577σy — a 15.5% gap. If a code governs your work, use the criterion that code names.

Can I use this for strain instead of stress?

Yes, with one substitution: enter εx, εy and half the engineering shear strain, γxy/2, in place of σx, σy and τxy. The transformation mathematics is identical because the strain tensor has the same structure. The outputs then read as principal strains and half the maximum engineering shear strain. The von Mises output will not be meaningful in that mode, since it is defined on stresses.

Why is Mohr's circle taught at all when the formula is short?

Because it makes several results obvious that the algebra hides: that the principal planes are always 90° apart, that the maximum-shear planes are always 45° from them, that the average normal stress never changes with rotation, and that a state with equal normal stresses and no shear has no preferred direction at all. Those relationships come up constantly in failure analysis, and reasoning about them geometrically is faster and less error-prone than re-deriving them.

What does a zero-radius circle mean physically?

It means the stress state is hydrostatic in the plane: the same normal stress in every direction and no shear on any plane. A thin plate loaded equally in both in-plane directions is the usual example. Because ductile yielding is driven by shear, a purely hydrostatic state produces no yielding no matter how large the stress — which is why the von Mises expression subtracts the σ₁σ₂ cross term and returns σ, not 2σ, for that case.

References

  • Mechanics of Materials, 10th ed. — plane stress transformation and Mohr's circle — R. C. Hibbeler, Pearson
  • Advanced Mechanics of Materials, 6th ed. — A. P. Boresi and R. J. Schmidt, Wiley
  • Roark's Formulas for Stress and Strain, 9th ed. — W. C. Young, R. G. Budynas and A. M. Sadegh, McGraw-Hill
  • FE Reference Handbook — Mechanics of Materials section — National Council of Examiners for Engineering and Surveying (NCEES)