JP3673689B2 - 歯車動的性能の評価システム及び評価方法 - Google Patents

Description

【０００１】
【発明の属する技術分野】
本発明は、歯車動的性能の評価システム及び評価方法に係り、特に、離散ウェーブレット変換を用いた歯車動的性能の評価システム及び評価方法に関する。
【０００２】
【従来の技術】
道路を行き交う多くの自動車や自動二輪車、オフィス革命をもたらしたＯＡ機器、産業用ロボットなど、あらゆる分野に歯車が利用されているが、歯車は装置の中に組み込まれており、その姿を直接見ることはほとんどない。歯車は、２軸の間に動力を伝えたいとき、また速度を変えたいとき、大きなトルクを得たいときなどに用いられている重要な機械要素の一つである。歯車工学には、歯形形状および歯形の機構学的問題に関する分野、動力伝達用として必要な強さを有するように設計する設計分野、さらに運転や振動・騒音などに関する分野などがある。歯車における最近の課題は、歯車装置の大きさに対して従来よりも小型で負荷能力の大きいものを安価に作ることと、振動・騒音を低減することである。
【０００３】
ここで歯車の動的性能の評価が重要となる。一般に、歯車の動的性能とは、歯車装置の振動、騒音、歯元ひずみ等に関するものである。近年、歯車が高速・高負荷条件下で運転されるようになり、それに伴って歯車装置の振動・騒音への関心が高まっている。従来から歯車の材料としては鋼材が広く用いられているが、近年ではコスト的な面から焼結材の利用も広まっている。焼結歯車はその製造工程上、内部に気孔を有するため強度的には鋼製歯車に劣るが、制振性については優れているといわれている。
【０００４】
【発明が解決しようとする課題】
歯車の損傷は歯の折損だけではなく、歯面の接触の繰り返しによる疲労が原因で歯面がはく離するピッチングやスポーリング損傷がある。また、歯車材料自体の疲労ではないが長期の運転により潤滑油が劣化することが原因で十分な油膜が歯面間に形成されず、歯面間が金属接触することで焼付き（Ｓｃｕｆｆｉｎｇ、Ｓｃｏｒｉｎｇ）が発生する場合がある。歯車の疲労に伴う故障診断を行う際に、歯の折損を対象とする場合には歯元曲げき裂が生じると歯の剛性の値が大きく変化するので、歯車運転性能も大きく変化する。しかし、歯面の損傷の場合には、多少の歯面はく離が生じても歯の剛性に及ぼす影響は、歯の折損に比べ非常に小さい。そのため、ピッチングやスポーリングなどの歯面はく離の損傷を早期に診断するには、劣化状態をモニターする検出器や歯車運転性能波形の解析方法を工夫しなければならない。
【０００５】
従来より、歯車や軸受などの故障診断を時間・周波数解析を用いて行っている研究は国内外に存在する。しかしながら、歯車に関しては、歯の折損に関する故障診断を行っている研究はあるものの、歯面の損傷に関する研究はあまりなされていない。また、従来、歯車動的性能の評価方法としては、高速フーリエ変換などを用いての周波数解析によっていたが、動的性能の時間的周波数変化を捉えることができないため、歯車装置の故障の初期発見や歯車制振特性などを評価し難い状況にあった。
【０００６】
本発明は、以上の点に鑑み、連続ならびに離散ウェーブレット変換を用いて、動的性能の時間・周波数解析を行い、得られた時間的変化を統計的手法ならびに時系列的データ処理手法等により総合して評価することを目的とする。また、本発明では、焼結歯車と鋼製歯車の運転時の歯元ひずみ、歯車箱振動加速度、歯車箱近傍音圧等の各種パラメータについて、ウェーブレット変換を用いて時間・周波数解析し、歯車の制振特性を評価することを目的とする。
【０００７】
【課題を解決するための手段】
本発明において、具体的には、例えば、連続ならびに離散ウェーブレット変換により、歯車の歯元ひずみの時間的変化の波形を任意の周波数ごとに分解し、これからの歯車材料による歯元ひずみの周波数特性の差を明確化した。これにより焼結歯車が鋼製歯車より制振特性に優れていることがわかった。本発明は、歯車材料による制振特性の差の把握、歯車装置の故障の初期発見、故障診断、歯車の損傷・析れの把握等に有用である。
【０００８】
本発明の第１の解決手段によると、
少なくとも１組の歯車がトルクの負荷を与えられた状態でかみあわされた被測定歯車と、
前記被測定歯車を駆動する駆動部と、
前記被測定歯車の所定の測定値を測定する測定部と、
前記測定部により測定された測定値を解析する解析部と
を備え、
前記解析部は、
特定の時間の測定値に基づき、離散ウェーブレット変換により、特定の時間及び特定の振動周波数を変数として計算し、ウェーブレット成分を低周波成分及び高周波成分に分解して求め、前記歯車対の動的性能を評価するようにした歯車動的性能の評価システム装置
を提供する。
【０００９】
また、本発明の第２の解決手段によると、
少なくとも１組の歯車がトルクの負荷を与えられた状態でかみあわされた被測定歯車の歯車動的性能の評価方法であって、
前記被測定歯車を駆動し、測定された前記被測定歯車の所定の測定値を入力する入力ステップと、
入力された測定値を解析する解析ステップと
を含み、
前記解析ステップは、
特定の時間の測定値に基づき、離散ウェーブレット変換により、特定の時間及び特定の振動周波数を変数として計算し、ウェーブレット成分を低周波成分及び高周波成分に分解して求め、被測定歯車の動的性能を評価するようにした歯車動的性能の評価方法
を提供する。
【００１０】
【発明の実施の形態】
（１）歯車のかみあい
図１に、歯車のかみあいの模式図を示す。
歯車は、歯がかみあって動力を伝える機械要素である。歯車のかみあい歯数の状態を表す指標としてかみあい率ε_ａがある。かみあい率が１．ｘｘｘの場合には、歯ｂに着目すると、かみあい始めには歯ａとｂの２対の歯がかみあう。かみあいが進むと歯ｂのみでかみあい、かみあい終わりには歯ｂとｃの２対の歯でかみあう。かみあい率が２．ｘｘｘの場合には３対と２対の歯のかみあいが繰り返されることとなる。このような歯のかみあいを繰り返すことにより歯車は動力を伝えることができ、歯数ｚの歯車が一定の回転速度ｎ（ｒｐｍ）で回転する場合には、次式で示すかみあい周波数ｆ_ｚあるいはかみあい周期Ｔ_ｚで歯がかみあう。
ｆ_ｚ＝１／Ｔ_ｚ＝ｚｎ／６０
【００１１】
図２に、一般化した歯車対の振動モデルの説明図を示す。歯車の軸のねじり剛性が歯の剛性に比べて充分に大きいときには、軸のねじり剛性を無視することができ、トルクＴ_Ｑｉを伝える一対の歯車の回転方向の運動方程式は、歯をばねとして歯車本体の慣性質量を作用線方向に換算した等価慣性質量とする１自由度の非線形運動方程式として表すことができる。
その運動方程式を式に示す。
Ｍｘ”＋Ｄｘ’＋ΣＫ（ｔ，ｘ）ｘ＝Ｗ＋ΣＫ（ｔ，ｘ）ｅ（ｔ，ｘ）
【００１２】
ここで、ｘは作用線上の相対変位、Ｍは等価慣性質量、Ｄは減衰係数、ΣＫ（ｔ，ｘ）はかみあっている歯の対の剛性、Ｗは歯面法線に働く静荷重、ｅ（ｔ，ｘ）は歯形誤差である。この式の強制外力はΣＫ（ｔ，ｘ）ｅ（ｔ，ｘ）である。歯の対の剛性と強制外力はかみあい周期で周期的に変化することとなる。
【００１３】
図３に、歯車の運転性能測定についての説明図を示す。
実際に歯車の運転性能を測定する場合、図に示すように歯車箱４の振動や騒音、歯４２の歯元にひずみゲージ４３を貼り付けて歯元ひずみ（歯元応力）が測定されている。歯元ひずみは、歯車軸を含む歯車の軸系のねじり振動の挙動を間接的に表している。また、歯車箱の振動は歯車のねじり振動が軸、軸受を伝ぱして起こるものであり、その歯車箱の振動が空気を振動させることにより騒音が発生する。
【００１４】
以上のことから歯のかみあいによって起こる歯車の運転性能波形は、かみあい周波数あるいはかみあい周期に大きく依存していることがわかる。
【００１５】
（２）評価システム
図４に、本発明に係る歯車動的性能の評価システムの構成図を示す。
試験に用いた評価システムは、動力循環式歯車試験機で歯数比と軸間距離を等しくする２組の歯車対（即ち、試験歯車対１、動力伝達用歯車対２）のそれぞれをねじり軸で結んでいる。一方の歯車側のねじり軸４にはトルク負荷用カップリング５が取り付けられており、重錘式レバーを用いてトルクを負荷した状態でカップリングを締結することにより歯車に負荷を加えることができる。コップ式無段変速機６を介して三相誘導電動機７の動力が動力伝達用歯車対２のやまば歯車を駆動し、入力動力は歯車列内を循環する。
【００１６】
この例では、大歯車を被動側、小歯車を駆動側として、単位歯幅当たり歯面法線荷重Ｐ_ｎ／ｂ＝１４２Ｎ／ｍｍ（小歯車の負荷トルク３０Ｎｍ）を歯車対に負荷し、大歯車回転速度ｎ_２を１６００ｒｐｍから１００００ｒｐｍの範囲で２００ｒｐｍ毎に変化させて、歯車の動的性能を測定した。小歯車の歯車軸に取り付けたトルク負荷用カップリング５を用いて、重錘式レバーによりトルクを負荷することで歯車対に負荷を加えることができる。潤滑油はＥＰ３０９０ギヤ油で、油温を３１３±５Ｋに制御した。
【００１７】
歯車箱振動加速度については、試験歯車対１の歯車箱側面中央に取り付けた圧電型振動加速度ピックアップ８（例えば、応答周波数１Ｈｚ〜２５ｋＨｚ）によりチャージアンプ１６を介して測定した。歯車近傍音圧については歯車箱側面より法線方向に３００ｍｍ離れた位置にコンデンサ型マイクロフォン９（例えば、応答周波数５Ｈｚ〜１２．５ｋＨｚ）を固定し、騒音計１４及びペンレコーダ１５により測定した。大歯車の歯の圧縮側歯元の最弱断面位置付近にひずみゲージ（例えば、ゲージ抵抗１２０Ω、ゲージ長０．３ｍｍ）を貼付し（図３参照）、一組のブリッジを組み、スリップリング１０、ストレインアンプ１１（例えば、応答周波数ＤＣ〜１００ｋＨｚ）を介して歯元ひずみを測定した。両歯車軸に取り付けたスリット板とフォトセンサ１２、１３により、歯車の歯の位置及び回転数を測定した。これらの測定結果は、このフォトセンサ１２、１３の同期信号とともに、データレコーダ１７（例えば、４ｃｈ）に記録され、信号解析には特定の歯がかみ合うときのデータを、一例として、８回加算平均したものを利用した。さらに、データレコーダ１７の出力をコンピュータ１８で解析した。
【００１８】
ここで、試験歯車対１について説明する。図５は、試験歯車対の諸元を示す説明図である。
試験歯車対は、小歯車（ピニオン）と大歯車（ギヤ）を含み、試験歯車は大歯車である。ここでは、一例として、試験歯車には焼結歯車と鋼製歯車を用いた。相手側小歯車には鋼製浸炭硬化歯車を用いた。ヤング率およびポアソン比は焼結歯車で１５２ＧＰａ、０．２５、鋼製歯車で２０６ＧＰａ、０．３０である。焼結歯車には、１２７〜１７５μｍの粉末粒子径を有する合金鋼粉を用い、その鋼粉に黒鉛とステアリン酸亜鉛を混合した後、圧粉成形した。そして、窒素ガス雰囲気中で焼結し、ホブによる歯切りと歯面研削の後、イオン窒化を施した。鋼製歯車の素材にはＳＣＭ４４０鋼を用い、ホブによる歯切りの後、高周波焼入れし、歯面研削を行った。
【００１９】
図では、小歯車（ピニオン）及び大歯車（ギヤ）に関して、それぞれ、以下のパラメータが示される。すなわち、モジュール、基準圧力角、歯数、転位係数、歯先円直径、中心間距離、歯幅、かみあい率、精度、歯面仕上げ、研削である。
【００２０】
（３）ウェーブレット変換による時間・周波数解析
（３−１）連続ウェーブレット変換
ウェーブレット変換（ＷＴ：Ｗａｖｅｌｅｔ Ｔｒａｎｓｆｏｒｍ）は、時間的にも周波数的にも局在したアナライジングウェーブレットと呼ばれる関数ψ（ｔ）の相似変形と平行移動を利用した時間・周波数２次元解析であり、入力信号ｆ（ｔ）のＷＴは（１）式で定義される。このｆ（ｔ）が時間的変化する（時刻ｔ）各種センサによる測定値を表す。ここで、ψ^＊（ｔ）はψ（ｔ）の複素共役で、ａ、ｂはそれぞれ周波数と時間に関するパラメータである。パラメータａは、例えば、周波数０から、測定に使用される各種センサの最高応答周波数まで適宜変化させることができる。パラメータｂは、測定する時間の長さ分変化させることができる。また、ウェーブレット逆変換は（２）式で表される。ただし、ψ（ｔ）は、（３）式のアドミッシブル条件を満たされなければならない。ここで、ωは角周波数、ψ（ω）の＾（ハット）はフーリエ変換である。
【００２１】
【数２】

