CC BY
50Electronic Journal «Technical acoustics» http://webcenter.ru/~eeaa/ejta
2 (2002) 6.1-6.15
A. G. TROSHIN
"Izmeron Ltd" Tool; Plant, Laboratory of linear measurements and electronic devices.
Novgorodskaya street 13, 193144, Saint-Petersburg, Russia, e-mail: eeaa@online.ru
Structural energy flows in water filled pipes: implementation of measurement technique based on PVDF reusable strain gauges*
Received 17.05.2002, published 01.07.2002
The purpose of this work is an implementation the newly developed measurement technique for structural energy estimation in-situ. The technique based on reusable PVDF strain gauges in measurement of structural energy flows transferred by longitudinal and shear forces, bending and torsional moments in water filled pipes are presented and discussed. Developed technique was implemented for determination of the efficiency of noise control means in the pipeline with the test rig «Flow». The technique also applied in order to determine the components of vibrational energy flows generated by pump incorporated into the rig «Impedance». Measurements were performed when real pumps have been run. The relationships between four structural components in pipe walls and water ones transferred with plane waves in the pipe channel were determined. Degree of uncertainty in proposed technique was revealed using error analyses taking into account the real data about coherence function and measured phase angle between force/moment and linear/angular velocity. The reliability of preformed measurements at tonal components of structural energy spectra proofed employing an error analyses and energy balance calculation in pipe cross section. The possibilities and limitations of the technique were addressed.
1. INTRODUCTION
Experimental determination of vibrational energy flows is a powerful tool in the modern applied acoustics. All energy transfer paths are forming the total vibrational or sound field in transportation (car, ship aircraft) compartments and therefore the lack of data about just single energy transfer channel makes the situation uncertain. User should has at hand entire information concerning all transfer paths such as machine resilient supports, pipeline hangers and pipe walls, air medium and plane waves propagation in liquid filled pipes. The appropriate experimental technique should be applied at each individual case. Energy flow approach opens a number of additional options for designer during development of noise control strategy. Having at hand data about energy flows it is possible to judge in straightforward manner, which transfer path dominants when vibrational energy transmits from machine (source) to foundation (receiving structure).
J. Verheij [1], R. Pinnington and R. White [2] and others proposed the technique for estimation of energy flows propagated via machinery supports and pipe hangers. Acoustic
*This paper had no linguistic support from English speaking expert. Editorial staff is not responsible for lingual mistakes in the article
power radiated in the form of air borne noise can be also measured by means of commercially available measurement equipment and appropriate software. Double or triple hydrophone methods which are resemble to Fahy [3] and Chung [4] technique for air borne noise can be also successfully applied in measurement of energy component propagated throughout the liquid channel of energy transmission in water filled pipes.
Methods for measurement of energy flow in pipe wall structure were proposed by J. Verheij [1] and G. Pavic [5]. The methods were based on conventional transducers array (accelerometers) and finite difference technique. The methods could be implemented in specific cases when structure (beam or rod) straight enough because the measurement of longitudinal component of energy flow requires a long measurement distance due to high speed of longitudinal waves. A measurement of bending wave component using 4-point technique can be performed in narrow frequency range due to its high sensitivity to phase mismatch between measuring chains. The chain implies here a measurement channel including transducer, preamplifier and FFT analyzer. The two or three point technique are disregarded the effect of evanescent waves whose contribution could add 3 or 6 dB into the total measurement uncertainty.
Developing a new generation of electromechanical transducer — reusable PVDF strain gauges — reveals the number of advantages for measuring the vibrational energy flows, vibration induced strains, forces and moments in pipe and shaft structures by measuring surface strain. The PVDF strain gauges produce the electrical charge which is directly proportional to dynamic strain and therefore it is possible to use them jointly with traditional measurement equipment, namely, charge preamplifiers and modern FFT analyzer if their channels are equipped according to ICP standard. The PVDF strain gauges also make possible the measurement of extremely small magnitudes of strain about 10-5 (|im/m according to
denotation accepted by Russian standard). In fact it can be explained that PVDF film more elastic than wire and quartz. Therefore the film is able to produce more powerful signal. Despite lower piezoelectric modulii than ceramic and quartz the cross coupling mechanic to electric efficiency of PVDF film is more than traditional peizoelectrically active material. This phenomenon makes film very useful in designing the transducers, actuators and other practical applications.
The reusable PVDF strain gauges allow to measure strain and force directly which make possible to reduce the order of derivative approximated by FD and give an opportunity to alleviate drawbacks of parametric methods with conventional wire strain gauges. The latter become poor for measurement of small dynamic (vibrational) strain. The advantage of reusable PVDF strain gauges is their high sensitivity to strain and possibility to calibrate and use them after having at hand the sensitivity vs. frequency as for usual vibration transducer. The formulae for measurement of individual components in total energy flow become simpler in comparison with traditional FD approach. This circumstance gives an additional advantage for overcoming the specific error inherent to FD technique for example narrow frequency range and long measurement distance. All techniques discussed earlier (see for example [6] and [7]) were based on single used PVDF gauges and they had implemented in so-called «academic structures» - free-free beams and beams with anechoic ends. This paper contains experimental results and their error analyses obtained with real pipeline configurations and running pumps. This configuration allows revealing the advantages and drawbacks of proposed technique.
2. METHODS FORMULATION AND TRANSDUCER ARRAYS
Measurement technique for vibrational energy propagated via pipe wall structure based on mathematics model, which describes a structure oscillations by use of Euler-Bernoully model. In other words the beam layer become plane and parallel each other due to shear force action (there is no shear strain dependence on beam thickness and rotational inertia of beam element also negligible). The used symbols and notations are given in Appendix A.
Error of the technique at upper cut-off frequency due to implementation of above model can be calculated by formula (1) obtained by Nikolay A. Kuznetzov (Krylov Institute).
8 =
1 -
P h,
high
ES
1-
E
<pG
+1
+ n fh
high
P h
ES
1+
E
<pG
(1)
6(1 + v )(1 + n2)2
n ■■
Dn
(7 + 6v)(1 + n2)2 + (20 + 12v)n2
- diameter ratio; v - Poisson ratio of pipe material.
- form quotient for annular cross-section of the pipe;
Douter
Actually, the equation (1) obtained as a difference between bending wave number derived according to rigorous solution for equation of beam motion taking into account rotational inertia and thickness shearing of cross section and Euler-Bernoully’s model.
The plots of the error due to implementation of the theory as a function of pipe geometry and frequency are shown in Fig. 1.
Another assumption is that the measurement part of the pipe should be straight and uniform. The interaction between liquid medium and pipe structure supposed to be negligible employing the mathematics model in measurement technique. Expression (1) demonstrates that technique based on PVDF application has also high frequency limitation due to model involved in the technique.
2
2
Fig. 1. Error due to technical model vs. frequency and parameters of pipe cross section
(—) - inner diameter 50 mm, wall thickness 5 mm; ( ) - inner diameter 100 mm, wall
thickness 5 mm; (.....) - inner diameter 150 mm, wall thickness 5 mm;
( - inner diameter 200 mm, wall thickness 5 mm
Forces and moments acting in a beam cross section can be expressed by well-known equations:
Ff■ (2)
M'= GJos dx ’ (3)
d M w
F --df ■ «
Mv=-EIbend • (5)
The partial derivatives in the equations (2)-(5) can be extracted using data about surface
strains (f) > S x (f )> ¿3* (/)> S4 x (f )and disPlacements of beam structure £ly (f), ^ y (f),
Zy (f), ^4y (f):
F (f ) = ES
^1, (f ) + ^2x (f )'
2
v
(6)
FT
K(f) = -D=‘L(eJ.f)-£,-(/))■ (7)
outer
FT
F (f ) = - A T - u (f ) ~e‘- (f ) - ^• (f ) + e*- (f ))■ (8)
bend outer
M„(f ) = GB1 (y (f )-fey (f )-ily (f )+¿4, (f ))• (9)
tors outer
For linear and rotational vibration velocities the formulae are:
v (f ) = LftMfl ■ (10)
2jrn
v, ( f ) = i^f) ■ (11)
jrn
/ (f ) = i-(f )—¿2-(f ^ ; (12)
nouter
,(f ) =f2f ■ ^ ) = fUf) ■ (n)
j®Douter j®Douter
d(p Ç, -(p2 ; dp = (f) -^2y (f ) -^3y (f) + ^4y (f))
d x Aors ’ d x Ator,A»*r
(14)
Arrangements of the strain gauges for measurement of the strain e1x (f), e2 x (f), e3 x (f), S4x(f) and accelerometers with double integration their signals for measuring the displacements |1 (f), ^2 (f), ^3y (f), ^4y(f) are shown from Fig. 2 to Fig. 4 including block diagrams of measurement set-ups.
Fig. 2. Arrangement of PVDF gauges at pipeline for measurement of flow components due to
longitudinal force and bending moment
1 - charge preamplifiers; 2 - PVDF strain gauges; 3 - accelerometers; 4 - FFT analyzer
Fig. 3. Arrangement of PVDF gauges at pipeline for measurement of flow components due to
shear force
1 - charge preamplifiers; 2 - PVDF strain gauges; 3 - accelerometers; 4 - FFT analyzer
Fig. 4. Arrangement of accelerometers at pipeline for measurement of flow components due to
torsional moment (Verheij’s method)
1 - charge preamplifiers; 2 - accelerometers; 3 - active summation and subtracting device; 4 -
FFT analyzer
Expressions for determination of the vibrational energy radiated due to harmonic forces and moments are:
WPx (f ) = iRefc; (f ) < (f )}; W, (f ) = |Re{Fz (f ) • v;(f )};
WMv (f ) = |ReK (f ) V; (f )};
wbending (f ) = WK (f ) + Wm,„ (f );
WMif (f ) = |ReM,(f ) •<P; (f )}.
(15)
(16)
(17)
(18) (19)
Substituting the forces/moments and linear/angular velocities to Eq. (15)—( 19) and applying Fourier transform it is possible to derive the formulae for determination of the energy flows in terms of cross spectra as it given further on:
For component of longitudinal force
ES
.............. (20)
where Aex = £ix +£2* and = H* +<*.
For component of bending moment
WM (f ) = E1bend {im Ga , (f )}
(21)
where = £1x -^2x and V = <L -%2x •
There are two advantages of technique employing the equation (9) and (10). One of them an opportunity to make measurement both longitudinal and bending moment components simultaneously. Another one the lack of lower and upper frequency limitation due to finite difference method because of direct measurement of longitudinal force and bending moment. For shear force component
WQ:(f) = D EIT m{,mG£I,(/)}• (22)
outer bend where Aez = £ - £x - £x + £4x .
For component of torsional moment
G j
WM' (f) = Jm G^(f)’ (23)
outer tor
where X1 = %1y + E)1 - instantaneous sum of vibrational acceleration from opposite positioned accelerometers at the points 1 and 2, X2 = %3y + E,3y - instantaneous sum of vibrational
acceleration from opposite positioned accelerometers at the points 3 and 4, respectively.
The analog devices have to be used for measurement performing with PVDF strain gauges and accelerometers. The devices provide more precise result vs. post processing of averaged spectral functions as it shown in [8]. The analog device for simultaneous summing and subtracting of time domain signal was designed and tested by author at Chalmers University of Technology [9]. The spectra of vibrational and hydrodynamic noise generated by pumps comprise powerful tonal components which usually concentrated in frequency band from 5 to 400 Hz with 90% of their energy. Therefore a frequency analysis of vibrational power is limited here by upper frequency 400 Hz in this particular case. Such upper frequency limit has to be chosen also due to limitation of Euler-Bernoully model at high frequencies.
3. ENERGY FLOWS ALONG RIG PIPELINES
3.1 MEASUREMENT OF ENERGY FLOWS WITH «FLOW» RIG
The purpose of this experiment to determine an efficiency of applied noise control means for pipelines using parameter of vibrational energy. A scheme of the measurement part on the rig is depicted in Fig. 5.
Fig. 5. Arrangement of measurement point for monitoring of structural energy flows
Four cross sections of the pipeline were measured during experiments with the running pump. Vibrational energy transferred by plane waves in water medium was also measured using Fahy’s technique with two hydrophones inserted into water feeding and intake pipes
[3].
The measurement results are shown from Fig. 6 to Fig. 8. Energy was measured at point 1,
3 and 4. The energy emitted into the resilient mount at point 2 was estimated using complex transmissibility technique [2] taking into account all six degree of freedom, which means that rotational stiffness values of isolator were measured and substituted into equations of energy estimation. The energy balance was calculated using equation (24) and final result is shown in Fig 7.
Vibrational energy dB ref. 1e-12 Watt
Frequency (Hz)
Fig. 6. Energy balance in pipe cross section Eq. (24) (—) - energy input (—) - energy
transferred
Here we assumed that water borne energy losses is negligible (there were no any special means applied) thus only structural components were considered.
11/input | jjrinput | tjjinput ^ jjrtranst , jjrtranst , jjrtranst , tjr total .x
Wlong + Wbend + Wtors — Wlong +W bend + Wtors + Whanger , (24)
We can see that inequality (24) becomes true over all frequency range however it is impossible to get the exact energy balance due to transformation of structural component of the energy to the air-borne noise and bending waves with polarization out-of-plane. The air borne noise and bending waves in another plane were not monitored during this experiment.
The apparent violation of energy conservation law at frequencies about of 70 and 260 Hz (see Fig. 6) is due to the back flow from water tank. Hanger substantially reduces the forward waves and flow due to the injection branch of the pipe becomes dominant.
As an example how to implement the proposed technique for estimating the efficiency of structural vibration reduction by flexible hanger and hose in energy terms is depicted in Fig. 7 and Fig. 8. Figures 7 and 8 show that hanger more efficient here comparing to the flexible hose.
Vibrational energy dB ref.1e-12
Frequency (Hz)
Fig. 7. Structural energy reduction by flexible hose,
(—) - power input, (- - -) - power transfer
Vibrational eneigy dB ref. 1e-12 Watt
Frequency (Hz)
Fig. 8. Structural energy reduction by hanger, (—) - power input, (- - - ) - power transfer
3.2 MEASUREMENT OF ENERGY FLOWS WITH «IMPEDANCE» RIG
This experiment was made to reveal the relationship between different components of energy flow generated by pump in situ in order to describe it as a source of water and structure borne sound in energy terms. A transducer arrangement is depicted in Fig. 9.
O
Fig. 9. Transducers arrangement at the pipe of “Impedance” rig
In this specific case the components of shear and bending moment are almost equal each other in a wide frequency range. Hence this fact makes possible the usage of far field approximation when total bending energy flow can be represented in form of double component of the bending moment (see Fig. 10)
Vibrational energy dB ref.1e-12 Watt
Frequency (Hz)
Fig. 10. Shearing force (-------) and bending moment ( — ) components
This circumstance could make measurement procedure easier and there is no reason to use the shearing force measurement arrangement. Figure 11 demonstrates that energy transferred via water of 10 dB higher than total structural component at frequency which proportional to product of two values, namely, the pump rotor number of revolutions and number of propeller blades.
Vibrational energy dB ref.1e-12 Watt
Frequency (Hz)
Fig. 11. Structure (—) and water borne (-------------) components
At other frequencies vibrational energy transferred by structural components only. The longitudinal component dominates in total structural energy flow due to the rigid connection of the pump and attached pipelines. (see. Fig. 12).
Vibrational energy dB ref.1e-12 Watt
Frequency (Hz)
Fig. 12. Three structural components: ( — ) - longitudinal, (..................) - bending
and (-------) - torsional
4. ERROR ANALYSES USING ACTUAL MEASURED DATA
Here we shall consider two types of error: one due to low signal to noise ratio or low value of coherence function between force/moment and linear/angular velocities. The second type is error due to ratio of phase mismatch between measurement chains and measured phase angle between force and linear velocity (moment and angular velocity). If the measurement error can exceed more than 100% let us put into consideration the term - “measurement uncertainty”.
The uncertainty of first type has a random nature and depends on the coherence function and phase angle. Seybert derived the expression for the random error [10] in the case when two microphones method is implemented when sound intensity measurement technique proposed by Fahy [3] is used. Seybert’s approach is universal and expression for random error can be applied for all cases where imaginary part of cross spectrum is used as a primary data in vibrational energy estimation. The formulae for bias (instrumentation) error are the function of phase mismatch between force/moment and linear/rotational acceleration channels. Here we assumed that real measured phase angle or its smoothed estimation has to be taken into account. The error can be written in the form
0
bias
Aa
tana
(25)
Another type of error is random one, which can be expressed by Seybert’s equation:
random
1 I 1
nd
fPv
+
/i 2 \ 1-Y2
n.1/2
I Fv
2Y
Fv
(tana)2
(26)
Here and after the a angle is actual one obtained during measurement as a function of frequency. Usually there are numerous singular frequency points where phase angle becomes small as well as coherence function. It means that energy transferred by corresponding waves is small due to lightly dampening of the structure or high reflections occur. The low coherence function shows that components of the wave motion namely, wave magnitudes, are also small.
Vibrational energy dB ref. ie-12 Watt Fig* I3* Structural energy
spectrum and its bias error
Bias error, dB
A. - I.-! i L
JlL , J.11
50
100
150 200 250
Frequency (Hz)
300
350
400
1
5
0
Vibrational energy dB ref. 1e-12 Watt
Random error, dB
Frequency (Hz)
Fig. 14. Structural energy spectrum and its random error
(500 independent samples)
If the structure is lightly damped the uncertainty of measurement will rapidly increase. The results of the estimation are depicted in Fig. 13 and Fig. 14. It can be seen that bias error does not exceed ±2 dB. The uncertainty of measurement rapidly increases at frequencies where signal to noise ratio is small (see Fig. 14) In other words the random error dominates.
The uncertainty can reach values of 10-30 dB at frequencies where energy values smaller
—12
than 30 dB referring to the threshold level 10 Watt. It should be noted that real values of vibrational energy that can be estimated correctly (with acceptable uncertainty ±3 dB) have to
_9
be more than 10 Watt (30 dB) when PVDF strain gauges sensitivity are of 30 mV per micro strain. There is only one way for error reduction. It is increasing a number of samples, which limited by observation time and time interval where vibration of machine can be considered as steady state random process. A logical connection between Seybert’s principle in intensity measurements and uncertainty principle in quantum mechanics could be made. Both of them do not allow getting the desirable accuracy by measurer will.
DISCUSSION
Technique involved reusable PVDF strain gauges was developed and discussed in the paper. In contrast to FD technique the proposed one has number of advantages that allow to make correct estimation of structure borne vibrational energy in situ with running pumps and real pipe configurations. High frequency limitation due to involved mathematical model has calculated in each individual case. The error analyses with actual measured data about coherence function and phase angle of cross spectra allow to determine the low tolerable limit
—9
of energy estimation. It is approximately 10 Watt involving this measurement technique and attainable PVDF gauge sensitivity. The sensitivity can be increased by use of ICP electronic circuit. The preliminary experiments have shown the attainable sensitivity with ICP about of 300—500 mV per micro strain.
It is shown that pump as a source of vibrational activity can be described by use of waterborne and structure-borne components and sometimes it is not correct to consider both transfer channels separately or skip one of them from transfer path analyses. This
phenomenon was investigated with the «Impedance» rig. The relationships between four structural components and water ones transferred by plane waves in the pipe channel were determined and discussed. In the case with «Impedance» rig a structural component transferred mostly by longitudinal waves in the pipe walls whilst water borne sound dominates at several frequency. Using the developed technique an efficiency of noise control means applied to «Flow» rig in energy terms was estimated. The efficiency of resilient hanger becomes more than efficiency of flexible hose in situ. The energy balance in pipe cross section and error analyses shows a high confidence level of obtained data. Moreover it can be declared that total error in vibrational energy estimation at the tonal components of energy spectrum does not exceed ±3 dB where coherence and phase angle between force/velocity, moment/angular velocity are high enough.
ACKNOWLEDGEMENTS
I would like to express my gratitude to Urij I. Kuznetzov from Krylov Institute for his help in measurement of water borne energy component propagated via pipe medium. His comments, friend attitude and personal participation during experimental testing are gratefully acknowledged. Prof. Alexei S. Nikiforov’s valuable comments and Dr. Igor Grushetsky’s support during preparation of the manuscript are very appreciated.
REFERENCES
[1] J. W. Verheij. Multi path sound transfer from resiliently mounted shipboard machinery. Teschnisch Physische Dienst TNO-TH Delft, 1982.
[2] R. J. Pinnington, R. G. White. Power flow through machine isolators to resonant and non resonant beams. J of Sound and Vib. 75 (1981) 179-197.
[3] Fahy F. J. Measurement of acoustic intensity using the cross-spectral density of two microphone signals. JASA 62(4) 1057-9, 1977.
[4] Chung J. V. Cross-spectral method of measurement acoustic intensity without error caused by instrument phase mismatch JASA 64(6), 1613-1616, 1978.
[5] G. Pavic. Technique for determination of vibration transmission mechanism in structures. Dissertation. Institute of sound and Vibration Research, Faculty of Engineering and Applied Science, University of Southampton, 1976.
[6] A. G. Troshin and V. I. Popkov. Measurement of Vibration Power flow in rod structures by using piezo-electric film sensors. 4-th International Congress on Structural Intensity Technique, Senlis, France, 1993, pp. 169-174.
[7] C. R. Helkyard, B. R. Mace. A wave approach to structural intensity in beams. 4-th International congress on intensity technique, Senlis, France, 1993, August 31—September 2, pp. 183-191.
[8] A. G. Troshin, M. A. Sanderson and L. Ivarsson. Vibration Isolation of Structure borne power transmission by fans in Building: a theoretical and experimental examination using Mobility and Complex Transmissibility methods, Report S 97-06, Gothenburg, December 1997, Department of Applied Acoustics, Chalmers University of Technology.
[9] A. G. Troshin. Summation-subtraction device for a six degrees of freedom of motion transducer comprised of six linear accelerometers Gothenburg, September 1997, Report S 97-
05, Release Note Department of Applied Acoustics, Chalmers University of Technology.
[10] A. F. Seybert. Statistic error in acoustic intensity measurement. J of Sound and Vib. 75 (1981)585-595.
APPENDIX A: SYMBOLS AND NOTATION
E Young modulus
G Shearing modulus
V Poisson ratio
P Material density of pipeline or shaft
$ , V , Vibrational displacement velocity and acceleration
F Longitudinal force
Fz Shear force
K Bending moment
Torsional moment
£m, Dynamic strain at point m along axis i
Douter Outer pipe diameter
Dmner Inner pipe diameter
bbend Moment of inertia for ring cross section refer to axes y and z
J tors Polar moment of inertia of annular cross section
A Distance between accelerometers for measuring a component of
tors torsional moment
A , Distance between PVDF strain gauges for measuring a component
bending of shear force
Yfv Coherence function between force and velocity
6 Fv Phase angle between force and velocity
Aa Channels phase mismatch
nd Number of samples
Gab (f) Complex valued cross spectrum between processes A(t) u B(t)
^bias Bias measurement error
^ rand Random error
W Vibrational power
LW Vibrational power level
Im, Re Real and imaginary part of complex value
S Annular cross section area
(O Angular frequency
f Frequency
