An Introduction to Uniform Torsion in Cylindrical Shafts

The Fundamental Engineering Mechanics Series – Part 3

by Dr Seán Carroll
|
Published: 15th September 2022
|
Tutorial

Welcome to part 3 of the Fundamental Engineering Mechanics tutorial series. This part will focus on torsion in circular shafts and the stresses and strains it induces. For the previous instalment in the series, go here.

In the previous tutorial in this series, we covered free-body diagrams and how to calculate reaction forces. It may seem like quite a jump to now be exploring torison in part 3. This is because, before torsion, I would usually suggest covering truss analysis and shear force and bending moment analysis for beams and frames. But I’ve already covered these in multiple tutorials, so there’s no need to cover them again as part of this series. Feel free to cover those before diving into Torsion with me here.

Remember that you can access any video and downloadable versions of the notes by joining the accompanying course below.

1.0 Introducing torsion

As mentioned above, in previous tutorials, we’ve discussed trusses and observed that the members that make up a truss experience only axial force. We could say that the action experienced by these members is axial force. We’ve also discussed structural elements whose primary actions were shear forces and bending moments (beams and frames).

We should clarify here that what distinguishes an axial force from a shear force is simply the axis along which the force is applied. Axial forces have their line of action parallel to the longitudinal axis of the member, while shear forces have a line of action perpendicular to this axis.

When discussing bending moments applied to 2D structures in the x-y plane, we didn’t labour the point at the time, but the moments considered were all about an axis perpendicular to the x-y plane (page or screen) – the z-axis, see Fig 1 below.

Intro to torsion 1 | DegreeTutors.com
Fig 1. Beam subject to shearing forces (blue), axial force (green) and bending moment (pink). All of these actions have been discussed in previous tutorials.

In this tutorial, we introduce another action – torsion. This is simply a bending moment applied about a different axis – the longitudinal axis, Fig 2. Torsion will cause twisting about the longitudinal axis and is a very common form of loading that we need to be able to analyse. We often refer to torsional moments as torques (we can use the terms torsional moment and torque interchangeably).

Intro to torsion 2 | DegreeTutors.com
Fig 2. Beam subject to a torsional moment about the longitudinal x-axis.

We can represent torsion (or any moment) with the curved arrow you’ll be familiar with from previous tutorials or with a double-headed straight arrow. The double-headed arrow is typically more convenient for representing torsion applied to 2D structures. To determine the sense of rotation indicated by a double-headed arrow, we use the right hand screw rule.

According to this sign convention, if you point the thumb or your right hand in the direction of the double-headed arrow and make a fist, your fingers will curve around and point in the direction of the positive sense of rotation for the torque. This will be an important convention to understand and keep in mind.

2.0 Torsional Deformation

Our starting point here will be to explore the concept of strain as it applies to circular bars and to derive an equation that relates the strain to the angle of twist in the bar. First, let’s state some assumptions that will apply to our discussion:

  • the bar or shaft is circular in cross-section and prismatic, which means it has the same cross-section shape throughout its length.
  • the angle of twist between both ends of the bar, which we’ll denote as \phi, remains small.
  • the applied torque is constant throughout the length of the bar – we therefore refer to the bar as being in pure torsion.
  • angles of twist are measured in radians, recall that \pi\: \text{radians}=180^\circ.

Consider the bar with cross-sectional radius r and length L, pictured below. It’s fixed against torsional rotation at end p. The applied torque at q, T, induces an angle of twist \phi such that point q rotates around to point q' at the free end of the bar. We also note that the angle of twist varies linearly between x=0 and x=L.

Intro to torsion 3 | DegreeTutors.com
Fig 3. Prismatic circular bar subject to torque T, showing resulting angular deformation \phi(x).

Like all good engineering derivations, we now consider a small segment of the bar with length \mathrm{d}x, pictured below. We note that the right cross-section rotates relative to the left cross-section by an amount \mathrm{d}\phi. As a result, the original element abcd becomes ab'c'd with b rotating to b' and c rotating to c'.

