JP3724829B2 - Vibration control method for elastic structure - Google Patents

Vibration control method for elastic structure Download PDF

Info

Publication number
JP3724829B2
JP3724829B2 JP31356994A JP31356994A JP3724829B2 JP 3724829 B2 JP3724829 B2 JP 3724829B2 JP 31356994 A JP31356994 A JP 31356994A JP 31356994 A JP31356994 A JP 31356994A JP 3724829 B2 JP3724829 B2 JP 3724829B2
Authority
JP
Japan
Prior art keywords
vibration
experimental
actuator
mode
control
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 - Lifetime
Application number
JP31356994A
Other languages
Japanese (ja)
Other versions
JPH08171403A (en
Inventor
一登 背戸
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.)
Oiles Corp
Original Assignee
Oiles Corp
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 Oiles Corp filed Critical Oiles Corp
Priority to JP31356994A priority Critical patent/JP3724829B2/en
Publication of JPH08171403A publication Critical patent/JPH08171403A/en
Application granted granted Critical
Publication of JP3724829B2 publication Critical patent/JP3724829B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Landscapes

  • Building Environments (AREA)
  • Buildings Adapted To Withstand Abnormal External Influences (AREA)
  • Feedback Control In General (AREA)

Description

【0001】
【産業上の利用分野】
本発明は、弾性構造物のモデル化方法及び振動制御方法に係り、特に、超高層ビルなどの弾性構造物の制振設計に際して有効な弾性構造物のモデル化方法及び振動制御方法に関するものである。
【0002】
【従来の技術】
科学技術の発達に伴い、高層ビルや橋の主塔などの構造物の一層の高層化が可能となったが、風に対する揺れを制御する課題が生じてきた。そこで、質量の慣性力を利用して制振力を得るアクティブマスダンパー方式の弾性構造物の振動制御方法が注目され、実際に高層ビルへの適用例も多く見られるようになってきた。しかし、最近発表された1000m級の超高層ビル未来構想のような高層化が更に進めば、風による揺れは長周期になり、アクティブマスダンパー方式では、充分な制御力を得ることは困難になることが予想される。そこで、これに代わる振動制御方法として、構造物を並立させ、各構造物間に複数のアクチュエータを配置して、各構造物の相互作用により振動を同時に制御するビル連結方式という方法が提案されている(特願平5-141917号等)。
【0003】
詳述すると、この振動制御方法は、対向する弾性構造物間の相互作用を利用して、両者の間に配置されたアクチュエータによって構造物の振動をアクティブにかつ同時に制御する方法である。この方式をアクティブマスダンパー方式と比較すると、次のような特徴がある。
【0004】
(1)制振力を各構造物の反力で受けるので、静的、動的制振力が得られる。長周期の揺れに対する制振力の問題はない。
【0005】
(2)補助質量が不要であるので、制振装置が軽量化され、その設置空間が要らない。
【0006】
(3)アクチュエータによる結合部が両ビル間を結ぶ通路などに有効利用できる。
【0007】
図10は、この方式を適用した並立する2つの弾性構造物即ち超高層ビルを示す。各ビルa,b単体では強風に対して揺れて居住性が甚だしく損なわれるが、超高層ビルa,b間に幾つかの制御装置(アクチュエータ)cを配してビルa,b相互に揺れを止め合うならば、この問題は解消される。制御装置cの内部を通路に利用することも可能である。これを積極的に利用するならば、立体的な交通システムが形成され、便利で居心地の良い超高層ビル群が構成されよう。
【0008】
【発明が解決しようとする課題】
ところで、従来では、制御系が設計しやすい集中定数系に近い2つの弾性構造物を並立させ、制振対象を1自由度系としてビル連結方式の制振効果を調べ、さらにこれを2自由度系に拡張し、数値解析と実験によりその有効性を確認している。
【0009】
しかしながら、実際の構造物は分布定数系であり、分布定数系を振動制御する有効な制御系設計法がないために結果の一般性に問題が生ずる。
【0010】
そこで、上記課題を解決すべく本発明は創案されたものであり、その目的は、実際の弾性構造物を最適に集中定数系に低次元モデル化し、このモデルに基づく最適な振動制御を可能とする弾性構造物のモデル化方法及び振動制御方法を提供することにある。
【0011】
【課題を解決するための手段】
上記目的を達成するために本発明は、並立される弾性構造物をアクチュエータで互いに連結して、該アクチュエータを介して力を及ぼし合うことにより振動を制御するに際して、分布定数系特性を有する柔軟な板構造物を並立させ、それら各板構造物の振動モード形を解析し、低次元化のために、制御対象とする振動モードの次数を求めてそれ以上の高次を無視すると共に、無視した最低次の高次モードの振動の節を求め、その求めた節に実験アクチュエータを各々配置すると共に上記板構造物の間に挟んで、これにより上記弾性構造物を仮想した実験モデルを作成し、上記実験モデルの各板構造物を、上記最低次の高次モードの振動の節を質点として集中定数系にモデル化し、かつ上記実験アクチュエータをモデル化して、力学モデルを作成し、その力学モデルから算出された状態方程式を基に、上記実験モデルの振動を制御するための上記実験アクチュエータの制御量を求めて、これに基づいて上記弾性構造物の制振系を最適化するものである。
【0013】
【作用】
前者によれば、弾性構造物は分布定数系構造物として仮想され、分布定数系モデルとして扱われる。
【0014】
また、後者によれば、分布定数系モデルに基づいて制振系が最適化される。
【0015】
【実施例】
以下本発明の好適実施例を添付図面に基づいて詳述する。
【0016】
図1は、図10で示した並立する弾性構造物を分布定数系板構造物にモデル化した図である。図示するように、モデルM(実験モデル)は、並立する柔軟な2枚の板構造物1,2と、その間に挟まれた2つのアクチュエータ(実験アクチュエータ)3,4より成っている。板構造物1,2は矩形状とされ、それぞれ高さがL=1000mm、幅がW=70mmと等しくされると共に、厚さがt1 = 8mm、t2 =12mmと異なっている。また板構造物1,2は典型的な分布定数系特性を有する。実際の超高層ビルは分布定数系に近い弾性構造物であるために、これを振動モデルたる板構造物1,2として仮想することにより、実際の弾性構造物を最適にモデル化することができる。なおアクチュエータ3,4の詳細については後述する。
【0017】
ところで、このような分布定数系特性を有する板構造物を制振対象に用いた場合、課題となるのが分布定数系から集中定数系への低次元化の方法である。そこで以下においては、並立する2つの板構造物の1次と2次の振動モードを制御対象と決め、この並立する板構造物を2つの単体構造物と考え、それぞれについて2自由度系低次元化物理モデルを作成する。特にその際には、本発明者らが提案する「不可制御、不可観測性を活用した低次元化モデル作成法」を用いている(参考文献;背戸一登、振動のアクティブ制御、日本音響学会誌第47巻 9号,pp. 668-677 ,1991年)。
【0018】
制御系設計にはLQ制御理論を用いるのが一般的であるが、適切な低次元化モデル作成法が存在しないため、しばしばスピルオーバと呼ばれる不安定振動の発生が問題になっている。そこで、提案されたこの手法は、モード形の節が構造的に持つ不可制御、不可観測性の性質を活用して、低次元化のために無視した高次モードのスピルオーバ防止を考慮して、指定した点を質点とする集中定数系モデルを作成する手法である。実験モード解析と感度解析を用いて固有モード行列を修正し、質量やばね定数などの物理量を得ている。以下に、低次元化モデル作成法の手順を簡単に述べる。
【0019】
(1) 実験モード解析法やFEM などによって各板構造物の振動モード形を解析する。
【0020】
(2) 低次元化のために無視したい最低次の高次モードの振動の節を集中定数系モデルの質点に指定する。
【0021】
(3) その点で制御対象モードのモード成分を読取り、仮のモード行列を作る。
【0022】
(4) 感度解析法を応用したモード行列の修正を行い、それを用いて集中定数系の物理量を求める。
【0023】
図2には、実験モード解析法により求めた各板構造物1,2の1次から4次までの振動モード形を示す。ここでは2次モードまでが制振対象モードであるから、3次モードの2つの節を質点にとって、図に示すような2自由度系モデルを設定する。そして、そこにアクチュエータを取り付ければ3次モードは励起されず、3次モードの制御スピルオーバは防止される。また、4次モードの節にセンサを取り付ければ4次モードの観測スピルオーバが起こらない。
【0024】
ここでは、図2に示す2つの質量とばねをもつ物理モデルを仮定し、1次と2次のモーダルパラメータから2自由度系の物理モデルを作成する。ここで、モード質量を単位行列になるように正規化した固有モードを正規化固有モードΦと呼び、固有振動数を対角要素とする周波数行列をΩとおくと、質量行列M、剛性行列Kは以下のようになる。
【0025】
【数1】

