• 検索結果がありません。

平時・出水時において 橋脚の定点で観測

N/A
N/A
Protected

Academic year: 2022

シェア "平時・出水時において 橋脚の定点で観測"

Copied!
3
0
0

読み込み中.... (全文を見る)

全文

(1)時系列振動解析を用いた河道内状態推定に関する研究 群馬大学大学院 学生会員 藤田智之 群馬大学工学部建設工学科 正会員 松本健作 群馬大学工学部建設工学科 非会員 神戸大典 株式会社数理設計研究所 非会員. 菅正信. ・ 研究目的 出水中は大量の激流が流れており、河道内に立ち入って計測等を行うことが非常に困難であるが、河 川工学上、出水中における河道内の現象を把握することは非常に重要であると考えられる。 本研究は時系列振動解析技術を、流体・河床および河川構造物からなる連成系に適用し、橋脚その他 河川構造物に設置した振動計によって記録される加速度信号から出水中における種々の水理諸量を推定 し、防災の観点からみた河川環境の健全性を診断するシステムの開発を試みるものである。. ・ 研究内容 本研究では自動制御の分野で用いられる定常時系列の概念を 物理モデルの作成. 平時・出水時において 橋脚の定点で観測. 運動方程式. 用いて線形系の応答の予測を行い、その結果から系のパラメー タを推定する手法を用いる。河川構造物の動力学モデルを作成. 連続時間システム. 観測結果より得られた自己 相関関数でYule-Walker方 程式を解く. インパルス応答 離散時間モデル. し、そこから離散時間モデルの理論的導出を行う。そして力学 パラメータと時系列パラメータの対応関係導出を行う。. 伝達関数の完成. ARパラメータの推定結果. 力学パラメータと時 系列パラメータの対 応関係を導出. Yule-Walker 方程式を経由することにより自己回帰・移動平均. 橋脚の 固有振動数 減衰定数. さらに、観測した加速度信号の自己相関関数を用いて. と. ARパラメータの関係式 を求める。. 過程の係数を推定する。この係数と構造物の固有振動数と減衰 定数の間に成立する関係式を用いることにより、橋梁の振動観 測データからリアルタイムで水位や流量、河道内地形などの水. 微動の観測データから 直接、構造系の振動特 性を推定. 理諸量を直接推定しようとするものである。ここでは最初の段 階として振動特性と水位との相関性ついて検討を行った。. ・ 現在までの成果 ① 条件の異なる観測結果より求められる自己相関関数とパワースペクトルの比較 栃木県葉鹿橋において任意に設置箇所を選び、加速度センサーを順次取り付けて、それぞれの箇 所について、横断、流下、鉛直方向と三つの方向を設定し、それぞれの加速度のデータを得た。下 図に示すのは橋梁の同じ地点において観測条件によって最も顕著に違いの現れた結果である。また 観測は7月11日、8月9日のものであり、水位は国土交通省から公表になっているデータよりそ れぞれ1.25m、3.85mとなっている。7月11日は平時の状態、8月9日は台風 10 号によ る出水時の状態でそれぞれ同じ地点で観測を行った。図1は平時(7月11日) 、出水時(8月9日) の観測データで、上から横断、流下鉛直方向を表しており、原信号とその自己相関関数を用いて雑 音を除去したもののそれぞれの比較である。そのサンプリング数はおよそ 1190Hzである。また図 2は原信号にフーリエ変換を施しそれをパワースペクトルに変換したものの平時と出水時の比較で ある。特に相違が見られたのが流下方向であり平時、出水時によって自己相関ではっきりと違いが 見られた。またその周期は T=0.04(s)となっている。. キーワード:時系列振動解析、河道、状態推定 連絡先:〒376-0041 群馬県桐生市天神町 1-5-1 群馬大学工学部建設工学科 TEL 0277-30-1640.

