水平板を用いた消波機構における指向性 アクチュエータの境界要素法による性能解析

Slides:



Advertisements
Similar presentations
第 7 週目: 周波数伝達関数とボード線図 周波数伝達関数 ボード線図 TUT, System & Control laboratory 1/16.
Advertisements

定在波型熱音響エンジンにおける 臨界温度比推定のための適応制御系の 安定性に関する実験と理論の比較 長岡技術科学大学
24 両端単純支持梁に対する外乱抑制制御系の製作
プロセス制御工学 3.伝達関数と過渡応答 京都大学  加納 学.
スペクトル法による数値計算の原理 -一次元線形・非線形移流問題の場合-
フィードバック制御に基づく 定在波型熱音響エンジンにおける 自励発振条件の特徴付け
超磁歪アクチュエータを用いた キャビテーション発生機構における 機械的特性の解析
プロセス制御工学 6.PID制御 京都大学  加納 学.
底面形状を考慮したH∞制御に基づく水槽内の定在波抑制制御
ー 第1日目 ー 確率過程について 抵抗の熱雑音の測定実験
ー 第3日目 ー ねじれ型振動子のブラウン運動の測定
不安定な補償器を用いた 低剛性・高慣性比の 二慣性ねじり振動系における 外乱抑制制御性能の改善
電気回路学Ⅱ エネルギーインテリジェンスコース 5セメ 山田 博仁.
2.伝送線路の基礎 2.1 分布定数線路 2.1.1 伝送線路と分布定数線路 集中定数回路:fが低い場合に適用
機械創造工学課程 08104288 鈴木翔 担当教員 小林泰秀 准教授
水平板を用いた消波機構における指向性 アクチュエータの境界要素法による性能解析
25 ロバスト制御に基づく柔軟ベルト駆動二慣性系の外乱抑制制御 機械創造工学課程 西村光博 担当教員 小林泰秀 准教授
定在波型熱音響エンジンの共鳴現象に 対するフィードバック制御の効果 長岡技術科学大学 ☆角島 悠太 小林 泰秀, 山田 昇.
学籍番号:   氏名:新保尚敬  指導教員:小林泰秀 准教授
28 PICマイコンを用いた能動騒音制御系の制御性能
閉ループ系を安定限界に保持する 適応制御に基づく定在波型 熱音響エンジンの定常発振制御
プロセス制御工学 7.多変数プロセスの制御 京都大学  加納 学.
31 ループ管熱音響システムにおける管内圧力の可視化 長岡技術科学大学 機械創造工学課程 梅本康平 担当教員 小林泰秀 准教授
従動側角速度フィードバック による不安定化に基づく 二慣性系の外乱抑制性能の改善
中性子干渉実験 2008/3/10 A4SB2068 鈴木 善明.
機械創造工学課程 西久保智昭 担当教員 小林泰秀 准教授
制御系における指向性アクチュエータの効果
ベクトル線図 周波数応答 G(jw) (– < w < ) を複素平面内に描いたものが、ベクトル線図である。
川崎浩司:沿岸域工学,コロナ社 第2章(pp.12-22)
1-P-6 パラボラ反射板を用いたアクティブマイクロフォンによる方向推定
安東 正樹池本尚史,小林洸,坪野公夫 (東京大学 理学系研究科)
逐次伝達法による 散乱波の解析 G05MM050 本多哲也.
基本システムのボード線図 ボード線図による基本システムの同定
高次システムのボード線図 周波数応答によるシステムの同定
高分解能ビーム軌道傾きモニターの設計開発
ー 第3日目 ー ねじれ型振動子のブラウン運動の測定
26ロバスト制御に基づく片持ち梁の外乱抑制制御系の設計
電気回路学Ⅱ コミュニケーションネットワークコース 5セメ 山田 博仁.
C4 能動騒音制御を用いたループ管熱音響冷却機の製作
22 物理パラメータに陽に依存する補償器を用いた低剛性二慣性系の速度制御実験 高山誠 指導教員 小林泰秀
7.一次元ダクトの消音制御系における低コスト化
30 両端単純支持梁に対する外乱抑制制御系の製作 機械創造工学課程 11307489 古澤大輔 担当教員 小林泰秀 准教授
電磁気学Ⅱ Electromagnetics Ⅱ 8/11講義分 点電荷による電磁波の放射 山田 博仁.
両端単純支持梁の フィードフォワード外乱抑制制御系における 指向性アクチュエータの効果
9. ナイキスト線図と安定余裕 教科書 7.2, 7.3.
21 柔軟片持ち梁の振動制御における移動可能なアクチュエータの製作
1-1-6 ロバスト能動騒音制御に基づく ループ管熱音響システムにおける 定在波抑制制御の効果
熱音響コアが多段接続された 電力フィードバック進行波型熱音響発電機の 発振条件及び実験
ナイキストの安定判別に基づく熱音響システムの自励発振解析における発振余裕と 定常発振状態における圧力振幅の関係
電力フィードバック回路の調整による 熱音響発電機の発振余裕の最大化
フィードバック制御に基づく 熱音響発電システムの検討
低剛性・高慣性比の二慣性系の 外乱抑制制御問題に対して 任意の制御性能を達成する 不安定な補償器
落下水膜の振動特性に関する実験的研究 3m 理工学研究科   中村 亮.
34 PICマイコンを用いた能動騒音制御系の製作
振動体の振幅を一定とする 振動発電機負荷のフィードバック制御系の 安定性解析 長岡技術科学大学 ○ 永井 和貴 稲田 千翔之 小林 泰秀
タンク内圧力の変動を考慮した コンプレッサーの能動騒音制御
高慣性比二慣性系の外乱抑制問題に対する 慣性比 の解析解に基づく 補償器の構成
PI補償器の出力を時変係数とする 定常発振制御系の安定性解析
第 5 章 :周波数応答 5.1 周波数応答と伝達関数 周波数伝達関数,ゲイン,位相 キーワード : 5.2 ベクトル軌跡 ベクトル軌跡
音響伝達特性を用いたシングルチャネル音源方向推定
物理学実験 II ブラウン運動 ー 第2日目 ー 電気力学結合系の特性評価 物理学実験II (ブラウン運動) 説明資料.
KAGRA用防振装置のプレアイソレータの性能測定
31 ループ管熱音響システムにおける管内圧力計測系の製作 機械創造工学課程 梅本康平 担当教員 小林泰秀 准教授
23 造波機構における水位計の製作 1 はじめに 4 再現性の低下要因の実験 水位計の再現性の向上を目的としている. 2 実験装置
振動体の振幅を一定とする 振動発電機負荷のフィードバック制御 長岡技術科学大学 ○ 永井 和貴 齋藤 浄 小林 泰秀
長岡技術科学大学 大学院 工学研究科 機械創造工学専攻 髙山 誠 指導教員 小林 泰秀 准教授
電気回路学Ⅱ 通信工学コース 5セメ 山田 博仁.
臨界温度比推定のために熱音響エンジンを 定常発振させる時変ゲインを用いた 定エネルギー制御系の安定性解析
各種荷重を受ける 中空押出形成材の構造最適化
北大MMCセミナー 第94回 附属社会創造数学センター主催 Date: 2019年1月25日(金) 16:30~18:00
Presentation transcript:

