 Open Access
 Total Downloads : 448
 Authors : Mohammed Hjaji, Abdurahman Hassen, Hasan Nagiar, Ezedine Allaboudi
 Paper ID : IJERTV5IS090300
 Volume & Issue : Volume 05, Issue 09 (September 2016)
 Published (First Online): 15092016
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Exact Finite Element for Torsional Vibration of Shafts under Harmonic Torsion
Mohammed Hjaji1, Abdurahman Hassen2, Hasan Nagiar3 and Ezedine Allaboudi4 Applied Mechanics Division, Mechanical and Industrial Engineering Department University of Tripoli
Tripoli, Libya
Abstract The present paper investigate the dynamic analysis for the torsional vibration of shafts subjected to various harmonic twisting moments. The governing torsional vibration equation and related boundary conditions for the shafts are derived using Hamiltons variational principle. The exact closedform solutions for the shafts having a cantilevered and simply supported boundary conditions and under various harmonic torsional loads are determined. The closed form solution is consequently used to develop a family of exact shape functions which exactly satisfy the exact homogeneous solution of the governing torsional equation. A superconvergent twonoded finite beam element is then formulated based on the exact shape functions. The proposed beam element developed involves no special discretization errors normally encountered in conventional finite element formulations and provide results in excellent agreement with a minimal number of degrees of freedom. The present solutions are shown to successfully capture the static and steady state dynamic responses of shafts. Both solutions are also able to predict the natural torsional frequencies and related modeshapes of the shafts. The validity and the accuracy of the present exact closedform and finite beam element solutions are achieved throughout the numerical examples presented and compared with wellestablished ABAQUS finite beam solution and other exact solution.
Keywords Exact finite element, closed form solution, torsional response, harmonic torsional loading, steady state response

INTRODUCTION AND OBJECTIVE
Shafts and beams of circular crosssections are important structural elements of rotating and nonrotating machines and systems, mainly used to transmit torque and rotation. These shafts are often idealized as simple rotating and nonrotating beam models such as helicopter blades, aircraft rotary wings and gas turbine blades, etc. In such applications, these shafts are subjected to dynamic torsional loadings caused by machinery, environmental shocks, and dynamic exciting torque as repeated harmonic twisting moments in the long range duration, Therefore, the torsional vibration of shafts under harmonic torsional cyclic loading produce stresses reversal that affect the structural machine integrity and the life of components and then are prone to fatigue failure. Under harmonic twisting moments, the steady state dynamic component of the total response is sustained for a long time and is thus of great importance in fatigue design of such shafts. While the transient component of total response which is introduced only at the beginning of the excitation tends to dampen out quickly and is thus of a very little importance in assessing the fatigue life of the shaft. Thus, the objective of this paper is to develop an exact closed form solution and accurate and efficient finite element solution which captures and isolates the steady state torsional dynamic response of shafts subjected to general harmonic twisting moments. The present solutions developed in this study are also capable accurately to capture the quasistatic response and predict the eigenfrequencies and eigenmodes of the shaft.

