JP2022024573A - Device and method of abnormality detection for rolling bearings - Google Patents

Device and method of abnormality detection for rolling bearings Download PDF

Info

Publication number
JP2022024573A
JP2022024573A JP2020127243A JP2020127243A JP2022024573A JP 2022024573 A JP2022024573 A JP 2022024573A JP 2020127243 A JP2020127243 A JP 2020127243A JP 2020127243 A JP2020127243 A JP 2020127243A JP 2022024573 A JP2022024573 A JP 2022024573A
Authority
JP
Japan
Prior art keywords
vibration
rotation speed
frequency
rolling bearing
abnormality
Prior art date
Legal status (The legal status 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 status listed.)
Granted
Application number
JP2020127243A
Other languages
Japanese (ja)
Other versions
JP7394031B2 (en
Inventor
豊美 吉田
Toyomi Yoshida
真 辺見
Makoto Henmi
靖 早坂
Yasushi Hayasaka
貴之 小山
Takayuki Koyama
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Industrial Products Ltd
Original Assignee
Hitachi Industrial Products 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 Hitachi Industrial Products Ltd filed Critical Hitachi Industrial Products Ltd
Priority to JP2020127243A priority Critical patent/JP7394031B2/en
Publication of JP2022024573A publication Critical patent/JP2022024573A/en
Application granted granted Critical
Publication of JP7394031B2 publication Critical patent/JP7394031B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

To provide a device and a method of abnormality detection for rolling bearing capable of calculating the rotational speed of a rotary machine which is supported by rolling bearings and whose rotational speed changes.SOLUTION: An abnormality detection device comprises: a vibration detector 15 that detects vibration of rolling bearings 20 supporting a rotary machine 30 whose rotational speed changes; a data acquisition unit 11 that acquires a waveform in a time domain of the vibration detected by the vibration detector 15; a rotational speed calculation unit 12 that finds the envelope of the waveform in the time domain, converts the waveform of the envelope into intensity distribution data in a frequency domain, and finds the rotational speed of the rotary machine 30 using the intensity distribution data in the frequency domain; and an abnormality specifying unit 13 that specifies abnormal components of the rolling bearings 20 using the rotational speed found by the rotational speed calculation unit 12, and the characteristic frequency of the rolling bearings 20. The rotational speed calculation unit 12 finds, from the intensity distribution data in the frequency domain, vibration intensities at a plurality of frequencies proportional to the reference frequency in a predetermined range, and finds the rotational speed using the found intensities.SELECTED DRAWING: Figure 1

Description

本発明は、転がり軸受の異常検出装置と異常検出方法に関し、特に、回転速度を求めて転がり軸受の異常を検出する装置と方法に関する。 The present invention relates to an abnormality detecting device and a method for detecting an abnormality in a rolling bearing, and more particularly to an apparatus and a method for detecting an abnormality in a rolling bearing by obtaining a rotation speed.

転がり軸受の異常(例えば、損傷や欠陥など)は、回転速度と振動に基づいて検出することができる。回転速度が変化する回転機械を転がり軸受が支持している場合では、転がり軸受の異常の検出に際して、計測した振動データを分析するためには、回転機械の回転速度も同時に得る必要がある。しかし、回転機械の回転速度を計測して収集するには、データロガーのチャンネルを1つ増やしてデータロガーから回転速度の情報を取得したり、回転機械に回転速度センサを設置したりする必要があるため、使用する装置や素子が増えるという課題がある。 Anomalies in rolling bearings (eg, damage or defects) can be detected based on rotational speed and vibration. When a rolling bearing supports a rotating machine whose rotation speed changes, it is necessary to obtain the rotation speed of the rotating machine at the same time in order to analyze the measured vibration data when detecting an abnormality in the rolling bearing. However, in order to measure and collect the rotation speed of the rotating machine, it is necessary to increase the channel of the data logger by one to acquire the rotation speed information from the data logger, or to install the rotation speed sensor in the rotating machine. Therefore, there is a problem that the number of devices and elements used increases.

このような課題を解決する技術として、例えば、特許文献1に記載された転がり軸受装置がある。特許文献1に記載された転がり軸受装置は、転がり軸受で検出した振動波形を利用して回転数を算出する装置であり、固定輪に取り付けられている振動検出器からの検出信号に基づいて軸受本体の振動及び回転数を算出する演算部を備える。演算部は、検出信号から得られる振動の周波数成分と、軸受本体の固有の振動回数とに基づいて軸受本体の回転数を算出する。 As a technique for solving such a problem, for example, there is a rolling bearing device described in Patent Document 1. The rolling bearing device described in Patent Document 1 is a device that calculates the number of rotations by using the vibration waveform detected by the rolling bearing, and is a bearing based on a detection signal from a vibration detector attached to a fixed wheel. It is equipped with a calculation unit that calculates the vibration and rotation speed of the main body. The calculation unit calculates the rotation speed of the bearing body based on the frequency component of the vibration obtained from the detection signal and the natural number of vibrations of the bearing body.

特開2017-133580号公報Japanese Unexamined Patent Publication No. 2017-133580

特許文献1に記載された技術などの従来の技術では、データロガーから回転速度の情報を取得したり、回転機械に回転速度センサを設置したりしなくても、計測した振動波形を用いて回転数を算出することで、回転機械の回転速度を推定することができる。しかし、特許文献1に記載された技術のような振動波形の周波数分析を用いる回転速度の推定では、アンバランスが小さく振動が小さい回転軸の場合や、回転速度が筐体などの構造物の固有振動の周波数に近い場合などは、回転速度に対応する周波数における振動の強度が必ずしも最大にならず、回転速度を正確に推定できないことがある。 In the conventional technology such as the technology described in Patent Document 1, the measured vibration waveform is used for rotation without acquiring the rotation speed information from the data logger or installing the rotation speed sensor in the rotating machine. By calculating the number, the rotation speed of the rotating machine can be estimated. However, in the estimation of the rotation speed using the frequency analysis of the vibration waveform as in the technique described in Patent Document 1, in the case of a rotation axis with a small imbalance and a small vibration, or when the rotation speed is unique to a structure such as a housing. When the frequency is close to the vibration frequency, the vibration intensity at the frequency corresponding to the rotation speed is not always maximized, and the rotation speed may not be estimated accurately.

本発明の目的は、転がり軸受で支持されて回転速度が変化する回転機械の回転速度を算出できる、転がり軸受の異常検出装置と異常検出方法を提供することである。 An object of the present invention is to provide an abnormality detection device for a rolling bearing and an abnormality detection method capable of calculating the rotation speed of a rotating machine supported by a rolling bearing and whose rotation speed changes.

本発明による異常検出装置は、回転速度が変化する回転機械を支持する転がり軸受に設置され、前記転がり軸受の振動を検出する振動検出器と、前記振動検出器が検出した前記振動の時間領域の波形を取得するデータ取得部と、時間領域の前記波形の包絡線を求め、高速フーリエ変換処理により時間領域の前記包絡線の波形を周波数領域の強度分布データに変換し、周波数領域の前記強度分布データを用いて前記回転機械の回転速度を求める回転速度算出部と、前記回転速度算出部が求めた前記回転速度と、前記転がり軸受の振動の要因を示す特徴周波数とを用いて、前記転がり軸受の異常部品を特定する異常特定部とを備える。前記回転速度算出部は、周波数領域の前記強度分布データについて、予め定めた範囲内の周波数である基準周波数に比例する複数の周波数における振動の強度を求め、求めた前記強度を用いて前記回転速度を求める。 The abnormality detection device according to the present invention is installed in a rolling bearing that supports a rotating machine whose rotation speed changes, and has a vibration detector that detects the vibration of the rolling bearing and a time region of the vibration detected by the vibration detector. The data acquisition unit for acquiring the waveform and the envelope of the waveform in the time region are obtained, and the waveform of the envelope in the time region is converted into the intensity distribution data in the frequency region by high-speed Fourier conversion processing, and the intensity distribution in the frequency region is obtained. The rolling bearing using a rotation speed calculation unit that obtains the rotation speed of the rotating machine using data, the rotation speed obtained by the rotation speed calculation unit, and a characteristic frequency indicating a factor of vibration of the rolling bearing. It is provided with an abnormality specifying part that identifies an abnormal part of the above. The rotation speed calculation unit obtains the vibration intensity at a plurality of frequencies proportional to the reference frequency, which is a frequency within a predetermined range, with respect to the intensity distribution data in the frequency domain, and uses the obtained intensity to obtain the rotation speed. Ask for.

本発明による異常検出方法は、異常検出装置が、回転速度が変化する回転機械を支持する転がり軸受の振動を検出する振動検出器が検出した前記振動の時間領域の波形を取得するデータ取得工程と、異常検出装置が、時間領域の前記波形の包絡線を求め、高速フーリエ変換処理により時間領域の前記包絡線の波形を周波数領域の強度分布データに変換し、周波数領域の前記強度分布データを用いて前記回転機械の回転速度を求める回転速度算出工程と、異常検出装置が、前記回転速度算出工程で求めた前記回転速度と、前記転がり軸受の振動の要因を示す特徴周波数とを用いて、前記転がり軸受の異常部品を特定する異常特定工程とを有する。前記回転速度算出工程では、異常検出装置が、周波数領域の前記強度分布データについて、予め定めた範囲内の周波数である基準周波数に比例する複数の周波数における振動の強度を求め、求めた前記強度を用いて前記回転速度を求める。 The abnormality detection method according to the present invention includes a data acquisition step of acquiring a waveform in the time region of the vibration detected by the vibration detector that detects the vibration of the rolling bearing that supports the rotating machine whose rotation speed changes. , The anomaly detection device obtains the envelope of the waveform in the time region, converts the waveform of the envelope in the time region into the intensity distribution data in the frequency region by high-speed Fourier conversion processing, and uses the intensity distribution data in the frequency region. The rotation speed calculation step for obtaining the rotation speed of the rotary machine, the rotation speed obtained by the abnormality detection device in the rotation speed calculation step, and the characteristic frequency indicating the cause of the vibration of the rolling bearing are used. It has an abnormality identification process for identifying an abnormality part of a rolling bearing. In the rotation speed calculation step, the abnormality detecting device obtains the vibration intensity at a plurality of frequencies proportional to the reference frequency, which is a frequency within a predetermined range, with respect to the intensity distribution data in the frequency domain, and obtains the obtained intensity. The rotation speed is obtained by using.

本発明によると、転がり軸受で支持されて回転速度が変化する回転機械の回転速度を算出できる、転がり軸受の異常検出装置と異常検出方法を提供することができる。 According to the present invention, it is possible to provide an abnormality detection device for a rolling bearing and an abnormality detection method capable of calculating the rotation speed of a rotating machine supported by a rolling bearing and whose rotation speed changes.

本発明の実施例1による、転がり軸受の異常検出装置を示す模式図である。It is a schematic diagram which shows the abnormality detection apparatus of a rolling bearing according to Example 1 of this invention. 転がり軸受の構成を示す模式図である。It is a schematic diagram which shows the structure of a rolling bearing. 転がり軸受を、回転機械の軸に直交する方向から見た断面の模式図である。It is a schematic diagram of the cross section of the rolling bearing seen from the direction orthogonal to the axis of the rotating machine. 本発明の実施例1による異常検出装置が行う、転がり軸受の異常部品を特定する処理の例を示すフローチャートである。It is a flowchart which shows the example of the process for identifying the abnormal part of a rolling bearing performed by the abnormality detection apparatus according to Example 1 of this invention. 図4のステップS5で異常検出装置が行う、回転機械の回転速度の算出の手順の一例を示すフローチャートである。It is a flowchart which shows an example of the procedure of calculation of the rotation speed of a rotating machine performed by an abnormality detection device in step S5 of FIG. 図4のステップS3で異常検出装置が行う、振動の時間領域の波形に包絡線処理を行って包絡線波形を求める処理の手順の一例を示すフローチャートである。It is a flowchart which shows an example of the procedure of the process of performing the envelope process on the waveform of the time domain of vibration, and obtaining the envelope waveform, which is performed by the abnormality detection apparatus in step S3 of FIG. 図4のステップS6で異常検出装置が行う、転がり軸受の異常部品を特定する処理の手順の一例を示すフローチャートである。It is a flowchart which shows an example of the procedure of the process of identifying the abnormal part of a rolling bearing performed by the abnormality detection apparatus in step S6 of FIG. 図4のステップS5で異常検出装置が行う、回転機械の回転速度の算出の手順の別の例を示すフローチャートである。It is a flowchart which shows another example of the procedure of calculation of the rotation speed of a rotating machine performed by an abnormality detection device in step S5 of FIG. 図4のステップS5で異常検出装置が行う、回転機械の回転速度の算出の手順のさらに別の例を示すフローチャートである。It is a flowchart which shows still another example of the procedure of calculation of the rotation speed of a rotating machine performed by an abnormality detection apparatus in step S5 of FIG. 転がり軸受の振動の時間領域の波形の例を示す図である。It is a figure which shows the example of the waveform of the time domain of the vibration of a rolling bearing. 転がり軸受の振動の時間領域の包絡線波形の例を示す図である。It is a figure which shows the example of the envelope waveform of the time domain of the vibration of a rolling bearing. 転がり軸受の振動の時間領域の包絡線波形の周波数(強度)分析の例を示す図である。It is a figure which shows the example of the frequency (intensity) analysis of the envelope waveform of the time domain of the vibration of a rolling bearing.

本発明による、転がり軸受の異常検出装置と異常検出方法は、転がり軸受で支持されて回転速度が変化する回転機械の回転速度を転がり軸受の振動を利用して算出し、算出した回転速度を用いて転がり軸受の異常部品を特定することで、転がり軸受の異常を検出する。本発明による、転がり軸受の異常検出装置と異常検出方法では、転がり軸受の異常検出に用いる回転速度を、次のようにして算出する。すなわち、転がり軸受の振動の波形に包絡線処理を行った後で高速フーリエ変換(以下、「FFT」とも記す)処理を行った周波数領域の強度分布データにおいて、基準周波数を設定し、この基準周波数に比例する複数の周波数における振動の強度を求め、求めた複数の周波数における振動の強度の和を求める。そして、複数の基準周波数についてこの振動の強度の和を求め、振動の強度の和が最大となる基準周波数を、回転機械の回転速度として算出する。求めた複数の周波数における振動の強度の平均が最大となる基準周波数を、回転機械の回転速度として算出することもでき、求めた複数の周波数における振動の強度の2乗平均(2乗平均平方根)が最大となる基準周波数を、回転機械の回転速度として算出することもできる。 In the abnormality detection device and abnormality detection method for rolling bearings according to the present invention, the rotation speed of a rotating machine supported by a rolling bearing and whose rotation speed changes is calculated by using the vibration of the rolling bearing, and the calculated rotation speed is used. By identifying the abnormal parts of the rolling bearing, the abnormality of the rolling bearing is detected. In the rolling bearing abnormality detection device and the abnormality detection method according to the present invention, the rotation speed used for the abnormality detection of the rolling bearing is calculated as follows. That is, a reference frequency is set in the intensity distribution data in the frequency domain in which the waveform of the vibration of the rolling bearing is subjected to the envelope processing and then the high-speed Fourier conversion (hereinafter, also referred to as “FFT”) processing, and this reference frequency is set. The intensity of vibration at a plurality of frequencies proportional to is obtained, and the sum of the vibration intensities at the obtained multiple frequencies is obtained. Then, the sum of the vibration intensities is obtained for a plurality of reference frequencies, and the reference frequency at which the sum of the vibration intensities is maximized is calculated as the rotation speed of the rotating machine. The reference frequency at which the average of the vibration intensities at the obtained multiple frequencies is maximized can be calculated as the rotation speed of the rotating machine, and the squared average of the vibration intensities at the obtained multiple frequencies (root mean square). The reference frequency at which is the maximum can also be calculated as the rotation speed of the rotating machine.

以下、本発明の実施例による、転がり軸受の異常検出装置と異常検出方法を、図面を用いて説明する。 Hereinafter, an abnormality detection device for rolling bearings and an abnormality detection method according to an embodiment of the present invention will be described with reference to the drawings.

図1は、本発明の実施例1による、転がり軸受の異常検出装置10を示す模式図である。図1には、異常検出装置10が異常を検出する転がり軸受20と、転がり軸受20が支持する回転機械30も示している。 FIG. 1 is a schematic view showing an abnormality detection device 10 for rolling bearings according to the first embodiment of the present invention. FIG. 1 also shows a rolling bearing 20 in which the abnormality detecting device 10 detects an abnormality, and a rotary machine 30 supported by the rolling bearing 20.

異常検出装置10は、データ取得部11と、回転速度算出部12と、異常特定部13と、出力部14を備え、データロガー16に接続される。また、異常検出装置10は、転がり軸受20に設置された振動検出器15と、不要な振動成分を除去するためのバンドパスフィルタ18(「BPF」とも記す)を備える。 The abnormality detection device 10 includes a data acquisition unit 11, a rotation speed calculation unit 12, an abnormality identification unit 13, and an output unit 14, and is connected to the data logger 16. Further, the abnormality detection device 10 includes a vibration detector 15 installed in the rolling bearing 20 and a bandpass filter 18 (also referred to as “BPF”) for removing unnecessary vibration components.

データ取得部11は、振動検出器15が検出した振動の波形を取得する。回転速度算出部12は、データ取得部11が取得した波形から回転機械30の回転速度を求める。異常特定部13は、回転速度算出部12が求めた回転速度を用いて、転がり軸受20の異常部品を特定する。出力部14は、表示画面や、外部装置へのインターフェースで構成することができ、異常検出の結果を出力する。 The data acquisition unit 11 acquires the vibration waveform detected by the vibration detector 15. The rotation speed calculation unit 12 obtains the rotation speed of the rotating machine 30 from the waveform acquired by the data acquisition unit 11. The abnormality specifying unit 13 identifies an abnormal component of the rolling bearing 20 by using the rotation speed obtained by the rotation speed calculation unit 12. The output unit 14 can be configured with a display screen or an interface to an external device, and outputs an abnormality detection result.

異常検出装置10は、パソコンなどのコンピュータ17で操作することができる装置である。異常検出装置10は、コンピュータ17で構成することもできる。 The abnormality detection device 10 is a device that can be operated by a computer 17 such as a personal computer. The abnormality detection device 10 can also be configured by a computer 17.

振動検出器15は、例えば加速度センサで構成されており、転がり軸受20の振動を検出する。 The vibration detector 15 is composed of, for example, an acceleration sensor, and detects the vibration of the rolling bearing 20.

データロガー16は、振動検出器15が検出した転がり軸受20の振動のデータ、すなわち、転がり軸受20の振動の時間領域の波形のデータを保存する保存装置である。 The data logger 16 is a storage device that stores the vibration data of the rolling bearing 20 detected by the vibration detector 15, that is, the waveform data of the vibration of the rolling bearing 20 in the time domain.

異常検出装置10のデータ取得部11は、データロガー16から、転がり軸受20の振動の時間領域の波形を取得する。 The data acquisition unit 11 of the abnormality detection device 10 acquires the waveform of the vibration of the rolling bearing 20 in the time domain from the data logger 16.

転がり軸受20は、回転機械30を支持する。 The rolling bearing 20 supports the rotary machine 30.

回転機械30は、回転する軸31を備え、軸31が転がり軸受20で支持されている。回転機械30は、運転中に回転速度(すなわち、軸31の回転速度)が変化する。回転機械30の例は、電動機、遠心圧縮機、及びポンプである。 The rotary machine 30 includes a rotating shaft 31, and the shaft 31 is supported by a rolling bearing 20. The rotation speed of the rotary machine 30 (that is, the rotation speed of the shaft 31) changes during operation. Examples of the rotary machine 30 are a motor, a centrifugal compressor, and a pump.

図2は、転がり軸受20の構成を示す模式図であり、転がり軸受20を、回転機械30の軸31に平行な方向から見た模式図である。転がり軸受20は、内輪21と、外輪22と、複数の転動体23と、保持器24を備える。内輪21は、回転機械30の軸31に固定されている環状部材であり、軸31とともに回転する。外輪22は、転がり軸受20のハウジングに固定されている環状部材であり、内輪21と同心円状に配置されている。複数の転動体23は、内輪21と外輪22との間の空間に配置された部材であり、自転しながら内輪21と同じ周方向に回転(公転)する。保持器24は、複数の転動体23を、互いの周方向の相対位置を保つように保持する部材である。 FIG. 2 is a schematic view showing the configuration of the rolling bearing 20, and is a schematic view of the rolling bearing 20 as viewed from a direction parallel to the axis 31 of the rotary machine 30. The rolling bearing 20 includes an inner ring 21, an outer ring 22, a plurality of rolling elements 23, and a cage 24. The inner ring 21 is an annular member fixed to the shaft 31 of the rotating machine 30 and rotates together with the shaft 31. The outer ring 22 is an annular member fixed to the housing of the rolling bearing 20, and is arranged concentrically with the inner ring 21. The plurality of rolling elements 23 are members arranged in the space between the inner ring 21 and the outer ring 22, and rotate (revolve) in the same circumferential direction as the inner ring 21 while rotating. The cage 24 is a member that holds a plurality of rolling elements 23 so as to maintain relative positions in the circumferential direction with each other.

図2には、転動体23の径dと、ピッチ円径D(転動体23の中心を通る円の直径)も示している。 FIG. 2 also shows the diameter d of the rolling element 23 and the pitch circle diameter D (the diameter of the circle passing through the center of the rolling element 23).

図3は、転がり軸受20を、回転機械30の軸31に直交する方向から見た断面の模式図である。図3には、転がり軸受20の接触角αを示している。接触角αは、転動体23の転走面(転動体23と接する内輪21と外輪22の面)と転動体23との間にかかる荷重の方向25と、転がり軸受20の径方向26(軸31に直交する方向)との間の角度である。 FIG. 3 is a schematic cross-sectional view of the rolling bearing 20 as viewed from a direction orthogonal to the axis 31 of the rotary machine 30. FIG. 3 shows the contact angle α of the rolling bearing 20. The contact angle α is the direction 25 of the load applied between the rolling surface of the rolling element 23 (the surface of the inner ring 21 and the outer ring 22 in contact with the rolling element 23) and the rolling element 23, and the radial direction 26 (axis) of the rolling bearing 20. The angle between the direction perpendicular to 31).

