Axisymmetric torsion of a bar

You can use a 1D torsional model if you compute a “torsional constant” following this link. The 1D torsional model is quite simple (only a function of rotation of the cross section).

Additionally, a 3D implementation exists here at this link