WO2021117752A1 - Procédé de surveillance d'état de palier à roulement et dispositif de surveillance d'état de palier à roulement - Google Patents

Procédé de surveillance d'état de palier à roulement et dispositif de surveillance d'état de palier à roulement Download PDF

Info

Publication number
WO2021117752A1
WO2021117752A1 PCT/JP2020/045772 JP2020045772W WO2021117752A1 WO 2021117752 A1 WO2021117752 A1 WO 2021117752A1 JP 2020045772 W JP2020045772 W JP 2020045772W WO 2021117752 A1 WO2021117752 A1 WO 2021117752A1
Authority
WO
WIPO (PCT)
Prior art keywords
rolling bearing
regression model
damage
rolling
generating
Prior art date
Application number
PCT/JP2020/045772
Other languages
English (en)
Japanese (ja)
Inventor
正嗣 北井
Original Assignee
Ntn株式会社
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
Priority claimed from JP2019223524A external-priority patent/JP2021092455A/ja
Priority claimed from JP2019226287A external-priority patent/JP2021096102A/ja
Application filed by Ntn株式会社 filed Critical Ntn株式会社
Publication of WO2021117752A1 publication Critical patent/WO2021117752A1/fr

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16—ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16C—SHAFTS; FLEXIBLE SHAFTS; ELEMENTS OR CRANKSHAFT MECHANISMS; ROTARY BODIES OTHER THAN GEARING ELEMENTS; BEARINGS
    • F16C33/00—Parts of bearings; Special methods for making bearings or parts thereof
    • F16C33/30—Parts of ball or roller bearings
    • F16C33/58—Raceways; Race rings
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00—Testing of machine parts
    • G01M13/04—Bearings
    • G01M13/045—Acoustic or vibration analysis