図4は、本発明の実施例1による異常検出装置10が行う、転がり軸受20の異常部品を特定する処理の例を示すフローチャートである。 FIG. 4 is a flowchart showing an example of a process for identifying an abnormal component of the rolling bearing 20 performed by the abnormality detecting device 10 according to the first embodiment of the present invention.

ステップS1で、異常検出装置10が起動され、処理が開始される。 In step S1, the abnormality detection device 10 is started and the process is started.

ステップS2で、異常検出装置10のデータ取得部11は、データロガー16から、転がり軸受20の振動の時間領域の波形を取得する。 In step S2, the data acquisition unit 11 of the abnormality detection device 10 acquires the waveform of the vibration of the rolling bearing 20 in the time domain from the data logger 16.

図10は、ステップS2で取得した、転がり軸受20の振動の時間領域の波形の例を示す図である。図10に示すように、転がり軸受20には、回転機械30の回転(転がり軸受20の内輪21の回転)に応じて、周期的な振動(加速度の変化)が発生している。 FIG. 10 is a diagram showing an example of the waveform in the time domain of the vibration of the rolling bearing 20 acquired in step S2. As shown in FIG. 10, in the rolling bearing 20, periodic vibration (change in acceleration) is generated in accordance with the rotation of the rotating machine 30 (rotation of the inner ring 21 of the rolling bearing 20).

