# An Introduction to Uniform Torsion in Cylindrical Shafts

### The Fundamental Engineering Mechanics Series – Part 3

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.

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).

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 , 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 .

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

Like all good engineering derivations, we now consider a small segment of the bar with length , pictured below. We note that the right cross-section rotates relative to the left cross-section by an amount . As a result, the original element becomes with rotating to and rotating to .

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 and remain unchanged. From this, we can conclude that the deformation of the element is completely defined by the decrease in the angle . We recognise this as a shear deformation, , 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.

We can see from Fig 4(a) above that the magnitude of the shear strain is given by,

(1)

We note that for small angles we can replace the distance with , yielding,

(2)

The quantity is the **rate of twist or twist per unit length** and is replaced by ,

(3)

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

(4)

Note that for elements of length experiencing pure torsion (constant torque magnitude throughout the bar), and **for pure torsion only** we have,

(5)

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 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.

## 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, . 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 and strain are related by the shear modulus of elasticity, ,

(6)

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

(7)

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

(8)

where 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.

**3.1 Relating shear stress with torque**

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

(9)

Replacing with the expression we derived previously yields,

(10)

Note that the summation of all of the moments , acting on all elements of area , equals the applied torque, . We can perform this summation by integrating over the full cross-sectional area,

(11)

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

(12)

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

(13)

Therefore we have,

(14)

Rearranging, we obtain the torsion formula,

(15)

If we recall our expression from above, where was the rate of twist or twist per unit length, we can combine this with the torsion formula to obtain,

(16)

In the above equation, the quantity is known as the torsional rigidity and is analogous to the flexural rigidity 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 and has the following parameter values:

- diameter
- length
- shear modulus

Use this information to:

- Determine the maximum shear stress that develops in the bar
- Determine the angle of twist between the bar ends
- If the maximum allowable shear stress is and the maximum allowable angle of twist is , what is the maximum allowable torque?

**Solution**

(1) The maximum shear stress is obtained using,

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.

1 2 3 4 5 6 7 8 9 10 11 12 13 |
import math #Basic math functionality import numpy as np #Numpy for working with arrays G = 75*10**9 #(N/m^2) L = 2 #(m) T = 800 #(Nm) d = 0.05 #(m) r = d/2 #(m) Ip = math.pi*d**4/32 #(m^4) tau_max = T*r/Ip #(N/m^2) print(f'The maximum shear stress is {round(tau_max/10**6,1)} N/mm^2') |

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

(2) The angle of twist is obtained using,

1 2 3 4 5 |
phi_radians = T*L/(G*Ip) #(rads) phi_deg = phi_radians*180/math.pi #(deg) print(f'The angle of twist between ends is {round(phi_radians,3)} radians or {round(phi_deg,2)} degrees') |

*-> 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, . In this case, we can solve for the maximum torque by rearranging the torsion formula as follows,

1 2 3 4 5 |
tau_max = 35 #(N/mm^2) T_max_1 = (tau_max*10**6)*Ip/r print(f'The maximum torque limited by shear stress is {round(T_max_1,1)} Nm') |

*-> 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,

1 2 3 4 5 6 |
phi_max_deg = 1.75 #(degrees) phi_max_rad = phi_max_deg*math.pi/180 #(rads) T_max_2 = phi_max_rad*G*Ip/L print(f'- The maximum torque limited by angle of twist is {round(T_max_2,1)} Nm') |

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

We can see that the maximum allowable torque is 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.

1 2 3 4 5 6 7 8 9 |
#Programmatically test which is the limiting factor if(T_max_1<T_max_2): print(f'(c) The maximum allowable torque is {round(T_max_1,1)} Nm, limited by shear stress') elif(T_max_1>T_max_2): print(f'(c) The maximum allowable torque is {round(T_max_2,1)} Nm, limited by angle of twist') else: print(f'(c) The allowable torque is {round(T_max_1,1)} Nm and simultaneously limited by shear stress and angle of twist') |

*-> 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 while the hollow shaft has an inner diameter of . Both shafts experience a torque . 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 . Therefore we simply need to compare the ratio of values to obtain a ratio for shear stress for the shafts. For the solid shaft,

For the hollow shaft,

1 2 3 4 5 6 7 8 9 10 11 |
D_outer = 1 #Value representing outer diameter D_ratio = 0.75 #Ratio of diameters, D_inner/D_outer Ip_solid = (math.pi*D_outer**4)/2 Ip_hollow = Ip_solid - (math.pi*(D_ratio*D_outer)**4)/2 tau_ratio = (1/Ip_hollow)/(1/Ip_solid) print(f'(a) The ratio of maximum shear stress in the hollow shaft to that on the solid shaft is {round(tau_ratio,2)}') print(f'In other words the hollow shaft experiences a {round(100*(tau_ratio-1))}% larger maximum stress for the same applied torque') |