(2) 原信号(上)と自己相関関数(下) 方 向. 信 号. 横 断. 元 信 号. 平 時 (7 月 1 1 日 ). パワースペクトル比較. 出 水 時 (8 月 9 日 ). 平時(7月11日) 横 断. 自 己 相 関. 流 下. 元 信 号. 流 下. 自 己 相 関. 鉛 直. 出水時(8月9日). 鉛 直. 元 信 号. 自 己 相 関. 図 1. 図2. ② 物理モデルの作成から力学パラメータと時系列パラメータの対応関係導出まで ①の結果をうけて振動に関連した出水によるモードの変化を確認したため河川 m. 構造物の動力学モデルを作成した。まず河川構造物のモデルを簡略化するため、流. x. 下方向、鉛直方向に分け考えることにした。力学モデルは地盤、橋脚からなる 1 自. f c. k. 由度系とし、左図のように表現した。左図は鉛直方向の図である。ここで記号 m,c,k は質量、粘性減衰係数、剛性を意味した力学パラメータと呼ばれるものである。 出水中の河川構造物は要素が複雑であるため振動系への洗掘力が未知である。よ. って本研究では未知である橋脚に対する入力f(t)を白色雑音とし、橋梁に取り付けた加速度計から観測 される振動加速度を出力とした加工モデルの定式化を行った。同定対象の物理法則を表す微分方程式を ラプラス変換し連続時間システムの伝達関数 H(s)を求める。伝達関数 H(s)を逆ラプラス変換して連続時 間システムのインパルス応答をh(t)を求め、それを離散時間化し離散時間システムのインパルス応答h (n)を求める。h(n)を Z 変換して離散時間システムの伝達関数 H(z)を求める。ここで求めた伝達関数が ARMA モデルの伝達関数 H(z)= B(z) / A(z) となる。以上を踏まえ本研究で求められたのが①式である。 H. =. ( z ). 1 −. 1 − a 1 Z. b − 1. 1. Z −. ―――①. − 1. a. Z. 2. −. 2. ここで AR パラメータと力学パラメータ m,k,c の関係は以下のようになる。 a a. 1 2. = =. φ. ( 1 ) 1. − φ. + ( 1 ). 1. φ. 2. •. φ. φ. ( 1 ). ( 1. ). =. e. ( 1. ). =. e. 1. −. c. −. c. ( 1 ) 2. φ. 2. ±. 2. c 2. . 2. c 2. −. 4. mk. −. 4. mk. m. m. T. T. ・ おわりに 物理モデルの作成から力学パラメータと時系列パラメータの対応関係導出まで求めることができた。 これに関する精度検証を AIC(赤池情報基準)を計算し行いたい。平時・出水時ともに様々な振動特性を示 す事から、水位(流量)の他にある風速・気温・湿度など振動に影響を与えうる因子が多く存在するも のと思われる。それら河川構造物の振動に影響する雑音を消去し、出水のみによる橋梁への振動特性へ の関連性を検討していきたい。現在、振動特性より水位の推定を行っているがこの振動特性より次は流 量を推定し、構造物の状態診断の検討を進めていきたい。その結果に関しては講演時に発表したい。 <参考文献> 1)中溝高好:信号解析とシステム同定:コロナ社:1988 年3月 2)大類重範:ディジタル信号処理:日本理工会出版:2001 年 12 月.

(3) 3)土岐憲三・佐藤忠信・生木泰秀:自己回帰・移動平均法による構造物の動的パラメータの推定:京 都大学防災研究所年報第 21 号B-2:昭 53 年 4 月.

(4)

参照

関連したドキュメント

本研究では.観測点を定めると,それぞれが受けもつ 領域が汚染発生分布とは独立に定まると

フィリピン付近の地震による津波を過去日本で、観測さ れたものは第

  阪神高速道路では,ASR橋脚の  劣化度の判定は,ASR構造物の維 持管理マニュアル 1)  に基づき,外観

次に、表3に天文観測から観測ロケットへの一般的な要求をまとめた。最も重要な要求は指 向精度要求である。ほとんどの天文観測において、特定天体

3.1 捕獲した Malware の種類と感染ポート 収集した Malware を ClamAV でウィルススキャン した結果,表 1 に示すように,ほとんどがボット関 連の

1.研究成果の概要 地下1000メートルの東京大学宇宙線研究所神岡地下観測所内にラドン観測網を設

153

観測画像修復における因果および半因果モデルの比較 98