ステップS3で、異常検出装置10の回転速度算出部12は、ステップS2で取得した振動の波形に包絡線処理を行い、振動の時間領域の波形の包絡線を求める。すなわち、異常検出装置10は、振動の時間領域の波形の輪郭をなぞる曲線(包絡線)を求める。 In step S3, the rotation speed calculation unit 12 of the abnormality detection device 10 performs envelope processing on the vibration waveform acquired in step S2, and obtains the envelope of the waveform in the vibration time domain. That is, the abnormality detecting device 10 obtains a curve (envelope) that traces the contour of the waveform in the time domain of vibration.

図11は、ステップS3で求めた、転がり軸受20の振動の時間領域の包絡線波形の例を示す図である。図11は、転がり軸受20に発生した周期的な振動(加速度の変化)の波形の包絡線を示している。 FIG. 11 is a diagram showing an example of the envelope waveform in the time domain of the vibration of the rolling bearing 20 obtained in step S3. FIG. 11 shows the envelope of the waveform of the periodic vibration (change in acceleration) generated in the rolling bearing 20.

ステップS4で、異常検出装置10の回転速度算出部12は、ステップS3で求めた時間領域の包絡線波形に高速フーリエ変換処理(FFT処理)を行い、時間領域の包絡線波形を周波数領域の強度分布データに変換する。 In step S4, the rotation speed calculation unit 12 of the abnormality detection device 10 performs a fast Fourier transform process (FFT process) on the time domain envelope waveform obtained in step S3, and converts the time domain envelope waveform into the frequency domain intensity. Convert to distribution data.

図12は、ステップS4で求めた、転がり軸受20の振動の時間領域の包絡線波形の周波数(強度)分析の例を示す図である。図12に示すように、振動の時間領域の包絡線波形の周波数(強度)分析では、いくつかの周波数において、振動の強度が大きくなっている(すなわち、ピークが発生している)ことがわかる。 FIG. 12 is a diagram showing an example of frequency (intensity) analysis of the envelope waveform in the time domain of the vibration of the rolling bearing 20 obtained in step S4. As shown in FIG. 12, the frequency (intensity) analysis of the envelope waveform in the time domain of vibration shows that the vibration intensity is high (that is, a peak is generated) at some frequencies. ..

ステップS5で、異常検出装置10の回転速度算出部12は、ステップS4で得た時間領域の包絡線波形の周波数(強度)分析を用いて、回転機械30の回転速度を算出する。 In step S5, the rotation speed calculation unit 12 of the abnormality detecting device 10 calculates the rotation speed of the rotating machine 30 by using the frequency (intensity) analysis of the envelope waveform in the time domain obtained in step S4.

ステップS6で、異常検出装置10の異常特定部13は、ステップS5で算出した回転機械30の回転速度と、転がり軸受20の特徴周波数を用いて、転がり軸受20の異常部品を特定する。 In step S6, the abnormality specifying unit 13 of the abnormality detecting device 10 identifies an abnormal component of the rolling bearing 20 by using the rotation speed of the rotating machine 30 calculated in step S5 and the characteristic frequency of the rolling bearing 20.

ステップS7で、異常検出装置10は、転がり軸受20の異常検出を終了して停止する。 In step S7, the abnormality detection device 10 ends the abnormality detection of the rolling bearing 20 and stops.

