Torsion of a Shaft with Circular Cross-Section

Let us consider a shaft with circular cross-section of radius R. Length of shaft is L (see figure).

Torsion of a Shaft with Circular Cross-Section

Select the coordinate system with the z-axis directed along the axis of the shaft, and the coordinate z=0 located at the left edge of the shaft.
The shaft is subjected to the externally applied torque τ. The torque is applied at the right end of the shaft, the left end of the shaft is rigidly clamped.

Torsion of a Shaft with Circular Cross-Section, the finite element model with applied loads and restraints

The finite element model with applied loads and restraints

Let us use the following initial data: length L of the shaft is 0.6 m, radius of cross-section R of the shaft is 0.02 m, the magnitude of the applied torque τ is 100 N-m.
Material characteristics: E = 2.1E+011 Pa, ν = 0.28.
To find the angle of twist, let us use the following relation:
,
where ϕ0 – angle of twist of the cross-section z=0, G=E/2(1+ν) – shear modulus, Jp=πR4 / 2 – polar moment of inertia of the circular cross-section.
Since by formulation of the task, the left end of the shaft is clamped, ϕ0=0. Then, at a distance z=0.5L from the clamped edge of the shaft, the angle of twist ϕ is given by the formula:

Thus, ϕ0.5L= 1.4551E-003 rad.
After carrying out calculation with the help of AutoFEM, the following results are obtained:

Table 1.Parameters of finite element mesh

Finite Element Type

Number of nodes

Number of Finite Elements

quadratic tetrahedron

3270

11237

Absolute value of displacement (at z=0.5L) Δu = 2.9488E-005 m. (mean at Sensor 1-4).

Table 2. Result "Angle of twist" *

Numerical Solution
Angle of twist Ψ = arcsin(Δu / R) ,rad

Analytical Solution
Angle of twist ϕ, rad

Error δ = 100%* | Ψ - ϕ | / | ϕ |

1.4744E-003

1.4551E-003

1.33

Torsion of a Shaft with Circular Cross-Section, Result "Displacement, magnitude"

Conclusions:

The relative error of the numerical solution compared to the analytical solution is 1.33%.

*The results of numerical tests depend on the finite element mesh and may differ slightly from those given in the table.

 

Read more about AutoFEM Static Analysis

autofem.com

Return to contents