| THEORY & FORMULAE |
Solid cylinderical shafts involved in transmission of power by means of gears, flywheel, pulleys are subjected not only to torsion but also bending. All rotating shafts even in the absence of external load, deflect during rotation. The combined weight of the shaft and the load can cause deflection that will create resonant vibration at certain speed known as the Critical Speed.
One of the key methods for calculating critical speed is the Dunkerley Equation. In 1894, Stanley Dunkerley published a study "On the Whirling and Vibration of Shafts," Phil. Trans. R. Soc., London. The first sentence of his paper reads, "It is well known that every shaft, however nearly balanced, when driven at a particular speed, bends, and, unless the amount of deflection be limited, might even break, although at higher speeds the shaft again runs true. This particular speed or 'critical speed' depends on the manner in which the shaft is supported, its size and modulus of elasticity, and the sizes, weights, and positions of any pulleys it carries". Dunkerley equation is an approximation to the first natural frequency of vibration of the system, which is assumed to be nearly equal to the critical speed of rotation.
The method consists of reducing the actual system of multi-rotors to a series of single rotor system. The critical speed of a single rotor system can be calculated by a direct formula (i.e. natural frequency of a pinned-pinned beam with intermediate load). Finally these individual critical speeds are combined to obtain the critical speed of the system as given in the equations below:
    
where
     L = length of shaft
     ai = location of rotor i from left end of shaft
     mi = mass of rotor i
     n = total number of rotors, including shaft; i = 1, ...n
     ωc = critical speed (cycle/sec[Hz], revolution/minute)
     ωi = critical speed for rotor i
     δi = static deflection at at position of rotor i, due to rotor i only
     g = acceleration due to gravity: 9.81 m/s2, 32.2 ft/s2
     E = modulus of elasticity of shaft material
     I = area moment of inertia of cross-sectional area about axis through centroid
The shaft is also considered a rotor with mass ms concentrated at its mid-point.
◊ Use link
EXAMPLE Of Input/Output
to demo data entry expectations and results; you may edit & use it as starting point