JP3686871B2 - 歯付ベルトの騒音予測方法及び装置、並びにプログラム - Google Patents
歯付ベルトの騒音予測方法及び装置、並びにプログラム Download PDFInfo
- Publication number
- JP3686871B2 JP3686871B2 JP2002050863A JP2002050863A JP3686871B2 JP 3686871 B2 JP3686871 B2 JP 3686871B2 JP 2002050863 A JP2002050863 A JP 2002050863A JP 2002050863 A JP2002050863 A JP 2002050863A JP 3686871 B2 JP3686871 B2 JP 3686871B2
- Authority
- JP
- Japan
- Prior art keywords
- belt
- displacement
- toothed
- calculating
- pulley
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims description 64
- 238000006073 displacement reaction Methods 0.000 claims description 109
- 230000005855 radiation Effects 0.000 claims description 64
- 238000004364 calculation method Methods 0.000 claims description 51
- 238000004458 analytical method Methods 0.000 claims description 16
- 230000004044 response Effects 0.000 claims description 16
- 230000008859 change Effects 0.000 claims description 14
- 230000005284 excitation Effects 0.000 claims description 13
- 230000000694 effects Effects 0.000 claims description 12
- 230000007246 mechanism Effects 0.000 description 11
- 238000012360 testing method Methods 0.000 description 9
- 239000011295 pitch Substances 0.000 description 8
- 238000010586 diagram Methods 0.000 description 6
- 230000009467 reduction Effects 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000007613 environmental effect Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10T—TECHNICAL SUBJECTS COVERED BY FORMER US CLASSIFICATION
- Y10T74/00—Machine element or mechanism
- Y10T74/18—Mechanical movements
- Y10T74/18568—Reciprocating or oscillating to or from alternating rotary
- Y10T74/18832—Reciprocating or oscillating to or from alternating rotary including flexible drive connector [e.g., belt, chain, strand, etc.]
- Y10T74/18848—Reciprocating or oscillating to or from alternating rotary including flexible drive connector [e.g., belt, chain, strand, etc.] with pulley
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Description
【発明の属する技術分野】
本発明は、歯付ベルトの騒音予測方法及び装置、並びにプログラムに関するものである。
【0002】
【従来の技術】
歯付ベルトの騒音は、従来から問題となっており、従来より実験的に歯付ベルトの発生騒音の特徴を調べ、実験結果をベースに騒音発生メカニズムや低騒音化を検討する研究が行われていた。なお、従来の実験結果の多くは、歯付ベルトの騒音メカニズムは、プーリ歯頂部とベルト歯溝部の衝突による噛合衝突音であるとしていた。
【0003】
【発明の解決しようとする課題】
さらに、歯付ベルトの振動、騒音に関係するパラメータとして、ベルト歯形形状、ベルト材料構成と材料物性、プーリ歯形形状、プーリ回転数及び張力(負荷)が考えられるが、従来の技術は実験的検討が主体であるため、これらのパラメータ変更に伴う騒音の予測評価をするには、ベルト、プーリの製作、試験が必要であり、大変な時間と費用を要するという問題があった。
【0004】
また、実験結果は、実験方法や装置等によって少しずつ異なる可能性があるため、再現性、実験精度の面で問題があった。更に、上述した歯付ベルトのパラメータの全てについて実験を行うのは困難であり、騒音発生メカニズムや騒音低減方法等についての体系的な整理をするのは難しいという問題があった。一方で、歯付ベルトの騒音を純解析的に予測する技術は、従来全く見られなかった。
【0005】
本発明は、上記問題点に鑑みてされたものであり、ベルト噛合部において、プーリ及びベルトの多角形効果、噛合部でのベルトの変形という要因によって生じるベルトの変位がベルトを加振し、ベルトの弦振動によって騒音が放射されるという歯付ベルトの騒音発生メカニズムに基づいて、歯付ベルトの騒音を純解析的に予測することができる歯付ベルトの騒音予測方法、装置、並びにプログラムを提供するものである。
【0006】
【課題を解決するための手段】
前記課題を解決するために本発明の請求項1に記載の歯付ベルトの騒音予測方法は、歯付ベルトの両端を歯付プーリで支持されたベルトスパンのベルト噛み合い時に、歯付プーリの多角形効果及び歯付ベルトのベルト噛合い点近傍でのベルトの変形により生じるベルト走行方向に対して直角方向のベルト変位に起因するベルトの弦振動により発射される騒音を予測する歯付ベルトの騒音予測方法であって、ベルト噛合い点近傍で生じるベルト走行方向に対して直角方向のベルト変位を計算するステップと、前記ベルト変位を、プーリ回転数とプーリ歯数等の諸元より、時間の関数に変換して、変動速度を計算するステップと、前記変動速度及び前記ベルトの寸法から、前記ベルトの弦振動による放射音圧レベルを計算するステップとを備えることを特徴とする。
【0007】
本発明の請求項5に記載の歯付ベルトの騒音予測装置は、歯付ベルトの両端を歯付プーリで支持されたベルトスパンのベルト噛み合い時に、歯付プーリの多角形効果及び歯付ベルトのベルト噛合い点近傍でのベルトの変形により生じるベルト走行方向に対して直角方向のベルト変位に起因するベルトの弦振動により発射される騒音を予測する歯付ベルトの騒音予測装置であって、ベルト噛合い点近傍で生じるベルト走行方向に対して直角方向のベルト変位を計算するベルト変位計算手段と、前記ベルト変位を、プーリ回転数とプーリ歯数等の諸元より、時間の関数に変換して、変動速度を計算する変動速度計算手段と、前記変動速度及び前記ベルトの寸法から、前記ベルトの弦振動による放射音圧レベルを計算する放射音圧レベル計算手段とを備えることを特徴とする。
【0008】
本発明の請求項9に記載の歯付ベルトの騒音予測プログラムは、コンピュータを、請求項5のような歯付ベルトの騒音予測装置として機能させるためのプログラムである。
【0009】
本発明の請求項2に記載の歯付ベルトの騒音予測方法は、歯付ベルトの両端を歯付プーリで支持されたベルトスパンのベルト噛み合い時に、歯付プーリの多角形効果及び歯付ベルトのベルト噛合い点近傍でのベルトの変形により生じるベルト走行方向に対して直角方向のベルト変位に起因するベルトの弦振動により発射される騒音を予測する歯付ベルトの騒音予測方法であって、ベルト噛合い点近傍で生じるベルト走行方向に対して直角方向のベルト変位を計算するステップと、前記ベルト変位を、プーリ回転数とプーリ歯数等の諸元より時間の関数に変換し、時間、或いは周波数の関数としての変位、或いは変動速度を計算するステップと、前記変換された変位、或いは変動速度を、前記ベルトスパンの噛合端の加振力として、ベルト弦振動の応答変位、もしくは振動モード、及びベルト弦振動の速度を計算するステップと、前記ベルト弦振動の応答変位、もしくは振動モード、及び前記ベルトの寸法から、前記ベルトの弦振動による音響放射効率を計算するステップと、前記ベルト弦振動の速度、及び前記音響放射効率から、放射音圧レベルまたは音響放射パワーレベルを計算するステップとを備えることを特徴とする。
【0010】
本発明の請求項6に記載の歯付ベルトの騒音予測装置は、歯付ベルトの両端を歯付プーリで支持されたベルトスパンのベルト噛み合い時に、歯付プーリの多角形効果及び歯付ベルトのベルト噛合い点近傍でのベルトの変形により生じるベルト走行方向に対して直角方向のベルト変位に起因するベルトの弦振動により発射される騒音を予測する歯付ベルトの騒音予測装置であって、ベルト噛合い点近傍で生じるベルト走行方向に対して直角方向のベルト変位を計算するベルト変位計算手段と、前記ベルト変位を、プーリ回転数とプーリ歯数等の諸元より時間の関数に変換し、時間、或いは周波数の関数としての変位、或いは変動速度を計算する変動速度計算手段と、前記変換された変位、或いは変動速度を、前記ベルトスパンの噛合端の加振力として、ベルト弦振動の応答変位、もしくは振動モード、及びベルト弦振動の速度を計算するベルト弦振動計算手段と、前記ベルト弦振動の応答変位、もしくは振動モード、及び前記ベルトの寸法から、前記ベルトの弦振動による音響放射効率を計算する音響放射効率計算手段と、前記ベルト弦振動の速度、及び前記音響放射効率から、放射音圧レベルまたは音響放射パワーレベルを計算する放射音圧レベルまたは音響放射パワーレベル計算手段とを備えることを特徴とする。
【0011】
本発明の請求項10に記載の歯付ベルトの騒音予測プログラムは、コンピュータを、請求項6のような歯付ベルトの騒音予測装置として機能させるためのプログラムである。
【0012】
請求項3,7,及び11に記載の歯付ベルトの騒音予測方法及び装置、並びにプログラムは、前記ベルト変位が、有限要素法解析を適用して計算されることを特徴とする。
【0013】
請求項4,8,及び12に記載の歯付ベルトの騒音予測方法及び装置、並びにプログラムは、前記有限要素法解析が、前記ベルト変位を、ベルトに固定した点の噛み合い進行に伴う噛合部までのベルト走行方向に対して直角方向の変化として計算することを特徴とする。
【0014】
上記記載の歯付ベルトの騒音予測方法及び装置、並びにプログラムによると、有限要素法解析を適用することによって、噛み合い点付近のベルト変位、即ち、ベルトに固定した点(例えば歯山中央の心線位置、或いは歯間中央の心線位置)の噛み合い進行に伴う噛合部までのベルト走行方向に対して直角方向の変化を計算する。そして、この変位をもとに、両端を歯付プーリで支持されたベルト表面の放射音圧レベル或いは音響放射パワーレベルを計算することができる。
【0015】
従って、ベルト噛合部で生じるベルト変位に起因するベルトの弦振動によって騒音が発射されるという歯付ベルトの騒音メカニズムに基づいて、歯付ベルトの発生騒音を体系的に整理することができる。そして、歯付ベルトの発生騒音を事前に予測することができることにより、歯付ベルトの低騒音化が可能となり、環境騒音の改善に貢献することができる。また、歯付ベルトの発生騒音の事前予測が可能となることにより、発生騒音の評価や低騒音化対策のための試験が不要となる、或いは簡素化されることで人的、物的な資源の節約となる。
【0016】
【発明の実施の形態】
以下、図面を参照しつつ、本発明における歯付ベルトの騒音予測方法及び装置、並びにプログラムについて説明する。
【0017】
まず、本実施形態が適用される歯付ベルト1及び歯付ベルトに噛合する歯付プーリ2,3,4の伝動機構について、図5及び図6に基づいて説明する。図5は、本実施形態に係る歯付ベルトと歯付プーリの伝動機構の概略図である。図6は、本実施形態に係る歯付ベルトと駆動歯付プーリの断面略図である。
【0018】
図5に示すように、歯付ベルト1は、回転軸方向に沿った歯部を外周部にそれぞれ有する駆動歯付プーリ2と、アイドラ歯付プーリ3と、従動歯付プーリ4の3つの歯付プーリと噛合している。ここで、駆動歯付プーリ2と、アイドラ歯付プーリ3と、従動歯付プーリ4とでベルトスパンを構成している。また、駆動歯付プーリ2にはモータMが備えられており、モータMを駆動することにより、歯付ベルト1は、駆動歯付プーリ2によって図中反時計回りに駆動させられ、アイドラ歯付プーリ3及び従動歯付プーリ4を図中反時計回りに回転させる。
【0019】
図6に示すように、歯付ベルト1の内側面には、その幅方向に沿って多数の円弧歯形の歯部1a、および歯底部1bが形成されている。また、駆動歯付プーリ2のプーリ中心を原点として、ベルト走行方向に対して直角方向、即ち原点から歯付ベルト1と駆動歯付プーリ2とが噛み合う定位置に向かってy軸とし、ベルト走行方向に対して平行方向、即ち原点から歯付ベルト1の走行方向と逆に向かってx軸として座標軸を採っている。従って、歯付ベルト1は−x方向に進行する。なお、歯付ベルト1は、円弧形のベルト歯部1aを有しているが、歯付ベルトの歯型として、台形他様々な形状が考えられる。
【0020】
[ 実施形態1]
次に、第1の実施形態における歯付ベルトの騒音予測装置10について、図2に基づいて説明する。図2は、第1の実施形態に係る歯付ベルトの騒音予測装置のブロック線図である。
【0021】
第1の実施形態においては、歯付ベルトの騒音予測装置10は、ベルト変位計算部11と、時間間隔計算部12と、変位、或いは変動速度計算部13と、ベルト弦振動計算部14と、音響放射効率計算部15と、放射音圧レベルまたは放射パワーレベル計算部16とを有している。
【0022】
まず、歯付ベルト1の噛み合い1ピッチ、或いは複数ピッチ分をnステップに分割し、n個、或いは(n+1)個の噛み合いの各ステップにおいて、噛み合い進行に伴うベルトに固定した点のy座標位置の変化、或いは噛み合い進行に伴う噛み合い点付近に固定されたx座標に対するベルトのy座標位置の変化がベルト変位計算部11によって求められる。ここで、nは正の整数である。
【0023】
ここで、噛み合い点の各ステップに対するベルト変位を、有限要素法解析により、ベルトに固定した心線位置の噛み合い進行に伴うy座標位置の変化として求める場合、プーリ中心を原点として図6に示す座標軸を採り、張り側ベルトの他端は、張り側張力を保持しながら心線方向(−x方向)に移動させ、y方向はフリーにし、緩み側ベルトの他端は、緩み側張力を保持しながら心線方向に移動させるような有限要素法解析条件を用いる。有限要素法解析結果の例として図7に示すような噛み合いステップと歯付ベルト1のy座標位置の関係が得られる。
【0024】
また、歯付ベルトの騒音予測装置10の操作者が、キーボードなどの入力装置から解析対象となる駆動歯付プーリ2のプーリ径D(m)、歯ピッチp(m)或いはプーリ歯数Zを入力し、また、駆動歯付プーリ2の回転数N(rpm)が入力され、時間間隔演算部12によって、1ステップ間の時間間隔Δt(sec)が求められる。具体的には、1ステップ間の時間間隔Δtは次式により計算される。
Δt=(60×p)/(N×π×D×n)=60/(N×Z×n)
【0025】
次に、変位、或いは変動速度計算部13によって、ベルト変位計算部11から入力されたベルト変位Δy(m)と、時間間隔演算部12から入力された1ステップ間の時間間隔Δtとにより、y座標位置の変化が、具体的な時間の関数として求められる。また、y方向の変動速度v(m/sec)が次式により時間の関数として計算される。
v=Δy/Δt
なお、変位、或いは変動速度計算部13においては、ベルト変位Δyをフーリエ変換により変位、或いは変動速度を周波数の関数として求めることもできる。
【0026】
次に、ベルト弦振動計算部14によって、変位、或いは変動速度計算部13で求めたy方向変位をベルトスパンの噛合端の境界条件(加振力)としてベルト弦振動の応答変位とそのモードを求める。演算方法としては、ベルトを図8に示すようにマスとバネ、並びにダッシュポットからなる質点系に置き換え、以下の式によりベルト各点の応答変位、モードを求める。なお、演算方法として有限要素法解析等を用いても良い。
【0027】
【数1】
【0028】
次に、音響放射効率計算部15によって、ベルト弦振動計算部14で求めた、加振角周波数ωに対するベルト各点(質点)の変位である振動モード{u}と、ベルト寸法(スパン長さl(m)、幅b(m))とから、ベルト振動による音響放射効率を計算する。弦振動の1次の固有値以下の周波数を対象とする場合、音響放射効率σは、例えばベルトを等価な面積を有する自由空間に置かれた半径a(m)の円形ピストンと仮定することにより、以下の式により求める。なお、厳密にベルト振動の音響放射効率σを求めるため、上記方法に替えて、境界要素法等を用いて解析的に求めたり、実験的に求めたりする方法を用いても良い。
【0029】
【数2】
【0030】
そして、放射音圧レベルまたは音響放射パワーレベル計算部16によって、ベルト弦振動計算部14で求めたベルト弦振動の速度と、音響放射効率計算部15で求めた音響放射効率から、音響放射パワーレベル或いは放射音圧レベルを以下の式により計算して、出力する。
【0031】
【数3】
【0032】
ここで、図2に示されている歯付ベルトの騒音予測装置10の各部11〜16は、例えば汎用のパーソナルコンピュータによって構成されている。かかるパーソナルコンピュータには、CPU、ROM、RAM、ハードディスク、FDやCDの駆動装置などのハードウェアが収納されており、ハードディスクには、パーソナルコンピュータを歯付ベルトの騒音予測プログラム(このプログラムは、CD−ROMやFD、MOなどのリムーバブルな記録媒体に記録しておくことにより、様々なコンピュータにインストールすることが可能である)を含む各種のソフトウェアが記録されている。そして、これらのハードウェアおよびソフトウェアが組み合わされることによって、上述の各部11〜16が構築されている。更に具体的には、ベルト変位計算部11は、市販の有限要素法解析プログラムが記憶されたものから構成される。
【0033】
また、歯付ベルトの騒音予測装置10の各部11〜16で得られるデータは、図示しないディスプレイに表示されたり、プリンタで印刷されたりすることで、歯付ベルトの騒音予測装置10の操作者に通知される。
【0034】
次に、図1を参照しつつ、第1の実施形態における歯付ベルトの騒音予測方法の手順について説明する。図1は、第1の実施形態における歯付ベルトの騒音予測方法の手順について説明したフローチャートである。
【0035】
ステップS1では、歯付ベルトの噛み合い1ピッチ、或いは複数ピッチ分をnステップに分割し、n個、或いは(n+1)個の噛み合いの各ステップにおいて、噛み合い進行に伴うベルトに固定した点のy座標位置の変化、或いは噛み合い進行に伴う噛み合い点付近に固定されたx座標に対するベルトのy座標位置の変化を求める。ここで、nは正の整数である。
【0036】
また、ステップS2において、歯付ベルトの騒音予測装置10の操作者が、キーボードなどの入力装置から解析対象となる駆動歯付プーリ2のプーリ径D、歯ピッチp或いはプーリ歯数Z、また、駆動歯付プーリ2の回転数Nが入力され、これらにより、上記1ステップ間の時間間隔Δtを求める。
【0037】
そして、ステップS3において、ステップS1で求めたy座標位置の変位ΔyとステップS2で求められた時間間隔Δtを用いて、各時間ステップに対するベルトのy座標の変位Δy、或いはy方向速度v=Δy/Δtを時間の関数として求める。また、これらをフーリエ変換することにより変位Δy、或いはy方向速度vは周波数の関数としても求められる。
【0038】
次に、ステップS4において、ステップS3で時間、または周波数の関数として求められたベルトの変位Δy、或いはy方向速度vをベルトスパンの噛合端の境界条件(加振力)として与え、ベルト弦振動の応答変位、或いは応答速度とそのモードを求める。
【0039】
次に、ステップS5において、ステップS4で求められた加振力の周波数、及びベルト弦振動の応答モードを用いてベルトの音響放射効率を計算する。
【0040】
そして、ステップS6において、ステップS4で求められたベルト弦振動の速度と振動モード、及びステップS5で求められた音響放射効率を用いて音響放射パワーレベル、或いは放射音圧レベルを計算して、出力する。
【0041】
ここで、前述と同様に、各ステップS1〜S6は、例えば汎用のパーソナルコンピュータを用いて行われ、特にステップS1は市販の有限要素法解析プログラムを用いることができる。また各ステップS1〜S6で得られるデータは、図示しないディスプレイに表示されたり、プリンタで印刷されることで、歯付ベルトの騒音予測装置10の操作者に通知される。
【0042】
このように、第1の実施形態に係る歯付ベルトの騒音予測方法及び装置、並びにプログラムによれば、有限要素法解析を適用することによって、歯付ベルト1と駆動プーリ2との噛み合い点付近のベルト変位を計算し、この変位をもとに弦振動モデルの加振力を導入することができる。そして、得られた加振力を弦振動モデルに適用し、弦振動速度を算出した上で、歯付ベルト1の音響放射効率を算出し、最終的にベルト表面の放射音圧レベルを計算することができる。
【0043】
[ 実施形態2]
次に、第2の実施形態における歯付ベルトの騒音予測装置20について、図4に基づいて説明する。図4は、第2の実施形態に係る歯付ベルトの騒音予測装置のブロック線図である。尚、第1の実施形態に係る歯付ベルトの騒音予測装置10と同一の部材については、同一の符号を付記してその説明を省略する。
【0044】
第2の実施形態においては、歯付ベルトの騒音予測装置20は、ベルト変位計算部11と、時間間隔計算部12と、変動速度計算部17と、放射音圧レベル計算部18とを有している。
【0045】
ベルト変位計算部11と時間間隔計算部12は、第1の実施形態に係る歯付ベルトの騒音予測装置10と同様である。
【0046】
変動速度計算部17によって、ベルト変位計算部11から入力されたベルト変位Δyと、時間間隔演算部12から入力された1ステップ間の時間間隔Δtとにより、変動速度v=Δy/Δtが求められる。
【0047】
放射音圧レベル計算部18によって、変動速度計算部17において噛合点で求められた変動速度vを用いて、以下の式により放射音圧レベルが計算されて出力される。
【0048】
【数4】
【0049】
ここで、第1の実施形態と同様に、図4に示されている歯付ベルトの騒音予測装置20の各部11、12及び17、18は、例えば汎用のパーソナルコンピュータによって構成されている。かかるパーソナルコンピュータには、CPU、ROM、RAM、ハードディスク、FDやCDの駆動装置などのハードウェアが収納されており、ハードディスクには、パーソナルコンピュータを歯付ベルトの騒音予測プログラム(このプログラムは、CD−ROMやFD、MOなどのリムーバブルな記録媒体に記録しておくことにより、様々なコンピュータにインストールすることが可能である)を含む各種のソフトウェアが記録されている。そして、これらのハードウェアおよびソフトウェアが組み合わされることによって、上述の各部11、12及び17、18が構築されている。更に具体的には、ベルト変位計算部11は、市販の有限要素法解析プログラムが記憶されたものから構成される。
【0050】
また、歯付ベルトの騒音予測装置20の各部11、12及び17、18で得られるデータは、図示しないディスプレイに表示されたり、プリンタで印刷されたりすることで、歯付ベルトの騒音予測装置20の操作者に通知される。
【0051】
次に、図3を参照しつつ、第2の実施形態における歯付ベルトの騒音予測方法の手順について説明する。図3は、第2の実施形態における歯付ベルトの騒音予測方法の手順について説明したフローチャートである。尚、第1の実施形態に係る歯付ベルトの騒音予測方法と同一のステップについては、同一の符号を付記してその説明を省略する。
【0052】
ステップS1とステップS2は、第1の実施形態に係る歯付ベルトの騒音予測方法と同様である。
【0053】
ステップS7において、ステップS1で求めたy座標位置の変位ΔyとステップS2で求められた時間間隔Δtを用いて、各時間ステップに対するベルトのy座標の変動速度v=Δy/Δtを時間の関数として求める。
【0054】
ステップS8において、ステップS7で求められた変動速度vを用いて放射音圧レベルを計算して出力する。
【0055】
ここで、前述と同様に、各ステップS1、S2及びS7、S8は、例えば汎用のパーソナルコンピュータを用いて行われ、特にステップS1は市販の有限要素法解析プログラムを用いることができる。また各ステップS1、S2及びS7、S8で得られるデータは、図示しないディスプレイに表示されたり、プリンタで印刷されることで、歯付ベルトの騒音予測装置20の操作者に通知される。
【0056】
このように、第2の実施形態に係る歯付ベルトの騒音予測方法及び装置、並びにプログラムによれば、有限要素法解析を適用することによって、歯付ベルト1と駆動歯付プーリ2との噛み合い点付近のベルト変位を計算し、この変位をもとに簡易にベルト表面の放射音圧レベルを計算することができる。また、ベルト寸法(スパン長、幅)、ベルトの材料物性(単位長さあたりの質量、引っ張り剛性)が殆ど同じであり、且つ加振周波数が同じ場合には、ベルトの応答量は加振力にほぼ比例することから、歯形形状の変化等による発生騒音の変化を予測する場合には有効である。
【0057】
【実施例】
次に、上記のような歯付ベルト1の騒音予測方法を検証するために行った試験について説明する。図9は、試験装置に用いた歯付ベルトと歯付プーリの伝動機構の概略図である。
【0058】
図9に示すように、歯付ベルト1は、回転軸方向に沿った歯部を外周部にそれぞれ有する駆動歯付プーリ2と、アイドラ歯付プーリ3と、従動歯付プーリ4の3つの歯付プーリと噛合している。また、駆動歯付プーリ2とアイドラ歯付プーリ3との間の距離が104mm、アイドラ歯付プーリ3と従動歯付プーリ4との間の距離が300mm、従動歯付プーリ4と駆動歯付プーリ2との間の距離が440mmとなり、駆動歯付プーリ2の接触角が148°となるように構成されている。更に、駆動歯付プーリ2と、アイドラ歯付プーリ3との間にはマイクロホン5が設けられ、実際の騒音を測定することができるように構成されている。
【0059】
また、試験では、駆動機器騒音の影響を除くために計測室外に設けられたモータにより駆動歯付プーリ2を駆動している。歯付ベルト1にかかる初期張力は従動歯付プーリ4の従動軸を移動することで、負荷は従動軸に連結された負荷懸架装置によって所定の値に設定されている。
【0060】
実験条件は、ベルト歯ピッチ:8mm、ベルト幅:19mm、駆動歯付プーリ:ピッチ8mm×歯数24個、張り側張力:43.5Kgf、緩み側張力:19.6Kgf、有効張力:23.9Kgf、接触角:148°、回転数:780rpmである。
【0061】
上記の実験条件に基づいて、本発明による歯付ベルトの騒音予測方法を用いて、試験対象ベルトとプーリの歯高さ、歯幅をそれぞれ僅かづつ変化させて計8ケースの組合せについての試験を行った。また、有限要素法による解析もこれらに対応させた。
【0062】
本発明による歯付ベルトの騒音予測方法による予測結果(計算値)とマイクロホン5による計測結果(実験値)の比較を、図10(a),(b)に示している。図10(a)は、歯高さ方向の隙間(+は隙間が空き、−は圧縮となる)を横軸にとり、放射音をプロットしたものである。図10(b)は、歯幅方向隙間(バックラッシュ)を横軸にとり、放射音をプロットしたものである。なお、放射音は、歯付ベルトの噛み合い周波数での音圧レベルを採用している。
【0063】
これらのグラフから、本発明による歯付ベルトの騒音予測方法による予測結果(計算値)と計測結果(実験値)とは、実験値にバラツキが見られるものの、計算値は平均的には2〜3dB程度の精度で実験値と対応しており、略一致していることが分かる。よって、本発明による歯付ベルトの騒音予測方法が信頼でき、実用可能であることが明らかとなった。
【0064】
なお、本発明に係る歯付ベルトの騒音予測方法及び装置、並びにプログラムの実施形態は、前記実施形態に限定されるものではなく、特許請求の範囲に記載した限りにおいてさまざまな設計変更が可能である。
【0065】
【発明の効果】
本発明の歯付ベルトの騒音予測方法及び装置、並びにプログラムによると、有限要素法解析を適用することによって、噛み合い点付近に設けた定位置のベルト変位を計算し、この変位をもとに、両端を歯付プーリで支持されたベルト表面の放射音圧レベル或いは音響放射パワーレベルを計算することができる。
【0066】
従って、ベルト噛合部で生じるベルト変位に起因するベルトの弦振動によって騒音が発射されるという歯付ベルトの騒音メカニズムに基づいて、歯付ベルトの発生騒音を体系的に整理することができる。そして、歯付ベルトの発生騒音を事前に予測することができることにより、歯付ベルトの低騒音化が可能となり、環境騒音の改善に貢献することができる。また、歯付ベルトの発生騒音の事前予測が可能となることにより、発生騒音の評価や低騒音化対策のための試験が不要となる、或いは簡素化されることで人的、物的な資源の節約となる。
【図面の簡単な説明】
【図1】第1の実施形態における歯付ベルトの騒音予測方法の手順について説明したフローチャートである。
【図2】第1の実施形態に係る歯付ベルトの騒音予測装置のブロック線図である。
【図3】第2の実施形態における歯付ベルトの騒音予測方法の手順について説明したフローチャートである。
【図4】第2の実施形態に係る歯付ベルトの騒音予測装置のブロック線図である。
【図5】本実施形態に係る歯付ベルトと歯付プーリの伝動機構の概略図である。
【図6】本実施形態に係る歯付ベルトと駆動歯付プーリの断面略図である。
【図7】噛み合いステップと歯付ベルトのy座標位置の関係を示した図である。
【図8】ベルトをマスとバネとダッシュポットからなる質点系に置き換えた図である。
【図9】試験装置に用いた歯付ベルトと歯付プーリの伝動機構の概略図である。
【図10】本発明による歯付ベルトの騒音予測結果と試験結果の比較であり、(a)は横軸に歯高さ方向の隙間を、(b)は横軸に歯幅方向の隙間をとったものである。
【符号の説明】
1 歯付ベルト
2 駆動歯付プーリ
3 アイドラ歯付プーリ
4 従動歯付プーリ
10 歯付プーリの騒音予測装置
11 ベルト変位計算部(ベルト変位計算手段)
13 変位、或いは変動速度計算部(変位、或いは変動速度計算手段)
14 ベルト弦振動計算部(ベルト弦振動計算手段)
15 音響放射効率計算部(音響放射効率計算手段)
16 放射音圧レベルまたは音響放射パワーレベル計算部(放射音圧レベルまたは音響放射パワーレベル計算手段)
17 変動速度計算部(変動速度計算手段)
18 放射音圧レベル計算部(放射音圧レベル計算手段)
20 歯付プーリの騒音予測装置
Claims (12)
- 歯付ベルトの両端を歯付プーリで支持されたベルトスパンのベルト噛み合い時に、歯付プーリの多角形効果及び歯付ベルトのベルト噛合い点近傍でのベルトの変形により生じるベルト走行方向に対して直角方向のベルト変位に起因するベルトの弦振動により発射される騒音を予測する歯付ベルトの騒音予測方法であって、
ベルト噛合い点近傍で生じるベルト走行方向に対して直角方向のベルト変位を計算するステップと、
前記ベルト変位を、プーリ回転数とプーリ歯数等の諸元より、時間の関数に変換し、変動速度を計算するステップと、
前記変動速度及び前記ベルトの寸法から、前記ベルトの弦振動による放射音圧レベルを計算するステップとを備えることを特徴とする歯付ベルトの騒音予測方法。 - 歯付ベルトの両端を歯付プーリで支持されたベルトスパンのベルト噛み合い時に、歯付プーリの多角形効果及び歯付ベルトのベルト噛合い点近傍でのベルトの変形により生じるベルト走行方向に対して直角方向のベルト変位に起因するベルトの弦振動により発射される騒音を予測する歯付ベルトの騒音予測方法であって、ベルト噛合い点近傍で生じるベルト走行方向に対して直角方向のベルト変位を計算するステップと、
前記ベルト変位を、プーリ回転数とプーリ歯数等の諸元より時間の関数に変換し、時間、或いは周波数の関数としての変位、或いは変動速度を計算するステップと、
前記変換された変位、或いは変動速度を、前記ベルトスパンの噛合端の加振力として、ベルト弦振動の応答変位、もしくは振動モード、及びベルト弦振動の速度を計算するステップと、
前記ベルト弦振動の応答変位、もしくは振動モード、及び前記ベルトの寸法から、前記ベルトの弦振動による音響放射効率を計算するステップと、
前記ベルト弦振動の速度、及び前記音響放射効率から、放射音圧レベルまたは音響放射パワーレベルを計算するステップとを備えることを特徴とする歯付ベルトの騒音予測方法。 - 前記ベルト変位は、有限要素法解析を適用して計算されることを特徴とする請求項1または2に記載の歯付ベルトの騒音予測方法。
- 前記有限要素法解析は、前記ベルト変位を、ベルトに固定した点の噛み合い進行に伴う噛合部までのベルト走行方向に対して直角方向の変化として計算することを特徴とする請求項3に記載の歯付ベルトの騒音予測方法。
- 歯付ベルトの両端を歯付プーリで支持されたベルトスパンのベルト噛み合い時に、歯付プーリの多角形効果及び歯付ベルトのベルト噛合い点近傍でのベルトの変形により生じるベルト走行方向に対して直角方向のベルト変位に起因するベルトの弦振動により発射される騒音を予測する歯付ベルトの騒音予測装置であって、
ベルト噛合い点近傍で生じるベルト走行方向に対して直角方向のベルト変位を計算するベルト変位計算手段と、
前記ベルト変位を、プーリ回転数とプーリ歯数等の諸元より、時間の関数に変換し、変動速度を計算する変動速度計算手段と、
前記変動速度及び前記ベルトの寸法から、前記ベルトの弦振動による放射音圧レベルを計算する放射音圧レベル計算手段とを備えることを特徴とする歯付ベルトの騒音予測装置。 - 歯付ベルトの両端を歯付プーリで支持されたベルトスパンのベルト噛み合い時に、歯付プーリの多角形効果及び歯付ベルトのベルト噛合い点近傍でのベルトの変形により生じるベルト走行方向に対して直角方向のベルト変位に起因するベルトの弦振動により発射される騒音を予測する歯付ベルトの騒音予測装置であって、
ベルト噛合い点近傍で生じるベルト走行方向に対して直角方向のベルト変位を計算するベルト変位計算手段と、
前記ベルト変位を、プーリ回転数とプーリ歯数等の諸元より、時間の関数に変換し、時間、或いは周波数の関数としての変位、或いは変動速度を計算する変位、或いは変動速度計算手段と、前記変換された変位、或いは変動速度を、前記ベルトスパンの噛合端の加振力として、ベルト弦振動の応答変位、もしくは振動モード、及びベルト弦振動の速度を計算するベルト弦振動計算手段と、
前記ベルト弦振動の応答変位、もしくは振動モード、及び前記ベルトの寸法から、前記ベルトの弦振動による音響放射効率を計算する音響放射効率計算手段と、
前記ベルト弦振動の速度、及び前記音響放射効率から、放射音圧レベルまたは音響放射パワーレベルを計算する放射音圧レベルまたは音響放射パワーレベル計算手段とを備えることを特徴とする歯付ベルトの騒音予測装置。 - 前記ベルト変位計算手段は、有限要素法解析を適用して前記ベルト変位を計算することを特徴とする請求項5または6に記載の歯付ベルトの騒音予測装置。
- 前記有限要素法解析は、前記ベルト変位を、ベルトに固定した点の噛み合い進行に伴う噛合部までのベルト走行方向に対して直角方向の変化として計算することを特徴とする請求項7に記載の歯付ベルトの騒音予測装置。
- 歯付ベルトの両端を歯付プーリで支持されたベルトスパンのベルト噛み合い時に、歯付プーリの多角形効果及び歯付ベルトのベルト噛合い点近傍でのベルトの変形により生じるベルト走行方向に対して直角方向のベルト変位に起因するベルトの弦振動により発射される騒音を予測する歯付ベルトの騒音予測プログラムであって、
ベルト噛合い点近傍で生じるベルト走行方向に対して直角方向のベルト変位を計算するベルト変位計算手段、
前記ベルト変位を、プーリ回転数とプーリ歯数等の諸元より、時間の関数に変換して、変動速度を計算する変動速度計算手段、
前記変動速度及び前記ベルトの寸法から、前記ベルトの弦振動による放射音圧レベルを計算する放射音圧レベル計算手段、としてコンピュータを機能させるための歯付ベルトの騒音予測プログラム。 - 歯付ベルトの両端を歯付プーリで支持されたベルトスパンのベルト噛み合い時に、歯付プーリの多角形効果及び歯付ベルトのベルト噛合い点近傍でのベルトの変形により生じるベルト走行方向に対して直角方向のベルト変位に起因するベルトの弦振動により発射される騒音を予測する歯付ベルトの騒音予測プログラムであって、
ベルト噛合い点近傍で生じるベルト走行方向に対して直角方向のベルト変位を計算するベルト変位計算手段、
前記ベルト変位を、プーリ回転数とプーリ歯数等の諸元より、時間の関数に変換し、時間、或いは周波数の関数としての変位、或いは変動速度を計算する変位、或いは変動速度計算手段、
前記変換された変位、或いは変動速度を、前記ベルトスパンの噛合端の加振力として、ベルト弦振動の応答変位、もしくは振動モード、及びベルト弦振動の速度を計算するベルト弦振動計算手段、
前記ベルト弦振動の応答変位、もしくは振動モード、及び前記ベルトの寸法から、前記ベルトの弦振動による音響放射効率を計算する音響放射効率計算手段、
前記ベルト弦振動の速度、及び前記音響放射効率から、放射音圧レベルまたは音響放射パワーレベルを計算する放射音圧レベルまたは音響放射パワーレベル計算手段、としてコンピュータを機能させるための歯付ベルトの騒音予測プログラム。 - 前記ベルト変位計算手段は、有限要素法解析を適用して前記ベルト変位を計算することを特徴とする請求項9または10に記載の歯付ベルトの騒音予測プログラム。
- 前記有限要素法解析は、前記ベルト変位を、ベルトに固定した点の噛み合い進行に伴う噛合部までのベルト走行方向に対して直角方向の変化として計算することを特徴とする請求項11に記載の歯付ベルトの騒音予測プログラム。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002050863A JP3686871B2 (ja) | 2002-02-27 | 2002-02-27 | 歯付ベルトの騒音予測方法及び装置、並びにプログラム |
| US10/374,843 US6882945B2 (en) | 2002-02-27 | 2003-02-26 | Method, apparatus, and program for estimating noise generation for a synchronous belt |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002050863A JP3686871B2 (ja) | 2002-02-27 | 2002-02-27 | 歯付ベルトの騒音予測方法及び装置、並びにプログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2003247916A JP2003247916A (ja) | 2003-09-05 |
| JP3686871B2 true JP3686871B2 (ja) | 2005-08-24 |
Family
ID=28662983
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002050863A Expired - Fee Related JP3686871B2 (ja) | 2002-02-27 | 2002-02-27 | 歯付ベルトの騒音予測方法及び装置、並びにプログラム |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US6882945B2 (ja) |
| JP (1) | JP3686871B2 (ja) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4839242B2 (ja) * | 2007-02-15 | 2011-12-21 | 株式会社リコー | 転写ベルト駆動機構のモデル化方法及びそのモデル化方法を用いた制御系の設計支援方法 |
| JP4995687B2 (ja) * | 2007-10-11 | 2012-08-08 | 三ツ星ベルト株式会社 | ベルト鳴きの評価方法及び装置、並びにプログラム |
| JP5676244B2 (ja) * | 2010-12-28 | 2015-02-25 | 株式会社ブリヂストン | 騒音予測装置及び騒音予測方法 |
| JP6165043B2 (ja) * | 2013-12-13 | 2017-07-19 | 三菱重工業株式会社 | 騒音評価装置及び騒音評価方法 |
| CN106154216B (zh) * | 2015-04-03 | 2019-05-07 | 北京智谷睿拓技术服务有限公司 | 牙齿定位方法和设备 |
| CN106154215B (zh) * | 2015-04-03 | 2019-05-07 | 北京智谷睿拓技术服务有限公司 | 牙齿定位方法和设备 |
| CN113190912B (zh) * | 2021-03-12 | 2022-10-28 | 西安交通大学 | 舰艇振动壳体的辐射噪声预测方法 |
| CN113492406B (zh) * | 2021-07-19 | 2022-07-12 | 中国人民解放军92578部队 | 一种舱室噪声声压自动检测方法 |
Family Cites Families (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3881364A (en) * | 1974-03-22 | 1975-05-06 | William S Rouverol | Low noise gearing |
| JPS5724997Y2 (ja) * | 1978-10-30 | 1982-05-31 | ||
| US4416649A (en) * | 1979-09-14 | 1983-11-22 | Uniroyal, Inc. | Reduced noise positive drive power transmission system |
| US4650445A (en) * | 1985-02-20 | 1987-03-17 | Borg-Warner Automotive, Inc. | Chain-belt |
| US4824424A (en) * | 1987-08-24 | 1989-04-25 | Fuji Jukogyo Kabushiki Kaisha | Belt for a belt drive device |
| US5405299A (en) * | 1992-12-29 | 1995-04-11 | Bridgestone Corporation | Toothed belt and toothed belt-adapted pulley |
| IT1265351B1 (it) * | 1993-11-26 | 1996-11-22 | Dayco Pti Spa | Trasmissione a cinghia dentata per il comando di alberi di equilibratura dinamica in motori endotermici |
| JP2960648B2 (ja) * | 1994-06-13 | 1999-10-12 | バンドー化学株式会社 | ベルトの荷重分担予測方法 |
| US5554807A (en) * | 1995-04-06 | 1996-09-10 | Tranquilla; Michael N. | Belt damping |
| JPH0949786A (ja) * | 1995-08-08 | 1997-02-18 | Toyota Central Res & Dev Lab Inc | 歯付きベルト伝動機構の騒音評価方法と騒音評価装置 |
| JPH09125558A (ja) * | 1995-08-31 | 1997-05-13 | Nippon Steel Corp | 音響減衰特性に優れた金属薄板構造体 |
| JPH09133592A (ja) * | 1995-11-08 | 1997-05-20 | Bando Chem Ind Ltd | 歯付ベルトの寿命予測方法及び寿命予測装置 |
| US5798456A (en) * | 1996-11-19 | 1998-08-25 | Unisys Corp. | Predicting behavior of synchronous flexible webs |
| US5940788A (en) * | 1997-05-02 | 1999-08-17 | Ford Global Technologies, Inc. | Method and system for designing vehicle door seals based on predicted sound transmission characteristics |
| GB0016259D0 (en) * | 2000-07-03 | 2000-08-23 | Lms International | A computer-aided engineering method for predicting the acoustic signature of vibrating structures using discrete mode |
| US20020101135A1 (en) * | 2000-07-28 | 2002-08-01 | Marco Giovanardi | Method and device for noise damping |
| JP3698625B2 (ja) * | 2000-09-08 | 2005-09-21 | バンドー化学株式会社 | 伝動ベルト |
| US6568280B1 (en) * | 2000-11-06 | 2003-05-27 | Mechanical Dynamics, Inc. | Method of computing belt tension distribution of metal-pushing V-belt for CVT |
-
2002
- 2002-02-27 JP JP2002050863A patent/JP3686871B2/ja not_active Expired - Fee Related
-
2003
- 2003-02-26 US US10/374,843 patent/US6882945B2/en not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| US6882945B2 (en) | 2005-04-19 |
| JP2003247916A (ja) | 2003-09-05 |
| US20030225545A1 (en) | 2003-12-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3686871B2 (ja) | 歯付ベルトの騒音予測方法及び装置、並びにプログラム | |
| Lyon | Machinery noise and diagnostics | |
| Andersson et al. | A dynamic model to determine vibrations in involute helical gears | |
| Yu et al. | A new dynamic model of a cylindrical gear pair with localized spalling defects | |
| He et al. | Novel mathematical modelling method for meshing impact of helical gear | |
| US6789035B1 (en) | Method of analyzing vibration in a transmission belt and an apparatus and program useable to perform the method | |
| KR20150080525A (ko) | 벨트 장력 산출 프로그램 및 벨트 고유 주파수 산출 프로그램, 그리고 이들의 방법 및 장치 | |
| Bortnowski et al. | Modelling of transverse vibration of conveyor belt in aspect of the trough angle | |
| Tanaka et al. | Vibration and sound-radiation analysis for designing a low-noise gearbox with a multi-stage helical gear system | |
| JPWO2017145850A1 (ja) | 検査装置、検査方法、及び、検査プログラム | |
| JP5792705B2 (ja) | 歪分布測定装置、歪分布測定方法、および、歪分布測定プログラム | |
| Nagai et al. | Experiments and analysis on chaotic vibrations of a shallow cylindrical shell-panel | |
| JP2009294729A (ja) | 歯車設計支援方法、歯車設計支援装置、歯車設計支援プログラムおよび記憶媒体 | |
| JP2005156394A (ja) | 伝動ベルトの振動解析方法及び装置、並びにプログラム | |
| RU2728584C1 (ru) | Способ исключения резонансных режимов в фрикционных дисках гидромеханической трансмиссии транспортной машины | |
| JP2003049905A (ja) | 有限要素法解析を用いた伝動ベルトの荷重分担予測方法及び装置、並びにプログラム | |
| Kagotani et al. | Transmission error in synchronous belt with resonance under installation tension | |
| CN119004921A (zh) | 齿轮振动噪声评价方法及装置、存储介质、计算机设备 | |
| JP5538942B2 (ja) | 放射音予測方法 | |
| CN116163914B (zh) | 一种空调器压缩机预警方法、装置和系统 | |
| Savant et al. | Numerical investigation of annoying CVT belt noise | |
| JP5020860B2 (ja) | 振動低減部材 | |
| JP2001311456A (ja) | 伝動ベルト駆動系の振動解析方法及び設計支援方法 | |
| Chen et al. | Noise Modeling of Synchronous Belts | |
| JP2002257797A (ja) | 軸受損傷評価装置及び軸受損傷評価方法及び軸受損傷評価プログラム及びこのプログラムを記録した記憶媒体 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20040802 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20040817 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20041013 |
|
| 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: 20050524 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050606 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 3686871 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090610 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090610 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100610 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100610 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110610 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120610 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130610 Year of fee payment: 8 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |
