Torsional Natural Frequencies PDF
Torsional Natural Frequencies PDF
by
Qingyu Wang
Mechanical Engineer
Elliott Group
Jeannette, Pennsylvania
Troy D. Feese
Senior Project Engineer
Engineering Dynamics Inc.
San Antonio, Texas
and
Brian C. Pettinato
Manager, Product Development
Elliott Group
Jeannette, Pennsylvania
Normalized Amplitude
0
Figure 2. Lumped Inertia-Stiffness (Mass-Elastic) System.
A typical representation of SM is the interference or Campbell Feese and Hill, (2009) show a measured AF, and the AF
diagram as shown in Figure 5. for common equipment trains is between 10 and 50 (damping
ratio 0.01 and 0.05). The system gain can be calculated using
Equation 5 for different speeds (different SMs), and different
Compressor Speed
2X Line Frequency
Line Frequency
damping ratios/AFs, as shown in Figure 6.
Motor Speed
25000
50
22962 =0.01
40
20000 2*RS
30
=0.0167
|G (i )|
15000
Frequency (cpm)
20
=0.025
=0.05
10
10000 1*RS
0
0.8 0.85 0.9 0.95 1 1.05 1.1 1.15 1.2
n
5000
865 The numerical results in Table 1 show that for the API
0
0 2000 4000 6000 8000
recommended separation margin of 10%, amplification would
Speed (rpm) be approximately 5 for most cases, and amplitudes would
almost double if the separation margin was 5% instead of the
Figure 5. Interference Diagram. recommended 10%. Low separation of 1% or less will
increase the amplitudes significantly.
Calculation of the AF is not a required part of the API 617
torsional analysis nor is it related to the SM requirement (10%). Table 1. Magnification at Different SMs and AFs.
However, SM and AF are related to vibration amplitude, which
can be demonstrated using a simple model. Damping Separation Margin (SM)
Ratio 0% 1% 5% 10% 15%
A single-degree-of-freedom mass-damper-spring system is
depicted having equation of motion as shown in Equation 4, 0.05 AF = 10 10 7 5 3
and the gain/vibration amplitude as shown in Equation 5 (see 0.025 20 19 9 5 4
Meirovitch, 1986 for details). 0.0167 30 26 10 5 4
0.01 50 36 10 5 4
x(t ) 2 n x (t ) n2 x(t ) n2 Ae it (4)
Uncertainties in Modeling
where n is the natural frequency, and is the damping ratio.
If the lumped inertia-stiffness model truly represents the
x(t ) A G(i ) ei (t ) ,where equipment train, then the calculated TNFs will match the real
systems TNFs. All factors that cause a difference between the
1
G(i ) calculated TNFs and the true TNFs of the system are
1 / 2 /
(5) considered uncertainties in this paper. There are two types of
2 2 2
n n uncertainties in the modeling: errors and variations.
Influence of Uncertainties S.O.1-7 (M-G-C) and S.O.9 (T-C-M, i.e. no gear) are used
for the second study, estimating the worst case scenario.
With the estimation of the uncertainties for each Assume the uncertainties of all inertias are within 5%, and
component, it is possible to evaluate uncertainties in the uncertainties of all stiffness are within 5% except for the
predicted/calculated TNFs. coupling stiffness. Two sets of uncertainties of the coupling
stiffness are used, 15% and 25% (see Table 2 for reference),
Motor-gear-compressor trains are a common torsionally and the results are shown in Table 5.
concerned category. This paper conducted a study on TNFs
using seven shop orders of motor-gear-compressor trains whose The first five columns in Table 5 are the same as in
TNFs have all been measured (all with partial or no load). Table 4. The last two columns are the differences (in absolute
Some parameters are listed in Table 3. Two studies were values) between column 3 and the calculated TNFs with the
performed with these shop orders. assumed uncertainties, i.e., they are the uncertainties of the
predictions. The prediction uncertainties can be considered as
First Study - Sensitivity the overall sensitivity.
A sensitivity study was performed by varying each Since Table 5 shows the worst case scenario, it means if
parameter sequentially, such as the inertia of the motor or the the coupling stiffness has accuracy within 15%, and all other
stiffness of a coupling, to see the change in the calculated stiffness and inertias are within 5% of the real values, then
TNFs. The purpose is to find whether discrepancies between predicted TNFs should be within 10% of the true values. The
prediction and measurement are related to parameter variation. prediction uncertainty approaches 15%, when the coupling
stiffness accuracy is within 25%.
The simulation is done by applying 1% variation to the
inertia of the motor, gear, compressor, and the stiffness of the To investigate the discrepancies, some of the torsional
two couplings. The results are shown in Table 4, where I1-4 analyses were thoroughly checked to make sure there were no
stands for the inertia of the motor/gear/compressor/turbine, and human errors present. For some of the shop orders, the motor
K1-3 represents the torsional stiffness of the couplings. The and gear vendors were requested to confirm/verify their
predicted and measured values are listed for the coupling numbers. One coupling vendor even performed a static test to
modes of each shop order. In some cases, the 2 nd TNF could verify their predicted coupling stiffness. All checks indicated
not be measured. no major errors/mistakes, yet the differences are more than
10% from the measurements for some shop orders. One
The sensitivity study results, such as -0.35% and 0.00% in explanation would be the variation-type uncertainties for the
the table for S.O.1 in column I1, are calculated as: use +1% couplings, especially since all tests were performed with only
inertia for the motor and keep all other inertias and stiffness partial or no load. Another explanation would be that there
unchanged, calculate the TNFs (new TNFs), and the difference were still human errors, since some vendors may not be fully
between the new TNFs and the predicted for 1st and 2nd TNFs aware of the torsional modeling techniques.
Shop
Start Year Driver Rated Power (HP) Driver speed (RPM) Gear Ratio Train Type1
Order
S.O.1 2010 11000 1000-1550 8.03 M-G-C
S.O.2 2009 4000 1200-1600 6.55 M-G-C
S.O.3 2009 4000 1200-1600 8.19 M-G-C
S.O.4 2009 9300 1450-1900 5.32 M-G-C
S.O.5 2009 13000 1400-1900 4.93 M-G-C
S.O.6 2009 12000 1400-1900 2.86 M-G-C
S.O.7 2009 1770 1400-1880 3.35 M-G-C
S.O.8 2005 1350 1040-1559 7.5588 M-G-C
S.O.9 2003 16200 3420-3636 1 T-C-M
S.O.10 2003 16200 3420-3636 1 T-C-C-M
S.O.11 1971 39793 3465-4851 1 T-C-C-C
S.O.12 1979 3485 5266-7899 1.4544 T-G-C
S.O.13 1977 6000 1800 4.2644 M-G-C
1 Trains are listed from driver to driven side. M means motor if on driver side (such as M-G-C), and means motor/generator if on
driven side (such as T-C-M/G). G means gear, C means compressor, and T means turbine.
2 This train has a variable speed drive, and its internal low stiffness rather than the couplings dominates the first critical.
70000
60000
50000
Static Torque (lb-in)
40000
30000
20000
10000
0
-0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1
-10000
Angle (deg)
The compressor system consisted of the following: As shown, the system is most sensitive to motor and
flywheel inertia as well as motor shaft extension stiffness.
There is some sensitivity to the coupling stiffness and no
Induction motor, rated 400 HP (298 kW) at 595 RPM
sensitivity to motor core stiffness, and mass-elastic properties
Disc pack coupling - hub on motor side and adapter on
of the compressor within a 30% to 300% range of the provided
compressor side values.
Flywheel mounted to compressor crankshaft
Reciprocating compressor with two throws To explain this, the percent strain energy and kinetic
energy for the first torsional mode were calculated. Results
A torsional analysis was performed in the design stage. show that 65% of the strain energy is in the motor shaft
The coupling size was selected based on proper service factor. extension and 30% is in the coupling. Therefore, these are the
The flywheel was sized to tune the first TNF between most sensitive components in terms of stiffness. Results show
compressor orders. For this case, the first TNF was predicted that the motor core and flywheel each have approximately 47%
to be 55 Hz, which would be between 5 and 6 compressor kinetic energy for the first torsional mode. This makes both of
speed. The second TNF was predicted to be above 30 these components the most sensitive from an inertia standpoint.
running speed and therefore not considered a problem.
The flywheel inertia calculation is straight forward and
Because of previous concerns with coupling stiffness should therefore be accurate. The compressor inertia and
values provided by the manufacturer and the criticality of stiffness were shown not to affect the first TNF. So for this
tuning the first TNF between harmonics to avoid high torsional system the uncertainty is most likely with the motor inertia and
vibration, measurements were recommended. A strain gage stiffness. Unfortunately, since the system was installed, there
telemetry system was mounted to the motor shaft extension. was no practical way to verify motor dimensions, weights, etc.
The first TNF was determined to be 59 Hz instead of 55 Hz as
predicted. The coupling stiffness could also affect the results, but
would need to have a larger percentage of error (80% for the
The percent error is approximately 7%, which is not that coupling versus only 25% for the motor). It was unknown if
far from the expected range of 5%. However, the first TNF the coupling manufacturer had actual factory test data to
is now only 1% from the 6 harmonic and electrical frequency substantiate the provided value.
of 60 Hz. This is much less than the API recommended
separation margin of 10%. For this particular system, it would This particular system had a heavy shrink fit on the
not be possible to achieve a 10% separation margin from both coupling hub, small diameter motor shaft, and no keyway. It
5 and 6 compressor harmonics even if the actual TNF had is interesting to note that if the 1/3 rule is omitted, then the
been 55 Hz as predicted. calculated first TNF would match the measured value.
Ignoring the penetration effect would be equivalent to modeling
The compressor system was tested over all load conditions a stepped shaft.
and found to still have acceptable stress levels in the motor
shaft and dynamic torque in the coupling. However, the
safety factor was reduced to less than two, which is normally
Multiple disc pack coupling failures were caused by high Based on this simple assessment, it would seem that when
torsional vibration due to insufficient separation margin from a the coupling spacer is torsionally much softer than the disc
significant compressor harmonic. All of the compressors were pack, uncertainty of the overall coupling stiffness will be
driven by 8000 HP (5966 kW) synchronous motors at 720 relatively low (example shown in Figure 12). However, when
RPM. The same flexible disc couplings were used on all units the spacer is torsionally stiff relative to the disc pack and no
between the motor shaft and compressor flywheel. hub (flanged directly to flywheel therefore no 1/3 rule),
uncertainty could be much higher than normal as shown in
Units that experience failures had a Siemens motor driving Figure 16. In the second example, the disc packs control the
an Ariel KBV/6 compressor. The first TNF was predicted to overall torsional stiffness of the coupling and are the most
be 3203 CPM (53.4 Hz). Measurements showed the actual difficult to predict. Due to the high torsional stiffness of the
frequency to be 58.5 60.5 Hz on the three units tested. larger coupling, it may be difficult for the manufacturer to test.
Analysis error was approximately 9% to 12% and the
separation margin was non-existent.
The target for the first TNF was 54 Hz, which is half-way With the engine damper functioning properly, operating
between 4 (48 Hz) and 5 (60 Hz) compressor harmonics. throughout the entire speed range does not pose a problem for
The second TNF was measured at 80 Hz, which is 6.7 the system; therefore, it is very important to check the viscous
compressor speed. The separation margins of 10% for the dampers on a periodic basis. When the dampers are located
first TNF and 5% for the second TNF are considered inside the engine frames, they can be subjected to elevated
acceptable. By properly tuning the TNFs, the units have been temperatures and sometimes overlooked during routine
running satisfactorily. maintenance (Feese and Hill, 2009).
-25.00%
Table 11 shows a variation of 33.8 Hz 35.6 Hz for the TNFs
first TNF. The average TNF was 34.9 Hz with a standard
deviation of 1.8% and having total variation of 2% to 3%. Figure 17. Discrepancies between Prediction and Measurement
Table 11. Comparison of TNFs Measured The distribution of the discrepancies is summarized in
Table 12, where for instance, 90% of all data points for the first
Unit Frequency (Hz) Speed (RPM)
TNF are within 15%.
Station H
1 34.5 460
2 34.5 460 Table 12. Percentage of Discrepancies
3 35.3 470
4 33.8 450 TNF within within within within within
Station D 25% 20% 15% 10% 5%
1 35.6 475
1st 100% 97.5% 90% 70% 30%
2 35.2 469
nd
3 35.3 471 2 100% 100% 93% 87% 53%
Because the fan normally operates 1000 to 1200 RPM, with the mechanical engineers that this change was being made
which was the speed range where excessive dynamic torque and did not realize how it could impact the mass-elastic
occurred, this is believed to be the reason for the coupling torsional system.
failure. For example, the VFD excitation was approximately
5% of the full load torque (FLT) and at torsional resonance the American Petroleum Institute (API) recommends that a
dynamic torque is amplified by a factor of 30. Therefore, the torsional analysis be performed in the design stage to prevent
maximum alternating torque was approximately 150% of the failures. A separation margin (SM) of at least 10% between
transmitted torque, which exceeded the rating of the coupling. the torsional natural frequencies and the excitation frequencies
is recommended to avoid running at a torsional resonance
Due to the large diameter and weight involved, the inertia unless shown to be safe. Many times satisfying the 10% SM
of the fan is many times greater than the inertia of the motor. is impractical for VFD motor systems that operate over a large
For the first torsional mode, the motor is typically near an speed range.
anti-node and acts like a torsional pendulum. The fan is
usually near the node and acts as an anchor. The VFD infers A torsional analysis of the fan system was never performed
load changes by monitoring motor current, which could also with either motor. After the coupling failure, the motor inertia
contain variation from the first TNF. In a torsionally stiff, values were compared. It was found that the replacement
lightly damped system, the first TNF is very sensitive to any motor had a much lower inertia (WR2) value than the original
harmonic excitation or sudden speed adjustments from the VFD motor. Reducing the motor inertia caused the first TNF of the
motor (Feese and Maxfield, 2008). system, which was originally below the minimum speed, to
increase into the normal operating speed range.
After further discussion with plant personnel, it was
determined that the fan was originally driven by another motor Since switching back to the original motor was not an
from a different manufacturer. Motor repairs were needed and option; a temporary solution was recommended where the
would take longer than acceptable. Therefore, an alternate running speed should be limited to a maximum motor speed of
motor was acquired and installed. The new motor from a 1000 RPM (VFD frequency of 50 Hz) to avoid exciting the first
different manufacturer was similar in electrical performance, TNF at 58 Hz. This provided a SM of 13% between the VFD
but was vastly different in physical size and inertia. excitation frequency and the first TNF of the system.
Unfortunately, the electrical engineers did not communicate
During the first event, the compressor unit tripped on high In general, induction motors, gearboxes, and centrifugal
vibration due to a failure of the low speed coupling. A picture compressors do not normally produce significant torsional
of the failed diaphragm coupling is shown in Figure 24. The excitation. For these types of rotating systems, the transmitted
coupling manufacturer concluded that the failure was due to torque should be smooth (less than 10% dynamic torque);
torsional fatigue (Corcoran, et al., 2010). Furthermore, the however, as noted in several reference papers, VFDs have been
manufacturer confirmed that the coupling was properly sized shown to create significant torsional excitation in systems and
for this application, speed, and load. Metallurgical study caused torsional vibration and failures of the system.
indicated fatigue cracks. All of these coupling failures
exhibited characteristics typical of high torsional vibration. Various tests were performed by sweeping the speed of the
VFD motor. The first TNF was measured at 13.5 Hz. This
compared to the predicted TNF of 12 Hz or 720 CPM;
therefore, the actual first TNF was 12.5% higher than predicted.
Shunt calibration is usually satisfactory for calibration of Figure A.4. Encoder on End of Compressor Crankshaft.
the strain gage system. For this method, it is assumed that the
strain gages have been properly affixed to the shaft surface and An alternate way to install an encoder is shown in Figure
orientated. If so, it is generally thought that the shunt A.5. In this case, a jack shaft from the front of the engine was
calibration will be within 3% of the actual value, which is more driving a cooling fan, which prevented installing the encoder
than adequate for most trouble-shooting field studies. directly to the engine crankshaft; therefore, a special mounting
bracket was designed with a spring to apply a constant force to
By applying a known unbalance in the bridge circuit, the the encoder wheel to prevent slippage relative to the jack shaft.
telemetry system can be calibrated with the shunt resistor. With The measurements were adjusted by the diameter ratio of the
some telemetry systems, the shunt resistor may be built-in. jack shaft and the rider wheel.
Laser Vibrometer
Electrical Measurements
API Standard 617, 1973, Centrifugal Compressors for General Fraunhofer ITWM, A Contactless Torque Sensor for Online
Refinery Services, Third Edition, American Petroleum Monitoring of Torsional Oscillations, 2007.
Institute, Washington, D.C.
Holset Flexible Couplings Catalog, 510/3.90/F.H., Application
API Standard 617, 2002, Axial and Centrifugal Compressor Information Type PM Industrial, Cincinnati, Ohio.
and Expander-Compressors for Petroleum, Chemical, and
Gas Refinery Services, Seventh Edition, American Hudson, J., Feese, T., Torsional Vibration A Segment of
Petroleum Institute, Washington, D.C. API 684, Proceedings of the 35th Turbomachinery
Symposium, Texas A&M University, College Station,
API Standard 671, Special Purpose Couplings for Texas, 2006.
Petroleum, Chemical and Gas Industry Services,
September 2010. Ker Wilson, W., Practical Solution of Torsional Vibration
Problems, Volume 1, New York, New York: John Wiley
API Standard 684, API Standard Paragraphs Rotordynamic and Sons, 3rd Ed., 1956.
Tutorial: Lateral Critical Speeds, Unbalance Response,
Stability, Train Torsionals, and Rotor Balancing, 2nd Kerkman, R.J., Theisen, J. and Shah, Kirti, 2008, PWM
Edition, August 2005. Inverters Producing Torsional Components in AC Motors,
Petroleum and Chemical Industry Technical Conference,
Bosin, D., Ehrich, R., and Stark, M., 1999, Torsional PCIC 2008. 55th IEEE.
Instabilities of Motor Driven Turbomachinery,
Turbomachinery International, pp 18-20. Kocur, J.A. and Muench, M.G., Impact of Electrical Noise on
the Torsional Response of VFD Compressor Trains,
Calistrat, M. M., and Leaseburge, G. G., 1972, Torsional Proceedings of the First Middle East Turbomachinery
Stiffness of Interference Fit Connections, ASME Paper Symposium, 2011.
72-PTG-37.
Meirovitch, L., 1986, Elements of Vibration Analysis, New
Corcoran, Kocur, and Mitsingas, Preventing Undetected Train York, New York: McGraw-Hill Inc.
Torsional Oscillations, Proceedings of the 39th
Turbomachinery Symposium, Houston, Texas 2010. Mondy, R. E. and Mirro, J., 1982, The Calculation and
Verification of Torsional Natural Frequencies for
Dawson, B. and Davies, M., 1975, An Improved Holzer Turbomachinery Equipment Strings, Proceedings of the
Procedure for Torsional Vibration Analysis, The Journal Twelfth Turbomachinery Symposium, Turbomachinery
of Mechanical Engineering Science, 17, 1, pp 26-30. Laboratory, Texas A&M University, College Station,
Texas, pp. 151-156.
De Choudhury, P., Torsional Analysis of Unit 59-K-400A
Comparative Study with Comments, Elliott Company Murray, B.D., Howes, B.C., Zacharias, V., Chui, J.,
Inter-Office Letter, April 3, 1979. Sensitivity of Torsional Analyses to Uncertainty in
System Mass-Elastic Properties, Presented at the
Ewins, D.J., 2000, Modal Testing: Theroy, Practice and International Pipeline Conference, Calgary, Alberta, June
Application, Second Edition, Research Studies Press Ltd. 1996.
Feese, T.D., and Hill, C.H., 2009, Preventing Torsional NEMA Standards Publication MG 1-2011, Motors and
Vibration Problems in Reciprocating Machinery, Generators, Rosslyn, Virginia, 2011.
Proceedings of the Thirty-Eighth Turbomachinery
Symposium, Turbomachinery Laboratory, Texas A&M Nestorides, E.J., A Handbook on Torsional Vibration,
University, College Station, Texas, pp. 213-238. Cambridge, 1958.