【００２２】
（１）式の計算には台形公式を利用し、その結果はエネルギー強度について次式で示す方法で表した。ここで、ＲｅとＩｍは、それぞれ実数部と虚数部であるウェーブレットマップを表示する際には、マップ中のウェーブレットの最大強度を１００として表示した。
【００２３】
【数３】

（３−２）離散ウェーブレット変換
【００２４】
（１）式の値を平面上の点（ｂ、１／ａ）にプロットしたウェーブレットマップは、信号の性質を知るのに便利であるが、演算を行うには（１）式の右辺の積分がそれほど簡単ではないこと、さらに（Ｗ_ψｆ）（ｂ，ａ）には多くの情報が重複していることなどのため、必ずしも効率的ではない。このことから、（１）式で示した連続ウェーブレット変換の座標を２つの整数ｊとｋによって（２^−ｊｋ、２^−ｊ）として離散化することにより、ある関数ｆ（ｔ）に関する離散ウェーブレット変換は（１）、（２）式より（Ｗ_ψｆ）（ｂ，ａ）をｄ_ｋ ^（ｊ）とおくと（６）式で与えられる。ここで、ｊはレベルと呼ばれ、周波数に関するパラメータである。また、ｋは時間に関するパラメータである。これらパラメータはａ、ｂと同様に変化させる。上述で示したような（ｂ、ａ）の値を任意にとり、変換を行う方法を連続ウェーブレット変換、後半で示したように（ｂ、ａ）を（２^−ｊｋ、２^−ｊ）に置き換えて、離散化させて変換を行う方法を離散ウェーブレット変換という。
【００２５】
また、離散ウェーブレット逆変換は（７）式で定義される。ここで、ｇ_ｊ（ｔ）は、（８）式のように定義した。また、ｆ（ｔ）の任意のレベルでの関数をｆ_ｊ（ｔ）とすると、（９）式と表すことができる。
【００２６】
【数４】

