Dynamics: Torsional Vibration

Torsional vibration is the oscillatory twisting of shafts, gears, and rotors in mechanical systems, a phenomenon distinct from lateral bending or axial vibration. While often invisible, its effects can be catastrophic, leading to sudden fatigue failures in power transmission systems like engine crankshafts, turbine-generator sets, and propeller driveshafts. Understanding its principles is essential for designing reliable rotating machinery, predicting resonant conditions, and implementing effective damping or detuning strategies to prevent costly downtime and dangerous fractures.

Fundamentals of Stiffness and Inertia

At the heart of any vibrational analysis are two core properties: stiffness and inertia. In torsional systems, torsional stiffness ($k_t$) quantifies a shaft's resistance to twist. It is defined as the torque required to produce a unit angular deflection. For a uniform, solid, circular shaft of length $L$, shear modulus $G$, and polar second moment of area $J$, the stiffness is given by $k_t=\frac{GJ}{L}$. A higher stiffness means the shaft twists less for a given applied torque.

The counterpart to stiffness is mass moment of inertia ($I$), which quantifies a rotating mass's resistance to angular acceleration. It depends on the mass's geometry and its distribution relative to the axis of rotation. For a simple disc, $I=\frac{1}{2}m{r}^2$ for a solid cylinder. Inertia stores kinetic energy, while stiffness stores potential strain energy; their interaction governs the system's vibratory behavior. You must accurately calculate or obtain these values for each rotating component in your system as the first step in any analysis.

Equation of Motion and Natural Frequency

The simplest torsional system is a single disc with mass moment of inertia $I$ attached to a massless shaft of torsional stiffness $k_t$, fixed at the far end. When disturbed, it exhibits free vibration. Applying Newton's second law for rotation, the sum of the inertial torque ($I\ddot{\theta }$) and the restoring spring torque ($k_t\theta$) equals zero. This yields the fundamental equation of motion for torsional vibration:

$$I\ddot{\theta }+k_t\theta =0$$

This is a standard second-order linear differential equation. Its solution describes simple harmonic motion: $\theta (t)=\Theta \cos ({\omega }_{n}t+\varphi )$, where $\Theta$ is the amplitude and $\varphi$ the phase angle. The system's natural frequency of torsional systems (${\omega }_n$) in radians per second is derived from the equation's characteristic equation:

$${\omega }_n=\sqrt{\frac{k_t}{I}}$$

In Hertz (cycles per second), the natural frequency is $f_n=\frac{1}{2\pi }\sqrt{\frac{k_t}{I}}$. This frequency depends solely on the system's physical properties ($k_t$ and $I$), not on the initial disturbance. A stiffer shaft or a smaller inertia will result in a higher natural frequency.

Analyzing Multi-Disc Torsional Systems

Real machinery often involves multiple rotating inertias connected by shafts, forming a multi-disc system. A two-disc system connected by a shaft is a fundamental building block. Here, both ends are free to rotate, and the shaft twists between them. The system has one natural frequency (for a two-disc, one-shaft system) corresponding to a mode shape where the discs oscillate in opposition. The natural frequency is calculated using the concept of an equivalent system. One common method uses the formula:

$${\omega }_n=\sqrt{k_t\left(\frac{1}{I_1}+\frac{1}{I_2}\right)}$$

where $I_1$ and $I_2$ are the inertias of the two discs. For systems with three or more discs, you must set up and solve the equations of motion for each inertia, leading to a multi-degree-of-freedom eigenvalue problem. This yields multiple natural frequencies, each with a distinct mode shape showing the relative twist amplitudes of each disc. Solving these typically requires matrix methods: $[-{\omega }^2[I]+[K]]\{\Theta \}=\{0\}$, where $[I]$ is the inertia matrix and $[K]$ is the stiffness matrix.

Modeling Geared Systems and Equivalent Inertia

Geared systems introduce a kinematic complication because different shafts rotate at different speeds. To analyze the torsional vibration of the entire system on a single reference shaft (usually the prime mover shaft), you must calculate equivalent inertia. The key principle is the conservation of kinetic energy. A mass moment of inertia $I_2$ on a shaft rotating at a speed $N_2$, connected via gears with a speed ratio $n=N_1/N_2$ to a reference shaft (speed $N_1$), has an equivalent inertia on the reference shaft of $I_{2,eq}=I_2\cdot n^2$.