LITERATURE REVIEW
Several studies have been conducted on the exact solution and finite element for torsional vibration analysis of shafts under various dynamic loads and boundary conditions. Among them, [1] derived the exact expressions for torsional frequencies and mode shapes for the free torsional vibrations of circular shafts and piping systems constrained by unsymmetrical torsional springs and carrying multiple unequal rotational masses at arbitrary locations on the shaft. Reference [2] developed a boundary element solution for the general linear elastic nonuniform torsion problem of homogeneous or composite prismatic bars of arbitrary cross section subjected to various twisting moment and linear torsional boundary conditions. In machinery systems, torsional vibration problems are important research topic of rotorshaft systems. Therefore, a lot of researchers have studied the free vibration frequencies of rotating shafts and beams. Rotating shaft is differed from the nonrotating one in having an additional centrifugal force and Coriolis effects on its dynamic behavior. For instance, [3] derived the finite element stiffness and mass matrices for both open and closed section shafts, based on the solution of the static equations, and thus obtained an approximate solution for the natural frequencies. Reference [4] introduced the classical Euler Bernoulli beam theory to present an analytical solution to investigate the vibrational behavior of a beam rotating with constant speed about its longitudinal axis for all boundary conditions. Reference [5] presented an analytical solution for free flexural vibration of a spinning finite Timoshenko beam subjected to a moving load for the general boundary conditions. Their formulation accounts for the rotary inertia and shear deformation effects. Reference [6] studied the dynamic response of a rotating shaft subject to axial force and moving loads is analyzed by using Timoshenko beam theory and the assumed mode method. Reference [7] presented an analytical approach to study the torsional vibrations of the drive shafts and mechanisms. Reference [8] studied the free vibrations and stability of internally damped rotating shafts with general boundary conditions. Reference [9] derived the exact governing equations for linear vibration of a rotating Timoshenko beam by using dAlembert principle and the virtual work principle. Based on the assumptions that the beam is linear elastic and the steady state axial strain is small, the formulation is captured the effect of Coriolis force on the natural frequency of the rotating Timoshenko beam. Reference [10] investigated the free vibrations analysis of a rotating shaft with nonlinearities in curvature and inertia.
Moreover, the finite element method is used for the torsional vibration analysis of solid and hollow shafts and
beams in determining the natural frequencies and mode shapes. This includes the work of [11] developed a finite element solution for determining the free vibration characteristics of rotating uniform Timoshenko beams. His
in which up (x,t) , vp (x,t) and wp (x,t) are the
longitudinal, lateral and transverse displacements of point p(x, y, z) along the principle axes ( X ,Y ,Z ) , (x,t) is
formulation incorporated the effects of shear deformation and rotary inertia on the natural frequencies of the rotating beams. Reference [12] derived explicit expressions for the finite
the torsional displacement, and point p along the principal axes.
y, z are the coordinates of
element mass and stiffness matrices using consistent mass formulation for the vibration of a rotating tapered beam. Reference [13] developed a three nodal CÂ° Timoshenko finite beam element to evaluate the natural whirling speeds of a rotating shaft with various end boundary conditions. Reference [14] presented a simple spinning composite shaft model based on a firstorder shear deformable beam theory. The finite element model is used to analyze the critical speeds, natural frequencies and related mode shapes of the composite shaft system. A feature common to the above finite element studies is use of approximate shape functions
Mx(,t)
O
Z
mx(x,t) Mx(L,t)
X
L
+(x,t)
Y O
involving spatial discretization errors, and thus requiring fine meshes to converge to the actual solution. In contrast, the present study avoids discretization errors by formulating exact shape functions which exactly satisfy the homogeneous solution of the dynamic governing torsional equation. Another commonality between the above studies is the fact they focus on extracting the free vibration characteristics
z P(y,z)
r
r
wp
vp
r
( , )
y
Z
including predicting the natural frequencies and mode shapes. In contrast, the present solution aims at directly extracting the steady state torsional dynamic response without the need to extracting the natural frequencies and mode shapes.

KINEMATICS FUNCTIONS
Fig.1: Coordinate system and displacements
StrainDisplacement Relations
Based on the assumption of small displacements, the nonzero shear strains are given by:
A straight uniform circular closed crosssection beam of
vp up z z
(4)
length L has is shown in Fig. (1). The shaft is referenced to a righthanded rectangular coordinates system ( X ,Y , Z ) ,
xy x y x
where the axis X is the longitudinal axis of the
wp up y y
(5)
beam, Y and Z are the principal axes of the crosssection. The present theoretical formulation is based on the following main assumptions:

The formulation is applicable to shafts having circular solid/hollow crosssections,

Crosssection is assumed to remain perfectly rigid in its own plane throughout deformation,

The material is assumed to remain linearly elastic throughout deformation,

Displacements, strains and rotations are assumed small, and

Damping is neglected in the formulation.