*-> (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…

1 2 3 4 |
print(f'(b) The ratio of angle of twist in the hollow shaft to that on the solid shaft is also {round(tau_ratio,2)}') print(f'Again, the hollow shaft experiences {round(100*(tau_ratio-1))}% more rotation for the same applied torque') |

*-> (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.

1 2 3 4 5 6 7 |
A_solid = math.pi*D_outer**2 A_hollow = A_solid - math.pi*(D_ratio*D_outer)**2 weight_ratio = A_hollow/A_solid print(f'(c) The ratio of weight of the hollow shaft to that of the solid shaft is {weight_ratio}') print(f'The hollow shaft is {round(100*(1-weight_ratio))}% lighter than the equivalent diameter solid shaft') |

*-> (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.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 |
#Outer diameter values are constant D_outer = 1 Ip_solid = (math.pi*D_outer**4)/2 A_solid = math.pi*D_outer**2 #Containers to hold our calculated ratios tr = np.array([]) wr = np.array([]) R = np.arange(0.1,0.81,0.01) #A selection of diameter ratios for r in R: D_ratio = r #Ratio of diamters, D_inner/D_outer Ip_hollow = Ip_solid - (math.pi*(D_ratio*D_outer)**4)/2 A_hollow = A_solid - math.pi*(D_ratio*D_outer)**2 tau_ratio = (1/Ip_hollow)/(1/Ip_solid) #Shear stress ratio tr = np.append(tr, tau_ratio) #Save for later weight_ratio = A_hollow/A_solid #Weight ratio wr = np.append(wr, weight_ratio) #Save for later |

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

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 |
import matplotlib.pyplot as plt #Plotting functionality fig, ax1 = plt.subplots(figsize=(12,6)) ax2 = ax1.twinx() ax1.plot(R,tr,'b') ax2.plot(R,wr,'r') #Housekeeping ax1.set_xlabel('Diameter ratio (Di/Do)') ax1.set_ylabel('Stress ratio', color='b') ax2.set_ylabel('Weight ratio',color='r') ax1.set_title('Diameter ratio versus stress and weight ratios') ax1.set_xlim([0.1,1]) plt.grid() plt.show() |

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 without exceeding an allowable shear stress of or an allowable rate of twist of . The shear modulus of elasticity is . A hollow or solid shaft may be used. To assess the suitability of each option:

- determine the required diameter of the solid shaft,
- determine the required outer diameter of the hollow shaft, , if the wall thickness is of the outer diameter
- determine the diameter ratio between shafts, 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,

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

(17)

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

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 |
#CONSTANTS T = 5000 #(Nm) tau_max = 50*10**6 #(N/m^2) theta_max = 0.45 #(deg/m) theta_max_rad = theta_max*math.pi/180 #(rads) G = 75*10**9 #(N/m^2) #SOLID SHAFT - LIMITED BY ALLOWABLE SHEAR STRESS r = ((2*T)/(math.pi*tau_max))**(1/3) d_solid_1 = 2*r #(m) Solid diameter limited by shear stress #SOLID SHAFT - LIMITED BY ALLOWABLE RATE OF TWIST r = ((2*T)/(math.pi*G*theta_max_rad))**(1/4) d_solid_2 = 2*r #(m) Solid diameter limited by rate of twist #Determine which condition governs if(d_solid_1>d_solid_2): d_solid = d_solid_1 print(f'(a). The req. dia. of the solid shaft is {round(d_solid_1,3)} m, limited by allowable shear stress') elif(d_solid_1<d_solid_2): d_solid = d_solid_2 print(f'(a). The req. dia. of the solid shaft is {round(d_solid_2,3)} m, limited by allowable rate of twist') else: d_solid = d_solid_1 print(f'(a). The req. dia. of the solid shaft is {round(d_solid_1,3)} m, simultaneously limited by allowable shear stress and rate of twist') |

*-> (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, . Therefore,

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 |
#HOLLOW SHAFT - LIMITED BY ALLOWABLE SHEAR STRESS d_hollow_1 = ((0.5*T)/(0.058*tau_max))**(1/3) #(m) Hollow outer diameter limited by shear stress #HOLLOW SHAFT - LIMITED BY ALLOWABLE RATE OF TWIST d_hollow_2 = (T/(G*theta_max_rad*0.058))**(1/4) #Determine which condition governs if(d_hollow_1>d_hollow_2): d_hollow = d_hollow_1 print(f'(b). The req. outer dia. of the hollow shaft is {round(d_hollow_1,3)} m, limited by allowable shear stress') elif(d_solid_1<d_solid_2): d_hollow = d_hollow_2 print(f'(b). The req. outer dia. of the hollow shaft is {round(d_hollow_2,3)} m, limited by allowable rate of twist') else: d_hollow = d_hollow_1 print(f'(b). The req. outer dia. of the hollow shaft is {round(d_hollow_1,3)} m, simultaneously limited by allowable shear stress and rate of twist') |

*-> (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.

1 2 3 4 5 6 7 8 |
diam_ratio = d_hollow/d_solid #(diameter ratio) A_solid = math.pi*d_solid**2/4 #(m^2) Cross-sectional area of solid shaft d_i = 0.8*d_hollow #(m) inner diameter of hollow bar A_hollow = math.pi*d_hollow**2/4 - math.pi*d_i**2/4 #(m^2) Cross-sectional area of hollow shaft weight_ratio = A_hollow/A_solid #Weight ratio (=area ratio) print(f'c. The dia. ratio of the hollow to solid shafts is {round(diam_ratio,2)} while the weight ratio is {round(weight_ratio,2)}') |

*-> 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 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, , and the angle of twist per unit length ,

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

- We derived an equation that relates the rate of twist per unit length, , with the applied torque, and torsional rigidity, .

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.

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.