Because we’ve made the initial assumption that the angle of rotation between the ends of the bar is small, we can say that the lengths of the sides ab' and dc' remain unchanged. From this, we can conclude that the deformation of the element abcd is completely defined by the decrease in the angle \angle bad. We recognise this as a shear deformation, \gamma, and conclude that the element is in pure shear which means it is only experiencing shear strain, as opposed to normal strains that would be induced by axial forces or bending moments.

Intro to torsion 4 | DegreeTutors.com
Fig 4. (a) Prismatic circular bar in a state of pure torsion. (b) Sub-length \mathrm{d}x experiencing an angle of twist \mathrm{d}\phi giving rise to a shear deformation \gamma_{\text{max}} on the outer surface of the bar.

We can see from Fig 4(a) above that the magnitude of the shear strain \gamma_{\text{max}} is given by,

(1)   \begin{equation*}\gamma_{\text{max}} = \frac{bb'}{\mathrm{d}x}\end{equation*}

We note that for small angles we can replace the distance bb' with r\:\mathrm{d}\phi, yielding,

(2)   \begin{equation*}\gamma_{\text{max}} = \frac{r\:\mathrm{d}\phi}{\mathrm{d}x}\end{equation*}

The quantity \mathrm{d}\phi/\mathrm{d}x is the rate of twist or twist per unit length and is replaced by \theta,

(3)   \begin{equation*}\theta = \frac{\mathrm{d}\phi}{\mathrm{d}x}\end{equation*}

As a result, we obtain the equation for the shear strain at radius r,

(4)   \begin{equation*}\boxed{\gamma_{\text{max}} = r\theta}\end{equation*}

Note that for elements of length L experiencing pure torsion (constant torque magnitude throughout the bar), \theta = \phi/L and for pure torsion only we have,

(5)   \begin{equation*}\boxed{\gamma_{\text{max}} = \frac{r\phi}{L}}\end{equation*}

We can see from the preceding equations that for any given cross-section the shear strain varies linearly with the radius. For hollow circular tubes, for example, we can obtain the shear strain magnitude at the inner and outer walls simply by substituting the relevant value for the radius, (see below). Therefore, the max in the subscript \gamma_{\text{max}} refers to the fact that this shear strain was evaluated for the outer surface of the bar, where the radius and therefore shear strain are a maximum.

Intro to torsion 4 | DegreeTutors.com
Fig 5. Shear strain \gamma varies linearly with radius.

3.0 Torsion in linearly elastic circular bars

So, at this point, we have an equation that relates the angle of twist with the associated shear strain. Next, we’d like to tie this into the shear stress, \tau. If the bar is made of a linearly elastic material, we know the shear stress will be linearly proportional to the shear strain. In the case of shear (as opposed to tension and compression), the stress (\tau) and strain (\gamma) are related by the shear modulus of elasticity, G,

(6)   \begin{equation*}\tau = G\gamma\end{equation*}

We can therefore state the shear stress, \tau at a radius r as,

(7)   \begin{equation*}\boxed{\tau = G r \theta}\end{equation*}

If we wish to determine the shear stress at some internal radius \rho, we could do so as follows,

(8)   \begin{equation*}\tau = \frac{\rho}{r}\tau_{\text{max}}\end{equation*}

where \tau_{\text{max}} is the maximum shear stress on the outer surface.

Figure 6(a) below shows a bar in pure torsion and an enlarged stress element with the directions of the resulting shear stresses at the surface indicated. In image (b), we can see the shear stresses and corresponding shear deformation on the element. We again see in image (c) that the shear stress varies linearly with radius.

Intro to torsion 6 | DegreeTutors.com
Fig 6. (a) Bar in pure torsion, (b) resulting stresses and strains reflected on an element taken from the bar, (c) cross-section of the bar showing linear variation of shear stress with radius.

3.1 Relating shear stress with torque

Now we’ll derive the torsion formula that relates the applied torque T with the shear stress induced, \tau. We start by considering the stressed cross-section and an element of area within this cross-section \mathrm{d}A (see Fig 7 below). This infinitesimal area is located at radius \rho and experiences a shear stress, \tau. Therefore we can calculate the force experienced by \mathrm{d}A as \tau \mathrm{d}A. This force generates a moment about the axis of the bar, \mathrm{d}M.

(9)   \begin{equation*}\mathrm{d}M = \tau \rho \mathrm{d}A\end{equation*}

Replacing \tau with the expression we derived previously yields,

(10)   \begin{equation*}\mathrm{d}M = \frac{\tau_{\text{max}}}{r} \rho^2 \mathrm{d}A\end{equation*}

Note that the summation of all of the moments \mathrm{d}M, acting on all elements of area \mathrm{d}A, equals the applied torque, T. We can perform this summation by integrating \mathrm{d}M over the full cross-sectional area,

(11)   \begin{equation*}T = \int_A \mathrm{d}M = \frac{\tau_{\text{max}}}{r}\int_A \rho^2 \mathrm{d}A\end{equation*}

We now define the polar moment of inertia of the circular cross-section, I_P,

(12)   \begin{equation*}I_P = \int_A \rho^2 \mathrm{d}A\end{equation*}

This is a property of the cross-section shape (analogous to the second moment of area, I, which we encounter when considering flexural rigidity, EI in bending). For a circle with diameter d, the polar moment of inertia is given by,

(13)   \begin{equation*}I_P = \frac{\pi d^4}{32}\end{equation*}

Therefore we have,

(14)   \begin{equation*}T = \frac{\tau_{\text{max}}}{r}I_P\end{equation*}

Rearranging, we obtain the torsion formula,

(15)   \begin{equation*}\boxed{\tau_{\text{max}} = \frac{T r}{I_p}}\end{equation*}

Intro to torsion 7 | DegreeTutors.com
Fig 7. Cross-section of a bar in pure torsion showing an element of area \mathrm{d}A experiencing a torsional stress \tau.

If we recall our expression from above, \tau_{\text{max}} = Gr\theta where \theta was the rate of twist or twist per unit length, we can combine this with the torsion formula to obtain,

(16)   \begin{equation*}\theta = \frac{T}{GI_P}\end{equation*}

In the above equation, the quantity GI_P is known as the torsional rigidity and is analogous to the flexural rigidity EI for bending.

Limitations

Our discussion and derivations so far have assumed the material remains linearly elastic and at no time is the applied torque such that the proportional limit for the material is exceeded. Our derivations are also only valid in the absence of stress concentrations that result from discontinuities in geometry or at points of load application. Next, we can consider some numerical examples to demonstrate these equations in action.

3.2 Example 1

Consider a solid steel bar with a circular cross-section. The bar is subject to an applied torque T=800\:Nm and has the following parameter values:

  • diameter d = 50\:mm
  • length L=2.0\:m
  • shear modulus G = 75\:kN/mm^2

Use this information to:

  1. Determine the maximum shear stress that develops in the bar
  2. Determine the angle of twist between the bar ends
  3. If the maximum allowable shear stress is 35 \:N/mm^2 and the maximum allowable angle of twist is 1.75^\circ, what is the maximum allowable torque?

Solution

(1) The maximum shear stress is obtained using,

    \begin{equation*}\tau_{\text{max}} = \frac{T r}{I_p}\end{equation*}

We’ll use some simple Python to perform the calculations but you can easily work this out with a pen and paper if you prefer.

-> The maximum shear stress is 32.6 N/mm^2

(2) The angle of twist is obtained using,

    \begin{equation*}\phi = \frac{TL}{GI_P}\end{equation*}

-> The angle of twist between ends is 0.035 radians or 1.99 degrees

(3) Let’s first assume the torque is limited by the maximum allowable shear stress, \tau_{\text{max}}=35\:N/mm^2. In this case, we can solve for the maximum torque by rearranging the torsion formula as follows,

    \begin{equation*}T = \frac{\tau_{\text{max}}I_p}{r}\end{equation*}

-> The maximum torque limited by shear stress is 859.0 Nm

Now we can test what the maximum torque would be if the angle of twist were the limiting criteria using,

    \begin{equation*}T = \frac{\phi G I_p}{L}\end{equation*}

-> The maximum torque limited by angle of twist is 702.8 Nm

We can see that the maximum allowable torque is 702.8\:Nm and is limited by the allowable angle of twist specified. However, keeping with the Python theme here and so that we have a reusable code that yields the final answer without our intervention at the end, we’ll write some code to output the final conclusion.

-> The maximum allowable torque is 702.8 Nm, limited by angle of twist

3.3 Example 2

Consider two shafts, one hollow and one solid. Both have an outer diameter D while the hollow shaft has an inner diameter of 0.75D. Both shafts experience a torque T. Our task is to analyse the relative efficiency of each shaft at transmitting torque. We’ll do this by establishing the ratio of the following properties for each shaft:

  • shear stress
  • angle of twist
  • weight

Solution

(a) If the torque and outer radius are the same, the shear stress is proportional to 1/I_P. Therefore we simply need to compare the ratio of I_P values to obtain a ratio for shear stress for the shafts. For the solid shaft,

    \begin{equation*}I_P = \frac{\pi D^4}{2}\end{equation*}

For the hollow shaft,

    \begin{equation*}I_P = \frac{\pi D^4}{2} - \frac{\pi (0.75D)^4}{2}\end{equation*}

-> (a) The ratio of maximum shear stress in the hollow shaft to that on the solid shaft is 1.46
-> (b) In other words the hollow shaft experiences a 46% larger maximum stress for the same applied torque

(b) Because both shafts are the same length and made of the same material, the angles of twist are also proportional to the inverse of the polar moment of inertia, therefore…

-> (b) The ratio of angle of twist in the hollow shaft to that on the solid shaft is also 1.46
-> Again, the hollow shaft experiences 46% more rotation for the same applied torque

(c) The weights of the shafts are proportional to their cross-sectional areas.

-> (c) The ratio of weight of the hollow shaft to that of the solid shaft is 0.4375
-> The hollow shaft is 56% lighter than the equivalent diameter solid shaft

From this analysis, we can see the inherent efficiency offered by a hollow shaft. Although it develops a larger magnitude shear stress under the same applied torque and a larger rotation, it is significantly lighter. So, if the material the shaft is made from can withstand the larger stress and if the larger angle of twist can be accommodated, it makes more sense to select a hollow shaft.

We can get a better sense of how these relationships change by plotting the ratios for different values of diameter ratio. Since we’ve already written the code to perform the calculation for one diameter ratio, all we need to do is wrap this in a for loop and save the values calculated on each iteration.

Now that we have all of the values calculated and saved, we can plot them so see the behaviour more clearly.

Intro to torsion 8 | DegreeTutors.com
Gig 8. Plot of stress ratio and weight ratio against ratio of inner and outer shaft diameters.

The blue line shows us how the hollow tube generates more stress when compared to the solid tube, as the inner diameter approaches the outer diameter and the walls of the tube get thinner. At the same time, we can see from the red line how the weight of the hollow tube progressively reduces. Note that the cross-over point here does not have any special significance as we have two different y-scales.

3.4 Example 3

A circular steel shaft must transmit a torque of 5000\:Nm without exceeding an allowable shear stress of 50\:N/mm^2 or an allowable rate of twist of 0.45^\circ/m. The shear modulus of elasticity is 75\: GPa. A hollow or solid shaft may be used. To assess the suitability of each option:

  1. determine the required diameter of the solid shaft, d
  2. determine the required outer diameter of the hollow shaft, d_o, if the wall thickness is 10\% of the outer diameter
  3. determine the diameter ratio between shafts, d_o/d and the ratio of weights of the hollow and solid shafts

Solution

For the solid and hollow shaft, the diameter may be limited by either the allowable shear stress or allowable rate of twist. We need to determine diameters based on both criteria and select the largest as the minimum diameter for each shaft. To evaluate a diameter limited by the allowable stress we will use,

    \begin{equation*}\tau_{\text{max}} = \frac{T r}{I_p}\end{equation*}

To evaluate a diameter limited by the allowable rate of twist we will use,

(17)   \begin{equation*}\theta = \frac{T}{GI_p}\end{equation*}

In both cases, we can rearrange and solve to determine the diameter.

-> (a). The req. dia. of the solid shaft is 0.096 m, limited by allowable rate of twist

For the hollow shaft we note the inner diameter, d_i = 0.8\times d_o. Therefore,

    \begin{align*}I_P &= \frac{\pi}{32}\left(d_o^4 - (0.8d_o)^4\right)\\I_P &= \frac{\pi}{32}\left(0.5904 \:d_o^4\right)\\I_P &= 0.058\:d_o^4\end{align*}

-> (b) The req. outer dia. of the hollow shaft is 0.11 m, limited by allowable rate of twist

Finally, we can determine the diameter and weight ratios as follows.

-> The dia. ratio of the hollow to solid shafts is 1.14 while the weight ratio is 0.47

From this we conclude that although the hollow shaft would have a 14\% larger diameter, it would be less than half the weight of a solid shaft.

4.0 Summary

Let’s quickly recap what we established in this tutorial:

  • Torque or torsion is just a bending moment about the longitudinal axis
  • We established a relationship between shear strain on the surface of a circular shaft, \gamma_{max}, and the angle of twist per unit length \theta,

    \begin{equation*}\gamma_{\text{max}} = r\theta\end{equation*}

  • We derived the torsion formula, which relates the shear stress on the surface of a circular shaft with the applied Torque, T, shaft radius, r, and polar moment of inertia, I_p,

    \begin{equation*}\tau_{\text{max}} = \frac{T r}{I_p}\end{equation*}

  • We derived an equation that relates the rate of twist per unit length, \theta, with the applied torque, T and torsional rigidity, GI_p.

Everything we discussed assumed that:

  • We were dealing with prismatic circular shafts
  • The shaft material remains linearly elastic
  • The shafts were in pure torsion
  • The rotational deformation of the shaft remained small

In the next tutorial in the series, we’ll expand on our discussion of torsion and consider non-uniform torsion in circular shafts. If you want to be notified when the next instalment is published, join the free Fundamentals of Mechanics course by signing up using the link at the top of this tutorial.

Author

Dr Seán Carroll

BEng (Hons), MSc, PhD, CEng MIEI, FHEA

Hi, I’m Seán, the founder of DegreeTutors.com. I hope you found this tutorial helpful. After spending 10 years as a university lecturer in structural engineering, I started DegreeTutors.com to help more people understand engineering and get as much enjoyment from studying it as I do. Feel free to get in touch or follow DegreeTutors on any of the social accounts.


If you want full access to the complete DegreeTutors course (and code!) library, for one low annual fee, consider subscribing to the DegreeTutors All Access Membership. You can read more about the benefits of Membership here or skim the headlines and subscribe directly below. If you’re a student at a college or university, become a DegreeTutors student member to access a selection of premium courses, completely free.

Access Every Course and Tool

$279 / YEAR

(less than a coffee per week ☕)

Price Guarantee – avoid future price rises as we grow
Over 740 lectures & over 140 hours of HD video content
 Unlimited access to online analysis toolboxes. 🛠️
Packed development roadmap of courses & tools 🏗️
Access to DegreeTutors:Labs member’s hub 🧪
Unlimited streaming of all current and future courses
Early access to new courses
Access all downloads, pdf guides & Python codes
Course completion certificates
Priority Q&A support