【００２７】
ここで、｛ｃ_ｋ ^（ｊ）｝は任意の数列であり、φ（ｔ）はスケーリング関数と呼ばれ、マザーウェーブレット関数ψ（ｔ）とともに（１０）、（１１）式のトゥースケール関係を満たす。ここで、｛ｐ_ｋ｝と｛ｑ_ｋ｝はトゥースケール数列である。ウェーブレット変換を用いて関数ｆ（ｔ）を分解するアルゴリズムと再構成するアルゴリズムは、（１２）式と（１３）式でそれぞれ与えられる。ここで、｛ａ_ｋ｝と｛ｂ_ｋ｝は分解アルゴリズムで使用される数列である。あるレベルでの関数ｇ_ｊ（ｔ）とｆ_ｊ（ｔ）は（８）、（９）式から計算することができ、それら関数の間には（１４）式の関係が成り立つ。
【００２８】
【数５】

【００２９】
与えられた離散データを分解するために、レベル０における数列｛ｃ_ｋ ^（０）｝は（１６）、（１７）式から求めた。数列｛β_ｋ｝は｜ｋ｜＞７のとき、０として近似した。
【００３０】
【数６】

【００３１】
本発明では、アナライジングウェーブレットとして、一例として、連続変換にはＧａｂｏｒ関数を、離散変換にはＨａａｒ関数を用いた。以下に、Ｇａｂｏｒ関数とＨａａｒ関数について示す。
Ｇａｂｏｒ関数は次式で与えられる。
【００３２】
【数７】

【００３３】
ここでω_ｐは中心周波数、ｒは周波数領域での局在の幅を決めるパラメータである。図６には、Ｇａｂｏｒ関数の時間波形図と周波数特性図を示す。Ｇａｂｏｒ関数は時間的にも周波数的にも局在性の良い関数であることがわかる。
つぎに、Ｈａａｒ関数は次のように定義される。
【００３４】
【数８】

【００３５】
図７に、Ｈａａｒ関数の時間波形図を示す。ここではψ（２^ｊｔ−ｋ）の一例を示す。この関数ではレベルｊでの矩形波の波長が１／２^ｊでｋ（ｋ＝０、１、…２^ｊ−１）によって位置が決まる。周波数パラメータａは２^−ｊ、時間パラメータｂは２^−ｊｋと置くため、ａとｂとの間にはｂ＝ａｋという関係が成り立つ。このため時間ｂと周波数１／ａとは反比例の関係になり、同時に細かくとることはできない。つまり、周波数を細かくとると時間の分解能が悪くなり、逆に、時間を細かくとると周波数の分解能が悪くなる。
【００３６】
この他に、スケーリング関数に４階のカーディナルＢスプライン関数Ｎ_４（ｔ）を用いることができる。ｍ階のカーディナルＢスプライン関数Ｎ_ｍ（ｔ）は次式で与えられる。
【００３７】
【数９】

【００３８】
図８は、４階（＝４）のスプライン関数Ｎ_４（ｔ）を用いたφ（ｔ）とψ（ｔ）を示す図である。スプライン関数Ｎ_４（ｔ）が、スケーリング関数φ（ｔ）自身である。Ｎ_４（ｔ）を用いたマザーウェーブレット関数ψ（ｔ）は、滑らかで対称性を有していることがわかる。なお、与えらた関数をＮ_４（ｔ）に基づいて分解・再構成する際に必要な（１２）、（１３）式 中の数列｛ｐ_ｋ｝、｛ｑ_ｋ｝、｛ａ_ｋ｝ 、｛ｂ_ｋ｝の値は、一例として、文献、「榊原進, ウェーヴレットビキナーズガイド, (1998), 東京電機大学出版局」、及び、「C. K. Chui著, 桜井・新井訳, ウェーブレット入門, (1994), 東京電機大学出版局」を参照した。
【００３９】
つぎに、図９に、離散ウェーブレット変換による評価方法のフローチャートを示す。図示のように、処理が開始されるとまず測定値ｆ（ｔ）を入力する（Ｓ１０）。測定値ｆ（ｔ）は、上述のように、振動加速度ピックアップ８、マイクロフォン９、ひずみゲージ４３、フォトセンサ１２、１３等の各種センサからの信号の時間的変化量である。つぎに、上述した式に従いｄ_ｋ ^（ｊ）とｃ_ｋ ^（ｊ）を計算する（Ｓ１３）。つぎにｆ_ｊ（ｔ）、ｇ_ｊ（ｔ）を計算する（Ｓ１５）。つぎに、計算されたｆ_ｊ（ｔ）、ｇ_ｊ（ｔ）を表示して（Ｓ１７）、処理を完了する。表示方法としては、表やグラフ等の適宜の手段をとることができる。
【００４０】
（４）ウェーブレット変換による解析例
図１０は、ＷＴによる解析例を示す図である。解析した信号は以下のように表され、その信号波形は図（ａ）に示すとおりである。
【００４１】
【数１０】

