Rectangular Shaft in Torsion (Saint-Venant)

stressmass-cost

Verified build 12 relations · 2 identities proven · 5 modeling steps · 3 parity samples

Twist a solid circular shaft and nothing surprises you: the shear stress grows straight out from the axis and peaks all the way around the rim (that is the torsion shaft page). Twist a solid rectangular bar and two things happen that no first-year intuition predicts — the peak shear stress sits at the midpoint of the long side, and the corners carry exactly zero stress. This page is about those two facts, the coefficients that quantify them, and why a rectangular shaft is an inefficient way to carry torque.

Convention (fixed for the whole page): aa is the long side, bb the short side, so the aspect ratio a/b1a/b \ge 1. (Sources label these differently — some swap aa and bb, some use 2a×2b2a\times 2b half-sides — so pin it once and read everything below against it.)

τmax=Tc1ab2,θ=Tc2ab3G\tau_{max} = \frac{T}{c_1\,a\,b^2}, \qquad \theta' = \frac{T}{c_2\,a\,b^3\,G}

Here τmax\tau_{max} is the peak shear stress (at the long-side midpoint), θ\theta' is the twist per unit length, and c1(a/b)c_1(a/b), c2(a/b)c_2(a/b) are two dimensionless coefficients that depend only on the shape ratio a/ba/b, not on the absolute size. They come from the classical Saint-Venant solution and are read from a cited table (below).

Why the corners are dead and the long-side middle is hottest

The clean way to see it is Prandtl’s membrane analogy (Timoshenko & Goodier §107): imagine a soap film stretched over a hole the shape of the cross-section and pushed up by a uniform pressure. The height of that film is the torsion stress function, and the shear stress at any point is the slope of the film, running along its contour lines. Two consequences fall straight out:

Drive the aspect ratio in the widget and watch the shear humps in the cross-section: the long-side hump is always the tallest, the corners always dead.

One material property, and the efficiency trap

Look at what each material property touches. GG (the shear modulus) sets the twist θ\theta' and nothing else; σy\sigma_y sets the margin; and τmax\tau_{max} — geometry and torque alone — is blind to the material. Swap steel for aluminium (G7926G \approx 79 \to 26 GPa) and the same bar twists about three times as much per metre at an identical peak stress. That is the same stiffness-versus-strength split the round torsion shaft makes.

Now the comparison that makes the page’s point. Replace the rectangle with a solid round shaft of the same cross-sectional area (same material, same weight per metre): πreq2=ab\pi r_{eq}^2 = a\,b. That round shaft carries the same torque at a lower peak stress and a smaller twist — the rectangle loses on both counts, and loses by more as it gets flatter. A square section carries about 36 % more peak stress than the equal-area circle; a 2 : 1 rectangle, about 62 % more. The corners are dead weight and the flats are underworked: that is why round shafts win, and why a square shaft is a bad deal. The widget shows the equal-area circle as a dashed overlay with its own stress readout.

The coefficients, and where they stop

c1c_1 and c2c_2 are cited data: the classical Saint-Venant coefficients as tabulated by Timoshenko & Goodier (§109). The widget looks them up by a/ba/b and interpolates linearly between the published rows. As the section flattens toward an infinitely thin strip both coefficients climb to the same limit,

c1=c2=13(a/b),c_1 = c_2 = \tfrac{1}{3} \qquad (a/b \to \infty),

which is the exact membrane solution for an infinite strip — the consistency check the build’s physics test pins the whole table against. Two boundaries are enforced, each a global refusal rather than a fabricated number:

Every authored coefficient is cross-checked two independent ways — against the exact Fourier-series solution and against Roark’s separately published closed form — and the results are on the /verification/ page, which states exactly what is machine-proven here (the lookup, the interpolation, the round-shaft algebra, the two refusals) and what rests on citation (the tabulated coefficients themselves). For the thin-walled cousin of this problem — where the section is a closed tube and the shear runs uniformly around the wall — see thin-tube torsion.

Try it

Material

T3, bare flat sheet 0.010-0.128 in. thick, AMS 4037 / AMS-QQ-A-250/4 (MIL-HDBK-5J Table 3.2.3.0(b1), p. 3-71)