水平板を用いた消波機構における指向性 アクチュエータの境界要素法による性能解析 0804   水平板を用いた消波機構における指向性  アクチュエータの境界要素法による性能解析 Analysis of directional wave actuator in active water-wave control system using horizontal plates based on boundary element method ○ 山本 龍彦(長岡技科大) 小林 泰秀(長岡技科大)

背景 性能改善 消波機構の用途:防波堤 外乱 外乱 反射 反射 透過 外乱 水平板タイプ 堤防タイプ (潜堤) ・消波機構の用途として防波堤が考えられる. ・堤防タイプがあり,外乱を100%反射することができるが左と右の水の循環に問題がある ・一方外乱を100%反射することはできないが,左と右の水の循環を考え,環境に配慮した堤防として,水平板タイプが提案されている ・水平板タイプの性能改善を目的に,鉛直方向に振動させるタイプが提案されている. <質問> ・防波堤に設置する場合に水平板の大きさはどのくらいになると考えられるか?また巨大になる場合に,それを動かすことは現実的に可能か? 外乱 鉛直方向に振動する水平板に対してロバスト制御系を設計し,外乱抑制を行う制御系を構成した.            [廣田,2010]

背景 従来研究 外乱 実験では示されているが,解析では示されていない. そのため,今回は解析を行う. Reference sensor  指向性アクチュエータを用いて外乱抑制を行う 制御系を構成した.                                                 [吉岡,2011] 外乱 Control actuator1 Control actuator2 指向性アクチュエータの効果 ・港湾内の水が港湾の形状によって共振することが知られている. ・従来研究として,この港湾内の共振を抑えるために,アクチュエータを用いた消波機構の導入を考えている. ・外乱をリファレンスセンサで測定し,鉛直方向に振動するアクチュエータ1つで消波する.右側(下流域)ではこれにより消波できる しかし,アクチュエータは左側(上流側)にも波を発生させてしまい,この波がリファレンスセンサで検出されてしまう. この問題を解決するために,制御用アクチュエータ2を導入に,リファレンスセンサに直接伝播する波を相殺するように駆動させることが提案されている. ・ここでは,このアクチュエータを指向性アクチュエータと呼び,制御用アクチュエータで生成した波がリファレンスセンサで検出されるまでの時間を延ばすことができる. ・この効果は,周波数応答結果で確認することができ,この時 制御用アクチュエータ1からリファレンスセンサまでの伝達関数である,フィードバックパス伝達関数に位相遅れが生じる. ・この研究は実験的に示されており,解析的には示されていない.そのため,本研究では解析的に検証する.  制御用アクチュエータからリファレンスセンサまでの距離が長い時に,良い閉ループ性能が得られる. 周波数応答  制御用アクチュエータ1からリファレンスセンサまでの伝達関数 (フィードバックパス伝達関数)に位相遅れが生じる. 実験では示されているが,解析では示されていない. そのため,今回は解析を行う.