転がり軸受20の特徴周波数とは、転がり軸受20の幾何学的寸法と回転機械30の回転速度(回転周波数)から決まる値であり、転がり軸受20の異常な振動の要因を示す周波数である。特徴周波数は、例えば、転動体23の公転周波数(保持器24の回転周波数)FTF、転動体23の自転周波数BSF、外輪転動体通過周波数BPFO、及び内輪転動体通過周波数BPFIである。これらの特徴周波数は、以下の式で計算することができる。
転動体23の公転周波数(保持器24の回転周波数)FTF
FTF=1/2*(1-d/D*cosα)*fs (1)
転動体23の自転周波数BSF
BSF=D/(2d)*[1-(d/D*cosα)]*fs (2)
外輪転動体通過周波数(外輪22の一地点を転動体23が通過していく周波数)BPFO
BPFO=z/2*(1-d/D*cosα)*fs (3)
内輪転動体通過周波数(内輪21の一地点を転動体23が通過していく周波数)BPFI
BPFI=z/2*(1+d/D*cosα)*fs (4)
但し、式(1)~(4)において、
dは転動体23の径(mm)、
Dはピッチ円径(mm)
αは転がり軸受20の接触角(ラジアン)
zは転動体23の数
fsは回転機械30の回転周波数(Hz)
である。
The characteristic frequency of the rolling bearing 20 is a value determined by the geometrical dimensions of the rolling bearing 20 and the rotation speed (rotational frequency) of the rotating machine 30, and is a frequency indicating a factor of abnormal vibration of the rolling bearing 20. The characteristic frequencies are, for example, the revolution frequency of the rolling element 23 (rotational frequency of the cage 24) FTF, the rotation frequency BSF of the rolling element 23, the outer ring rolling element passing frequency BPFO, and the inner ring rolling element passing frequency BPFI. These characteristic frequencies can be calculated by the following formula.
Revolution frequency of rolling element 23 (rotation frequency of cage 24) FTF
FTF = 1/2 * (1-d / D * cosα) * fs (1)
Rotation frequency BSF of rolling element 23
BSF = D / (2d) * [1- (d / D * cosα) 2 ] * fs (2)
Outer ring rolling element passing frequency (frequency at which the rolling element 23 passes through one point on the outer ring 22) BPFO
BPFO = z / 2 * (1-d / D * cosα) * fs (3)
Inner ring rolling element passing frequency (frequency at which the rolling element 23 passes through one point of the inner ring 21) BPFI
BPFI = z / 2 * (1 + d / D * cosα) * fs (4)
However, in the formulas (1) to (4),
d is the diameter (mm) of the rolling element 23,
D is the pitch circle diameter (mm)
α is the contact angle (radian) of the rolling bearing 20.
z is the number fs of the rolling elements 23 and is the rotation frequency (Hz) of the rotating machine 30.
Is.

転がり軸受20に損傷や欠陥などの異常があると、転がり軸受20の振動の周波数領域の波形には、異常の部位に関連した特徴周波数とその高調波の周波数にピークが発生する。すなわち、保持器24に異常がある場合は、式(1)で表される周波数とその高調波の周波数にピークが発生する。転動体23に異常がある場合は、式(2)の2倍の周波数とその高調波の周波数にピークが発生する。外輪22に異常がある場合は、式(3)の周波数とその高調波の周波数にピークが発生する。内輪21に異常がある場合は、式(4)の周波数とその高調波の周波数にピークが発生する。 When the rolling bearing 20 has an abnormality such as damage or a defect, a peak occurs in the characteristic frequency related to the abnormal portion and the frequency of its harmonics in the waveform of the vibration frequency region of the rolling bearing 20. That is, when there is an abnormality in the cage 24, a peak occurs in the frequency represented by the equation (1) and the frequency of its harmonics. If there is an abnormality in the rolling element 23, a peak occurs at twice the frequency of Eq. (2) and the frequency of its harmonics. When the outer ring 22 has an abnormality, a peak occurs in the frequency of the equation (3) and the frequency of its harmonics. When the inner ring 21 has an abnormality, a peak occurs in the frequency of the equation (4) and the frequency of its harmonics.

なお、これらのピークは、転がり軸受20に損傷や欠陥などの異常がほとんど無い場合でも、微弱なものを確認できることが多い。このため、これらのピークを与える周波数(すなわち、特徴周波数とその高調波の周波数)は、後述するように、回転機械30の回転速度の算出に用いることができる。 It should be noted that these peaks can often be confirmed to be weak even when there are almost no abnormalities such as damage or defects in the rolling bearing 20. Therefore, the frequency giving these peaks (that is, the characteristic frequency and the frequency of its harmonics) can be used for calculating the rotation speed of the rotating machine 30 as described later.

図5は、図4のステップS5で異常検出装置10の回転速度算出部12が行う、回転機械30の回転速度の算出の手順の一例を示すフローチャートである。転がり軸受20の特徴周波数を用いて異常部品を特定するには、式(1)~(4)に示したように、回転機械30の回転速度(回転周波数fs)が既知であることが前提である。本実施例では、ステップS4で得た周波数領域の強度分布データについて、基準周波数を予め定めた範囲内の周波数として設定し、基準周波数に比例する複数の周波数における振動の強度を求め、求めた複数の周波数における振動の強度の和を求める。複数の基準周波数についてこの振動の強度の和を求め、振動の強度の和が最大となる基準周波数を、回転機械30の回転速度として求める。 FIG. 5 is a flowchart showing an example of a procedure for calculating the rotation speed of the rotary machine 30 performed by the rotation speed calculation unit 12 of the abnormality detection device 10 in step S5 of FIG. In order to identify abnormal parts using the characteristic frequency of the rolling bearing 20, it is premised that the rotation speed (rotation frequency fs) of the rotary machine 30 is known as shown in the equations (1) to (4). be. In this embodiment, for the intensity distribution data in the frequency domain obtained in step S4, the reference frequency is set as a frequency within a predetermined range, and the vibration intensities at a plurality of frequencies proportional to the reference frequency are obtained and obtained. Find the sum of the vibration intensities at the frequency of. The sum of the vibration intensities is obtained for a plurality of reference frequencies, and the reference frequency at which the sum of the vibration intensities is maximized is obtained as the rotation speed of the rotating machine 30.

ステップS501で、異常検出装置10は、使用する変数であるWmax、SPD、及びispに初期値を設定する。Wmaxは、基準周波数に比例する複数の周波数における振動の強度の和の最大値を表す。SPDは、回転機械30の回転速度を表す。ispは、基準周波数に比例する複数の周波数における振動の強度の和が最大となるときのループカウンタiの値である。本実施例では、一例として、Wmax=0.00、SPD=0、isp=0と設定した。 In step S501, the abnormality detection device 10 sets initial values for the variables Wmax, SPD, and ISP to be used. Wmax represents the maximum value of the sum of the vibration intensities at a plurality of frequencies proportional to the reference frequency. SPD represents the rotation speed of the rotating machine 30. isp is the value of the loop counter i when the sum of the vibration intensities at a plurality of frequencies proportional to the reference frequency is maximized. In this embodiment, as an example, Wmax = 0.00, SPD = 0, and ISP = 0 are set.

ステップS502で、異常検出装置10は、基準周波数fs(i)の範囲を設定する。基準周波数fs(i)の最小値をf1とし、基準周波数fs(i)の最大値をfNとすると、基準周波数fs(i)の範囲は、f1からfNである。基準周波数fs(i)の範囲(最小値f1と最大値fN)は、例えば、回転機械30の運転周波数の範囲などに基づいて予め定めることができる。基準周波数fs(i)の範囲は、異常検出装置10の操作者が最小値f1と最大値fNを異常検出装置10に入力することで、異常検出装置10に設定される。 In step S502, the abnormality detection device 10 sets the range of the reference frequency fs (i). Assuming that the minimum value of the reference frequency fs (i) is f1 and the maximum value of the reference frequency fs (i) is fN, the range of the reference frequency fs (i) is f1 to fN. The range of the reference frequency fs (i) (minimum value f1 and maximum value fN) can be predetermined based on, for example, the range of the operating frequency of the rotary machine 30. The range of the reference frequency fs (i) is set in the abnormality detection device 10 by the operator of the abnormality detection device 10 inputting the minimum value f1 and the maximum value fN into the abnormality detection device 10.

ステップS503で、異常検出装置10は、ループカウンタiの値を1に設定する。 In step S503, the abnormality detection device 10 sets the value of the loop counter i to 1.

ステップS504で、異常検出装置10は、基準周波数fs(i)にf1を設定する。 In step S504, the abnormality detection device 10 sets f1 to the reference frequency fs (i).

ステップS505で、異常検出装置10は、基準周波数fs(i)を変化させる周波数刻みΔfを設定する。周波数刻みΔfは、基準周波数fs(i)がステップS502で設定した範囲内で変化するように任意に定めることができ、例えば、得られる回転速度の精度と計算に要する時間を考慮して定めることができる。周波数刻みΔfは、異常検出装置10の操作者が異常検出装置10に入力することで、異常検出装置10に設定される。 In step S505, the abnormality detection device 10 sets the frequency step Δf that changes the reference frequency fs (i). The frequency step Δf can be arbitrarily set so that the reference frequency fs (i) changes within the range set in step S502. For example, the frequency step Δf is set in consideration of the accuracy of the obtained rotation speed and the time required for calculation. Can be done. The frequency step Δf is set in the abnormality detection device 10 by being input to the abnormality detection device 10 by the operator of the abnormality detection device 10.

基準周波数fs(i)は、転がり軸受20の特徴周波数とその高調波の周波数を用いて定めることもできる。すなわち、転がり軸受20の異常部位が予想できる場合は、基準周波数fs(i)を、その異常部位に関連する特徴周波数から決められる周波数(その高調波も含む)としてもよい。転がり軸受20の異常部位が予想できない場合は、基準周波数fs(i)を、転がり軸受20の全ての特徴周波数から決められる周波数(その高調波も含む)としてもよい。 The reference frequency fs (i) can also be determined by using the characteristic frequency of the rolling bearing 20 and the frequency of its harmonics. That is, if an abnormal portion of the rolling bearing 20 can be predicted, the reference frequency fs (i) may be set to a frequency (including its harmonics) determined from the characteristic frequency related to the abnormal portion. If the abnormal portion of the rolling bearing 20 cannot be predicted, the reference frequency fs (i) may be set to a frequency (including its harmonics) determined from all the characteristic frequencies of the rolling bearing 20.

ステップS506で、異常検出装置10は、周波数領域の強度分布データについて、基準周波数fs(i)に比例する複数の周波数における振動の強度を求める。このとき、基準周波数fs(i)に比例する複数の周波数は、Δf×比例倍数の範囲で最大となる周波数とする。また、異常(損傷)部位があらかじめ特定できる機械の場合は、特定の損傷周波数(特徴周波数から決まる周波数)としてもよい。異常(損傷)部位があらかじめ特定できない機械の場合は、すべての損傷周波数(特徴周波数から決まる周波数)としてもよい。これらの周波数は、その高周波を含めてもよい。振動の強度は、例えば、振幅の大きさから求めることができる。振動の強度を求める周波数の数は、例えば、特徴周波数の数を基に定めることができる。本実施例では、特徴周波数の数が4個(式(1)~(4))であるので、振動の強度を求める周波数の数は、最低4個であり、n次の高調波まで考慮するのであれば4n個とすることができる(nの最大値は、例えば4である)。 In step S506, the abnormality detection device 10 obtains the vibration intensity at a plurality of frequencies proportional to the reference frequency fs (i) for the intensity distribution data in the frequency domain. At this time, the plurality of frequencies proportional to the reference frequency fs (i) are set to be the maximum frequencies in the range of Δf × proportional multiple. Further, in the case of a machine in which an abnormal (damaged) part can be specified in advance, a specific damage frequency (frequency determined from the characteristic frequency) may be used. In the case of a machine in which the abnormal (damaged) part cannot be identified in advance, all damaged frequencies (frequency determined from the characteristic frequency) may be used. These frequencies may include the high frequency. The vibration intensity can be obtained from, for example, the magnitude of the amplitude. The number of frequencies for which the vibration intensity is obtained can be determined, for example, based on the number of characteristic frequencies. In this embodiment, since the number of characteristic frequencies is four (equations (1) to (4)), the number of frequencies for which the vibration intensity is obtained is at least four, and even the nth harmonic is considered. If, the number can be 4n (the maximum value of n is, for example, 4).

ステップS507で、異常検出装置10は、ステップS506で求めた複数の周波数における振動の強度の和W(i)を求める。 In step S507, the abnormality detection device 10 obtains the sum W (i) of the vibration intensities at the plurality of frequencies obtained in step S506.

ステップS508で、異常検出装置10は、ステップS507で求めた振動の強度の和W(i)とWmaxの大きさを比較する。振動の強度の和W(i)がWmaxより大きい場合は、ステップS509に進み、W(i)がWmax以下である場合は、ステップS510に進む。 In step S508, the abnormality detection device 10 compares the sum W (i) and the magnitude of Wmax of the vibration intensities obtained in step S507. If the sum W (i) of the vibration intensities is larger than Wmax, the process proceeds to step S509, and if W (i) is Wmax or less, the process proceeds to step S510.

ステップS509で、異常検出装置10は、WmaxにW(i)の値を代入し、ispにiを設定した後、ステップS510に進む。 In step S509, the abnormality detection device 10 substitutes the value of W (i) for Wmax, sets i for ISP, and then proceeds to step S510.

ステップS510で、異常検出装置10は、ループカウンタiの値に1を加える。 In step S510, the abnormality detection device 10 adds 1 to the value of the loop counter i.

ステップS511で、異常検出装置10は、基準周波数fs(i)にfs(i-1)+Δfを設定する。すなわち、基準周波数fs(i)は、Δfだけ増加するように変化する。 In step S511, the abnormality detection device 10 sets the reference frequency fs (i) to fs (i-1) + Δf. That is, the reference frequency fs (i) changes so as to increase by Δf.

ステップS512で、異常検出装置10は、基準周波数fs(i)が基準周波数fs(i)の範囲の最大値fN以下であるか判定する。基準周波数fs(i)が最大値fN以下の場合は、異常検出装置10は、ステップS506からステップS511までの処理を繰り返す。ステップS511で基準周波数fs(i)をΔfだけ変化させているので、ステップS506では、異常検出装置10は、Δfだけ変化させた基準周波数fs(i)に比例する複数の周波数における振動の強度を求めることになる。基準周波数fs(i)が最大値fNより大きい場合は、ステップS513に進む。 In step S512, the abnormality detection device 10 determines whether the reference frequency fs (i) is equal to or less than the maximum value fN in the range of the reference frequency fs (i). When the reference frequency fs (i) is equal to or less than the maximum value fN, the abnormality detection device 10 repeats the processes from step S506 to step S511. Since the reference frequency fs (i) is changed by Δf in step S511, in step S506, the abnormality detection device 10 determines the vibration intensity at a plurality of frequencies proportional to the reference frequency fs (i) changed by Δf. You will be asked. If the reference frequency fs (i) is larger than the maximum value fN, the process proceeds to step S513.

ステップS513で、異常検出装置10は、ispの値がゼロか判定する。ispの値がゼロの場合は、ステップS505からステップS512までの処理を繰り返す。ispの値がゼロでない場合は、ステップS514に進む。 In step S513, the abnormality detection device 10 determines whether the value of ISP is zero. If the value of ISP is zero, the processes from step S505 to step S512 are repeated. If the value of ISP is not zero, the process proceeds to step S514.

ステップS514で、異常検出装置10は、基準周波数fs(isp)を回転機械30の回転速度SPDとする。回転機械30の回転速度SPDは、このようにして算出することで求められる。 In step S514, the abnormality detection device 10 sets the reference frequency fs (isp) to the rotation speed SPD of the rotary machine 30. The rotation speed SPD of the rotary machine 30 is obtained by calculating in this way.

この回転速度SPDの算出には、基準周波数fs(i)に比例する周波数以外の周波数における振動の強度は、通常、大きくならないことを利用している。転がり軸受20の振動は、回転機械30の回転速度SPDに比例した周波数で発生し、この周波数以外の周波数ではほとんど発生せず、発生したとしても小さな振動である(なお、比例の倍数は、振動の原因により異なる)。従って、基準周波数fs(i)が回転速度SPDと一致した場合には、基準周波数fs(i)に比例する周波数における振動の強度が大きくなる。但し、どの原因(軸受の場合は、損傷や欠陥の位置)によって振動が発生するのか不明であるので、振動を起こす可能性のある複数の周波数における振動の強度の和W(i)を求めて、和W(i)が最大となる基準周波数fs(isp)を回転機械30の回転速度SPDとして求めることにしている。 The calculation of the rotation speed SPD utilizes the fact that the vibration intensity at a frequency other than the frequency proportional to the reference frequency fs (i) usually does not increase. The vibration of the rolling bearing 20 occurs at a frequency proportional to the rotation speed SPD of the rotating machine 30, and hardly occurs at frequencies other than this frequency, and even if it occurs, it is a small vibration (note that a multiple of the proportional vibration is vibration). Depends on the cause of). Therefore, when the reference frequency fs (i) coincides with the rotation speed SPD, the vibration intensity at a frequency proportional to the reference frequency fs (i) increases. However, since it is unclear which cause (in the case of bearings, the position of damage or defect) causes vibration, the sum of the vibration intensities W (i) at multiple frequencies that may cause vibration is obtained. , The reference frequency fs (isp) at which the sum W (i) is maximized is determined as the rotation speed SPD of the rotating machine 30.

図6は、図4のステップS3で異常検出装置10の回転速度算出部12が行う、振動の時間領域の波形に包絡線処理を行って包絡線波形を求める処理の手順の一例を示すフローチャートである。 FIG. 6 is a flowchart showing an example of a procedure of processing for obtaining an envelope waveform by performing an envelope process on a waveform in a vibration time domain, which is performed by the rotation speed calculation unit 12 of the abnormality detection device 10 in step S3 of FIG. be.

転がり軸受20に異常があると、振動の時間領域の波形として、高周波数の振動が周期的に変動した成分を含む波形が得られる。この時間領域の波形に対し、単に高速フーリエ変換処理(FFT処理)を行うと、高い周波数成分のスペクトルが得られるのみで、高周波数の振動が周期的に変動するスペクトルは得られない。このため、周期的に変動する高周波数の振動の周期性を抽出するためには、振動の波形に包絡線処理を行ってからFFT処理を行うのが有効である。 If there is an abnormality in the rolling bearing 20, a waveform including a component in which high-frequency vibration periodically fluctuates can be obtained as a waveform in the time domain of vibration. When the fast Fourier transform process (FFT process) is simply performed on the waveform in this time domain, only the spectrum of the high frequency component can be obtained, and the spectrum in which the high frequency vibration fluctuates periodically cannot be obtained. Therefore, in order to extract the periodicity of high-frequency vibration that fluctuates periodically, it is effective to perform envelope processing on the vibration waveform and then perform FFT processing.

ステップS301で、異常検出装置10は、振動の時間領域の波形(図10)にFFT処理を行い、振動波形のデータを時間領域から周波数領域に変換する。 In step S301, the abnormality detection device 10 performs FFT processing on the waveform (FIG. 10) in the time domain of vibration, and converts the vibration waveform data from the time domain to the frequency domain.

