JP3823397B2 - アクティブ制振方法 - Google Patents
アクティブ制振方法 Download PDFInfo
- Publication number
- JP3823397B2 JP3823397B2 JP28693996A JP28693996A JP3823397B2 JP 3823397 B2 JP3823397 B2 JP 3823397B2 JP 28693996 A JP28693996 A JP 28693996A JP 28693996 A JP28693996 A JP 28693996A JP 3823397 B2 JP3823397 B2 JP 3823397B2
- Authority
- JP
- Japan
- Prior art keywords
- frequency
- vibration
- signal component
- mode
- actuator
- 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
Images
Landscapes
- Feedback Control In General (AREA)
- Vibration Prevention Devices (AREA)
Description
【0001】
【発明の属する技術分野】
本発明は、対象物の振動を打ち消すように積極的に加振を行うアクティブ制振の方法に係り、特に、制振帯域外に生じる暴れを解消するアクティブ制振方法に関するものである。
【0002】
【従来の技術】
建造物、機械、移動体、流体等の対象物に好ましくない振動が加わるか又は生じている場合に、この対象物の振動を打ち消すように積極的に加振を行い静止を図る、いわゆるアクティブ制振の技術が知られている。
【0003】
アクティブ制振にあっては、対象物に、振動を検出するセンサと対象物を加振するアクチュエータとを取り付け、センサの検出信号に応じてアクチュエータを駆動するようになっている。
【0004】
対象物の振動はその固有の振動特性により顕著に振動する振動周波数が飛び飛びに存在し、これらの周波数における振動を周波数の低い方より順に第1次,第2次,…第m次モード(mは整数)と呼ぶ。一般に周波数の低い振動モードが大きく作用するので、主として周波数の低い振動モードについて制振を行うことにより、対象物をほぼ静止させることができる。
【0005】
一方、アクチュエータにも駆動制御が可能な周波数に限りがあるため、上記対象物の振動特性に照らし合わせ、アクチュエータの能力を予め設定することが重要である。
【0006】
従って、対象物の振動モードのうち周波数の低い振動モード、即ち第1次モードより必要な第m次モードまでが含まれる振動周波数の範囲において、アクチュエータ駆動による制振が有効であればよいことになる。
【0007】
従来のアクティブ制振方法では、このような周波数の低い低次振動モードに対応するべく、センサの検出信号を低域通過フィルタに通し、この低域信号成分に応じてアクチュエータを駆動制御するようになっている。
【0008】
このような低域通過フィルタにはデジタルフィルタが一般に用いられる。
【0009】
低域通過フィルタの周波数特性は、第m次モードの周波数以下では利得が1、第m+1次モードの周波数以上では利得が0であることを理想とする。しかし、実際には急峻な遮断特性を得ようとすると位相が大きくずれたり、利得が平坦にならなかったりするので、第m次モードの周波数より下から第m+1次モードの周波数より上までに亘ってなだらかに利得が小さくなるような遮断特性を用いている(図9参照)。
【0010】
【発明が解決しようとする課題】
従来技術にあっては、アクティブ制振のためのアクチュエータの駆動制御は、設定された振動周波数の範囲(以下、これを制振帯域という)において最も有効に働くようになっており、その半面で、制振帯域外については必ずしも制振に有効でない。例えば、第2次モードまでを対象とするシステムでは第3次モードから上には有効でない。第3次モードまでを対象に広げたものは、今度は第4次モードから上には有効でない。即ち、第m次モードを対象とするシステムでは第m+1次モード(これを打ち切りモードという)から上には有効でない。
【0011】
従って、図7に示されるように、制振を行わない場合(無制御;破線)に、各次モードの周波数毎に同じ大きさの応答振幅があるものとしたとき、制振を行った場合(制御;実線)には、各次モードの周波数における応答振幅は第m次モードの周波数ωm までについては確かに小さくなっている。
【0012】
これに対し、低域通過フィルタの周波数特性は、第m次モードの周波数以下から第m+1次モードの周波数以上までに亘るなだらかな遮断特性となっているので、アクチュエータの駆動制御の基となる信号に高い周波数の成分が含まれることになる。アクチュエータは、これらの高周波数成分を含む駆動を強いられ、不安定になる。このため、図7に示されるように、制振を行った場合(制御;実線)において、打ち切りモードである第m+1次モードの周波数ωm+1 で、かえって振動が激しくなるという不具合(これをスピルオーバという)が生じる。
【0013】
このように従来のアクティブ制振方法は、第m+1次モードにスピルオーバが生じ、しかも、第m次モードの周波数以下を極力抑えようとすると、ますますスピルオーバが激しくなるというジレンマがある。
【0014】
そこで、本発明の目的は、上記課題を解決し、図8に示されるように、第m+1次モードの周波数ωm+1 において応答振幅を無制御のときと同じ程度まで抑えること、即ち、制振帯域外に生じる暴れを解消するアクティブ制振方法を提供することにある。
【0015】
【課題を解決するための手段】
上記目的を達成するために本発明は、対象物の振動を検出し、その検出信号に応じてアクチュエータを駆動することにより、対象物の振動を打ち消すアクティブ制振の方法において、上記アクチュエータ駆動による制振が有効な振動周波数の範囲全域で上記検出信号を遮断し、それより高い周波数では通過させる高域通過フィルタを用い、この高域通過フィルタを現在時間までの検出信号に適用して所定時間前の高域信号成分を抽出し、この高域信号成分の波形を現在時間まで外挿し、この外挿による現在時間の高域信号成分を上記検出信号から差し引くことにより上記アクチュエータ駆動による制振が有効な低域信号成分を検出し、その低域信号成分に応じて上記アクチュエータを駆動するものである。
【0016】
上記高域通過フィルタは、上記アクチュエータ駆動による制振が有効な振動周波数の範囲に上記対象物の振動モードのうち第1次モードより第m次モードまでが含まれるとき、第m次モードの周波数以下では利得が0、第m+1次モードの周波数以上では利得が1の周波数特性を有し、上記検出信号から第m+1次モードの周波数以上の高域信号成分を抽出する。
【0017】
上記高域通過フィルタは、予めその周波数特性をフーリエ逆変換して得られる時間関数と過去の所定時間t−T(tは現時刻を表す)より未来の所定時間t+Tまでの有限時間とを設定してなり、この時間関数と検出信号とを過去の所定時間t−2Tより現在時間まで畳み込むことにより、過去の所定時間t−Tにおける高域信号成分を出力する。
【0018】
上記高域信号成分の外挿波形は第m+1次モードの周波数を持つ正弦波である。
【0019】
【発明の実施の形態】
以下、本発明の一実施形態を添付図面に基づいて詳述する。
【0020】
図1に示されたモデルは、比較的柔軟な構造物を対象物としたものである。対象物1は細い柔軟な梁である。対象物1の一端は壁2に固定され、他端は自由となっている。この自由端を加振点3とする。この加振点3はこれから制振しようとする振動をこのモデルに与えるために設けた。対象物1の適宜な位置にアクチュエータが取り付けられており、その位置を制御点4という。また、振動を検出するセンサが取り付けられており、その位置をセンシング点5という。このモデルでは対象物1は壁2から加振点3までの150cmの長さがある。ここで用いられるセンサは速度センサであり、壁2から104cmの位置に制御点4及び速度センシング点5が設けられている。なお、対象物1は厚さ2mm、ヤング率2.0×1011N/m2 、ポアソン比0.3、各振動モードに減衰比0.01を設定してある。また、制御点4、加振点3には力発生用ボイスコイル、センサによる質量175g,100gがそれぞれ付加される。
【0021】
図2に示されるように、アクティブ制振装置は、対象物である柔軟構造体21に取り付けられたセンサ22と、センサ22の検出信号を所定周波数でサンプリングしてデジタルデータ列に変換するA/D変換器23と、そのデジタルデータ列を処理してアクチュエータ駆動指令値列を生成するパソコン等の計算機24と、そのアクチュエータ駆動指令値列を上記サンプリング周波数でアナログ波形信号に変換するD/A変換器25と、このアナログ波形信号で駆動制御され、柔軟構造体21を加振するアクチュエータ26とからなる。センサ22は柔軟構造体21の振動を検出するものであり、その検出量は変位、速度、加速度のいずれでもよい。これに対し、アクチュエータ26は力、圧力等を働かせるものである。アクチュエータ26は、柔軟構造体21外の図示されない固定点で支持されたものでもよいし、柔軟構造体21上でマスを移動させ反力を取り出すようなものでもよい。
【0022】
本発明に用いるために計算機24に組み込まれた高域通過フィルタの周波数特性を図3に示す。図示されるように、対象物の振動モードの第m次モードの周波数ωm より低周波では利得が0、かつ位相遅れが0であり、第m+1次モードの周波数ωm+1 より高周波では利得が1、かつ位相遅れが0である。周波数ωm からωm+1 までの領域における利得の変化は任意である。
【0023】
この特性を周波数伝達関数F(ω)で示すと、
F(ω)=0 0<ω≦ωm
F(ω)=1 ωm+1 ≦ω
であり、図4に示されるように、この高域通過フィルタ41において、センサの検出信号をy(t)、フィルタ出力をy1 (t)とし、積分範囲をt−∞≦τ≦t+∞とする代わりに、t−T≦τ≦t+Tの範囲のみで計算すると、
【0024】
【数1】
【0025】
となる。このフィルタ出力y1 (t)が高域信号成分である。
【0026】
ここで、h(t)は、周波数伝達関数F(ω)のフーリエ逆変換として定まる時間関数即ちインパルス応答関数である。即ち、
【0027】
【数2】
【0028】
である。この式(2)の計算は複雑かつ計算ステップ数が膨大であるが、フィルタ固有のものであるから、予め大型計算機等の上記計算機24とは別の計算機で計算しておき、その結果を計算機24に組み込み、計算機24では式(1)を計算するようにしてもよい。
【0029】
インパルス応答関数h(t)の一例を図5に示す。インパルス応答関数h(t)は、t=0で最大値を示し、他にもいくつかのピークがある。図には描ききれないが、インパルス応答関数h(t)は、過去無限大時間から未来無限大時間まで存在する。しかし、実用上は、このような無限大時間のデータは取り扱えないので、過去の所定時間−Tより未来の所定時間Tまでの有限時間を設定する。所定時間Tの決め方としては、例えばt=0における最大値に対しh(t)の絶対値が数%の値となるピークまでを含むものとする。この範囲外は無限大時間まで存在しても値が小さいので無視し得る。
【0030】
式(1)のようにインパルス応答関数h(t)と検出信号y(t)とを−T≦t≦Tの範囲で畳み込むことによりフィルタ出力y1 (t)が得られる。有限時間を設定してあるので、計算ステップ数が少なく、従って、パソコン等の計算機24によるリアルタイム処理が可能である。しかし、積分範囲に未来の時間が含まれており、未来の時間の入力信号は未知であるため、このままでは実用できない。
【0031】
そこで、式(1)の適用を時間Tだけ過去にシフトする。即ち、インパルス応答関数と検出信号とを過去の所定時間t−2Tより現在時間まで畳み込むことにより、過去の所定時間t−Tにおけるフィルタ出力を計算することにする。
【0032】
【数3】
【0033】
この式(3)によるフィルタ出力y1 (t−T)は、過去の所定時間t−Tにおける高域信号成分である。図6に示されるように、現在時間tまでの検出信号y(t)は既知であり、この検出信号y(t)を用いて、時間t−Tまでのフィルタ出力61を計算することは可能である。このようにして過去の所定時間t−Tにおける高域信号成分y1 (t−T)が求まる。
【0034】
しかし、求めたい値は現在時間のものである。そこで、フィルタ出力61の時間的変化を調べることにより、その概ねの周期、振幅が判る。これに従って、フィルタ出力62を現在時間まで外挿し、現在時間の値を得る。具体的には、y1 (τ)(τ≦t−T)を周波数ωm+1 を持つ正弦波又は余弦波等の調和関数で表す。この調和関数を時間Tだけ延長することにより、y1 (t)を容易に得ることができる。これは、カーブフィットと呼ばれる手法であり、特に柔軟な構造物にあっては実際の高域信号成分と殆ど同じものが得られる。
【0035】
この外挿による現在時間の高域信号成分を検出信号から差し引くことにより、低域信号成分のみを検出することができる。即ち、式(4)により低域信号成分ylow (t)が求まる。
【0036】
【数4】
【0037】
図2のアクティブ制振装置では、センサ22の検出信号はA/D変換器23により所定周波数でサンプリングしてデジタルデータ列yi に変換されているので、計算機24では、
【0038】
【数5】
【0039】
を計算することになる。式(5)のAjは、式(3)及びその外挿部分に関して予めオフラインで求めた定数である。従って、デジタルデータ列yi から簡単に低域信号成分のみからなるデジタルデータ列ylow,i を計算できる。
【0040】
計算機24では、低域信号成分のみからなるデジタルデータ列ylow,i に基づき、対象物である柔軟構造体21の振動を打ち消すように各サンプリング時間毎のアクチュエータの指令値が計算される。このようにして得られたアクチュエータ駆動指令値列がD/A変換器25によりアナログ波形信号に変換される。このアナログ波形信号でアクチュエータを駆動する。
【0041】
アクチュエータ駆動指令値の計算の基になった低域信号成分ylow (t)は、図3の特性を有する高域通過フィルタの出力を検出信号から差し引いたものであるから、図9に実線で示した特性を有する低域通過フィルタの出力に等しい。従って、計算機24で行われるアクチュエータ駆動指令値の計算には第m+1次モードの周波数以上の成分が関与しなくなり、アクチュエータはその能力で十分対応しきれる駆動がなされるようになり、スピルオーバが生じない。これにより、図8に示されるように、第m+1次モードの周波数ωm+1 において応答振幅が無制御のときと同じ程度まで抑えることができる。また、制御ゲインを上げてもスピルオーバが小さくなる。
【0042】
図1に示されたモデルに対し、第2次モードまでを制振帯域としたときの結果を説明する。図10及び図11は、横軸に加振点3に加えた加振周波数、縦軸に(センシング点5における速度振幅)/(加振点3に加えた加振力振幅)を取ってある。図10は従来技術、図11は本発明の結果を示す。
【0043】
図10に示されるように、無制御のとき、第1次モード,第2次モード,第3次モードの振動ピークが存在することが判る。第2次モードまでを対象に従来技術によるアクティブ制振を行うと、第3次モードでは振動が15dB(6倍)も強調されてしまう。
【0044】
図11に示されるように、第2次モードまでを対象に本発明によるアクティブ制振を行うと、第1次モード,第2次モードの振動が小さくなることは勿論、第3次モードは無制御のときと同程度にとどまっている。
【0045】
【発明の効果】
本発明は次の如き優れた効果を発揮する。
【0046】
(1)第m+1次モードのスピルオーバがなくなり、対象物が目論見通り制振できるようになる。
【0047】
(2)アクチュエータが高い周波数に追従しなくてもよいので、簡素で安価なものを使用することができる。
【図面の簡単な説明】
【図1】本発明の一実施形態を示す振動体のモデル図である。
【図2】本発明に係るアクティブ制振装置のブロック図である。
【図3】本発明に用いる高域通過フィルタの周波数特性図である。
【図4】図3の特性を持つ高域通過フィルタのブロック図である。
【図5】図3の特性を持つ高域通過フィルタのインパルス応答関数の波形図である。
【図6】本発明の高域通過フィルタの出力算出方法を示す時間波形図である。
【図7】従来のアクティブ制振方法による応答振幅特性図である。
【図8】本発明の目的とする応答振幅特性図である。
【図9】従来及び本発明の目的とする低域通過フィルタの周波数特性図である。
【図10】モデル実験の結果を示す従来のアクティブ制振の周波数特性図である。
【図11】モデル実験の結果を示す本発明のアクティブ制振の周波数特性図である。
【符号の説明】
1 対象物
2 壁
3 加振点
4 制御点
5 センシング点
41 高域通過フィルタ
【発明の属する技術分野】
本発明は、対象物の振動を打ち消すように積極的に加振を行うアクティブ制振の方法に係り、特に、制振帯域外に生じる暴れを解消するアクティブ制振方法に関するものである。
【0002】
【従来の技術】
建造物、機械、移動体、流体等の対象物に好ましくない振動が加わるか又は生じている場合に、この対象物の振動を打ち消すように積極的に加振を行い静止を図る、いわゆるアクティブ制振の技術が知られている。
【0003】
アクティブ制振にあっては、対象物に、振動を検出するセンサと対象物を加振するアクチュエータとを取り付け、センサの検出信号に応じてアクチュエータを駆動するようになっている。
【0004】
対象物の振動はその固有の振動特性により顕著に振動する振動周波数が飛び飛びに存在し、これらの周波数における振動を周波数の低い方より順に第1次,第2次,…第m次モード(mは整数)と呼ぶ。一般に周波数の低い振動モードが大きく作用するので、主として周波数の低い振動モードについて制振を行うことにより、対象物をほぼ静止させることができる。
【0005】
一方、アクチュエータにも駆動制御が可能な周波数に限りがあるため、上記対象物の振動特性に照らし合わせ、アクチュエータの能力を予め設定することが重要である。
【0006】
従って、対象物の振動モードのうち周波数の低い振動モード、即ち第1次モードより必要な第m次モードまでが含まれる振動周波数の範囲において、アクチュエータ駆動による制振が有効であればよいことになる。
【0007】
従来のアクティブ制振方法では、このような周波数の低い低次振動モードに対応するべく、センサの検出信号を低域通過フィルタに通し、この低域信号成分に応じてアクチュエータを駆動制御するようになっている。
【0008】
このような低域通過フィルタにはデジタルフィルタが一般に用いられる。
【0009】
低域通過フィルタの周波数特性は、第m次モードの周波数以下では利得が1、第m+1次モードの周波数以上では利得が0であることを理想とする。しかし、実際には急峻な遮断特性を得ようとすると位相が大きくずれたり、利得が平坦にならなかったりするので、第m次モードの周波数より下から第m+1次モードの周波数より上までに亘ってなだらかに利得が小さくなるような遮断特性を用いている(図9参照)。
【0010】
【発明が解決しようとする課題】
従来技術にあっては、アクティブ制振のためのアクチュエータの駆動制御は、設定された振動周波数の範囲(以下、これを制振帯域という)において最も有効に働くようになっており、その半面で、制振帯域外については必ずしも制振に有効でない。例えば、第2次モードまでを対象とするシステムでは第3次モードから上には有効でない。第3次モードまでを対象に広げたものは、今度は第4次モードから上には有効でない。即ち、第m次モードを対象とするシステムでは第m+1次モード(これを打ち切りモードという)から上には有効でない。
【0011】
従って、図7に示されるように、制振を行わない場合(無制御;破線)に、各次モードの周波数毎に同じ大きさの応答振幅があるものとしたとき、制振を行った場合(制御;実線)には、各次モードの周波数における応答振幅は第m次モードの周波数ωm までについては確かに小さくなっている。
【0012】
これに対し、低域通過フィルタの周波数特性は、第m次モードの周波数以下から第m+1次モードの周波数以上までに亘るなだらかな遮断特性となっているので、アクチュエータの駆動制御の基となる信号に高い周波数の成分が含まれることになる。アクチュエータは、これらの高周波数成分を含む駆動を強いられ、不安定になる。このため、図7に示されるように、制振を行った場合(制御;実線)において、打ち切りモードである第m+1次モードの周波数ωm+1 で、かえって振動が激しくなるという不具合(これをスピルオーバという)が生じる。
【0013】
このように従来のアクティブ制振方法は、第m+1次モードにスピルオーバが生じ、しかも、第m次モードの周波数以下を極力抑えようとすると、ますますスピルオーバが激しくなるというジレンマがある。
【0014】
そこで、本発明の目的は、上記課題を解決し、図8に示されるように、第m+1次モードの周波数ωm+1 において応答振幅を無制御のときと同じ程度まで抑えること、即ち、制振帯域外に生じる暴れを解消するアクティブ制振方法を提供することにある。
【0015】
【課題を解決するための手段】
上記目的を達成するために本発明は、対象物の振動を検出し、その検出信号に応じてアクチュエータを駆動することにより、対象物の振動を打ち消すアクティブ制振の方法において、上記アクチュエータ駆動による制振が有効な振動周波数の範囲全域で上記検出信号を遮断し、それより高い周波数では通過させる高域通過フィルタを用い、この高域通過フィルタを現在時間までの検出信号に適用して所定時間前の高域信号成分を抽出し、この高域信号成分の波形を現在時間まで外挿し、この外挿による現在時間の高域信号成分を上記検出信号から差し引くことにより上記アクチュエータ駆動による制振が有効な低域信号成分を検出し、その低域信号成分に応じて上記アクチュエータを駆動するものである。
【0016】
上記高域通過フィルタは、上記アクチュエータ駆動による制振が有効な振動周波数の範囲に上記対象物の振動モードのうち第1次モードより第m次モードまでが含まれるとき、第m次モードの周波数以下では利得が0、第m+1次モードの周波数以上では利得が1の周波数特性を有し、上記検出信号から第m+1次モードの周波数以上の高域信号成分を抽出する。
【0017】
上記高域通過フィルタは、予めその周波数特性をフーリエ逆変換して得られる時間関数と過去の所定時間t−T(tは現時刻を表す)より未来の所定時間t+Tまでの有限時間とを設定してなり、この時間関数と検出信号とを過去の所定時間t−2Tより現在時間まで畳み込むことにより、過去の所定時間t−Tにおける高域信号成分を出力する。
【0018】
上記高域信号成分の外挿波形は第m+1次モードの周波数を持つ正弦波である。
【0019】
【発明の実施の形態】
以下、本発明の一実施形態を添付図面に基づいて詳述する。
【0020】
図1に示されたモデルは、比較的柔軟な構造物を対象物としたものである。対象物1は細い柔軟な梁である。対象物1の一端は壁2に固定され、他端は自由となっている。この自由端を加振点3とする。この加振点3はこれから制振しようとする振動をこのモデルに与えるために設けた。対象物1の適宜な位置にアクチュエータが取り付けられており、その位置を制御点4という。また、振動を検出するセンサが取り付けられており、その位置をセンシング点5という。このモデルでは対象物1は壁2から加振点3までの150cmの長さがある。ここで用いられるセンサは速度センサであり、壁2から104cmの位置に制御点4及び速度センシング点5が設けられている。なお、対象物1は厚さ2mm、ヤング率2.0×1011N/m2 、ポアソン比0.3、各振動モードに減衰比0.01を設定してある。また、制御点4、加振点3には力発生用ボイスコイル、センサによる質量175g,100gがそれぞれ付加される。
【0021】
図2に示されるように、アクティブ制振装置は、対象物である柔軟構造体21に取り付けられたセンサ22と、センサ22の検出信号を所定周波数でサンプリングしてデジタルデータ列に変換するA/D変換器23と、そのデジタルデータ列を処理してアクチュエータ駆動指令値列を生成するパソコン等の計算機24と、そのアクチュエータ駆動指令値列を上記サンプリング周波数でアナログ波形信号に変換するD/A変換器25と、このアナログ波形信号で駆動制御され、柔軟構造体21を加振するアクチュエータ26とからなる。センサ22は柔軟構造体21の振動を検出するものであり、その検出量は変位、速度、加速度のいずれでもよい。これに対し、アクチュエータ26は力、圧力等を働かせるものである。アクチュエータ26は、柔軟構造体21外の図示されない固定点で支持されたものでもよいし、柔軟構造体21上でマスを移動させ反力を取り出すようなものでもよい。
【0022】
本発明に用いるために計算機24に組み込まれた高域通過フィルタの周波数特性を図3に示す。図示されるように、対象物の振動モードの第m次モードの周波数ωm より低周波では利得が0、かつ位相遅れが0であり、第m+1次モードの周波数ωm+1 より高周波では利得が1、かつ位相遅れが0である。周波数ωm からωm+1 までの領域における利得の変化は任意である。
【0023】
この特性を周波数伝達関数F(ω)で示すと、
F(ω)=0 0<ω≦ωm
F(ω)=1 ωm+1 ≦ω
であり、図4に示されるように、この高域通過フィルタ41において、センサの検出信号をy(t)、フィルタ出力をy1 (t)とし、積分範囲をt−∞≦τ≦t+∞とする代わりに、t−T≦τ≦t+Tの範囲のみで計算すると、
【0024】
【数1】
【0025】
となる。このフィルタ出力y1 (t)が高域信号成分である。
【0026】
ここで、h(t)は、周波数伝達関数F(ω)のフーリエ逆変換として定まる時間関数即ちインパルス応答関数である。即ち、
【0027】
【数2】
【0028】
である。この式(2)の計算は複雑かつ計算ステップ数が膨大であるが、フィルタ固有のものであるから、予め大型計算機等の上記計算機24とは別の計算機で計算しておき、その結果を計算機24に組み込み、計算機24では式(1)を計算するようにしてもよい。
【0029】
インパルス応答関数h(t)の一例を図5に示す。インパルス応答関数h(t)は、t=0で最大値を示し、他にもいくつかのピークがある。図には描ききれないが、インパルス応答関数h(t)は、過去無限大時間から未来無限大時間まで存在する。しかし、実用上は、このような無限大時間のデータは取り扱えないので、過去の所定時間−Tより未来の所定時間Tまでの有限時間を設定する。所定時間Tの決め方としては、例えばt=0における最大値に対しh(t)の絶対値が数%の値となるピークまでを含むものとする。この範囲外は無限大時間まで存在しても値が小さいので無視し得る。
【0030】
式(1)のようにインパルス応答関数h(t)と検出信号y(t)とを−T≦t≦Tの範囲で畳み込むことによりフィルタ出力y1 (t)が得られる。有限時間を設定してあるので、計算ステップ数が少なく、従って、パソコン等の計算機24によるリアルタイム処理が可能である。しかし、積分範囲に未来の時間が含まれており、未来の時間の入力信号は未知であるため、このままでは実用できない。
【0031】
そこで、式(1)の適用を時間Tだけ過去にシフトする。即ち、インパルス応答関数と検出信号とを過去の所定時間t−2Tより現在時間まで畳み込むことにより、過去の所定時間t−Tにおけるフィルタ出力を計算することにする。
【0032】
【数3】
【0033】
この式(3)によるフィルタ出力y1 (t−T)は、過去の所定時間t−Tにおける高域信号成分である。図6に示されるように、現在時間tまでの検出信号y(t)は既知であり、この検出信号y(t)を用いて、時間t−Tまでのフィルタ出力61を計算することは可能である。このようにして過去の所定時間t−Tにおける高域信号成分y1 (t−T)が求まる。
【0034】
しかし、求めたい値は現在時間のものである。そこで、フィルタ出力61の時間的変化を調べることにより、その概ねの周期、振幅が判る。これに従って、フィルタ出力62を現在時間まで外挿し、現在時間の値を得る。具体的には、y1 (τ)(τ≦t−T)を周波数ωm+1 を持つ正弦波又は余弦波等の調和関数で表す。この調和関数を時間Tだけ延長することにより、y1 (t)を容易に得ることができる。これは、カーブフィットと呼ばれる手法であり、特に柔軟な構造物にあっては実際の高域信号成分と殆ど同じものが得られる。
【0035】
この外挿による現在時間の高域信号成分を検出信号から差し引くことにより、低域信号成分のみを検出することができる。即ち、式(4)により低域信号成分ylow (t)が求まる。
【0036】
【数4】
【0037】
図2のアクティブ制振装置では、センサ22の検出信号はA/D変換器23により所定周波数でサンプリングしてデジタルデータ列yi に変換されているので、計算機24では、
【0038】
【数5】
【0039】
を計算することになる。式(5)のAjは、式(3)及びその外挿部分に関して予めオフラインで求めた定数である。従って、デジタルデータ列yi から簡単に低域信号成分のみからなるデジタルデータ列ylow,i を計算できる。
【0040】
計算機24では、低域信号成分のみからなるデジタルデータ列ylow,i に基づき、対象物である柔軟構造体21の振動を打ち消すように各サンプリング時間毎のアクチュエータの指令値が計算される。このようにして得られたアクチュエータ駆動指令値列がD/A変換器25によりアナログ波形信号に変換される。このアナログ波形信号でアクチュエータを駆動する。
【0041】
アクチュエータ駆動指令値の計算の基になった低域信号成分ylow (t)は、図3の特性を有する高域通過フィルタの出力を検出信号から差し引いたものであるから、図9に実線で示した特性を有する低域通過フィルタの出力に等しい。従って、計算機24で行われるアクチュエータ駆動指令値の計算には第m+1次モードの周波数以上の成分が関与しなくなり、アクチュエータはその能力で十分対応しきれる駆動がなされるようになり、スピルオーバが生じない。これにより、図8に示されるように、第m+1次モードの周波数ωm+1 において応答振幅が無制御のときと同じ程度まで抑えることができる。また、制御ゲインを上げてもスピルオーバが小さくなる。
【0042】
図1に示されたモデルに対し、第2次モードまでを制振帯域としたときの結果を説明する。図10及び図11は、横軸に加振点3に加えた加振周波数、縦軸に(センシング点5における速度振幅)/(加振点3に加えた加振力振幅)を取ってある。図10は従来技術、図11は本発明の結果を示す。
【0043】
図10に示されるように、無制御のとき、第1次モード,第2次モード,第3次モードの振動ピークが存在することが判る。第2次モードまでを対象に従来技術によるアクティブ制振を行うと、第3次モードでは振動が15dB(6倍)も強調されてしまう。
【0044】
図11に示されるように、第2次モードまでを対象に本発明によるアクティブ制振を行うと、第1次モード,第2次モードの振動が小さくなることは勿論、第3次モードは無制御のときと同程度にとどまっている。
【0045】
【発明の効果】
本発明は次の如き優れた効果を発揮する。
【0046】
(1)第m+1次モードのスピルオーバがなくなり、対象物が目論見通り制振できるようになる。
【0047】
(2)アクチュエータが高い周波数に追従しなくてもよいので、簡素で安価なものを使用することができる。
【図面の簡単な説明】
【図1】本発明の一実施形態を示す振動体のモデル図である。
【図2】本発明に係るアクティブ制振装置のブロック図である。
【図3】本発明に用いる高域通過フィルタの周波数特性図である。
【図4】図3の特性を持つ高域通過フィルタのブロック図である。
【図5】図3の特性を持つ高域通過フィルタのインパルス応答関数の波形図である。
【図6】本発明の高域通過フィルタの出力算出方法を示す時間波形図である。
【図7】従来のアクティブ制振方法による応答振幅特性図である。
【図8】本発明の目的とする応答振幅特性図である。
【図9】従来及び本発明の目的とする低域通過フィルタの周波数特性図である。
【図10】モデル実験の結果を示す従来のアクティブ制振の周波数特性図である。
【図11】モデル実験の結果を示す本発明のアクティブ制振の周波数特性図である。
【符号の説明】
1 対象物
2 壁
3 加振点
4 制御点
5 センシング点
41 高域通過フィルタ
Claims (4)
- 対象物の振動を検出し、その検出信号に応じてアクチュエータを駆動することにより、対象物の振動を打ち消すアクティブ制振の方法において、上記アクチュエータ駆動による制振が有効な振動周波数の範囲全域で上記検出信号を遮断し、それより高い周波数では通過させる高域通過フィルタを用い、この高域通過フィルタを現在時間までの検出信号に適用して所定時間前の高域信号成分を抽出し、この高域信号成分の波形を現在時間まで外挿し、この外挿による現在時間の高域信号成分を上記検出信号から差し引くことにより上記アクチュエータ駆動による制振が有効な低域信号成分を検出し、その低域信号成分に応じて上記アクチュエータを駆動することを特徴とするアクティブ制振方法。
- 上記高域通過フィルタは、上記アクチュエータ駆動による制振が有効な振動周波数の範囲に上記対象物の振動モードのうち第1次モードより第m次モードまでが含まれるとき、第m次モードの周波数以下では利得が0、第m+1次モードの周波数以上では利得が1の周波数特性を有し、上記検出信号から第m+1次モードの周波数以上の高域信号成分を抽出することを特徴とする請求項1記載のアクティブ制振方法。
- 上記高域通過フィルタは、予めその周波数特性をフーリエ逆変換して得られる時間関数と過去の所定時間t−T(tは現時刻を表す)より未来の所定時間t+Tまでの有限時間とを設定してなり、この時間関数と検出信号とを過去の所定時間t−2Tより現在時間まで畳み込むことにより、過去の所定時間t−Tにおける高域信号成分を出力することを特徴とする請求項1又は2記載のアクティブ制振方法。
- 上記高域信号成分の外挿波形は第m+1次モードの周波数を持つ正弦波であることを特徴とする請求項1〜3いずれか記載のアクティブ制振方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP28693996A JP3823397B2 (ja) | 1996-10-29 | 1996-10-29 | アクティブ制振方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP28693996A JP3823397B2 (ja) | 1996-10-29 | 1996-10-29 | アクティブ制振方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH10133746A JPH10133746A (ja) | 1998-05-22 |
JP3823397B2 true JP3823397B2 (ja) | 2006-09-20 |
Family
ID=17710917
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP28693996A Expired - Fee Related JP3823397B2 (ja) | 1996-10-29 | 1996-10-29 | アクティブ制振方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3823397B2 (ja) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4616520B2 (ja) * | 2001-06-29 | 2011-01-19 | 日立オートモティブシステムズ株式会社 | 制振装置 |
JP4228389B2 (ja) * | 2003-06-30 | 2009-02-25 | 株式会社安川電機 | 位置決め制御装置 |
JP4713355B2 (ja) * | 2006-02-02 | 2011-06-29 | 日立建機株式会社 | アクティブ振動制御装置 |
EP2706034B1 (de) | 2012-09-10 | 2015-11-04 | Integrated Dynamics Engineering GmbH | Aktiver Tilger für tieffrequent schwingende Strukturen |
-
1996
- 1996-10-29 JP JP28693996A patent/JP3823397B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JPH10133746A (ja) | 1998-05-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR101628119B1 (ko) | 소음제어시스템 및 그 방법 | |
JP5113145B2 (ja) | 並列適応フィルター構成を用いるアクティブノイズコントロールのためのシステム | |
JP3451891B2 (ja) | 能動型振動制御装置 | |
JP5691804B2 (ja) | マイクロホンアレイ装置及び音信号処理プログラム | |
US5319715A (en) | Noise sound controller | |
JP3486673B2 (ja) | 能動振動減衰電磁加振機及び能動振動減衰制御方法 | |
JP3823397B2 (ja) | アクティブ制振方法 | |
Halkyard et al. | Feedforward adaptive control of flexural vibration in a beam using wave amplitudes | |
JP3732227B2 (ja) | 繰り返し事象を制御する適応制御システム | |
JPS63208913A (ja) | 振動制御装置 | |
JPH0844377A (ja) | 周期性信号の適応制御方法 | |
JP2005257919A (ja) | 像振れ補正装置 | |
JPH06235439A (ja) | 振動制御方法および装置 | |
JP2013087857A (ja) | 制振システムおよび制振方法 | |
JPS6353617A (ja) | アクテイブ振動絶縁方法 | |
JP4760175B2 (ja) | 像振れ補正装置 | |
JPH08261277A (ja) | 車両の振動低減装置 | |
EP0598120B1 (en) | Noise control apparatus | |
JP4392700B2 (ja) | 板状体の振動抑制方法 | |
JPH0580777A (ja) | 車室内騒音の能動消音装置 | |
JP2841585B2 (ja) | 車室内騒音の低減装置 | |
JP3371477B2 (ja) | 車両用能動型振動制御装置及び車両用能動型騒音制御装置 | |
JP7288841B2 (ja) | センサシステム、および該センサシステムを備える除振装置 | |
JP3265793B2 (ja) | 能動型振動制御装置及び車両用能動型振動制御装置 | |
JP4713355B2 (ja) | アクティブ振動制御装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
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: 20060606 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20060619 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090707 Year of fee payment: 3 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
LAPS | Cancellation because of no payment of annual fees |