WO2004013589A1 - 振動検出方法及びその装置 - Google Patents

振動検出方法及びその装置 Download PDF

Info

Publication number
WO2004013589A1
WO2004013589A1 PCT/JP2002/007949 JP0207949W WO2004013589A1 WO 2004013589 A1 WO2004013589 A1 WO 2004013589A1 JP 0207949 W JP0207949 W JP 0207949W WO 2004013589 A1 WO2004013589 A1 WO 2004013589A1
Authority
WO
WIPO (PCT)
Prior art keywords
vibration
threshold
value
data
time
Prior art date
Application number
PCT/JP2002/007949
Other languages
English (en)
French (fr)
Inventor
Yoichi Ageishi
Hidetoshi Yamaguchi
Keisuke Suzuki
Original Assignee
Synchro Co., Ltd.
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Synchro Co., Ltd. filed Critical Synchro Co., Ltd.
Priority to PCT/JP2002/007949 priority Critical patent/WO2004013589A1/ja
Priority to AU2002323803A priority patent/AU2002323803A1/en
Publication of WO2004013589A1 publication Critical patent/WO2004013589A1/ja

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/04Bearings
    • G01M13/045Acoustic or vibration analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H1/00Measuring characteristics of vibrations in solids by using direct conduction to the detector
    • G01H1/003Measuring characteristics of vibrations in solids by using direct conduction to the detector of rotating machines
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/02Gearings; Transmission mechanisms
    • G01M13/028Acoustic or vibration analysis