Definitions

  • the present invention relates to a rolling bearing condition monitoring method and a rolling bearing condition monitoring device.
  • Rolling bearings that are not easy to replace and are used under relatively low speed conditions, such as spindle bearings for wind power generators, are often used continuously even if initial damage occurs. On the other hand, if the size of the damaged portion increases and the displacement between the inner and outer rings of the bearing increases and the displacement exceeds the permissible range, the device itself including the rolling bearing may be destroyed.
  • Patent Document 1 determines the relationship between a state variable including output data of a sensor that detects the state of an industrial machine or an ambient environment and a failure occurrence, that is, a failure condition. learn. As a result, the machine learning device can output the presence / absence of failure or the degree of failure in response to the input state variable.
  • Patent Document 1 describes the remaining life of a rolling bearing from the time of damage diagnosis until the damaged part progresses to the allowable limit size based on the effective vibration value obtained by the vibration sensor. The theoretical formula to predict is described.
  • Patent Document 1 the degree of failure is exemplified as an index value, and a configuration for directly estimating the size of the damaged portion, which indicates the state of damage, is not shown.
  • Patent Document 2 since the remaining life of the actually damaged rolling bearing is not measured, the accuracy of the above theoretical formula has not been sufficiently confirmed. Therefore, there is no guarantee that the remaining life can be accurately estimated, especially when it is predicted that it is difficult to express by a simple theoretical formula such that the form of damage progress changes as the operation time increases. I am concerned.
  • a main object of the present invention is to provide a condition monitoring method and a condition monitoring device capable of accurately estimating the size of a damaged part of a rolling bearing or the remaining life of a rolling bearing.
  • the method of monitoring the state of a rolling bearing is a step of generating a regression model by machine learning based on the vibration data of the first rolling bearing and a step of generating a regression model by using the regression model of the second rolling bearing. It includes a step to monitor the condition.
  • Each of the first rolling bearing and the second rolling bearing is arranged between the rotating wheel, the stationary wheel, and the rotating wheel and the stationary wheel, and moves on the raceway surface of the stationary wheel as the rotating wheel rotates. It has a plurality of rolling elements.
  • the step of generating the regression model is the vibration data of the first rolling bearing and the damage of the first rolling bearing.
  • the step of monitoring the state of the second rolling bearing is a step of acquiring the vibration data of the second rolling bearing and a feature quantity vector of the second rolling bearing composed of the feature quantity of the vibration data of the second rolling bearing. It includes a step of generating and a step of estimating the damage size of the second rolling bearing from the feature vector of the second rolling bearing using a regression model.
  • the present invention it is possible to provide a rolling bearing condition monitoring method and a condition monitoring device capable of accurately estimating the size of a damaged portion of a rolling bearing or the remaining life of a rolling bearing.
  • FIG. 1 is a diagram showing a cross-sectional view of a rolling bearing monitored by the condition monitoring device according to the first embodiment of the present invention and a block diagram of the condition monitoring device.
  • a rolling bearing is formed by a roller bearing in which the outer ring is a rotating wheel is typically described, but the scope of application of the present invention is limited to such a bearing condition monitoring device.
  • the rolling bearing to be monitored may have an inner ring of a rotating wheel, a ball bearing, or the like.
  • the rolling bearing 10 includes an inner ring 12, an outer ring 16, and a plurality of rolling elements 18.
  • the inner ring 12 is fitted onto the non-rotating shaft body 14.
  • the outer ring 16 is provided on the outer peripheral side of the inner ring 12 and rotates integrally with a rotating body (not shown).
  • Each of the plurality of rolling elements 18 is a cylindrical "roller", and is interposed between the inner ring 12 and the outer ring 16 while being held at equal intervals with the adjacent rolling elements by a cage (not shown).
  • the inner ring 12 receives a radial load from a plurality of rolling elements 18 that are passing through the load range.
  • a load region is formed on the upper side in the vertical direction (on the paper surface) with respect to the center of rotation. Then, in the normal state in which the raceway surface (outer peripheral surface of the inner ring 12) 20 is not damaged, the inner ring 12 is located at the center of the load region located in the vertically upward direction of the inner ring central axis among the plurality of rolling elements 18. Receives maximum load from passing rolling elements.
  • Such rolling bearings 10 may continue to be used even after initial damage such as small-scale peeling if they are used under relatively low speed conditions and / or if they are not easy to replace. There are many. On the other hand, if the size of the damaged portion is increased and the displacement between the inner and outer rings of the bearing is increased and the displacement exceeds the allowable range, the device including the rolling bearing 10 may be destroyed.
  • the rolling bearing 10 in order to use the rolling bearing 10 to the limit of the safe range, it is necessary to grasp the size of the damaged part. However, it is difficult to measure the size of the damaged part during use. Therefore, it is desirable to estimate the size of the damaged portion by measuring measurement data that reflects the size of the damaged portion, which can be easily measured.
  • condition monitoring device 100 of the rolling bearing 10 accurately estimates the size of the damaged portion of the rolling bearing 10 from the measurement data that can be easily measured.
  • the condition monitoring device 100 is configured to accurately estimate the circumferential length of the damaged portion from the measurement data of the vibration acceleration of the rolling bearing 10 by using a regression model by machine learning. ..
  • condition monitoring device 100 for monitoring the state of the rolling bearing 10 includes an acceleration sensor 102 and a diagnostic unit 106.
  • the x direction corresponds to the circumferential direction of the rolling bearing 10
  • the y direction corresponds to the axial direction of the rolling bearing 10.
  • the X direction and the Y direction correspond to the measurement directions of the acceleration sensor 102.
  • the acceleration sensor 102 is fixedly installed in a place where the vibration acceleration of the bearing in the bearing box or the outer diameter surface of the equipment can be measured.
  • the Y direction corresponds to the vertical direction (or radial direction) of the outer peripheral surface (raceway surface 20) of the inner ring 12 of the rolling bearing 10.
  • the X direction corresponds to the horizontal direction (or horizontal direction) at the installation position of the acceleration sensor 102.
  • the acceleration sensor 102 is a sensor for detecting the acceleration of vibration of the rolling bearing 10 in the vertical direction (radial direction).
  • the acceleration of this vibration is also referred to as “vibration acceleration” or simply “acceleration”.
  • the acceleration sensor 102 may be configured to detect the vibration acceleration in the horizontal direction (horizontal direction).
  • the diagnosis unit 106 monitors and diagnoses the state of the rolling bearing 10 based on the detected value of the acceleration sensor 102. Specifically, the diagnostic unit 106 estimates the size of the damaged portion of the rolling bearing 10 based on the detected value of the acceleration sensor 102. Hereinafter, the size of the damaged portion of the rolling bearing 10 will be described.
  • (Explanation of damage size) 2 to 4 are views showing a stage in which the size of the damaged portion of the rolling bearing 10 progresses.
  • the load range of the rolling bearing 10 is enlarged and shown.
  • the y direction corresponds to the axial direction of the rolling bearing 10
  • the x direction corresponds to the circumferential direction of the rolling bearing 10.
  • FIG. 2 shows an example in which the initial damage D1 occurs on the raceway surface 20 of the inner ring 12.
  • stage 0 the state before the initial damage to the raceway surface 20
  • stage 1 the state in which the raceway surface 20 is initially damaged
  • the size of the damaged portion refers to the circumferential length of the damaged portion formed on the raceway surface of the inner ring 12 (stationary wheel).
  • the size of the damaged portion indicates the maximum value of the length of the damaged portion in the circumferential direction.
  • the size of such a damaged portion is also simply referred to as “damage size”. That is, the damage size of the initial damage D1 is sz1 in FIG.
  • FIG. 3 is a diagram showing a state in which the initial damage generated on the raceway surface 20 is progressing in the axial direction. Since the portion adjacent to the initial damage D1 as shown in FIG. 2 in the circumferential direction (x direction in FIG. 3) is the portion where the rolling element comes into contact with the rolling element before or after passing through the damaged portion. , The rolling element and the raceway surface 20 of the inner ring 12 come into contact with each other on the entire surface. Therefore, the stress distribution on the contact surface becomes uniform. On the other hand, the portion adjacent to the initial damage D1 in the axial direction (y direction in FIG. 3) is a portion that comes into contact with the rolling element while the rolling element is passing through the damaged portion.
  • the initial damage D1 basically progresses along the axial direction. Damage in which such initial damage spreads mainly in the axial direction is referred to as "damage D2".
  • the damage size of damage D2 is sz2 in FIG.
  • stage 2 a state in which the initial damage mainly progresses in the axial direction, that is, a state in which the damage D2 occurs (FIG. 3) is referred to as “stage 2”.
  • FIG. 4 is a diagram showing a state in which the damage is progressing in the circumferential direction.
  • the damage D2 shown in FIG. 3 expands over the entire axial direction, the damage then expands in the circumferential direction. After such damage spreads over the entire axial direction, the damage that spreads in the circumferential direction is referred to as "damage D3".
  • the damage size of damage D3 is sz3 in FIG.
  • stage 3 a state in which the damage spreads in the entire axial direction and then progresses in the circumferential direction.
  • FIG. 5 shows the experimental results of measuring the change in damage size after a slight initial damage D1 occurs on a part of the raceway surface 20 in the plurality of rolling bearings 10.
  • the damage size of each rolling bearing 10 is determined by, for example, periodically stopping the device including a plurality of rolling bearings 10, photographing the stationary wheels of each rolling bearing 10 with a camera, and performing image processing on the photographed image. Can be done.
  • FIG. 5 shows the operating time of the rolling bearing 10 after the initial damage occurs on the raceway surface 20, and the vertical axis shows the damage size.
  • the data before the initial damage occurred is omitted.
  • FIG. 5 shows the experimental results until the initial damage progresses and the rolling bearing 10 reaches the usage limit.
  • the reference point of the usage limit is, for example, the time when the amount of displacement between the inner and outer rings of the bearing due to damage reaches the upper limit of the allowable range of the clearance in the equipment.
  • FIG. 6 is a diagram for explaining the time course of damage size in FIG.
  • FIG. 6 illustrates one of the plurality of measurement data of FIG.
  • the times t1 and t2 in FIG. 6 indicate the times when the stage 2 and the stage 3 were entered, respectively.
  • the time t3 in FIG. 6 indicates the time when the rolling bearing 10 reaches the usage limit, that is, the time when the damage size reaches the upper limit of the allowable range.
  • the damage size sz2 indicating the circumferential length of the damage D2 does not change significantly from the damage size sz1 of the damage D1. Therefore, as shown in FIG. 6, the propagation speed of the damage size in the circumferential direction is slow and slow.
  • stage 3 when the damage spreads over the entire axial direction, the damage further spreads in the circumferential direction (stage 3).
  • the rolling element 18 passing through the damage is in a no-load state in which a load is hardly received from the raceway surface. Therefore, the vibration of the rolling element 18 becomes intense, and as shown in FIG. 6, the propagation speed of the damage size in the circumferential direction becomes faster than that of the stage 2.
  • the reference point for the usage limit of the rolling bearing 10 is, for example, the time when the amount of displacement between the inner and outer rings of the bearing due to damage becomes the upper limit of the allowable range of clearance in the equipment. Further, it may be a time point when the damage size becomes a constant multiple of the interval of the rolling elements 18. With this setting, for example, in a certain rolling bearing 10, when the damage size becomes 1.5 to 2 times the distance between the rolling elements, the two rolling elements 18 passing through the damaged part receive almost no load. As a result, it is possible to detect the time when the vibration of the rolling bearing 10 becomes too large.
  • the rolling bearing 10 (hereinafter, also referred to as “learning sample”) for generating a regression model by machine learning is determined.
  • the remaining one rolling bearing 10 was used as a "test sample” for estimating the damage size using the regression model.
  • the learning sample corresponds to one embodiment of the "first rolling bearing”.
  • the test sample corresponds to one embodiment of "second rolling bearing”.
  • the diagnostic unit 106 separates the vibration acceleration data of the learning sample by frequency filtering. Specifically, the vibration acceleration data is separated into data corresponding to 20 to 1000 Hz, 1000 to 5000 Hz, 5000 to 20000 Hz by a low-pass filter, a band-pass filter, and a high-pass filter.
  • the diagnostic unit 106 divides the vibration acceleration data separated for each frequency into a plurality of segments, and calculates a feature quantity vector composed of a combination of a plurality of feature quantities for each segment.
  • FIG. 7 is a diagram for explaining a segment and a feature amount vector.
  • the measurement data of the time length T1 is divided into segments of the time length T2 shorter than the time length T1.
  • T2 is set to a length that is a constant multiple of the time required for one rotation of the rotation axis.
  • T2 is set at 0.6 seconds, which corresponds to the time required for the rotation axis to rotate five times.
  • a feature vector is generated for each of the further divided data.
  • the feature quantity vector treats a plurality of feature quantities as a set of vectors.
  • the feature amount is, for example, the level and / or shape of the waveform such as effective value, maximum value, crest factor, kurtosis, and skewness for at least one parameter of time, frequency, and skewness. It can be a parameter representing the characteristics of.
  • the feature amount may be processed by various frequency filters. Further, the feature amount extracted from the image data obtained by the time frequency analysis from the vibration data may be used.
  • a feature vector for generating a regression model by machine learning is generated.
  • the vibration acceleration data of the test sample is also converted into a feature vector by the same preprocessing and then input to the regression model.
  • the diagnostic unit 106 obtains the damage size in each segment based on the damage size data measured in the training sample. Specifically, the damage size at the time corresponding to each segment is obtained. As a result, it is possible to obtain a combination of the feature amount vector and the damage size corresponding to each segment. That is, time-dependent data showing the relationship between the feature vector and the damage size at each stage in which the damage progresses can be obtained.
  • FIG. 8 is a diagram for explaining the generation of a regression model by machine learning and the estimation of the damage size using the regression model.
  • the diagnostic unit 106 inputs the feature vector of the learning sample for each segment, outputs the damage size, and generates a regression model (first regression model) by machine learning.
  • a regression model for example, regression methods such as support vector machine, random forest, kernel ridge regression, deep learning regression, and combinations thereof can be used.
  • the first regression model is performed from the initial stage of damage (that is, stage 2) in which the damage progress rate is slow, the final stage of damage (that is, stage 3) in which the damage progress rate is fast, and after the initial damage occurs. Generated for each of the whole (ie, stages 2 and 3) up to the operable limit.
  • the diagnostic unit 106 generates a first regression model in which the damage size can be automatically estimated by inputting the feature vector (see FIG. 8).
  • the diagnostic unit 106 holds the first regression model in a storage unit (not shown) so that it can be read out and used at any time.
  • the corresponding feature vector of the test sample is applied to the first regression models of the above three types, the initial stage of damage, the final stage of damage, and the whole. This gives an estimate of the damage size of the test sample as the output of the first regression model.
  • FIG. 9 is a diagram for explaining the evaluation of the estimation accuracy of the first regression model.
  • the root mean square error square root and the coefficient of determination of the damage size estimated by the first regression model according to the first embodiment and the actually measured damage size are shown.
  • the root mean square error and the coefficient of determination when a multiple regression model is used instead of the regression model by machine learning are also shown.
  • the root mean square error and the coefficient of determination are calculated for each of the three types of first regression models described above. That is, the statistical value is calculated when the first regression model is generated in each of the initial stage of damage (stage 2), the final stage of damage (stage 3), and all of them (stage 2 and stage 3). There is.
  • the "root mean square error square root” is the square root obtained by averaging the squares of the error between the measured value of the damage size and the estimated value. That is, the root mean square error square root indicates that the smaller the value, the higher the accuracy of estimating the damage size.
  • the value of the root mean square error square root is smaller in the embodiment than in the comparative example. That is, it is shown that the estimation accuracy of the embodiment is higher than that of the comparative example. In addition, it is considered that the estimation accuracy is improved because the value of the root mean square error is greatly reduced, especially at the end of the damage.
  • the "coefficient of determination” is an index showing how well the estimated value obtained by regression analysis matches the actual measured value. Specifically, the closer the coefficient of determination is to 1, the closer the estimated value and the measured value are.
  • the coefficient of determination value of the embodiment is closer to 1 than that of the comparative example. That is, it is shown that the estimation accuracy of the embodiment is higher than that of the comparative example. In addition, it is considered that the estimation accuracy is improved because the value of the coefficient of determination increases significantly, especially in the final stage of injury.
  • the damage size can be estimated with higher accuracy than the comparative example using the multiple regression analysis.
  • FIG. 10 is an example of a flowchart for explaining the process of generating a regression model by machine learning.
  • the flowchart shown in FIG. 10 is executed by the diagnostic unit 106.
  • the diagnostic unit 106 acquires the vibration acceleration data and the damage size in the training sample. Specifically, for example, the diagnostic unit 106 receives vibration acceleration data from the acceleration sensor 102 wirelessly or by wire. Further, for example, the diagnostic unit 106 receives the damage size data wirelessly or by wire from the device in which the damage size value at each measurement time measured by the user is input.
  • step S02 the diagnostic unit 106 divides the vibration data into a plurality of segments.
  • the segment is set to, for example, a length that is a constant multiple of the time required for one rotation of the rotation axis.
  • step S03 the diagnostic unit 106 generates a feature amount vector for each segment.
  • the feature amount vector treats a plurality of feature amounts of vibration data as a set of vectors.
  • step S04 the diagnostic unit 106 obtains the corresponding damage size for each segment.
  • step S05 the diagnostic unit 106 uses machine learning to generate a damage size regression model (first regression model) that inputs the feature vector and outputs the damage size, and ends the process.
  • the diagnostic unit 106 is configured to perform machine learning about the combination of the feature vector and the damage size in each segment of the whole (that is, stages 2 and 3) from the occurrence of the initial damage to the operable limit. Will be done.
  • FIG. 11 is an example of a flowchart for explaining the process of estimating the damage size using the generated regression model.
  • the flowchart shown in FIG. 11 is executed by the diagnostic unit 106.
  • step S21 the diagnostic unit 106 acquires vibration data in the test sample. Specifically, for example, the diagnostic unit 106 receives vibration acceleration data from the acceleration sensor 102 wirelessly or by wire.
  • step S22 the diagnostic unit 106 divides the vibration data into segments.
  • the segment is set to, for example, a length that is a constant multiple of the time required for one rotation of the rotation axis.
  • step S23 the diagnostic unit 106 generates a feature amount vector for each segment.
  • step S24 the diagnostic unit 106 inputs the feature quantity vector of the test sample into the regression model (first regression model) generated by using machine learning for the relationship between the vibration data of the training sample and the damage size, and the test sample. Get an estimate of the damage size of and finish the process.
  • This regression model is, for example, by performing machine learning about the combination of the feature vector and the damage size in each segment of the whole (that is, stages 2 and 3) from the occurrence of the above initial damage to the operable limit. It was generated.
  • the condition monitoring device 100 may be configured to notify the user of the rolling bearing 10 of the estimated damage size by using voice, visual display, or the like.
  • the damage size is used as an index of the damage status of the rolling bearing 10.
  • the rated operating torque, power consumption, the rate of wear debris mixed in grease, etc. can be used as an index of the damage status instead of the damage size. ..
  • the condition monitoring device 100 generates a damage size regression model based on machine learning based on the relationship between the vibration data of the rolling bearing 10 for learning and the damage size, and is used for evaluation. By inputting the vibration data of the rolling bearing 10 into the regression model, the damage size is estimated. With such a configuration, the condition monitoring device 100 can estimate the size of the damaged part of the rolling bearing more accurately than the estimation accuracy by the conventional multiple regression analysis based on the vibration data that can be easily measured. it can.
  • the rolling bearing condition monitoring device 100 described above can be applied to various mechanical devices, and is particularly suitable for condition monitoring of the spindle bearing of a wind power generator. That is, the spindle bearing of the wind power generator is not easy to replace, is used under relatively low speed conditions, and is often continuously used even if the bearing is damaged. For the spindle bearings of such wind power generation equipment, it is an issue to clarify the bearing replacement time due to damage.
  • FIG. 12 is a diagram schematically showing the configuration of a wind power generation device to which the rolling bearing condition monitoring device 100 according to the first embodiment is applied.
  • the wind power generator 210 includes a spindle 220, a blade 230, a speed increaser 240, a generator 250, a spindle bearing (hereinafter, simply referred to as “bearing”) 260, and an acceleration sensor. It includes 270 and a data processing device 280.
  • the speed increaser 240, the generator 250, the bearing 260, the acceleration sensor 270 and the data processing device 280 are housed in the nacelle 290, and the nacelle 290 is supported by the tower 300.
  • the spindle 220 enters the nacelle 290, is connected to the input shaft of the speed increaser 240, and is rotatably supported by the bearing 260. Then, the main shaft 220 transmits the rotational torque generated by the blade 230 that has received the wind power to the input shaft of the speed increaser 240.
  • the blade 230 is provided at the tip of the main shaft 220, converts wind power into rotational torque, and transmits the wind power to the main shaft 220.
  • the speed increaser 240 is provided between the spindle 220 and the generator 250, and increases the rotational speed of the spindle 220 to output to the generator 250.
  • the speed increaser 240 is configured by a gear speed increase mechanism including a planetary gear, an intermediate shaft, a high speed shaft, and the like.
  • a plurality of bearings for rotatably supporting a plurality of shafts are also provided in the speed increaser 240.
  • the generator 250 is connected to the output shaft of the speed increaser 240 and generates electricity by the rotational torque received from the speed increaser 240.
  • the generator 250 is composed of, for example, an induction generator, but the type of the generator 250 is not limited to this.
  • a bearing that rotatably supports the rotor is also provided in the generator 250.
  • the bearing 260 is fixed in the nacelle 290 and rotatably supports the spindle 220.
  • the bearing 260 is a rolling bearing, and is a bearing to be monitored by the condition monitoring device 100 according to this embodiment.
  • the bearing 260 is different from the rolling bearing 10 described below in FIG. 1 in that the inner ring is a rotating wheel and the outer ring is a stationary wheel.
  • the condition monitoring device 100 according to the above embodiment is described in this way. It is also applicable to bearing 260.
  • the outer ring is a stationary ring, the load region is formed on the outer ring below the central axis in the vertical direction, and the initial damage occurs on the raceway surface of the inner peripheral surface of the outer ring.
  • the acceleration sensor 270 is a sensor for detecting the acceleration in the vertical direction (or radial direction) of the bearing 260.
  • the acceleration sensor 270 is fixed at a place where the vibration acceleration of the bearing on the outer diameter surface of the equipment or the bearing box can be measured, for example. Further, the acceleration sensor 270 may be a sensor for detecting the vibration acceleration in the horizontal direction (or horizontal direction) of the bearing 260.
  • the data processing device 280 is provided in the nacelle 290 and receives the detected value of the acceleration sensor 270. Then, the data processing device 280 monitors the state of the bearing 260 according to a preset program. Specifically, the data processing device 280 detects the vibration acceleration of the bearing 260 based on the detection value of the acceleration sensor 270. The data processing device 280 generates a feature vector based on the vibration acceleration. The data processing device 280 is configured to be capable of generating a regression model by machine learning that inputs the feature amount vector and outputs the damage size based on the feature amount vector and the damage size. The data processing device 280 is configured so that the damage size can be estimated by inputting the feature amount vector based on the regression model.
  • the data processing device 280 realizes the function of the diagnostic unit 106 (FIG. 1) described above. Further, the data processing device 280 and the acceleration sensor 270 constitute the above-mentioned condition monitoring device 100 (FIG. 1).
  • condition monitoring device 100 can estimate the damage size of the rolling bearing 10 with high accuracy based on the vibration acceleration of the rolling bearing 10 by using the regression model by machine learning. .. Therefore, it is possible to provide a condition monitoring method and a condition monitoring device capable of accurately estimating the size of the damaged portion of the rolling bearing.
  • condition monitoring device 100 In the second embodiment, a method of estimating the remaining life of the rolling bearing 10 based on the vibration acceleration of the rolling bearing 10 will be described using a regression model by machine learning. Since the condition monitoring device 100 according to the second embodiment has the same basic configuration as the condition monitoring device 100 shown in FIG. 1, the description will not be repeated.
  • the remaining one rolling bearing 10 was used as a "test sample” for estimating the remaining life using the regression model.
  • the learning sample corresponds to one embodiment of the "first rolling bearing”.
  • the test sample corresponds to one embodiment of "second rolling bearing”.
  • the diagnostic unit 106 separates the vibration acceleration data of the learning sample by frequency filtering. Specifically, the vibration acceleration data is separated into data corresponding to 20 to 1000 Hz, 1000 to 5000 Hz, 5000 to 20000 Hz by a low-pass filter, a band-pass filter, and a high-pass filter.
  • the diagnostic unit 106 divides the vibration acceleration data separated for each frequency into a plurality of segments, and calculates a feature quantity vector composed of a combination of a plurality of feature quantities for each segment.
  • the measurement data of the time length T1 is divided into segments of the time length T2 shorter than the time length T1.
  • T2 is set to a length that is a constant multiple of the time required for one rotation of the rotation axis.
  • T2 is set at 0.6 seconds, which corresponds to the time required for the rotation axis to rotate five times.
  • the feature quantity vector treats a plurality of feature quantities as a set of vectors.
  • the feature amount is, for example, the level and / or shape of the waveform such as effective value, maximum value, crest factor, kurtosis, and skewness for at least one parameter of time, frequency, and skewness. It can be a parameter representing the characteristics of.
  • the feature amount may be processed by various frequency filters. Further, the feature amount extracted from the image data obtained by the time frequency analysis from the vibration data may be used.
  • a feature vector for generating a regression model by machine learning is generated.
  • the vibration acceleration data of the test sample is also converted into a feature vector by the same preprocessing and then input to the regression model.
  • the diagnostic unit 106 obtains the damage size in each segment based on the damage size data measured in the training sample. Specifically, the damage size at the time corresponding to each segment is obtained. This damage size value can be used as an index of the reference of the usable limit of the rolling bearing 10 described above.
  • the diagnostic unit 106 obtains the remaining life in each segment.
  • the remaining life is the time from the segment to reach the standard of the usable limit.
  • the time until the reference of the usable limit is reached is the time until the damage size exceeds a constant multiple of the interval of the rolling elements 18 in the rolling bearing 10 in this experiment.
  • the remaining life is calculated by subtracting the operating time from the start of use of the rolling bearing 10 to the corresponding segment from the total operating time from the start of use of the rolling bearing 10 to the reaching of the limit standard in which the rolling bearing 10 can be used. Is sought after.
  • FIG. 13 is a diagram for explaining the generation of a regression model by machine learning and the estimation of the remaining life using the regression model.
  • the diagnostic unit 106 inputs the feature vector of the learning sample for each segment, outputs the remaining life, and generates a regression model (second regression model) by machine learning.
  • machine learning for example, regression methods such as support vector machine, random forest, kernel ridge regression, deep learning regression, and combinations thereof can be used.
  • the second regression model is used from the initial stage of damage (that is, stage 2) in which the damage progress rate is slow, the final stage of damage (that is, stage 3) in which the damage progress rate is fast, and after the initial damage occurs. Generated for each of the whole (ie, stages 2 and 3) up to the operable limit.
  • the diagnostic unit 106 generates a second regression model in which the remaining life can be automatically estimated by inputting the feature vector (see FIG. 13).
  • the diagnostic unit 106 holds the second regression model in a storage unit (not shown) so that it can be read out and used at any time.
  • the vibration acceleration data of the test sample is converted into a feature vector by the same preprocessing as the vibration acceleration data of the training sample.
  • the corresponding feature vector of the test sample is applied to the above three types of second regression models: initial damage, final damage, and the whole.
  • initial damage initial damage
  • final damage final damage
  • whole the whole of second regression models
  • FIG. 14 is a diagram for explaining the evaluation of the estimation accuracy of the regression model.
  • the root mean square error square root and the coefficient of determination of the residual life estimated by the regression model (second regression model) according to the second embodiment and the measured residual life are shown.
  • the actually measured remaining life or the measured value of the remaining life indicates the remaining life calculated from the total operating time that has reached the standard of the limit that can be actually used.
  • the root mean square error and the coefficient of determination when a multiple regression model is used instead of the regression model by machine learning are also shown.
  • the root mean square error and the coefficient of determination are calculated for each of the three types of second regression models described above. That is, the statistical value is calculated when a second regression model is generated in each of the initial stage of damage (stage 2), the final stage of damage (stage 3), and all of them (stage 2 and stage 3). There is.
  • the "root mean square error square root” is the square root obtained by averaging the square of the error between the measured value and the estimated value of the remaining life. That is, the root mean square error square root indicates that the smaller the value, the higher the accuracy of estimating the remaining life.
  • the value of the root mean square error square root is smaller in the embodiment than in the comparative example. That is, it is shown that the estimation accuracy of the embodiment is higher than that of the comparative example. In addition, it is considered that the estimation accuracy is improved because the value of the root mean square error is greatly reduced, especially at the end of the damage.
  • the "coefficient of determination” is an index showing how well the estimated value obtained by regression analysis matches the actual measured value. Specifically, the closer the coefficient of determination is to 1, the closer the estimated value and the measured value are.
  • the coefficient of determination value of the embodiment is closer to 1 than that of the comparative example. That is, it is shown that the estimation accuracy of the embodiment is higher than that of the comparative example. In addition, it is considered that the estimation accuracy is improved because the value of the coefficient of determination increases significantly, especially at the end of injury.
  • the remaining life can be estimated with higher accuracy than the comparative example using the multiple regression analysis.
  • FIG. 15 is an example of a flowchart for explaining the process of generating a regression model by machine learning.
  • the flowchart shown in FIG. 15 is executed by the diagnostic unit 106.
  • the diagnostic unit 106 acquires the vibration acceleration data in the learning sample and the total operating time when the bearing 10 reaches the usable limit reference. Specifically, for example, the diagnostic unit 106 receives vibration acceleration data from the acceleration sensor 102 wirelessly or by wire. Further, for example, the diagnostic unit 106 receives the data of the total operating time wirelessly or by wire from the device in which the total operating time is input by the user.
  • step S32 the diagnostic unit 106 divides the vibration data into a plurality of segments.
  • the segment is set to, for example, a length that is a constant multiple of the time required for one rotation of the rotation axis.
  • step S33 the diagnostic unit 106 generates a feature amount vector for each segment.
  • the feature amount vector treats a plurality of feature amounts of vibration data as a set of vectors.
  • step S34 the diagnostic unit 106 obtains the remaining life corresponding to each segment. Specifically, the diagnostic unit 106 calculates the remaining life by subtracting the operation time corresponding to the segment from the total operation time.
  • step S35 the diagnostic unit 106 uses machine learning to generate a regression model of the remaining life (second regression model) that inputs the feature vector and outputs the remaining life, and ends the process.
  • the diagnostic unit 106 is configured to perform machine learning about the combination of the feature vector and the remaining life in each segment of the whole (that is, stages 2 and 3) from the occurrence of the initial damage to the operable limit. Will be done.
  • FIG. 16 is an example of a flowchart for explaining a process of estimating the remaining life using the generated regression model.
  • the flowchart shown in FIG. 16 is executed by the diagnostic unit 106.
  • step S41 the diagnostic unit 106 acquires vibration data in the test sample. Specifically, for example, the diagnostic unit 106 receives vibration acceleration data from the acceleration sensor 102 wirelessly or by wire.
  • step S42 the diagnostic unit 106 divides the vibration data into segments.
  • the segment is set to, for example, a length that is a constant multiple of the time required for one rotation of the rotation axis.
  • step S43 the diagnostic unit 106 generates a feature amount vector for each segment.
  • step S44 the diagnostic unit 106 inputs the feature quantity vector of the test sample into the regression model (second regression model) generated by using machine learning for the relationship between the vibration data of the training sample and the remaining life, and the test sample.
  • the estimated value of the remaining life of is obtained, and the process is terminated.
  • the second regression model is, for example, machine learning about the combination of the feature vector and the remaining life in each segment of the whole (that is, stages 2 and 3) from the occurrence of the above initial damage to the operable limit. It is generated by doing.
  • the condition monitoring device 100 may be configured to notify the user of the rolling bearing 10 of the estimated remaining life by using voice, visual display, or the like.
  • the condition monitoring device 100 generates a regression model of the remaining life based on machine learning based on the relationship between the vibration data of the rolling bearing 10 for learning and the remaining life, and is used for evaluation. By inputting the vibration data of the rolling bearing 10 into the regression model, the remaining life is estimated. With such a configuration, the condition monitoring device 100 can estimate the remaining life of the rolling bearing more accurately than the estimation accuracy by the conventional multiple regression analysis based on the vibration data that can be easily measured.
  • the rolling bearing condition monitoring device 100 described above can be applied to various mechanical devices, and is particularly suitable for condition monitoring of the spindle bearing of a wind power generator. That is, the spindle bearing of the wind power generator is not easy to replace, is used under relatively low speed conditions, and is often continuously used even if the bearing is damaged. For the spindle bearings of such wind power generation equipment, it is an issue to clarify the bearing replacement time due to damage.
  • FIG. 17 is a diagram schematically showing the configuration of a wind power generation device to which the rolling bearing condition monitoring device 100 according to the second embodiment is applied.
  • the wind power generation device 210 according to the second embodiment replaces the data processing device 280 in the wind power generation device 210 shown in FIG. 12 with the data processing device 280A.
  • the description of the common parts with the first embodiment will not be repeated.
  • the data processing device 280A is provided in the nacelle 290 and receives the detected value of the acceleration sensor 270. Then, the data processing device 280A monitors the state of the bearing 260 according to a preset program. Specifically, the data processing device 280A detects the vibration acceleration of the bearing 260 based on the detection value of the acceleration sensor 270. The data processing device 280A generates a feature vector based on the vibration acceleration.
  • the data processing device 280 is configured to be capable of generating a regression model by machine learning based on the feature amount vector and the remaining life, with the feature amount vector as an input and the remaining life as an output. The data processing device 280 is configured so that the remaining life can be estimated by inputting the feature amount vector based on the regression model.
  • the data processing device 280A realizes the function of the diagnostic unit 106 (FIG. 1). Further, the data processing device 280A and the acceleration sensor 270 constitute the condition monitoring device according to the second embodiment.
  • condition monitoring device 100 can estimate the remaining life of the rolling bearing 10 with high accuracy based on the vibration acceleration of the rolling bearing 10 by using the regression model by machine learning. .. Therefore, it is possible to provide a condition monitoring method and a condition monitoring device capable of accurately estimating the remaining life of the rolling bearing.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Mechanical Engineering (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

La présente invention comprend : une étape consistant à générer des vecteurs de quantité de caractéristiques d'un premier palier à roulement, qui sont formés de quantités de caractéristiques de données de vibration sur le premier palier à roulement ; une étape consistant à effectuer un apprentissage machine pour les vecteurs de quantité de caractéristiques du premier palier à roulement et des tailles d'endommagement du premier palier à roulement pour générer un premier modèle de régression pour estimer les tailles d'endommagement à partir des quantités de caractéristiques ; une étape consistant à générer un vecteur de quantité de caractéristiques d'un second palier à roulement, qui est formé d'une quantité de caractéristiques de données de vibration sur le second palier à roulement ; et une étape consistant à estimer, à l'aide du premier modèle de régression, une taille d'endommagement du second palier à roulement à partir du vecteur de quantité de caractéristiques du second palier à roulement.
PCT/JP2020/045772 2019-12-11 2020-12-09 Procédé de surveillance d'état de palier à roulement et dispositif de surveillance d'état de palier à roulement WO2021117752A1 (fr)

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2019-223524 2019-12-11
JP2019223524A JP2021092455A (ja) 2019-12-11 2019-12-11 転がり軸受の状態監視方法及び転がり軸受の状態監視装置
JP2019226287A JP2021096102A (ja) 2019-12-16 2019-12-16 転がり軸受の状態監視方法及び転がり軸受の状態監視装置
JP2019-226287 2019-12-16

Publications (1)

Publication Number Publication Date
WO2021117752A1 true WO2021117752A1 (fr) 2021-06-17

Family

ID=76329932

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2020/045772 WO2021117752A1 (fr) 2019-12-11 2020-12-09 Procédé de surveillance d'état de palier à roulement et dispositif de surveillance d'état de palier à roulement

Country Status (1)

Country Link
WO (1) WO2021117752A1 (fr)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114674562A (zh) * 2022-03-22 2022-06-28 中车大连机车研究所有限公司 一种考虑寿命监测条件的轨道交通圆锥滚子轴承寿命预测方法
WO2024150672A1 (fr) * 2023-01-11 2024-07-18 日本精工株式会社 Procédé d'analyse de développement d'écaillage de bague de chemin de roulement de palier, dispositif d'analyse de développement d'écaillage, et programme

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140365411A1 (en) * 2013-06-05 2014-12-11 The Trustees Of Columbia University In The City Of New York Monitoring Health of Dynamic System Using Speaker Recognition Techniques
CN107144428A (zh) * 2017-03-17 2017-09-08 北京交通大学 一种基于故障诊断的轨道交通车辆轴承剩余寿命预测方法
WO2017159784A1 (fr) * 2016-03-17 2017-09-21 Ntn株式会社 Système de surveillance de conditions et dispositif de génération d'énergie éolien
WO2019026980A1 (fr) * 2017-08-04 2019-02-07 新日鐵住金株式会社 Dispositif de traitement d'informations, procédé de traitement d'informations et programme
JP2019045434A (ja) * 2017-09-07 2019-03-22 理研計器株式会社 ガス分析方法およびガス分析装置
JP2019045241A (ja) * 2017-08-31 2019-03-22 Ntn株式会社 転がり軸受の状態監視方法および状態監視装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140365411A1 (en) * 2013-06-05 2014-12-11 The Trustees Of Columbia University In The City Of New York Monitoring Health of Dynamic System Using Speaker Recognition Techniques
WO2017159784A1 (fr) * 2016-03-17 2017-09-21 Ntn株式会社 Système de surveillance de conditions et dispositif de génération d'énergie éolien
CN107144428A (zh) * 2017-03-17 2017-09-08 北京交通大学 一种基于故障诊断的轨道交通车辆轴承剩余寿命预测方法
WO2019026980A1 (fr) * 2017-08-04 2019-02-07 新日鐵住金株式会社 Dispositif de traitement d'informations, procédé de traitement d'informations et programme
JP2019045241A (ja) * 2017-08-31 2019-03-22 Ntn株式会社 転がり軸受の状態監視方法および状態監視装置
JP2019045434A (ja) * 2017-09-07 2019-03-22 理研計器株式会社 ガス分析方法およびガス分析装置

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114674562A (zh) * 2022-03-22 2022-06-28 中车大连机车研究所有限公司 一种考虑寿命监测条件的轨道交通圆锥滚子轴承寿命预测方法
CN114674562B (zh) * 2022-03-22 2022-10-25 中车大连机车研究所有限公司 一种考虑寿命监测条件的轨道交通圆锥滚子轴承寿命预测方法
WO2024150672A1 (fr) * 2023-01-11 2024-07-18 日本精工株式会社 Procédé d'analyse de développement d'écaillage de bague de chemin de roulement de palier, dispositif d'analyse de développement d'écaillage, et programme

Similar Documents

Publication Publication Date Title
JP6196093B2 (ja) 軸受装置の振動解析方法、軸受装置の振動解析装置、および転がり軸受の状態監視装置
JP6553970B2 (ja) 異常診断装置およびセンサ外れ検知方法
JP6605865B2 (ja) 転がり軸受の状態監視装置及び転がり軸受の状態監視方法
JP2018141804A (ja) 回転機械の残存耐用年数を測定する方法
WO2021117752A1 (fr) Procédé de surveillance d'état de palier à roulement et dispositif de surveillance d'état de palier à roulement
WO2016042983A1 (fr) Dispositif de diagnostic d'anomalie de palier à roulement, générateur d'énergie éolienne, et procédé de diagnostic d'anomalie de palier à roulement
WO2011081085A1 (fr) Dispositif de diagnostic d'anomalie destiné à un roulement, éolienne électrique et système de diagnostic d'anomalie
JP6650030B2 (ja) 転がり軸受疲労状態予測装置及び転がり軸受疲労状態予測方法
JP7027782B2 (ja) 転がり軸受の異常診断装置
JP6958068B2 (ja) 回転機械設備の異常診断システム及び異常診断方法
WO2017145687A1 (fr) Dispositif de diagnostic d'anomalie et procédé de diagnostic d'anomalie
JP6374234B2 (ja) 状態監視システム及びそれを備えた風力発電システム
JP2021096102A (ja) 転がり軸受の状態監視方法及び転がり軸受の状態監視装置
JPWO2020040280A1 (ja) 転がり軸受の異常診断方法及び異常診断装置
JP6654002B2 (ja) 転がり軸受の異常判定しきい値の設定方法、及び転がり軸受の状態監視方法
JP6192414B2 (ja) 転がり軸受の状態監視装置
JP6192413B2 (ja) 軸受装置の振動解析方法、軸受装置の振動解析装置、および転がり軸受の状態監視装置
JP6686625B2 (ja) 転がり軸受診断装置
JP7113668B2 (ja) 状態監視システム
JP6824076B2 (ja) 状態監視システムおよび風力発電装置
JP7098399B2 (ja) 状態監視装置および状態監視方法
JP6639266B2 (ja) 異常診断装置および異常診断方法
JP2021092455A (ja) 転がり軸受の状態監視方法及び転がり軸受の状態監視装置
JP6639265B2 (ja) 異常診断装置および異常診断方法
JP2015175828A (ja) 転動装置の状態監視装置およびそれを備える風力発電設備

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 20899063

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 20899063

Country of ref document: EP

Kind code of ref document: A1