This invention relates to proximity sensors, such as eddy current sensors, and to the processing of signals from such sensors. In particular it relates to the application of such sensors to the measurement of position in relation to objects such as turbine blades and the like.
Gas turbine engines employ sets of turbine blades mounted on rotatable shafts, typically with one or more sets of blades acting as a compressor, feeding air into a combustion chamber, and one or more sets located behind the combustion chamber, acting to power the compressor. Typically there may be between 20 and 200 blades forming a compressor or combustion turbine, and shafts (or spools) may rotate typically at 10,000 revolutions per minute (RPM). To achieve reasonable efficiency it is beneficial in many circumstances for the outside edges of the turbine blades to be within a certain distance from the cowling surrounding the turbine, this being typically between 0.5 and 5 mm. Measurement of this distance while the engine is running is therefore very useful in that the result can be used when adjusting the blade to cowling distance to achieve optimum efficiency. The measurement can also be very useful in detecting (and, if required, compensating for) turbine blade “growth”, which occurs due to the centripetal forces acting on the blades when rotating at high speed, or in detecting overall changes in turbine tip to casing clearance. Defective or incorrectly fitted blades may also be detected, should they unilaterally change height.
Eddy current sensors are commonly in use for making the measurements described above. U.S. Pat. No. 5,942,893 describes such an eddy current sensor that may be used in measuring turbine blade clearance. In use, the sensor is located in the cowling so that it is in close proximity to the outer edge of the turbine blades. The sensors contain one or more coils. The sensors contain a means to create a magnetic field. The field extends into the region through which the objects to be sensed pass. The means to create the field can be either a permanent magnet (for example a rare-earth magnet), or a current passing through one or more of the coils. The field may be either DC or AC. A permanent magnet or a DC current creates a non-uniform DC field. Alternatively, an AC current creates an AC field. Eddy currents are induced in a conducting object as it experiences changes in the magnetic field that surrounds it. In the case of a DC field, the object experiences field changes as it moves through the non-uniform excitation field. The faster the object moves, the greater the rate of change of field, hence higher currents are induced at higher speeds. In the case of an AC illuminating field, the object experiences a changing field whether it is moving or not. Usually the AC frequency is chosen such that it is much greater than the blade passing rate. In this case, the effect of movement through the magnetic field is small compared to the effect of the changing current in the coil.
The eddy currents flowing in the conducting object result in a secondary field. The secondary field can be sensed by various means, such as by detecting the induced voltage in a separate detection coil, or by sensing impedance or apparent inductance changes in the drive coil.
Factors other than distance also affect the measurements made, important ones being changes in the temperatures of the sensor and the blades. Gas turbine engines are subject to extremes of temperature, both hot and cold. It is important therefore to reduce the effects of temperature on the distance reading in such applications. U.S. Pat. No. 4,716,366 describes a differencing technique to reduce temperature effects. U.S. Pat. No. 4,893,079 describes method where temperature correction is performed based on measurements of coil resistance. U.S. Pat. Nos. 7,324,908 and 4,970,670 describes methods wherein a temperature sensor is used to generate a correction factor. U.S. Pat. Nos. 6,479,990, 7,162,384 and 7,323,868 describe alternative techniques.
Other types of proximity sensor exist. These include capacitive types, such as those produced by Tyco Thermal Controls LLC, which measure the capacitance between a sensor and the object being measured, and use the principle that the capacitance will vary inversely with the separation distance.
Yet other types of proximity sensor employ antennas operating at high frequencies to create an electromagnetic field, and measure a disturbance to the field caused by an object entering the field. Such a sensor is described in International patent application, “Rotor Blade Sensor”, publication No. WO2009/034305.
According to a first aspect of the present invention there is provided a proximity sensor system comprising a transducer for detecting proximity to an object, means for receiving a signal from the transducer, and means for processing the received signal, the processing means being adapted to extract, from the received signal, a signal having a form related to the proximity of the object to the transducer, wherein the processing means is arranged to fit a curve to the extracted signal, the curve being chosen to approximate to a form of the signal, and to extract parameters pertaining to a width measurement and a height measurement from the fitted curve, with the parameters providing an indication as to the proximity distance of the object to the transducer.
By fitting a curve to the measured data, and taking measurements from the curve it has been found that a more consistent measurement can be obtained, that still provides the required accuracy of measurement of distance, but which is less susceptible to noise.
Also, under many circumstances (such as when measurements of high speed objects are made), when the received signal is digitised, the number of samples making up the extracted signal is limited, and hence measuring parameters directly from it may lead to significant quantisation errors in both temporal (or equivalently angular) measurements, and in any amplitude measurements made of the signal. The present invention, by fitting a curve to the extracted signal (and extracting desired temporal and/or amplitude parameters from the fitted curve) significantly reduces the effect of such quantisation errors.
The means for processing the received signal may also be arranged to carry out an averaging process whereby two or more extracted signals are averaged together in some manner. This may be done by combining pulse data before the curve fitting is carried out, the curve fitting then being done on combined data. The pulses being combined, comprising two or more pulses, may be time aligned by, for example, subtracting an appropriate time offset from the pulses. A curve fitting process as described herein may then be applied to the combined data. The pulses may relate to different, e.g. consecutive, objects passing the sensor, or alternatively the pulses may relate to multiple instances of the same object passing the sensor.
Averaging of data, such as extracted parameters, following the curve fitting on individual pulses may also be done, but overall requires more computation.
The transducer may optionally be an eddy current sensor. The eddy current sensor may be a single coil sensor, driven with an AC or DC signal, wherein impedance or inductance changes are detected as the coil is driven. Alternatively the eddy current sensor may have one or more coils configured as receive coils, with a DC magnetic field being produced with a permanent magnet. Alternatively, the eddy current sensor may have multiple coils, with one or more driven coils and one or more receive coils.
Alternatively, the transducer may be a capacitive sensor, such as that of the type mentioned above.
Alternatively, the transducer may be a radio frequency sensor that utilises an electromagnetic antenna, such as that of the type mentioned above.
When the transducer is an eddy current sensor, the drive signal may be an AC (alternating current) signal, or a DC (direct current) signal. If a DC current is used then the voltage level of returned signal will vary directly with the distance of the coil to the object, within an operational range, and will increase as the speed increases. If an AC signal is used then the envelope of the signal will be proportional to this distance within the operational range. The operational range is typically between 0.25 mm and 10 mm, with better performance generally being obtained with ranges up to around 7 mm.
If AC is used, then a signal may be extracted by demodulating the signal received from the coil. The demodulation may be a simple envelope detection, or may use a phase coherent approach. If the latter is used then the processing means is preferably arranged to receive a version of the signal sent to the coil, as well as signal received from it. This can act as a clock reference for the coherent demodulation process. The AC signal may be a sine wave. The frequency of the AC signal may be chosen to be compatible with the various components used in the system. Preferably the frequency is high enough to allow several, such as at least twenty, more preferably at least fifty, and still more preferably at least two hundred cycles to occur during the passing of each object of interest, such as a turbine blade. For good results the demodulated data should have enough bandwidth to accurately represent the change in signal as a function of position as the blade passes the sensor. From Nyquist theorem at least two samples per cycle of output bandwidth are required. In practice, because of the desire to filter out the strong signal at 2× the drive frequency which the demodulation process produces from the signal, but preserve the shape of the demodulated pulse, the operation frequency should be rather higher than this. As an example consider an engine with a blade frequency of 10 kHz. Suppose also that a low-pass −3 dB point of around 200 kHz is wanted in the demodulated data and that the 2× drive frequency must be attenuated by 80 dB using a 4th order Bessel filter. In this case a drive frequency of at least 1.5 MHz should be used.
The processing means may be a combination of analogue and digital circuitry. The analogue circuitry may comprise buffer, filters or demodulators, or may comprise merely the input of an analogue to digital converter (ADC). The digital processor may comprise a microprocessor, along with associated memory. The digital processor may have an output to a display, to another indicator (such as an audible alert) or to another system, to which the result of the distance measurement or any of the measured parameters may be passed.
The distance measurement may be compared against a threshold, and an output made from the processor if the separation distance is outside of a predetermined range when compared against the threshold.
The curve used to model the received signal is preferably chosen as being broadly similar in form to the extracted signal. The closer the similarity the better the performance will be in terms of providing an accurate measurement distance, but the skilled person will be aware that some curve forms may involve too much computation to be technically or financially practical on a given system, and so will choose an appropriate curve form based on these factors.
The process of choosing the curve to approximate to the form of the signal is advantageously done in a design or commission phase, wherein test measurements may be taken to ascertain the basic form of the returned signal, and a curve found that is then used in subsequent operations.
The curve may be a function derived from an empirical model of a system on which the present invention is to be used. The empirical model function may be based upon measured signals taken from a real system under test conditions. The function may be produced by, for example, using a look-up table, or by fitting a spline to the measured signals, or by using a Taylor series representation of the measured signal, or by any other suitable means.
The extracted signal may comprise a pulse waveform, with the pulses corresponding to the passing of objects past the transducer. In a turbine application the pulses may correspond to the turbine blade tip passing the transducer. Some turbine blades have a blade tip that is arranged to interface to the tips of the adjacent blades, so that a continuous shroud is formed at the edge. In this case the joints between the tips will still provide a discontinuity as far as generated eddy currents are concerned, and so will still provide a pulse output.
Each pulse corresponds to the passing of an object. It is therefore preferable to isolate each pulse and perform a curve fit on the single pulse. This will provide a measure relating to each individual object. The curve fitting process carried out by the processing means may be carried out using any standard method. For example, the curve may be fitted to the extracted signal using a least squares technique.
The parameters taken from the fitted curve may comprise for example the pulse height, the width at half height (or at any other desired fraction of the pulse height, such as the top or bottom 10%), and the position in time, or other suitable unit, of the pulse peak.
The measured parameters may be converted to a separation distance using, for example, a lookup table. This may be stored in the memory of the processor, and may be initially generated in a calibration step. Other methods of producing the look-up table may be used, such as by mathematical or computer modelling or other analysis. Multiple look up tables may be used, each associated with a different measured parameter.
The proximity distance values may be used as taken from the look-up table, or may instead be processed by, for example averaging the values to reduce the effects of random noise. Systematic errors, such as caused by temperature effects, may also be reduced by calculating correction factors from the averaged values, and applying the correction factor to instantaneous estimates or measurements taken.
The means for providing a drive signal to a coil in the transducer may comprise an oscillator, along with suitable amplification means. The oscillator and amplifier may be located at a suitable distance from the transducer to keep it clear of excess thermal or mechanical stresses, and may be connected to the transducer using wires. The means for receiving the signal from a coil in the transducer (which may be the same coil as is being driven by the driving means) may again comprise wiring from the transducer to the processing means. The driving means and the processing means may be conveniently co-located.
According to a second aspect of the invention there is provided a method for determining distance from a transducer to an object of interest, the transducer being a transducer for detecting proximity, comprising the steps of:
Where the curve has a width, such as, for example, a pulse signal, then the pulse width may advantageously be the parameter extracted. Other parameters that may be extracted include a peak height of the curve, and a position of the peak centre of the curve.
The curve used to approximate to the extracted signal may be derived from a known curve, such as those described above, or may be an empirically derived curve, formed from test data taken from measurements on a representative system. The measurements may be, for example, averaged data from many runs of the representative system, or may be taken from controlled, calibration runs.
The invention will now be described, by way of example only, with respect to the following Figures, of which:
The drive/tuning impedance (25) is in series with the transducer coil (23). Changes to the transducer coil (23) impedance or inductance caused by electrically conductive objects such as turbine blades coming into close proximity will affect the voltage or current (or both), and these changes can be detected by looking at the voltage across the drive/tuning impedance (25). Preferably the drive/tuning impedance (25) has very similar impedance characteristics to that of transducer coil (23). The impedance 25 may therefore be realised by using an identical component to that of transducer coil (23).
The signal at the output of the drive/tuning impedance, as well as feeding the transducer coil (23) via a coaxial cable (co-ax) (26), also feeds a screen (27) of the co-ax, via screen drive amplifier (34). The screen is driven so that the capacitance between the coil connection and the screen is effectively reduced to near zero. This is mainly to prevent this capacitance (which is in parallel with the sense coil) acting to lower the self-resonant frequency of the coil (23). To avoid the influence of this capacitance on the measurements the sensor should be operated well below this frequency. Capacitive effects are particularly undesirable in applications such as gas turbine engines because in the engine conductive carbon deposits can occur which have a large influence on the varying capacitance between the sense coil and the turbine blade. For this reason it is also desirable to have the ‘earthy’ end of the coil closest to the blade.
As said above, the voltage across the drive/tuning impedance (25) is used to detect the presence or otherwise of the object being in close proximity to the coil (23). The voltage at the output end of impedance (25) is applied to the positive input of a difference amplifier (28). The negative input of the difference amplifier (28) is fed by the voltage on the other side of the drive impedance (25), after being attenuated in attenuator (24). This subtracts off the input signal from the amplifier (21), offsetting it and improving the system dynamic range and increasing sensitivity. The attenuator (24) is adjustable and is used to provide some system adjustment to the amount subtracted. The output from the difference amplifier (28) is low pass filtered in filter (29) and then provided as a first input to demodulator (30).
The second path from the amplifier (21) feeds a second input to demodulator (30), this effectively acting as a phase reference signal. The demodulator (30) comprises a mixer circuit. The output of demodulator (30) is filtered in high, and low-pass filters (31,32) to remove unwanted frequency elements such as the phase reference signal, and any DC component, before being digitised in ADC (33). Subsequent processing is done digitally in processor (22).
The processing described in relation to
A person of ordinary skill in the art would realise that the circuitry discussed in relation to
The Anderson squared function shown is:
The dipole-dipole function is based on a point dipole representation of both the eddy current coil and the elongate elements of the test jig and is given by:
The Gaussian model is given by:
The Lorentzian is given by:
In all of the above, the parameter a is the peak height of the pulse, the parameter b is the width of the pulse, the parameter c is the time of the peak centre. The variable, t is time relative to an arbitrary start time. For turbine blade measurement applications the dimensions of parameters b and c may be converted to angular or distance units for convenience, given knowledge of the rotation speed of the turbine. In particular, the width parameter in dimensions of distance (e.g. around the circumference of a turbine) can be computed using the formula w=sb where s is the speed of the passing elements (e.g. turbine blade tips) and w is the pulse width in dimensions of distance (e.g. units of mm). The signal received from the sensor may be used to compute the speed. In a turbine application, given knowledge of the number of blades Nb on the rotor and of the diameter d (in metres) of the rotor, by measuring the frequency fb (in Hz) with which pulses occur, the speed, s (in metres/second), of the blade tips can be calculated using s=πdfb/Nb. Alternatively, measuring time T (in seconds) for n blades to pass, use s=πdn/(NT).
Of course, the clearance distance between the coil and the element is not related to the peak centre parameter, c; however, it is used as an indicator of which particular element (e.g. turbine blade) is currently being measured.
Data have been measured using simulated turbine blades on a test rig, using a single coil eddy current sensor. Calculations on the measured data showed that the Anderson squared function was found to produce the best fit, in terms of minimal error, defined as the mean squared difference between the measured data and the model. It can also be seen by eye that the Anderson squared function appears to most closely resemble the pulse.
Thus, for this example the Anderson squared function was chosen. The skilled person will understand that other applications, such as systems employing different coil types, different demodulation systems or different element shapes etc may find different functions, including ones not described above, or an empirically derived model function (as described above), have a better fit.
The pulse waveform shown in
An embodiment of the invention records the pulse data using circuitry described in relation to
b is a graph showing in a first plot the relationship between w, the width of the curve fitted to the pulse (shown with circular markers), and in a second plot the half-height width of the directly measured pulse (square markers), against clearance distance. The inset graph shows values at small clearance distances. The two plots are not identical because, owing to the form of the model function, the paramater b in the Anderson squared model represents the pulse width at 0.512 of the pulse amplitude. It has been found that this measurement (pulse width) is much less dependent upon the temperature of the coil and element, and so is useful at the high temperatures frequently encountered in turbines. As clearance distance increases the pulse signal becomes smaller and inevitably more noisy. This can be seen in
Once values for a and w have been obtained, they can be used to estimate the clearance between sensor coil and the passing element. This may be done for example either by reference to a predetermined calibration curve, or by interpolation from a look-up table. These may be themselves produced by, for example, making measurements in controlled conditions where the proximity distance is known, and recording several different measurements at different proximity distances.
The shape of the signal pulses depends on many factors. In addition to the clearance between sensor coil and the passing element, these factors include (but are not limited to) the temperature of the blade and the temperature of the sensor. The amplitude parameter, a, is strongly dependent on clearance, but is also significantly affected by blade temperature, and by resistance changes (which may be thermally induced) in the coil. Consequently thermal variations in the engine can introduce errors in the clearance distance estimated by using amplitude parameter, a, alone. In contrast, the width parameter, w, is only weakly affected by coil and blade temperature. It is more selectively dependent on the clearance. However the dependence is weaker, so random noise affects an estimate of clearance distance estimated using width parameter, w, more significantly than when using amplitude parameter, a.
Another embodiment makes use of the characteristic of the measurement of w being much less affected by changes in coil resistance (but generally having a poorer signal to noise ratio), and the characteristic of the measurement of a generally having a good signal to noise ratio (but relatively poor resilience to changes in coil resistivity, and hence temperature), as shown below.
The technique uses an instantaneous value of a to derive an estimate of the clearance, and then corrects that estimate using average values of w and a.
In more detail, curve fitting is used as described herein, and values for a and b for each pulse, i.e. each passing of the element being detected, are extracted from the fitted curves. Width parameter, b, is converted to dimensions of distance and subsequently denoted w. A look-up table is then used to produce a first estimate of clearance from the value of a, using e.g. a graph similar to that of
Given a set N of instantaneous pulse amplitude and width measurements ai and wi i.e. measurements derived from each individual passing element, average values of each are produced:
The values ā and
Suppose temperature effects cause ĥa, the clearance estimated using parameter a only, to be incorrect so that
ĥa=khactual
where k is an as yet unknown factor and hactual is the actual clearance. The same temperature effects cause negligible error on ĥw, the clearance estimated using parameter w only. However, ĥw is subject to random noise with standard deviation σhw. If an estimate of clearance ĥ
The average of k is then approximated by
The values of ĥā and ĥ
The value for N, i.e. the number of pulses over which the data is averaged, is preferably chosen to span a time that is short enough to enable the averages to track temperature variations of the objects being detected and of the sensing coil. N should also be large enough to ensure that the noise in ĥ
The examples above all use an eddy current sensor as the transducer. It will be appreciated that any form of sensor that produces a signal having characteristics that change according to the proximity between it and an object being measured, and which produces an output that may be modelled either with an analytic function or an empirically derived model, may be used. The examples also show a pulse signal as being that modelled. Other transducers, or other methods of processing the transducer outputs, may present a signal of a form different from a pulse as shown. It may be, for example, a bipolar pulse, a phase change, or a transition from one steady state level to another, such as a rising or falling edge. The present invention may be used with any such signal that may be modelled either with an analytic function or an empirically derived model. Of course, with signal shapes other than pulses, then different signal characteristics, such as rise time, fall time etc. may be those that are taken from any model to measure the proximity data.
The above examples have been disclosed for illustrative purposes, and those skilled in the art will appreciate that various modifications, additions and substitutions are possible, without departing from the scope of the invention as disclosed in the accompanying claims.
Number | Date | Country | Kind |
---|---|---|---|
1102542.6 | Feb 2011 | GB | national |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/GB12/00142 | 2/13/2012 | WO | 00 | 7/29/2013 |