For example, if a gear on the reference shaft (1) has 20 teeth driving a gear with 40 teeth on shaft 2, the speed ratio $n=N_1/N_2=40/20=2$. An inertia $I_2=5{\text{ kg-m}}^2$ on shaft 2 would appear as $I_{2,eq}=5\times (2{)}^2=20{\text{ kg-m}}^2$ on shaft 1. Shaft stiffnesses must also be referred to the reference shaft using the square of the speed ratio. Once all inertias and stiffnesses are reduced to a common shaft, you can model the system as a linear chain of discs and springs and calculate its natural frequencies using multi-disc methods.

Applications to Engine Crankshaft Vibration Analysis

The analysis of engine crankshaft vibration is a critical practical application. A crankshaft is a complex, multi-inertia system subjected to periodic torque pulses from cylinder firing. These pulses contain energy at the engine firing frequency and its multiples (harmonics). If any harmonic excites a torsional natural frequency of the crankshaft system, severe resonant vibrations occur.

The analysis involves modeling the crankshaft, flywheel, damper, and driven load (e.g., generator, propeller) as a series of concentrated inertias connected by stiffnesses representing crankshaft throws and shaft sections. Dampers (viscous or tuned) are often included in the model. The goals are to: 1) Calculate the system's natural frequencies and mode shapes, 2) Map the exciting torque harmonics across the engine's operating speed range (a Campbell diagram), and 3) Identify critical speeds where resonance occurs. The designer then aims to shift natural frequencies away from major excitation lines or incorporate a torsional vibration damper to attenuate resonant responses, ensuring the crankshaft operates safely below its endurance limit.

Common Pitfalls

1. Ignoring Gyroscopic Effects in Geared Systems: When reducing geared systems, a common mistake is to forget that the direction of rotation matters if the axis of vibration is also precessing. For pure torsional analysis about a fixed axis, the equivalent inertia method using $n^2$ is correct. However, if lateral or gyroscopic coupling is significant, a more complex model is needed.
2. Incorrectly Calculating Equivalent Stiffness: For shafts in series (e.g., a crankshaft with multiple throws), the equivalent torsional stiffness is not simply the sum. Torsional stiffnesses in series combine like springs in series: $\frac{1}{k_{eq}}=\frac{1}{k_1}+\frac{1}{k_2}+...$. Using the sum will overestimate stiffness and thus over-predict the natural frequency.
3. Confusing Mass Moment of Inertia with Polar Second Moment of Area: A critical error is using $J$ (polar second moment of area, with units m⁴) in place of $I$ (mass moment of inertia, with units kg-m²) in the natural frequency formula. Remember: $k_t$ uses $J$ (a geometric property), while the equation of motion uses $I$ (a mass distribution property).
4. Neglecting Damping in Resonance Assessments: While the undamped natural frequency is vital for initial design, completely ignoring damping when predicting resonant amplitude can be misleading. Damping determines how severe the resonance peak will be. In systems like crankshafts with dampers, evaluating the system's response requires including damping in a forced vibration analysis.

Summary

• Torsional vibration analysis centers on the interplay between torsional stiffness ($k_t$), which resists twist, and mass moment of inertia ($I$), which resists angular acceleration.
• The fundamental equation of motion $I\ddot{\theta }+k_t\theta =0$ leads to a natural frequency of ${\omega }_n=\sqrt{k_t/I}$ for a single inertia system.
• Multi-disc systems require setting up equations for each inertia, leading to multiple natural frequencies and mode shapes solvable via matrix eigenvalue problems.
• In geared systems, all parameters must be reduced to a common reference shaft using kinetic energy equivalence, where equivalent inertia is found using the square of the speed ratio ($I_{eq}=I\cdot n^2$).
• Engine crankshaft vibration analysis applies these principles to model the shaft assembly, predict natural frequencies, and avoid resonant conditions induced by periodic firing pulses, often using Campbell diagrams and necessitating dampers.