◀ Back to Note

Tony's Intertwingled Memex • Reference Manual

Circular Drive Shaft Torsional Shear Stress (Mechanical Engineering) Solution

Formulas & Definitions

A mechanical design engineer at the UNLV Howard R. Hughes College of Engineering is analyzing the torsional load capacity of a solid alloy steel drive shaft powering an industrial centrifugal water chiller.

The drive shaft has a circular cross-section of radius R=3 cmR = 3\text{ cm}R=3 cm. When transmitted engine torque twists the shaft, internal elastic torsional shear stress τ\tauτ develops. According to the mechanics of materials elastic torsion formula, shear stress varies directly with radial distance from the central rotational axis, modeled in Cartesian coordinates as:

τ(x,y)=12x2+y2[Newtons per square centimeter, N/cm2]\tau(x, y) = 12\sqrt{x^2 + y^2} \quad [\text{Newtons per square centimeter, N/cm}^2]τ(x,y)=12x2+y2​[Newtons per square centimeter, N/cm2]
where x2+y2≤9x^2 + y^2 \le 9x2+y2≤9.

Step-by-Step Solution:

  1. Define Polar Domain of Integration:
    State the Circle Polar Radius Identity:

    x2+y2≤R2  ⟹  r2≤32  ⟹  0≤r≤3 cmx^2 + y^2 \le R^2 \implies r^2 \le 3^2 \implies 0 \le r \le 3\text{ cm}x2+y2≤R2⟹r2≤32⟹0≤r≤3 cm
    Because the shaft is a complete solid cylinder, the angular sweep traverses the full circle:
    0≤θ≤2π rad0 \le \theta \le 2\pi\text{ rad}0≤θ≤2π rad
    Thus, the integration domain is D={(r,θ)∣0≤r≤3, 0≤θ≤2π}D = \{(r, \theta) \mid 0 \le r \le 3, \, 0 \le \theta \le 2\pi\}D={(r,θ)∣0≤r≤3,0≤θ≤2π}.

  2. Transform Integrand and Area Differential:
    State the Polar Coordinate Radial Distance Identity:

    x2+y2=r\sqrt{x^2 + y^2} = rx2+y2​=r
    Substitute into the shear stress function:
    τ(r)=12r N/cm2\tau(r) = 12r\text{ N/cm}^2τ(r)=12r N/cm2
    The internal moment arm multiplying the shear stress is rrr. Incorporating the Polar Area Element Formula with Jacobian (dA=r dr dθdA = r \, dr \, d\thetadA=rdrdθ):
    dT=(r)⋅τ(r) dA=(r)(12r)(r dr dθ)=12r3 dr dθdT = (r) \cdot \tau(r) \, dA = (r)(12r)(r \, dr \, d\theta) = 12r^3 \, dr \, d\thetadT=(r)⋅τ(r)dA=(r)(12r)(rdrdθ)=12r3drdθ
    Formulate the double integral for total torque TTT:
    T=∬D12r3 dr dθ=∫02π∫0312r3 dr dθT = \iint_D 12r^3 \, dr \, d\theta = \int_{0}^{2\pi} \int_{0}^{3} 12r^3 \, dr \, d\thetaT=∬D​12r3drdθ=∫02π​∫03​12r3drdθ

  3. Evaluate the Double Integral:
    Apply Fubini's Separability Theorem in Polar Coordinates (defined in Separability in Polar Integrals):

    T=(∫02πdθ)(∫0312r3 dr)T = \left( \int_{0}^{2\pi} d\theta \right) \left( \int_{0}^{3} 12r^3 \, dr \right)T=(∫02π​dθ)(∫03​12r3dr)

    Evaluate the radial integral:

    ∫0312r3 dr=[3r4]03=3(3)4−0=3(81)=243 N⋅cm\int_{0}^{3} 12r^3 \, dr = \left[ 3r^4 \right]_{0}^{3} = 3(3)^4 - 0 = 3(81) = 243\text{ N}\cdot\text{cm}∫03​12r3dr=[3r4]03​=3(3)4−0=3(81)=243 N⋅cm

    Evaluate the angular integral:

    ∫02πdθ=[θ]02π=2π rad\int_{0}^{2\pi} d\theta = [\theta]_0^{2\pi} = 2\pi\text{ rad}∫02π​dθ=[θ]02π​=2π rad

    Multiply and convert units:

    T=2π×243=486π N⋅cm≈1526.814 N⋅cmT = 2\pi \times 243 = 486\pi\text{ N}\cdot\text{cm} \approx 1526.814\text{ N}\cdot\text{cm}T=2π×243=486π N⋅cm≈1526.814 N⋅cm
    Convert to standard SI units (1 m=100 cm1\text{ m} = 100\text{ cm}1 m=100 cm):
    T=486π100 N⋅m=4.86π N⋅m≈15.268 N⋅mT = \frac{486\pi}{100}\text{ N}\cdot\text{m} = 4.86\pi\text{ N}\cdot\text{m} \approx 15.268\text{ N}\cdot\text{m}T=100486π​ N⋅m=4.86π N⋅m≈15.268 N⋅m

Meaning of the Answer & Real-Life Application:

  • Quantitative Interpretation: The total internal elastic resisting torque generated across the cross-section is 486π N⋅cm≈15.27 N⋅m486\pi\text{ N}\cdot\text{cm} \approx 15.27\text{ N}\cdot\text{m}486π N⋅cm≈15.27 N⋅m.
  • Engineering Application: In mechanical driveline design, integrating shear stress over polar coordinates confirms the classical torsion formula T=τmax⁡JRT = \frac{\tau_{\max} J}{R}T=Rτmax​J​. For a solid circular shaft of radius R=3 cmR = 3\text{ cm}R=3 cm, the polar moment of inertia is J=π2R4=π2(81)=40.5π cm4J = \frac{\pi}{2} R^4 = \frac{\pi}{2}(81) = 40.5\pi\text{ cm}^4J=2π​R4=2π​(81)=40.5π cm4. At the outer boundary (r=3 cmr = 3\text{ cm}r=3 cm), τmax⁡=12(3)=36 N/cm2\tau_{\max} = 12(3) = 36\text{ N/cm}^2τmax​=12(3)=36 N/cm2, giving T=(36)(40.5π)3=486π N⋅cmT = \frac{(36)(40.5\pi)}{3} = 486\pi\text{ N}\cdot\text{cm}T=3(36)(40.5π)​=486π N⋅cm. Knowing the exact resisting torque enables engineers to size shaft diameters to avoid plastic deformation under dynamic starting torque surges.