A shaft is modelled using three-dimensional solid finite elements. The shear-deformation and rotary inertia effects are automatically included through the three-dimensional elasticity formulation. The formulation allows warping of plane cross-sections and takes care of gyroscopic effect. Unlike a beam element model, the present model allows the actual rotor geometry to be modelled. Shafts with complicated geometry can be modelled provided that the shaft cross-section has two axes of symmetry with equal or unequal second moment of areas. The acceleration of a point on the shaft is determined in inertial and rotating frames. It is found that the finite element formulation becomes much simpler in a rotating frame of reference that rotates about the centre-line of the bearings with an angular velocity equal to the shafts spin speed. The finite element formulation in the above frame is ideally suited to non-circular shafts with solid or hollow, prismatic or tapered sections and continuous or abrupt change in cross-sections. The shaft and the disc can be modelled using the same types of element and this makes it possible to take into account the flexibility of the disc. The formulation also allows edge cracks to be modelled. A two-dimensional model of shaft disc systems executing synchronous whirl on isotropic bearings is presented. The application of the two-dimensional formulation is limited but it reduces the number of degrees of freedom. The three-dimensional solid and two-dimensional plane stress finite element models are extensively validated using standard available results.
LARGE EDDY SIMULATION OF TURBULENT INCOMPRESSIBLE FLUID FLOWS BY A NINE-NODES CONTROL VOLUME-FINITE ELEMENT METHOD