解析方法 境界要素法 有限要素法 境界要素法 本研究では を使用する. 特徴 特徴 今回の研究では  2次元水槽を用いて,水中の障害物における波の反射率を境界要素法で解析した.              [R. Shin,et al, 2005]  2次元で非線形な非定常波の自由水面を有限要素法で解析した.             [G.X.Wu & R.Eatock Taylor ,1994]  3次元タンク内のスロッシング波を有限要素法で解析した.           [G.X.Wu, et al. 1998] 特徴 ・ 計算領域の境界から解析を行う. ・ 境界のみを解析するので,有限要素法よりも解析時間が短い. ・ 境界の情報だけを知りたい時に  使用される. ・解析方法として,有限要素法と境界要素法が考えられる. ・有限要素法では,非定常波の自由水面やスロッシング波の解析に使用されている. ・FEMの特徴として,計算領域を三角形要素に分割して解析を行い,3次元の複雑な構造物を対象とする時に使用される. ・境界要素法では,水中の障害物における反射率の解析に使用されている. ・BEMの特徴として,計算領域の境界から解析を行い,境界の情報だけを知りたい時に使用される. 今回の研究では, ・水面(境界)の情報だけを知りたい ・解析時間を短くしたい この2点から,境界要素法を用いて解析を行う. <質問> ・スロッシング波の解析は,境界要素法でも可能か? ・スロッシング波の消波には,導入できる可能性はあるのか? ・ 計算領域を三角形要素に分割して解析を行う. ・ 3次元の複雑な構造物における流体の解析に使用される. 今回の研究では ・水面(境界)の情報だけを知りたい. ・解析時間を短くしたい. 本研究では 境界要素法 を使用する.