The formulation is captured only the steady state dynamic response.
xz x z x
Torsional Displacement Function
The shaft is assumed to be subjected to harmonic torsion;
(i) distributed twisting moment mx (x,t) mx (x) eit acting along the shaft axis, and (ii) concentrated twisting moments M (xe ,t) Mx (xe ) eit at the shaft both ends (i.e., e0,L ), as shown in Fig. (1). The harmonic twisting
moments can be written by:
mx x,t , Mx xe ,t mx x, Mx xe eit (6)
where is the frequency of the applied twisting moments,
According to the assumptions described above, the
i 1
is the value of the imaginary,
mx (x,t)
is the
displacement functions of an arbitrary point
harmonic distributed torsion, Mx (xe ,t) is the concentrated
p(x, y, z) located on the shaft crosssection, as shown in Figure (1), can be expressed as:
harmonic twisting moment applied at both ends of the given shaft. Under the effect of the applied harmonic twisting
up x,t 0
vp x,tr sin z (x,t)
(1)
(2)
moments, the angular displacement x,t
harmonic, i.e.,
is assumed to be
wp x,t r cos y (x,t)
(3)
x,t x xeit
(7)
in which x x represents the space function of the steady state torsional response. In line with the main objective of this paper focusing on steady state response, the torsional displacement function postulated in equation (7) neglect the transient component of the response.


TORSIONAL VIBRATION EQUATION
Exact Homogeneous Solution for Torsional Displacement
The exact homogeneous solution of the torsional equation is obtained by (a) setting the right hand side of the equation
(12) equal to zero, i.e., mx x 0 , and (b) assuming the
torsional displacement to take the following exponential form:
The variational form of the Hamiltonian functional H is
xh x Ai emi x
(14)
taken to be stationary, i.e.,
From equation (14) by substituting into homogeneous form of
H
t2 T U W dt 0
t1
t1
for x,t1 x,t2 0
(8)
equation (12), the exact homogeneous solution for torsional displacement is written in matrix form as:
A
xh x
em1x
em2x
1
E(x)
A
(15)
in which the integration is performed between arbitrary time
12 A2 21
12 21
limits t1 and t 2 , where T is the variation of kinetic energy,
defined by:
where
m1,2 i
G ,
E(x)
em1x
em2x ,
T L
u u v v w w dAdx
(9)
12 12
A A A are unknown integration constants.
0 A p p p p p p
12 1 2 12
U is the variation of internal strain energy, defined by:
Equation (15) governs the steady state dynamic solution for torsional vibration of shafts under harmonic twisting
L
L
A
A
U 0
G xy xy G xz xz dAdx
(10)
moments.
Particular Solution for Torsional Displacement
and W is the variation of the virtual work done for the shaft
For a shaft subjected to harmonic distributed twisting
subjected to given applied torsions, given as:
moment
mx x ei t =mx ei t , the corresponding particular
W
m x,t x,tdxM
x,t x,t
(11)
solution p x for the governing torsional vibration equation
L
L
0 x
x xxe
(12) is obtained by:
where is the density of the shaft material, and G is the modulus of rigidity of the shaft.
xp x mx
J 2
(16)
From Equations (17) and by substituting into energy
expressions (911), the resulting equations into Hamiltons equation (8), performing integration by parts, the governing equation for torsional vibration and related boundary
The complete steady state exact solution for torsional vibration response is obtained by adding the homogeneous solution (15) to the particular solution (16), gives:
conditions of shafts under harmonic twisting moments are given in terms of variable x as:
x x em1x em2x
A
mx
2
(17)
GJ xx J 2x x mx x
(12)
12
21
J
and the boundary conditions are:
Equation (17) represents the complete steady state solution for torsional vibration of shafts under harmonic torsional
L L L
loading, where the integration constants A can be
x x0 0 or GJx x0 Mx x0
(13)
21
determined from the relevant boundary conditions of the
where J is the torsional coefficient defined by;
problem.
J Iz I y A
( y2 z2 ) dA .
Exact Solutions for Cantilevered and Simplysupported Shafts under Harmonic Torsion

GENERAL EXACT CLOSED FORM SOLUTION The general exact closed form solution x x of equation
(12) consists of two parts, homogeneous solution xh x and
particular solutionxp x .
Shafts having cantilevered and simplysupported boundary conditions subjected to harmonic twisting moments are shown in Fig (2). The boundary conditions of the cantilevered shaft at both ends are given as:
x (x) x0 0 and GJx xxL Mx (x)xL
while for simplysupported shaft are given by:
x x0 x L0 0 .
Z displacements N , the vector of integration constants
mx(x,t) Mx(L,t)
A
21
is expressed in terms of the nodal torsional
O X displacements
N 1
x1
x2 1 by imposing the
conditions xh 0x1 andxh Lx2 (Fig. 3), yielding:
L

Cantilever shaft
xh 0
L
p>L
Z N 1
xh
1
mx(x,t)
E 0
(20)
O X E L
12
A1
G A1
12 22
L
where
N 1
x1
x2 1 is the vector of the nodal

Simplysupported shaft
Fig. (2): Shafts with cantilevered and simplysupported boundary conditions under harmonic torsions
torsional displacements.
From equation (20) and by substituting into equation (17), the following expression is obtained as:
Substituting the above boundary conditions into equation (17), the general exact closed form solutions governing the torsional steady state dynamic responses for cantilevered and simplysupported shafts under harmonic twisting moments are represented, respectively, by the following relations:
xh x
Hx N
12 1
12 1
Node 1
Z
x1
(21)
mx
1 1
1 J 2
xc x em1x em2x
m L m L
12 m1e 1
m2e 2
22 Mx L
mx
GJ
21
Node 2 Le
J 2
Twonoded Beam element
X
mx
1
(18)
Fig. (3): A proposed twonoded finite beam element
1
1 1
J 2
in which
Hx
E(x) 12G is the matrix of exact
x em1x
em2x
12
xs 12 em1L
em2L
mx
shape functions for the steady state dynamic torsional
22
response. It is noted that the exact shape functions obtained in
m
J 2
21
equation (21) exactly satisfy the exact homogeneous solution of the governing torsional equation. The exact shape
x
J 2
(19)


FINITE ELEMENT FORMULATION
In this section, a new twonoded finite beam element is developed for dynamic torsional analysis of structural members with circular crosssections under various harmonic twisting moments. Figure (3) shows the proposed twonoded finite beam element with two degrees of freedom per
functions developed in this study are depends on the length of the beam element, the properties of the crosssection and the exciting frequency of the torsional loading.
Matrix Formulation for Torsional Response
By using equations (17) with torsional displacement equation (15), and by substituting into the energy expressions (911), the following variations of energy expressions are given in terms of nodal torsional displacement as:
T 2 LJ Hx Hx dx (22)
L
L
element. A family of exact shape functions which exactly
0 N
21 12 N
1
satisfy the homogeneous solution of the dynamic torsional
U
GJ
H x H x
dx (23)
equation is employed to formulate the exact stiffness and
0 N
21 12 N 1
L
L
mass matrices and the load potential vector.
W
m x Hx
dxM x Hx L
(24)
Formulating Exact Torsional Shape Functions
To relate the exact homogeneous solution of the torsional
0 x N 12 1
x N 12 10
displacement function xh x to the nodal torsional
From equations (2224), by substituting into variational form of Hamiltons principle equation (8), leads to typical finite beam element model as:
solutions available in the literature. In the finite element Abaqus model, a twonoded B31 beam element having six degrees of freedom per element (i.e., three translations and three rotations) is used for comparison.
K
2 M
F
(25)
Example (1): Cantilever Hollow Shaft under Concentrated
e
e
N 1
e 1
Harmonic Twisting Moment
A 3000mm cantilever hollow shaft subjected to a
where the beam element stiffness matrix Ke is given by:
concentrated harmonic twisting moment
Mx 12.0eit kNm
K
GJH x H x dx
(26)
applied at the free end of the cantilever is shown in Fig. (4).
L
L
e 0 21 12
The mass stiffness matrix Me of the beam element is:
The shaft has an outer diameter of 100mm and inner diameter of 80mm and made of steel material with the following properties: the modulus of elasticity E 200GPa , shear
modulus G 70GPa and the density 7800kg / m3 . The
M
L2 JH x H x dx
(27)
purpose of this example is to assess the accuracy and validity
e
0 21 12
of the results obtained from the present finite element
The potential energy load vector Fe1 of the element is:
formulation. It is required to investigate the following:

The static torsional response of the shaft using very low
F
m xHx
dxM
xHx L
(28)
exciting frequency 0.011 , where the first natural
L
L
e 1
0 x
1 x
10
frequency of the given shaft is
1 261.7Hz ,
The above expressions for stiffness, mass and potential load vector formulated for twonoded finite beam element using the exact torsional shape functions developed to investigate the steady state dynamic torsional response of shafts under various harmonic torsional loadings.


NUMERICAL EXAMPLES AND DISCUSSIONS The present solutions (exact closed form and finite beam
element solutions) developed in this paper governing the torsional vibration of shafts under various harmonic twisting moments can be used to capture the following analyses:

The steady state dynamic response for the shafts with solid and hollow crosssections and under harmonic twisting moment at exciting frequency ,

The quasistatic response for the shafts under harmonic twisting moments by using very low exciting frequency compared to the first torsional natural frequency of the shaft,

Extracting the natural torsional frequencies and related mode shapes of the shaft under harmonic twisting moments.
To demonstrate the validity, accuracy and applicability of the exact closed form solution and finite beam element formulation developed, three examples are presented in this study. These examples investigate the torsional vibration of the shaft under various harmonic twisting moments and various boundary conditions. The present finite beam element developed is based on the exact torsional shapes functions which exactly satisfy the homogeneous form of the governing torsional vibration equation of the shaft. Due to this treatment, the mesh discretization errors induced in the classical finite element solutions using polynomial interpolation shape functions are eliminated. As a result, it is observed that, the results obtained based on a single finite beam element exactly matched with the corresponding results based on the exact closedform solutions developed in this study up to five significant digits. The numerical nodal results obtained from the present finite beam element are compared with the established finite beam element Abaqus and exact