ステップS302で、異常検出装置10は、バンドパスフィルタ18(「BPF」とも記す)を用いて、周波数領域に変換した振動の波形から不要な振動成分を除去する。回転1次成分や低周波の振動成分は、高周波数の振動の波形を抽出する際にノイズとなる。このため、バンドパスフィルタ18を用いて、これらの不要な振動成分を除去し、必要な周波数範囲の振動成分を取り出す。バンドパスフィルタ18が取り出す周波数成分は、転がり軸受20の構成や、回転機械30の運転周波数の範囲に基づいて定めることができる。例えば、損傷した転がり軸受20の構成部品は、その固有振動の周波数で振動するので、バンドパスフィルタ18は、この固有振動の周波数成分を取り出せるように設定される。 In step S302, the abnormality detection device 10 uses a bandpass filter 18 (also referred to as “BPF”) to remove unnecessary vibration components from the vibration waveform converted into the frequency domain. The rotation primary component and the low frequency vibration component become noise when extracting the waveform of the high frequency vibration. Therefore, the bandpass filter 18 is used to remove these unnecessary vibration components and take out the vibration components in the required frequency range. The frequency component taken out by the bandpass filter 18 can be determined based on the configuration of the rolling bearing 20 and the range of the operating frequency of the rotary machine 30. For example, since the component of the damaged rolling bearing 20 vibrates at the frequency of its natural vibration, the bandpass filter 18 is set so that the frequency component of this natural vibration can be taken out.

ステップS303で、異常検出装置10は、バンドパスフィルタ18で不要な振動成分を除去した振動の波形(周波数領域の波形)に対し、逆高速フーリエ変換処理(IFFT処理)を行い、振動波形のデータを周波数領域から時間領域に戻す。 In step S303, the abnormality detection device 10 performs inverse high-speed Fourier conversion processing (IFFT processing) on the vibration waveform (waveform in the frequency domain) from which unnecessary vibration components have been removed by the bandpass filter 18, and the vibration waveform data. Is returned from the frequency domain to the time domain.

ステップS304で、異常検出装置10は、時間領域に戻した振動の波形に包絡線処理を行い、振動の時間領域の包絡線波形(図11)を得る。包絡線処理は、振動の時間領域の波形の輪郭をなぞる曲線(包絡線)を求める処理である。包絡線処理には、波形の信号にヒルベルト変換を施す方法など、任意の既存の方法を用いることができる。 In step S304, the abnormality detection device 10 performs envelope processing on the vibration waveform returned to the time domain, and obtains an envelope waveform (FIG. 11) in the vibration time domain. The envelope processing is a process for obtaining a curve (envelope) that traces the contour of the waveform in the time domain of vibration. For the envelope processing, any existing method such as a method of applying a Hilbert transform to a waveform signal can be used.

以上のようにして振動の波形に包絡線処理を行うと、周期的に変動する高周波数の振動の周期性を抽出することが可能になるため、転がり軸受20の異常の特定に用いる特徴周波数を抽出することができる。 By performing envelope processing on the vibration waveform as described above, it is possible to extract the periodicity of high-frequency vibration that fluctuates periodically. Therefore, the characteristic frequency used to identify the abnormality of the rolling bearing 20 is selected. Can be extracted.

図7は、図4のステップS6で異常検出装置10が行う、転がり軸受20の異常部品を特定する処理の手順の一例を示すフローチャートである。 FIG. 7 is a flowchart showing an example of a procedure for identifying an abnormal component of the rolling bearing 20 performed by the abnormality detecting device 10 in step S6 of FIG.

ステップS601で、異常検出装置10の異常特定部13は、図4のステップS4で得られた、振動の周波数領域の包絡線波形に、予め任意に定めた閾値より大きい振幅が含まれているかを調べる。振動の周波数領域の強度分布データに、この閾値より大きい振幅が含まれている場合は、転がり軸受20に異常があるとして、ステップS602に進む。振動の周波数領域の強度分布データに、この閾値より大きい振幅が含まれていない場合は、転がり軸受20に異常がないとして、ステップS603に進む。 In step S601, the abnormality identification unit 13 of the abnormality detection device 10 determines whether the envelope waveform in the vibration frequency region obtained in step S4 of FIG. 4 contains an amplitude larger than a predetermined threshold value. investigate. If the intensity distribution data in the vibration frequency region contains an amplitude larger than this threshold value, it is considered that the rolling bearing 20 has an abnormality, and the process proceeds to step S602. If the intensity distribution data in the vibration frequency region does not include an amplitude larger than this threshold value, it is assumed that the rolling bearing 20 has no abnormality, and the process proceeds to step S603.

ステップS602で、異常検出装置10の異常特定部13は、特徴周波数を用いて転がり軸受20の異常部品を特定する。特徴周波数は、図4のステップS5で算出した回転機械30の回転速度SPDを式(1)~(4)の回転周波数fsとすることで、式(1)~(4)から求めることができる。異常検出装置10は、ステップS601で求めた、閾値より大きい振幅について、この振幅のピークを与える周波数を特定する。そして、異常検出装置10は、ピークを与える周波数が、式(1)~(4)で与えられる特徴周波数やその高調波の周波数のうち、どの周波数に対応するか調べる。異常検出装置10は、ピークを与える周波数に対応した特徴周波数(特徴周波数の高調波も含む)から、転がり軸受20の異常部品が保持器24、転動体23、外輪22、及び内輪21のうちのどれであるかを特定する。 In step S602, the abnormality specifying unit 13 of the abnormality detecting device 10 identifies the abnormal component of the rolling bearing 20 by using the characteristic frequency. The characteristic frequency can be obtained from the equations (1) to (4) by setting the rotational speed SPD of the rotary machine 30 calculated in step S5 of FIG. 4 to the rotational frequency fs of the equations (1) to (4). .. The anomaly detection device 10 specifies a frequency that gives a peak of this amplitude with respect to an amplitude larger than the threshold value obtained in step S601. Then, the abnormality detecting device 10 examines which frequency corresponds to the characteristic frequency given by the equations (1) to (4) and the frequency of the harmonic thereof. In the abnormality detection device 10, the abnormal parts of the rolling bearing 20 are included in the cage 24, the rolling element 23, the outer ring 22, and the inner ring 21 from the characteristic frequency (including the harmonic of the characteristic frequency) corresponding to the frequency giving the peak. Identify which one.

ステップS603で、異常検出装置10の出力部14は、転がり軸受20に異常があるかないかという検出結果を出力する。異常検出装置10は、転がり軸受20に異常がある場合には、ステップS602で特定した異常部品が何かを出力することもできる。異常検出装置10は、検出結果や特定した異常部品を、異常検出装置10の表示画面に出力したり、コンピュータ17や外部の表示装置、記憶装置、または印刷装置に出力したりすることができる。 In step S603, the output unit 14 of the abnormality detection device 10 outputs a detection result as to whether or not there is an abnormality in the rolling bearing 20. When the rolling bearing 20 has an abnormality, the abnormality detecting device 10 can also output what the abnormal component identified in step S602 is. The abnormality detection device 10 can output the detection result and the identified abnormal component to the display screen of the abnormality detection device 10, or output to the computer 17, an external display device, a storage device, or a printing device.

図5に示したフローチャートでは、ステップS507において、異常検出装置10は、ステップS506で求めた複数の周波数における振動の強度の和W(i)を求め、ステップS514において、この和W(i)が最大となる基準周波数fs(isp)を回転機械30の回転速度SPDとして算出する。本実施例による異常検出装置10は、複数の周波数における振動の強度の和W(i)の他に、振動の強度の平均W1(i)や振動の強度の2乗平均W2(i)を用いて、回転機械30の回転速度SPDを算出することもできる。 In the flowchart shown in FIG. 5, in step S507, the abnormality detecting device 10 obtains the sum W (i) of the vibration intensities at the plurality of frequencies obtained in step S506, and in step S514, the sum W (i) is calculated. The maximum reference frequency fs (isp) is calculated as the rotation speed SPD of the rotating machine 30. The abnormality detecting device 10 according to the present embodiment uses the average W1 (i) of the vibration intensity and the squared average W2 (i) of the vibration intensity in addition to the sum W (i) of the vibration intensities at a plurality of frequencies. It is also possible to calculate the rotation speed SPD of the rotating machine 30.

上述したように、基準周波数fs(i)が回転速度SPDと一致した場合には、基準周波数fs(i)に比例する周波数における振動の強度が大きくなる。振動の発生原因が不明であるので、振動を起こす可能性のある複数の周波数における振動の強度の平均W1(i)や振動の強度の2乗平均W2(i)を求めて、平均W1(i)や2乗平均W2(i)が最大となる基準周波数fs(isp)を回転機械30の回転速度SPDとして求めることもできる。 As described above, when the reference frequency fs (i) coincides with the rotation speed SPD, the vibration intensity at a frequency proportional to the reference frequency fs (i) increases. Since the cause of the vibration is unknown, the average W1 (i) of the vibration intensity and the squared average W2 (i) of the vibration intensity at a plurality of frequencies that may cause the vibration are obtained, and the average W1 (i) is obtained. ) And the reference frequency fs (isp) at which the squared average W2 (i) is maximized can be obtained as the rotation speed SPD of the rotating machine 30.