Figure 0003724829
【0026】
ここで、3次モードの節である点20と点34の1次の固有モード成分をφ 11,φ21、2次の成分をφ12,φ22とすると、仮の正規モード行列は次のようになる。
【0027】
【数2】
Figure 0003724829
【0028】
よって、式(1) は次式で表すことができる。
【0029】
【数3】
Figure 0003724829
【0030】
これより質量集中の条件を満たすためには、次の拘束条件を満足しなければならない。
【0031】
【数4】
Figure 0003724829
【0032】
そこで、上式を物理モデルの変換誤差と考え、次の誤差関数ε1 を定義する。
【0033】
【数5】
Figure 0003724829
【0034】
この誤差関数ε1 を零にするような正規化固有モードΦの修正を行うために、Φに対する誤差関数の感度行列を次のようにおく。
【0035】
【数6】
Figure 0003724829
【0036】
ここで、φ11、φ21、φ12、φ22について修正量をそれぞれδφ11、δφ21、δφ12、δφ22とすると、誤差関数ε1 を零とする修正量は次式で求められる。
【0037】
【数7】
Figure 0003724829
【0038】
上式に最小ノルム解を使った一般化逆行列を適用すると、修正量は次式で求めることができる。
【0039】
【数8】
Figure 0003724829
【0040】
この修正量を用いて誤差関数を零に収束させれば、式(5) を満足する固有モードが得られる。
【0041】
こうして得た固有モードを式(1),(2) に代入して、質量行列M、剛性行列Kを求めるわけだが、このような値をもつ2自由度系物理モデルは、自由端に仮想のばねを取り付けることによって実現が可能となる。
【0042】
以下に得られた各物理定数を示す。
【0043】
11=2.0Kg ,M12=1.5Kg
21=1.1Kg ,M22=0.8Kg
10=-18532N/m ,K11=44904N/m,K12=40177N/m
20=-3060N/m,K21=7233N/m ,K22=6752N/m
さて、以上により作成した集中定数系モデルに基づき状態方程式を導入する。図3に2つの2自由度系に低次元化された制御対象と2つのハイブリッド型のアクチュエータからなる力学モデルを示す。図において、K10、K20は低次元化の際に考慮した仮想のばねである。また、用いたハイブリッドアクチュエータの減衰係数は次のようになっている。
【0044】
1 =2.2Ns/m ,C2 =1.9Ns/m
この力学モデルの運動方程式を求め、状態変数ベクトルを次のように定義すると、
【0045】
【数9】
Figure 0003724829
【0046】
状態方程式は次のようになる。
【0047】
【数10】
Figure 0003724829
【0048】
ただし、
【0049】
【数11】
Figure 0003724829
【0050】
c1,Kc2:各々のアクチュエータの制御力fc と制御量uの関係を示す力変換係数
次に、制御系の設計にはLQ制御理論を用いる。設計パラメータは、次に示す線形二次形式の評価関数Jに与える重み係数行列Q、Rである。ただし、Qは状態ベクトルXに、Rは制御量uに掛かる重み係数行列である。
【0051】
【数12】
Figure 0003724829
【0052】
LQ制御理論に基づけば、この評価関数Jを最小にする制御量uは次のように定式化されており、リカッチ方程式を解くことによりフィードバックゲインKを決定することができる。
【0053】
【数13】
Figure 0003724829
【0054】
ここに、Pは次のリカッチ方程式の解である。
【0055】
【数14】
Figure 0003724829
【0056】
図4は、重み係数行列Qの違いによる、各構造物の周波数応答を比較したシミュレーション結果を示す。ここで、細い実線は非制御時、中太実線は変位項に次のような重みを掛けた場合、
【0057】
【数15】
Figure 0003724829
【0058】
また、太い実線は速度項に次のような重みを掛けた場合を示す。
【0059】
【数16】
Figure 0003724829
【0060】
同図(a)は、各構造物の点20を加振したことによる、各構造物1の点35、(b)は構造物2の点35で、各々観測された周波数応答である。見られるように、構造物1の1次と2次の共振ピークは、5.8Hz と37Hzに現れ、構造物2のそれは、3.2Hz と21Hzに現れている。これらの共振ピークは、変位に重みを掛けた場合には各構造物共に1次のピークは良く抑制されるが、2次はあまり抑制されない。一方、速度に重みを掛けることによって、各構造物共に1次、2次の共振ピークが良く抑制され、良好な振動制御結果が得られることが分る。
【0061】
次に、各質点のインパルス応答をシミュレーションによって比較したものを図5に示す。ここに、細線は非制御時の応答、太い実線は上記の重みを掛けた場合の制御時の応答を示す。図5(a)に見られるように、変位に重みを掛けた場合には、明かに2次モードの振動が良く抑制されないままに持続している。一方、速度項に重みを掛けた場合のインパルス応答を図5(b)に示す。制御を掛けることにより全ての質点において1次と2次モードの振動が速やかに収束していることが確認される。
【0062】
また、非制御時にも振動が収束していることから、ここで提案するハイブリッド型のアクチュエータは、アクティブ制御が使用できない場合にも、減衰力によりある程度の信頼性を確保できることが確認できる。
【0063】
次に、上述のシミュレーション結果を確認するために行った制御実験の結果を示す。
【0064】
先ず図6に実験装置の構成を示す。実験装置は4つの変位センサ、パーソナルコンピュータ、A/D 変換器、D/A 変換器、二つの電磁アクチュエータ3,4、および駆動回路から構成されている。アクチュエータ3,4は、図7に示すコイルと永久磁石からなる非接触の電磁力タイプであり、コイルを巻いたアルミケースと永久磁石との間の渦電流損失によりパッシブ制御力も得ることのできるハイブリッド型のものである。
【0065】
各板構造物1,2の絶対変位は、変位センサにより検出され、A/D 変換器を介してコンピュータに入力される。そしてコンピュータ内で、入力された変位信号の差分をとることで速度信号を求め、8個の状態量を得ている。また、シミュレーションにより得られたフィードバックゲインを実験でも直接使用して制御系を決定し、制御実験を行っている。
【0066】
前述の通り、アクチュエータを3次モードの節に設置することで3次モードによる制御スピルオーバを、センサを4次モードの節に設置することで4次モードによる観測スピルオーバを防いでいる。ちなみに、アクチュエータは図2に示す20点と34点付近に、センサは25点と36点付近にそれぞれ設置している。
【0067】
なお、ここで扱ったような、一様断面、同一高さの構造物では、振動モード形はほぼ類似しており、振動の節も同じ様な位置に現れるので、アクチュエータを同一高さにとる手法にはかなり一般性があると考えている。
【0068】
図8には、図5(b)のシミュレーションに対応する、速度に重みを掛けた場合のインパルス応答の実験結果を示す。細線が非制御時、太線が制御時の応答であり、これもシミュレーション結果と実験結果は極めてよく一致している。このことは、ここで扱うような分布定数系特性を有する弾性構造物の振動制御における、本モデル化方法の有効性が改めて確認されたことになる。
【0069】
また、いずれも制御結果においても、全ての質点においてスピルオーバを起こすことなく、1次と2次の振動が速やかに収束していることが確認される。これより、3次と4次のスピルオーバの防止においても、本方法が有効であることが確認できる。また、非制御時の応答よりアクチュエータの信頼性が確認できる。
【0070】
最後に、各構造物の上部2つのセンサ測定された周波数応答の実験結果を図9に示す。これは速度に重みを掛けた図4のシミュレーション結果に対応するものであるが、これらの実験結果もシミュレーション結果と良く一致しており、各共振ピークが制御によって十分に抑制されている。このことは、各振動モードがアクティブにダンピングを得たことを示すものであり、超々高層ビルのダンピング不足による風励振の問題は、このような振動制御方法で解消できることを示している。
【0071】
このように、アクチュエータを含む制振系が最適化され、超々高層ビルなどで課題となる風による揺れ制御に、かかる弾性構造物のモデル化方法及び振動制御方法が有効であることが確認された。
【0072】
尚、かかる方法は、高層ビルに限らず宇宙構造物、弾性車両、建築機械等、並立する構造物が存在する多くの対象に対して適用可能である。
【0073】
【発明の効果】
本発明は次の如き優れた効果を発揮する。
【0074】
(1)実際の弾性構造物を最適にモデル化することができる。
【0075】
(2)最適なモデルに基づく最適な振動制御が可能になる。
【図面の簡単な説明】
【図1】図10の並立する弾性構造物をモデル化した斜視図である。
【図2】各板構造物の1次から4次までの振動モード形と集中定数系モデルを示し、(a)は板構造物1の場合、(b)は板構造物2の場合である。
【図3】2つのアクチュエータを有する並立構造物の力学モデルを示す。
【図4】周波数応答のシミュレーション結果を示し、(a)は板構造物1の場合、(b)は板構造物2の場合である。
【図5】インパルス応答のシミュレーション結果を示し、(a)は変位に重みを掛けた場合、(b)は速度に重みを掛けた場合である。
【図6】実験装置の構成図である。
【図7】アクチュエータを示す分解斜視図である。
【図8】速度に重みを掛けたときのインパルス応答の実験結果を示し、(a)は板構造物1の場合、(b)は板構造物2の場合である。
【図9】周波数応答の実験結果を示し、(a)は板構造物1の場合、(b)は板構造物2の場合である。
【図10】並立する弾性構造物即ち超高層ビルを示す斜視図である。
【符号の説明】
1,2 板構造物
a,b 弾性構造物
c アクチュエータ[0001]
[Industrial application fields]
The present invention relates to an elastic structure modeling method and a vibration control method, and more particularly to an elastic structure modeling method and a vibration control method that are effective in damping control of an elastic structure such as a high-rise building. .
[0002]
[Prior art]
With the development of science and technology, it has become possible to further increase the height of structures such as high-rise buildings and main towers of bridges, but there has been a problem of controlling fluctuations in the wind. Accordingly, attention has been paid to a vibration control method of an elastic structure of an active mass damper system that obtains a damping force by using an inertial force of mass, and many examples of application to a high-rise building have come to be seen in practice. However, if the high-rise building such as the recently announced 1000m-class skyscraper future concept is further advanced, the vibration caused by the wind will become a long cycle, and it will be difficult to obtain sufficient control power with the active mass damper system. It is expected that. Therefore, as an alternative vibration control method, a method called a building connection method has been proposed in which structures are arranged side by side, a plurality of actuators are arranged between the structures, and vibration is controlled simultaneously by the interaction of the structures. (Japanese Patent Application No.5-141917 etc.).
[0003]
More specifically, this vibration control method is a method of actively and simultaneously controlling the vibration of a structure by an actuator disposed between the two using an interaction between opposing elastic structures. When this method is compared with the active mass damper method, it has the following characteristics.
[0004]
(1) Since the damping force is received by the reaction force of each structure, static and dynamic damping force can be obtained. There is no problem of damping force against long-period shaking.
[0005]
(2) Since no auxiliary mass is required, the vibration damping device is reduced in weight, and its installation space is not required.
[0006]
(3) The connecting portion by the actuator can be used effectively for a passage connecting the two buildings.
[0007]
FIG. 10 shows two side-by-side elastic structures or skyscrapers to which this scheme is applied. Each building a and b alone is shaken by strong winds and the habitability is severely impaired. However, several control devices (actuators) c are arranged between the high-rise buildings a and b, and the buildings a and b shake each other. If we stop, this problem is solved. It is also possible to use the inside of the control device c as a passage. If this is used positively, a three-dimensional transportation system will be formed, and a convenient and comfortable skyscraper group will be constructed.
[0008]
[Problems to be solved by the invention]
By the way, in the past, two elastic structures close to a lumped parameter system, which is easy to design the control system, are arranged side by side, and the vibration control effect of the building connection system is investigated with the object of vibration control as a one-degree-of-freedom system. The system has been expanded to confirm its effectiveness through numerical analysis and experiments.
[0009]
However, the actual structure is a distributed parameter system, and there is no effective control system design method for controlling the vibration of the distributed parameter system.
[0010]
Accordingly, the present invention has been created to solve the above-mentioned problems, and the purpose of the present invention is to make an actual elastic structure into a low-dimensional model optimally in a lumped parameter system, and to enable optimal vibration control based on this model. Another object is to provide a modeling method and a vibration control method for an elastic structure.
[0011]
[Means for Solving the Problems]
In order to achieve the above object, the present invention provides a flexible constant distribution characteristic when controlling vibrations by connecting parallel elastic structures with an actuator and applying force via the actuator. The plate structures are arranged side by side, the vibration mode shape of each plate structure is analyzed, and the order of the vibration mode to be controlled is obtained and the higher order is ignored and ignored for the purpose of reducing the order. Determine the vibration node of the lowest order higher order mode, place the experimental actuators in the obtained nodes and sandwich them between the plate structures, thereby creating an experimental model that virtualizes the elastic structure, Each plate structure of the above experimental model is modeled as a lumped parameter system with the vibration node of the lowest order higher order mode as the mass point, and the experimental actuator is modeled to create a dynamic model. Then, based on the equation of state calculated from the dynamic model, obtain the control amount of the experimental actuator to control the vibration of the experimental model, and optimize the damping system of the elastic structure based on this To do .
[0013]
[Action]
According to the former, the elastic structure is hypothesized as a distributed constant system structure and is treated as a distributed constant system model.
[0014]
Further, according to the latter, the vibration damping system is optimized based on the distributed parameter system model.
[0015]
【Example】
Hereinafter, preferred embodiments of the present invention will be described in detail with reference to the accompanying drawings.
[0016]
FIG. 1 is a diagram in which the side-by-side elastic structure shown in FIG. 10 is modeled as a distributed constant plate structure. As shown in the figure, the model M (experimental model) includes two flexible plate structures 1 and 2 arranged side by side, and two actuators (experimental actuators) 3 and 4 sandwiched therebetween. The plate structures 1 and 2 have a rectangular shape, each having a height equal to L = 1000 mm, a width equal to W = 70 mm, and a thickness different from t 1 = 8 mm and t 2 = 12 mm. The plate structures 1 and 2 have typical distributed constant system characteristics. Since an actual skyscraper is an elastic structure close to a distributed constant system, an actual elastic structure can be optimally modeled by virtualizing it as plate structures 1 and 2 that are vibration models. . The details of the actuators 3 and 4 will be described later.
[0017]
By the way, when a plate structure having such a distributed constant system characteristic is used as a vibration control target, a problem is how to reduce the order from a distributed constant system to a lumped constant system. Therefore, in the following, the primary and secondary vibration modes of two parallel plate structures are determined as control objects, and the parallel plate structures are considered as two single structures, and each has a two-degree-of-freedom low-dimensional system. Create a physical model. In particular, in this case, the “method of creating a low-dimensional model utilizing uncontrollability and unobservability” proposed by the present inventors is used (reference: Kazuto Sado, active control of vibration, Acoustical Society of Japan 47, 9 (pp. 668-677, 1991).
[0018]
In general, LQ control theory is used for control system design. However, since there is no appropriate low-dimensional model creation method, generation of unstable vibration called spillover is a problem. Therefore, this proposed method takes advantage of the uncontrollable and unobservable nature of the mode-shaped nodes structurally, taking into account the prevention of spillover of higher-order modes ignored for lower dimensions This is a technique for creating a lumped parameter model with a specified point as a mass point. The experimental mode analysis and sensitivity analysis are used to correct the eigenmode matrix to obtain physical quantities such as mass and spring constant. In the following, the procedure for creating a low-dimensional model is briefly described.
[0019]
(1) Analyze the vibration mode shape of each plate structure by the experimental mode analysis method or FEM.
[0020]
(2) Specify the vibration node of the lowest order higher order mode that should be ignored for reduction of dimensions as the mass point of the lumped parameter model.
[0021]
(3) At that point, the mode component of the controlled object mode is read and a temporary mode matrix is created.
[0022]
(4) The mode matrix is modified by applying the sensitivity analysis method, and the physical quantity of the lumped parameter system is obtained using it.
[0023]
FIG. 2 shows the first to fourth vibration mode shapes of the plate structures 1 and 2 obtained by the experimental mode analysis method. Here, up to the second-order mode is the vibration suppression target mode, and the two-degree-of-freedom system model as shown in the figure is set with the two nodes of the third-order mode as the mass points. If an actuator is attached to the actuator, the third-order mode is not excited and control spillover of the third-order mode is prevented. In addition, if a sensor is attached to the fourth-order mode node, the fourth-order mode observation spillover will not occur.
[0024]
Here, a physical model having two masses and a spring shown in FIG. 2 is assumed, and a two-degree-of-freedom physical model is created from first-order and second-order modal parameters. Here, the eigenmode that is normalized so that the mode mass becomes a unit matrix is called a normalized eigenmode Φ, and when the frequency matrix having the natural frequency as a diagonal element is denoted by Ω, the mass matrix M and the stiffness matrix K Is as follows.
[0025]
[Expression 1]
Figure 0003724829
[0026]
Here, assuming that the primary eigenmode components of points 20 and 34 which are nodes of the third-order mode are φ 11 and φ 21 , and the second-order components are φ 12 and φ 22 , the temporary normal mode matrix is It becomes like this.
[0027]
[Expression 2]
Figure 0003724829
[0028]
Therefore, Formula (1) can be expressed by the following formula.
[0029]
[Equation 3]
Figure 0003724829
[0030]
In order to satisfy the mass concentration condition, the following constraint condition must be satisfied.
[0031]
[Expression 4]
Figure 0003724829
[0032]
Therefore, considering the above equation as a conversion error of the physical model, the following error function ε 1 is defined.
[0033]
[Equation 5]
Figure 0003724829
[0034]
In order to correct the normalized eigenmode Φ so that the error function ε 1 becomes zero, the sensitivity matrix of the error function with respect to Φ is set as follows.
[0035]
[Formula 6]
Figure 0003724829
[0036]
Here, assuming that the correction amounts for φ 11 , φ 21 , φ 12 , and φ 22 are δφ 11 , δφ 21 , δφ 12 , and δφ 22 , the correction amounts that make the error function ε 1 zero are obtained by the following equations.
[0037]
[Expression 7]
Figure 0003724829
[0038]
Applying the generalized inverse matrix using the minimum norm solution to the above equation, the correction amount can be obtained by the following equation.
[0039]
[Equation 8]
Figure 0003724829
[0040]
If the error function is converged to zero using this correction amount, an eigenmode satisfying Equation (5) is obtained.
[0041]
Substituting the eigenmodes obtained in this way into equations (1) and (2) to find the mass matrix M and stiffness matrix K. The two-degree-of-freedom physical model with these values Realization is possible by attaching a spring.
[0042]
Each physical constant obtained is shown below.
[0043]
M 11 = 2.0Kg, M 12 = 1.5Kg
M 21 = 1.1Kg, M 22 = 0.8Kg
K 10 = -18532N / m, K 11 = 44904N / m, K 12 = 40177N / m
K 20 = -3060N / m, K 21 = 7233N / m, K 22 = 6752N / m
Now, the equation of state is introduced based on the lumped parameter system model created as described above. FIG. 3 shows a dynamic model composed of a controlled object reduced to two two-degree-of-freedom systems and two hybrid actuators. In the figure, K 10 and K 20 are virtual springs that are taken into consideration when the dimensions are reduced. The damping coefficient of the hybrid actuator used is as follows.
[0044]
C 1 = 2.2Ns / m, C 2 = 1.9Ns / m
Obtaining the equation of motion of this dynamic model and defining the state variable vector as
[0045]
[Equation 9]
Figure 0003724829
[0046]
The equation of state is as follows.
[0047]
[Expression 10]
Figure 0003724829
[0048]
However,
[0049]
[Expression 11]
Figure 0003724829
[0050]
K c1 , K c2 : Force conversion coefficients indicating the relationship between the control force f c of each actuator and the control amount u Next, LQ control theory is used for designing the control system. The design parameters are weight coefficient matrices Q and R given to the evaluation function J in the linear quadratic form shown below. Here, Q is a state vector X, and R is a weighting coefficient matrix applied to the control amount u.
[0051]
[Expression 12]
Figure 0003724829
[0052]
Based on the LQ control theory, the control amount u that minimizes the evaluation function J is formulated as follows, and the feedback gain K can be determined by solving the Riccati equation.
[0053]
[Formula 13]
Figure 0003724829
[0054]
Here, P is a solution of the following Riccati equation.
[0055]
[Expression 14]
Figure 0003724829
[0056]
FIG. 4 shows a simulation result in which the frequency response of each structure is compared due to the difference in the weighting coefficient matrix Q. Here, when the thin solid line is not controlled, the thick solid line is multiplied by the following weight on the displacement term:
[0057]
[Expression 15]
Figure 0003724829
[0058]
A thick solid line indicates a case where the following weight is applied to the speed term.
[0059]
[Expression 16]
Figure 0003724829
[0060]
FIG. 6A shows the frequency response observed at the point 35 of each structure 1 and FIG. 5B shows the frequency response observed at the point 35 of the structure 2 due to the vibration of the point 20 of each structure. As can be seen, the primary and secondary resonance peaks of structure 1 appear at 5.8 Hz and 37 Hz, and that of structure 2 appears at 3.2 Hz and 21 Hz. As for these resonance peaks, when the weight is applied to the displacement, the primary peak is well suppressed in each structure, but the secondary is not much suppressed. On the other hand, it can be seen that by applying a weight to the speed, the primary and secondary resonance peaks of each structure are well suppressed and a good vibration control result can be obtained.
[0061]
Next, FIG. 5 shows a comparison of impulse responses of the respective mass points by simulation. Here, the thin line shows the response at the time of non-control, and the thick solid line shows the response at the time of control when the above weight is applied. As seen in FIG. 5 (a), when the displacement is weighted, the vibration in the secondary mode is clearly not suppressed well. On the other hand, FIG. 5B shows the impulse response when the speed term is weighted. By applying the control, it is confirmed that the vibrations of the primary and secondary modes are quickly converged at all mass points.
[0062]
In addition, since the vibration converges even during non-control, it can be confirmed that the hybrid actuator proposed here can secure a certain degree of reliability by the damping force even when active control cannot be used.
[0063]
Next, the result of the control experiment performed in order to confirm the above-mentioned simulation result is shown.
[0064]
First, FIG. 6 shows the configuration of the experimental apparatus. The experimental device consists of four displacement sensors, a personal computer, an A / D converter, a D / A converter, two electromagnetic actuators 3 and 4, and a drive circuit. The actuators 3 and 4 are non-contact electromagnetic force types composed of a coil and a permanent magnet shown in FIG. 7, and a hybrid capable of obtaining a passive control force by eddy current loss between the aluminum case wound with the coil and the permanent magnet. Of the type.
[0065]
The absolute displacement of each plate structure 1 and 2 is detected by a displacement sensor and input to a computer via an A / D converter. In the computer, the speed signal is obtained by taking the difference between the input displacement signals, and eight state quantities are obtained. In addition, the feedback gain obtained by the simulation is directly used in the experiment to determine the control system and conduct the control experiment.
[0066]
As described above, the control spillover by the third-order mode is prevented by installing the actuator in the third-order mode section, and the observation spillover by the fourth-order mode is prevented by installing the sensor in the fourth-order mode section. Incidentally, the actuators are installed near the 20 and 34 points shown in FIG. 2, and the sensors are installed near the 25 and 36 points, respectively.
[0067]
In the structure with the same cross section and the same height as described here, the vibration mode shape is almost similar, and the vibration nodes appear at the same position. I think the method is quite general.
[0068]
FIG. 8 shows the experimental results of the impulse response when the speed is weighted, corresponding to the simulation of FIG. The thin line is the non-control response and the thick line is the control response. The simulation results and the experimental results also agree very well. This confirms the effectiveness of the present modeling method in vibration control of an elastic structure having distributed constant system characteristics as dealt with here.
[0069]
Also, in both control results, it is confirmed that the primary and secondary vibrations converge quickly without causing spillover at all mass points. From this, it can be confirmed that the present method is also effective in preventing third-order and fourth-order spillover. In addition, the reliability of the actuator can be confirmed from the response during non-control.
[0070]
Finally, the experimental results of the measured frequency response of the top two sensors of each structure are shown in FIG. This corresponds to the simulation result of FIG. 4 in which the speed is weighted, but these experimental results also agree well with the simulation results, and each resonance peak is sufficiently suppressed by the control. This indicates that each vibration mode has actively obtained damping, and that the problem of wind excitation due to insufficient damping of a super high-rise building can be solved by such a vibration control method.
[0071]
In this way, the damping system including the actuator has been optimized, and it has been confirmed that the modeling method and the vibration control method of such an elastic structure are effective for the vibration control due to wind, which is a problem in ultra-high-rise buildings and the like. .
[0072]
In addition, this method is applicable not only to a high-rise building but to many objects in which a parallel structure exists such as a space structure, an elastic vehicle, and a construction machine.
[0073]
【The invention's effect】
The present invention exhibits the following excellent effects.
[0074]
(1) An actual elastic structure can be modeled optimally.
[0075]
(2) Optimal vibration control based on the optimal model becomes possible.
[Brief description of the drawings]
1 is a perspective view modeling the side-by-side elastic structures of FIG. 10;
FIGS. 2A and 2B show vibration mode shapes from the first to fourth orders and a lumped parameter system model of each plate structure, where FIG. 2A shows the case of the plate structure 1 and FIG. 2B shows the case of the plate structure 2; .
FIG. 3 shows a dynamic model of a side-by-side structure with two actuators.
4A and 4B show simulation results of frequency response, where FIG. 4A shows the case of the plate structure 1 and FIG. 4B shows the case of the plate structure 2;
FIGS. 5A and 5B show simulation results of impulse responses, in which FIG. 5A shows a case where weight is applied to displacement, and FIG. 5B shows a case where weight is applied to velocity.
FIG. 6 is a configuration diagram of an experimental apparatus.
FIG. 7 is an exploded perspective view showing an actuator.
FIGS. 8A and 8B show experimental results of impulse responses when weight is applied to the speed, where FIG. 8A shows the case of the plate structure 1 and FIG. 8B shows the case of the plate structure 2;
9A and 9B show experimental results of frequency response, where FIG. 9A shows the case of the plate structure 1, and FIG. 9B shows the case of the plate structure 2. FIG.
FIG. 10 is a perspective view showing a side-by-side elastic structure or skyscraper.
[Explanation of symbols]
1, 2 Plate structure a, b Elastic structure c Actuator

Claims (1)

並立される弾性構造物をアクチュエータで互いに連結して、該アクチュエータを介して力を及ぼし合うことにより振動を制御するに際して、
分布定数系特性を有する柔軟な板構造物を並立させ、それら各板構造物の振動モード形を解析し、低次元化のために、制御対象とする振動モードの次数を求めてそれ以上の高次を無視すると共に、無視した最低次の高次モードの振動の節を求め、その求めた節に実験アクチュエータを各々配置すると共に上記板構造物の間に挟んで、これにより上記弾性構造物を仮想した実験モデルを作成し、
上記実験モデルの各板構造物を、上記最低次の高次モードの振動の節を質点として集中定数系にモデル化し、かつ上記実験アクチュエータをモデル化して、力学モデルを作成し、
その力学モデルから算出された状態方程式を基に、上記実験モデルの振動を制御するための上記実験アクチュエータの制御量を求めて、これに基づいて上記弾性構造物の制振系を最適化することを特徴とする弾性構造物の振動制御方法。
When controlling vibrations by connecting parallel elastic structures to each other with an actuator and applying a force through the actuator,
Flexible plate structures with distributed constant system characteristics are arranged side by side, the vibration mode shapes of each plate structure are analyzed, and the order of the vibration mode to be controlled is obtained for further reduction in order to reduce the dimension. Ignore the next and determine the lowest order higher order mode vibration node, and place the experimental actuators in the determined nodes and sandwich them between the plate structures, so that the elastic structure Create a virtual experimental model,
Each plate structure of the experimental model is modeled as a lumped parameter system with the vibration node of the lowest order higher order mode as a mass point, and the experimental actuator is modeled to create a dynamic model,
Based on the equation of state calculated from the dynamic model, obtain the control amount of the experimental actuator for controlling the vibration of the experimental model, and optimize the damping system of the elastic structure based on this A vibration control method for an elastic structure characterized by the above.
JP31356994A 1994-12-16 1994-12-16 Vibration control method for elastic structure Expired - Lifetime JP3724829B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP31356994A JP3724829B2 (en) 1994-12-16 1994-12-16 Vibration control method for elastic structure

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP31356994A JP3724829B2 (en) 1994-12-16 1994-12-16 Vibration control method for elastic structure

Publications (2)

Publication Number Publication Date
JPH08171403A JPH08171403A (en) 1996-07-02
JP3724829B2 true JP3724829B2 (en) 2005-12-07

Family

ID=18042895

Family Applications (1)

Application Number Title Priority Date Filing Date
JP31356994A Expired - Lifetime JP3724829B2 (en) 1994-12-16 1994-12-16 Vibration control method for elastic structure

Country Status (1)

Country Link
JP (1) JP3724829B2 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4587011B2 (en) * 2000-11-24 2010-11-24 清水建設株式会社 Active seismic isolation control method and control device
CN102411380B (en) * 2011-11-30 2013-09-25 华南理工大学 Multi-channel control device and method of bending and torsion vibration of flexible hinged plate
CN106017839B (en) * 2016-07-18 2019-05-14 华南理工大学 Based on flexible, hinged plate benging and twisting vibration detection control apparatus and method

Also Published As

Publication number Publication date
JPH08171403A (en) 1996-07-02

Similar Documents

Publication Publication Date Title
Inman Active modal control for smart structures
Kapuria et al. Active vibration control of piezoelectric laminated beams with electroded actuators and sensors using an efficient finite element involving an electric node
Bandyopadhyay et al. Modeling, control and implementation of smart structures: a FEM-state space approach
Balamurugan et al. Finite element formulation and active vibration control study on beams using smart constrained layer damping (SCLD) treatment
Aldraihem et al. Optimal size and location of piezoelectric actuator/sensors: practical considerations
Holterman Vibration control of high-precision machines with active structural elements
Liu et al. Piezoelectric actuator placement optimization and active vibration control of a membrane structure
US5265387A (en) Vibration suppressing structure
Ramegowda et al. A novel coupling algorithm for the electric field–structure interaction using a transformation method between solid and shell elements in a thin piezoelectric bimorph plate analysis
Henrioulle et al. Experimental validation of a collocated PVDF volume velocity sensor/actuator pair
Kermani et al. Parameter selection and control design for vibration suppression using piezoelectric transducers
Chen et al. Structural analysis and optimal design of a dynamic absorbing beam
Ostrowski et al. Semi‐active modal control of structures with lockable joints: General methodology and applications
JP3724829B2 (en) Vibration control method for elastic structure
Baz et al. Active modal control of vortex-induced vibrations of a flexible cylinder
Piefort et al. Finite element modeling of piezoelectric structures
Strassberger et al. Active noise reduction by structural control using piezo-electric actuators
Kusculuoglu et al. Finite element formulation for composite plates with piezoceramic layers for optimal vibration control applications
Li et al. Suppression of thermal postbuckling and nonlinear panel flutter motions using piezoelectric actuators
Gardonio et al. Comparison of smart panels for tonal and broadband vibration and sound transmission active control
Barhoush et al. Vortex shedding analysis by finite elements
JPH03250165A (en) Hybrid dynamic vibration reducer
Adhikari et al. Modal space sliding‐mode control of structures
Tairidis et al. Classical and soft robust active control of smart beams
JP2732681B2 (en) Optimal active dynamic vibration absorber control system for multi-degree-of-freedom buildings

Legal Events

Date Code Title Description
A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20050627

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20050920

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090930

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090930

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100930

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100930

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110930

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110930

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120930

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120930

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130930

Year of fee payment: 8

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

EXPY Cancellation because of completion of term