The dynamic response of the shaft under harmonic torsion at exciting frequency 1.401 .
Z
Mx(3m,t)=12.0eit kNm
X
Y
L=3000mm
L=3000mm
80mm
100mm
Fig. (4): A cantilever hollow shaft under end harmonic twisting moment
In the Abaqus finite element model, the shaft is divided into 80 beam B31 element along the longitudinal axis of the shaft. In other words, the model has 486 degrees of freedom in order to achieve the required accuracy in this example. In constraint, the results obtained from the present finite element dveloped use only one twonoded beam element with two degree of freedom to attain the exact solution results.
Static Torsional Response
In order to approach the static response, the exciting frequency should be taken significantly lower than the first natural frequency of the cantilever hollow shaft. The static results for maximum torsional displacement xmax(x2 ) of the shaft under harmonic torsional loading at exciting frequency 0.011 are given in Table I. It is observed that the results for the nodal torsional displacement x2 obtained from the present finite element based on one twonoded beam element (i.e., 2 dof) provide excellent agreement with Abaqus finite element model based on 80 beam B31 elements (i.e., 486 dof). The present finite element formulation based on exact shape functions demonstrated that the new twonoded beam element provides excellent agreement with present exact closed form solution and Abaqus beam model by keeping the number of degrees of freedom a minimum.
Dynamic Torsional Response
The maximum steady state torsional displacement 2 (at the cantilever free end) of the cantilever shaft under the given harmonic twisting moment Mx 12.0eit kNm with exciting frequency 1.401 366.4Hz is provided in Table I. The maximum torsional displacement results based on the formulations developed in this study are compared with those based on Abaqus finite beam B31 element solution. The torsional displacement results obtained from the present exact solution and finite element formulation using a single beam element with 2 dof are found exactly identical to Abaqus beam model based on 80 beam elements with 486 dof in order to achieve the solution accuracy.
Extracting of Natural Torsional Frequencies
Under the given harmonic torsional loading, the natural torsional frequencies are extracted from the multiple steady state torsional dynamic analyses in which the exciting frequency of the torsional loading varying from nearly zero to 1600Hz. Fig. 6 shows the results for the nodal (maximum) torsional displacement 2 at the cantilever end against the forcing frequency . The torsional natural frequencies are determined at the peaks of the torsional displacementfrequency diagram as observed in Fig. 6. Peaks on the diagram indicate the resonance and are then identify the torsional natural frequencies of the given cantilever shaft.
Maximum torsional function x max (rad)
Maximum torsional function x max (rad)
0.200
Table I: Static and dynamic results for cantilever hollow shaft under end harmonic torsion
0.160
0.120
1=195.6Hz
2=586.9Hz
Present FE
Abaqus FE
3=978.2Hz
TYPE OF RESPONSE
Variable
Present FE
(2 dof)
Abaqus FE
(486 dof)
Exact Solution
Static
0.011
xmax (103 rad)
80.74
80.75
80.74
Dynamic
1.401
xmax (103rad)
50.53
50.53
50.53
TYPE OF RESPONSE
Variable
Present FE
(2 dof)
Abaqus FE
(486 dof)
Exact Solution
Static
0.011
xmax (103 rad)
80.74
80.75
80.74
Dynamic
1.401
xmax (103rad)
50.53
50.53
50.53
0.080
0.040
4=1369.5Hz
Example (2): Cantilever Shaft under Distributed Torsion
0.000
0 200 400 600 800 1000 1200 1400 1600
Exciting frequency (Hz)
A 4000mm cantilever shaft having solid crosssection of radius 40mm subjected to uniform distributed harmonic twisting moment m(x,t) 4000 eit Nm is considered as
shown in Fig. 5. The material of the cantilever is made from steel with the following mechanical properties; modulus of elasticity E 200GPa , modulus of rigidity G 70GPa ,
and the material density 7850kg / m3 . The objective of this example is to:

Extract the natural torsional frequencies from the dynamic response of the cantilever shaft under the given harmonic torsion, and

