【机械类毕业论文中英文对照文献翻译】近似熵作为诊断工具为机器健康监测
收藏
资源目录
压缩包内文档预览：(预览前20页/共24页)
编号：77693893
类型：共享资源
大小：1.35MB
格式：RAR
上传时间：2020-05-07
上传人：柒哥
认证信息
个人认证
杨**（实名认证）
湖南
IP属地：湖南
6
积分
- 关 键 词：
-
机械类毕业论文中英文对照文献翻译
机械类
毕业论文
中英文
对照
文献
翻译
近似
作为
诊断
工具
机器
健康
监测
- 资源描述：
-
【机械类毕业论文中英文对照文献翻译】近似熵作为诊断工具为机器健康监测,机械类毕业论文中英文对照文献翻译,机械类,毕业论文,中英文,对照,文献,翻译,近似,作为,诊断,工具,机器,健康,监测
- 内容简介：
-
Mechanical SystemsandSignal ProcessingMechanical Systems and Signal Processing 21 (2007) 824839Approximate Entropy as a diagnostic tool for machinehealth monitoringRuqiang Yan, Robert X. Gao?Department of Mechanical and Industrial Engineering, University of Massachusetts, Amherst, MA 01003, USAReceived 7 April 2005; received in revised form 16 February 2006; accepted 22 February 2006Available online 4 April 2006AbstractThis paper presents a new approach to machine health monitoring based on the Approximate Entropy (ApEn), which is astatistical measure that quantifies the regularity of a time series, such as vibration signals measured from an electricalmotor or a rolling bearing. As the working condition of a machine system deteriorates due to the initiation and/orprogression of structural defects, the number of frequency components contained in the vibration signal will increase,resulting in a decrease in its regularity and an increase in its corresponding ApEn value. After introducing the theoreticalframework, numerical simulation of an analytic signal is presented that establishes a quantitative relationship between theseverity of signal degradation and the ApEn values. The results of the simulation are then verified experimentally, throughvibration measurement on a realistic bearing test bed. The study has shown that ApEn can effectively characterise theseverity of structural defect, with good computational efficiency and high robustness.r 2006 Elsevier Ltd. All rights reserved.Keywords: Approximate Entropy; Machine health monitoring; Computational efficiency1. IntroductionThe past decades have seen increasing attention from the research community worldwide on machinecondition monitoring and health diagnosis. The efforts have fuelled the continued advancement of sensing aswell as signal-processing technologies. In addition to commonly used time (statistical 14) and frequency(spectral 2,5,6) domain techniques, advanced signal processing techniques, such as blind source separation,and wavelet transforms, have been investigated as new tools for fault detection in mechanical systems 710.Due to instantaneous variations in friction, damping, or loading conditions, machine systems are oftencharacterised by non-linear behaviours. Therefore, techniques for non-linear parameter estimation provide agood alternative to extracting defect-related features hidden in the measured signals that may not beeffectively identified using other techniques.A number of non-linear parameter identification techniques have been investigated, such as CorrelationDimension 11 and Complexity Measure 12,13, for patient sleep evaluation 14, gearbox tooth defectARTICLE IN PRESS/locate/jnlabr/ymssp0888-3270/$-see front matter r 2006 Elsevier Ltd. All rights reserved.doi:10.1016/j.ymssp.2006.02.009?Corresponding author. Tel.: +14135450868; fax: +14135451027.E-mail address: gao (R.X. Gao).diagnosis 15, and rolling bearing defect detection 16,17. Reliable estimation of the Correlation Dimension ofa time series requires that a large quantity of data points be calculated 18,19, leading to lengthycomputational time and insuitability for on-line, real-time applications. The Complexity Measure, incomparison, is computationally more efficient. The complexity of a signal can be described by two measures:the LempelZiv Complexity 12,20 and Approximate Entropy (ApEn) 13,21. The LempelZiv Complexityindicates the degree of regularity of a time series, and has shown to be effective in quantifying the relationshipbetween brain activity patterns and depth of anesthesia (DOA) in patients 22. As an example, theelectroencephalographs (EEG) of test patients have shown a higher Complexity value in the awake state thanin the sleep state. Furthermore, studies on the hearts conducting system have shown that important indicatorsof the heart rhythm regularity, such as Sinus rhythm (SR), ventricular tachycardia (VT), and ventricularfibrillation (VF) 23, could be effectively detected and separated through Complexity thresholding, whereasother techniques, such as Fourier transform-based VF-filtering, has shown difficulty. Analysis of Complexityvalues resulting from vibration measurement on a large rotating machine 24 has shown that the inceptionand growth of faults in the machine could be correlated with changes in the Complexity value. Another studyon defects detection in rolling bearings has established quantitative relationship between the Complexity valueand the defect size, and consequently the defect severity level 25.Evaluation of the Complexity values of a time-series is based on a coarse-graining operation by which thedata sequence is transformed into a pattern of only a few symbols, e.g. 0 and 1, and involves data sequencecomparison and number counting in only one dimension 20. In comparison, Approximate Entropyexpresses the regularity of a time series in multiple dimensions, and contains more time-related information21. This makes it an attractive tool for monitoring system dynamics, as information on the temporalprogression of a defect is valuable not only for diagnosing the current machine condition, but also inpredicting its future behaviour. In a study on turbo generator, it was found that the ApEn values haveincreased significantly after looseness of the bearing bushing was identified, when compared with normaloperation conditions 26. In the field of biomedical engineering, studies on foetal heart rates usingcardiotocography (CTG) have shown correlation between the increase of ApEn values and pathologicalconditions 27. Similar studies on EEG signals have found that an increase of the ApEn values can be used todetect and characterise eplieptic activities 28. In another study, different sleep states were successfullyidentified by variations of the ApEn values of EEG signals 29. In the field of physiological data analysisfor plants, ApEn values were found to be clearly indicative of changes in the sap flow dynamics of twotropical tree species under water deficit 30. In addition, ApEn has been applied to the detection of humanmotion, in which static background, motion envents, and temporal clutter could be distinguished from eachother, based on their ApEn values 31. Furthermore, ApEn has shown to be effective as a feature detector, foridentifying potential areas of interest within an image for target detection and image compression 32,33.Recently, the application of ApEn has been extended to assess subtle, exploitable changes in financial markettime series 34.This paper investigates the utility of ApEn as a measure for characterising machine operationconditions, with the goal to introduce a quantitative diagnostic tool for machine health monitoring. Afterpresenting the theoretical background, the effects of several parameters that are of key relevance to the ApEncalculation, e.g. data length, sampling rate, dimension, and tolerance value, are discussed. Numericalsimulation of an analytic signal is then presented to establish a quantitative relationship between the signalregularity (which is affected by the degradation of structural health due to defect propagation) and the ApEnvalues. Simulation results are then experimentally verified, using vibration signals measured on a realisticbearing test bed. Good agreement between the simulation and experimentation was found, confirming theeffectiveness of ApEn in characterising the severity of structural defect, and its computational efficiency andhigh robustness.2. Theoretical backgroundFor a time series S containing N data points x(1), x(2),y,x(N), its regularity can be measuredby ApEn in a multiple dimensional space 13, in which a series of vectors are constructed andARTICLE IN PRESSR. Yan, R.X. Gao / Mechanical Systems and Signal Processing 21 (2007) 824839825expressed asX1 fx1;x2;.;xmgX2 fx2;x3;.;xm 1g.XN ? m 1 fxN ? m 1;xN ? m 2;.;xNg:(1)In Eq. (1), each of the vectors is composed of m consecutive and discrete data points of the time series S13,21. The distance dXi; Xj between two vectors X(i) and X(j) can be defined as the maximum differencein their respective corresponding elements:dXi;Xj maxk1;2;:;mjxi k ? 1 ? xj k ? 1j,(2)where i 1;2;.;N ? m 1, j 1;2;.;N ? m 1, and N is the number of data points contained in thetime series.For each vector X(i), a measure that describes the similarity between the vector X(i) and all other vectorsXj; j 1;2;.;N ? m 1; jai can be constructed asCmir 1N ? m ? 1XjaiYfr ? dXi;Xj?g;j 1;2;.;N ? m 1,(3)whereYfxg 1;xX0;0;xo0:(4)The symbol r in Eq. (3) represents a predetermined tolerance value, defined asr k ? stdS,(5)where k is a constant (k40), and stdd represents the standard deviation of the time series. By definingfmr 1N ? m 1XilnCmir?;i 1;2;.;N ? m 1(6)the ApEn value of the time series can be calculated asApEnm;r limN!1fmr ? fm1r?.(7)For practical applications, a finite time series consisting N data points is used to estimate the ApEn value ofthe time series, which is defined asApEnm;r;N fmr ? fm1r.(8)Eq. (8) indicates the similarity among the reconstructed vectors within the time series, when the dimensionof the vectors has increased from m to m+1 21. This affects the regularity of the time series being analysedand consequently, the associated ApEn value. The greater the regularity is, the lower the ApEn value. As anexample, a periodic time series containing only a single frequency component will have a low ApEn value(which approaches zero), due to the high regularity of the signal. In contrast, complex time series containingmultiple frequency components will have high ApEn values, due to a low level of regularity.To calculate the ApEn value of a given time series S, data points contained within the series are firstrearranged into a series of m and an m+1 dimensional vectors, respectively, as illustrated in Fig. 1. Thedistances between two corresponding data points from each vector are then calculated. Based on Eq. (3), thesimilarity measure for the reconstructed vectors under a given r value is obtained, for each dimension m andm+1. Subsequently, the ApEn value for the time series S with N data points is calculated using Eq. (8). Toensure consistency of the ApEn calculation, a minimum data length N, as well as appropriate dimension m andtolerance r need to be determined. In the presented study, the effects of data length, dimension, and toleranceon the ApEn value calculation have been investigated, based on realistic vibration signals measured from abearing test bed.ARTICLE IN PRESSR. Yan, R.X. Gao / Mechanical Systems and Signal Processing 21 (2007) 8248398262.1. Effect of data lengthFig. 2 illustrates a vibration signal measured from an undamaged, healthy bearing under a radial load of9kN and a sampling rate of 10kHz. The power spectrum density of the signal indicates the existence of fourmajor frequency components: (1) bearing unbalance frequency fu, which occurs when the gravitational centreof the bearing does not coincide with its rotational center, as illustrated in Fig. 3(a); (2) misalignmentARTICLE IN PRESSFig. 1. Major steps in the Approximate Entropy calculation.2000300040005000010006000-0.4-0.200.20.40.6-0.6Sampling Data PointsSignal (V)05010015020000.020.040.060.080.10FrequencyPSD (Watts/Hz)fufmfBPFO2% fBPFOFig. 2. Vibration signal from an undamaged, healthy bearing.R. Yan, R.X. Gao / Mechanical Systems and Signal Processing 21 (2007) 824839827frequency fm, which occurs when the two raceways of the bearing (inner and outer) fall out of the same plane,resulting in a raceway axis that is no longer parallel to the axis of the rotating shaft, as illustrated in Fig. 3(b);(3) ball passing frequency in the outer raceway fBPFO, which is related to the periodic passing of the rollingballs over a fixed reference position; and (4) its harmonic, with the numerical value determined by (2?fBPFO).The relationship between the ApEn value and the data length taken for the calculation is illustrated in Fig. 4,where six vibration signals measured on the same bearing are comparatively displayed, under a sampling rateof 6, 10, 14, 18, 22 and 26kHz, respectively. It is seen that when the data length is greater than 1000 datapoints, variation of the ApEn values with respect to each specific sampling rate becomes insignificant.Furthermore, the difference between the ApEn values caused by the different sampling rates approaches aconstant and small value. For example, the difference of ApEn values between 6 and 26kHz is 0.104; as thesampling rate increases to over 14kHz, the difference quickly diminishes to below 0.01. Such result isconsistent with the theoretical analysis presented in 13,21, where a data length between 750 and 5000 datapoints was considered sufficient for achieving stable and consistent ApEn values.Once the data length is fixed, a higher sampling rate translates into a shorter data sampling period, and viceversa. A shorter sampling period means a shorter observation span of the time series, and subsequently, aARTICLE IN PRESSFig. 3. Models of unbalance and misalignment of a rolling bearing.Number of Data Points0100020003000400050006000ApEn0.050.100.150.200.250.30fs=6 kHzfs=10 kHzfs=14 kHzfs=18 kHzfs=22 kHzfs=26 kHz0 SamplingRate (kHz)610141822260.1420.0840.0540.0450.0410.038ApEnDifference0.0580.0300.0090.0040.003Fig. 4. Effect of data length on the ApEn values.R. Yan, R.X. Gao / Mechanical Systems and Signal Processing 21 (2007) 824839828reduced amount of temporal information to be included for the ApEn calculation, leading to a lower ApEnvalue. Such a trend is noted in Fig. 4, where the highest sampling rate fs 26kHz has resulted in the lowestApEn value (e.g. ApEn 0.038 for 1000 data points). To further quantify the effect of data sampling rate,several vibration signals from a healthy bearing were measured under different sampling rates. The ApEnvalues calculated from these signals are shown in Fig. 5, in which the dimension m (m 2, 3,y,10) serves as aparameter. It is seen that the ApEn value decreases with an increasing sampling rate. When the sampling ratehas increased beyond a certain threshold (e.g. 22kHz), the ApEn values approach a steady-state value (around0.04), with the fluctuations less than 4%.2.2. Effect of toleranceAnother parameter affecting the ApEn value calculation is the tolerance r. As defined in Eq. (5), selection ofr depends on the k value for a specific time series S that contains a fixed number of data points. The effect of kon the ApEn value calculation is shown in Fig. 6, where nine different dimensions m were evaluated to verifythe consistency of the ApEn calculation results.ARTICLE IN PRESS812162000.4000.100.200.30Sampling Rate (KHz)ApEn2824m = 2m = 3m = 4m = 5m = 6m = 7m = 8m = 9m = 104Fig. 5. Effect of sampling rate on the ApEn values.ApEnm = 2m = 3m = 4m = 5m = 6m = 7m = 8m = 9m = 10k Value00.250.300.350.400.050.100.150.200.2500.450.50Fig. 6. Effect of k on the ApEn values.R. Yan, R.X. Gao / Mechanical Systems and Signal Processing 21 (2007) 824839829As the k value increases, variability in ApEn decreases continually for the corresponding dimension m, atdifferent rates. The difference in ApEn associated with the different m values becomes negligibly small, as kapproaches 0.4. For example, at k 0:4, variation of the ApEn due to the different m values is less than 4%.Such a variation was considered acceptable for the presented study, thus the tolerance value r correspondingk 0:4 was chosen for the subsequent ApEn calculations.2.3. Effect of dimensionAs indicated in Eq. (1), increasing the dimension m will lead to more temporal information on the signal tobe included in the ApEn calculation. While improving the overall calculation accuracy, a higher m value leadsto higher computational cost. This is confirmed experimentally in Fig. 7, where an increase of the dimension mfrom 2 to 10 has resulted in an increase of over three times for the computation time (from 0.5 to 2.05s). Theevaluation was performed for k40:4, given that the ApEn values did not change appreciably beyond thisvalue, as shown in Fig. 6. Since the corresponding ApEn values changed less than 6.6% when m increased from2 to 10, a low m value is generally preferred for the ApEn calculation. Thus, m 2 was chosen for thepresented study.3. Formulation of the test signalTo quantitatively characterise a machines health state and its degradation (as the severity of the structuraldefects increases), a test signal is needed to provide a reference base. Such a test signal was formulated for thepresented study, based on vibration signals measured from a rolling bearing, as illustrated in Fig. 2. Throughan inverse Fourier transform, four major frequency components (out of a total of 213 8192 frequencycomponents) identified in the bearing signal spectrum (at 10, 20, 41, and 82Hz, respectively) were combinedwith an additive noise e(t) to construct the test signal S(t). Together, the four components represented over91% of the energy content of the original signal. The test signal is expressed asSt Xt et.(9)In Eq. (9), Xt represents the reconstructed signal from the four frequency components of the originalbearing signal, and is defined byXt 0:016 sin2p10t 0:112 sin2p20t 0:325 sin2p41t 0:029 sin2p82t.(10)ARTICLE IN PRESS35812.500.51.01.5m ValueComputation Time(s)1092.0611247mTime(s)23456789100.500.681.001.191.352.051.841.551.68Fig. 7. Dimension m vs. computation time t.R. Yan, R.X. Gao / Mechanical Systems and Signal Processing 21 (2007) 824839830The noise component e(t) is defined by the signal-to-noise ratio (SNR) 35 asSNR 10log10Xt?2=Met?2=L,(11)whereXt?2PXt2,et?2Pet2, and M and L are the length of X(t) and e(t), respectively.The constructed test signal is shown side-by-side with the measured bearing vibration signal in Fig. 8. Usingthis test signal, the ApEn values corresponding to different SNRs were calculated, as listed in Table 1.It is noted that the ApEn value of the test signal has increased by over eight times when the SNR decreasedfrom 100 to 0dB. The decrease of the SNR is reflected in the waveform becoming increasingly noisier. Thisis analogous to the deterioration of a machine system where defects initiate and propagate through thestructure. Furthermore, since the noise signal was added to the test signal with an increasing bandwidth, as theSNR ratio was decreased, it corresponds to an increasing number of frequency components contained inthe signal, as the column Spectrum in Table 1 indicates. This simulation result confirms that the ApEn valueprovides a quantitative measure for characterising the degradation of a dynamic signal, which can berepresented by the deterioration of a machines health condition.ARTICLE IN PRESS2000300040005000010006000-0.4-0.200.20.40.6-0.6Sampling Data PointsSignal (V)-0.4-0.200.20.40.6-0.6Signal (V)2000300040005000010006000Sampling Data Points(a)(b)Fig. 8. Test signal formulated as a reference base for quantitative ApEn evaluation.Table 1ApEn value of test signal under various SNRR. Yan, R.X. Gao / Mechanical Systems and Signal Processing 21 (2007) 8248398314. Experimental verificationTo investigate the suitability of Approximate Entropy for machine health monitoring and diagnosis,experimental studies were conducted on a custom-designed bearing test bed, as shown in Fig. 9. Vibrationsignals from five ball bearings (i.e. four NU205 bearings from SKF 36 and one 1100 KR bearing fromTimken 37) with pre-seeded structural defects were measured under a range of rotational speeds and loadingconditions. Since the vibration data contain physical information on the working condition (or health state) ofthe bearing being monitored, evaluation of their corresponding ApEn values would reveal the physicalcondition of the bearing. The bearing test system consisted of a DC motor (Lesson Cl6D34FT18), twosupporting pillow blocks (SKF 209-112), the test bearing with supporting housing, a hydraulic cylinder (Miller4Z644B), a hydraulic pump (Enerpac P392), and an optical encoder for rotational speed measurement. Thehydraulic cylinder was able to apply a wide range of loads to the bearing in the radial direction. Motor speedwas adjusted by using an integrated speed controller. An accelerometer (PCB 621B40) was installed on the testbearing housing for vibration measurement. To differentiate the effect of various parameters on ApEncalculation, experiments were conducted using bearings with different structural defect severities, underdifferent rotational speeds and loading conditions, as described in the following.4.1. Effect of defect severityShown in Fig. 10 are four vibration signals acquired from four test bearings (SKF NU205 36), under fourdifferent defect severity conditions: healthy (undamaged), light damage (0.1mm hole), medium damage(0.5mm hole), and severe damage (1mm hole). The signals were obtained under a rotational speed of 600rpm anda hydraulic pipeline pressure of 500 PSI, which corresponds to a radial load of 9.163kN applied to the test bearing.The difference among the four bearing conditions is not readily recognisable from the time domain displays.However, the corresponding ApEn values calculated have revealed a consistent increasing trend, as the bearinghealth condition deteriorates. Table 2 illustrates the ApEn values for each of the four bearing defect severityconditions, calculated out of the average of five repetitive test runs. The ApEn values are stable for all the fourconditions, as indicated by the associated standard deviations.ARTICLE IN PRESSFig. 9. Experiment test bed configuration.R. Yan, R.X. Gao / Mechanical Systems and Signal Processing 21 (2007) 824839832As shown in Table 2, the ApEn values of the defective bearings are consistently higher than that of thehealthy bearing. The larger the defect size is, the higher the ApEn value. This can be explained by the fact thatvibration signals from a defective bearing contain more frequency components, as indicated in Table 1,leading to a lower regularity of the signal pattern and subsequently, a higher ApEn value. Realistically, anynon-ideal operation conditions would introduce additional frequency components to the bearing spectrum,which in turn increases its ApEn values. Examples of such non-ideal conditions include unbalance, asillustrated in Fig. 3(a), or misalignment, as illustrated in Fig. 3(b). As the defect size increases, defect-inducedvibration would become stronger than that caused by a smaller defect, and more impacts would be generatedwhen the rolling elements roll into and out of the defect, leading to a longer period of defect-related vibration.As a result, more frequency components will be contained within the frequency spectrum, and the ApEn valuewill become higher.To verify the robustness of the ApEn as a measure for defect severity classification, further tests wereconducted for the same set of test bearings, under different loads (e.g. 7.331kN or 400 PSI) and rotationalspeeds (e.g. 1200rpm). The results listed in Table 3 indicate a similar trend of the ApEn value distribution tothose listed in Table 2.ARTICLE IN PRESSTable 2Approximate Entropy values for different bearing defect severity levels: test 1 (radial load: 9.163kN; speed: 600rpm)Approximate entropyHealthy bearingDefective bearing (defect size: mm)ApEnAvg.0.0830.3160.4280.512Std. dev.0.0300.0300.0400.030-0.4-0.200.20.4020406080100Time (ms)Signal (V)-0.4-0.200.20.4020406080100Time (ms)Signal (V)-0.4-0.200.20.4020406080100Time (ms)Signal (V)-0.4-0.200.20.4020406080100Time (ms)Signal (V)(a)(b)(c)(d)Fig. 10. Bearing vibration signals for four different health conditions.Table 3Approximate Entropy values for different bearing defect severity levels: test 2 (radial load: 7.331kN; speed: 1200rpm)Approximate entropyHealthy bearingDefective bearing (defect size: mm)ApEnAvg.0.2250.4810.5870.676Std. dev.0.0410.0390.0310.032R. Yan, R.X. Gao / Mechanical Systems and Signal Processing 21 (2007) 824839833Fig. 11 illustrates the increase of the ApEn value as the defect size grows, where the error bars indicate thevariability of the ApEn value calculations. It is evident that ApEn provides a quantitative measure forclassifying the defect severity levels.4.2. Effect of rotational speedThe effect of bearing rotational speed on the ApEn values was investigated, with the mean value andstandard deviations for each operating condition given in Table 4 and Fig. 12. When the speed increases, theApEn values increase as well, for both the healthy and defective bearings. This can be explained by the factthat the frequency band of the bearing vibration signal broadens as the rotational speed increases. Also, therate of new vibration pattern generation will increase, resulting in an increase of the ApEn value. Under thesame speed, vibration from a defective bearing contains more frequency components due to periodicinteractions between the rolling balls and the structural defect, thus resulting in a higher ApEn value. Thisexplains why the curve of the ApEn value for the defective bearing lies above the curve for the healthy bearing.4.3. Effect of loadThe effect of loading on the ApEn values of the bearing is illustrated in Table 5 and Fig. 13, respectively.Radial loads were supplied to the test bearing from a hydraulic pump system on the bearing test bed. Theresults indicate that ApEn values did not change appreciably within the range of loads tested, for either thehealthy or defective bearings. This is because load changes do not affect the frequency composition of thevibration signals, thus no ApEn variations will incur.ARTICLE IN PRESSACBA) healthy - new bearing;B) light defect - one 0.1 mm diameter hole;C) medium defect - one 0.5 mm diameter hole;D) severe defect - one 1 mm diameter hole.Defect SizeApEn0.150.300.450.600.9000.759.163 KN, 600 rpm7.331 KN, 1200 rpmDFig. 11. Approximate Entropy values for different defect sizes.Table 4Approximate Entropy under different rotational speeds (Radial load: 9163N)Speed (rpm)Healthy0.1mm defectAvg.Std. dev.Avg.Std. dev.3000.0440.0100.2250.0306000.0830.0100.3160.0209000.1430.0200.4220.02012000.2260.0400.4820.040R. Yan, R.X. Gao / Mechanical Systems and Signal Processing 21 (2007) 8248398344.4. Life cycle testTo evaluate the suitability of the ApEn value as a measure for machine health indicator, a run-to-failureexperiment was conducted on a deep groove ball bearing of 52mm outer diameter (type 1100 KR) 37, underARTICLE IN PRESSSpeed (rpm)030060090012001500ApEn0.10.20.30.40.50.60Light defect - one 0.1 mm diameter holeHealthy - new bearingFig. 12. Effect of rotational speed on the ApEn values (radial load: 9163N).Table 5Approximate Entropy under different radial loads (Speed: 600rpm)Radial load (kN)Healthy0.1mm defectAvg.Std. dev.Avg.Std. dev.3.6650.0830.0100.3200.0307.3310.0900.0300.3160.0209.1630.0830.0200.3160.02014.6610.0810.0100.3140.02018.3270.0810.0100.3130.020Radial Load (kN)05101520ApEn0.10.20.30.40.50.60Light defect - one 0.1 mm diameter holeHealthy - new bearingFig. 13. Effect of radial loads on the ApEn value (Speed: 600rpm).R. Yan, R.X. Gao / Mechanical Systems and Signal Processing 21 (2007) 824839835a radial load of 5498N. This experiment was continually run under a rotational speed of 2000rpm. Tosimulate a structural defect, a 0.27mm-wide groove was introduced across the outer raceway, as shown on theright side of Fig. 14. After running continually for about 2.34 million revolutions, the defect size was measuredto have enlarged to 5.5mm, as shown on the left side of Fig. 14. Upon reaching approximately 2.7 millionrevolutions, the defect had propagated throughout the entire raceway and rendered the bearing practicallynon-functional. Vibration signals were taken during the experiment at an interval of every 7min, representingthe process of defect propagation.For each vibration signal obtained from the life cycle test, 40 data segments of 1000 data points were chosento calculate the ApEn values. The mean values were used to represent the actual ApEn value of the respectivesignal segment. The relationship between the ApEn value and revolution is shown in Fig. 15, in which threedefect propagation phases (I through III) were specified. The size of the bearing defect was measured in eachof the phases. It can be seen that the ApEn value has shown a generally increasing trend as the generalcondition of the bearing deteriorates (from Phase I to Phase III). Such a trend is consistent with the trend ofthe vibration amplitude measured from the test bearing, as illustrated in Fig. 15.To enable performance comparison, the LempelZiv complexity 12,20 was also calculated for the signalsobtained during the bearing life cycle test. Fig. 16 illustrates comparatively the percentage increase of theLempelZiv complexity and ApEn values as a function of the bearing revolution. While both have shown anincreasing trend corresponding the deterioration of the bearings health status, the ApEn values have shown alarger percentage of increase than that of the LempelZiv complexity, for each phase transition. This isquantitatively illustrated in Table 6, and can be attributed to the fact that more time-related information isinherently contained within the ApEn calculation procedure than the complexity formulation. As a result,ARTICLE IN PRESSFig. 14. Test bearing with a 0.27mm outer raceway initial scratch.Vibration Amplitude (V)1.01.52.02.500.50.600.300.400.500.55Revolution (million)ApEn3.0Phase IPhase IIPhase IIITest Condition:Radial Load - 5.498 KNShaft Speed - 2000 rpm0.450.351.500.300.701.101.300.900.50ApEnVibration AmplitudeFig. 15. ApEn and vibration magnitude values at various phases of the bearing life cycle test (radial load: 5.498kN, speed: 2000rpm).R. Yan, R.X. Gao / Mechanical Systems and Signal Processing 21 (2007) 824839836ApEn presents an inherently more effective measure for differentiating the working status of a machine beingmonitored than the LempelZiv complexity.5. ConclusionsApproximate Entropy has shown to present an effective measure for indicating the health state of a machinesystem. Experimental studies on a bearing test bed have shown that deterioration of a test bearing under aconstant speed can be effectively identified through the increase in the ApEn values (e.g. 276% increase when0.1mm diameter hole was introduced into the outer raceway of an initially healthy, undamaged bearing). TheApEn values have also shown to be sensitive to increase of the machines speed. For a healthy bearing, theApEn value increased by about 513% when the bearing rotating speed increased from 300 to 1200rpm, whilefor a defective bearing, ApEn value increased by about 214%. The change of bearing load has shown noappreciable effect on the ApEn values. From a run-to-failure test, ApEn has shown to be an effective measurefor characterising the degradation of the spindle-bearing system.To investigate the computational efficiency of the ApEn algorithm, experiments were conducted on aLaptop computer with a 2.0GHz CPU and 512MB memory. On average, it took the algorithm 0.75s tocomplete a calculation for 1000 data points. Such a result indicates that ApEn is generally appropriate for on-line applications.Research is being continued to analyse vibration signals from different types of defects and on differenttypes of bearings, to further validate the broad applicability of this technique for machine system healthmonitoring and diagnosis. Furthermore, the integration of ApEn with recurrent neural networks is beinginvestigated to realise an effective machine remaining life prognosis system.ARTICLE IN PRESS00.51.01.52.02.53.001020304050ComplexityApEnPhase IPhase IIPhase IIIPercentage of Increase (%)Revolution (million)Fig. 16. Percentage of increase for ApEn and complexity value as a function of the bearing revolution (radial load: 5.498kN, speed:2000rpm).Table 6Approximate Entropy and complexity value at each phase of the bearing life cycle testPhaseDefect size (mm)ApEnComplexityAvg.Percentage of increaseAvg.Percentage of increaseI0.270.37411.20.50210.9II5.500.4160.557III12.000.56049.70.66632.7R. Yan, R.X. Gao / Mechanical Systems and Signal Processing 21 (2007) 824839837AcknowledgmentThe authors gratefully acknowledge funding provided to this research by the National Science Foundationunder award # DMI-0218161. Experimental support from the SKF and Timken companies are appreciated.References1 D. Dyer, R.M. Stewart, Detection of rolling element bearing damage by statistical vibration analysis, Journal of Mechanical Design100 (1978) 229235.2 N. Tandon, A comparison of some vibration parameters for the condition monitoring of rolling element bearings, Journal of theInternational Measurement Confederation 12 (3) (1994) 285289.3 H.R. Martin, F. Honarvar, Application of statistical moments to bearing failure detection, Applied Acoustics 44 (1995) 6777.4 F. Honarvar, H.R. Martin, New statistical moments for diagnostics of rolling element bearings, Journal of Manufacturing Scienceand Engineering 119 (1997) 425432.5 A. Sturm, S. Billhardt, Envelope curve analysis of machines with rolling-element bearings, in: Proceedings of IFAC/IMACSSymposium on Fault Detection, Supervision and Safety for Technical ProcessesSAFEPROCESS91, Baden-Baden, Germany,September 1013, 1991, pp. 529533.6 R. Jones, Enveloping for bearing analysis, Sound and Vibration 30 (1996) 1015.7 G. Gelle, J. Galy, G. Delaunay, Blind source separation: a tool for system monitoring and fault detection? in: Proceedings of IFAC/IMACS Symposium on Fault Detection, Supervision and Safety for Technical ProcessesSAFEPROCESS 2000, Budapest,Hungary, 1416 June, 2000, pp. 705710.8 S. Petropol, S. Lesecq, A. Barraud, Fault detection using on-line wavelet analysis: application induction motor, in: Proceedings ofIFAC/IMACS Symposium on Fault Detection, Supervision and Safety for Technical ProcessesSAFEPROCESS 2000, Budapest,Hungary, 1416 June, 2000, pp. 11251130.9 C. Wang, R. Gao, Wavelet transform with spectral post-processing for enhanced feature extraction, IEEE Transactions onInstrumentation and Measurement 52 (4) (2003) 12961301.10 Z.K. Peng, F.L. Chu, Application of the wavelet transform in machine condition monitoring and fault diagnostics: a review withbibliography, Mechanical Systems and Signal Processing 18 (2) (2004) 199221.11 D. Rulle, Deterministic chaos: the science and the fiction, Proceedings of the Royal Society of London A 427 (1990) 241248.12 A. Lempel, J. Ziv, On the complexity of finite sequences, IEEE Transactions on Information Theory 22 (1) (1976) 7581.13 S.M. Pincus, Approximate Entropy as a measure of system complexity, Proceedings of the National Academy of Sciences of theUnited States of America 88 (6) (1991) 22972301.14 T. Kobayashi, S. Madokoro, Y. Wada, Human sleep EEG analysis u
- 温馨提示:
1: 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
人人文库网所有资源均是用户自行上传分享，仅供网友学习交流，未经上传用户书面授权，请勿作他用。
川公网安备: 51019002004831号