Bound properties of 2024-T3 aluminum sheet (bare)
G4 Msitypicalmil-hdbk-5j
sigma_y47 ksidesign min.mil-hdbk-5j
rho0.1 lb/inch**3typicalmil-hdbk-5j
Inputs
N·m
Aspect ratio (long / short)
Stress coefficient c1(a/b)
Stiffness coefficient c2(a/b)
Cross-sectional area
Max shear stress (mid of long side)
Twist rate (per length)
Angle of twist (over L)
Equal-area round radius
Round-shaft shear stress
Round-shaft twist rate
Stress penalty (rect / round)
Twist penalty (rect / round)
Safety factor (shear yield)
Bar mass
kg

4 materials in the database are not listed here: no published value in our cited sources for every property this THING needs.

Materials modeled here: 2024-T3 aluminum sheet (bare) 304 stainless steel 6061-T6 aluminum 7075-T6 aluminum AISI 1045 medium-carbon steel AISI 4340 low-alloy steel (Ni-Cr-Mo) ASTM A36 structural steel (hot-rolled) C26000 Cartridge Brass (70/30) Ti-6Al-4V

Governing relations

a/b    (long÷short)a/b \;\;(\text{long} \div \text{short})

Assumes: a is the long side and b the short side, so the ratio a/b is at least 1; the coefficients c1, c2 depend on this ratio alone, not on the absolute size · Valid while: a is meant to be the LONG side, so a/b must be at least 1. This bar has b longer than a — swap your labels (put the longer dimension in a) and the coefficients apply. The published coefficient table stops at a/b = 10; beyond that there is no tabulated value and interpolation would invent one. For a very thin strip use the closed-form limit c1 = c2 = 1/3 discussed on this page.

Source: Timoshenko, S. P., & Goodier, J. N., Theory of Elasticity, 3rd ed., McGraw-Hill, 1970 — §107 (membrane analogy) and §108-109 (torsion of a bar of rectangular cross section: the coefficient table for maximum stress and torsional rigidity, and the thin-strip limit 1/3).

τmax=Tc1ab2\tau_{max} = \dfrac{T}{c_1\,a\,b^{2}}

Assumes: Saint-Venant pure torsion of a prismatic bar; no warping restraint, no axial or bending load; linear elastic; the peak sits at the midpoint of the long side, the corners carry zero stress · Valid while: Peak shear stress exceeds the shear yield strength (σ_y/2 by the maximum-shear-stress criterion) — the long-side midpoint has yielded and every elastic number here stops being the truth.

Source: Timoshenko, S. P., & Goodier, J. N., Theory of Elasticity, 3rd ed., McGraw-Hill, 1970 — §107 (membrane analogy) and §108-109 (torsion of a bar of rectangular cross section: the coefficient table for maximum stress and torsional rigidity, and the thin-strip limit 1/3).

θ=Tc2ab3G\theta' = \dfrac{T}{c_2\,a\,b^{3}\,G}

Assumes: torsional rigidity of the section is G·c2·a·b^3; c2 is the cited stiffness coefficient; uniform section and torque along the length; twist per unit length is constant

Source: Timoshenko, S. P., & Goodier, J. N., Theory of Elasticity, 3rd ed., McGraw-Hill, 1970 — §107 (membrane analogy) and §108-109 (torsion of a bar of rectangular cross section: the coefficient table for maximum stress and torsional rigidity, and the thin-strip limit 1/3).

θ=θL\theta = \theta' L

Assumes: the twist rate is uniform along the bar, so the total angle is simply theta' times length

Source: Gere, J. M., & Goodno, B. J., Mechanics of Materials, 9th ed., Cengage, 2018 — §3.3-3.5 (circular shafts in torsion: tau = 16T/(πd^3), theta = TL/(GJ), J = πd^4/32), used here for the equal-area round-shaft comparison.

A=abA = a\,b

Assumes: solid rectangular section

Source: Gere, J. M., & Goodno, B. J., Mechanics of Materials, 9th ed., Cengage, 2018 — §3.3-3.5 (circular shafts in torsion: tau = 16T/(πd^3), theta = TL/(GJ), J = πd^4/32), used here for the equal-area round-shaft comparison.

req=abπr_{eq} = \sqrt{\dfrac{a\,b}{\pi}}

Assumes: the round shaft compared against carries the SAME cross-sectional area (same material used)

Source: Gere, J. M., & Goodno, B. J., Mechanics of Materials, 9th ed., Cengage, 2018 — §3.3-3.5 (circular shafts in torsion: tau = 16T/(πd^3), theta = TL/(GJ), J = πd^4/32), used here for the equal-area round-shaft comparison.