Determine the steady state torsional mode shapes of the cantilever shaft corresponded to the natural torsional frequencies.
The cantilever shaft in this example is modelled in Abaqus finite element using 100 beam B31 elements with six degrees of freedom per node along the longitudinal axis of the shaft, i.e., a total of 606 degrees of freedom is used in order to yield the accuracy of this problem. While the present finite element uses a single twonoded beam element with 2 dof to approach the corresponding results obtained from exact closed form solution.
Z
Fig. 6: Natural torsional frequencies of the cantilever under harmonic torsion
The first four torsional natural frequencies extracted from the steady state dynamic responses at the peaks are given in Table II. The values of the torsional natural frequencies obtained from the finite beam element developed in this study using a single twonoded beam element (2 dof) are presented and compared with the corresponding results obtained from the exact closed form solution developed and Abaqus solution using 100 beam B31 element (606 dof). It is noted from the results that, the present finite beam element solution with a minimum degrees of freedom exhibit excellent agreement when compared with other solutions, the exact closedform solution and Abaqus finite beam B31 element model having a large number of degrees of freedom. Thus, the present finite beam element model is able to capture the eigenfrequencies of the given cantilever shaft.
Table II: Natural torsional frequencies (in Hz) for the cantilever shaft under harmonic twisting moment
Number
Exact Solution
Abaqus FE (606 dof)
Present FE (2 dof)
1
195.6
195.6
195.6
2
586.9
586.7
586.9
3
978.2
978.0
978.2
4
1369.6
1369.3
1369.5
Number
Exact Solution
Abaqus FE (606 dof)
Present FE (2 dof)
1
195.6
195.6
195.6
2
586.9
586.7
586.9
3
978.2
978.0
978.2
4
1369.6
1369.3
1369.5
mx=4.0eit kNm/m
O X
L=4000mm
80mm
Fig. 5: A cantilever shaft under harmonic distributed twisting moment
Steady State Torsional Mode Shapes
The first five steady state torsional mode shapes for the dynamic vibration response of the cantilever shaft under the given distributed harmonic torsion
mx (x,t) 4000eit Nm
is shown in Fig. 6. The normalized
steady state torsional modes (2 xmax ) obtained using the
present finite beam element are plotted for the first four torsional exciting frequencies:
1 195.6Hz , 2 586.9Hz , 3 978.2Hz and
4 1369.5Hz , respectively.
Normalized Torsional function (2/x max)
Normalized Torsional function (2/x max)
1.20
0.80
Fig. 7: A clampedclamped shaft under various harmonic twisting moments
In order to establish the validity and accuracy of the present finite element based ontwonoded beam element, the nodal degrees of freedom results for static torsional response and steady state torsional dynamic response are obtained and compared against the results based on established Abaqus finite beam element. Under the present finite element solution, only four twonoded beam elements with a total of 5 degrees of freedom are used while in Abaqus finite element
0.40
0.00
0.40
0.80
1.20
1=195.6Hz
4=1369.5Hz
3
3
2=586.9Hz =978.2Hz
model, the shaft is consisted of twohundred twonoded beam B31 elements with a total of 1206 degrees of freedom along the shaft longitudinal axis to achieve the convergence.
Static Torsional Response
The static results for the nodal torsional displacement
xi (for i 1,2,3,4,5) are plotted against the shaft longitudinal
0.0 0.8 1.6 2.4 3.2 4.0
Shaft axis x (m)
Fig. 7: Normalized torsional mode shapes for cantilever shaft under harmonic distributed torsion
Example (3): ClampedClamped Shaft Verification of The Present Finite Beam Element
This example is presented to exhibit the ability of the present finite element developed in the paper to achieve the required accuracy for the static and dynamic responses for the given problem by comparing the present results with those based on the Abaqus finite beam model solution. A clamped clamped shaft under various concentrated and distributed harmonic twisting moments is considered as shown in Fig. 8. The shaft has a circular solid crosssection of radius 50mm and 8000mm span, while the mechanical material properties are given as: E 210 GPa , G 72 GPa and 7850kg / m3 .
It is required to assess the efficiency and accuracy of the present finite element formulation in evaluating the nodal torsional degrees of freedom (xi for i 1, 2,3, 4,5) for the following analyses:

The static torsional analysis of the clampedclamped shaft under the given harmonic twisting moments at very low exciting frequency in order to capture the static response.

The torsional dynamic analysis of the clampedclamped shaft under the given harmonic twisting moments at exciting frequency 50Hz .
axis x as illustrated in Fig. 8a, in which the static torsional response is approached by using a very small exciting frequency . The static results shows that, the nodal
torsional degrees of freedom obtained from the present finite beam element formulation having 5 degrees of freedom coincide on the corresponding results obtained from Abaqus beam model having 1206 dof and then provide an excellent agreement.
Dynamic Torsional Response
Fig. 7b shows the nodal degrees of freedom results for the steady state torsional response plotted against the shaft coordinate axis x . It is observed that, the developed finite beam element results based on four twonoded beam elements with 5 degrees of freedom shows again an excellent agreement with those results based on Abaqus finite element model using 200 beam B31 elements with 1206 degrees of freedom. This is a natural outcome of the fact that the present finite beam element formulation is based on the exact shape functions which exactly satisfy the homogeneous solution of the governing torsional vibration equation. This treatment eliminates the discretization errors occurred in the conventional finite element formulations which based on approximate interpolation shape functions.
Nodal Torsional displacement i (rad)
Nodal Torsional displacement i (rad)
0.020
0.016
0.012