目的 指向性アクチュエータの効果を境界要素法で検証する. 1.実験と解析による周波数応答結果を比較する.  実験と解析,無指向性と指向性アクチュエータの比較を行う. ・目的として,指向性アクチュエータの効果を境界要素法で検証する ・BEMの解析と実験結果の比較を行う.具体的には両者の周波数応答を示し,ゲイン線図と位相線図で比較する. ・指向性アクチュエータの効果を解析的に検証する.具体的には,フィードバックパス伝達関数の位相遅れの増加分を検証する. 2.指向性アクチュエータの効果を解析的に検証する.  指向性アクチュエータによるフィードバックパス伝達関数の位相遅れの 増加分を検証する.

実験装置の概要 ・港湾内の共振現象を模擬するために,左側のアクチュエータで外乱を生成し,右端のエラーセンサで消波することを考えます. ・外乱をリファレンスセンサで測定し,この2つの制御用アクチュエータで消波を行います. ・この時,水槽の左側を上流方向,右側を下流方向とします. ・制御用アクチュエータの駆動方法は2種類あり,制御用アクチュエータ1のみを駆動させる場合をCase Aとします. ・Case Aに対して,リファレンスセンサに直接伝播する波を相殺するために,むだ時間を設定した制御用アクチュエータ2を設定しめす.これがCase Bになります. ・ここで制御用アクチュエータ2で設定するむだ時間は,波の速度から計算することができます. 音波の場合では速度が周波数に依存しないため,全周波数帯域で指向性のある波を生成できます. しかし,水の波の場合では速度が周波数に依存してしまうため,1つの周波数に対して指向性のある波を生成することができます. ・本研究ではこの周波数を3次共振周波数と設定します. <質問> むだ時間を設定する時に,3次共振周波数を選んだ理由はなぜか? 上流方向 下流方向

境界要素法による数値解析 境界要素法による数値モデル 要素の幅:ΔS=5として,実験装置を数値モデル化した. Control actuator1 先ほど説明した実験装置を数値モデル化した結果です. ・水槽の外側・3つのアクチュエータをそれぞれ境界として,設定しました. Disturbance actuator Control actuator2 境界要素法による数値モデル

解析モデル ・ラプラス方程式 ・自由表面の設定 ・水槽の水底・壁面,水平板(静止)の境界条件 ・水平板(振動)の境界条件 o x z ・ラプラス方程式 ・自由表面の設定 →エネルギー散逸を考慮した係数αを追加 :速度ポテンシャル ・水槽の水底・壁面,水平板(静止)の境界条件 ・数値モデルの境界条件を設定します. ・境界内の条件をラプラス方程式と仮定します ・自由表面は,エネルギー散逸を考慮したαを設定しました ・水槽の水底・壁面,また静止している水平板の境界条件をこのように設定しました ・水平板(振動)の境界条件 ・速度ポテンシャルと位置zにおける変位の設定 ・周波数応答関数の定義 →水平板の変位から位置xにおける水面の変位まで周波数応答G(x, jω)を定義  周波数応答関数を求めるために,    を数値的に求める.