τ=2Tπreq3\tau_{\circ} = \dfrac{2\,T}{\pi\,r_{eq}^{3}}

Assumes: solid circular shaft carrying the same torque; tau = 16T/(πd^3) = 2T/(πr^3)

Source: Gere, J. M., & Goodno, B. J., Mechanics of Materials, 9th ed., Cengage, 2018 — §3.3-3.5 (circular shafts in torsion: tau = 16T/(πd^3), theta = TL/(GJ), J = πd^4/32), used here for the equal-area round-shaft comparison.

θ=2Tπreq4G\theta'_{\circ} = \dfrac{2\,T}{\pi\,r_{eq}^{4}\,G}

Assumes: solid circular shaft, polar second moment J = πr^4/2, so theta' = T/(G·J) = 2T/(π r^4 G)

Source: Gere, J. M., & Goodno, B. J., Mechanics of Materials, 9th ed., Cengage, 2018 — §3.3-3.5 (circular shafts in torsion: tau = 16T/(πd^3), theta = TL/(GJ), J = πd^4/32), used here for the equal-area round-shaft comparison.

ητ=τmaxτ\eta_\tau = \dfrac{\tau_{max}}{\tau_{\circ}}

Assumes: ratio of the rectangle's peak stress to the equal-area round shaft's; greater than 1 means the rectangle is worse, and it worsens as the section elongates

Source: Gere, J. M., & Goodno, B. J., Mechanics of Materials, 9th ed., Cengage, 2018 — §3.3-3.5 (circular shafts in torsion: tau = 16T/(πd^3), theta = TL/(GJ), J = πd^4/32), used here for the equal-area round-shaft comparison.

ηθ=θθ\eta_\theta = \dfrac{\theta'}{\theta'_{\circ}}

Assumes: ratio of the rectangle's twist rate to the equal-area round shaft's; also greater than 1

Source: Gere, J. M., & Goodno, B. J., Mechanics of Materials, 9th ed., Cengage, 2018 — §3.3-3.5 (circular shafts in torsion: tau = 16T/(πd^3), theta = TL/(GJ), J = πd^4/32), used here for the equal-area round-shaft comparison.

SF=σy/2τmax\mathrm{SF} = \dfrac{\sigma_y/2}{\tau_{max}}

Assumes: maximum-shear-stress (Tresca) criterion — shear yield at σ_y/2; conservative next to the distortion-energy value 0.577 σ_y

Source: Budynas, R. G., & Nisbett, J. K., Shigley's Mechanical Engineering Design, 10th ed., McGraw-Hill, 2015 — §5-4 (maximum-shear-stress theory, S_sy = 0.5 S_y).

m=ρabLm = \rho\,a\,b\,L

Assumes: prismatic solid bar, uniform density

Source: Gere, J. M., & Goodno, B. J., Mechanics of Materials, 9th ed., Cengage, 2018 — §3.3-3.5 (circular shafts in torsion: tau = 16T/(πd^3), theta = TL/(GJ), J = πd^4/32), used here for the equal-area round-shaft comparison.

Derivation

Steps marked modeling step are where physics enters by citation — every other line is machine-proven to follow from them, and the cited models are independently re-derived in the test pipeline where possible. See what is and isn't machine-verified.

τmax=Tab2c1\tau_{max} = \frac{T}{a b^{2} c_{1}}

1. Prandtl's membrane analogy: the torsion stress function is the height of a soap film blown over the cross-section, and the shear stress is the film's SLOPE, running parallel to the boundary. The film is steepest where the wall is closest to the centre — the midpoint of the LONG side — so the peak stress sits there, not at the corners, where a convex corner forces the slope (and the stress) to zero. Timoshenko tabulates the coefficient c1(a/b); this is cited data, not something the build derives. — modeling: membrane analogy (Timoshenko §107); c1(a/b) from the §109 table modeling step

thetap=TGab3c2thetap = \frac{T}{G a b^{3} c_{2}}

2. The same solution gives the torsional rigidity as G·c2·a·b^3, so the twist per unit length is theta' = T/(c2·a·b^3·G). G is the ONLY material property that touches this page: swap steel for aluminium (G ≈ 26 vs 79 GPa) and the same bar twists about three times as much per metre at an identical peak stress — because tau_max carries no material property at all. — modeling: torsional rigidity G·c2·a·b^3; c2(a/b) from the §109 table modeling step