図8は、図4のステップS5で異常検出装置10の回転速度算出部12が行う、回転機械30の回転速度の算出の手順の別の例を示すフローチャートである。図8に示すフローチャートは、図5に示したフローチャートと、ステップS517からステップS519までの処理が異なる。図8に示すフローチャートでは、ステップS506で求めた複数の周波数における振動の強度の平均W1(i)が最大となる基準周波数fs(isp)を、回転速度SPDとして算出する。 FIG. 8 is a flowchart showing another example of the procedure for calculating the rotation speed of the rotary machine 30 performed by the rotation speed calculation unit 12 of the abnormality detection device 10 in step S5 of FIG. The flowchart shown in FIG. 8 is different from the flowchart shown in FIG. 5 in the processes from step S517 to step S519. In the flowchart shown in FIG. 8, the reference frequency fs (isp) at which the average W1 (i) of the vibration intensities at the plurality of frequencies obtained in step S506 is maximized is calculated as the rotation speed SPD.

ステップS517で、異常検出装置10は、ステップS506で求めた複数の周波数における振動の強度の平均W1(i)を求める。 In step S517, the abnormality detection device 10 obtains the average W1 (i) of the vibration intensities at the plurality of frequencies obtained in step S506.

ステップS518で、異常検出装置10は、ステップS517で求めた振動の強度の平均W1(i)とWmaxの大きさを比較する。 In step S518, the abnormality detection device 10 compares the magnitudes of W1 (i) and Wmax of the vibration intensities obtained in step S517.

ステップS509で、異常検出装置10は、WmaxにW1(i)の値を代入し、ispにiを設定する。 In step S509, the abnormality detection device 10 substitutes the value of W1 (i) for Wmax and sets i for ISP.

図9は、図4のステップS5で異常検出装置10の回転速度算出部12が行う、回転機械30の回転速度の算出の手順のさらに別の例を示すフローチャートである。図9に示すフローチャートは、図5に示したフローチャートと、ステップS527からステップS529までの処理が異なる。図9に示すフローチャートでは、ステップS506で求めた複数の周波数における振動の強度の2乗平均W2(i)が最大となる基準周波数fs(isp)を、回転速度SPDとして算出する。 FIG. 9 is a flowchart showing still another example of the procedure for calculating the rotation speed of the rotary machine 30 performed by the rotation speed calculation unit 12 of the abnormality detection device 10 in step S5 of FIG. The flowchart shown in FIG. 9 is different from the flowchart shown in FIG. 5 in the processes from step S527 to step S529. In the flowchart shown in FIG. 9, the reference frequency fs (isp) at which the squared average W2 (i) of the vibration intensities at the plurality of frequencies obtained in step S506 is maximized is calculated as the rotation speed SPD.

ステップS527で、異常検出装置10は、ステップS506で求めた複数の周波数における振動の強度の2乗平均W2(i)を求める。 In step S527, the abnormality detection device 10 obtains the squared average W2 (i) of the vibration intensities at the plurality of frequencies obtained in step S506.

ステップS528で、異常検出装置10は、ステップS527で求めた振動の強度の2乗平均W2(i)とWmaxの大きさを比較する。 In step S528, the abnormality detection device 10 compares the squared average W2 (i) of the vibration intensity obtained in step S527 with the magnitude of Wmax.

ステップS529で、異常検出装置10は、WmaxにW2(i)の値を代入し、ispにiを設定する。 In step S529, the abnormality detection device 10 substitutes the value of W2 (i) for Wmax and sets i for ISP.

本実施例による異常検出装置10は、以上に説明したように、回転速度が変化する回転機械30の回転速度を、転がり軸受20の振動データを利用して算出することができる。したがって、データロガー16から回転速度の情報を取得したり、回転機械30に回転速度センサを設置したりしなくても、異常検出装置10で回転機械30の回転速度を正確に算出することができ、コストダウンや作業効率と作業精度の向上を図ることができる。 As described above, the abnormality detection device 10 according to the present embodiment can calculate the rotation speed of the rotating machine 30 whose rotation speed changes by using the vibration data of the rolling bearing 20. Therefore, the rotation speed of the rotation machine 30 can be accurately calculated by the abnormality detecting device 10 without acquiring the rotation speed information from the data logger 16 or installing the rotation speed sensor in the rotation machine 30. , Cost reduction and improvement of work efficiency and work accuracy can be achieved.

なお、本発明は、上記の実施例に限定されるものではなく、様々な変形が可能である。例えば、上記の実施例は、本発明を分かりやすく説明するために詳細に説明したものであり、本発明は、必ずしも説明した全ての構成を備える態様に限定されるものではない。また、ある実施例の構成の一部を他の実施例の構成に置き換えることが可能である。また、ある実施例の構成に他の実施例の構成を加えることも可能である。また、各実施例の構成の一部について、削除したり、他の構成を追加・置換したりすることが可能である。 The present invention is not limited to the above embodiment, and various modifications are possible. For example, the above-described embodiment has been described in detail in order to explain the present invention in an easy-to-understand manner, and the present invention is not necessarily limited to the embodiment including all the described configurations. Further, it is possible to replace a part of the configuration of one embodiment with the configuration of another embodiment. It is also possible to add the configuration of another embodiment to the configuration of one embodiment. Further, it is possible to delete a part of the configuration of each embodiment and add / replace another configuration.

10…異常検出装置、11…データ取得部、12…回転速度算出部、13…異常特定部、14…出力部、15…振動検出器、16…データロガー、17…コンピュータ、18…バンドパスフィルタ、20…転がり軸受、21…内輪、22…外輪、23…転動体、24…保持器、25…荷重の方向、26…径方向、30…回転機械、31…軸。 10 ... Abnormality detection device, 11 ... Data acquisition unit, 12 ... Rotation speed calculation unit, 13 ... Abnormality identification unit, 14 ... Output unit, 15 ... Vibration detector, 16 ... Data logger, 17 ... Computer, 18 ... Band path filter , 20 ... rolling bearing, 21 ... inner ring, 22 ... outer ring, 23 ... rolling element, 24 ... cage, 25 ... load direction, 26 ... radial direction, 30 ... rotating machine, 31 ... shaft.

Claims (10)