アクチュエータの駆動方法 Case A 制御用アクチュエータ1のみを使用する(v(t)=0). 無指向性アクチュエータ 制御用アクチュエータ1のみを使用する(v(t)=0).  制御用アクチュエータ1,2を使用し,任意の周波数で 指向性の波を生成する(v(t)=-u(t-τ)).  τはむだ時間であり,3次共振周波数の波の速度Cpと 制御用アクチュエータ間の距離dより求める. τ=d/Cp=0.356[s] Case B       指向性アクチュエータ  全周波数帯域で指向性の波を生成する      (v(jω)=-Θ(jω)u(jω), Θ(jω)=exp(-jkLd)).  Ldは制御用アクチュエータ1,2間の距離,kは波数であり, 角周波数ωと波の速度C(ω)より求める. k=ω/C(ω)       Case C       厳密な指向性アクチュエータ ・CaseCの説明(式などを加える) →実現させることは難しいが,全周波数で単方向波を生成できるアクチュエータを仮想的に考える. ・比較のために(制アク1→壁との距離の2倍)を考慮したCaseA→CaseA2とする ・CaseA2とCaseCが非常に合っていることがわかる.これより,CaseCが制アク1→壁→リファレンスセンサに検出されたことを意味する      Case A2        むだ時間を考慮した    無指向性アクチュエータ  Case Aにむだ時間(exp(jkl))を加える.  kは波数,lは制御用アクチュエータ1からエラーセンサ側の壁までの距離の2倍とする. l=315×2=630 mm Control actuator1 むだ時間 Ref. sensor Case A Wall Case A2

周波数応答結果(実験と解析の比較) おおむね一致 Case AよりもCase Bの位相が遅れている 解析結果 実験結果 ・実験と解析による周波数応答結果を比較します. ・図の説明:実験結果と解析結果,横軸・縦軸の説明,Case AとBの説明,縦軸の説明(共振周波数) ・実験結果について,2次共振周波数より低周波側ではノイズの影響があるため,解析との比較が難しい ・実験結果について,3次共振周波数より高周波側では細かな波が発生し,共振のピークが確認できない ・そのため,2次と3次共振周波数付近を比較する. ・ゲイン線図では,2次と3次共振周波数が一致していること,この共振周波数付近のゲイン特性が似ていることがわかる ・位相線図では,設定した3次共振周波数付近でCase AよりもCase Bの位相が遅れていることが確認できる <質問> ・2次と3次共振周波数周辺以外で,解析結果は使えないのか? ⇒実験結果から,3次共振周波数より高周波側で,ゲインが小さい周波数が確認できる.解析結果も同様に確認できる これより少ないが確認できる部分もある. Case AよりもCase Bの位相が遅れている 解析結果 実験結果

指向性アクチュエータの効果の検証 Case Cについて 位相特性を比較すると, Case A2とCase Cが ほぼ一致している.  位相特性を比較すると,  Case A2とCase Cが    ほぼ一致している. ・Case A2とCase Cを比較する. ・位相線図に注目すると,両者がほぼ一致している. ・Case Cは,制御用アクチュエータ1からリファレンスセンサに直接伝播するのではなく,壁で反射されてからリファレンスセンサで検出されていることがわかる. Case Cについて Control actuator1 Control actuator2 Ref. sensor Wall

指向性アクチュエータの効果の検証 指向性アクチュエータとして利用できる ・CaseBとCを比較する. ・CaseBとCを比較し,2.0Hzまで一致していることがわかる.これより,設定した3次共振周波数付近では指向性の波を生成できていることがわかる. (補足情報) それ以降は,直接上流側へ放射される進行波を相殺できなくなり、CaseCに対して位相が進んだ。 位相が進んだ理由としては,CaseBで直接上流側へ放射される進行波を相殺できなくなったためと考えられる. 指向性アクチュエータとして利用できる

結言 1.実験と解析による周波数応答結果を比較した. 2.指向性アクチュエータがフィードバックパス 伝達関数に与える効果を検証した.   実験結果と解析結果を比較し,共振周波数の一部と 共振周波数付近のゲイン特性がおおむね一致した. 2.指向性アクチュエータがフィードバックパス 伝達関数に与える効果を検証した. ・指向性アクチュエータは,むだ時間を設定した周波数付近で厳密な指向性アクチュエータとおおむね一致した. ・指向性アクチュエータを用いることで,フィードバックパス伝達関数における位相遅れが増加する.その増加分は,制御用アクチュエータ1とエラーセンサ側の壁までの距離の2倍に相当する. 以上より,指向性アクチュエータの効果を解析により示した.