θ=Lthetap\theta = L thetap

3. The twist rate is uniform along a prismatic bar, so the total angle of twist is just the rate times the length. — kinematics: uniform twist along the length modeling step

req=abπr_{eq} = \frac{\sqrt{a} \sqrt{b}}{\sqrt{\pi}}

4. Now the fair comparison. Replace the rectangle with a solid ROUND shaft of the SAME cross-sectional area (same material, same weight per length): π·r_eq^2 = a·b fixes its radius. — equal-area substitution

τround=2Tπreq3\tau_{round} = \frac{2 T}{\pi r_{eq}^{3}}

5. That round shaft carries the same torque with the familiar tau = 16T/(πd^3) = 2T/(π r^3). The rectangle's tau_max is LARGER (their ratio eta_tau exceeds 1), and the penalty grows as the section flattens — the corners are dead weight that carries almost no shear. This is why a square shaft is a bad deal and a thin strip is a terrible one. — circular torsion formula, equal area

ητ=τmaxτround\eta_{\tau} = \frac{\tau_{max}}{\tau_{round}}

6. Collecting the comparison: eta_tau is the stress penalty (about 1.36 for a square, rising with a/b) and eta_theta the twist penalty. Both exceed 1 for every rectangle — the round section is the efficient one. — definition: efficiency ratios vs the equal-area circle modeling step

SF=σy2τmaxSF = \frac{\sigma_{y}}{2 \tau_{max}}

7. Yield in shear arrives at σ_y/2 by the maximum-shear-stress (Tresca) criterion. Note again what each property touches: G set the twist, σ_y sets this margin, and tau_max — geometry and load alone — is blind to the material. — maximum-shear-stress criterion modeling step

How it fails

The widget answers one narrow, elastic question — given this rectangular section and this torque, what is the peak shear stress and the twist? — and reports a static shear-yield margin on that peak. Most of the ways a bar like this actually gets into trouble live in what the page deliberately does not model.

  • Fixed-Fixed Torsion Shaft (Interior Torque)

    A solid circular shaft built into a wall at BOTH ends, with a torque applied at an interior station. Equilibrium gives one equation for the two wall reaction torques; the missing equation is compatibility — the twist at the load point is single-valued — and the build solves the coupled 2×2 system exactly. The larger reaction lands on the SHORTER segment, and the material cancels out.

    • stress
    • mass-cost
  • Shaft in Torsion (Solid, Circular)

    The power-transmission workhorse: twist a solid circular bar and shear stress winds around it. Three material properties drive three different outputs — stiffness (G) sets the twist, strength (σ_y) sets the margin, and the stress itself doesn't care what the shaft is made of at all.

    • stress
    • torque-power
    • mass-cost
  • Shaft under Combined Bending + Torsion

    Real shafts never get to choose: the belt that twists them also bends them. Bending and torsion land on the same surface element, Mohr's circle finds the worst plane, and two failure criteria — Tresca and von Mises — disagree by up to 15 % about how bad it is.

    • stress
    • mass-cost
  • Thin-Walled Tube in Torsion (Bredt)

    Why driveshafts, airframes, and bike frames are closed tubes: in torsion, what matters is not how much metal you have but how much AREA the wall encloses. Bredt's shear flow makes any closed section solvable with two knobs — and the isoperimetric inequality polices which sections can exist at all.

    • stress
    • mass-cost
  • Cantilever Beam (End Load)

    A beam fixed at one end, loaded at the other — the fruit-fly of structures. One widget shows why stiffness (E) and strength (σ_y) are independent axes: swap steel for titanium and deflection goes UP while the safety factor also goes up.

    • stress
    • mass-cost
  • Composite Bar (Core + Sleeve)

    A solid core inside a concentric sleeve, bonded between rigid end plates and pushed by a centric axial load. The two materials must stretch together, so the load splits in proportion to each member's axial stiffness A·E — and the build solves that coupled 2×2 share exactly. Swap the sleeve's metal and watch the load migrate to the stiffer member.

    • stress
    • mass-cost

Chains with

Outputs whose SI dimension and quantity kind match another THING's input — the only wires the planner's connectionLegal accepts (invariant 2, computed at build time, not hand-listed). Wire these on the chaining demo.

+ 25 more THINGs its outputs can legally feed (showing the first 8 in course order).

Sources