【００４２】
ｆ（ｔ）は振幅が時間的に変化する周期１２８μｓ（周波数７．８１ｋＨｚ）の正弦波形であり、ｔ＝１０２４μｓからはその正弦波形に周期３２μｓ（周波数３１．２５ｋＨｚ）の余弦波が加わったものである。図（ｂ）はｆ（ｔ）の信号に対して、Ｇａｂｏｒ関数を用いたウェーブレットマップを示す。図のスケールがエネルギー強度を示し、濃い部分が強いことを示す。この図より、ｆ（ｔ）には時間的に強度が変化する周波数７．８１ｋＨｚの成分に加え、ｔ＝１０２４μｓからは周波数３１．２５ｋＨｚの成分の強度が現われていることがわかる。
【００４３】
次に、上述の信号に対してＨａａｒ関数を用いて離散ウェーブレット変換を行う。ここでサンプリング時間をｔ_ｓ、サンプリング点数をＮ_ｓとすると、レベルｊでの解析周波数ｆ_ｊはｆ_ｊ＝２^ｊ／（Ｎ_Ｓ ｔ_Ｓ）で近似できる。図（ｃ）（ｄ）には、ｆ（ｔ）を周波数１５．６３ｋＨｚを境に高周波領域と低周波領域に分割・再構成した波形をそれぞれ示す。周波数１５．６３ｋＨｚ以下では、時間的に振幅の変化する周波数７．８１ｋＨｚの波形が近似的に再構成されており、また高周波領域にはｔ＝１０２４μｓから周波数３１．２５ｋＨｚの波形が再構成されている。以上の結果から、ＷＴにより信号の時間に対する周波数成分の変化を捉えることができ、任意の周波数で波形を分解・再構成できることが確認できた。
【００４４】
（５） 制振特性の評価
つぎに、制振特性の評価について述べる。
【００４５】
（５−１） 従来の評価法
まず、従来の歯車動荷重評価について説明する。
図１１は、大歯車回転速度ｎ_２に伴う動荷重比σ_ｄｍａｘ／σ_ｓｍａｘの変化を示す図である。σ_ｄｍａｘは各回転速度での動的歯元応力の最大値であり、σ_ｓｍａｘは静的歯元応力の最大値である。ここでσ_ｓｍａｘはｎ_２＝６ｒｐｍでの歯元応力である。ｆ_ｚはかみあい周波数で、ｆ_ｅは歯車対のねじり振動の固有振動数であり、その値は約２．５６ｋＨｚであった。焼結歯車および鋼製歯車ともにｎ_２の増加に伴い動荷重比は変動しており、ｎ_２＝４２００ｒｐｍ（ｆ_ｚ＝１．１２ｋＨｚ）、６０００ｒｐｍ（ｆ_ｚ＝１．６０ｋＨｚ）、９０００ｒｐｍ（ｆ_ｚ＝２．４０ｋＨｚ）で極大となっている。ｎ_２＝９０００ｒｐｍ、４２００ｒｐｍではかみあい周波数がｆ_ｅあるいはその１／２に近い周波数であったため、歯車対が共振し動荷重比が高くなったものと考えられる。また、歯車箱の固有振動数は約５ｋＨｚであったため、ｎ_２＝９０００ｒｐｍ、６０００ｒｐｍ、４２００ｒｐｍのかみあい周波数が歯車箱の固有振動数のほぼ１／２、１／３、１／４に対応して、歯車箱全体が共振し動荷重比が高くなったものと考えられる。同じ大歯車回転速度に対して、焼結歯車の動荷重比のほうが、鋼製歯車よりも概して小さい傾向にある。
【００４６】
（５−２）連続ウェーブレット変換による評価法
つぎに、連続ウェーブレット変換による評価について述べる。
図１２は、ウェーブレットマップの図である。図（ａ）は鋼製歯車、図（ｂ）は焼結歯車についての図である。この図は、動荷重比が最大となったｎ_２＝９０００ｒｐｍでの歯元ひずみ波形のウェーブレットマップを示す。このマップでは、解析した時間・周波数領域での強度の最大値を１００として正規化してある。ｎ_２＝９０００ｒｐｍにおける１対の歯のかみあう時間は約０．５２ｍｓであり、この時間を半周期とする波形の周波数は約０．９６ｋＨｚである。そのため、低周波領域での強度が高いところは歯のかみあいの進行に伴う荷重移動に起因し、高周波領域の強度が高いところは歯車軸を含む歯のねじり振動に起因していると考えられる。両歯車ともに、ピークになる時間はピッチ点付近であるが、焼結歯車の方が非常に小さい強度を示している。
【００４７】
（５−３）離散ウェーブレット変換による評価
つぎに、離散ウェーブレット変換による評価について述べる。連続ウェーブレット変換により得られたウェーブレットマップでは、各周波数における波形の形については直感的にわかりづらい。すなわち、連続ウェーブレット変換では、エネルギーの大きさが得られるものの、波形の時間的変化が得られない。そこで離散ウェーブレット変換を用いて、歯元ひずみ波形を特定の周波数で分解して考察する。
【００４８】
図１３に、歯元ひずみの分解波形図を示す。この図は、ｎ_２＝９０００ｒｐｍの歯元ひずみ波形を離散ウェーブレット変換し、レベルｊごとに再構成した波形ｇ_ｊ（ｔ）と、そのレベルｊ未満で再構成した次式で表される波形ｆ_ｊ（ｔ）を示す。
【００４９】
ここで、ｆ_ｊ（ｔ）は、低周波数成分であり、図中上側の波形を示す。ｇ_ｊ（ｔ）は、レベルｊの周波数における振動を表す高周波数成分であり、図中下側の波形を示す。ｇ_ｊ（ｔ）は分解されたウェーブレット成分である。これら離散値は、次式の関係にあり、上述のようにコンピュータにより高速に計算することができる。
ｆ_ｊ（ｔ）＝ｆ_ｊ−１（ｔ）＋ｇ_ｊ−１（ｔ）
【００５０】
例えば、この図のｆ_３（ｔ）はほぼ荷重移動による成分のみの歯元ひずみ変動を表していると考えられるが、両歯車ともに大差はない。しかし、レベル５、６、７の高周波成分になると、鋼製歯車の方が大きな振幅を示している。
【００５１】
次に高周波領域全体を考察するために、特定の周波数を境に分解・再構成する。
図１４に、大歯車回転速度と周波数推移についての図を示す。この図は、各大歯車回転速度における鋼製歯車のウェーブレットマップにおいて、最大の強度（ピーク値）における周波数と、その次に高い強度における周波数を示す。図中の黒色の曲線はｆ_ｚまたはｆ_ｚの整数倍を表し、灰色の曲線は１／（２εＴ_ｚ）の周波数を表す。ここでＴ_ｚ（＝１／ｆ_ｚ）はかみあい周期で、εはかみあい率である。歯元ひずみ波形の荷重移動に起因する成分と歯車のねじり振動に起因する成分を分離するためには、付近で歯元ひずみを分解するのが妥当であると考えられる。
【００５２】
図１５に、ｎ_２＝９０００ｒｐｍにおける歯元ひずみのｆ_ｚ以上での再構成波形およびε_ｎ−ε_ｎ＝６の波形図を示す。この図は、ｆ_ｚ＝９０００ｒｐｍにおける歯元ひずみについて、ｆ_ｚより高い周波数領域で再構成した波形と、ｎ_２＝９０００ｒｐｍの波形からｎ_２＝６ｒｐｍで測定した静的歯元ひずみ波形を引算した波形ε_ｎ−ε_ｎ＝６を示す。ε_ｎとε_ｎ＝６では波形自体の時間のオーダは異なるが、Ｔ_ｚで各波形を正規化してε_ｎ−ε_ｎ＝６を求めた。静的な歯元ひずみ変化を差し引いたε_ｎ−ε_ｎ＝６の波形は近似的に歯車のねじり振動を表していると考えられる。各歯車において２つの波形は類似していることがわかる。つまり、離散ウェーブレット変換を用いてｆ_ｚより高い周波数成分を取り出すことで、歯元ひずみ波形に含まれるねじり振動に起因する成分を取り出すことができるといえる。両歯車について比較すると、焼結歯車ではかみあい中の振幅の変化が少ないが、鋼製歯車では振幅が大きく変化し、ピッチ点付近で大きな振幅を示している。
【００５３】
図１６は、大歯車回転速度に伴う歯元ひずみの変化を示す図である。この図は、ｆ_ｚを境に歯元ひずみの波形を高周波領域と低周波領域に分解・再構成した波形の最大値と最小値との差を回転速度毎に示す。低周波領域については、両歯車ともに大きな変動はなく、ほぼ一定の値であり、概して、焼結歯車のほうがひずみの値が大きい。これはヤング率とポアソン比が鋼に比べ小さいことに由来すると考えられる。回転数の増加に伴う高周波領域の変動については、図１１の動荷重比に見られた周期的な変動と定性的に同じである。その変動の振幅については、概して焼結歯車のほうが小さい。
【００５４】
図１７（ａ）は、大歯車回転速度に伴う歯車箱振動加速度の変化を示す図である。図１７（ｂ）は、大歯車回転速度に伴う歯車箱近傍音圧の変化の図を示す図である。これらの図は、同様に、歯車箱振動加速度と歯車箱近傍音圧についても、離散ウェーブレット変換を用いてｆ_ｚ付近を境に各波形を分解・再構成し、その最大値と最小値との差を各回転数毎に示す。ｆ_ｚ以下の低周波領域では両歯車ともに大差はないが、高周波領域では鋼製歯車のほうが大きな値を示している。歯車のかみあいによって生じるねじり振動が、軸と軸受を伝ぱして歯車箱を振動させ、その歯車箱の振動が空気中に放射され騒音となる。このため、歯車のねじり振動を抑えている焼結歯車のほうが、歯車箱振動加速度と歯車箱音圧においても、ｆ_ｚ以上の高周波領域でともに低い値を示しているといえる。
【００５５】
（６）歯車運転性能評価
つぎに、歯車運転性能評価について説明する。ＷＴを用いた歯車の運転性能評価を行うために、異なる歯車材質を用いてそれら歯車の性能を測定する。歯車材質として、鋼の原料粉（クロム・モリブデン鋼に相当）から作られた焼結材と一般的に歯車に使用されている鋼材（ＪＩＳ：ＳＣＭ４４０）を用いた。焼結材は、その内部に気孔を有しているために、鋼に比べ制振性があといわれる材料である。歯車諸元は上述の図に示したとおりで、大歯車の材質を焼結材と鋼材とし、相手小歯車には鋼材を用いた。歯車運転性能の測定方法は、上述のように、大歯車回転速度を１６００ｒｐｍから１００００ｒｐｍまで２００ｒｐｍ毎に変化させ、歯元ひずみ、歯車箱振動加速度、歯車箱近傍音圧を測定した。また、動荷重を求める際に動的な歯元応力の最大値σ_ｄｍａｘと静的な歯元応力の最大値σ_ｓｍａｘが必要となる。
そこで、ｎ_２＝６ｒｐｍで測定した歯元ひずみから静的歯元応力を求めた。
【００５６】
（６―１） 従来の性能評価法
大歯車回転速度ｎ_２の増加に伴う動荷重比σ_ｄｍａｘ／σ_ｓｍａｘの変化は上述と同様である。
図１８は、（ａ）大歯車回転速度ｎ_２の増加に伴う歯車箱振動加速度レベルＶＡＬと、（ｂ）歯車箱近傍音圧レベルＳＰＬの変化を示す図である。振動加速度レベルＶＡＬおよび音圧レベルＳＰＬは、次の式で定義した。
ＶＡＬ＝２０ｌｏｇ ａ_ｒｍｓ／ａ_０
ＳＰＬ＝２０ｌｏｇ ｐ_ｒｍｓ／ｐ_０
ここで、ａ_ｒｍｓは大歯車１回転中の振動加速度波形の実効値であり、ａ_０は１ｘ１０^−２ｍ／ｓ^２とした。また、ｐ_ｒｍｓは大歯車１回転中の音圧波形の実効値であり、ｐ_０は２ｘ１０^−５Ｐａとした。焼結歯車と鋼製歯車ともに、ＶＡＬとＳＰＬはｎ_２の増加とともに変動しながら増加する傾向にある。動荷重比が極大となるｎ_２＝４２００ｒｐｍ、６０００ｒｐｍ、９０００ｒｐｍにおいて、ＶＡＬとＳＰＬは極大値を示す傾向にある。また、鋼製歯車よりも焼結歯車のほうがＶＡＬとＳＰＬのレベルの値は概して小さい。したがって、焼結歯車のほうが歯車のねじり振動に関する制振特性が優れている。すなわち、歯車対の回転方向振動を強制振動の運動方程式で表したときの減衰係数は焼結歯車のほうが小さいために、歯車箱の振動と音も焼結歯車のほうが鋼製歯車よりも低下したといえる。
【００５７】
（５−２） 連続ウェーブレット変換による性能評価法
図１９は、ｎ_２＝１８００ｒｐｍにおける歯元ひずみ波形ε、そのＷＴマップとＦＦＴによる周波数分析結果を焼結歯車および鋼製歯車について示す図である。連続ウェーブレット変換には、ガボール関数をマザーウェーブレットとして用いた。ＷＴマップとＦＦＴの結果は、最大値を１００として正規化されている。両歯車ともに、ｎ_２＝１８００ｒｐｍにおける１対の歯のかみあう時間は、２．６ｍｓであり、その時間を半周期とする波形の周波数は１９２Ｈｚである。そのため、ＷＴマップにおいてかみあい周波数ｆ_ｚ＝４８０Ｈｚより低周波の領域で強度が高いところは歯のかみあいの進行に伴う荷重移動に起因しており、ｆ_ｚより高周波の領域で強度が高いところは歯車軸を含む歯車全体のねじり振動に起因していると考えられる。鋼製歯車に比べ、焼結歯車ではｆ_ｚより高周波の領域での強度は相対的に低い。この傾向は、ほかの回転速度においても定性的に同じであった。上記のことから、図１１に示した動荷重比に関して焼結歯車のほうが鋼製歯車に比べ低い理由は、歯車材に由来する歯車軸を含む歯車のねじり振動が小さかったことが考えられる。
【００５８】
次ぎに、各大歯車回転速度におけるＷＴマップにおいて強度が最大となる周波数とその次に大きい強度となる周波数を調べた。
図２０は、鋼製歯車における大歯車回転速度ｎ_２増加に伴う各領域でのＷＴマップ中の周波数変化を示す図である。ｆ_ｚを境にして、それよりも高い周波数と低い周波数でＷＴの強度は極大値を示していることがわかる。特に、ｆ_ｚより低い周波数ｆ_ｌｏｗは、次式で示すようにｎ_２の増加に伴って増す傾向にある。
ｆ_ｌｏｗ＝（１／２ε_ａ）（ｚ_２ｎ_２／６０）＝ｆ_ｚ／（２ε_ａ）
ここで、ε_ａはかみあい率で、本歯車対の場合１．２４６である。ｚ_２は大歯車の歯数で、１６枚である。したがって、ｆ_ｌｏｗ＝０．４ｆ_ｚとなる。図に示した大歯車回転速度に対する周波数の変化は焼結歯車においても同様であり、ｆ_ｚより低い周波数ｆ_ｌｏｗは０．４ ｆ_ｚで増加していた。このような結果から、歯元ひずみは、かみあい周波数ｆ_ｚより低い周波成分（荷重移動に依存）と高い周波数成分（歯車のねじり振動に依存）で構成されていることがわかる。
【００５９】
（５−３） 離散ウェーブレット変換による性能評価法
そこで、離散ウェーブレット変換を用いて歯元ひずみ波形を分解し、ｆ_ｚを境にして再構成することで、歯元ひずみを利用して歯車材の違いによる歯車制振特性の詳細な検討を行う。離散ウェーブレット変換に用いたスケーリング関数とマザーウェーブレット関数は、４階のＢスプライン関数に基づいている。
【００６０】
図２１は、焼結歯車と鋼製歯車について、大歯車回転速度ｎ_２＝１８００ｒｐｍにおける歯元ひずみを離散ウェーブレット変換により分解した結果を示す図である。焼結歯車と鋼製歯車の各レベルにおけるｇ_ｊを比べると、特にレベル−３とレベル−４では振幅の値が焼結歯車のほうが明らかに小さい。また、かみあい周波数ｆ_ｚ以下の周波数成分の波形にほぼ相当するｆ_−６の波形を比べると、その波形の最大値は焼結歯車のほうが大きい。焼結材のヤング率とポアソン比は１５２ＧＰａと０．２５であり、鋼のそれらは２０６ＧＰａと０．３であるので、同じ荷重が歯に加わった場合、大きいひずみが生じるのは焼結材である。したがって、荷重の移動に依存するｆ_ｚ以下の周波数成分の歯元ひずみ波形は焼結歯車のほうが大きくなることが理解できる。
【００６１】
図２２は、ｎ_２＝１８００ｒｐｍにおける焼結歯車の歯元ひずみ波形をｆ_ｚを境にして再構成した波形と、ｎ_２＝１８００ｒｐｍの歯元ひずみ波形からｎ_２＝６ｒｐｍで測定した静的歯元ひずみ波形を引算した波形ε_ｎ２−ε_６ｒｐｍを示す図である。１８００ｒｐｍと６ｒｐｍでは、歯元ひずみ波形の時間のオーダは異なるが、かみあい周期Ｔ_ｚで各ひずみ波形の時間を正規化してε_ｎ２−ε_６ｒｐｍを求めた。静的な歯元ひずみを差し引いたε_ｎ２−ε_６ｒｐｍの波形は近似的に歯車のねじり振動に起因する波形変化を表していると考えられる。ｆ_ｚ以上の波形とε_ｎ２−ε_６ｒｐｍの波形を比べると、実際に歯がかみあっている時間範囲では両者の波形の変動の様子は類似していることがわかる。したがって、離散ウェーブレット変換を用いて再構成したｆ_ｚ以上の波形は、歯車のねじり振動に起因する歯元ひずみ波形を精度よく表しているといえる。
【００６２】
離散ウェーブレット変換を用いて解析した結果から、歯車のねじり振動が支配的であるｆ_ｚ以上の周波数成分の歯元ひずみ波形は、歯車材の減衰係数に依存しており、ｆ_ｚ以下の荷重移動に伴う歯元ひずみ波形は歯車材のヤング率に依存していることがわかる。
【００６３】
（７）新しい運転性能評価
図２３は、従来の動荷重比、ｆ_ｚ以下で再構成した歯元ひずみ波形の最大値、ｆ_ｚ以上で再構成した歯元ひずみ波形の振幅とＷＴを用いて求めた動荷重比のそれぞれの変化を大歯車回転速度の増加に対して示す図である。ｆ_ｚ以上で再構成した各回転速度における歯元ひずみ波形の振幅は、回転速度の増加とともに変動しており、その変動の様子は従来の動荷重比の変動と似ている。なお、当然のことであるがｎ_２＝６ｒｐｍでの値はほぼ０に近い。ｆ_ｚ以下で再構成した各回転速度における歯元ひずみ波形の最大値は、ｎ_２＝６ｒｐｍでの値よりもわずかに小さい傾向にあるが、回転速度に関わらずほぼ一定の値であるといえる。このことは、ｆ_ｚを境にして再構成したｆ_ｚ以下の歯元ひずみ波形の最大値は、静的な歯元ひずみの最大値にほぼ相当すると考えられる。そこで、各回転速度毎に元の歯元ひずみ波形の最大値を再構成したｆ_ｚ以下の歯元ひずみの最大値で除して動荷重比を求めた。回転速度に対するＷＴによる新しい動荷重比の変動の傾向と従来の動荷重比の傾向は概して同じであることがわかる。しかし、従来の動荷重比に比べＷＴによる動荷重比の値のほうが高い。歯のかみあいにおいて、２対かみあいから１対かみあいへ移行するときの歯元ひずみ波形の変化は急激で、非常に高周波の成分に相当する。したがって、再構成したｆ_ｚ以下の歯元ひずみ波形に急激な変化は含まれないことが原因であると考えられる。しかしながら、多くの歯車装置は非常に低速で回転させることが困難である。そのため、設計段階で設定した安全率（動荷重比）が実際に製作された歯車装置に対し妥当であるかを判断する基準にここで提案したＷＴによる新しい動荷重比評価法を利用することができると考えられる。
【００６４】
図２４は、焼結歯車に関するｎ_２＝１８００ｒｐｍにおける歯車箱振動加速度波形の離散ウェーブレット変換解析結果を示す図である。大歯車１回転（３０ｍｓ）中の振動波形は低周波成分と高周波成分から構成されていることが理解できる。また、図２５は、歯元ひずみ波形と同様に、（ａ）振動加速度波形と（ｂ）歯車箱近傍音圧波形を、ｆ_ｚを境に再構成し、回転速度の増加に対する各波形の実効値の変化を示す図である。焼結歯車ならびに鋼製歯車に関する振動加速度波形と音圧波形ともにｆ_ｚ以下で再構成された波形の実行値は回転速度ｎ_２の２乗に比例して増加している。すなわち、両歯車ともに同じトルクが負荷されているので、荷重移動に依存するｆ_ｚ以下の振動と音圧波形は歯車回転エネルギーで決まるといえる。したがって、材質が異なる焼結歯車と鋼製歯車ともにほぼ同じ値を示している。歯車のねじり振動に起因するｆ_ｚ以上の振動と音圧の実効値は４２００ｒｐｍ、６０００ｒｐｍ、９０００ｒｐｍで極大を示しており、元の振動と音圧の実効値の変化の様子と同じである。このことは、各波形の実効値の変動は歯車のねじり振動に依存してことを明確に示している。歯車箱振動加速度は、歯のかみあいによって生じるねじり振動が軸と軸受を伝ぱして歯車箱を振動させることにより起こる。そして、歯車箱振動が空気中に放射され音圧（騒音）が発生する。そのために、歯車のねじり振動が抑えられている焼結歯車のほうが、歯車箱振動と音圧に関しても鋼製歯車より低い値を示していることが理解できる。
【００６５】
【発明の効果】
１）本発明の評価システムによると、歯元ひずみ波形において、歯車のねじり振動に起因する歯元ひずみは、かみあい周波数ｆ_ｚよりも高い周波数成分で構成されていたということを評価することができる。また、その歯元ひずみの高周波成分は、焼結歯車のほうが鋼製歯車よりも小さな値を示していることから、焼結歯車は制振特性において優れているといえ、歯車運転時の振動・騒音の原因となる歯車のねじり振動を低く抑えていることがわかった。
２）本発明によると、焼結歯車は鋼製歯車に比べ制振特性に優れるため、歯車装置に焼結歯車を適用することは振動・騒音を低減するという点で有効であることが評価することができる。
３）本発明によると、歯元ひずみのウェーブレット変換を行うことで、従来の歯元ひずみを用いた歯車動荷重評価では得ることができない歯車動的性能評価を行うことができた。
【００６６】
さらに、本発明の歯車運転性能評価システム及び評価方法については、歯車装置の検査に応用できると考えられる。また、本発明によると、かみあい周波数ｆ_ｚを境に歯車装置の性能波形を再構成し、歯車荷重伝達に関する波形と歯車ねじり振動に関する波形から、装置の設計段階で設定した安全率が妥当であるかを、新しい動荷重評価法を用いて判定することが可能である。また、本発明によると、材質の違いによる減衰係数の違いは、主にかみあい周波数ｆ_ｚ以上の周波数成分の波形に反映され、相対的な減衰係数の差も判別することができる。
【図面の簡単な説明】
【図１】歯車のかみあいの模式図。
【図２】一般化した歯車対の振動モデルの説明図。
【図３】歯車の運転性能測定についての説明図。
【図４】本発明に係る歯車動的性能の評価システムの構成図。
【図５】試験歯車対の諸元を示す説明図。
【図６】Ｇａｂｏｒ関数の時間波形図と周波数特性図。
【図７】Ｈａａｒ関数の時間波形図。
【図８】４階のスプライン関数Ｎ_４（ｔ）を用いたφ（ｔ）とψ（ｔ）を示す図。
【図９】離散ウェーブレット変換による評価方法のフローチャートを示す図。
【図１０】ＷＴによる解析例を示す図。
【図１１】大歯車回転速度ｎ_２に伴う動荷重比σ_ｄｍａｘ／σ_ｓｍａｘの変化を示す図。
【図１２】ウェーブレットマップを示す図。図（ａ）は鋼製歯車、図（ｂ）は焼結歯車についての図。
【図１３】歯元ひずみの分解波形図。
【図１４】大歯車回転速度と周波数推移についての図。
【図１５】ｎ_２＝９０００ｒｐｍにおける歯元ひずみのｆ_ｚ以上での再構成波形およびε_ｎ−ε_ｎ＝６の波形図。
【図１６】大歯車回転速度に伴う歯元ひずみの変化を示す図。
【図１７】（ａ）は、大歯車回転速度に伴う歯車箱振動加速度の変化を示す図、及び（ｂ）は、大歯車回転速度に伴う歯車箱近傍音圧の変化の図を示す図。
【図１８】（ａ）は、大歯車回転速度ｎ_２の増加に伴う歯車箱振動加速度レベルＶＡＬと、（ｂ）歯車箱近傍音圧レベルＳＰＬの変化を示す図。
【図１９】ｎ_２＝１８００ｒｐｍにおける歯元ひずみ波形ε、そのＷＴマップとＦＦＴによる周波数分析結果を焼結歯車および鋼製歯車について示す図。
【図２０】鋼製歯車における大歯車回転速度ｎ_２増加に伴う各領域でのＷＴマップ中の周波数変化を示す図。
【図２１】焼結歯車と鋼製歯車について、大歯車回転速度ｎ_２＝１８００ｒｐｍにおける歯元ひずみを離散ウェーブレット変換により分解した結果を示す図。
【図２２】ｎ_２＝１８００ｒｐｍにおける焼結歯車の歯元ひずみ波形をｆ_ｚを境にして再構成した波形と、ｎ_２＝１８００ｒｐｍの歯元ひずみ波形からｎ_２＝６ｒｐｍで測定した静的歯元ひずみ波形を引算した波形ε_ｎ２−ε_６ｒｐｍを示す図。
【図２３】従来の動荷重比、ｆ_ｚ以下で再構成した歯元ひずみ波形の最大値、ｆ_ｚ以上で再構成した歯元ひずみ波形の振幅とＷＴを用いて求めた動荷重比のそれぞれの変化を大歯車回転速度の増加に対して示す図。
【図２４】焼結歯車に関するｎ_２＝１８００ｒｐｍにおける歯車箱振動加速度波形の離散ウェーブレット変換解析結果を示す図。
【図２５】歯元ひずみ波形と同様に、（ａ）振動加速度波形と（ｂ）歯車箱近傍音圧波形を、ｆ_ｚを境に再構成し、回転速度の増加に対する各波形の実効値の変化を示す図。
【符号の説明】
１ 試験歯車対
２ 動力伝達用歯車対
３、４ ねじり軸
５ トルク負荷用カップリング
６ コップ式無段変速機
７ 三相誘導電動機
８ 圧電型振動加速度ピックアップ
９ マイクロフォン
１０ スリップリング
１１ ストレインアンプ
１２，１３ フォトセンサ
１４ 騒音計
１５ ペンレコーダ
１７ データレコーダ
１８ コンピュータ