Comparison with established Abaqus finite beam element and exact solutions available in the literature demonstrates the validity and accuracy of the present exact closed form solution and finite element formulation.


0.008
0.004
0.000
(a)
Nodal Torsional displacement i (rad)
Nodal Torsional displacement i (rad)
0.000
0.002
0.004
0.006
0.008
0.010
(b)
Abaqus Solution Present FE Solution
0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0
Shaft axis x (m)
Abaqus Solution Present FE Solution
Abaqus Solution Present FE Solution
0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0
Shaft axis x (m)
REFERENCES

C. K. Rao, Torsional frequencies and mode shapes of generally constrained shafts and piping, Journal of Sound and Vibration, 125(1), (1988), 115121.

E. J. Sapountzakis and V. G. Mokos, Warping shear stresses in non uniform torsion by BEM, Computational Mechanics, 30(2), (2003), 131142.

D.V. Mallick and R. Dungar, Dynamic characteristics of core wall structures subjected to torsion and bending. The Structural Engineer, 55, (1977), 251261.

J. W. Zu and R. P. Han, Natural Frequencies and Normal Modes of a Spinning Timoshenko Beam With General Boundary Conditions, Journal of Applied Mechanics 59(2) , (1992), 197204.

H. F. Bauer, Vibration of a Rotating Uniform Beam, Part 1: Orientation in the Axis of Rotation, Journal of Sound and Vibration 72, (1981), 177189.

H. P. Lee, dynamic response of a rotating Timoshenko shaft subject to axial force and moving loads, Journal of Sound and Vibration, 1995, 181 (1), 169177.

K. Koser and F. Pasin, Torsional vibrations of the drive shafts of mechanisms, Journal of Sound and Vibration, 199, (1997), 559565.

J. Melanson and J.W. Zu, Free vibration and stability analysis of internally damped rotating shafts with general boundary conditions, Journal of Vibration and Acoustics, 120(3), (1998), 776783.

S. C. Lin and K. M. Hsiao, Vibration Analysis of A Rotating Timoshenko Beam, Journal of Sound and Vibration 240(2), (2001), 303322.

S. A. Hosseini and S.E. Khadem, Free vibrations analysis of a rotating
Fig. 8: Nodal torsional displacements for Static and dynamic responses of clampedclamped shaft under various harmonic torsions


SUMMARY AND CONCLUSION
From the numerical results conducted throughout this study, the following concluding remarks are made:

The dynamic equation of motion for torsional vibration and related boundary conditions for shafts under various harmonic twisting moments are derived via Hamiltons principle.

Exact closedform solutions of steady state torsional response of shafts are obtained for cantilever and simply supported shafts.

The exact closedform solution derived is successfully used to formulate a family of exact shape functions which based on the homogeneous solution of the governing torsional equation.

The exact shape functions are used to formulate a super convergent finite beam element for the shafts. The proposed beam element has a two nodes and two degrees of freedom.

The present exact closed form solution and finite element formulation developed in this study are able to efficiently capture the quasistatic and steady state response of beams under harmonic torsional loading. It is also capable of extracting the eigenfrequencies and eigen modes.

The new beam element involves no discretization errors and generally provides excellent results compared with Abaqus finite element solution while keeping the number of degrees of freedom a minimum.
shaft with nonlinearities in curvature and inertia, Mechanism and Machine Theory, 44, (2009), 272288.

T. Yokoyama, Free vibration characteristics of rotating Timoshenko beams, International Journal of Mechanical Sciences, 30(10), (1988), 743755.

Y. A. Khulief, Vibration frequencies of a rotating tapered beam with end mass, Journal of Sound and Vibration 134, (1989), 8797.

L. W. Chen and D. M. Ku, Finite element analysis of natural whirl speeds of rotating shafts, Computer and Structures, 40, (1991), 741 747.

M. Y. Cang, J. K. Chen, and C. Y. Chang, A simple spinning laminated composite shaft model, International Journal of Solids and Structures, 41(34), (2004), 637662.