Definitions

  • the present invention relates to an inspection method and an inspection device for inspecting whether an abnormal portion such as a void or a crack exists inside a concrete structure.
  • a mechanism is used to repeatedly apply a mechanical impact to the wall of a structure while moving in a tunnel, and analyze the impact sound of the impact to identify the position and size of the abnormal part at high speed.
  • the present invention relates to a suitable vibration inspection method and apparatus.
  • the present invention can be applied to an inspection method and an apparatus for determining whether or not there is an abnormality in a rotating shaft, a bearing, an impeller, or the like of a rotating device such as a pump motor.
  • the present invention relates to a vibration inspection method and apparatus suitable for analyzing the waveform of background noise and vibration generated by an object and specifying the type of abnormality with high accuracy. Background technology '''
  • the frequency analysis method converts the amplitude of the waveform into a spectral component using a method such as FFT (Fast Fourier Transform), and plots it on the frequency axis so that it can be observed as a component intensity distribution table.
  • FFT Fast Fourier Transform
  • a detection window for detecting a range in which the spectrum component intensity of the frequency changes is represented by four numerical values of an upper limit value and a lower limit value of the frequency, and an upper limit value and a lower limit value of the spectrum intensity.
  • the spectrum component of the new data is compared with the numerical value of the above detection window, and whether or not the spectrum component exists in the range is determined, and the abnormal state and the normal state are determined based on the identification result.
  • the time axis analysis method includes a method of comparing the effective value of the waveform, a method of comparing the peak value of the waveform, and a response waveform for one impact such as a hitting analysis using a strike.
  • a method of comparing the shapes of the two is known.
  • the measurement data contains a lot of noise
  • the measurement data is passed through a frequency filter to minimize the effect of noise, and the components in a specific frequency band are attenuated to reduce the noise.
  • a method of analyzing the waveform data is described. This method is effective when the frequency component of the noise is biased to a specific frequency range and is significantly different from the frequency range of the required signal frequency component.
  • a frequency filter that removes a specific frequency band is prepared.
  • the method of comparing the response waveform shape for one impact is a method in which reflected waves due to structural deformation are likely to appear and the effect can be faithfully reflected.
  • the noise component that changes randomly is averaged to reduce the average intensity, and the signal component appears repeatedly, so the fact that the intensity does not change even after averaging is used.
  • the vibration value after installation is significantly larger than the vibration value at the time of shipping inspection of the product to be inspected due to the difference in the installation location such as the rigidity of the floor on which the inspection target is installed is reduced May be strongly different.
  • the conventional method of comparing the peak values with one threshold value even a normal target may become abnormal. Therefore, it was necessary to observe for a sufficiently long time at the installation site to set a threshold with a margin for the peak value of the normal vibration waveform.
  • the measurement becomes difficult there is a problem that setting of the value becomes extremely difficult unless the initial deterioration is detected.
  • a method of comparing the shape of the response waveform for one impact is a method in which reflected waves due to structural deformation are likely to appear in a laboratory or the like, and the effect can be faithfully reflected.
  • the subject of concrete structures is located outside a noisy environment such as a tunnel or a bridge, and is susceptible to noise.
  • the method of analyzing each impact makes the structure surface dense. In order to perform detailed and efficient scrutiny, it is necessary to increase the number of hits per unit time, and to increase the hitting capacity, shorten the hitting interval, and at the same time, increase the computational capacity required for analysis. There o
  • a first object of the present invention is to provide a vibration inspection method and a vibration inspection method which have a high detection rate of abnormalities and abnormalities and are excellent in distinguishability from noise included in vibration data.
  • a second object of the present invention is to provide a vibration inspection method and an apparatus for easily confirming the results of the vibration inspection.
  • a third object of the present invention is to provide a vibration inspection method and apparatus capable of automatically judging the presence or absence of an abnormality or deformation.
  • a fourth object of the present invention is to provide a vibration inspection method and a vibration inspection method capable of easily setting a judgment criterion for automatically judging abnormalities or abnormalities.
  • a fifth object of the present invention is to provide a vibration inspection method and a vibration inspection method capable of giving a continuous impact with a simple mechanism in a sound inspection of a concrete structure or the like and determining a plurality of response waveforms at once. To be.
  • the above purpose is to set two or more threshold values in advance, measure the vibration data of the inspection object for a certain period of time, and convert the measured vibration data into digital waveforms. This is achieved by calculating the time that exceeds the threshold value, setting the time as the threshold excess time, and comparing the magnitude of this value to determine the degree of abnormality of the inspection object.
  • a value smaller than the maximum peak value of the vibration waveform data of a normal inspection object is set as the first threshold, and a second threshold is set to a value equal to or less than the above threshold. Then, by comparing the magnitude of the time exceeding the threshold, the degree of abnormality of the inspection object is determined.
  • the method is characterized in that matching is determined.
  • a value smaller than the maximum peak value of the vibration waveform data of a normal inspection object is set as the first threshold value, and a second threshold value is set to a value equal to or less than the above threshold value.
  • the threshold is always determined to have a finite value, and the degree of abnormality is determined.
  • the above purpose is to continuously apply impulse-like vibrations to the structure from an external vibration source when the inspection object does not have a vibration source such as a power section inside, Alternatively, it is achieved by measuring an attenuation component of a passing vibration waveform.
  • the above objective is to calculate the threshold excess time for the plurality of thresholds and then take this logarithm, and use the logarithmic trend of the threshold excess time as a logarithmic slope to calculate the normal part at the same threshold. This can be achieved by using the difference between the logarithmic slope due to the tendency of the logarithmic threshold excess time and the logarithmic slope due to the logarithmic threshold excess time of the abnormal part to determine the degree of deterioration of the abnormal part. Disclosure of the invention
  • a vibration inspection device that detects abnormalities in an inspection object
  • two or more threshold values are set in advance, vibration data of the inspection object is measured for a fixed time, and vibration waveform data obtained by digitizing the measured vibration data is displayed.
  • This is a vibration inspection method that calculates a time when a threshold value is exceeded, defines the time as exceeding the threshold value, and determines the degree of abnormality of the inspection object by comparing the magnitude of the value.
  • a value smaller than the maximum peak value of the vibration waveform data of the normal inspection object is set as one threshold value, and a second threshold value is set to a value equal to or less than the above threshold value.
  • the logarithm is taken, and the tendency of the logarithmic threshold excess time is defined as a logarithmic slope, and the logarithmic threshold of the normal part at the same threshold is calculated.
  • a vibration inspection method characterized in that the degree of deterioration of an abnormal portion is determined by using a difference between a logarithmic slope due to a tendency of a value excess time and a logarithmic slope due to a time exceeding a logarithmic threshold of an abnormal portion.
  • a vibration detector a measuring device that replaces the vibration data from the vibration detector with digital information and generates vibration waveform data, and an absolute value of the vibration waveform data from the measuring device, and a preset threshold value
  • a determination device that calculates the time or number of times exceeding the threshold value and compares it with a reference value to determine the degree of abnormality
  • a display that displays data from the determination device as an abnormality degree or waveform.
  • a vibration inspection device BRIEF DESCRIPTION OF THE FIGURES
  • FIG. 1 is a configuration diagram of a vibration inspection device according to one embodiment of the present invention.
  • FIG. 2 is a functional configuration diagram of the vibration inspection device shown in FIG.
  • FIG. 3 is a flowchart showing a vibration analysis procedure of the vibration inspection device shown in FIG.
  • Figure 4 shows (a) the observed vibration data (vibration waveform), (b) the converted waveform with the frequency selection range set to the high frequency band, and (c) the converted waveform with the frequency selection range set to the low frequency band. (D) This is a screen on which threshold values are set.
  • Figure 5 shows (a) the vibration data (vibration waveform) of the deteriorated pump A, (b) the vibration waveform of the normal pump B, and (c) the waveform obtained by converting the vibration waveform of the pump A in the high frequency region.
  • Pump B vibration waveform converted in the high frequency range (e) Pump A squared waveform screen with threshold set, and ( ⁇ ) Pump ⁇ squared waveform This is a screen on which values are set.
  • Figure 6 is a screen showing (a) a waveform obtained by expanding the vibration waveform of the deteriorated pump A over time and (b) a waveform obtained by expanding the vibration waveform of the normal pump B.
  • Figure 7 shows (a) a partially enlarged waveform of the vibration waveform of the deteriorated pump A, (b) a noise waveform, and (c) a waveform in which noise is superimposed on the vibration waveform.
  • FIG. 8 is a table showing the values obtained by comparing the absolute value data with the two threshold values and calculating the threshold excess time and the ratio of the slope of the same section in the pumps A and B.
  • Figure 9 is a table in which the slope is calculated for each two thresholds, and the degree of deterioration is calculated based on the difference between the slopes.
  • FIG. 1 is a configuration diagram of a vibration inspection device according to one embodiment of the present invention.
  • the rotating device 5 is operated, and if there is an abnormality in the rotating device, the vibration generated from the abnormal location is propagated to the surface of the device such as the bearing portion 6, and the vibration is measured, and the degree of the abnormality is measured. This is an example of determining.
  • the vibration inspection device includes a vibration detector (sensor) 1 and a vibration device 2 that sets the vibration waveform data from the vibration data obtained by taking in the vibration data of the vibration detector 1 and measuring the vibration data. It comprises a determination device 3 for determining the degree of abnormality from the measurement result, and a display device 4 for displaying the measurement result, the determination result, and the like.
  • the vibration detector (sensor) 1 is mounted on the rotating device bearing 6, and measures vibration data transmitted from inside.
  • the measuring device 2 replaces the vibration data observed by the vibration detector (sensor) 1 with digital information to obtain a vibration waveform.
  • the judging device 3 converts the vibration waveform data into an absolute value, compares the absolute value with a preset threshold value, calculates the time or number of times the threshold value is exceeded, and compares the calculated time or number with a reference value to determine the degree of abnormality.
  • the display device 4 displays the degree of abnormality and the waveform.
  • FIG. 2 is a functional block diagram of the vibration inspection device shown in FIG.
  • the detection signal (vibration data) of the vibration detector (sensor) 1 is first input to the amplifier 21 of the measuring device 2, and the analog-to-digital conversion process is performed by the AD converter 22. Then, the data is stored in the memory 124 as vibration waveform data d 1 cut out at regular time intervals by the data cut-out function 23.
  • the judgment device 3 calls the vibration waveform data d1 stored in the memory 24 and inputs it to the signal processing unit 31.
  • the vibration waveform data d1 is accumulated in the vibration waveform data file 32 so that it can be used as many times as necessary for subsequent arithmetic functions and output to a display device.
  • the frequency band is selected using the filter parameter 34 set each time.
  • the selected frequency band is sent to the absolute value function 35 as the filter-converted data d2, and becomes the absolute value data d3.
  • Absolute value data d 3 is compared with threshold value parameter 3 7 by threshold value comparison function 3 6, and the total time during which the absolute value data exceeds the threshold value is calculated. It is passed to the judgment unit 38 as d4.
  • the judging unit 38 compares the preset reference value with the above threshold excess time data d4, and accumulates the judgment result data d5 in the management data file 39.
  • the display device 4 calls the threshold excess time data d 4 for the necessary period set separately by the rough setting means 41 from the judgment result data file, and displays the time-series graph via the time-series display 4 2. Display or call the judgment result data d5 and display the management table via the management table display 43.
  • the vibration waveform data d1 at the required time is called from the vibration waveform data file 32 and displayed via the waveform display 44.
  • FIG. 3 is a flowchart of a processing procedure showing a vibration analysis procedure in the vibration inspection apparatus having the above-described configuration.
  • the vibration measured by the vibration detector (sensor) 1 is taken into the amplifier / AD converter as vibration data (vibration waveform signal).
  • step S1 the amplitude processing and AD com- plex processing are performed, and the analog signal is converted into a digital signal. It is converted into a signal and stored in the memory as vibration waveform data D1.
  • the measurement data is cut out in step S1 according to the measurement time set in advance in step S2 so that the comparison for a fixed time can always be performed.
  • step S1 the next measurement instruction signal (described later) is available.
  • the type of abnormalities to be determined next is selected by the step S3 determination control function.
  • the step S 4 parameter setting function selects one filter and one parameter from a preset filter and parameter group, and passes it to the digital filter processing in step S 5.
  • the measurement data is converted to filter conversion data d2 by convolution (inner product operation) filtering so that signals in a specific frequency band can pass easily. You.
  • the filter-to-converted data d2 is sent in the absolute value conversion processing in step S6 to become absolute value data d3.
  • step S3 the type of abnormality to be determined by the determination control function is selected, and in response to the selection of the determination type, step S7: in the threshold setting function, a pair of preset threshold parameters is set. (2) thresholds are selected and passed to step S8, threshold excess time calculation process.
  • step S3 the above process is repeated to select the type of abnormality to be determined each time, and a threshold is set for all types of abnormality N for one vibration waveform data d1 measured by one data acquisition. The above process is repeated until the calculation of the value excess time is completed.
  • step S 9 2 N pieces of threshold excess time data are generated.
  • step S 9 2 N pieces of data are used as N pieces of vector data, and N kinds of judgment reference vectors are used. Calculate the judgment result by comparing with and record it in the judgment result data file.
  • step S3 after all the processing for one vibration waveform data d1 is completed, after a certain time has elapsed by the timer for measuring the measurement cycle time, the process proceeds to step S1. The next measurement is instructed, and a new measurement is performed in step S1 based on the measurement instruction signal.
  • the step S10 display control has screen display parameter input means, and when it is started, a signal is sent to the step S11 graph setting, step S12 management table setting, step S13 waveform setting. Forward to indicate start of display.
  • step S11 graph setting and step S12 management table setting and step S13 waveform setting data file reading information and display conditions are created using preset values via step S10 display control. I do.
  • the S14 graph creation reads the information of the data file according to the read information generated in the step S11 graph setting and arranges the screen.
  • step S15 the management table is prepared by reading the information of the data file based on the read information generated in step S11 management table setting and arranging the screen.
  • step S16 waveform creation the information in the data file is read and the screen is arranged according to the read information generated in step S13 waveform setting.
  • the step S17 graph display displays the screen arrangement information of the step S14 graph generation.
  • the step S18 management table display displays screen arrangement information for creating the step S15 management table.
  • Step 19 The waveform display displays the screen array information for creating the step S16 waveform.
  • the abnormalities and deformations of the inspection object are determined based on the plurality of thresholds and the time when the absolute value of the vibration waveform is compared with a plurality of thresholds L and values. In this case, it is determined that the vector formed by the plurality of thresholds deviates from the vector in the normal state and approximates the vector for abnormality determination. If the rate at which the amplitude waveform contains an attenuated waveform that attenuates over time is large, it can be determined as if it were a betattle. Hereinafter, this method will be described in detail together with the observed waveform.
  • the vibration inspection method if the object to be inspected is a rotating device, the vibration data generated when the abnormal part makes rotational contact with another structure or slides is measured as sound or the amplitude of acceleration, and is used as vibration waveform data. .
  • many abnormalities of rotating equipment are mostly caused by bearing steel ball scratches inside the mechanism.
  • a sound sensor such as a noise meter as a vibration detector.
  • an acceleration sensor is attached directly to the mechanism as a vibration detector and installed, and L and weak vibrations that cannot be captured as sound are measured as needed.
  • the vibration waveform of the impact vibration generated by other mechanisms passing through or reflected at the deformed part of the structure and returning to the surface is measured. It is operationally advantageous to avoid inefficient contact type measurement and to perform non-contact measurement because a wide and efficient inspection is required. For this reason, in general, it is possible to obtain the vibration of the required amplitude by adjusting the impact force.
  • Figure 4 (b) shows the time waveform resulting from selecting the frequency band in step S5.
  • the purpose of frequency selection is mainly to remove noise and extract a specific frequency signal, and a filter used at that time is often a high-order filter having a sharp frequency characteristic.
  • the purpose of noise removal is not the purpose but to estimate the cause of the abnormality.
  • the passband is set to 12 KHz to 15 KHz
  • this method uses 10 KHz to 4 OKHz ( )).
  • Fig. 4 (b) is an example of setting from 10KHz to 4OKHz (upper limit)
  • Fig. 4 (c) is a waveform set from 50KHz to 250KHz.
  • Step S4 The parameter setting is to set the pass band of the filter for each cause to be determined, and the frequency band is divided as follows. For example, when vibrations occur due to minute scratches on bearings or high-speed rotation of hydraulic impellers, ultrasonic vibration is generated when the structure undergoes shear failure. However, when the ultrasonic vibration propagates through the mechanism, it spreads to the vibration in the lower audible region, and the component in the narrow ultrasonic region is lost. In this case, if the filter is set to 2 OKHz or more, the diffused audible vibration can be captured.
  • the attenuation component of the waveform in the high frequency band of 1 kHz or more, it is possible to determine and detect the presence of a relatively large flaw on the rotating shaft rolling surface member.
  • the presence of eccentricity of the rotating shaft, load imbalance of the rotating shaft, and the like can be determined by extracting the attenuation component of the waveform in a low frequency band of 1 KHz or less.
  • step S3 determination control and searching for the parameter is added, different models can be inspected continuously.
  • the vibration waveform data is amplitude data of both polarities, but when comparing by a threshold value, processing is easier with unipolarity. For this reason, the step S6 absolute value processing may not be necessary, but the processing is easier and the amount of information is doubled compared to the case without the absolute value processing, and the improvement in accuracy can be expected.
  • Fig. 4 (d) is a waveform obtained by squaring Fig. 4 (b) and making it an absolute value.
  • a horizontal line is drawn, and the horizontal line illustrates the threshold value.
  • the calculation of the threshold excess time is performed by comparing the waveform with the threshold, recording the coordinates at which the waveform is greater than the threshold and the coordinates at which the waveform is smaller than the threshold, and taking the total length of time.
  • the sum of the time lengths may be set as the threshold excess time.
  • the observation waveform extraction length varies depending on the model to be inspected, and the above threshold excess time is determined by the ratio of the observation time corresponding to the standard time length. It is desirable to normalize the excess time.
  • the number of times the threshold has been exceeded can be easily obtained from the coordinates where the threshold has been exceeded, and a dog-like waveform having a relatively large peak value such as cavitation observation of a water treatment pump can be obtained. It has the effect of comparing the number of cavitation pulses when mixing.
  • a plurality of threshold values are used, but one threshold value is used for one calculation, and the same processing is performed after calculating the threshold excess time and changing the threshold value.
  • Fig. 5 shows a comparison between the waveform of a pump of the same model that has deteriorated and that of a normal pump.
  • Fig. 5 (a) shows the vibration waveform data of pump A several years after operation
  • Fig. 5 (b) shows the waveform immediately after operation.
  • 7 shows the vibration waveform data of the normal pump B.
  • Figures 5 (c) and 5 (d) show the filter-converted waveforms converted by setting the frequency band to high frequency, respectively.
  • the amplitude of the deteriorated waveform of Pump A is more than twice as large as the amplitude of Pump B.
  • Fig. 5 (e) and (f) show the absolute value data obtained by squaring the filter-converted waveform.
  • the above difference is further expanded to 6 times or more.
  • the table shown in Fig. 8 shows the values obtained by comparing the above absolute value data with the two threshold values and calculating the threshold excess time and the ratio of the slope of the same section for pumps A and B. Waveforms that appear to have a difference of about 6 times when compared to the absolute value data show a very large difference of 20 to 700 times when compared using threshold values.
  • the reason for this is that when the time axis is enlarged and displayed as shown in Figs. 6 (a) and (b), the deteriorated pump A has multiple intermittent amplitude peaks larger than the pump B amplitude.
  • the decay waveform is such that the amplitude gradually decays in the time after each peak.
  • Pump B has small peaks but no attenuating shape, and many sporadic noise-like waveforms. Due to such a difference in the tendency, the waveform of pump A, which has a large amount of damping component, has both the excess time of threshold 1 and the excess time of threshold 2 larger than that of pump B. The tendency becomes remarkable. Conversely, when the threshold value is reduced, the difference between the two becomes smaller.
  • the threshold excess time be less than one-hundredth of the observation time for a normal pump B waveform. From the experience that the difference between the users is rapidly reduced, if the threshold overtime is set to be more than one-hundredth of the observation time, a sufficient difference can be observed between Pump A and Pump B, and a good judgment can be made. be able to. Also, if the threshold is set to a large value, the irregularities in the waveform of the pump B will often exceed the threshold only for noise-like waveforms, which is not preferable for comparison. For example, even if the sampling period is 500,000 times every second, if the observation time is 1 second and the threshold excess time is 1 / 10,000, the data exceeds the threshold only 5 times.
  • step S9 what is the normal condition of pump B. If the condition is the same in the environment, the degree of deterioration is determined by comparing directly with the threshold excess time shown in Table 1. Is also good. However, the elasticity of the floor differs between the first floor and the second floor of the building, and there is a difference even if the installation work cannot be completed, which may result in a difference in the propagation characteristics inside the pump. In such a case, it is necessary to set a threshold value in accordance with the amplitude of the pump separately for each site, and there is a problem that the criterion cannot be determined.
  • the slope is calculated for every two thresholds, and the degree of deterioration is determined based on the difference between the slopes. If the threshold exceeded time obtained from the observation is determined to be 1/100, 1/100, 1/100 of the observation time, the threshold is determined.
  • the determination threshold value is determined based on the rule when the determination is made.
  • the threshold value is determined based on the height of the peak of the normal waveform, and the waveform exceeding the threshold value initially set due to effects such as lubrication due to the driving experience of the pump itself is determined at a later date.
  • Time to exceed threshold without leaving is determined based on the height of the peak of the normal waveform, and the waveform exceeding the threshold value initially set due to effects such as lubrication due to the driving experience of the pump itself is determined at a later date.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本発明は建造物内に発生したクラックや回転機器の劣化を観測波形から安定かつ精度良く検知する振動検出方法を提供することにある。このため、本発明は、検査対象物である回転機器5の劣化を検出する振動検査装置において、軸受け部6に設置した振動検出器(センサ)1と、振動検出器(センサ)1で検出した振動データをデジタル化する計測装置2と、計測装置2より出力する振動波形データを判定する判定装置と、判定装置からの出力を表示する表示装置とを備え、予め2つ以上のしきい値を設定し、検査対象物の振動データを一定時間分計測し、計測した振動データをデジタル化した振動波形データがしきい値を超過する時間を算出し、しきい値超過時間とし、この値の大小を比較することにより検査対象物の異状の度合いを判定することを特徴としている。

Description

振動検出方法及びその装置 技術分野
本発明は、 コンクリート構造物内部の空洞ゃクラックなどの異状部が存在する か否かを検査する検査方法およびその装置に関するものである。 特に、 トンネル 内を移動しながら機械的打撃を繰り返し構造物壁面に与える機構を用いるように し、 打撃の反響音を解析して、 異状部の位置、 大きさなどを高速度で特定するの に好適な振動検査方法およびその装置に関するものである。 また、 本発明は、 ポ ンプゃモーターなどの回転機器の回転軸、 '軸受け、 羽根車などに異状が存在する か否かを判定する検査方法およびその装置にも適用できる。 特に、 対象が発する 暗騒音、 振動の波形を解析し異状の種類などを高精度に特定するのに好適な振動 検査方法およびその装置に関するものである。 背景技術 '
検査対象の異状や変状を音響や振動を用 、て検査する従来の振動検査方法にお いては、 大別して周波数分析法と時間軸分析法とがある。 周波数分析法は F F T (高速フーリエ変換) などの方法により波形の振幅をスぺクトル成分に変換し、 それを周波数軸上にプロッ トして成分強度分布表として観測できるようにするも のである。 成分強度分布を正常、 異状の判定に用いるためには、 最初に正常なデ 一タスぺクトル強度を成分強度分布表上にプロッ 卜し、 異状時に正常データの分 布が偏ってずれる周波数に着目する。 そして、 その周波数のスぺクトル成分強度 が変化する範囲を検知する検知窓を周波数の上限値および下限値、 そしてスぺク トル強度の上限値および下限値の 4つの数値で表す。 新たなデータのスぺグトル 成分を上記の検知窓の数値と比較し、 その範囲にスぺクトル成分が存在するか否 かを識別し、 識別結果から異状と正常とを判別する方法である。
また、 時間軸分析法には、 波形の実効値を比較する方法と、 波形の波高値とを 比較する方法と、 打撃を用いた打音解析のように 1インパク 卜に対する応答波形 の形状で比較する方法などが知られている。 波形の実効値を比較する方法では、 計測データにノイズが多く混入する場合、 ノィズの影響を最小化するために計測 データを周波数フィルターに通し特定の周波数帯域の成分を減衰させてノイズの 少なくなつた波形データを解析する方法が一般的である。 この方法はノイズの周 波数成分が特定の周波数領域に偏っていて必要な信号の周波数成分の周波数領域 と大きく異なる場合には有効であり、 特定の周波数帯域を除去する周波数フィル ターを用意する。 次に 1インパク卜に対する応答波形の形状で比較する方法は構 造的な変状による反射波などが現れやすくその影響を忠実に反映できる方法であ しかし上述した従来の振動検査方法のうち周波数分析法では、 データに混入す るノィズが多い場合、 ノィズの影響を最小化するためには極端に長く観測し分析 するデータ長を大きくすることが必要となる。 こうすることでランダムに変化す るノイズ成分が平均化されて平均強度が小さくなり、 信号成分は繰り返し現れる ために平均化しても強度は変化しないことを利用している。 このような方法では 観測時間を長時間に設定する必要があることから、 観測時間が極端に長くなり、 点検周期が伸び巡回的に監視する回転機器台数を削減する必要が生じたり、 トン ネル内の移動点検では緻密な点検ができなくなるなどの問題がある。 あるいは解 析機器の演算能力の向上が必要になつて機器費用が高くなるなどの問題がある。 また時間軸分析法のうち実行値を比較する方法では、 ノィズが広い周波数帯域 に分散して現れたり、 時間によって様々な周波数成分に変化したりする。 このた め、 周波数成分が特定の周波数領域に偏らないことがあり、 必要な信号の周波数 成分の周波数領域に重なることが多く、 ノイズを除去する周波数フィルターを用 いることができないという問題がある。
次に波高値を比較する方法では、 検査対象の設置される床の剛性の低下など設 置場所の違いから検査対象である商品の出荷検査時の振動値に比べ設置後の振動 値とが大幅に強く異なることがある。 このため、 従来のように 1つのしきい値で 波高値を比較する方法では、 正常な対象でも異状となってしまうことがある。 そ のため設置した現場において正常な振動波形の波高値に余裕を持ったしきい値を 設定するために十分長時間の観測が必要となつた。 同時に劣化した振動波形の予 測が困難となるため、 初期の劣化を検知するような微妙なしき Ϊ、値の設定は極め て困難となる問題があった。
次に 1ィンパクトに対する応答波形の形状で比較する方法は、 実験室などにお いては構造的な変状による反射波などが現れやすくその影響を忠実に反映できる 方法である。 し力、し、 コンクリ一ト構造物の対象はトンネルや橋 など騷音の多 い屋外にありノイズの影響を受けやすいという問題と同時に、 1インパクトごと に解析を行う方法では構造物表面を緻密にかつ効率よく精査しょうとすると単位 時間当たりの打撃数を増加する必要があり、 打撃能力を高め打撃間隔を短縮し、 それにあわせて解析に必要な演算能力も高くする必要が生じるなどの問題がある o
そこで、
本発明の第 1の目的は、 異状や変状の検出率が高く、 振動データに含まれる雑 音との識別性に優れた振動検査方法およびその装置を提供するにある。
本発明の第 2の目的は、 振動検査の結果を容易に確認できる振動検査法及びそ の装置を提供するにある .
本発明の第 3の目的は、 異状や変状の有無を自動的に判別することが可能な振 動検査方法およびその装置を提供するにある。
本発明の第 4の目的は、 異状や変状の自動判定を行うための判定基準を容易に 設定できることが可能な振動検査方法およびその装置を提供するにある。
本発明の第 5の目的は、 コンクリート構造物などの打音検査において簡単な機 構で連続打撃を与え、 複数の応答波形を一度に判定することが可能な振動検査方 法およびその装置を提供するにある。
上記目的は、 予め 2つ以上のしきい値を設定し、 検査対象物の振動データを一 定時間分計測し、 計測した振動デ一タをデジ夕ル化した振動波形デー夕がしきい 値を超過する時間を算出ししきい値超過時間としこの値の大小を比較することに より検査対象物の異状の度合いを判定することで達成される。
好適には、 1つめのしきい値として正常な検査対象物の振動波形データの最大 ピ一ク値より小さい値を設定し、 上記のしきい値以下の値に 2つめのしきい値を 設定して、 しきい値超過時間の大小を比較することにより検査対象物の異状の度 合いを判定することを特徴とする。
さらに好適には 1つめのしきい値として正常な検査対象物の振動波形データの 最大ピーク値より小さい値を設定し、 上記のしきい値以下の値に 2つめのしきい 値を設定して、 常にしきい値超過時間比率が有限の値を持つようにして、 異状の 度合いを判定することを特徴とする
上記目的は検査対象物が内部に動力部などの振動発生源を持たない構造物にお いては、 外部の振動発生源より構造物にインパルス状の振動を連続して与え、 異 状部で反射あるいは通過する振動波形の減衰成分を計測することで達成される。 上記目的は、 前記複数のしきい値に対するしきい値超過時間を算出した後にこ の対数を取り、 対数化したしきい値超過時間の傾向を対数傾きとし、 同じしきい 値での正常部の対数しきい値超過時間の傾向による対数傾きと異状部の対数しき い値超過時間による対数傾きとの差異を利用して異状部の劣化の度合いを判定す ることで達成できる。 発明の開示
上記目的を達成するために本発明が採用した技術解決手段は、
検査対象物の異状を検出する振動検査装置において、 予め 2つ以上のしきい値 を設定し、 検査対象物の振動データを一定時間分計測し、 計測した振動データを デジタル化した振動波形データがしきい値を超過する時間を算出し、 しきい値超 過時間とし、 この値の大小を比較することにより検査対象物の異状の度合いを判 定することを特徵とする振動検査方法である。
また、 1つめのしきい値として正常な検査対象物の振動波形データの最大ピー ク値より小ざい値を設定し、 上記のしきい値以下の値に 2つめのしきい値を設定 して、 しきい値超過時間の値の大小を比較することにより異状の度合いを判定す ること.を特徴とする振動検査方法である。
また、 1つのしきい値として正常な検査対象物の振動波形データの最大ピーク 値より小さい値を設定し、 上記のしきい値以下の値に 2つめのしきい値を設定し て、 2つ以上のしきい値超過時間よりしきい値超過時間比率に置き換え、 常にし きい値超過時間比率が有限の値を持つようにして、 異状の^合いを判定すること T/JP2002/007949 を特徴とする振動検査方法である。
また、 検査対象物が内部に動力部などの振動発生源を持たない構造物において は、 外部の振動発生源より構造物にインパルス状の振動を連続して与え、 異状部 で反射ある L、は通過する振動波形の減衰成分を一定時間計測し区間ごとの判定を 行えるようにすることを特徴とする振動検査方法である。
また、 前記複数のしきい値に対するしきい値超過時間を算出した後にこの対数 を取り、 対数化したしきい値超過時間の傾向を対数傾きとし、 同じしきい値での 正常部の対数しきい値超過時間の傾向による対数傾きと異状部の対数しきい値超 過時間による対数傾きとの差異を利用して異状部の劣化の度合いを判定すること を特徴とする振動検査方法である。
また、 振動検出器と、 前記振動検出器からの振動データをデジタル情報に置き なおして振動波形データとする計測装置と、 前記計測装置からの振動波形データ を絶対値化し予め設定されたしきい値と比較し、 しきい値を超過する時間あるい は回数を算出し基準値と比較して異状の度合いを定める判定装置と、 前記判定装 置からのデータを異状の度合いや波形として表示する表示装置とから構成される ことを特徴とする振動検査装置である。 図面の簡単な説明
図 1は本発明の一実施形態に係る振動検査装置の構成図である。
図 2は図 1に示す振動検査装置の機能構成図である。
図 3は図 1に示す振動検査装置の振動解析手順を示す流れ図である。
図 4は (a ) 観測した振動データ (振動波形) と (b ) 周波数選択範囲を高周 波帯域に設定し変換した波形と (c ) 周波数選択範囲を低周波帯域に設定し変換 した波形と (d ) しきい値を設定した画面である。
図 5は (a ) 劣化の進んだポンプ Aの振動データ (振動波形) と (b ) 正常な ポンプ Bの振動波形と (c ) ポンプ Aの振動波形を高周波数領域で変換した波形 と (d ) ポンプ B振動波形を高周波数領域で変換した波形と (e ) ポンプ Aの 2 乗処理を行った波形にしきい値を設定した画面と (ί ) ポンプ Βの 2乗処理を行 つた波形にしきい値を設定した画面である。 図 6は (a ) 劣化の進んだポンプ Aの振動波形を時間の拡大を行った波形と ( b ) 正常なポンプ Bの振動波形を時間の拡大を行った波形とを示した画面である
図 7は (a ) 劣化の進んだポンプ Aの振動波形の一部を拡大し表示した波形と ( b ) ノイズ波形と (c ) 振動波形にノイズが重畳した波形である。
図 8は絶対値データと 2つのしきい値との比較を行いしきい値超過時間と、 ポ ンプ A、 Bで同じ区間の傾きの比とを算出した値を示す表である。
図 9はしきい値 2つごとに傾きを算出し、 その傾きの差異で劣化の度合いを算 出した表である。 発明を実施するための最良の形態
以下、 本発明の一実施形態を図面を参照して説明する。
図 1は、 本発明の一実施形態に係る振動検査装置の構成図である。 この実施形 態は、 回転機器 5を動作させ、 回転機器に異状があった場合、 異状個所より発生 した振動が軸受け部 6などの機器表面に伝播し現れた振動を計測し、 異状の度合 いを判定する例である。
この実施形態に係る振動検査装置は、 振動検出器 (センサ) 1とこの振動検出 器 1の振動データを取り込み計測した振動デ一夕から振動波形データを設定する 計測装置 2と、 計測装置 2の計測結果から異状の度合いを判定する判定装置 3と 、 計測結果や判定結果などを表示する表示装置 4とを備えて構成される。
振動検出器 (センサ) 1は回転機器軸受け部 6に装着され、 内部より伝播する 振動データを計測する。 計測装置 2は振動検出器 (センサ) 1により観測された 振動データをデジタル情報に置き直して振動波形デ一夕とする。 判定装置 3は上 記振動波形データを絶対値化し予め設定されたしきい値と比較し、 しきい値を超 過する時間あるいは回数を算出し基準値と比較して異状の度合いを定める。 表示 装置 4は異状の度合いや波形を表示する。
図 2は図 1に示す振動検査装置の機能プロック図である。
振動検出器 (センサ) 1の検出信号 (振動データ) は先ず計測装置 2の増幅器 2 1に入力され、 A Dコンバータ 2 2によってアナログ—ディジタル変換処理が行 われ、 ついでデータ切り出し機能 2 3にて一定時間間隔で切り出された振動波形 データ d 1としてメモリ一 2 4に蓄積される。
判定装置 3では、 メモリー 2 4に蓄積された振動波形データ d 1を呼び出し、 信号処理部 3 1に入力する。 この振動波形データ d 1は、 後続の演算機能や、 表 示機への出力などへ適宜何度でも利用できるようにするために振動波形デー夕フ アイル 3 2に蓄積される。 つぎに周波数帯域選別機能 3 3にて判定すべき異状の 種類ごとに、 毎回設定されるフィルターパラメータ 3 4を用いて周波数帯域を選 別する。 選別された周波数帯域はフィルタ一変換データ d 2として絶対値化機能 3 5に送られ、 絶対値データ d 3となる。 上記絶対値データ d 3はしきい値比較 機能 3 6にてしきい値パラメータ 3 7と比較され、 絶対値データがしきい値を超 えている時間の総和を算出し、 しきい値超過時間データ d 4として判定部 3 8に 渡す。 判定部 3 8は予め設定してある基準値と上記しきい値超過時間データ d 4 との比較を行い、 判定結果データ d 5を管理データファイル 3 9に蓄積する。 表示装置 4では判定結果データフアイルから別途ダラフ設定手段 4 1にて設定 される必要な期間分のしきい値超過時間デ一タ d 4を呼び出し、 時系列表示 4 2 を介して時系列グラフを表示し、 あるいは判定結果データ d 5を呼び出し管理表 表示 4 3を介して管理表を表示する。 また同様にじて必要な時点の振動波形デー タ d 1を振動波形デ一タファイル 3 2から呼び出して波形表示 4 4を介して表示 する。
図 3は、 上述した構成の振動検査装置における振動解析手順を示した処理手順 の流れ図である。 振動検出器 (センサ) 1で計測される振動は、 振動データ (振 動波形信号) として増幅/ A Dコンバータに取り込まれ、 先ずステップ S 1で增 幅処理および A Dコンパ一ト処理されアナログ信号からディジタル信号に変換さ れ、 振動波形データ D 1としてメモリに格納される。 このとき常に一定時間の比 較が行えるようにステップ S 2にて予め設定された計測時間に従ってステップ S 1で計測データが切り出される。 通常 A Dコンバート処理を行う市販の I Cチッ プなどのハードウヱァには入力データ数設定機能が付与されているので、 上記の 波形切り出し処理は設定された計測時間に匹敵するデータ数を上記 I Cチップに セッ 卜する方法を取る。 ステップ S 1では次回計測の指示信号 (後述) があるま „
PCT/JP2002/007949 では次の計測はおこなわないようにしメモリに格納された計測データは次回計測 データが入力されるときまで保持される。
つぎに判定すべき異状の種類が複数あるときには、 ステップ S 3判定制御機能 にて次回の判定すべき異状の種類を選択する。 この判定種類選択を受けて、 ステ ップ S 4パラメータ設定機能では予め設定されたフィルタ一パラメータ群より 1 つのフィルタ一パラメータを選択し、 ステップ S 5デジタルフィルター処理に渡 す。 そして I I R型デジタルフィルタ一処理の係数パターンとして使用され、 特 定の周波数帯域の信号が通過しやすいようにコンボリューシヨン (内積演算) フ イ リタリング処理により計測データがフィルター変換データ d 2に変換される。 つぎにフィルタ一変換データ d 2はステップ S 6の絶対値化処理にて送られ絶対 値データ d 3となる。 上記計測データをそのまま比較することも容易であるが好 適には絶対値化処理を行い振動波形データがしきい値を過ぎる確率を 2倍に上げ るようにし、 2つ以上のしきい値を超過した時間からしきい値超過時間比率を計 測できるようにすると精度の向上が期待できる。
上記ステップ S 3判定制御機能にて判定すべき異状の種類を選択し、 この判定 種類選択を受けて、 ステップ S 7しきい値設定機能では予め設定されしきい値パ ラメ一夕群より 1対 (2個) のしきい値を選択し、 ステップ S 8しきい値超過時 間算出処理に渡す。 ステップ S 8しきい値超過時間算出処理では絶対値データ d 3としきい値パラメータとを比較し絶対値データ d 3がしきい値を超えている時 間の総和を算出し、 2個のしきい値(thnl thn2)に対して 2個のしきい値超過時 間データ (d4nl d4n2) を算出する。 ステップ S 3では上記の処理を繰り返して その都度判定すべき異状の種類を選択し、 1回のデータ取得により計測された 1 つの振動波形データ d 1に対してすべての異状の種類 Nに関してしきい値超過時 間を算出し終えるまで、 上記処理を繰り返す。
以上の処理結果としてしきい値超過時間デ一夕が 2 N個生成され、 ステップ S 9判定処理においては 2 N個のデータを N個のべクトルデータとし、 N種類の判 定基準べク トルと比較して判定結果を算出し、 判定結果データファイルに記録す る。 ステップ S 3では 1つの振動波形データ d 1に対する全ての処理が終了した 後、 計測サイクル時間を計るタイマーによって一定時間経過後にステップ S 1に 次回計測の指示を行い、 ステップ S 1は計測指示信号により新たな計測を行う。 ステップ S 1 0表示制御は画面表示パラメ一夕入力手段を有し、 さらに起動さ れたときに、 ステップ S 1 1グラフ設定、 ステップ S 1 2管理表設定、 ステップ S 1 3波形設定に信号を送り、 表示の開始を示す。
ステップ S 1 1グラフ設定およびステップ S 1 2管理表設定およびステップ S 1 3波形設定ではステップ S 1 0表示制御を介して予め設定された値を利用して データフアイルの読取情報と表示条件を作成する。
ついで S 1 4グラフ作成はステップ S 1 1グラフ設定に生成された読取情報に 従いデ一タフアイルの情報を読み取り画面配列を行う。
ステップ S 1 5管理表作成はステップ S 1 1管理表設定に生成された読取情報 に徒いデ一タフアイルの情報を読み取り画面配列を行う。
ステップ S 1 6波形作成はステップ S 1 3波形設定に生成された読取情報に従 いデータフアイルの情報を読み取り画面配列を行う。
ステップ S 1 7グラフ表示はステップ S 1 4グラフ生成の画面配列情報を表示 する。
ステップ S 1 8管理表表示はステップ S 1 5管理表作成の画面配列情報を表示 する。
ステップ 1 9波形表示はステツプ S16波形作成の画面配列情報を表示する。 本発明において、 検査対象物の異状や変状は、 絶対値化した振動波形を複数の しき L、値で比較した場合、 複数のしき 、値超過時間で判定され、 ·異状や変状のあ る場合には上記複数のしきい値の作るベタトルが正常時のべク トルから乖離し、 異状判別用のべクトルに近似することで判定される。 振幅波形に時間的に減衰す る減衰波形の含まれる率が多ければベタトルのようにみなして判定できる。 以下 、 この方法を観測波形と共に詳述する。
振動検査方法において、 検査対象が回転機器であれば、 異状部分が他の構造物 と回転接触したり摺動する際に起きる振動データを音響もしくは加速度の振幅と して計測し振動波形データとする。 一般には回転機器の異状の多くは機構部内部 にあるベアリング鋼球傷が多く、 この検出では機構部の損傷を防ぐ程度の検査で は、 振動検出器として騷音計などの音響センサーで計測することが多い。 しかし 、 損傷に至るまでの寿命を測って保全計画を立てるような目的に対しては、 より 微小な傷が発生した時点で、 傷が転送面と衝突するときに起きる微弱な振動を計 測する必要があることから、 振動検出器として直接機構部に加速度センサ一を接 着して設置しておき、 音響として捉えられな L、微弱な振動を随時計測する。 またトンネルのようなコンクリート構造物の 查では、 他の機構より発せられ た打撃振動が構造物内の変状部分にて通過や反射し表面に戻ってきた振動波形を 計測するが、 対象壁面が広く効率の良い検査が求められるなどの理由から効率の 悪い接触式を避け、 非接触で計測するのが運用上有利である。 このため一般には 打撃力を調整して必要な振幅の振動を得られることから音響センサ一で計測され O o
以上のように振動検出器 (センサ) の種別で長所短所が異なるために増幅器入 力は両者とも交換して用いられるようにすることが望ましい。
振動検出器 (加速度センサ) の信号が増幅器及び ADコンバート処理を通して 得られる振動波形デ一夕 d 1は図 4 (a) に示す時間一加速度曲線として観測さ れ 。
ステップ S 5で周波数帯域を選択した結果の時間波形を図 4 (b) に示す。 従来は周波数選択の目的が主にノィズ除去と特定周波数信号の抽出とをかねた ものが多く、 そのときに用いられるフィルタ一は高次フィルターで周波数特性の 韋峻なものを用いることが多い。 しかし異状が上記に述べたような衝突に基づく ような場合その周波数帯域は広範になる。 また本方法においてはノイズ除去が目 的ではなく異状の原因を推定するにあるので、 例えば従来では通過帯域を 12K Hz〜l 5 KHzに設定した例に対して本方法では 10KHz〜4 OKHz (上 限) に設定するようにしている。 図 4 (b) は 10KHz〜4 OKHz (上限) と設定した例で、 図 4 (c) は 50KHz〜250 K H zと設定した波形である o
ステップ S 4パラメータ設定は判別する原因別にフィルターの通過領域を設定 するものであるが、 その周波数帯域の区分は次のようになされる。 例えばべァリ ングの微小傷発生や水理系の羽根車の高速回転などによって起きるキヤビテーシ ョン現象などの発生時には、 構造物がせん断破壊を起こす際に超音波振動が発生 し、 この超音波振動が機構部を伝播する際にさらに低い可聴領域の振動に拡散し 超音狭領域の成分は失われてゆく。 この場合フィルターを 2 O K H z以上に選択 設定を行つた場合、 上記拡散した可聴領域の振動を捉えることができる。
また 1 K H z以上の高い周波数帯域での波形の減衰成分を抽出するようにして 、 回転軸転がり面部材の比較的大きな傷の存在を判定し検出することができる。 あるいは 1 K H z以下の低い周波数帯域での波形の減衰成分を抽出するようにし て、 回転軸の偏芯、 回転軸の負荷アンバランスなどの存在を判定することができ る。
好適にはステップ S 3判定制御から別途検査対象機種を設定して該ノ、°ラメータ を検索するような処理機能を付加すると、 異機種を連続して点検することができ る。
振動波形データは両極性の振幅データであるが、 しきい値によって比較する場 合、 片側極性であるほうが処理が容易である。 このためにステップ S 6絶対値化 処理は無くても良いが、 有った方が無い場合に比べ処理が容易でかつ情報量が 2 倍になり精度の向上も期待できる。
好適には絶対値化処理を 2乗計算による絶対値化に変えて、 振幅の強弱をさら に大きく開くように強調し、 かつ振幅を物理的な比較に多く用いられるパワー ( =振幅2) として観測することができる。 図 4 ( d ) は図 4 ( b ) を 2乗して絶 対値化した波形である。
図 4 ( d ) には横線を入れてあり、 横線はしきい値を図示化したものである。 しきい値超過時間の算出は波形としきい値とを比較し、 波形がしきい値より大き くなる座標と小さくなる座標を記録し、 それらの時間長さ総計を取ったものであ る。
観測時の波形切り出しが常に一定時間に設定されていれば時間長さの総計をし きい値超過時間とすればよい。 .
好適には、 水処理ポンプのように大型で機種によつては極めて低周波領域での 判定が必要となる場合があり、 こういう場合にはある程度長時間の観測を行った 方が精度が向上するので、 検査対象機種により観測波形きり出し長さを変化し、 上記のしきい値超過時間は標準時間長さに対応する観測時間の比で上記しきい値 超過時間を正規化することが望まれる。
さらに好適には、 しきい値を超過した座標からしきい値超過回数も簡単に求め ることができ、 水処理ポンプのキヤビテーション観測のように比較的振幅のピー ク値の犬なる波形が混入する場合など、 キヤビテ一ションパルスの数量比較が行 えるなどの効果がある
ステップ S 8では複数のしきい値を用いるが、 一回の演算には 1つのしきい値 を用い、 しきい値超過時間を算出した後にしきい値を代えて同じ処理を行う。 図 5に同じ機種のポンプで劣化が進んだ物と正常品との波形比較を示す 図 5 ( a ) は稼動後数年を経過したポンプ Aの振動波形データを同図 (b ) は 稼動直後の正常なポンプ Bの振動波形データを示す。 目視で 1 . 3倍程度に差異 が有ると見えるが、 厳密な判定を行うには 2つの波形の振幅は大差が無いと考え るべきである。 図 5 ( c ) ( d ) はそれぞれ周波数帯域を高周波に設定して変換 したフィルタ一変換波形である。 ここで劣化が進んだポンプ Aの波形における振 幅がポンプ Bの振幅に対して 2倍強大きい差となっている。 次に図 5 ( e ) ( f ) にフィルタ一変換波形を 2乗した絶対値データを示すが、 ここでさらに上記の 差は 6倍以上に拡大している。
図 8に示す表に上記の絶対値データと 2つのしきい値との比較を行いしきい値 超過時間と、 ポンプ A、 Bで同じ区間の傾きの比とを算出した値を示す。 絶対値 データまでの比較では 6倍程度の差異に見えた波形が、 しきい値による比較を行 うと 2 0〜7 0 0倍と極めて大きな差異になる。
その理由は図 6 ( a ) ( b ) に示すように時間軸を拡大表示すると、 劣化の起 きているポンプ Aではポンプ Bの振幅より大きめの振幅ピークが間欠的に複数存 在し、 その各々のピーク以降の時間で振幅が徐々に減衰する減衰波形となつてい る。 またポンプ Bでは小さ目のピークがあるが減衰的な形状は見受けられず散発 的なノイズらしい波形が多い。 このような傾向の違いから、 減衰成分の多いボン プ Aの波形の方がポンプ Bよりもしきい値 1の超過時間もしきい値 2の超過時間 も共に多く、 かっしきい値を高くすればその傾向は顕著となる。 逆にしきい値の 値を小さく してゆくと両者の差異ほ次第に小さなものになる。 そのため好適には 正常なポンプ Bの波形でしきい値超過時間が観測時間の百分の一以下になると両 者の差異が急速に小さくなる経験から、 しきい地超過時間が観測時間の百分の一 以上になるように設定するとポンプ Aとポンプ Bとに十分なる差異が観測でき、 良好な判定を行うことができる。 またしきい値を課題に大きく設定するとポンプ Bの波形の中のィレギュラーな多くはノイズのような波形だけがしきい値を超過 するだけとなるので、 比較には好ましくない。 例えばサンプリング周期を 5 0 0 , 0 0 0回毎秒としても観測時間 1秒間のデータでしきい値超過時間が一万分の 一ではただの 5回だけデータがしきい値を超えたことになり、 ポンプ Bの様相を 正しく表しているとは言いがたい。 このため好適には観測時間とサンプリング周 期の積である観測データ数において 2桁以上のしきい値超過回数が得られるよう なしきい値を設定すれば良好な判定が確保できる。
このしきい値超過時間で比較する方法によれば図 7 ( a ) の減衰波形に同図 ( b ) のノイズが重畳し同図 (c ) のような波形になっても、 単に波形が細切れに なるだけでしきい値超過時間にはあまり変化が無く、 さらにこの波形が数百数千 個集まった振動波形データでは、 しきい値超過時間で考えれば、 それぞれの誤差 はキヤンセルされノイズ重畳前と大きな変化はなくなる。
ステップ S 9の判定においては、 ポンプ Bの正常な状態がどのような.環境にお いても同じであれば、 表 1に示すしきい値超過時間でそのまま比較して劣化の度 合いを決めても良い。 しかし設置場所が建物の 1階と 2階とでは床の弾性が異な り、 設置工事のでき不出来などでも違いが生じて.、 それがポンプ内部の伝播特性 の差異となることもある。 このような場合その現場ごとに別々にポンプの振幅に 合わせたしきい値設定が必要となってしまい、 判定基準は確定できないという問 題を生じる。
そこで好適には、 図 9に示す表に示すように、 しきい値 2つごとに傾きを算出 し、 その傾きの差異で劣化の度合いを判定することにし、 しきい値の規準を正常 波形の観測から得たしきい値超過時間が観測時間の 1 0 0 0 0分の 1、 1 0 0 0 分の 1、 1 0 0分の1、 となるようなしきい値に決定することにすると、 上記の ような問題を解決することができる。
さらに好適には、 正常なポンプ Bの 2点間傾きが時間 1 0倍あたり 1 0となる 区間で判定を行うといつた規則に基づいて判定しきい値を定めるようにすれば、
3 より正確な判定が可能となる。 具体的には、 正常波形のピークすれすれの高さで しきい値を決めてしまって、 後日になりボンプ自身の運転慣れによる潤滑化など の効果によつて当初設定したしきい値を越す波形が出ないでしきい値超過時間 =
0の状態となり比較がしにく くなるなどの弊害が解消できる。
以上、 本実施形態の一例を説明したが、 センサから検出した振動データの処理 、 判定装置はコンピュータを使用して行うが、 専用の回路で設定することも可能 であり、 また表示手段もプリンタ一など種々の表示装置を使用することができる 。 また、 本発明はその精神又は主要な特徴から逸脱することなく他の色々な形で 実施することができる。 また、 前述の実施例はあらゆる点で単なる例示に過ぎず 、 限定的に解釈してはならない。 産業上での利用可能性
本発明によれば、
異状や変状を高い検出率で検出できると共に、 雑音との識別性を高くすること ができまた、 振動検査の結果を容易に確認でき、 更に、 異状や変状の有無を自動 的に判別することが可能となり、 異状や変状の自動判定を行うための判定基準を 容易に設定することも可能である。 またコンクリ一卜構造物などの打音検査にお いて簡単な機構で連続打撃を与え、 複数の応答波形を一度に判定することも可能 となる。

Claims

• 請求の範囲
1 . 検査対象物の異状を検出する振動検査装置において、 予め 2つ以上のしき い値を設定し、 検査対象物の振動データを一定時間分計測し、 計測した振動デー タをデジタル化した振動波形データがしきい値を超過する時間を算出し、 しきい 値超過時間とし、 この値の大小を比較することにより検査対象物の異状の度合い を判定することを特徴とする振動検査方法。
2 . 請求項 1において、 1つめのしきい値として正常な検査対象物の振動波形 データの最大ピーク値より小さい値を設定し、 上記のしきい値以下の値に 2つめ のしきい値を設定して、 しき 、値超過時間の値の大小を比較することにより異状 の度合いを判定することを特徴とする振動検査方法。
3 . 請求項 1において、 1つのしきい値として正常な検査対象物の振動波形デ —夕の最大ピーク値より小さい値を設定し、 上記のしきい値以下の値に 2つめの しきい値を設定して、 2つ以上のしきい値超過時間よりしきい値超過時間比率に 置き換え、 常にしきい値超過時間比率が有限の値を持つようにして、 異状の度合 、を判定することを特徴とする振動検査方法。
4 . 請求項 1において、 検査対象物が内部に動力部などの振動発生源を持たな い構造物においては、 外部の振動発生源より構造物にィンパルス状の振動を連続 して与え、 異状部で反射あるいは通過する振動波形の減衰成分を一定時間計測し 区間ごとの判定を行えるようにすることを特徴とする振動検査方法。
5 . 請求項 1または請求項 2あるいは請求項 3のいずれかにおいて、 前記複数 のしき L、値に対するしき L、値超過時間を算出した後にこの対数を取り、 対数化し たしきい値超過時間の傾向を対数傾きとし、 同じしきい値での正常部の対数しき い値超過時間の傾向による対数傾きと異状部の対数しきい値超過時間による対数 傾きとの差異を利用して異状部の劣化の度合いを判定することを特徴とする振動 検査方法。
6 . 振動検出器と、 前記振動検出器からの振動データをデジタル情報に置きな おして振動波形デ一夕とする計測装置と、 前記計測装置からの振動波形デ一夕を 絶対値化し予め設定されたしきい値と比較し、 しきい値を超過する時間あるいは 回数を算出し基準値と比較して異状の度合いを定める判定装置と、 前記判定装置 からのデータを異状の度合いや波形として表示する表示装置とから構成されるこ とを特徴とする振動検査装置。
PCT/JP2002/007949 2002-08-05 2002-08-05 振動検出方法及びその装置 WO2004013589A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
PCT/JP2002/007949 WO2004013589A1 (ja) 2002-08-05 2002-08-05 振動検出方法及びその装置
AU2002323803A AU2002323803A1 (en) 2002-08-05 2002-08-05 Vibration detecting method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2002/007949 WO2004013589A1 (ja) 2002-08-05 2002-08-05 振動検出方法及びその装置

Publications (1)

Publication Number Publication Date
WO2004013589A1 true WO2004013589A1 (ja) 2004-02-12

Family

ID=31217255

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2002/007949 WO2004013589A1 (ja) 2002-08-05 2002-08-05 振動検出方法及びその装置

Country Status (2)

Country Link
AU (1) AU2002323803A1 (ja)
WO (1) WO2004013589A1 (ja)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102538954A (zh) * 2012-02-06 2012-07-04 成都阜特科技有限公司 基于直线电机的振动检定装置及其振动检定方法
CN103123279A (zh) * 2013-01-22 2013-05-29 宁波易荣机电科技有限公司 一种轴承振动自动测量仪
CN104165682A (zh) * 2014-08-05 2014-11-26 中国北车集团大连机车车辆有限公司 内燃机车通风机振动故障预警装置及方法
TWI686591B (zh) * 2018-11-29 2020-03-01 東訊股份有限公司 行動式振動檢測裝置及其檢測方法
CN112402736A (zh) * 2020-11-17 2021-02-26 杭州师范大学钱江学院 一种输液监视方法
US20220214311A1 (en) * 2019-04-25 2022-07-07 Thk Co., Ltd. Abnormality diagnosis system and abnormality diagnosis method

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10176948A (ja) * 1996-12-19 1998-06-30 Meidensha Corp 異常振動検出装置

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10176948A (ja) * 1996-12-19 1998-06-30 Meidensha Corp 異常振動検出装置

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102538954A (zh) * 2012-02-06 2012-07-04 成都阜特科技有限公司 基于直线电机的振动检定装置及其振动检定方法
CN102538954B (zh) * 2012-02-06 2015-04-08 成都阜特科技股份有限公司 基于直线电机的振动检定装置及其振动检定方法
CN103123279A (zh) * 2013-01-22 2013-05-29 宁波易荣机电科技有限公司 一种轴承振动自动测量仪
CN104165682A (zh) * 2014-08-05 2014-11-26 中国北车集团大连机车车辆有限公司 内燃机车通风机振动故障预警装置及方法
TWI686591B (zh) * 2018-11-29 2020-03-01 東訊股份有限公司 行動式振動檢測裝置及其檢測方法
US20220214311A1 (en) * 2019-04-25 2022-07-07 Thk Co., Ltd. Abnormality diagnosis system and abnormality diagnosis method
CN112402736A (zh) * 2020-11-17 2021-02-26 杭州师范大学钱江学院 一种输液监视方法

Also Published As

Publication number Publication date
AU2002323803A1 (en) 2004-02-23

Similar Documents

Publication Publication Date Title
CA2687785C (en) Parameter independent detection of rotating machinery faults
JP3321487B2 (ja) 機器/設備診断方法およびシステム
US20170241422A1 (en) Method and system for detection of faults in pump assembly via handheld communication device
CN105067239B (zh) 基于扫频激励振动的梁裂纹故障检测装置及方法
JP3046426B2 (ja) プラント機器の監視装置
Hadjileontiadis et al. Crack detection in beams using kurtosis
JP2008292288A (ja) 減速機の軸受診断装置
Klausen et al. Multi-band identification for enhancing bearing fault detection in variable speed conditions
TWI447371B (zh) 即時檢測系統和方法
JP6728808B2 (ja) 計測診断装置、及び計測診断方法
Nienhaus et al. Development of acoustic emission (AE) based defect parameters for slow rotating roller bearings
JP2002333437A (ja) 打撃検査装置
JP5133276B2 (ja) 構造物の欠陥検出方法
WO2004013589A1 (ja) 振動検出方法及びその装置
JPWO2004068078A1 (ja) 状態判定方法と状態予測方法及び装置
JP3920715B2 (ja) 振動信号の処理方法
KR20110132236A (ko) 베어링 상태 감시 방법 및 베어링 상태 감시 장치
JP6127283B2 (ja) コンクリート構造物の検査方法及びコンクリート構造物の検査装置
Jiang et al. Rolling bearing quality evaluation based on a morphological filter and a Kolmogorov complexity measure
JP4071161B2 (ja) 回転機械の劣化診断方法
JP2003202276A (ja) 機械設備又は機器の異常診断方法
JPH06323899A (ja) 低速回転機械の異常診断方法
JP2918683B2 (ja) 軸受疲労によるae信号の特定方法およびこの特定方法を用いた軸受破壊予知方法
JP2004279056A (ja) 回転機械の診断方法及びそのプログラム
Rodriguez et al. Comparative study between laser vibrometer and accelerometer measurements for mechanical fault detection of electric motors

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AU BR CA CN HU ID IL IN JP KR MX PL SG SK US VN

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR IE IT LU MC NL PT SE SK TR

121 Ep: the epo has been informed by wipo that ep was designated in this application
122 Ep: pct application non-entry in european phase
NENP Non-entry into the national phase

Ref country code: JP

WWW Wipo information: withdrawn in national office

Country of ref document: JP