Claims (7)

1. 少なくとも１組の歯車がトルクの負荷を与えられた状態でかみあわされた被測定歯車と、
   前記被測定歯車を駆動する駆動部と、
   前記被測定歯車の動的性能に関する所定の時間変化量である測定値ｆ（ｔ）を測定する測定部と、
   前記測定部により測定された測定値ｆ（ｔ）を、離散ウェーブレット変換により解析する解析部と
   を備え、
   前記解析部は、
   測定値ｆ（ｔ）に基づき、振動周波数及び時間に関する変数ｃ _ｋ ^（ｊ） 及びｄ _ｋ ^（ｊ） を次式（１６）、（１７）及び（１２）に基づき計算し、
   

   

   （ここで、ａ、ｂはそれぞれ周波数と時間に関するパラメータ、ｊはレベルと呼ばれ周波数に関するパラメータ、ｋは時間に関するパラメータである。）
   マザーウェーブレット関数ψ（ｔ）及びスケーリング関数φ（ｔ）として次式のｍ階カーディナルＢスプライン関数Ｎ _ｍ （ｔ）を用いて、次式（８）及び（９）に基づき離散ウェーブレット変換を行うことにより、ウェーブレット成分を低周波成分ｆ _ｊ （ｔ）及び高周波成分ｇ _ｊ （ｔ）に分解して計算し、
   

   

   低周波成分ｆ _ｊ （ｔ）及び高周波成分ｇ _ｊ （ｔ）を表示することにより、歯車対の動的性能を評価するようにした歯車動的性能の評価システム。