回転速度が変化する回転機械を支持する転がり軸受に設置され、前記転がり軸受の振動を検出する振動検出器と、
前記振動検出器が検出した前記振動の時間領域の波形を取得するデータ取得部と、
時間領域の前記波形の包絡線を求め、高速フーリエ変換処理により時間領域の前記包絡線の波形を周波数領域の強度分布データに変換し、周波数領域の前記強度分布データを用いて前記回転機械の回転速度を求める回転速度算出部と、
前記回転速度算出部が求めた前記回転速度と、前記転がり軸受の振動の要因を示す特徴周波数とを用いて、前記転がり軸受の異常部品を特定する異常特定部と、
を備え、
前記回転速度算出部は、周波数領域の前記強度分布データについて、予め定めた範囲内の周波数である基準周波数に比例する複数の周波数における振動の強度を求め、求めた前記強度を用いて前記回転速度を求める、
ことを特徴とする、転がり軸受の異常検出装置。
A vibration detector installed on a rolling bearing that supports a rotating machine whose rotation speed changes, and a vibration detector that detects the vibration of the rolling bearing.
A data acquisition unit that acquires a waveform in the time domain of the vibration detected by the vibration detector, and
The envelope of the waveform in the time domain is obtained, the waveform of the envelope in the time domain is converted into intensity distribution data in the frequency domain by a fast Fourier transform process, and the rotation of the rotating machine is performed using the intensity distribution data in the frequency domain. The rotation speed calculation unit that obtains the speed, and
Using the rotational speed obtained by the rotational speed calculation unit and the characteristic frequency indicating the cause of vibration of the rolling bearing, an abnormality specifying unit for identifying an abnormal component of the rolling bearing, and an abnormality specifying unit.
Equipped with
The rotation speed calculation unit obtains the vibration intensity at a plurality of frequencies proportional to the reference frequency, which is a frequency within a predetermined range, with respect to the intensity distribution data in the frequency domain, and uses the obtained intensity to obtain the rotation speed. Ask,
Anomaly detection device for rolling bearings.
前記回転速度算出部は、前記基準周波数を変化させ、変化させた前記基準周波数に比例する複数の周波数における振動の前記強度を求め、求めた前記強度を用いて前記回転速度を求める、
請求項1に記載の、転がり軸受の異常検出装置。
The rotation speed calculation unit changes the reference frequency, obtains the intensity of vibration at a plurality of frequencies proportional to the changed reference frequency, and obtains the rotation speed using the obtained intensity.
The abnormality detection device for rolling bearings according to claim 1.
前記回転速度算出部は、前記基準周波数に比例する複数の周波数における振動の強度の和を求め、前記和が最大となる前記基準周波数を前記回転速度として求める、
請求項2に記載の、転がり軸受の異常検出装置。
The rotation speed calculation unit obtains the sum of the vibration intensities at a plurality of frequencies proportional to the reference frequency, and obtains the reference frequency at which the sum is maximum as the rotation speed.
The abnormality detection device for rolling bearings according to claim 2.
前記回転速度算出部は、前記基準周波数に比例する複数の周波数における振動の強度の平均を求め、前記平均が最大となる前記基準周波数を前記回転速度として求める、
請求項2に記載の、転がり軸受の異常検出装置。
The rotation speed calculation unit obtains an average of vibration intensities at a plurality of frequencies proportional to the reference frequency, and obtains the reference frequency at which the average is maximum as the rotation speed.
The abnormality detection device for rolling bearings according to claim 2.
前記回転速度算出部は、前記基準周波数に比例する複数の周波数における振動の強度の2乗平均を求め、前記2乗平均が最大となる前記基準周波数を前記回転速度として求める、
請求項2に記載の、転がり軸受の異常検出装置。
The rotation speed calculation unit obtains the squared average of the vibration intensities at a plurality of frequencies proportional to the reference frequency, and obtains the reference frequency at which the squared average is maximum as the rotation speed.
The abnormality detection device for rolling bearings according to claim 2.
異常検出装置が、回転速度が変化する回転機械を支持する転がり軸受の振動を検出する振動検出器が検出した前記振動の時間領域の波形を取得するデータ取得工程と、
異常検出装置が、時間領域の前記波形の包絡線を求め、高速フーリエ変換処理により時間領域の前記包絡線の波形を周波数領域の強度分布データに変換し、周波数領域の前記強度分布データを用いて前記回転機械の回転速度を求める回転速度算出工程と、
異常検出装置が、前記回転速度算出工程で求めた前記回転速度と、前記転がり軸受の振動の要因を示す特徴周波数とを用いて、前記転がり軸受の異常部品を特定する異常特定工程と、
を有し、
前記回転速度算出工程では、異常検出装置が、周波数領域の前記強度分布データについて、予め定めた範囲内の周波数である基準周波数に比例する複数の周波数における振動の強度を求め、求めた前記強度を用いて前記回転速度を求める、
ことを特徴とする、転がり軸受の異常検出方法。
The data acquisition process of acquiring the waveform in the time domain of the vibration detected by the vibration detector that detects the vibration of the rolling bearing that supports the rotating machine whose rotation speed changes.
The anomaly detection device obtains the envelope of the waveform in the time domain, converts the waveform of the envelope in the time domain into the intensity distribution data in the frequency domain by the fast Fourier transform process, and uses the intensity distribution data in the frequency domain. The rotation speed calculation step for obtaining the rotation speed of the rotary machine, and
An abnormality specifying step of identifying an abnormal component of the rolling bearing by the abnormality detecting device using the rotational speed obtained in the rotational speed calculation step and a characteristic frequency indicating a factor of vibration of the rolling bearing.
Have,
In the rotation speed calculation step, the abnormality detecting device obtains the vibration intensity at a plurality of frequencies proportional to the reference frequency, which is a frequency within a predetermined range, with respect to the intensity distribution data in the frequency domain, and obtains the obtained intensity. The rotation speed is obtained by using.
A method for detecting anomalies in rolling bearings, which is characterized by the fact that.
前記回転速度算出工程では、異常検出装置が、前記基準周波数を変化させ、変化させた前記基準周波数に比例する複数の周波数における振動の前記強度を求め、求めた前記強度を用いて前記回転速度を求める、
請求項6に記載の、転がり軸受の異常検出方法。
In the rotation speed calculation step, the abnormality detecting device changes the reference frequency, obtains the intensity of vibration at a plurality of frequencies proportional to the changed reference frequency, and uses the obtained intensity to determine the rotation speed. Ask,
The method for detecting an abnormality in a rolling bearing according to claim 6.
前記回転速度算出工程では、異常検出装置が、前記基準周波数に比例する複数の周波数における振動の強度の和を求め、前記和が最大となる前記基準周波数を前記回転速度として求める、
請求項7に記載の、転がり軸受の異常検出方法。
In the rotation speed calculation step, the abnormality detecting device obtains the sum of the vibration intensities at a plurality of frequencies proportional to the reference frequency, and obtains the reference frequency at which the sum is maximum as the rotation speed.
The method for detecting an abnormality in a rolling bearing according to claim 7.
前記回転速度算出工程では、異常検出装置が、前記基準周波数に比例する複数の周波数における振動の強度の平均を求め、前記平均が最大となる前記基準周波数を前記回転速度として求める、
請求項7に記載の、転がり軸受の異常検出方法。
In the rotation speed calculation step, the abnormality detecting device obtains the average of the vibration intensities at a plurality of frequencies proportional to the reference frequency, and obtains the reference frequency at which the average is maximum as the rotation speed.
The method for detecting an abnormality in a rolling bearing according to claim 7.
前記回転速度算出工程では、異常検出装置が、前記基準周波数に比例する複数の周波数における振動の強度の2乗平均を求め、前記2乗平均が最大となる前記基準周波数を前記回転速度として求める、
請求項7に記載の、転がり軸受の異常検出方法。
In the rotation speed calculation step, the abnormality detecting device obtains the squared average of the vibration intensities at a plurality of frequencies proportional to the reference frequency, and obtains the reference frequency at which the squared average is maximum as the rotation speed.
The method for detecting an abnormality in a rolling bearing according to claim 7.
JP2020127243A 2020-07-28 2020-07-28 Abnormality detection device and abnormality detection method for rolling bearings Active JP7394031B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2020127243A JP7394031B2 (en) 2020-07-28 2020-07-28 Abnormality detection device and abnormality detection method for rolling bearings

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2020127243A JP7394031B2 (en) 2020-07-28 2020-07-28 Abnormality detection device and abnormality detection method for rolling bearings

Publications (2)

Publication Number Publication Date
JP2022024573A true JP2022024573A (en) 2022-02-09
JP7394031B2 JP7394031B2 (en) 2023-12-07

Family

ID=80265851

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020127243A Active JP7394031B2 (en) 2020-07-28 2020-07-28 Abnormality detection device and abnormality detection method for rolling bearings

Country Status (1)

Country Link
JP (1) JP7394031B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116625685A (en) * 2023-04-20 2023-08-22 广东海洋大学 Rolling bearing combination fault diagnosis method based on multi-domain feature construction
JP7389303B1 (en) * 2023-10-12 2023-11-29 旭精機工業株式会社 Data extraction device and abnormality monitoring device

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001208655A (en) 2000-01-28 2001-08-03 Rion Co Ltd Failure diagnostic method and its apparatus
JP3884895B2 (en) 2000-03-10 2007-02-21 株式会社東芝 Bearing diagnosis device for rotating machine
JP5375701B2 (en) 2010-03-23 2013-12-25 新日鐵住金株式会社 Method, apparatus and program for estimating rotational speed of rotating machine
EP3239667B1 (en) 2012-01-30 2018-11-28 S.P.M. Instrument AB An analysis system
JP6676982B2 (en) 2016-01-27 2020-04-08 株式会社ジェイテクト Rolling bearing device and method of detecting abnormality of rolling bearing

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116625685A (en) * 2023-04-20 2023-08-22 广东海洋大学 Rolling bearing combination fault diagnosis method based on multi-domain feature construction
CN116625685B (en) * 2023-04-20 2023-12-15 广东海洋大学 Rolling bearing combination fault diagnosis method based on multi-domain feature construction
JP7389303B1 (en) * 2023-10-12 2023-11-29 旭精機工業株式会社 Data extraction device and abnormality monitoring device

Also Published As

Publication number Publication date
JP7394031B2 (en) 2023-12-07

Similar Documents

Publication Publication Date Title
JP5565120B2 (en) High-frequency electromagnetic vibration component removal method and high-frequency electromagnetic vibration component removal device, rolling bearing diagnosis method and bearing diagnosis device for a rotating machine
Wang et al. Rotating speed isolation and its application to rolling element bearing fault diagnosis under large speed variation conditions
Singh et al. Motor current signature analysis for bearing fault detection in mechanical systems
JP5733913B2 (en) Abnormal diagnosis method for rotating machinery
KR101482509B1 (en) Diagnosis System and Method of Bearing Defect
RU2551388C2 (en) Monitoring of axial vibration for detection of shaft displacement in turbomachine units
JP6017649B2 (en) Abnormal diagnosis method for rotating machinery
JP5293300B2 (en) Vibration monitoring device and vibration monitoring method for rotating machine
JP2007108189A (en) Method and device for diagnosing abnormality in rotor
CZ306833B6 (en) A method of detecting and locating a partial rotor-stator contact during turbine operation
JP7394031B2 (en) Abnormality detection device and abnormality detection method for rolling bearings
US20200348207A1 (en) Method for Estimating Bearing Fault Severity for Induction Motors
WO2020189014A1 (en) Abnormality diagnosis device and abnormality diagnosis method
JP5958932B2 (en) Dynamic equipment condition monitoring system, method and program
KR102068077B1 (en) A Diagnosis Apparatus For Rotating Machinery Using Complex Signals
JP5828948B2 (en) Abnormal diagnosis method for rotating machinery
Kabul et al. An alternative approach for the detection of broken rotor bars and bearing faults of induction motor based on vibration signals
Salazar-Villanueva et al. Spectral analysis for identifying faults in induction motors by means of sound
JP3918931B2 (en) Rotating body abnormality diagnosis method and apparatus
WO2022152336A1 (en) A method for monitoring turbine blade vibration
JP6283591B2 (en) Automatic vibration diagnostic equipment for rotating machinery
Thanagasundram et al. Autoregressive based diagnostics scheme for detection of bearing faults
JP7040920B2 (en) Bearing condition monitoring device and abnormality diagnosis method
Roque et al. An approach to fault diagnosis of rolling bearings
Ertunç et al. Vibration analysis based localized bearing fault diagnosis under different load conditions

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230203

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20231018

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20231024

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20231107

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20231121

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20231127

R150 Certificate of patent or registration of utility model

Ref document number: 7394031

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150