Rotating machinery imbalance
Rotating machinery imbalance
(OP)
Forgive my multiple posts but I have been directed to a couple different forums to better answer my question.
I have an embedded DFT that identifies the velocity magnitude and phase of the heavy spot on rotors, from time-series accelerometer and key-phasor data.
I wish to validate and eventually certify these results.
Can someone on the forum refer me to an individual or lab that offers such services?
I have an embedded DFT that identifies the velocity magnitude and phase of the heavy spot on rotors, from time-series accelerometer and key-phasor data.
I wish to validate and eventually certify these results.
Can someone on the forum refer me to an individual or lab that offers such services?





RE: Rotating machinery imbalance
RE: Rotating machinery imbalance
Are you looking for certification of the sensors and instrumentation? Or for the amount of residual unbalance of the rotor?
Vibration levels can be, and often are related to unbalance. But are not automatically equivalent to unbalance. One way to Establish the relationship involves a survey made by placing trial masses on the rotor at multiple locations and recording the resulting 1X vibration and phase
RE: Rotating machinery imbalance
I have in-flight time-series accelerometer and key-phasor data collected from aircraft.
My embedded DFT is giving me a report as to above parameters. My immediate need is to validate the report provided by my embedded DFT algorithm (determine degree of accuracy).
RE: Rotating machinery imbalance
If you are balancing the rotor by itself, you can draw a correlation using influence coefficient. You can also move a known unbalance (inch-ounce) around the rotor and plot the magnitude response... it will have a steady component and a component which varies sinusoidally with position. The sinusoidally varying position is due to your known added unbalance and can be used as conversion to convert the dc component of velocity (unknown unbalance) into a known unbalance.
As far as validating FFT... still don't know what you're after.
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Rotating machinery imbalance
Or you could add a known imbalance to a rotor for which you have already have established the influence coefficients.
Sorry this doesn't actually answer your question
Incidentally in one of your previous posts you mention that you want to find a lab near Detroit to check this for you. I'd look in the yellow pages for an instrument calibration service.
Cheers
Greg Locock
New here? Try reading these, they might help FAQ731-376: Eng-Tips.com Forum Policies http://eng-tips.com/market.cfm?
RE: Rotating machinery imbalance
Here is a good calibration lab that may help:
http://wccl.com/contacts.htm
I assume you want to verify that instrument is workly correctly first and then get a certified electronic calibration (see lab link). One practical way to test end-to-end functions is to compare to another known instrument with measurments on a rotor kit or low cost balance simulator. You could then verify amplitude, speed/frequency, and phase by testing at different speeds and with different weight amounts and at different weight locations. I use a Sears variable speed grinder with a custom 6" diameter aluminum disk with several tapped holes for set-screw weights. You could beg, borrow, rent, or steal a balance analyzer or pay someone like me to do the functional tests and then go for lab calibration and certificate. You cannot get an unbalance weight display without first "calibrating" the machine with a known unbalance weight.
Walt
RE: Rotating machinery imbalance
We have talked in the past.
Yes this is correct. I have a commercial balancer (which may or not need calibration) with which to compare the results from my embedded algorithm, and the reports are very close, within 30 degrees and 0.02ips, repeatable with real-time acquisition from my desktop rotor simulator, and 10 degrees and similar phase discrepancy from previously acquired flight data. But I'd like to determine which is closer in order to characterize my algorithm.
Initially, I'd like to find someone whom I can email my acquisition data to, who is set up to analyze it.
RE: Rotating machinery imbalance
What is the format of the raw data? 30-degrees phase is not very close,unless it is a consistant shift that can be compensated by calculation. Does your vibration sensor or speed sensor have a phase lag that is different from the "commercial balancer"?
Walt
RE: Rotating machinery imbalance
Cheers
Greg Locock
New here? Try reading these, they might help FAQ731-376: Eng-Tips.com Forum Policies http://eng-tips.com/market.cfm?
RE: Rotating machinery imbalance
The 30 degree discrepancy appears to be quite repeatable. I suspect it has to do with the integrated Bessel filter in the MEMS accelerometer used to capture one data set, vs. the ICP accelerometer with hardware filtering used in the commercial balancer to capture the other data set.
A representative data set is attached.
The data parameters are:
sample rate = 1mS
accelerometer sensitivity = 1200mV/g
magnitude of the raw data is in volts.
Tell me the velocity in inches per second, and the phase in degrees.
RE: Rotating machinery imbalance
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Rotating machinery imbalance
Attempting again.
RE: Rotating machinery imbalance
We'll see if it remains this time.
RE: Rotating machinery imbalance
Hi Greg
I'd like to hear anything you can recall about this.
RE: Rotating machinery imbalance
The file cannot be downloaded because the path is on your computer and not accessible.
Walt
RE: Rotating machinery imbalance
Should be there now.
RE: Rotating machinery imbalance
Cheers
Greg Locock
New here? Try reading these, they might help FAQ731-376: Eng-Tips.com Forum Policies http://eng-tips.com/market.cfm?
RE: Rotating machinery imbalance
All we could derive from this would be the phase between the pulse and the same frequency of the data channel, it tells you nothing about the balance.
Cheers
Greg Locock
New here? Try reading these, they might help FAQ731-376: Eng-Tips.com Forum Policies http://eng-tips.com/market.cfm?
RE: Rotating machinery imbalance
Both channels have to be sampled simultaneously, therefore the time record should have the same number of samples for both channels (tacho and vibration). The tacho and accelerometer angular positions need to be known relative to the known heavy spot to verify that the vibration vector is correct.
Walt
RE: Rotating machinery imbalance
The enigmatic filename, CH13.csv was simply generated by my Tektronix scope (meaning channels 1 & 3).
My embedded algorithm counts complete revolutions in the keyphasor channel and disregards any data past the last full revolution, or 5000 samples, whichever happens first.
There are 16 full shaft revolutions represented in this data set. Each positive "pulse" in the keyphasor channel is the start of a new revolution (0 degrees).
As before, the first column is the keyphasor strobe data, and each sample in this column corellates directly to the sample in the second column, which is the raw acceleration data.
The sample rate is 1mS, the accelerometer sensitivity is 1200mV per g, the amplitudes for both channels are in volts.
Each pair of samples is digitized "simultaneously" (some insignificant latency applies), with a 12-bit A/D.
RE: Rotating machinery imbalance
Your data is very dirty, it would probably be a good idea to use a low pass filter on both channels.
I am a bit baffled, you appear to have an 8 bladed helicopter!
Cheers
Greg Locock
New here? Try reading these, they might help FAQ731-376: Eng-Tips.com Forum Policies http://eng-tips.com/market.cfm?
RE: Rotating machinery imbalance
Need velocity in IPS please.
Greg wrote: "energy correction"
Huh?
Greg wrote: "phase of 1.93 radians"
Close, 1.93rads = 110.5 degrees. I show 113 degrees.
Greg wrote: "Your data is very dirty"
Yes it is, and I'm not sure why. The MEMS accelerometer has an integral 2-pole Besel filter with nominal roll-off at 50Hz. Frequency of interest is 8.7Hz (fundamental shaft rotation freq), and I'm sampling at 1Khz. Don't care about the full spectrum. The big spike at 70Hz coincides with engine RPM (4200RPM). The other at 250Hz is a mystery at this point. That's well beyond the accelerometer roll-off. Might be some noise on my circuit board.
Greg wrote: "I am a bit baffled, you appear to have an 8 bladed helicopter!"
Nope, can only afford two blades :)
RE: Rotating machinery imbalance
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Rotating machinery imbalance
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Rotating machinery imbalance
The phase is very dependent on the details of how the fft is done. I used synchronous sampling at several different resolutions, the number I gave was the average for the different resolutions.
revs per frame/Vrms/phase
2.000000 0.054359 1.789612
3.000000 0.055126 1.868418
4.000000 0.057182 1.905492
5.000000 0.058639 1.937415
6.000000 0.059522 1.952189
7.000000 0.059790 1.955220
8.000000 0.059752 1.954846
9.000000 0.059860 1.958020
10.000000 0.060458 1.969961
11.000000 0.060310 1.959680
12.000000 0.059962 1.954206
13.000000 0.059732 1.923215
14.000000 0.059209 1.904673
15.000000 0.059458 1.885334
16.000000 0.059967 1.899490
Cheers
Greg Locock
New here? Try reading these, they might help FAQ731-376: Eng-Tips.com Forum Policies http://eng-tips.com/market.cfm?
RE: Rotating machinery imbalance
The pattern as Greg said has a lot more going on than simple once per revolution sinusoid. Additionally I don't think it's even periodic at once per revolution (there seem to be some erratic changes). One might think the correction would involve something other than simple balancing.
Since you mentioned 30 degree error. I noticed the width of the pulse is 0.014. The separation between pulses is 0.114. The width of the pulse is 44 degrees. Have you consistently worked with the rising edge, middle or trailing edge?
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Rotating machinery imbalance
Cheers
Greg Locock
New here? Try reading these, they might help FAQ731-376: Eng-Tips.com Forum Policies http://eng-tips.com/market.cfm?
RE: Rotating machinery imbalance
Slide 2 shows the magnitude and phase computed based on a rolling one-period window. There is plenty of variability. If the only thing present were harmonics of the fundamental, the result would be stable. It is far from periodic at rotating speed.
Slide 3 shows the magnitude and phase computed based on a rolling 10-period window. There is still some variability as the window slides.
My results roughly agree with Greg. I would eyeball from slide 3 the magnitude is around 0.6 and the phase is around –1.9 radians. I used a lagging phase convention, I assume Greg used a leading phase convention.
Below is the code FWIW. It is based on observation that one rotation corresdonds to 114 samples and the rising edge phase reference occurs at 0.181 +/- Integer*Rotation. I set Nperiods between 1 and 10 for slides 2 and 3 and respectively
CODE
' yrange contains the beginning of the range for values
' trange range contains the beginning of the range for time
' Constants for this problem
Const WindowLength As Integer = 114 ' Number of samples in on rotation period
Const dt As Double = 0.001
Const tref As Double = 0.181 ' Reference time at Rising Edge.
' Phase is 0 if the sinusoidal peak occurs at this time
' phase is positive if sinusoidal peak occurs slightly after this time (LAGGING PHASE CONVENTION)
Const Nperiods = 10 ' Number of Periods to use for averaging
Dim counter As Integer
Dim t As Double, y As Double
Dim a As Double ' cos term of fourier series
Dim b As Double ' sin term of fourier series
Dim pi As Double
Dim w As Double ' radian frequency corresponding to rotating speed
pi = Application.pi ' A constant
w = 2 * pi / (WindowLength * dt) ' Rot speed in rad/sec
a = 0 ' initialize cos term
b = 0 ' initialize sin term
For counter = 1 To WindowLength * Nperiods
t = trange.Offset(counter - 1, 0)
y = yrange.Offset(counter - 1, 0)
a = a + y * Cos(w * (t - tref))
b = b + y * Sin(w * (t - tref))
Next counter
' Find magnitude by combining cos and sine components
Mag = 2 * Sqr(a ^ 2 + b ^ 2) / WindowLength / Nperiods
' Find the phase..
' atan2 function takes argument x, y
' find angle as a LAG angle. Cos is 0 ref, Sin lags cos.
Dim Ang As Double
' Ang = Application.WorksheetFunction.Atan2(a, b)
End Function
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Rotating machinery imbalance
Notes
1 - It is in manual recalculation mode so you have to prese F9 before any recaclulation occurs.
2 - Note there are provisions for rescaling and browsing the data in the plot tab.
3 - The label "degrees" in the angle is wrong.. should be radians.
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Rotating machinery imbalance
The line of interest is the first order. You can see why I thought the rotor must have 8 blades.
Cheers
Greg Locock
New here? Try reading these, they might help FAQ731-376: Eng-Tips.com Forum Policies http://eng-tips.com/market.cfm?
RE: Rotating machinery imbalance
0.06 volts peak/0
Using 1.2 volts/g conversion, that equates to 0.05 g's pk/0
At a frequency of 8.77 hz, that corresponds to 0.35 ips pk/0
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Rotating machinery imbalance
Who woke you guys up? :)
The more I learn, the less I know. I have no expertise in this subject matter so I rely on others'.
I purchased the embedded DFT, and my objective was to get some corroboration on the result report, and to get a "tie breaker" between the report my algorythm is generating and that from a commercial balancer (which may well be out of calibration).
Greg's phase and velocity report is quite close to mine:
Greg's 0.34ips at 110.5 degrees
Mine 0.35ips at 113 degrees
Thank you Greg & Pete
As for Pete's coment on the periodicity, run a 7Hz to 10Hz bandpass filter on the accel data (2nd column) and the periodicity jumps right out.
ElectricPete wrote: "One might think the correction would involve something other than simple balancing"
Nope. The velocity and phase report specifies the location of a "dot" on a polar chart, specific to a particular aircraft. Striking a line from this dot to notes on the margins of the chart, point to the "corrective action". Apply the corrective action, then take more data. The next report will indicate the "quality" of the corrective action. If the velocity magnitude decreases, the dot is closer to the center of the chart. The center of the chart is perfect (un-achievable) balance. Repeat process until best balance is reached.
This is how every rotor/propeller balancer (icluding Chadwick-Hellmuth) I've used works.
RE: Rotating machinery imbalance
My comments were:
1 – The signal is not stable. Stable may be a vague term, but where I was headed is that it does not seem valuable as a calibration reference since the 1x component seems to vary over time.
2 – The entire waveform is not periodic at 1x (if it were the higher frequency content would not influence the computation of 1x in my algorithm as they did).
Further on discussion of 1, I performed filtering.
Slide 1 is original waveform and LP-filtered waveform.
I chose Kaiser Bessel 12hz cutoff, 5hz transisition width (about 2.5hz on each side), Allowable deviation at ends of transition width: 0.05, Resulting FIR Length: 253 Samples (output has half that many samples trimmed from each end). It does not preserve phase, but good enough to judge magnitude and phase stability imo.
Original and LP filtered waveform in slide 1.
Zoom-in on LP filtered waveform in slide 2.
Slide 3 is result of LP filtering. Since I know the periodicity of the signal of interest, I performed LP filtering by subtracting out a rolling average of length 115 samples, centered on the point of interest. You can already see substantial variation in the peak-to-peak values, more than just expected by result of sampling effects missing the true peak.
Slide 4 is sliding 1 period computation. The magnitude floats from 0.06 to 0.07 and phase floats –2.5 to –2.7 (again magnitude of phase no longer correct due to filtering, but variability is evident).
Slide 5 is zoom-in of the 1-period computation. Again you see magnitude vary 0.06 to 0.07. This is what I meant by not stable
Slides 6 and 7 are sliding 10-period computation. There is a variability 0.064 – 0.066 which has some sinusoidal 1x component in it.
Reviewing my algorithm:
1 – We do have very good knowledge of the fundamental period. There are 16 pulse rising edges and they are all separated by exactly 114 samples.
2 – The period is (114+1)*dt. 114 is the correct length for the sum since we don't want to count both endpoints (only one or the other). However I think I should have used 115 when computing w... i.e. should have used w = 2 * pi / ((WindowLength + 1) * dt).instead of w = 2 * pi / ((WindowLength) * dt)
3 – When I "corrected" w to be 2 * pi / ((WindowLength + 1) * dt) on slides 8 and 9, the once per revolution variability significantly unexpectedly increased. That was a surprise to me. I conclude for some reason using w = 2 * pi / ((WindowLength ) * dt) gives a better estimate than w = 2 * pi / ((WindowLength + 1) * dt), but I honestly don't know why. I would be interested to hear comments on why that is.
4 – Filtering seems to have slightly increased the magnitude unexpectedly – I don't know why that happened, but I'm not going to worry about that.
Regardless of those questions, I am left believing your signal has variability in the 1x component. I can understand you want to validate your algorithm using some kind of realistic signal, but I would think for calibration purposes you would want to use a signal with a repeatable (stable) 1x component.
I didn't say this is not the way to balance. My comment was that some other corrections beyond simple balancing may be appropriate.
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Rotating machinery imbalance
"Slide 3 is result of LP filtering..."
should have been
"Slide 3 is result of HP filtering..."
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Rotating machinery imbalance
There is a measure called 'coherence' that is an indication of the quality of your data, and would allow you to select the best compromise between resolution and total sample time.
It is defined as
sqrt((mod(Fxy))^2/(Fxx.*Fyy))
Where Fxy is the transfer function (FRF) between channels x and y
This I have always understood to be the ratio of the length of the vector sum of the averages to the total of the magnitudes of each vector, at each frequency.
If the signals are correlated you get 1, if not it drops towards zero. Essentially it is a measure of the phase stability of the transfer function at the frequency of interest.
Cheers
Greg Locock
New here? Try reading these, they might help FAQ731-376: Eng-Tips.com Forum Policies http://eng-tips.com/market.cfm?
RE: Rotating machinery imbalance
2.000000 0.054359 1.789612 0.941900
3.000000 0.055126 1.868418 0.978263
4.000000 0.057182 1.905492 0.989251
5.000000 0.058639 1.937415 0.995448
6.000000 0.059522 1.952189 0.996754
7.000000 0.059790 1.955220 0.995691
8.000000 0.059752 1.954846 0.997073
9.000000 0.059860 1.958020 0.997832
10.000000 0.060458 1.969961 0.999614
11.000000 0.060310 1.959680 0.999743
12.000000 0.059962 1.954206 0.999499
13.000000 0.059732 1.923215 0.999489
14.000000 0.059209 1.904673 0.999327
15.000000 0.059458 1.885334 0.999931
16.000000 0.059967 1.899490 1.000000
From this I'd say that 0.2-0.33 (ie 3-5 revs per frame) orders is a good compromise between enough resolution to get good phase stability, and enough averages.
Cheers
Greg Locock
New here? Try reading these, they might help FAQ731-376: Eng-Tips.com Forum Policies http://eng-tips.com/market.cfm?
RE: Rotating machinery imbalance
Also the integration from acceleration to velocity causes a 90 degree phase shift, so I'm not really working in velocity, I have rescaled the volts into velocity magnitude at that frequemcy, not the same thing.
Cheers
Greg Locock
New here? Try reading these, they might help FAQ731-376: Eng-Tips.com Forum Policies http://eng-tips.com/market.cfm?