2. 前記測定部で測定する所定測定値は、騒音に関するデータ、回転数、歯の位置、振動加速度及び歯元ひずみのいずれかひとつ又は複数であることを特徴とする請求項１に記載の歯車動的性能の評価システム。
3. 前記測定部で測定する所定測定値は、歯元ひずみであり、
   前記解析部は、振動特性の評価、運転性能の評価又は故障・損傷診断を行うことを特徴とする請求項１又は２に記載の歯車動的性能の評価システム。
4. 時間ｔにおける計測値ｆ（ｔ）に関する離散的ウェーブレット変換が、整数ｊを周波数に相当するパラメータ、整数ｋが時間に相当するパラメータ、ψはマザーウェーブレット関数としたとき、（Ｗ_ψｆ）（２^−ｊｋ，２^−ｊ）＝ｄ_ｋ ^（ｊ）とおくと、次式で与えられることを特徴とする請求項１乃至３いずれかに記載の歯車動的性能の評価システム。
   

   
5. ｆ_ｊ（ｔ）＝ｆ_ｊ−１（ｔ）＋ｇ_ｊ−１（ｔ） の関係となる低周波成分ｆ_ｊ（ｔ）、高周波成分ｇ_ｊ（ｔ）を求めることを特徴とする４又は５に記載の歯車動的性能の評価システム。
6. 少なくとも１組の歯車がトルクの負荷を与えられた状態でかみあわされた被測定歯車の歯車動的性能の評価方法であって、
   前記被測定歯車を駆動し、測定された前記被測定歯車の動的性能に関する所定の時間変化量である測定値ｆ（ｔ）を入力する入力ステップと、
   入力された測定値ｆ（ｔ）を、離散ウェーブレット変換により解析する解析ステップと
   を含み、
   前記解析ステップは、
   測定値ｆ（ｔ）に基づき、振動周波数及び時間に関する変数ｃ _ｋ ^（ｊ） 及びｄ _ｋ ^（ｊ） を次式（１６）、（１７）及び（１２）に基づき計算し、
   

   

   （ここで、ａ、ｂはそれぞれ周波数と時間に関するパラメータ、ｊはレベルと呼ばれ周波数に関するパラメータ、ｋは時間に関するパラメータである。）
   マザーウェーブレット関数ψ（ｔ）及びスケーリング関数φ（ｔ）として次式のｍ階カーディナルＢスプライン関数Ｎ _ｍ （ｔ）を用いて、次式（８）及び（９）に基づき離散ウェーブレット変換を行うことにより、ウェーブレット成分を低周波成分ｆ _ｊ （ｔ）及び高周波成分ｇ _ｊ （ｔ）に分解して計算し、
   

   

   低周波成分ｆ _ｊ （ｔ）及び高周波成分ｇ _ｊ （ｔ）を表示することにより、歯車対の動的性能を評価するようにした歯車動的性能の評価方法。
7. 前記測定部で測定する所定測定値は、騒音に関するデータ、回転数、歯の位置、振動加速度及び歯元ひずみのいずれかひとつ又は複数であることを特徴とする請求項６に記載の歯車動的性能の評価方法。

Images